CN110132261B - 一种基于数值拟合的高精度星上轨道预报方法 - Google Patents
一种基于数值拟合的高精度星上轨道预报方法 Download PDFInfo
- Publication number
- CN110132261B CN110132261B CN201811365986.7A CN201811365986A CN110132261B CN 110132261 B CN110132261 B CN 110132261B CN 201811365986 A CN201811365986 A CN 201811365986A CN 110132261 B CN110132261 B CN 110132261B
- Authority
- CN
- China
- Prior art keywords
- orbit
- forecasting
- satellite
- instantaneous
- precision
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/04—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
- G01C21/08—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于数值拟合的高精度星上轨道预报方法,首先获得包括最新太阳辐射流量和地磁指数在内的大气环境参数,卫星星体参数以及卫星最新的精密轨道根数,结合星体和空间环境参数,利用地面高精度数值法进行轨道预报,获得一组预报弧段内卫星瞬时轨道根数,然后转换为平均根数进行数值拟合,得到每个轨道根数的变化项系数,根据星上预报精度要求上注多组初轨进行分段外推,根据上注的初始轨道根数及变化项系数进行轨道预报,得到预报时刻卫星的平均根数,如果预报精度要求较高,注入多组初值进行分阶段预报,将预报时刻卫星的平均根数转换为瞬时根数,从而得到卫星的瞬时位置,本发明解决了现有技术中存在的星上轨道预报精度差的问题。
Description
技术领域
本发明属于航天测量与控制技术领域,具体涉及一种基于数值拟合的高精度星上轨道预报方法。
背景技术
由于受星载计算机等条件的限制,星上轨道预报通常采用简化的分析法模型,预报精度较差,对于低轨卫星,加之大气密度的影响,预报精度更差。
发明内容
本发明的目的是提供一种基于数值拟合的高精度星上轨道预报方法,解决了现有技术中存在的星上轨道预报精度差的问题。
本发明所采用的技术方案是,一种基于数值拟合的高精度星上轨道预报方法,具体按照以下步骤实施:
步骤1、获得包括最新太阳辐射流量和地磁指数在内的大气环境参数,卫星星体参数以及卫星最新的精密轨道根数,结合星体和空间环境参数,利用地面高精度数值法进行轨道预报,获得一组预报弧段内卫星瞬时轨道根数;
步骤2、将步骤1获得的预报弧段内的一组瞬时轨道根数转换为平均根数;
步骤3、用步骤2得到的预报弧段内的一组平均根数进行数值拟合,得到每个轨道根数的变化项系数;
步骤4、根据星上预报精度要求上注多组初轨进行分段外推,根据上注的初始轨道根数及变化项系数用步骤3进行轨道预报,得到预报时刻卫星的平均根数,如果预报精度要求较高,注入多组初值进行分阶段预报;
步骤5、将预报时刻卫星的平均根数转换为瞬时根数,从而得到卫星的瞬时位置。
本发明的特点还在于,
步骤1具体如下:
在地心天球坐标系中,卫星运动方程为:
根据公式(1)~公式(3)进行高精度数值法轨道预报,生成一组卫星的瞬时轨道根数。
其中,σ(t)为卫星t时刻瞬时根数,
σ=(a,e,i,Ω,ω,M) (5)
其中,a轨道半长径,Ω升交点经度,e轨道偏心率,ω近地点幅角,i轨道倾角,M为平近点角,
步骤3具体如下:
具体拟合方法为:设拟合数据序列包含k个数据,其中轨道半长轴分别为a1,a2,...,ak,则平均半长轴为:
轨道角速度n为:
对任意时刻的每个轨道根数,ti=[t0,tf],存在
f(ti)=ΦT(ti)P (8)
其中,i=1,2,…n,
系数向量P满足的线性方程组为:
ΦP=F (9)
P=[Pk0 Pk1 Pk1s Pk1c Pk2s Pk2c Pk3s Pk3c]T (12)
对上式中系数矩阵P进行估计求解,得到卫星每个根数的变化项系数。
步骤4中上注初轨数量为3组/天,即外推8小时换一组初轨,具体预报方法如下:
根据上式利用注入的初值和变化项系数进行轨道预报。
步骤5具体如下:
最后根据上式将预报结果转换为瞬时根数,即完成高精度星上轨道预报。
本发明的有益效果是,一种基于数值拟合的高精度星上轨道预报方法,先采用高精度数值法进行轨道预报,并将得到的一组瞬时根数转换为平均根数,基于这组平均根数进行多项式拟合,给出每个轨道根数的一阶和二阶变化项系数,再根据星上轨道预报精度要求,分几组初轨和拟合系数进行上注,最后根据预报时间判断用哪组初轨进行分段预报,从而达到提高星上轨道预报精度的目的。
附图说明
图1是本发明一种基于数值拟合的高精度星上轨道预报方法中外推24 小时的误差分析情况图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明所采用的技术方案是,一种基于数值拟合的高精度星上轨道预报方法,具体按照以下步骤实施:
步骤1、获得包括最新太阳辐射流量和地磁指数在内的大气环境参数,卫星星体参数以及卫星最新的精密轨道根数,结合星体和空间环境参数,利用地面高精度数值法进行轨道预报,获得一组预报弧段内卫星瞬时轨道根数, 具体如下:
在地心天球坐标系中,卫星运动方程为:
根据公式(1)~公式(3)进行高精度数值法轨道预报,生成一组卫星的瞬时轨道根数;
其中,σ(t)为卫星t时刻瞬时根数,
σ=(a,e,i,Ω,ω,M) (5)
其中,a轨道半长径,Ω升交点经度,e轨道偏心率,ω近地点幅角,i轨道倾角,M为平近点角,
步骤3、用步骤2得到的预报弧段内的一组平均根数进行数值拟合,得到每个轨道根数的变化项系数,具体如下:
具体拟合方法为:设拟合数据序列包含k个数据,其中轨道半长轴分别为a1,a2,...,ak,则平均半长轴为:
轨道角速度n为:
对任意时刻的每个轨道根数,ti=[t0,tf],存在
f(ti)=ΦT(ti)P (8)
其中,i=1,2,…n,
系数向量P满足的线性方程组为:
ΦP=F (9)
P=[Pk0 Pk1 Pk1s Pk1c Pk2s Pk2c Pk3s Pk3c]T (12)
对上式中系数矩阵P进行估计求解,得到卫星每个根数的变化项系数;
步骤4、根据星上预报精度要求上注多组初轨进行分段外推,根据上注的初始轨道根数及变化项系数用步骤3进行轨道预报,得到预报时刻卫星的平均根数,如果预报精度要求较高,注入多组初值进行分阶段预报,其中上注初轨数量为3组/天,即外推8小时换一组初轨,具体预报方法如下:
根据上式利用注入的初值和变化项系数进行轨道预报;
步骤5、将预报时刻卫星的平均根数转换为瞬时根数,从而得到卫星的瞬时位置,具体如下:
最后根据上式将预报结果转换为瞬时根数,即完成高精度星上轨道预报。
实施例
采取如表1所示的初始轨道数据:
表1:初始轨道数据
用高精度轨道预报软件预报两天的弹道,并转换为平根数,采用二次多项式函数拟合后得到每个轨道根数的变化项的系数,得到预报公式为:
选择弹道第一点为平根数注入需求点,相应的变率输出结果如表2所示:
表2:注入需求点与变率输出结果对应表
用拟合后的系数进行外推48小时,将外推结果与标准弹道比较,其位置误差最大为355m,主要是沿迹向,在径向和法线方向的误差较小,约为 20米;若外推24小时,其最大位置误差小于100m;若外推12小时约为20m;外推6小时约为10m。如图1所示,用简化的分析法模型进行预报24小时误差约为500m,由此可见,用该方法可大大提高了星上的轨道预报精度。
本发明一种基于数值拟合的高精度星上轨道预报方法,能够解决由于星上计算机条件限制采用简化分析法轨道预报精度较差的问题,使用地面高精度数值法轨道预报,并用预报弧段内轨道平均根数进行拟合轨道变化项,并按星上轨道预报的精度要求进行多组初轨预报的方法,达到提高星上轨道预报精度的目的。
Claims (3)
1.一种基于数值拟合的高精度星上轨道预报方法,其特征在于,具体按照以下步骤实施:
步骤1、获得包括最新太阳辐射流量和地磁指数在内的大气环境参数,卫星星体参数以及卫星最新的精密轨道根数,结合星体和空间环境参数,利用地面高精度数值法进行轨道预报,获得一组预报弧段内卫星瞬时轨道根数;
所述步骤1具体如下:
在地心天球坐标系中,卫星运动方程为:
根据公式(1)~公式(3)进行高精度数值法轨道预报,生成一组卫星的瞬时轨道根数;
步骤2、将步骤1获得的预报弧段内的一组瞬时轨道根数转换为平均根数;
其中,σ(t)为卫星t时刻瞬时根数,
σ=(a,e,i,Ω,ω,M) (5)
其中,a轨道半长径,Ω升交点经度,e轨道偏心率,ω近地点幅角,i轨道倾角,M为平近点角,
步骤3、用步骤2得到的预报弧段内的一组平均根数进行数值拟合,得到每个轨道根数的变化项系数;
所述步骤3具体如下:
具体拟合方法为:设拟合数据序列包含k个数据,其中轨道半长轴分别为a1,a2,...,ak,则平均半长轴为:
轨道角速度n为:
对任意时刻的每个轨道根数,ti=[t0,tf],存在
f(ti)=ΦT(ti)P (8)
其中,i=1,2,…n,
系数向量P满足的线性方程组为:
ΦP=F (9)
P=[Pk0 Pk1 Pk1s Pk1c Pk2s Pk2c Pk3s Pk3c]T (12)
对上式中系数矩阵P进行估计求解,得到卫星每个根数的变化项系数;
步骤4、根据星上预报精度要求上注多组初轨进行分段外推,根据上注的初始轨道根数及变化项系数用步骤3进行轨道预报,得到预报时刻卫星的平均根数,如果预报精度要求较高,注入多组初值进行分阶段预报;
步骤5、将预报时刻卫星的平均根数转换为瞬时根数,从而得到卫星的瞬时位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811365986.7A CN110132261B (zh) | 2018-11-16 | 2018-11-16 | 一种基于数值拟合的高精度星上轨道预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811365986.7A CN110132261B (zh) | 2018-11-16 | 2018-11-16 | 一种基于数值拟合的高精度星上轨道预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110132261A CN110132261A (zh) | 2019-08-16 |
CN110132261B true CN110132261B (zh) | 2023-03-14 |
Family
ID=67568181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811365986.7A Active CN110132261B (zh) | 2018-11-16 | 2018-11-16 | 一种基于数值拟合的高精度星上轨道预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110132261B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110816896B (zh) * | 2019-10-28 | 2021-06-11 | 中国空间技术研究院 | 一种卫星星上简易轨道外推方法 |
CN111060934B (zh) * | 2019-12-25 | 2023-08-04 | 中国西安卫星测控中心 | 低轨航天器短弧定轨瞬时根数精度评估方法 |
CN111854765B (zh) * | 2020-06-08 | 2022-04-26 | 中国人民解放军战略支援部队航天工程大学 | 一种中轨道导航卫星轨道长期预报方法 |
CN111953401B (zh) * | 2020-07-28 | 2022-06-07 | 中国西安卫星测控中心 | 一种微小卫星自主请求式轨道服务系统 |
CN112949067B (zh) * | 2021-03-05 | 2022-11-25 | 航天科工空间工程发展有限公司 | 一种星载空间目标轨道平滑方法 |
CN113297745B (zh) * | 2021-05-28 | 2022-09-02 | 中国人民解放军63921部队 | 一种基于短弧拟合位置的双弧段轨道改进方法 |
CN113447957B (zh) * | 2021-06-23 | 2022-11-25 | 北京航空航天大学 | 一种高精度实时卫星轨道瞬时半长轴的去振方法 |
CN113641949B (zh) * | 2021-08-05 | 2023-03-28 | 中国西安卫星测控中心 | 一种地球同步转移段轨道根数高精度拟合方法 |
CN114383619B (zh) * | 2021-12-07 | 2023-09-05 | 上海航天控制技术研究所 | 一种高精度轨道计算方法 |
CN114282166B (zh) * | 2021-12-22 | 2024-04-19 | 中国西安卫星测控中心 | 基于环境函数矩阵的卫星经度计算方法、变轨判别方法 |
CN115402539B (zh) * | 2022-11-02 | 2023-01-31 | 宁波天巡科技有限公司 | 一种基于轨道和空间环境数据的卫星变轨检测方法及系统 |
CN115828035B (zh) * | 2023-02-16 | 2023-04-21 | 北京控制工程研究所 | 静止轨道卫星太阳光压干扰力矩的在轨拟合估计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002249100A (ja) * | 2001-02-21 | 2002-09-03 | Mitsubishi Electric Corp | 衛星追尾システム及び衛星追尾方法 |
CN102591343A (zh) * | 2012-02-09 | 2012-07-18 | 航天东方红卫星有限公司 | 基于两行根数的卫星轨道维持控制方法 |
WO2012125293A2 (en) * | 2011-03-11 | 2012-09-20 | Sorce4 Llc. | Offline ephemeris prediction |
CN106959457A (zh) * | 2017-04-25 | 2017-07-18 | 清华大学 | 一种用于卫星导航的glonass历书参数估计方法 |
CN107797130A (zh) * | 2017-10-16 | 2018-03-13 | 中国西安卫星测控中心 | 低轨航天器多点多参数轨道上行数据计算方法 |
CN108761507A (zh) * | 2018-05-21 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 基于短弧定轨和预报的导航卫星轨道快速恢复方法 |
-
2018
- 2018-11-16 CN CN201811365986.7A patent/CN110132261B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002249100A (ja) * | 2001-02-21 | 2002-09-03 | Mitsubishi Electric Corp | 衛星追尾システム及び衛星追尾方法 |
WO2012125293A2 (en) * | 2011-03-11 | 2012-09-20 | Sorce4 Llc. | Offline ephemeris prediction |
CN102591343A (zh) * | 2012-02-09 | 2012-07-18 | 航天东方红卫星有限公司 | 基于两行根数的卫星轨道维持控制方法 |
CN106959457A (zh) * | 2017-04-25 | 2017-07-18 | 清华大学 | 一种用于卫星导航的glonass历书参数估计方法 |
CN107797130A (zh) * | 2017-10-16 | 2018-03-13 | 中国西安卫星测控中心 | 低轨航天器多点多参数轨道上行数据计算方法 |
CN108761507A (zh) * | 2018-05-21 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 基于短弧定轨和预报的导航卫星轨道快速恢复方法 |
Non-Patent Citations (2)
Title |
---|
基于GPS测量数据的卫星在轨轨道预报算法研究;刘燎等;《上海航天》;20170430(第02期);全文 * |
基于动力学轨道拟合的LEO卫星轨道预报精度分析;王友存等;《全球定位系统》;20180831(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110132261A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110132261B (zh) | 一种基于数值拟合的高精度星上轨道预报方法 | |
Müller et al. | High performance clocks and gravity field determination | |
CN107797130B (zh) | 低轨航天器多点多参数轨道上行数据计算方法 | |
CN107329146B (zh) | 一种导航卫星低轨监测星座的优化设计方法 | |
CN112257343B (zh) | 一种高精度地面轨迹重复轨道优化方法及系统 | |
CN102591343B (zh) | 基于两行根数的卫星轨道维持控制方法 | |
Canuto et al. | Spacecraft dynamics and control: the embedded model control approach | |
CN100501331C (zh) | 基于x射线脉冲星的导航卫星自主导航系统与方法 | |
CN102495950B (zh) | 一种适用于太阳同步轨道的倾角偏置量获取方法 | |
CN110068845B (zh) | 一种基于平根数理论确定卫星理论轨道的方法 | |
CN110609972B (zh) | 一种指定发射仰角的自由弹道构造方法 | |
CN109323698B (zh) | 一种空间目标陨落多模型跟踪引导方法 | |
CN110816896B (zh) | 一种卫星星上简易轨道外推方法 | |
Verma et al. | Electron density distribution and solar plasma correction of radio signals using MGS, MEX, and VEX spacecraft navigation data and its application to planetary ephemerides | |
CN102878997A (zh) | 一种大偏心率轨道的星上快速高精度外推方法 | |
CN104751012A (zh) | 沿飞行弹道的扰动引力快速逼近方法 | |
CN111238489B (zh) | 一种低轨卫星大气阻力摄动建模和计算方法 | |
CN104501804A (zh) | 一种基于gps测量数据的卫星在轨轨道预报方法 | |
CN112769466B (zh) | 低轨卫星星座构型保持方法 | |
CN112240941B (zh) | 一种重力卫星星载加速度计的相对标校方法及系统 | |
CN108875174A (zh) | 一种基于多段打靶法的不变拟周期轨道确定方法 | |
CN110059285B (zh) | 考虑j2项影响的导弹自由段弹道偏差解析预报方法 | |
Rutkowska et al. | SLR technique used for description of the Earth elasticity | |
CN114771873B (zh) | 一种超低轨卫星轨道自主精确维持方法 | |
CN103196451A (zh) | 一种基于Fisher信息矩阵的导航脉冲星选择方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |