两种观测技术综合精密定轨的探讨

利用T／P卫星的SLR和DORIS实测资料,对两种观测技术综合精密定轨中的加权及其对定轨影响的问题作了初步的探讨。根据所提出的一种经验性的加权方法进行了综合定轨计算,结果表明：对于两种不同技术的观测,相对权选取的恰当与否将影响综合定轨的精度;综合定轨的最优加权不仅依赖于观测资料的精度,还与观测资料的多少和几何分布有关;通过选用最优加权,可使得综合轨的精度优于仅用其中一种技术的定轨精度,综合定轨能有效地提高定轨精度。     

轨道改进中计算状态转移矩阵的分析方法

对当今人卫轨道改进问题,由于力学模型的复杂,精密星历和状态转移矩阵的计算均采用数值方法,这就需要积分两组常微分方程．本文针对状态转移矩阵在定轨中的作用,对定轨弧段不太长的情况,给出了状态转移矩阵的一种分析算法,从而避免数值求解两组常微分方程的问题,并以实际算例证实了这种算法的有效性     

采用重置参数的轨道改进算法

当使用精度差的初始根数作定轨计算时,被估值的模型参数会吸收初值中所含误差而偏离其合理数值(如CD约为2.2),使定轨计算过程的RMS已不再变化,但轨道收敛到与实际状态有偏离的轨道上。文中给出的算例采用重置被歪曲的估值模型参数方法,首先以TLE根数为初值用精密定轨程序解条件方程,然后以第一轮迭代计算结果作为初始根数并重置模型参数,再进行第二轮迭代计算,使定轨计算结果收敛到正确轨道上,文中还使用另一颗激光卫星的双行根数作初值验证了该方法的有效性。较好地解决了因初值不准所引起的定轨计算不收敛,或收敛到与实际状态有偏离的轨道上的问题。最终得出的RMS达到厘米级精度。文中图示了两次定轨计算的RMS变化曲线图、残差分布图,迭代过程的资料采用率及定轨计算结果。     

北斗卫星单系统精密定轨方法对比分析

提出了联合使用载波相位和相位平滑伪距实现北斗卫星双差动力法精密定轨,给出了北斗卫星非差动力法和双差动力法精密定轨的数据处理流程,分析了两种方法的异同.结合实测数据,对比了两种方法的实际定轨效果,结果表明:一定测站布局下,利用两种方法,GEO(Geostationary Earth Orbit Satellite)卫星3维精密定轨精度均能达到1 m左右量级,IGSO(Inclined Geosynchronous Earth Orbit Satellite)和MEO(Medium Earth Orbit Satellite)卫星优于0.5 m,3类卫星的径向定轨精度均优于10 cm.较之非差动力法,双差动力法对GEO卫星精密定轨精度具有一定的改善作用,两者在IGSO卫星精密定轨上效果基本相当,但在MEO卫星定轨上,非差动力法结果更优.     

人卫长弧定轨中的摄动计算问题

人卫精度定轨中摄动星历表的计算，常采用数值方法，对于长弧定轨，将会遇到数值稳定问题，有人提出采用Ｅｎｃｋｅ特别摄动法，但要保持数值稳定，关键在于坐标摄动法数值计算中参考轨道的选取，而不是用不用Ｅｎｃｋｅ变换，更有效的方法是一般数值方法中采用稳定化措施。     

运动学统计定轨法应用于CE-1落月轨迹确定

运动学统计定轨方法不同于常规动力学统计定轨方法和几何法,既不需要建立飞行器的动力学模型,又可以利用整个弧段的测量信息进行统计计算,所以采用运动学统计定轨的方法描述飞行轨迹特别适用于嫦娥二期的动力落月过程。对嫦娥一号(CE-1)落月过程,采用运动学统计定轨方法,联合USB和VLBI数据进行轨道计算并确定CE-1的落月轨迹以及撞月点的位置,将所得结果与动力学统计定轨结果进行比较。计算表明,运动学统计定轨方法精度可以满足工程要求,可以为探月工程二期软着陆阶段的导航提供一种较为可行的定轨方法。     

考虑地球扁率摄动影响的初轨计算方法

在二体问题意义下的短弧定轨,Laplace型方法是最主要最典型的一种初轨计算方法。若测角资料达到10^-4-10^-5精度(相当于2″—20″之间),那么要使定轨精度达到与其相应的程度,地球非球形引力位中的扁率项摄动应该考虑,在此前提下,同样可以采用相应的Laplace型定轨方法。即给出这种严格包含扁率摄动的初轨计算方法的原理和具体计算过程以及计算实例,除采用多资料定轨方法外,这种方法也是提高初轨计算精度的一种途径,它同样可用于多资料的情况,这种方法对于大扁率主天体(即中心天体)的卫星定轨将更有实用价值。     

