空间飞行器综合定轨与参数分析软件COMPASS的研发——动力学统计定轨程序compass的编写及调试

介绍了空间飞行器综合定轨与参数分析软件COMPASS的开发过程。软件的初级阶段目标是可以利用SLR观测对多颗激光卫星进行同时定轨、可以利用非差GPS观测对GPS星座进行同时定轨,并估计有意叉的地学参数。COMPASS的开发采取了由简到繁、循序渐进的技术策略,软件开发经历了这样几个主要过程：多星多技术定轨框架的建立。利用SLR观测确定GPS卫星的轨道,利用IGS的SP3轨道确定GPS星座的轨道,利用非差GPS伪距观测确定GPS星座的轨道,利用非差GPS伪距和相位观测确定GPS星座的轨道。激光卫星的定轨精度已经达到国际水平,可以用于提供国际服务（如IERSEOP;ILRS快速分析）;GPS定轨内符精度达到国际先进水平,平均外符精度好于30cm。     

面向导航增强的低轨卫星钟差确定及预报方法研究

低轨导航增强是未来导航发展的重要趋势,而高精度低轨卫星钟差是实现低轨导航增强的必要条件。基于Sentinel-6A卫星,对低轨卫星钟差特性进行了分析,给出了钟差确定方法及影响因素,介绍了顾及钟差特性的低轨卫星钟差预报方法。实验表明,低轨卫星钟差含有多个周期项,给低轨卫星建模和预报带来了困难。与使用运动学定轨模型相比,基于简化动力学的定轨模型可显著提升低轨卫星钟差精度;当基于运动学模型确定低轨卫星钟差时,相较于使用GPS单系统数据,多GNSS观测数据可提升低轨卫星钟差精度。研究表明,基于GPS和Galileo观测的Sentinel-6A卫星钟差精度相较于GPS单系统钟差精度改善了36%,同时,所使用的GNSS产品精度与低轨卫星钟差精度密切相关。利用顾及卫星钟差特性的低轨卫星钟差预报方法,当预报时长小于1 min,低轨卫星钟差预报精度（预报与解算值之差的RMSE）在0.1 ns之内,当预报时长小于5 min,预报精度在0.3 ns之内,随着预报时长的增长,预报精度显著下降。     

利用DORIS测轨系统实现高精度定轨

仿真研究DORIS测轨系统,重点考察了大气密度模型误差、测量精度、测轨网分布对定轨精度的影响。仿真结果表明,除了测轨网的地理分布,动力学模型中的大气模型误差对中低轨卫星精密定轨结果影响也较大。对ENVISAT卫星的DORIS实测数据进行了定轨分析,结果表明实测数据的定轨精度比仿真精度大约低1个量级不到。综合仿真结果和实测数据进行精度分析,推断对800km高度的太阳卫星轨道,采用8个DORIS信标站布设,24h定轨,定轨三维位置精度可以达到29cm(1σ),其中径向误差为3.4cm(1σ)。若采用30个DORIS信标站布设,定轨精度可提高30%。     

导航卫星的短弧运动学定轨

导航卫星在姿轨控和轨道恢复期间,由于观测数据有限,传统的统计定轨理论难以实现导航卫星精密定轨。本文尝试采用一种不依赖轨道动力学的、新的运动学定轨方法来处理短弧和复杂动力学过程中的定轨,提出了基于多项式拟合的短弧运动学定轨算法,并提出2种不同的实现方案。该算法充分利用了高采样率的测轨数据,减少了结果的噪声,其优点在于不需要长时间累积测轨数据,可以实现近实时快速计算,克服了动力学法定轨发散和单点定位无法获得卫星速度信息等不足。对COMPASS M-01导航卫星实测数据的处理表明,10min左右短弧运动学定轨的位置精度可以优于10m,速度精度优于4cm/s,满足了短弧跟踪条件下RDSS对卫星轨道精度的要求,实现了短弧跟踪条件下卫星精密定轨,但从轨道预报精度来看,该方法仅仅适用于短期预报。     

基于双星系统的卫星几何学精密定轨方法研究

