CN102878997A - 一种大偏心率轨道的星上快速高精度外推方法 - Google Patents
一种大偏心率轨道的星上快速高精度外推方法 Download PDFInfo
- Publication number
- CN102878997A CN102878997A CN2012104091566A CN201210409156A CN102878997A CN 102878997 A CN102878997 A CN 102878997A CN 2012104091566 A CN2012104091566 A CN 2012104091566A CN 201210409156 A CN201210409156 A CN 201210409156A CN 102878997 A CN102878997 A CN 102878997A
- Authority
- CN
- China
- Prior art keywords
- tau
- track
- chebyshev polynomials
- satellite
- orbit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种大偏心率轨道的星上快速高精度外推方法,火星探测器任务轨道为大偏心率(e>0.6)椭圆轨道,采用数值积分方法,要做到较高精度的轨道计算,需要考虑较高阶次的火星形状摄动以及其他各种摄动因素影响,从而导致计算量较大,若星上计算机复位或切机,数值积分方法将会中断;解析法只适用于偏心率较小(e<0.6)的椭圆轨道;利用精密轨道提供的探测器位置拟合切比雪夫多项式系数,将会带来较大的拟合误差。本发明通过引入精密轨道与二体轨道的位置差,利用位置差拟合切比雪夫多项式系数,可以很好地解决大偏心率椭圆轨道的星上轨道计算问题,计算量小,且精度高。
Description
技术领域
本发明涉及一种大偏心率轨道的星上方法,属于卫星轨道技术领域。
背景技术
航天器在空间的运动轨迹称为航天器的轨道,通过高精度的外推计算可以得到航天器在过去、当前和未来一段时间内任一时刻的运动状态。火星探测器处于星际飞行,飞行距离最大达到4亿公里,相比地球卫星和月球探测器,火星探测器面临星地距离远、时延大和长时间日凌等问题,火星探测器对GNC的自主性提出了更高要求。因此,星上高精度轨道计算和轨道外推非常重要。求解卫星轨道动力学方程主要有数值法和解析法两种方法。数值方法解算精度高,但存在计算量大以及星上计算机复位或切机时数值积分将会中断的弊端;由于火星探测器任务轨道为大偏心率(e>0.6)椭圆轨道,解析方法也不适用。目前有文献利用精密轨道提供的探测器位置拟合切比雪夫多项式系数,进而实现圆轨道的星上轨道计算,但该方法用于大偏心率椭圆轨道时将会带来较大的拟合误差。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提出一种大偏心率轨道的星上快速高精度外推方法,解决原来星上轨道外推算法拟合误差大的问题,从而提高星上轨道外推计算的精度。
本发明的技术解决方案是:一种大偏心率轨道的星上快速高精度外推方法,所述大偏心率指e>0.6的椭圆轨道,包括以下步骤:
(1)在地面求解大偏心率椭圆轨道的精密轨道与二体轨道的X,Y,Z三轴位置差,之后利用切比雪夫多项式拟合法拟合得到切比雪夫多项式系数;所述精密轨道是指卫星实际运行的轨道,由地面测轨提供,二体轨道是指卫星只受目标天体中心引力作用下的理想轨道,由理论计算得到;所述拟合过程求解如下:
(1.1)首先将时间t∈[t0,t0+Δt]变换成τ∈[-1,1],变换公式为 t∈[t0,t0+Δt] (1)
其中t0为开始历元,Δt为拟合时间区间的长度,τ为标准化的时间变量;
(1.2)将大偏心率椭圆轨道的精密轨道与二体轨道的三轴位置差的X,Y,Z分量采用如下切比雪夫多项式表示
其中公式(4)中n≥2;
其中公式(5)中n≥3;
(1.3)将卫星精密轨道和二体轨道的X轴位置差为观测值,则误差方程为
误差方程的矩阵展开式为
令
则有向量表达式
Vx=BM-fx (8)
为了求得唯一解,采用最小二乘法,使所求的切比雪夫多项式系数阵M满足Vx TVx=min,其中Vx T为Vx的转置矩阵,求解得到
M=(BTB)-1BTfx (9)
(3)卫星利用轨道初值计算二体轨道的三轴位置和三轴速度,利用切比雪夫多项式系数计算三轴位置差和三轴速度差,将其与对应时刻的二体轨道位置和速度相加,得到卫星的实际位置和速度。
本发明与现有技术相比的优点在于:
(1)本发明引入精密轨道与二体轨道的位置差,利用位置差拟合切比雪夫多项式系数,可以很好地解决大偏心率椭圆轨道的星上轨道计算问题,计算量小,且精度高;
(2)本发明本用于星上快速轨道计算,适用于大偏心率的轨道,算法简单,易操作,精度高,具有很强的实用性。
附图说明
图1为本发明的实现流程图;
图2为本发明实施例中利用精密轨道提供的探测器位置拟合切比雪夫多项式系数,进行环火大椭圆轨道的拟合结果;左图代表拟合的位置误差,右图代表拟合的速度误差;
图3为本发明实施例中利用精密轨道与二体轨道的位置差拟合切比雪夫多项式系数,进行环火大椭圆轨道的拟合结果;左图代表拟合的位置误差,右图代表拟合的速度误差。
具体实施方式
如图1所示,本发明具体实现包括以下步骤:
(1)在地面求解大偏心率椭圆轨道的精密轨道与二体轨道的X,Y,Z三轴位置差,之后利用切比雪夫多项式拟合法拟合得到切比雪夫多项式系数;
切比雪夫多项式拟合就是根据给定的数据拟合出一个函数,使其在给定点的函数值与给定值之间的方差和最小,且该函数是以切比雪夫多项式为基函数构成的函数。
所述精密轨道是指卫星实际运行的轨道,由地面测轨提供,二体轨道是指卫星只受目标天体中心引力作用下的理想轨道,由理论计算得到;二体轨道的初值与精密轨道的初值一致,而平近点角M的取值从精密轨道初始平近点角M0开始,以的步长增加,其中μm为火星的引力常数,a为精密轨道的半长轴,n0为轨道角速度,h为轨道角速度的整数倍。
所述拟合过程求解如下:
(1.1)首先将时间t∈[t0,t0+Δt]变换成τ∈[-1,1],变换公式为
其中t0为开始历元,Δt为拟合时间区间的长度,τ为标准化的时间变量;
(1.2)将大偏心率椭圆轨道的精密轨道与二体轨道的三轴位置差的X,Y,
Z分量采用如下切比雪夫多项式表示
其中n≥2;
其中n≥3;
(1.3)将卫星精密轨道和二体轨道的X轴位置差为观测值,则误差方程为
其中k=1,2,…,m为拟合时间区间内的观测值个数,Xk为对应k时刻的X轴位置差,Vx为相应的误差向量;
误差方程的矩阵展开式为
令
则有向量表达式
Vx=BM-fx
为了求得唯一解,采用最小二乘法,使所求的切比雪夫多项式系数阵M满足Vx TVx=min,其中Vx T为Vx的转置矩阵,求得法方程为
(BTB)M-BT fx=0
解得
M=(BTB)-1BT fx
(3)卫星利用轨道初值计算二体轨道的三轴位置和三轴速度,利用切比雪夫多项式系数计算三轴位置差和三轴速度差,将其与对应时刻的二体轨道位置和速度相加,得到卫星的实际位置和速度。
为了更清楚的表明本发明的优点,在此进行数学仿真,仿真条件:火星探测器半长轴9296.7公里,倾角87.6度,偏心率0.607403,升交点赤经73.5度,近地点幅角321.0度,真近点角0.0度。
图2为利用精密轨道提供的探测器位置拟合切比雪夫多项式系数,进行环火大椭圆轨道的拟合结果;图3为利用精密轨道与二体轨道的位置差拟合切比雪夫多项式系数,进行环火大椭圆轨道的拟合结果。可以看出,采用本发明方法后拟合精度得到显著提高。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (1)
1.一种大偏心率轨道的星上快速高精度外推方法,所述大偏心率指e>0.6的椭圆轨道,其特征在于:它包括以下步骤:
(1)在地面求解大偏心率椭圆轨道的精密轨道与二体轨道的X,Y,Z三轴位置差,之后利用切比雪夫多项式拟合法拟合得到切比雪夫多项式系数;所述精密轨道是指卫星实际运行的轨道,由地面测轨提供,二体轨道是指卫星只受目标天体中心引力作用下的理想轨道,由理论计算得到;所述拟合过程求解如下:
(1.1)首先将时间t∈[t0,t0+Δt]变换成τ∈[-1,1],变换公式为
其中t0为开始历元,Δt为拟合时间区间的长度,τ为标准化的时间变量;
(1.2)将大偏心率椭圆轨道的精密轨道与二体轨道的三轴位置差的X,Y,Z分量用切比雪夫多项式表示为
其中公式(4)中n≥2;
其中公式(5)中n≥3;
(1.3)将卫星精密轨道和二体轨道的X轴位置差为观测值,则误差方程为
其中k=1,2,…,m为拟合时间区间内的观测值个数,Xk为对应k时刻的X轴位置差,Vx为相应的误差向量;
误差方程的矩阵展开式为
令
则有向量表达式
VX=BM-fX (8)
为了求得唯一解,采用最小二乘法,使所求的切比雪夫多项式系数阵M满足VX TVX=min,其中Vx T为Vx的转置矩阵,求解得到
M=(BTB)-1BTfX (9)
(3)卫星利用轨道初值计算二体轨道的三轴位置和三轴速度,利用切比雪夫多项式系数计算三轴位置差和三轴速度差,将其与对应时刻的二体轨道位置和速度相加,得到卫星的实际位置和速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409156.6A CN102878997B (zh) | 2012-10-24 | 2012-10-24 | 一种大偏心率轨道的星上快速高精度外推方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409156.6A CN102878997B (zh) | 2012-10-24 | 2012-10-24 | 一种大偏心率轨道的星上快速高精度外推方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102878997A true CN102878997A (zh) | 2013-01-16 |
CN102878997B CN102878997B (zh) | 2015-11-25 |
Family
ID=47480417
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210409156.6A Active CN102878997B (zh) | 2012-10-24 | 2012-10-24 | 一种大偏心率轨道的星上快速高精度外推方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102878997B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514362A (zh) * | 2013-06-24 | 2014-01-15 | 中国电子科技集团公司第二十八研究所 | 基于模型误差补偿的两行根数生成方法 |
CN103744089A (zh) * | 2013-11-15 | 2014-04-23 | 福建纳威导航科技有限责任公司 | 非保守力作用下的导航卫星轨道确定改进的无奇点技术 |
CN103777208A (zh) * | 2013-11-06 | 2014-05-07 | 福建纳威导航科技有限责任公司 | 有势位力作用下的导航卫星轨道确定改进的无奇点技术 |
CN106092099A (zh) * | 2016-06-02 | 2016-11-09 | 哈尔滨工业大学 | 航天器相对位置增量定轨方法 |
CN110816896A (zh) * | 2019-10-28 | 2020-02-21 | 中国空间技术研究院 | 一种卫星星上简易轨道外推方法 |
CN111209523A (zh) * | 2020-01-06 | 2020-05-29 | 中国科学院紫金山天文台 | 适用于大偏心率轨道密集星历精密计算的快速处理方法 |
CN113641949A (zh) * | 2021-08-05 | 2021-11-12 | 中国西安卫星测控中心 | 一种地球同步转移段轨道根数高精度拟合方法 |
CN114383619A (zh) * | 2021-12-07 | 2022-04-22 | 上海航天控制技术研究所 | 一种高精度轨道计算方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102116630A (zh) * | 2009-12-31 | 2011-07-06 | 北京控制工程研究所 | 一种绕火星探测器的星上快速高精度确定方法 |
-
2012
- 2012-10-24 CN CN201210409156.6A patent/CN102878997B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102116630A (zh) * | 2009-12-31 | 2011-07-06 | 北京控制工程研究所 | 一种绕火星探测器的星上快速高精度确定方法 |
Non-Patent Citations (1)
Title |
---|
俞朔春等: "《基于GPS的静止轨道卫星自主定轨技术研究》", 《航天控制》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514362A (zh) * | 2013-06-24 | 2014-01-15 | 中国电子科技集团公司第二十八研究所 | 基于模型误差补偿的两行根数生成方法 |
CN103514362B (zh) * | 2013-06-24 | 2016-11-09 | 中国电子科技集团公司第二十八研究所 | 基于模型误差补偿的两行根数生成方法 |
CN103777208A (zh) * | 2013-11-06 | 2014-05-07 | 福建纳威导航科技有限责任公司 | 有势位力作用下的导航卫星轨道确定改进的无奇点技术 |
CN103744089A (zh) * | 2013-11-15 | 2014-04-23 | 福建纳威导航科技有限责任公司 | 非保守力作用下的导航卫星轨道确定改进的无奇点技术 |
CN106092099A (zh) * | 2016-06-02 | 2016-11-09 | 哈尔滨工业大学 | 航天器相对位置增量定轨方法 |
CN106092099B (zh) * | 2016-06-02 | 2018-10-02 | 哈尔滨工业大学 | 航天器相对位置增量定轨方法 |
CN110816896A (zh) * | 2019-10-28 | 2020-02-21 | 中国空间技术研究院 | 一种卫星星上简易轨道外推方法 |
CN111209523A (zh) * | 2020-01-06 | 2020-05-29 | 中国科学院紫金山天文台 | 适用于大偏心率轨道密集星历精密计算的快速处理方法 |
CN111209523B (zh) * | 2020-01-06 | 2020-12-29 | 中国科学院紫金山天文台 | 适用于大偏心率轨道密集星历精密计算的快速处理方法 |
WO2021139129A1 (zh) * | 2020-01-06 | 2021-07-15 | 中国科学院紫金山天文台 | 适用于大偏心率轨道密集星历精密计算的快速处理方法 |
US11319094B2 (en) * | 2020-01-06 | 2022-05-03 | Purple Mountain Observatory, Chinese Academy Of Sciences | Method for accurately and efficiently calculating dense ephemeris of high-eccentricity orbit |
CN113641949A (zh) * | 2021-08-05 | 2021-11-12 | 中国西安卫星测控中心 | 一种地球同步转移段轨道根数高精度拟合方法 |
CN114383619A (zh) * | 2021-12-07 | 2022-04-22 | 上海航天控制技术研究所 | 一种高精度轨道计算方法 |
CN114383619B (zh) * | 2021-12-07 | 2023-09-05 | 上海航天控制技术研究所 | 一种高精度轨道计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102878997B (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102878997B (zh) | 一种大偏心率轨道的星上快速高精度外推方法 | |
Visser et al. | Energy integral method for gravity field determination from satellite orbit coordinates | |
CN106697333B (zh) | 一种航天器轨道控制策略的鲁棒性分析方法 | |
CN102878995B (zh) | 一种静止轨道卫星自主导航方法 | |
CN110132261B (zh) | 一种基于数值拟合的高精度星上轨道预报方法 | |
CN109344449B (zh) | 航天器月地转移轨道逆向设计方法 | |
CN109032176B (zh) | 一种基于微分代数的地球同步轨道确定和参数确定方法 | |
CN104076819B (zh) | 一种圆参考轨道下卫星有界伴飞的边界控制方法 | |
Gerlach et al. | CHAMP gravity field recovery using the energy balance approach | |
CN101226062B (zh) | 一种星上实时计算环月轨道的方法 | |
CN102495950B (zh) | 一种适用于太阳同步轨道的倾角偏置量获取方法 | |
CN104848862A (zh) | 一种环火探测器精密同步定位守时方法及系统 | |
Abdelrahman et al. | Sigma-point Kalman filtering for spacecraft attitude and rate estimation using magnetometer measurements | |
CN104501804A (zh) | 一种基于gps测量数据的卫星在轨轨道预报方法 | |
CN107144283A (zh) | 一种用于深空探测器的高可观度光学脉冲星混合导航方法 | |
CN103645489A (zh) | 一种航天器gnss单天线定姿方法 | |
McLaughlin et al. | Drag coefficient estimation in orbit determination | |
De Marchi et al. | Constraining the Nordtvedt parameter with the BepiColombo Radioscience experiment | |
Somodi et al. | Application of numerical integration techniques for orbit determination of state-of-the-art LEO satellites | |
CN108959665B (zh) | 适用于低轨卫星的轨道预报误差经验模型生成方法及系统 | |
Bakhtiari et al. | Long-term effects of main-body’s obliquity on satellite formation perturbed by third-body gravity in elliptical and inclined orbit | |
Traub et al. | Assessment of the dependencies of realistic differential drag controlled in-plane reconfiguration maneuvers on relevant parameters | |
CN104932520A (zh) | 一种存在测量误差下的航天器雨滴形状悬停构型控制方法 | |
Zhang et al. | Robustness analysis method for orbit control | |
CN112327333A (zh) | 一种卫星位置计算方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |