SGP4/SDP4模型精度分析

本文基于最新发布的SGP4/SDP4(Simplified General Perturbation Version 4/Simplified Deep-space Perturbation Version 4)模型设计了一套定轨方案,从空间目标库中挑选出不同类型和轨道参数的1120个目标进行计算,定量给出了SGP4/SDP4模型处理不同类型空间目标的定轨预报精度.结果表明:近地目标定轨精度为百米量级;半同步和同步轨道定轨精度平均为0.7和1.9km.椭圆轨道目标的定轨精度与偏心率有关,除少数e>0.8的椭圆轨道目标,绝大多数椭圆轨道目标定轨误差均小于10km.用SGP4/SDP4模型对近地目标预报3天,半同步轨道预报30天,同步轨道预报15天,椭圆轨道预报1天,预报误差一般不超过40km.     

An Accuracy Analysis of the SGP4/SDP4 Model

Based on the latest release of the SGP4/SDP4 (Simplified General Perturbation Version 4/ Simplified Deep-space Perturbation Version 4) model, in this paper we have designed an orbit determination program. Through calculations for the 1120 objects with various types and orbital elements selected from the space objects database, we have obtained the accuracies of the orbit determination prediction dealt with various types of space objects by the SGP4/SDP4 model. The results show that the accuracies of the near-earth objects are in the order of magnitude of 100 meters; the averages of the orbit determination accuracies of the semi-synchronous and geosynchronous orbits are, respectively, 0.7 and 1.9 km. The orbit determination accuracies of the elliptical orbit objects are related to their eccentricities. Except for few elliptical orbit objects with e > 0.8, the orbit determination errors of the vast majority of the elliptical orbit objects are all less than 10 km. By using the SGP4/SDP4 model to make 3 days predictions for near-earth objects, 30 days for semi-synchronous orbit objects, 15 days for geosynchronous orbit objects and 1 day for elliptical orbit objects, the errors of prediction generally don’t exceed 40 km.     

太阳同步卫星的控制现状及发展趋势

在介绍太阳同步轨道相关知识的基础上,着重分析了该轨道的高度因素,给出了当前在轨太阳同步卫星轨道高度的上下限,介绍了几种典型的低轨道太阳同步卫星,同时给出了太阳同步卫星轨道控制方法的发展状况及趋势,最后综述了当前国内外太阳同步轨道的轨控模式与维持策略。     

同步卫星短弧定轨

同步卫星受到摄动力的影响,它的实际轨道有一点漂移．卫星需要不断的调轨调姿,以保证其正常运行．为了研究卫星在几小时,甚至更短的时间内的轨迹情况,采用短弧段定轨法．用动力学方法进行短弧定轨,分别研究1小时和15分钟定轨并进行比较,目的是为了在同步轨道卫星变轨后,能尽快地为卫星提供精密的预报轨道．此外,在系列短弧定轨后,得到精密轨道系列,为研究轨道变化的力学因素及研究短弧中卫星转发器时延变化规律等提供依据．     

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.     

GNSS精密轨道产品不连续性分析及评估

按照目前的国际规范, 高精度GNSS (Global Navigation Satellite System)轨道产品一般以天为周期进行发布, 提供给用户使用. 连续使用多天的产品存在不同天间的跳变问题. 利用德国地学研究中心(GFZ)、欧洲定轨中心(COD)、欧空局(ESA)、美国喷气试验室(JPL)以及上海天文台(SHA)共5个GNSS分析中心2013---2017年的轨道产品, 分析了轨道跳变的特性. 计算结果表明: GFZ、COD、ESA、SHA和JPL的3维轨道跳变平均分别为7.79cm、1.51cm、7.77cm、11.75cm和2.51cm. 轨道跳变序列的周期特性分析表明: 序列存在90d、120d、340d左右的显著周期项, 对应于海潮对地球自转的影响, 其振幅为数毫米至1cm左右. 表明精密轨道确定需要进一步精化该项模型; GPS的跳变序列还存在与卫星星座相关的175d和352d左右的交点年显著周期项. 此外, 针对COD产品外推轨道的分析, 验证了地球反照辐射压和太阳光压模型等动力学模型对轨道的差异.     

Discovery of Earth's quasi‐satellite