双星导航系统的建立，为我国实现自主的卫星天基测控技术提供了可能。同时由于其全天候、全天时、相对定轨精度高等特点可望成为一种低、中轨卫星精密定轨的有效手段。但由于双星系统中静止卫星数目仅有两颗，要得到低、中轨卫星的三维定位信息，还需要额外的第三维观测量。本文根据双星系统的现状，结合国内对低、中轨卫星精密定轨的要求，提出了一种基于双星系统的卫星几何学精密定轨的方法，概述了其基本原理及其组成，给出了卫星定轨的方法和数学模型，同时根据数据仿真分析了该定轨系统的定轨精度，为下一步工程实现提供了理论基础。     

广播星历旋转误差对低轨星载BDS-3/GPS实时精密定轨影响分析

基于LT-01A卫星星载BDS-3/GPS观测值进行了星载实时精密定轨研究,并重点分析了广播星历旋转误差对实时定轨精度的影响。通过赫尔默特转换评估了所选时段内GPS和BDS-3广播星历轨道旋转误差,显示BDS-3广播星历旋转误差可达-8.7 mas,平均量级较GPS大约2.5倍。BDS-3广播星历经旋转改正后,轨道切向、法向均方根(RMS)误差从25 cm左右提升至10 cm量级,提升幅度超过50%。因此,基于星载BDS-3以及BDS-3/GPS联合的实时定轨精度受BDS-3星历旋转误差影响严重,且主要作用于切向和法向。经过旋转改正后,单独BDS-3实时定轨在切向、法向、径向RMS分别为21.0 cm、10.7 cm及11.2 cm,其切向和法向精度比改正前分别提升15.0%和31.8%;BDS-3与GPS联合定轨进一步提升切向精度至19.4 cm。得益于BDS-3广播星历较高的精度,单BDS-3以及BDS-3/GPS联合的实时定轨在旋转改正前的三维RMS分别为31.9 cm和29.2 cm,较单GPS实时定轨分别提升9.1%和16.8%;添加旋转改正后,其定轨精度分别提升至26.7 ...     

双线程滑动窗口多系统GNSS超快精密定轨实现及评估

全球导航卫星系统（GNSS）超快精密定轨为GNSS实时应用提供了高精度空间基准。基于天地协同定位、导航与授时（PNT）网络服务中心实现了四系统GNSS卫星超快精密定轨,并对定轨结果进行精度评价。介绍了天地协同PNT网络的概念内涵以及网络服务中心部署的超快精密定轨软件架构和详细功能,并针对实时应用需求提出了一种双线程滑动窗口超快精密定轨策略。最后利用重叠弧段比较、与外部轨道产品比较以及卫星激光测距（SLR）检核3种方式对定轨结果进行了精度评价。结果表明,与武汉大学分析中心的最终事后精密轨道产品相比,四系统GNSS MEO卫星预报6 h弧段的径向均方根（RMS）误差整体在2~5 cm水平,BDS2 IGSO卫星最小一维RMS误差在10~15 cm水平;GPS和Galileo卫星的SLR检核残差均值在1~3 cm水平,标准差在3~6 cm水平,能够满足后续厘米级实时应用对空间基准的精度需求。     

一种数值定轨方案

本文给出一种有别于习惯采用的人造卫星精密定轨方法，该方法仍以分析定轨法的轨道根数作为待估状态量，但却采用数值法计算精度要求较高的卫星受摄星历，而用分析法去计算条件方程中所需的状态转移矩阵。这样既保持了清晰的轨道图象，又使计算卫星受摄星历和状态转移矩阵得以简化，相应的软件形式简单。只要定轨弧段Ｓ不太长（Ｓ〈１０＾３），就能在一定精度条件下明显提高卫星定轨的计算效率。     

低轨卫星对高轨卫星仅测角初轨计算方法

针对我国地面测站对高轨卫星监视能力缺乏的问题,提出一种低轨卫星对高轨卫星仅测角初轨计算方法。该算法引入天基跟踪坐标系,消除测距信息影响,建立仅测角观测方程;引入法向运动,增加摄动因素影响,建立扩展拉普拉斯动力学模型;推导分析观测模型和动力学模型的关系方程,将初轨计算问题转换为非线性方程求解问题,利用高斯全主元消去法完成方程求解。通过实战和仿真测角数据对方法进行检验,结果表明,该方法能利用仅测角数据对非合作目标进行初轨确定,精度在公里量级,可为我国地基监视系统提供补充参考。     

GPS单差观测量的运动学相对定轨研究

以双星编队飞行为工程背景,在两颗编队卫星均搭载GPS双频接收机和一颗搭载GPS双频接收机而另一颗搭载GPS单频接收机两种模式下,对相对定轨方法和数据处理流程进行了研究。运动学相对定轨采用最小二乘批处理方法进行轨道估计,定轨流程包括单点定位、动力学平滑、数据编辑、运动学绝对定轨和运动学相对定轨五个步骤。利用GRACE双星数据进行测试,结果通过KBR验证:利用消电离单差组合观测量,相对定轨精度达到1 cm左右;以双频GPS为参考星,双频消电离观测量和单频GRAPHIC观测量做单差,相对定轨精度优于18cm。     

星载双频GPS相位历元差缩减动力学定轨(英文)

研究了星载双频GPS相位历元差缩减动力学定轨方法,该方法既克服了相位历元差运动学方法在观测几何较差或数据不足情况下无法应用的缺点,又避免了相位非差动力学方法易受周跳和模糊度影响的缺点。历元差方法对相位周跳的影响不敏感,因此能够降低相位数据预处理中周跳探测的难度。在模型求解过程中,解决了高维矩阵的计算问题,将长弧段观测数据分成若干短弧段,相邻两个短弧段的连接处不做差分,仅在每个短弧段内部进行历元间差分处理。通过对GRACE卫星进行试算,并与GFZ事后科学轨道进行比较,结果表明相位历元差缩减动力学定轨精度在径向、沿迹方向和轨道法向分别可达1.92cm、3.83cm和3.80cm,三维位置精度可达5.76cm,该方法与相位非差缩减动力学定轨精度相当。     

基于改进UKF的小卫星磁测定轨方法

针对小卫星设计简单、功耗低的特点,提出一种利用三轴磁强计测量信息实现卫星磁测自主定轨的方法。为了提高自主定轨的计算效率,采用基于球体单形ε点采样的改进UKF(Unscented Kalman Filter)作为处理自主定轨的滤波算法。最后对这种定轨方法进行了仿真试验,结果表明,该滤波算法能有效保证收敛的速度和定轨精度。     

空间飞行器的多普勒实时自主定轨精度分析

空间飞行器需要实时的高精度轨道信息来完成对栽荷的指令操作和遥感数据的实时处理。除了星栽GPS技术,星载多普勒无线电定轨定位系统（DORIS,Doppler Orbitography and Radio—positioning Intergrated by Satellite）是仅有的有能力提供分米级精度的实时在轨轨道确定技术,它可通过测量星地相对多普勒频移,在星上完成实时定轨和预报,目前该技术已在国外多个卫星上实现,达到了较好的效果,而我国还没有建立这样实时自主定轨系统。为此,结合我国高分辨率空间对地观测系统的建设需求和我国航天器对实时自主定轨及其精度的要求,利用扩展卡尔曼滤波算法对多普勒测量进行了实时自主定轨仿真计算,分析了频率偏差估计与否、初轨误差、地面信标站地理分布以及观测精度等对实时自主定轨的滤波收敛时间和定轨精度的影响,为我国利用DORIS技术进行实时在轨轨道确定提供方案和软件原型。仿真计算表明,基于28个全球分布的地面站,对于高度为800km的卫星,在忽略其动力学模型误差的假设下,若初轨三维位置、速度误差分别为100m（或差至1km）、1m／S（1d）,2h后滤波可以达到稳定收敛,收敛后的实时定轨误差可以达到0．1m（1d）。滤波估计参数除了6个卫星轨道状态参数,还估计了地面信标相对于卫星超稳定振荡器的频率偏差;     

基于双星定位系统和星敏感器的低轨卫星自主定轨方法研究

