基于光学测角数据的风云四号同步轨道卫星精密定轨

针对风云四号同步卫星的精密定轨和精度评估需求，首先利用地面光学测角数据对FY-4A卫星进行精密定轨，定轨后方位角和高度角的残差rms分别为0.25＂和0.45＂。与基于测距数据的轨道相比，位置精度在有测角数据的弧段内小于50m。进一步联合测角数据和测距数据对FY-4A卫星进行联合定轨，定轨后轨道重叠精度优于15m。利用联合定轨结果评估了基于测距数据的实时轨道产品精度，可以明显发现轨道精度随着测距数据的积累而逐步提高。     

空间引力波探测器轨道确定及仿真分析

针对空间引力波探测器的精密轨道确定问题,选取LISA(Laser Interferometer Space Antenna)作为研究对象,建立了探测器的目标仿真环境,模拟生成美国深空网(Deep Space Network,DSN)和中国深空网(China Deep Space Network,CDSN)在不同测站下的测距测速数据,采用非线性加权最小二乘(Weighted Least Squares,WLS)和蒙特卡洛方法(Monte Carlo,MC)分析跟踪弧长、测轨数据类型、测站数量及其分布等因素对LISA探测器定轨精度的影响。仿真结果表明:①增加测站的定轨弧长,可有效提高轨道确定精度,当跟踪弧段达到20 d时,探测器位置确定精度可达92 m;②在现有测量精度条件下,测距信息对定轨精度的贡献较大,测距/测速组合定轨方式与仅利用测距、测速定轨相比,估计的位置精度分别提高了32.23%、99.52%;③采用多测控网定轨模式可以提高多站共视比率,有利于提高定轨精度和收敛速度,DSN/CDSN联合定轨比仅采用DSN定轨的平均位置精度提高43.73%。     

“龙江2号”月球轨道微卫星定轨分析

微纳卫星深空探测任务中,通常所分配的测控资源有限,因此有必要对有限测控资源条件下微纳卫星的定轨精度进行分析。以微纳卫星深空探测为背景,采用"龙江2号"微卫星的轨道测量数据对其定轨精度进行了分析。"龙江2号"微卫星只有USB轨道测量数据,且环月段测控资源相对紧张,每天有两站跟踪,共约3～4 h的轨道测量数据。首先介绍了"龙江2号"微卫星飞行任务及其飞行过程中影响测定轨的因素;其次给出了定轨的动力学模型,对微卫星地月转移段的定轨精度进行了分析;最后通过分析摄动力、动量轮卸载以及数据弧段长度的影响,给出了微卫星环月阶段所使用的定轨策略,并通过重叠弧段比较的模式,给出了微卫星环月段的定轨精度。研究结论可以为后续微纳卫星深空探测任务提供有益参考。     

对一种月球与火星探测多程微波测量链路定轨定位的数值模拟初步分析

轨道器精密定轨与着陆器的精确定位在深空探测任务中具有非常重要的科学意义。对一种月球与火星探测多程微波测量链路的定轨定位能力进行了初步仿真分析,推导了这种多程微波测量链路的测量模型,分析了该模型的优势。模拟仿真分析结果表明,此测量跟踪模式的数据具有提升轨道精度的潜在能力,并且同时求得着陆器的位置。定量分析表明,在考虑坐标系转换误差,重力场误差,行星历表误差以及星上转发误差的情况下,模拟1 mm/s的噪声,对于月球探测器来说,轨道器的定轨精度可达几米,着陆器的定位精度有望达到分米量级;对于火星探测器来说,轨道器的定轨精度可达到数10 m,着陆器的定位精度可达到几米。     

基于固定模糊度PPP的低轨卫星实时精密定轨

