劉慶博，任順清

（哈爾濱工業(yè)大學(xué) 空間控制與慣性技術(shù)研究中心，哈爾濱 150080）

空間目標(biāo)初始軌道確定吉布斯問題誤差分析

劉慶博，任順清

（哈爾濱工業(yè)大學(xué) 空間控制與慣性技術(shù)研究中心，哈爾濱 150080）

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

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

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

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

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

1 軌道根數(shù)與初始軌道確定

1.1 軌道根數(shù)

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

求解上述軌道根數(shù)的過程如下：

徑向速度大小vr為

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

圖1 軌道根數(shù)Fig.1 Orbital parameters

軌道傾角i：

升交點赤經(jīng)Ω：

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

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

近地點幅角ω：

真近點角θ：

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

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

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

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

其中，

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

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

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

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

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

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

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

NDr

2

將比角動量h對r1、r2和r3求導(dǎo)，則h的誤差與位置誤差δ1r、δr2和δr3的關(guān)系為

對位置矢量求導(dǎo)，得：

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

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

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

類似，略去詳細(xì)推導(dǎo)。

速度的靈敏度矩陣Av為

其中，

長半軸a的靈敏度矩陣為

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

升交點赤經(jīng)Ω的靈敏度矩陣為

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

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

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

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

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

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

利用STK軟件繪制軌道，設(shè)定軌道參數(shù)如下：

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

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

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

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

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

圖5 升交點赤經(jīng)Ω與位置誤差的關(guān)系Fig.5 Relationship between right ascension of the ascending node and position error

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

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

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

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

4 結(jié) 論

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

（References）：

[1] 劉林, 王海紅, 胡松杰. 衛(wèi)星定軌綜述[J]. 飛行器測控學(xué)報, 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] 李振華, 鮮勇, 劉炳琪, 等. 單站單圈衛(wèi)星初軌確定方法研究[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] 王歆. 一種初軌計算的穩(wěn)健方法[J]. 天文學(xué)報, 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接收機定位誤差對衛(wèi)星軌道根數(shù)的影響分析[J]. 測繪科學(xué), 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自適應(yīng)偽譜法的飛行器再入軌跡優(yōu)化與制導(dǎo)[J]. 中國慣性技術(shù)學(xué)報, 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

裝備預(yù)研基金項目（9140A09030313HT01121）；國家重大科學(xué)儀器設(shè)備開發(fā)專項（2013YQ310737）

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

聯(lián) 系 人：任順清（1967—），男，教授，博士生導(dǎo)師。E-mail: renshunqing@hit.edu.cn