利用双星定位系统的两个测距信息外加星敏感器的测角信息,基于信息融合和信息守恒理论,采用联合推广卡尔曼滤波算法,对低轨卫星进行自主定轨仿真研究。计算结果表明,该方法可以显著提高定轨精度,达到低轨卫星自主定轨的精度要求。     

卫星定轨综述

对卫星定轨中的热点问题作一综述,内容包括短弧定轨和长弧精密定轨。对于短弧定轨(通常指初轨),不仅从原理上分析其关键之处,还介绍如何将适用于无摄运动的特殊定轨方法推广用于同时考虑各种摄动的受摄运动,并可同时适用于圆锥曲线的各种情况：变化的椭圆轨道、抛物线轨道和双曲线轨道,后者在深空探测器的运动过程中将会遇到。对于长弧精密定轨,阐述当今提高精密定轨精度的力模型补偿问题以及几何定轨的提法,天地基网联合定轨和星-星相对测量的自主定轨问题。     

相对测量在编队卫星自主定轨中的应用

针对编队卫星自主定轨问题进行了研究,设计了一种完全不依赖于地面站和GPS系统的自主导航方案。利用星间测量信息进行卫星编队相对轨道状态的自主确定;并在利用磁强计进行卫星绝对轨道自主确定的基础上,引入星间测量信息提高绝对定轨精度;设计扩展卡尔曼滤波器进行卫星编队轨道状态估计,数学仿真结果验证了这种导航方案和算法的有效性。     

基于伪距信息的COMPASS-MEO导航卫星单星定轨分析

COMPASS-MEO导航卫星的伪距观测量包含星站距离、接收机钟差和卫星钟差以及各种噪声。本文首先分析了接收机钟差和卫星钟差在一定间隔内主要表现为线性变化的特性,可以考虑将接收机和卫星钟差作为测距的偏差和偏差变化率和轨道一起进行求解。其次,利用实测的MEO导航单星伪距数据,进行了数据预处理和定轨试验,分别对站间无时间同步信息和有时间同步信息两种情况进行定轨、残差分析和参数解算。使用重叠弧段、轨道预报和激光观测数据验证定轨结果的精度。结果表明,两种情况下的定轨结果无明显差别。该定轨方法可以运用于MEO导航单星精密轨道的计算。     

用于低轨卫星精密轨道确定的简化动力学法

全球定位系统(GPS)的测量量可促成一种简化动力学法，用于近地卫星的近厘末级定轨。利用这种方法，不同观测时刻的卫星状态间的转移由正规的动力学模型和由连续GPS载波相位数据以运动学法导出的卫星位置变化来完成。动力学信息和运动学信息间的相对权重可自由变化。协方差分析表明：当观测几何差而动力学模型好时，模型支配着状态转移的计算；动力模型差而观测几何好时，载波相位统治着状态转移的计算；而当运动学和动力学信息间没有明显优劣时，简化动力学联合法可显著提高轨道精度。本文给出简化动力学法中准最佳的择权准则，并研究了求解精度对权的敏感性。     

EUVE卫星的GPS测轨

1992年夏,随着极紫外探测卫星(EUVE)一起发射了一台单频GPS接收机。EUVE卫星轨道低,并在星体上直接安装了二个天线,可用来研究误差源,以及未来卫星可达到的定轨精度。利用GDPSY—OASLSⅡ软件处理了几段GPS教据。从EUVE来的数据和从全球GPS接机网来的数据相结合,消除了SA的影响。研究了几种消除单频数据中电离层延时误差以及利用简化动力学法求解轨道的方法。通过内符合检验并与哥达德航天中心产生的GTDS轨道解作比对,鉴定了GPS测轨精度。利用GPS的测轨精密度和准确度在一米量级。     

空间飞行器综合定轨与参数分析软件COMPASS的研发——结构设计

在空间飞行器综合定轨与参数分析软件（COMPASS）需求分析的基础上,重点针对定轨软件的总体流程、多卫星多传感器动力学定轨的流程、全球导航卫星系统（GNSS）数据处理流程进行结构设计,给出了详细的过程和数据流程图。这一工作为COMPASS的编写奠定了基础。