Abstract— The newly discovered asteroid 2003 YN₁₀₇ is currently a quasi‐satellite of the Earth, making a satellite‐like orbit of high inclination with apparent period of one year. The term quasi‐satellite is used since these large orbits are not completely closed, but rather perturbed portions of the asteroid's orbit around the Sun. Due to its extremely Earth‐like orbit, this asteroid is influenced by Earth's gravity to remain within 0.1 AU of the Earth for approximately 10 years (1997 to 2006). Prior to this, it had been on a horseshoe orbit closely following Earth's orbit for several hundred years. It will re‐enter such an orbit, and make one final libration of 123 years, after which it will have a close interaction with the Earth and transition to a circulating orbit. Chaotic effects limit our ability to determine the origin or fate of this object.     

月球卫星轨道力学综述

月球探测器的运动通常可分为3个阶段，这3个阶段分别对应3种不同类型的轨道：近地停泊轨道、向月飞行的过渡轨道与环月飞行的月球卫星轨道。近地停泊轨道实为一种地球卫星轨道；过渡轨道则涉及不同的过渡方式(大推力或小推力等)；环月飞行的月球卫星轨道则与地球卫星轨道有很多不同之处，它决不是地球卫星轨道的简单克隆。针对这一点，全面阐述月球卫星的轨道力学问题，特别是环月飞行中的一些热点问题，如轨道摄动解的构造、近月点高度的下降及其涉及的卫星轨道寿命、各种特殊卫星(如太阳同步卫星和冻结轨道卫星等)的轨道特征、月球卫星定轨等。     

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

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

The Application of GIM in Precise Orbit Determination for LEO Satellites with Single-Frequency GPS Measurements

With the precise GPS ephemeris and clock error available, the iono- spheric delay is left as the dominant error source in the single-frequency GPS data. Thus, the removal of ionospheric effects is a ma jor prerequisite for an improved orbit reconstruction of LEO satellites based on the single-frequency GPS data. In this paper, the use of Global Ionospheric Maps (GIM) in kine- matic and dynamic orbit determinations for LEO satellites with single-frequency GPS pseudorange measurements is discussed first, and then, estimating the iono- spheric scale factor to remove the ionospheric effects from the C/A-code pseu- dorange measurements for both kinematic and dynamic orbit determinations is addressed. As it is known that the ionospheric delay of space-borne GPS sig- nals is strongly dependent on the orbit altitudes of LEO satellites, we select the real C/A-code pseudorange measurement data of the CHAMP, GRACE, TerraSAR-X and SAC-C satellites with altitudes between 300 km and 800 km as sample data in this paper. It is demonstrated that the approach to eliminating ionospheric effects in C/A-code pseudorange measurements by estimating the ionospheric scale factor is highly effective. Employing this approach, the accu- racy of both kinematic and dynamic orbits can be improved notably. Among those five LEO satellites, CHAMP with the lowest orbit altitude has the most remarkable improvements in orbit accuracy, which are 55.6% and 47.6% for kine- matic and dynamic orbits, respectively. SAC-C with the highest orbit altitude has the least improvements in orbit accuracy accordingly, which are 47.8% and 38.2%, respectively.     

GIM在LEO卫星单频GPS定轨中的应用

电离层延迟误差是单频GPS(Global Positioning System)数据最主要的误差源,为提高基于单频GPS数据的LEO(Low Earth Orbiting)卫星定轨精度,必须消除/减弱GPS观测数据中电离层延迟影响.研究了全球电离层模型GIM(Global IonosphericMaps)在基于单频GPS伪距数据的低轨卫星运动学和动力学定轨中的应用,并通过估算电离层尺度因子的方法消除C/A码伪距观测量中电离层延迟影响.由于LEO卫星星载GPS信号受电离层延迟影响与卫星轨道高度相关,选取了轨道高度在300～800 km的CHAMP(CHAllenging Mini-satellite Payload)、GRACE(Gravity Recovery AndClimate Experiment)、TerraSAR-X及SAC-C等LEO卫星C/A码伪距观测量作为试算数据.CHAMP等卫星实测数据计算结果表明:以JPL(Jet Propulsion Laboratory)发布的GIM模型作为背景模型,通过电离层比例因子法能很好地消除C/A码伪距观测量中电离层延迟影响,提高LEO卫星运动学和动力学定轨精度,其中,CHAMP卫星轨道最低,受电离层延迟影响最严重,定轨精度提高最显著,分别为55.6%和47.6%;SAC-C卫星轨道高度最高,受电离层延迟影响最小,相应的定轨精度提高幅度也最低,分别为47.8%和38.2%.     

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.     

Luni-solar perturbations of an Earth satellite

