本发明专利技术的一种适用于卫星轨道参数递推计算的去奇异方法，方便递推算法根据当前轨道平根数递推卫星轨道信息，预测指定时刻卫星的位置，包括如下步骤：输入参数无奇异处理；卫星轨道根数递推计算；无奇点变量还原。该方法可以将输入的卫星轨道平根数进行去奇异，避免当偏心率接近零时可能出现的奇异问题，使算法适合长时间、高精度的卫星位置预测，可以适用于卫星在轨的多种任务需求。星在轨的多种任务需求。星在轨的多种任务需求。
【技术实现步骤摘要】
适用于卫星轨道参数递推计算的去奇异方法
[0001]本专利技术涉及卫星轨道动力学中的卫星轨道计算领域，具体涉及卫星轨道参数递推领域。
技术介绍
[0002]近年来，对于人造卫星在指定时刻位置的预测，已经成为越来越重要的问题。高精度的轨道预报作为航天技术中一项重要技术，在卫星轨道设计和轨道优化中起着至关重要的作用，同时可以为天线指向、跟踪定位等卫星在轨任务提供可靠的轨道信息参考。传统的卫星轨道参数递推需要以轨道平根数作为输入，当轨道偏心率接近零时，计算过程中可能会出现奇异问题。
[0003]目前，应用于地面站的卫星轨道递推方法和高精度定轨算法一般由高性能计算机实现，其中包括了高精度积分算法和高精度动力学模型，对计算资源需求较多。卫星在轨运行时，卫星定轨的一种常见手段是由GNSS接收机接收其轨道数据。这种方法的缺点是：GNSS实时定轨存在误差，且易由于卫星指向、电子干扰或环境遮挡等情况造成导航星可见数量小于4颗，这种情况会影响卫星轨道数据的连续性和准确性。另一种手段是由地面上注轨道根数进行外推。这种方法的缺点是：目前的应用的外推模型一般过于简单，导致轨道数据精度变差，且一般不适用于长期预报。
[0004]中国专利技术专利“一种卫星位置获取方法和系统”(专利号：CN104459732A)中，介绍了一种GLONASS卫星位置解算方法。该方法可以获取观测时刻卫星的位置和速度，但是采用是一种数值积分方法，并且依赖于GLONASS接收机。
[0005]中国专利技术专利“一种基于拉格朗日插值及卡尔曼滤波的轨道数据插值方法”(专利号：CN103268407A)中，介绍了一种基于拉格朗日插值及卡尔曼滤波的轨道数据插值方法。该算法可以实现卫星位置的高精度实时仿真，但是占用计算资源较大，会增加星上计算机的负担。
[0006]李丹，于洋等人在“基于轨道根数的低轨卫星轨道预测算法”(见《光学精密工程》，2016年，10期)论文中提出了一种利用椭圆曲线来预测卫星轨道的方法，但是在求解过程中需要计算系数的偏微分。
[0007]中国专利技术专利“一种适用于圆轨道卫星的星上自主轨道外推方法”(专利号：CN103995800A)中，介绍了一种适用于圆轨道卫星的轨道递推的方法。但是该方法只考虑了摄动。
[0008]基于以上考虑，本专利技术公开的适用于卫星轨道参数递推计算的去奇异方法，其优势在于，使算法精度高，适用于长期递推，既能满足多种在轨任务的精度需求，又不会占用星上计算机过多资源。
技术实现思路
[0009]针对卫星在轨时进行轨道递推计算时，若轨道偏心率很小(近圆轨道，e≈0)时，计
算中出现会奇点的问题，本专利技术提出了适用于卫星轨道参数递推计算的去奇异方法，该方法可以将轨道平根数转化为无奇点根数进行卫星轨道递推计算。
[0010]为了实现该方法，采用如下的技术方案：
[0011]适用于卫星轨道参数递推计算的去奇异方法，所述的预报方法以星历时刻t0、t0时刻的轨道平根数[a0,e0,i0,Ω0,w0,M0]和目标时刻t1作为输入，对卫星轨道根数进行去奇异处理。计算结束后，将无奇点根数还原。其中，a0为t0时刻轨道半长轴，e0为t0时刻轨道偏心率，i0为t0时刻轨道倾角，Ω0为t0时刻升交点赤径，ω0为t0时刻近地点幅角，M0为t0时刻平近点角；
[0012]本专利技术步骤具体如下：
[0013]步骤S1，输入参数无奇异处理
[0014]当轨道偏心率很小(近圆轨道，e≈0)时，为了避免计算中出现奇点，引入3个无奇点变量ξ0、η0和λ0，即
[0015]ξ0＝e0cos(ω0)
[0016]η0＝
‑
e0sin(ω0)
[0017]λ0＝ω0+M0[0018]步骤S2，轨道递推计算
[0019]利用无奇点的输入参数，进行卫星轨道位置递推计算，得到t1时刻的卫星无奇点轨道根数[a
s
,i
s
,Ω
s
,ξ
s
,η
s
,λ
s
]。
[0020]ξ
s
,η
s
,λ
s
为递推后的无奇点变量。
[0021]步骤S3，无奇点变量还原
[0022]计算结束后将3个无奇点变量还原
[0023][0024]ω
s
＝arccos(ξ
s
/e
s
)
[0025]至此通过轨道递推方法得到t1时刻卫星轨道瞬根数[a
s
,e
s
,i
s
,Ω
s
,ω
s
,M
s
]。
[0026]其中，a
s
为t1时刻轨道半长轴，e
s
为t1时刻轨道偏心率，i
s
为t1时刻轨道倾角，Ω
s
为t1时刻升交点赤径，ω
s
为t1时刻近地点幅角，M
s
为t1时刻平近点角
[0027]以上为适用于卫星轨道参数递推计算的去奇异方法的
技术实现思路
，该方法适用于长期递推，既能满足多种在轨任务的精度需求，又不会占用星上计算机过多资源。
附图说明
[0028]图1为适用于卫星轨道参数递推计算的去奇异方法流程图。
[0029]图2为卫星的星上轨道参数自主预报示意图。
具体实施方式
[0030]下面结合附图对本专利技术的具体实施方式进行进一步的详细描述。
[0031]卫星的位置预报可以为卫星在轨任务提供可靠的轨道信息参考，卫星的位置可以用轨道根数表示，也可以用坐标表示。本专利技术包括以下步骤(如图1所示)：步骤S1，输入参数无奇异处理；步骤S2，卫星轨道根数递推计算；步骤S3，无奇点变量还原。
[0032]本专利技术的具体方案是根据给定的初始星历时刻t0和t0时刻的轨道平根数[a0,e0,i0,Ω0,ω0,M0]，在给定目标时刻t1的情况下，进行去奇异处理，方便递推算法通过数学计算递推得到卫星在t1时刻的轨道瞬根数。6个轨道根数的物理意义如下表所示。
[0033]a0e0i0Ω0ω0M0轨道半长轴轨道偏心率轨道倾角升交点赤径近地点幅角平近点角
[0034]以FY
‑
3D星在J2000.0坐标系下在轨运行实测数据作为参考(如图2所示)，进行仿真分析，在轨运行实测数据如下
[0035]北京时间2019年1月4日16时00分00秒FY
‑
3D星实测轨道数据
[0036][0037][0038]北京时间2019年1月5日16时00分00秒FY
‑
3D星实测轨道数据
[0039][0040]将北京时间2019年1月4日16时00分00秒FY
‑
3D星实测轨道平根数作为输入，引入3个无奇点变量，即
[0本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种适用于卫星轨道参数递推计算的去奇异方法，其特征在于，所述适用于卫星轨道参数递推计算的去奇异方法利用以星历时刻t0、以及t0时刻的轨道平根数[a0,e0,i0,Ω0,ω0,M0]，得到卫星在t1时刻在J2000.0坐标系下的轨道瞬根数[a
s
,e
s
,i
s
,Ω
s
,ω
s
,M
s
]；其中，a0为t0时刻轨道半长轴，e0为t0时刻轨道偏心率，i0为t0时刻轨道倾角，Ω0为t0时刻升交点赤径，ω0为t0时刻近地点幅角，M0为t0时刻平近点角；a
s
为t1时刻轨道半长轴，e
s
为t1时刻轨道偏心率，i
s
为t1时刻轨道倾角，Ω
s
为t1时刻升交点赤径，ω
s
为t1时刻近地点幅角，M
s
为t1时刻平近点角。2.如权利要求1所述的一种适用于卫星轨道参数递推计算的去奇异方法，其特征在于，所述适用于卫星轨道参数递推计算的去奇异方法包括以下步骤：步骤S1，输入参数无奇异处理；步骤S2，卫星轨道根数递推计算；步骤S3，无奇点变量还原。3.如权利要求2所述的一种适用于卫星轨道参数递推计算的去奇异方法，其特征在于，所述步骤S1具体包括以下步骤：轨道偏心率e≈0时，为了避免计算中...
【专利技术属性】
技术研发人员：刘伟亮，吕旺，宋玉亭，王豪，钱元元，余润清，肖灿明，俞航，顾军，张健军，沈松涛，辛阁霖，
申请(专利权)人：上海航天空间技术有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。