[发明专利]适用于卫星轨道参数递推计算的去奇异方法

权利要求书

1.一种适用于卫星轨道参数递推计算的去奇异方法，其特征在于，所述适用于卫星轨道参数递推计算的去奇异方法利用以星历时刻t₀、以及t₀时刻的轨道平根数[a₀,e₀,i₀,Ω₀,ω₀,M₀]，得到卫星在t₁时刻在J2000.0坐标系下的轨道瞬根数[a_s,e_s,i_s,Ω_s,ω_s,M_s]；

其中，a₀为t₀时刻轨道半长轴，e₀为t₀时刻轨道偏心率，i₀为t₀时刻轨道倾角，Ω₀为t₀时刻升交点赤径，ω₀为t₀时刻近地点幅角，M₀为t₀时刻平近点角；a_s为t₁时刻轨道半长轴，e_s为t₁时刻轨道偏心率，i_s为t₁时刻轨道倾角，Ω_s为t₁时刻升交点赤径，ω_s为t₁时刻近地点幅角，M_s为t₁时刻平近点角。

2.如权利要求1所述的一种适用于卫星轨道参数递推计算的去奇异方法，其特征在于，所述适用于卫星轨道参数递推计算的去奇异方法包括以下步骤：

步骤S1，输入参数无奇异处理；

步骤S2，卫星轨道根数递推计算；

步骤S3，无奇点变量还原。

3.如权利要求2所述的一种适用于卫星轨道参数递推计算的去奇异方法，其特征在于，所述步骤S1具体包括以下步骤：

轨道偏心率e≈0时，为了避免计算中出现奇点，转化为3个无奇点变量ξ₀、η₀和λ₀，即：

ξ₀＝e₀cos(ω₀)

η₀＝-e₀sin(ω₀)

λ₀＝ω₀+M₀。

4.如权利要求3所述的一种适用于卫星轨道参数递推计算的去奇异方法，其特征在于，所述步骤S2具体包括以下步骤：

利用无奇点的输入参数，进行卫星轨道位置递推计算，得到t₁时刻的卫星无奇点轨道根数[a_s,i_s,Ω_s,ξ_s,η_s,λ_s]，其中，ξ_s,η_s,λ_s为递推后的无奇点变量。

5.如权利要求4所述的一种适用于卫星轨道参数递推计算的去奇异方法，其特征在于，所述步骤S3中采用以下公式将3个无奇点变量还原：

ω_s＝arccos(ξ_s/e_s)

M_s＝mod(λ_s-ω_s,2π)

得到t₁时刻卫星轨道瞬根数[a_s,e_s,i_s,Ω_s,ω_s,M_s]。