Luni-solar perturbations of the orbit of an artificial Earth satellite are given by modifying the analytical theory of an artificial lunar satellite derived by the author in recent papers. Expressions for the first-order changes, both secular and periodic, in the elements of the geocentric Keplerian orbit of the earth satellite are given, the moon's geocentric orbit, including solar perturbations in it, being found by using Brown's lunar theory.The effects of Sun and Moon on the satellite orbit are described to a high order of accuracy so that the theory may be used for distant earth satellites.     

Properties of linearized mappings associated with periodic orbits in the three-body problem

In this paper three results on the linearized mapping associated with the plane three body problem near a periodic orbit are established. It is first shown that linear stability of such an orbit is independent of initial position on the orbit and of coordinate system. Second, the relation of Hénon connecting the rates of change of rotation angle and period on an isoenergetic family of periodic orbits is proved, together with a similar relation for families of orbits closing exactly in a rotating coordinate system. Finally, a condition for a critical orbit is given which is applicable to any family of periodic orbits.     

Optimal periodic relative orbit and rectilinear relative orbits with eccentric reference orbits

The problem of two-body linearized periodic relative orbits with eccentric reference orbits is studied in this paper. The periodic relative orbit in the target-orbital coordinate system can be used in fly-around and formation-flying orbit design. Based on the closed-form solutions to the Tschauner–Hempel equations, the initial condition for periodic relative orbits is obtained. Then the minimum-fuel periodic-orbit condition with a single impulse is analytically derived for given initial position and velocity vectors. When considering the initial coasting time, the impulse position of the global minimum-fuel periodic orbit is proved to be near to the perigee of the target and can be obtained by numerical optimization algorithms. Moreover, the condition for a special periodic orbit, i.e., the rectilinear relative orbit in the target-orbital frame, is obtained. Numerical simulations are used to demonstrate the efficacy of the method, and show the geometry of the periodic relative orbit and the rectilinear relative orbit.     

The Data Depth-weight-kernel Estimation of the Model Error of Satellite Orbit Determination

It is an objective fact that there exists error in the satellite dynamic model and it will be transferred to satellite orbit determination algorithm, forming a part of the connotative model error. Mixed with the systematic error and random error of the measurements, they form the unitive model error and badly restrict the precision of the orbit determination. We deduce in detail the equations of orbit improvement for a system with dynamic model error, construct the parametric model for the explicit part of the model and nonparametric model for the error that can not be explicitly described. We also construct the partially linear orbit determination model, estimate and fit the model error using a two-stage estimation and a kernel function estimation, and finally make the corresponding compensation in the orbit determination. Beginning from the data depth theory, a data depth weight kernel estimator for model error is proposed for the sake of promoting the steadiness of model error estimation. Simulation experiments of SBSS are performed. The results show clearly that the model error is one of the most important effects that will influence the precision of the orbit determination. The kernel function method can effectively estimate the model error, with the window width as a major restrict parameter. A data depth-weight-kernel estimation, however, can improve largely the robustness of the kernel function and therefore improve the precision of orbit determination.     

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

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

转发器式卫星测轨方法

提出了转发器式卫星测轨方法。发射信号和接收信号的不同组合,形成不同模式的转发器式卫星测轨方法,并给出了不同模式下归算转发器式卫星测轨的公式。自发自收模式下的转发器式卫星测轨方法的观测和计算结果表明,定轨观测残差小于9cm。用转发器式卫星测轨方法,不但能给出高精度时间比对结果,而且能给出高精度卫星轨道和卫星预报轨道。     

空间碎片轨道预报中动力学模型基准的选取

在不同的轨道预报场景中, 使用的动力学模型也不同. 例如, 在低轨空间碎片的预报中大气阻力是十分重要的摄动力, 而到了中高轨, 大气阻力就可以忽略不计. 如何为不同轨道类型的空间碎片选择最优(满足精度要求下的最简)动力学模型还没有系统、详尽的研究. 将对不同精度需求、不同轨道类型下的大批量轨道进行预报, 通过比较不同动力学模型下的预报结果, 给出各种预报场景的最优动力学模型建议. 可以为不同轨道类型的空间碎片在轨道预报时选择基准动力学模型提供参考或标准.     

从地面发射月球探测器的窗口选择

典型的月球探测器飞行轨道包括地球停泊轨道段、地月转移轨道段、月球卫星轨道段和着月轨道段。首先介绍了设计从地面发射月球探测器轨道典型的约束条件;然后,借助于二体假设,建立解析表达式,分析各种约束对窗口选择的影响,给出了各轨道段概略的飞行时间和粗窗口;最后,利用精确的探测器轨道动力学模型,计算精窗口,并给出了一则算例,所得结论可为月球探测器轨道发射、轨道设计提供依据。