[發明專利]一種衛星位置計算方法及裝置

在審

• 申請號: 202011161320.7
• 申請日: 2020-10-26
• 公開（公告）號: CN112327333A
• 公開（公告）日: 2021-02-05
• 發明（設計）人: 車征
• 申請（專利權）人: 陜西星邑空間技術有限公司
• 主分類號: G01S19/27
• 分類號: G01S19/27;G01S19/37;G01S19/21
• 代理公司: 北京挺立專利事務所(普通合伙) 11265
• 代理人: 張智銳

1.一種衛星位置計算方法，其特征在于，所述方法包括：

獲取預設時間段內至少一個衛星軌道參數的變化規律；

根據所述衛星軌道參數的變化規律得到與所述衛星軌道參數對應的擬合函數，所述擬合函數中包括至少一個待擬合變量；

根據所述衛星軌道參數的變化規律和預設的待擬合變量初值計算所述擬合函數中的待擬合變量，得到不同時間下的衛星軌道參數；

根據所述不同時間下的衛星軌道參數，得到不同時間下衛星的位置。

2.根據權利要求1所述的方法，其特征在于，所述衛星軌道參數包括：半長軸、長軸偏心率、短軸偏心率、軌道傾角、升交點赤經和相位。

3.根據權利要求1所述的方法，其特征在于，所述獲取預設時間段內至少一個衛星軌道參數的變化規律，包括：

獲取預設時間段內根據動力學模型計算得到的至少一個衛星軌道參數的變化規律。

4.根據權利要求2所述的方法，其特征在于，當所述衛星軌道參數為半長軸時，所述與所述衛星軌道參數對應的擬合函數為：其中A為所述半長軸，A₀、A₁和A₂分別為所述半長軸擬合函數中的待擬合變量；

當所述衛星軌道參數為長軸偏心率時，所述與所述衛星軌道參數對應的擬合函數為：其中e_x為所述長軸偏心率，e_x0、e_x1、e_x2、e_x3和e_x4分別為所述長軸偏心率擬合函數中的待擬合變量；

當所述衛星軌道參數為短軸偏心率時，所述與所述衛星軌道參數對應的擬合函數為：其中e_y為所述短軸偏心率，e_y0、e_y1、e_y2、e_y3和e_y4分別為所述短軸偏心率擬合函數中的待擬合變量；

當所述衛星軌道參數為軌道傾角時，所述與所述衛星軌道參數對應的擬合函數為：其中i為所述軌道傾角，i₀、i₁、i₂、i₃和i₄分別為所述軌道傾角擬合函數中的待擬合變量；

當所述衛星軌道參數為升交點赤經時，所述與所述衛星軌道參數對應的擬合函數為：其中Ω為所述升交點赤經，Ω₀、Ω₁、Ω₂和Ω₃分別為所述升交點赤經擬合函數中的待擬合變量；

當所述衛星軌道參數為相位時，所述與所述衛星軌道參數對應的擬合函數為：其中λ為所述相位，λ₀、λ₁、λ₂和λ₃分別為所述相位擬合函數中的待擬合變量。

5.根據權利要求1所述的方法，其特征在于，所述根據所述衛星軌道參數的變化規律和預設的待擬合變量初值計算所述擬合函數中的待擬合變量，包括：

根據預設原則設置所述擬合函數中的待擬合變量的初值，得到第一函數；

將所述衛星軌道參數的變化規律中的多個時間數值和多個時間數值對應的多個衛星軌道參數的數值代入所述擬合函數中，得到多個遞推函數；

根據所述第一函數和所述多個遞推函數計算所述擬合函數中的待擬合變量。

6.根據權利要求5所述的方法，其特征在于，所述根據所述第一函數和所述多個遞推函數計算所述擬合函數中的待擬合變量，包括：

將第一遞推函數與所述第一函數做差，得到關于所述待擬合變量的第一差值數組，利用所述第一差值數組修正所述待擬合變量的初值，得到第一待擬合變量，并將所述第一待擬合變量帶入所述擬合函數中，得到第二函數，將所述第二遞推函數與所述第二函數做差，并重復上述過程，直至第n差值數組中的任一數值不再變化或所述第n差值數組中的任一數值與0的差值小于預設閾值，則所述第n差值數組對應的第n待擬合變量即為所述擬合函數中的待擬合變量。