空間目標初始軌道確定吉布斯問題誤差分析

劉慶博，任順清

（哈爾濱工業大學 空間控制與慣性技術研究中心，哈爾濱 150080）

空間目標初始軌道確定吉布斯問題誤差分析

劉慶博，任順清

（哈爾濱工業大學 空間控制與慣性技術研究中心，哈爾濱 150080）

現代戰爭中的空間武器平臺捕獲敵方的空間目標，確定其運行軌道后實施軍事打擊，所以對軌道確定參數的誤差進行分析對于精確打擊具有非常重要的意義。為了準確分析吉布斯三位置矢量定軌法的軌道根數的確定精度，根據吉布斯方法的軌道確定模型，采用向量求導的方法，詳細推導了軌道根數對于觀測位置誤差的靈敏度矩陣，并給出了軌道根數誤差與觀測位置誤差之間的關系，明確初始軌道的確定精度及誤差傳播規律。最后對吉布斯方法及其誤差分析進行仿真分析驗證，驗證了方法的正確性。

初始軌道確定；吉布斯方法；軌道根數；靈敏度矩陣

在軌道力學中，對空間目標進行軌道確定包括兩個過程：利用短弧段觀測數據的初始軌道確定以及長弧段下的精密定軌[1]。初始軌道確定一般都是采用二體模型，經典的初始軌道確定方法都是利用測角數據來進行計算的，主要有Laplace法和Gauss法兩種。

長期以來，國內外的眾多學者都針對以上兩種方法進行了更加深入的研究，并提出了多種初始軌道確定的方法，解決了一些使用傳統方法會出現不收斂或不穩定等問題。文獻[2]基于單站單圈測量數據研究了低軌衛星初軌確定的方法，對遺傳算法和最小二乘法在不同觀測條件下確定的初軌信息以及軌道預測信息進行了分析和比較，確定了兩種方法的優缺點和各自的適用條件。文獻[3]采用最小一乘方法建立了一種初軌計算的穩健方法，將初軌計算問題轉換為線性規劃問題求解，并通過bootstrap方法給出估計精度，數值計算結果表明該方法穩健有效，并具有較高的崩潰點。文獻[4]將采用角度和角速率數據進行初始軌道確定與只采用角度數據進行了對比，采用角度和角速率數據有著更好的平穩性和可靠性。文獻[5]提出了一種采用單時刻 GPS測量的方法來對地球靜止衛星進行初始軌道的確定，并采用EKF對該算法進行了驗證。文獻[6]利用遺傳算法來對極短弧的光學觀測進行初始軌道確定，并利用蒙特卡洛仿真來驗證算法的可靠性。文獻[7]將天基角度測量與遺傳算法相結合來確定初始軌道，得到了滿意的結果，并且多種群遺傳算法（MPGA）可以有效解決角度測量方法的迭代過程中出現的問題。文獻[8]中根據開普勒運動定律推導出確定斜距的方法，5個數據一組進行初始軌道的確定。除此之外，還可以利用其他方法來確定初始軌道，如利用高斯混合模型求解[9]。以上是對初始軌道確定方法研究現狀的概述，對軌道確定方法進行誤差分析的相關文獻相對較少。文獻[10]中分析了星載GPS接收機的定位誤差對衛星軌道根數的影響，只推導了長半軸，軌道傾角和升交點赤經的誤差與位置、速度誤差的關系式。

本文將主要對吉布斯三位置矢量定軌法進行誤差分析。首先給出軌道根數的計算方法和吉布斯三位置矢量定軌法。利用3個位置矢量就能確定軌道根數，但目前沒有文獻給出位置矢量本身的測量誤差、夾角的大小等對于6個軌道根數確定誤差的影響程度。在整個軌道確定的過程中，初始軌道確定的精度大多數情況下無法滿足精度要求，需要進行精密定軌，而初始軌道確定的精度往往會影響精密定軌初值的選取，高精度的初始軌道的確定可以減少整個軌道確定的時間。所以本文將對吉布斯方法進行誤差分析，推導出6個軌道根數的靈敏度矩陣，明確誤差傳播特性，并對其進行仿真分析與驗證。