传统动力学定轨法受制于动力学模型精度，传统几何定轨法精度受限，只能达到亚米级，而基于精密单点定位（PPP）模式的几何定轨法一般采用浮点解，定轨精度及可靠性较基于双差模式的相对定位较差。为提高PPP模式低轨定轨的定位性能，利用中国区域内外的IGS测站计算出当前所有卫星的宽巷和窄巷相位小数偏差产品，对经过中国大陆区域上空的国产低轨卫星海洋二号（HY-2）和资源三号 (ZY-3) 卫星进行固定模糊度PPP的定轨解算，与事后精密轨道结果进行比较，分析其外符合精度。结果表明：仅利用约10min弧段的HY-2和ZY-3卫星数据，切向与径向的定轨精度可达2cm左右，法向为5cm左右，较浮点解定轨精度大幅提升。基于固定模糊度PPP的定轨方法能够满足厘米级的实时精密定轨。     

USB观测数据时标偏差影响分析与评估

USB系统是目前中国载人航天和月球探测任务的主要测控网.由于USB测量设备本身以及无线电信号传播媒介以及其他误差因素的影响,USB测量数据中包含了各种误差,需要在定轨时对观测数据进行误差修正.通常,例行的USB测量误差修正包括对流层折射修正、电离层延迟修正和通道延迟修正,但对定轨过程中可能影响观测数据计算精度的时标偏差并未作修正.针对USB测距测速观测数据,详细研究了观测数据时标偏差对观测值计算精度的影响,分析了误差影响特性,建立了相应的误差修正模型,并通过与卫星星历偏差对USB测距测速观测值计算精度影响特征的比较,发现时标偏差对测量的影响与轨道沿迹误差对观测计算值的影响等效,这使得在定轨过程中分离时标偏差的难度较大.提出了基于星载GPS定位数据分离时标偏差的方法,并利用某次任务的实测数据,分离出了该次任务中USB测量的时标偏差.最后针对目前USB数据时标偏差影响和观测误差量级相当的情况,建议将目前的观测时标精度提高到优于0.1ms的水平,使得时标偏差的影响降低到比观测误差小一个量级的水平.     

火星探测VLBI测定轨技术

我国将于2020年首次发射由环绕器和着陆巡视器组成的火星探测器,火星探测器的跟踪及精密测定轨是完成工程任务和科学探测的基础。火星探测器的跟踪和测定轨,目前主要采用基于地面无线电测量的测距、测速和甚长基线干涉VLBI测角3种手段。主要针对VLBI技术予以介绍,主要内容为:△DOR型VLBI技术在国内外的应用情况、火星探测器VLBI测定轨技术分析、基于同波束VLBI的火星车定位技术、火星探测器VLBI观测等。这些内容对我国的火星探测器测定轨有重要的应用价值。     

基于单频星载GPS数据的低轨卫星精密定轨

为满足搭载单频GPS接收机低轨卫星的精密定轨需求以及深化单频定轨研究,文中解决了单频星载GPS数据的周跳探测问题,并利用“海洋二号”（HY-2A）卫星及“资源三号”（ZY-3）卫星的单频星载GPS实测数据采用两种方法确定了二者的简化动力学轨道,并通过观测值残差分析、与双频精密轨道比较、激光测卫数据检核等方法对所得轨道精度进行评定。结果表明,在不考虑电离层延迟影响的情况下,HY-2A卫星定轨精度为2~3dm,ZY-3卫星为1m左右;而采用半和改正组合消除电离层延迟一阶项影响后,二者定轨精度均显著提高,HY-2A卫星三维精度提高至1dm左右,ZY-3卫星提高至1~2dm。文章的研究成果表明,搭载单频GPS接收机的低轨卫星也可获得厘米级的定轨精度。     

嫦娥系列探月卫星无线电科学实验简介

在"嫦娥1/2/3号"系列探测任务中开展了行星无线电科学探测实验,这些工作包括:使用天文VLBI技术对探测器进行工程和科学测轨、定位观测,开环和闭环测速测距观测,基于微波观测重构和优化月球重力场模型,通过重力异常揭示质量瘤和撞击盆地。星载或器载主/被动雷达设备还用于探测月壤和月球内部结构。在"嫦娥2号"任务中,实现了对拉格朗日平动点利萨如轨道飞行的测控,以及对图塔蒂斯小行星的飞掠探测。在"嫦娥3号"任务中,实现了多通道开环3向相位测距和多普勒测速技术。该月球无线电相位测距技术LRPR作为一种新的空间测量技术,可以用于测定台站的位置、潮汐、星体的自转特性。还可以与月球激光测距技术LLR融合,监测月球动力学运动变化,并有潜力用于未来的火星探测任务。     

基于Gooding算法的天基光学目标跟踪定轨

开展了基于Gooding算法的400km天基平台光学目标监测的轨道确定研究，当测量误差为3”和6”时，分别对800，1500及36000km轨道高度目标进行初始轨道确定及轨道改进分析.仿真结果表明，利用400km轨道高度平台对800~1500km轨道高度目标进行初定轨，测量数据误差为3”~6”时，4~15min弧段的初定轨精度约在10km量级，1~2min弧段的初定轨精度约在100km量级；15min初定轨弧段轨道改进后误差在100m量级，弧段小于10min时轨道改进误差精度在km量级.利用400km轨道高度平台对36000km轨道高度目标进行初定轨，测量数据误差为3”时，15~20min弧段的初定轨精度约在数十km量级，8~10min弧段的初定轨精度在100km量级；轨道改进后误差在km量级.测量数据误差为6”时，20min弧段初定轨精度在10km量级，8~15min弧段初定轨精度在100km量级，轨道改进后误差精度在10km量级.      

地基行星大气掩星观测资料流程和数据整理

在中国火星探测萤火一号(YH-1)计划中, 包括了地基掩星观测反演火星大气的科研任务. 观测资料整理是反演流程的第一步. 本文描述了地基火星大气掩星观测处理软件系统的观测数据流程和观测资料整理模块,并详细介绍了观测资料整理模块的结构和功能, 其中包括时间系统转换、历表插值、坐标系变换、信号时延改正以及掩星平面建立. 利用行星数据系统公布的火星快车无线电科学数据和由SPICE得到的地球、火星历表以及火星快车的轨道数据, 结合本文的算法, 得到了一些实验结果.      

Precision orbit determination performance for CryoSat-2

In this paper we discuss our efforts to perform precision orbit determination (POD) of CryoSat-2 which depends on Doppler and satellite laser ranging tracking data. A dynamic orbit model is set-up and the residuals between the model and the tracking data is evaluated. The average r.m.s. of the 10?s averaged Doppler tracking pass residuals is approximately 0.39?mm/s; and the average of the laser tracking pass residuals becomes 1.42?cm. There are a number of other tests to verify the quality of the orbit solution, we compare our computed orbits against three independent external trajectories provided by the CNES. The CNES products are part of the CryoSat-2 products distributed by ESA. The radial differences of our solution relative to the CNES precision orbits shows an average r.m.s. of 1.25?cm between Jun-2010 and Apr-2017. The SIRAL altimeter crossover difference statistics demonstrate that the quality of our orbit solution is comparable to that of the POE solution computed by the CNES. In this paper we will discuss three important changes in our POD activities that have brought the orbit performance to this level. The improvements concern the way we implement temporal gravity accelerations observed by GRACE; the implementation of ITRF2014 coordinates and velocities for the DORIS beacons and the SLR tracking sites. We also discuss an adjustment of the SLR retroreflector position within the satellite reference frame. An unexpected result is that we find a systematic difference between the median of the 10 s Doppler tracking residuals which displays a statistically significant pattern in the South Atlantic Anomaly (SSA) area where the median of the velocity residuals varies in the range of ?0.15 to +0.15?mm/s.     

Geostationary orbit determination using SATRE

A new strategy of precise orbit determination (POD) for GEO (Geostationary Earth Orbit) satellite using SATRE (SAtellite Time and Ranging Equipment) is presented. Two observation modes are proposed and different channels of the same instruments are used to construct different observation modes, one mode receiving time signals from their own station and the other mode receiving time signals from each other for two stations called pairs of combined observations. Using data from such a tracking network in China, the results for both modes are compared. The precise orbit determination for the Sino-1 satellite using the data from 6 June 2005 to 13 June 2005 has been carried out in this work. The RMS (Root-Mean-Square) of observing residuals for 3-day solutions with the former mode is better than 9.1 cm. The RMS of observing residuals for 3-day solutions with the latter mode is better than 4.8 cm, much better than the former mode. Orbital overlapping (3-day orbit solution with 1-day orbit overlap) tests show that the RMS of the orbit difference for the former mode is 0.16 m in the radial direction, 0.53 m in the along-track direction, 0.97 m in the cross-track direction and 1.12 m in the 3-dimension position and the RMS of the orbit difference for the latter mode is 0.36 m in the radial direction, 0.89 m in the along-track direction, 1.18 m in the cross-track direction and 1.52 m in the 3-dimension position, almost the same as the former mode. All the experiments indicate that a meter-level accuracy of orbit determination for geostationary satellite is achievable.     

