[發(fā)明專利]一種空間飛行器軌道特征點(diǎn)的解析計(jì)算方法

有效

• 申請?zhí)枺?/td>: 201810069338.0
• 申請日: 2018-01-24
• 公開（公告）號: CN108334683B
• 公開（公告）日: 2022-02-18
• 發(fā)明（設(shè)計(jì)）人: 孫遠(yuǎn)輝;趙倩
• 申請（專利權(quán)）人: 北京電子工程總體研究所
• 主分類號: G06F30/15
• 分類號: G06F30/15
• 代理公司: 北京正理專利代理有限公司 11257
• 代理人: 付生輝

1.一種空間飛行器軌道特征點(diǎn)的解析計(jì)算方法，其特征在于，包括：

變換當(dāng)前時間點(diǎn)的軌道參數(shù)，獲取第一軌道六根數(shù)；

根據(jù)飛行器當(dāng)前參數(shù)和二體理論下的軌道特征點(diǎn)計(jì)算得到飛行時間；

結(jié)合飛行時間，利用基于中間軌道理論的外推算法計(jì)算得到偽特征點(diǎn)的參數(shù)；

根據(jù)二體理論變換偽特征點(diǎn)的參數(shù)，獲取第二軌道六根數(shù)；

根據(jù)偽特征點(diǎn)的參數(shù)和第二軌道六根數(shù)，計(jì)算得到真特征點(diǎn)參數(shù)；

其中，

結(jié)合飛行時間，利用基于中間軌道理論的外推算法計(jì)算得到偽特征點(diǎn)的參數(shù)，其推算步驟包括：

S31：將當(dāng)前時刻J2000坐標(biāo)系的矢量x(t_i)，轉(zhuǎn)換為橢球坐標(biāo)系；

具體的，進(jìn)行如下轉(zhuǎn)換，

其中x(t_i)由x_i，y_i，z_i，和組成，并且

d＝(r_i²-c²)+δ(2z_i+δ)

r_i²＝x_i²+y_i²+z_i²

D²＝(ρ_i²+c²)(1-η_i²)

S32：計(jì)算Jacobi前三個常數(shù)(α₁，α₃，α₂)；

S33：計(jì)算F(ρ)和G(η)的四次式；

F(ρ)＝μ[c²p₀(1-S₀)+(ρ²+c²)(γ₀ρ²+2ρ-p₀)]

F(ρ)＝μγ₁(γρ²+2ρ-p)(ρ²+2A₁ρ+B₁)

G(η)＝μ[-p₀(1-S₀)+(1-η²)(p₀+2δη+c²γ₀η²)]

G(η)＝μS₁p₀(S+2Pη-η²)(1+P₁η-Q₁η²)

S34：初始化六個積分系數(shù)(R₁，R₂，R₃，N₁，N₂，N₃)；

計(jì)算R₁，R₂，R₃，需要A_k和W_k，其中A₀＝1，通過因數(shù)分解得到A₁及其它，W₀＝V₀＝W為真近地點(diǎn)角，W₁＝(W+eV₁)/p，V₁＝sinW；

N₂＝D₁[u+k₁T₂/2+k₁T₂/2+3k₁²T₄8+5k₁³T₆16]

計(jì)算N₁，N₂，N₃，需要C_1k，C_2k和T_k，其中

T₀＝u

T₁＝1-cosu

T_k＝[(k-1)T_k-2-cosusin^k-1u]/k,k＝2,...,6

S35：計(jì)算Jacobi后三個常數(shù)(β₁，β₂，β₃)

具體的，采用下式計(jì)算：

t_i+β₁＝R₁(ρ_i)+c²N₁(η_i)

β₂＝-α₂R₂(ρ_i)+α₂N₂(η_i)

β₃＝φ_i+c²α₃R₃(ρ_i)-α₃N₃(η_i)

使用初始環(huán)境和已初始化的系數(shù)，設(shè)置，

β₁＝-T時間

β₂＝ω近地點(diǎn)幅角

β₃＝Ω升交點(diǎn)赤經(jīng)

S36：使用Jacobi常數(shù)(α₁，α₃，α₂，β₁，β₂，β₃)替換運(yùn)動方程，求解給定時間t_f時的ρ_f，η_f和φ_f；

具體的，采用下式計(jì)算：

t_f+β₁＝R₁(ρ_f)+c²N₁(η_f)

β₂＝-α₂R₂(ρ_f)+α₂N₂(η_f)

β₃＝φ_i+c²α₃R₃(ρ_f)-α₃N₃(η_f)

第一個運(yùn)動方程是開普勒方程的通式；只要計(jì)算使用橢球坐標(biāo)系，相應(yīng)的值由下式計(jì)算：

其中

S37：將時間t_f時橢球坐標(biāo)系下的位置矢量X(t_f)，轉(zhuǎn)換為J2000坐標(biāo)系的矢量x(t_f)

其中

所述真特征點(diǎn)的參數(shù)求解模型為：

其中，

r＝a(1-emx cosΔE+emysinΔE)

其中，ΔE為偏近點(diǎn)角增量，滿足Kepler方程：