1 軌道根數與初始軌道確定

1.1 軌道根數

軌道根數是用來描述天體在其軌道運行狀態的一組參數。6個軌道根數如圖1所示，其中，h為比角動量的模，i為軌道傾角，Ω為升交點赤經(RA)，e為偏心率，ω為近地點幅角，θ為真近點角，a為長半軸。當已知位置矢量和速度矢量

求解上述軌道根數的過程如下：

徑向速度大小vr為

根據文獻[11]，直接給出6個軌道根數的算法，比角動量h：

圖1 軌道根數Fig.1 Orbital parameters

軌道傾角i：

升交點赤經Ω：

偏心率矢量e為地心O指向近地點的矢量，表達式為

偏心率e為偏心率矢量e的模。

近地點幅角ω：

真近點角θ：

長半軸a：式中：μ = GMe是地心引力常數，μ = 398600km3/s2，Me是地球質量。

下面針對a、i、Ω、e、ω、θ這6個軌道根數進行誤差分析。

1.2 初始軌道確定——吉布斯三位置矢量定軌法

如果測量出空間目標軌道上P1、P2、P3位置的位置矢量 r1、r2和r3后，便可用吉布斯三位置矢量定軌法來計算6個軌道根數。如圖2所示。

其中，

根據式(9)，角動量h可由r1、r2和r3通過如下關系得到：

圖2 吉布斯三位置矢量定軌Fig.2 Gibbs’ method for initial orbit determination

從式(13)可以看出，測量得到3個位置矢量，能計算過度矢量N、D、S及其相應的標量，再計算速度矢量，最后根據式(1)～(8)確定6個軌道根數。然而測量的矢量r1、r2和r3是有誤差的，而且它們之間的夾角有大有小，如何根據位置矢量的誤差來確定6個軌道根數的誤差就是下面將要研究的內容。

2 初始軌道確定方法的誤差分析

對吉布斯三位置矢量定軌法進行誤差分析，要根據6個軌道根數的計算公式寫成關于3個位置矢量r1、r2和的微分表達式，求取關于 r1、r2和r3的偏導，轉換成矩陣的形式，得到6個軌道根數a、i、Ω、e、ω、θ的靈敏度矩陣，并且h作為重要的中間變量同樣需要求取靈敏度矩陣。在求解過程中，主要利用向量對向量求偏導的相關計算公式以及如下幾個常用的求取偏導公式：

1）如果M為矢量M的模，m為另一矢量，則：

首先求取比角動量的模 h對位置矢量的靈敏度矩陣。根據軌道比角動量為常矢量，h=r×v，將第2個位置的速度v2的表達式代入可以得到

NDr

2

將比角動量h對r1、r2和r3求導，則h的誤差與位置誤差δ1r、δr2和δr3的關系為

對位置矢量求導，得：

根據式(15)，寫出矩陣：

式中，比角動量的靈敏度矩陣

為3×9維的矩陣，其中各個元素的計算可將公式(16)～(24)代入即可計算，它們均是矢量 r1、r2和r3及其標量的函數。計算比角動量的模的靈敏度矩陣為其他5個軌道根數的靈敏度矩陣求取方法

類似，略去詳細推導。

速度的靈敏度矩陣Av為

其中，

長半軸a的靈敏度矩陣為

軌道傾角i的靈敏度矩陣為

升交點赤經Ω的靈敏度矩陣為

偏心率e的靈敏度矩陣由偏心率矢量e的靈敏度矩陣得到，則：

真近點角θ的靈敏度矩陣Aθ為

求取近地點幅角的靈敏度矩陣前先求解出矢量B關于位置誤差的靈敏度矩陣ΑB為

繼而求解出近地點幅角的靈敏度矩陣Αω為