Designing planetary protection into the Mars Observer mission.

Planetary protection has been an important consideration during the process of designing the Mars Observer mission. It affected trajectory design of both the interplanetary transfer and the orbits at Mars; these in turn affected the observation strategies developed for the mission. The Project relied mainly on the strategy of collision avoidance to prevent contamination of Mars. Conservative estimates of spacecraft reliability and Martian atmosphere density were used to evaluate decisions concerning the interplanetary trajectory, the orbit insertion phase at Mars, and operations in orbit at Mars and afterwards. Changes in the trajectory design, especially in the orbit insertion phase, required a refinement of those estimates.     

Precise orbit determination using the batch filter based on particle filtering with genetic resampling approach

In this study, genetic resampling (GRS) approach is utilized for precise orbit determination (POD) using the batch filter based on particle filtering (PF). Two genetic operations, which are arithmetic crossover and residual mutation, are used for GRS of the batch filter based on PF (PF batch filter). For POD, Laser-ranging Precise Orbit Determination System (LPODS) and satellite laser ranging (SLR) observations of the CHAMP satellite are used. Monte Carlo trials for POD are performed by one hundred times. The characteristics of the POD results by PF batch filter with GRS are compared with those of a PF batch filter with minimum residual resampling (MRRS). The post-fit residual, 3D error by external orbit comparison, and POD repeatability are analyzed for orbit quality assessments. The POD results are externally checked by NASA JPL’s orbits using totally different software, measurements, and techniques. For post-fit residuals and 3D errors, both MRRS and GRS give accurate estimation results whose mean root mean square (RMS) values are at a level of 5 cm and 10–13 cm, respectively. The mean radial orbit errors of both methods are at a level of 5 cm. For POD repeatability represented as the standard deviations of post-fit residuals and 3D errors by repetitive PODs, however, GRS yields 25% and 13% more robust estimation results than MRRS for post-fit residual and 3D error, respectively. This study shows that PF batch filter with GRS approach using genetic operations is superior to PF batch filter with MRRS in terms of robustness in POD with SLR observations.     

Accuracy assessment of TDRSS-based TOPEX/Poseidon orbit determination

Orbit determination results are obtained for the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using a batch-least-squares estimator available in the Goddard Trajector Determination System (GTDS) to process Tracking and Data Relay Satellite (TDRS) System (TDRSS) measurements. The GTDS orbit solutions are compared with the definitive Precision Orbit Determination (POD) orbit solutions. The root-mean-square (RMS) solution difference in the radial component is 28 centimeters.     

基于机动力模型的GEO卫星恢复期间定轨

北斗卫星导航系统(BDS)中GEO卫星频繁的轨道机动对高精度、实时不间断的导 航服务需求提出了更高要求, 如何在短弧跟踪条件下提高GEO卫星轨道快速 恢复能力, 是提升导航系统服务精度的关键因素. 针对该问题, 本文提出了基 于机动力模型的动力学定轨方法, 尝试利用高精度的C波段转发式测距数据, 辅 以机动期间的遥测遥控信息建立机动力模型, 联合轨控前后的观测数据进行动 力学长弧定轨. 利用BDS中GEO卫星实测数据进行了定轨试验与分析, 结果表明, 恢复期间需要采用解算机动推力的定轨方法, 联合机动前、机动期间和机 动后4h数据定轨的轨道位置精度在20m量级, 径向精度优于2.5m. 该方 法克服了短弧跟踪条件下动力学法定轨和单点定位中的诸多问题, 提供了解决 GEO卫星机动后轨道快速恢复问题的技术方法.