GPS星历对LEO星载GPS精密定轨精度的影响

目前,越来越多的低轨卫星上都搭载了用于精密定轨的星载GPs接收机,星载GPS已成为低轨卫星精密定轨的主要手段之一.星载GPS精密定轨精度依赖于GPS星历及钟差精度.基于SHORDE-Ⅲ非差动力学定轨功能,以2005年8月1日至8月7日一周的GRACE卫星实测数据为例,采用事后精密轨道(igs)、快速轨道(igr)和超快速轨道(igu)三种GPS星历在同等条件下定轨,估计GPS星历精度对低轨卫星定轨精度的影响,实际计算结果表明igs和igr两类GPS星历定轨精度相当,约为9.5 cm,igu星历定轨精度略低于igs和igr星历,约为10.5cm:高频GPS卫星钟差数据对定轨精度会产生1-6cm影响.     

Lageos卫星精密定轨中“类阻力”摄动的计算

在Lageos卫星的精密定轨中引入了“类阻力”摄动加速度，以减小沿轨误差的影响。本文讨论了“类阻力”摄动的阻力系数CD的不同计算模式对定轨计算的影响。通过实例计算和分析比较，表明了类阻力系数CD在各个短弧上拟合的值变化较大。在一个月的长弧段内反CD作为常系数拟合求解仅是一种求平均值的近拟方法。类阻力系数CD应采用既包含长期变化又包含若干周期变化的计算模式比较合理，具体计算时该取多少项，应根据选用弧段的长短及精度要求来确定。     

星载GPS精密测轨研究及应用

星载ＧＰＳ定轨系统由于其全天侯、价格低、不受卫星高度的影响可达到米级、分米级乃至厘米级的测轨精度等特点已成为低轨道卫星精密定轨的要求,概述了星载ＧＰＳ系统的组成和定轨原理,给出了星载ＧＰＳ测轨的几种主要方法和数学模型,同时根据ＴＯＰＥＸ卫星星载ＧＰＳ实测数据分析了各种定轨方法的测轨精度以及影响定轨精度的各种因素。     

Homogeneous and Heterogeneous Observational Data Weighting Method of Combined Orbit Determination of Near-earth Satellites of Bi-satellite Positioning System and its Applications

In the light of the problem of amalgamation and processing of multisource observational data in the combined orbit determination of near-earth satellites of the bi-satellite positioning system, the optimal weighting method of the improved variance component estimation of the two-step systematic error correction of homogeneous observational data is proposed. Analyses show that the multi-source amalgamation measurement model of the heterogeneous observational data essentially is a multi-structure, multi-parameter non-linear regression model, and the optimal weighting method of the combination of model structure characteristic analysis and variance component estimation of the heterogeneous observational data is established. The realization algorithms of the optimal weighting and the combined orbit determination parameter estimation of the two sorts of observational data are designed, and the simulation experiments of the combined orbit determination are carried out by taking the distances among the two satellites and the backup satellite and the homogeneous observational data and the distance between the two satellites and the heterogeneous observational data of satellite sensor angle measurements as the examples. The results of theoretical analysis and simulation calculation show that for the combined orbit determination of homogeneous observational data, the accuracy of orbit determination obtained by adopting the variance component estimation method of the two-step systematic error correction can be more superior than that obtained by means of the traditional empirical weighting method. For the combined orbit determination of heterogeneous observational data, through the introduction of the weighting factor by which the model structure is characterized the accuracies of the combined orbit determination of the near-earth satellite and geostationary satellite are both improved to a certain extent in comparison with the mean weighting mode.     

基于CEI测量的GEO目标快速轨道恢复

连线干涉测量(Connected Element Interferometry, CEI)是一种全天时全天候的被动测角技术, 已用于空间目标的跟踪监视. 地球静止轨道(Geostationary Earth Orbit, GEO)卫星需要频繁机动以保持轨位或完成其他任务, 其机动后的快速轨道恢复能力对于监视预警极为重要. 针对基于CEI的GEO短弧定轨和预报, 分析了定轨算法的形亏和数亏, 在附加先验轨道约束的短弧定轨基础上, 提出了轨道半长轴初值的自适应优化方法. 利用亚太七号卫星的CEI仿真和实测数据进行了短弧定轨和预报, 实验结果表明, 采用优化后的半长轴初值, 30min短弧定轨和10min预报的卫星位置分量精度均优于4km, 能够满足非合作GEO目标机动后快速轨道恢复的需求.     

月球探测器过渡轨道的短弧定轨方法

论述的短弧定轨,是指在无先验信息情况下又避开多变元迭代的初轨计算方法,它需要相应的动力学问题有一能反映短弧内达到一定精度的近似分析解．探测器进入月球引力作用范围后接近月球时可以处理成相对月球的受摄二体问题,而在地球附近,则可处理成相对地球的受摄二体问题,但在整个过渡段的力模型只能处理成一个受摄的限制性三体问题．而限制性三体问题无分析解,即使在月球引力作用范围外,对于大推力脉冲式的过渡方式,相对地球的变化椭圆轨道的偏心率很大(超过Laplace极限),在考虑月球引力摄动时亦无法构造摄动分析解．就此问题,考虑在地球非球形引力(只包含J2项)和月球引力共同作用下,构造了探测器飞抵月球过渡轨道段的时间幂级数解,在此基础上给出一种受摄二体问题意义下的初轨计算方法,经数值验证,定轨方法有效,可供地面测控系统参考．     

TLE-Aided Orbit Determination Using Single-Station SLR Data

It is difficult to use the single-station satellite laser ranging (SLR) data for orbit determination, due to the singular geometrical distribution of the observations. The single-station data produced by performing the diffuse- reflection SLR on the earth-orbiting space debris are therefore ineffective for orbit improvement. To solve this problem, we propose an orbit determination method by using single-station SLR data in aid of the two-line element set (TLE). For verifying its feasibility, this method is implemented and applied to the orbit determination of the satellite Ajisai, using the single-station SLR data of five passes in one day and the corresponding TLE. And on this basis, the five-day orbit prediction is generated, the result indicates that the errors of predicted positions are less than 40 m. In addition, the potential application of this method in the orbit improvement of space debris is discussed.     

Experiments in orbit determination using numerical methods

The applicability of advanced numerical methods for the solution of the orbit determination problem is studied.The dynamics of the observed object is written as a system of integral equations. This system is solved numerically by representing the components of the force function as linear combinations of B-splines and by applying the multigrid technique. In an outer loop the orbit determination problem is solved using Newton's method.The method is suitable for both preliminary orbit determination and orbit improvement.Most of the programming and part of the analysis was done by G. Zwier, in partial fulfilment of the requirements for obtaining the degree of master in applied mathematics at the Twente University of Technology.     

The Effect of GPS Ephemeris on the Accuracy of Precise Orbit Determination of LEO Satellite-borne GPS

The satellite-borne GPS receivers dedicated to precise orbit determination are now being carried by more and more low earth orbit (LEO) satellites and the satellite-borne GPS has become one of the main means for the precise orbit determination of low earth orbit satellites. The accuracy of satellite-borne GPS precise orbit determination depends on the accuracies of the GPS ephemeris and the clock error. Based on the orbit determination function of SHORDEIII zero-difference dynamics and using the observational data obtained by the GRACE satellites for the week from 2005 August 1 to 7 as an example, three versions of GPS ephemerides (igs, igr and igu) are used to carry out orbit determination under the same conditions and to estimate the effect of the GPS ephemeris accuracy on the accuracy of orbit determination of low earth orbit satellites. Our calculated results show that the two ephemerides, igs and igr, are equivalent to each other in orbit determination accuracy (about 9.5 cm), while igu is slightly less accurate, at about 10.5 cm. The effect produced by the data of the high frequency GPS satellite clock error on the accuracy of orbit determination is 1–6 cm.     

A new method for the determination of an intermediate orbit using three positions of a small body on the celestial sphere

We propose a new method for the determination of the preliminary orbit of a small celestial body using three pairs of its angular coordinates in three moments of time. The method is based on the use of the intermediate orbit we constructed earlier using three position vectors and the corresponding time moments. This intermediate orbit accounts for the main part of the perturbations of the motion of the body under study. We compare the results obtained by the classical Lagrange-Gauss method, Herrick-Gibbs method, generalized Herrick-Gibbs method, and the new method by the examples of the determination of the orbit of the small planet 1566 Icarus. The comparison showed that the new method is a highly efficient tool for the study of perturbed motion. It is especially efficient when applied to high-precision observational data covering short arcs of the orbit.     

A method of astronomical autonomous orbit and attitude determinations for satellites

A method of realtime autonomous orbit determination for earth satellites using the extended Kalman filtering is proposed. The observed quantities are: the satellite-sun direction vector measured by a sun sensor, the satellite-earth and satellite-moon direction vectors measured by an ultraviolet sensor, and the geocentric distance measured by a radar altimeter. At the same time the satellite attitude to the earth is also determined. Results of our simulation of the autonomous orbit determination show that the precision of the orbit determinations is better than 200 m. The effects of the sampling period, orbital inclination, orbital eccentricity and orbital altitude on the precision of orbit determination are analyzed and compared, and certain principles helpful for improving the precision of orbit determination are suggested.