標量對位置矢量的靈敏度矩陣均為1×9維矩陣，假設各個位置矢量測量的不確定度相等，均為σ，并假設各次測量均為獨立測量，再結合不確定度的合成公式，即當所求取軌道根數的靈敏度矩陣為A時，那么軌道根數的不確定度為表示靈敏度矩陣的元素，就可求出6個軌道根數的不確定度。

3 初始軌道確定方法的仿真分析

利用STK軟件繪制軌道，設定軌道參數如下：

在繪制好的軌道上選取3個位置矢量（km）：

對吉布斯三位置矢量定軌法進行Matlab編程，經過計算可以得到6個軌道根數，長半軸軌道傾角偏心率 e= 0.100132，升交點赤經近地點幅角真近點角

把以上建立的運行軌道模型中的位置信息當作基準數據，在基準數據中加入0.1 km的位置隨機誤差，多次取值，求取6個軌道根數的不確定度，并在此基礎之上逐步增大選取的 3個位置矢量之間的夾角，最大到120°左右，觀察6個軌道根數不確定度的變化趨勢，繪（相對于位置矢量r）。制出變化曲線。令位置不確定度根據前面得到的6個軌道根數的靈敏度矩陣，計算6個軌道根數的不確定度，與由仿真數據計算得到的不確定度進行比較，仿真結果如圖3～8所示。

圖3 長半軸a與位置誤差的關系Fig.3 Relationship between semimajor axis and position error

圖4 軌道傾角i與位置誤差的關系Fig.4 Relationship between inclination and position error

圖5 升交點赤經Ω與位置誤差的關系Fig.5 Relationship between right ascension of the ascending node and position error

圖6 偏心率e與位置誤差的關系Fig.6 Relationship between eccentricity and position error

從上述仿真結果可以看到，利用吉布斯三位置矢量定軌法計算得到的軌道根數與預設的軌道根數接近，由實際數據計算得到的軌道根數的不確定度與由靈敏度矩陣計算得到的軌道根數的不確定度，兩者的結果十分接近，驗證了前面所求取的軌道根數的靈敏度矩陣的準確性。從軌道根數的不確定度的變化趨勢可以看到，增大選取位置矢量之間的夾角，6個軌道根數的不確定度在減小，所以盡量增大選取的3個位置矢量之間的夾角有利于減小軌道根數的誤差。但實際確定初始軌道時，增大3個位置的夾角可能受到所需確定時間的限制，可根據具體所需的軌道確定精度和效率進行取舍。

圖7 真近點角θ與位置誤差的關系Fig.7 Relationship between true anomaly and position error

圖8 近地點幅角ω與位置誤差的關系Fig.8 Relationship between argument of perigee and position error

4 結 論

針對吉布斯方法確定的6個軌道根數的誤差，推導了6個軌道根數誤差相對于位置矢量測量誤差的靈敏度矩陣，用Matlab軟件編程計算了6個軌道根數的不確定度；然后，利用STK軟件繪制軌道，驗證了空間目標初始軌道確定方法的有效性；最后，對選取位置矢量之間的夾角做出改變進行了仿真分析，繪制了軌道根數不確定度的變化曲線，將真實數據計算得到的不確定度與由靈敏度矩陣計算得到的不確定度進行比較，驗證了前面所求取的靈敏度矩陣的準確性。由仿真結果可知，為減小軌道根數的誤差，可以通過增大選取的3個位置矢量之間的夾角來實現。在實際的初始軌道確定中，可根據本文給出的6個軌道根數靈敏度矩陣可以計算出所求解的軌道的軌道根數誤差。

（References）：

[1] 劉林, 王海紅, 胡松杰. 衛星定軌綜述[J]. 飛行器測控學報, 2005, 24(2): 28-34. Liu Lin, Wang Hai-hong, Hu Song-jie. Summary on satellite orbit determination[J]. Journal of Spacecraft TT&C Technology, 2005, 24(2): 28-34.

[2] 李振華, 鮮勇, 劉炳琪, 等. 單站單圈衛星初軌確定方法研究[J]. 火力與指揮控制, 2014, 39(12): 123-126. Li Zhen-hua, Xian Yong, Liu Bing-qi, et al. Study on initial orbit determination with single lap radar station[J]. Fire Control & Command Control, 2014, 39(12): 123-126.

[3] 王歆. 一種初軌計算的穩健方法[J]. 天文學報, 2013, 54(3): 274-281. Wang Xin. A robust method of preliminary orbit determination[J]. Acta Astronomica Sinica, 2013, 54(3): 274-281.

[4] Demars K J, Jah M K, Schumacher P W. Initial Orbit Determination using Short-Arc Angle and Angle Rate Data[J]. IEEE Transactions on Aerospace & Electronic Systems, 2012, 48(3): 2628-2637.

[5] Ghangho K, Chongwon K, Changdon K. Coarse initial orbit determination for a geostationary satellite using single-epoch GPS measurements[J]. Sensors, 2015, 15(4): 7878-7897.

[6] Ansalone L, Curti F. A genetic algorithm for initial orbit determination from a too short arc optical observation[J]. Advances in Space Research, 2013, 52: 477-489.

[7] Liu L, Tang G, Hu S. Initial orbit determination based on sparse space-based angle measurement and genetic algorithm[C]//2013 IEEE International Conference on Information and Automation. 2013: 850-855.

[8] Vitarius P J, Hahs D, Gregory D A. A Keplerian approach to angles-only orbit determination[C]//Defense and Security Symposium. International Society for Optics and Photonics, 2006: 62200M/1-62200M/7.

[9] Demars K J, Jah M K. Probabilistic initial orbit determination using Gaussian mixture models[J]. Journal of Guidance Control & Dynamics, 2013, 36(5): 1324-1335.

[10] 初東, 王剛. 星載 GPS接收機定位誤差對衛星軌道根數的影響分析[J]. 測繪科學, 2014, 39(6): 31-34. Chu Dong, Wang Gang. Influence of spaceborne GPS receiver positioning errors on satellite’s orbital elements[J]. Science of Surveying and Mapping, 2014, 39(6): 31-34.

[11] Curtis H D. Orbital mechanics for engineering students [M]. Burlington: Butterworth-Heinemann, 2010.

[12] 夏紅偉, 李秋實, 李莉, 等. 基于hp自適應偽譜法的飛行器再入軌跡優化與制導[J]. 中國慣性技術學報, 2015, 23(6): 818-823. Xia Hong-wei1, Li Qiu-shi, Li Li, et al. Trajectory optimization and guidance for reentry craft based on hp-adaptive pseudospectral method[J]. Journal of Chinese Inertial Technology, 2015, 23(6): 818-823.

Error analysis of Gibbs’ method for initial orbit determination of space target

LIU Qing-bo, REN Shun-qing

(Space Control and Inertial Technology Research Center, Harbin Institute of Technology, Harbin 150080, China)

In modern war, the space weapon platform can capture enemy’s space target and then attack it after determining the orbit, so the error analysis of orbital parameters is of great significance. To precisely analyze the orbit-determination accuracy of orbital elements by Gibbs three-position vector method, the sensitivity matrices of orbit parameters versus the position errors are deduced by using the vector derivatives according to the model of Gibbs’ method. The precision of initial orbit determination and error propagation characteristics are determined by the relationship between the orbital parameter errors and position errors. Simulation results show the correctness of the Gibbs method and its error analysis.

initial orbit determination; Gibbs’ method; orbital parameter; sensitivity matrix

U666.1

A

1005-6734(2016)02-0263-06

10.13695/j.cnki.12-1222/o3.2016.02.023

2015-12-14；

2016-03-11

裝備預研基金項目（9140A09030313HT01121）；國家重大科學儀器設備開發專項（2013YQ310737）

劉慶博（1991—），男，博士研究生，從事慣性技術研究。E-mail: lqb9104@163.com

聯 系 人：任順清（1967—），男，教授，博士生導師。E-mail: renshunqing@hit.edu.cn