，，，

1. 中國(guó)科學(xué)院 國(guó)家空間科學(xué)中心，北京 100190 2. 中國(guó)科學(xué)院大學(xué)，北京 100190

在國(guó)內(nèi)外編目實(shí)踐中，一些突發(fā)的情況會(huì)對(duì)正常的編目工作產(chǎn)生不利的影響，例如太陽(yáng)爆發(fā)干擾了大氣，使得空間目標(biāo)的預(yù)報(bào)精度大大降低，也降低了對(duì)近地目標(biāo)的編目維持能力，由于極端空間環(huán)境條件下大氣模型精度有限，影響了大量的軌道預(yù)報(bào)精度，使得目標(biāo)丟失。一些突發(fā)的衛(wèi)星事件解體會(huì)在很短時(shí)間內(nèi)增加數(shù)以千計(jì)的可跟蹤目標(biāo)[1-3]。在這種情況下，設(shè)備監(jiān)測(cè)到的許多目標(biāo)可能沒(méi)有編目，非關(guān)聯(lián)軌道 (Uncorrelated Tracks，UCTs)數(shù)據(jù)和沒(méi)有編號(hào)的目標(biāo)大大增加，目標(biāo)編目數(shù)據(jù)處理中心的負(fù)擔(dān)立即增大。此外在目標(biāo)解體初期，目標(biāo)的空間密度還可能引起目標(biāo)混淆。這兩種情況必須在處理中心增大處理能力，以及增加額外的關(guān)聯(lián)過(guò)程來(lái)克服。大型衛(wèi)星星座的部署也會(huì)增加許多空間目標(biāo)。有些空間碎片可能已處于設(shè)備探測(cè)能力的邊緣，在觀測(cè)條件較好時(shí)，設(shè)備偶然能得到觀測(cè)數(shù)據(jù)，但是，觀測(cè)數(shù)據(jù)不足以維持該目標(biāo)的編目，因此目標(biāo)可能丟失。對(duì)于光學(xué)觀測(cè)設(shè)備，還有可見(jiàn)期的問(wèn)題[4-10]。

解決這些問(wèn)題除了增加探測(cè)網(wǎng)的探測(cè)能力及提高軌道預(yù)報(bào)精度外，還需要改進(jìn)關(guān)聯(lián)算法。采用基于協(xié)方差理論關(guān)聯(lián)算法可以提高短弧段雷達(dá)觀測(cè)數(shù)據(jù)和稀疏可見(jiàn)期數(shù)據(jù)的關(guān)聯(lián)水平。該計(jì)算方法旨在解決UCTs自動(dòng)關(guān)聯(lián)問(wèn)題，建立一個(gè)更完整的空間對(duì)象編目庫(kù)?；趨f(xié)方差的軌道關(guān)聯(lián)動(dòng)態(tài)算法是基于軌道不確定性協(xié)方差的動(dòng)態(tài)算法，利用軌道數(shù)據(jù)的協(xié)方差信息生成完整的軌道狀態(tài)和誤差分布函數(shù)，并用于統(tǒng)計(jì)關(guān)聯(lián)軌道。如果協(xié)方差準(zhǔn)確表示實(shí)際軌道誤差分布，那么這個(gè)算法是最優(yōu)的、最大化的概率關(guān)聯(lián)。協(xié)方差矩陣可以更有效統(tǒng)計(jì)和自動(dòng)方式關(guān)聯(lián)UCTs。基于協(xié)方差的關(guān)聯(lián)已經(jīng)在精密軌道攝動(dòng)的模擬中測(cè)試，可完成編目10萬(wàn)個(gè)空間目標(biāo)[11]。仿真結(jié)果表明，相對(duì)于之前固定門限方法的軌道關(guān)聯(lián)，基于協(xié)方差的軌道關(guān)聯(lián)可降低錯(cuò)誤關(guān)聯(lián)值達(dá)兩個(gè)數(shù)量級(jí)?；趨f(xié)方差的數(shù)據(jù)關(guān)聯(lián)方法將提高空間監(jiān)測(cè)網(wǎng)低軌道雷達(dá)監(jiān)測(cè)數(shù)據(jù)的編目分辨能力。

本文建立了線性和無(wú)跡卡爾曼濾波協(xié)方差生成和演化算法，分析了兩種協(xié)方差演化的精度，無(wú)跡卡爾曼濾波的協(xié)方差演化方法相對(duì)線性協(xié)方差演化方法有較高的演化精度。把協(xié)方差演化方法應(yīng)用到印度一箭104星中的20顆衛(wèi)星軌道關(guān)聯(lián)中，仿真結(jié)果表明采用基于協(xié)方差軌道關(guān)聯(lián)具有很好的關(guān)聯(lián)效果。

1 基于軌道協(xié)方差UCTs動(dòng)態(tài)關(guān)聯(lián)算法模擬分析方法

基于協(xié)方差的UCTs動(dòng)態(tài)關(guān)聯(lián)算法是把軌道協(xié)方差作為關(guān)聯(lián)時(shí)的影響因素，關(guān)聯(lián)量是復(fù)合等概率橢球。當(dāng)探測(cè)器的測(cè)量誤差是高斯分布并且軌道動(dòng)力學(xué)為線性時(shí)，協(xié)方差關(guān)聯(lián)是最佳關(guān)聯(lián)量。采用此關(guān)聯(lián)方法，相同坐標(biāo)系下兩衛(wèi)星的狀態(tài)量和協(xié)方差須預(yù)報(bào)到同一時(shí)刻。例如，可以從先前t1時(shí)刻的軌道和協(xié)方差預(yù)報(bào)到t2時(shí)刻軌道，兩軌道狀態(tài)差為:

X2,1=X2(t2)-X1(t2)

(1)

式中：X1(t2)和X2(t2)為空間目標(biāo)在t2時(shí)刻的軌道狀態(tài)，包括位置和速度。關(guān)聯(lián)過(guò)程中通過(guò)計(jì)算統(tǒng)計(jì)量k來(lái)衡量?jī)尚l(wèi)星的關(guān)聯(lián)程度。k屬于χ2統(tǒng)計(jì)，也被稱為馬氏距離：

k2=[X2,1]T[P1(t2)+P2(t2)]-1[X2,1]

(2)

式中：P1和P2分別為兩衛(wèi)星的軌道協(xié)方差矩陣。

為了實(shí)際模擬基于協(xié)方差的軌道關(guān)聯(lián)，需要模擬空間目標(biāo)監(jiān)測(cè)網(wǎng)的監(jiān)測(cè)模型和編目過(guò)程。決定從真實(shí)雷達(dá)測(cè)站出發(fā)，不同雷達(dá)測(cè)量數(shù)據(jù)隨機(jī)加入測(cè)量噪聲，模擬雷達(dá)測(cè)距及測(cè)角數(shù)據(jù)。定軌模塊對(duì)這些測(cè)量數(shù)據(jù)進(jìn)行初始軌道確定及最小二乘處理生成軌道及協(xié)方差數(shù)據(jù)。關(guān)聯(lián)模塊根據(jù)定軌目標(biāo)的協(xié)方差數(shù)據(jù)進(jìn)行演化實(shí)現(xiàn)與編目庫(kù)和非關(guān)聯(lián)目標(biāo)的關(guān)聯(lián)。

1.1 真實(shí)軌道

真實(shí)軌道是由數(shù)值軌道傳播器生成的。動(dòng)力學(xué)方程模型包括太陽(yáng)輻射光壓，空氣阻力密度模型，第三體的影響作用力和21×21階的重力場(chǎng)[12]。50個(gè)物體的軌道狀態(tài)取自軌道碎片編目。這些物體是從空間目標(biāo)目錄里半長(zhǎng)軸為6860～6890之間抽樣的。目標(biāo)的橫截面積設(shè)為常值，面積和質(zhì)量與碎片編目的值相同。設(shè)置每一個(gè)物體的阻力系數(shù)和光壓反射系數(shù)為默認(rèn)值?？臻g目標(biāo)軌道攝動(dòng)動(dòng)力學(xué)模型采用的是EGM2008的21階地球重力場(chǎng)模型，NRLMSISE2000大氣阻力模型，球形光壓模型。依據(jù)軌道攝動(dòng)方程演化計(jì)算理論上空間目標(biāo)的真實(shí)軌道[13]。

1.2 雷達(dá)探測(cè)器模擬

依據(jù)現(xiàn)有的雷達(dá)跟蹤測(cè)量站，模擬空間目標(biāo)的跟蹤測(cè)量數(shù)據(jù)。用國(guó)外公布的雷達(dá)參數(shù)對(duì)雷達(dá)跟蹤站進(jìn)行建模。表1所示這些參數(shù)來(lái)自美國(guó)巴利亞多(Vallado)[14]。對(duì)雷達(dá)的觀測(cè)數(shù)據(jù)進(jìn)行模擬，包括測(cè)距范圍、仰角、方位角、雷達(dá)高度、測(cè)距標(biāo)準(zhǔn)差、仰角標(biāo)準(zhǔn)差以及方位角標(biāo)準(zhǔn)差。在觀測(cè)數(shù)據(jù)中加入了高斯噪聲[15]。每一步都用標(biāo)準(zhǔn)差生成隨機(jī)偏差模擬高斯噪聲。采用這種標(biāo)準(zhǔn)差是因?yàn)楦鶕?jù)一段時(shí)間對(duì)雷達(dá)探測(cè)器的分析顯示日常偏差標(biāo)準(zhǔn)差與測(cè)量標(biāo)準(zhǔn)差大致相等。如果觀測(cè)數(shù)據(jù)超出了方位角和仰角的限制就要去掉[16]。假定某個(gè)空間目標(biāo)截面積為1 m2時(shí)觀測(cè)站最大測(cè)量距離為ρmax。雷達(dá)的返回值與1/ρ4成正比，其中ρ是測(cè)距，所以觀測(cè)量應(yīng)滿足：

(3)

式中：A為空間目標(biāo)截面積，兩次測(cè)量之間為10 s間隔時(shí)，觀測(cè)數(shù)據(jù)設(shè)置限制在2′弧段內(nèi)。每個(gè)雷達(dá)探測(cè)器被分配最大數(shù)量的可跟蹤目標(biāo)。

表1 雷達(dá)探測(cè)器參數(shù)

1.3 軌道確定協(xié)方差

為了測(cè)試k值是否可以用來(lái)關(guān)聯(lián)編目軌道，選取物體的一個(gè)子集部分構(gòu)成了模擬目錄。選用目錄里每個(gè)物體的第一軌道，使用測(cè)距、方位角、仰角三要素方法進(jìn)行初始軌道確定。軌道狀態(tài)估計(jì)作為最小二乘批處理程序的先驗(yàn)狀態(tài)。當(dāng)執(zhí)行最小二乘批差分校正時(shí)，每一個(gè)目錄物體所有前三天的觀測(cè)數(shù)據(jù)合并在一起。每個(gè)物體的面積和質(zhì)量都設(shè)為默認(rèn)值(1m2和100 kg)，但是要估計(jì)阻力系數(shù)及光壓反射參數(shù)的修正值[17-18]。在軌道計(jì)算中用數(shù)值計(jì)算軌道的傳播。在未來(lái)大于三天的時(shí)間里，物體估計(jì)狀態(tài)隨時(shí)間傳播且儲(chǔ)存為星歷表、用來(lái)做為一條空間物體目錄。50個(gè)空間物體只有30個(gè)創(chuàng)建了編目星歷，其他物體的星歷未知。

在先驗(yàn)?zāi)夸泟?chuàng)建后，所有的真實(shí)軌道在未來(lái)24 h傳播。用真實(shí)軌道可以形成觀測(cè)數(shù)據(jù)，且已知物體的每個(gè)軌道都被設(shè)置了編號(hào)。最小二乘差分校正用來(lái)獲得每個(gè)軌道的狀態(tài)和協(xié)方差估計(jì)。由于這些軌道相當(dāng)?shù)亩?120 s)，這些單個(gè)軌道的估計(jì)有很大的不確定性，這在協(xié)方差上得到了反映。

1.4 數(shù)值積分模型的軌道協(xié)方差演化

1.4.1 誤差演化線性化方法

(4)

(5)

對(duì)式(5)微分：

(6)

(7)

e(t)=Φ(t,ti)ei+wd(t,ti)

(8)

其中，狀態(tài)轉(zhuǎn)換矩陣Φ(t,ti)滿足：

(9)

線性協(xié)方差轉(zhuǎn)換矩陣Φ(t,ti)計(jì)算過(guò)程中忽略了高階項(xiàng)的影響，只考慮了一階項(xiàng)。線性協(xié)方差演化是一種分析方法，其優(yōu)點(diǎn)是計(jì)算簡(jiǎn)化高效。但是對(duì)于非線性程度高、協(xié)方差演化周期長(zhǎng)或初始協(xié)方差較大情況其計(jì)算精度較低[19]，不適用于非高斯協(xié)方差演化。

1.4.2 UKF協(xié)方差演化方法

UKF是基于確定性樣點(diǎn)計(jì)算的濾波算法，以Unscented變換為基礎(chǔ)，采用卡爾曼線性濾波框架和確定性采樣方式進(jìn)行濾波。UKF采用一組有效的采樣點(diǎn)來(lái)將狀態(tài)近似為高斯隨機(jī)變量，捕捉高斯隨機(jī)變量的均值和協(xié)方差。當(dāng)狀態(tài)量通過(guò)實(shí)際的非線性系統(tǒng)之后，后驗(yàn)均值和協(xié)方差可以精確到三階(泰勒展開(kāi))。由于不需要推導(dǎo)和計(jì)算Jacobian矩陣，基于UT的UKF算法更便于實(shí)現(xiàn)。

(1)構(gòu)造Sigma點(diǎn)

(10)

(2)利用Sigma點(diǎn)進(jìn)行非線性傳播

根據(jù)計(jì)算所得Sigma點(diǎn)直接按照離散化非線性模型進(jìn)行狀態(tài)和觀測(cè)預(yù)測(cè)，對(duì)應(yīng)產(chǎn)生的Sigma樣本點(diǎn)為：

χk(i)=f[χk-1(i),k]

(11)

式中：χk(i)(i=0，…,2n)為預(yù)測(cè)的對(duì)稱分布Sigma點(diǎn)。

(3)根據(jù)樣本點(diǎn)計(jì)算狀態(tài)均值和方差

(12)

(13)

從上面的3個(gè)步驟中可以看出，Unscented變換不需要計(jì)算非線性函數(shù)的雅可比矩陣，變換中僅涉及函數(shù)計(jì)算，即使系統(tǒng)的模型比較復(fù)雜的情況下，算法實(shí)現(xiàn)的難度也并不會(huì)增加。

1.5 UCTs關(guān)聯(lián)

處理這些軌道的第一步是比較星歷看軌道是否屬于已經(jīng)存在的編目物體，然后更新編目星歷。對(duì)于與已經(jīng)編目在案的軌道相符的軌道，采用軌道狀態(tài)值和協(xié)方差計(jì)算關(guān)聯(lián)量k值。如果k值小于限定閾值，則認(rèn)為該軌道匹配且軌道的標(biāo)識(shí)或者標(biāo)簽保持不變。如果k值大于限定閾值，則更改標(biāo)簽以表明它是不相關(guān)軌道即UCTs。然后就是把UCTs與所有編目里的物體進(jìn)行比較。為了減少比較數(shù)量，采用預(yù)測(cè)濾波的方法只選擇軌道根數(shù)大致相當(dāng)?shù)男菤v。如果比較值k小于限定閾值，則該軌道被標(biāo)記為屬于物體編目。接著，所有被分配到編目里的新軌道與先前確定的軌道結(jié)合形成新的星歷。這些綜合的觀測(cè)數(shù)據(jù)用來(lái)估計(jì)一個(gè)更新的編目軌道并創(chuàng)立一個(gè)新的、為期三天的預(yù)報(bào)編目星歷。下一步就是UCTs與UCTs之間的關(guān)聯(lián)，所有的UCTs形成一個(gè)候選軌道數(shù)據(jù)庫(kù)，不能與編目庫(kù)關(guān)聯(lián)的軌道再與UCTs候選軌道關(guān)聯(lián)。

2 模擬分析

2.1 協(xié)方差演化分析

低軌道目標(biāo)分別選取400 km軌道高度的近圓軌道目標(biāo)國(guó)際空間站和1 500 km軌道高度的近圓軌道目標(biāo)Ajisai衛(wèi)星。對(duì)國(guó)際空間站和Ajisai衛(wèi)星設(shè)定初值協(xié)方差，采用線性協(xié)方差演化方法和UKF協(xié)方差演化方法計(jì)算演化時(shí)間為一天時(shí)，空間站和衛(wèi)星的軌道協(xié)方差在RSW坐標(biāo)系下的分量。RSW坐標(biāo)系的R軸(徑向)是從地球指向航天器方向的矢量，S軸(橫向)在軌道平面內(nèi)垂直于R軸，W軸(垂跡)垂直于軌道面。把演化結(jié)果同STK軟件中高精度軌道(HPOP)協(xié)方差演化結(jié)果比較，分析兩種方法協(xié)方差演化精度。軌道攝動(dòng)力主要考慮地球引力、大氣阻力、光壓攝動(dòng)力及第三體引力。地球引力采用的是21階的EGM2008模型，大氣阻力采用的是NRLMSISE2000模型，大氣阻力系數(shù)Cd=2，光壓系數(shù)Cr=1.8。

假設(shè)在J2000坐標(biāo)系下空間站位置三方向初值誤差均為10 m，速度三方向誤差均為0.01 m/s?？臻g站協(xié)方差演化時(shí)長(zhǎng)為一天，選取演化時(shí)長(zhǎng)中的100個(gè)點(diǎn)，比較線性方法和UKF方法演化的100個(gè)點(diǎn)與HPOP協(xié)方差演化的誤差。圖1中前3個(gè)圖采用是兩種方法與HPOP對(duì)R、S、W三個(gè)方向的演化誤差進(jìn)行比較，最后1個(gè)圖是采用STK中的HPOP對(duì)空間站軌道協(xié)方差的演化結(jié)果。

相比較發(fā)現(xiàn)，線性方法在S方向一天的協(xié)方差演化誤差是米量級(jí)，UKF方法演化誤差是分米量級(jí)。UKF協(xié)方差演化誤差精度高于線性方法1個(gè)量級(jí)。

采用同樣方法對(duì)Ajisai衛(wèi)星軌道協(xié)方差演化，結(jié)果如圖2所示，相同的初值協(xié)方差A(yù)jisai衛(wèi)星通過(guò)HPOP方法協(xié)方差演化一天的值小于空間站的演化值。這主要是空間站所受的大氣阻力大于Ajisai衛(wèi)星。對(duì)于Ajisai衛(wèi)星S方向演化一天的誤差，線性方法是米量級(jí)，UKF方法演化誤差是厘米量級(jí)。UKF誤差精度高于線性方法兩個(gè)量級(jí)。對(duì)于低軌道空間目標(biāo)，軌道高度越高，UKF協(xié)方差演化誤差精度越高。

從圖1和圖2可以看出，對(duì)于空間站和Ajisai衛(wèi)星，通過(guò)UKF方法演化的R方向協(xié)方差誤差大于線性方法演化誤差。這是由于軌道協(xié)方差隨時(shí)間演化呈現(xiàn)非高斯性引起的。采用蒙特卡洛方法可以驗(yàn)證當(dāng)軌道協(xié)方差演化時(shí)長(zhǎng)為一天時(shí)，空間目標(biāo)的軌道協(xié)方差不再是高斯分布的橢球，而是呈月牙狀。月牙的彎曲方向就是在徑向(R方向)。圖3是采用蒙特卡洛方法演化的軌道協(xié)方差分布圖。對(duì)Ajisai衛(wèi)星初始根數(shù)加上位置標(biāo)準(zhǔn)差10 m、速度標(biāo)準(zhǔn)差0.01 m的隨機(jī)誤差，構(gòu)成10 000個(gè)采樣點(diǎn)，含隨機(jī)誤差的采樣點(diǎn)在與HPOP相同模型攝動(dòng)力作用下向前傳播一天，把演化后的10 000個(gè)采樣點(diǎn)與標(biāo)準(zhǔn)精密軌道相差得到蒙特卡洛誤差采樣點(diǎn)分布。從圖3可以看出蒙特卡洛采樣點(diǎn)分布不再是標(biāo)準(zhǔn)橢球形，采樣點(diǎn)分布在R方向呈現(xiàn)彎曲的非高斯性分布。這種非高斯性彎曲分布導(dǎo)致R方向誤差協(xié)方差值變大，而UKF方法可以比較明顯地表示這種協(xié)方差演化特性。

2.2 UCTs關(guān)聯(lián)模擬

選取2017年2月印度發(fā)射的一箭104星中的20個(gè)目標(biāo)衛(wèi)星做UCTs關(guān)聯(lián)，衛(wèi)星軌道半長(zhǎng)軸6 860～6 890 km，軌道傾角97.5°，衛(wèi)星軌道區(qū)域分布比較集中。各顆衛(wèi)星之間具有相似的軌道類型，在與箭體分離后一段時(shí)間各個(gè)衛(wèi)星在同一軌道面上，軌道之間只存在較小的相位差異，采用之前半長(zhǎng)軸、軌道傾角篩選再進(jìn)行空間目標(biāo)編目存在困難，適用于驗(yàn)證基于協(xié)方差算法軌道關(guān)聯(lián)。

每顆衛(wèi)星選取2～4個(gè)不同時(shí)間點(diǎn)的軌道星歷，軌道間隔小于兩天。采用美國(guó)巴利亞多(Vallado)提供的雷達(dá)參數(shù)模擬20個(gè)目標(biāo)衛(wèi)星的雷達(dá)觀測(cè)數(shù)據(jù)。并用HerrickGibbs方法對(duì)短弧段各個(gè)衛(wèi)星初始軌道進(jìn)行定軌，用最小二乘法對(duì)軌道進(jìn)行精密定軌，生成衛(wèi)星在某點(diǎn)的軌道根數(shù)及協(xié)方差信息。應(yīng)用軌道預(yù)報(bào)方法對(duì)軌道根數(shù)進(jìn)行數(shù)值預(yù)報(bào)，通過(guò)線性協(xié)方差演化方法和UKF協(xié)方差演化方法對(duì)衛(wèi)星初始協(xié)方差進(jìn)行演化。

20個(gè)UCTs模擬真實(shí)軌道是基于space-track提供的兩行軌道根數(shù)，由于公布的兩行軌道根數(shù)不包含軌道協(xié)方差信息，因此假設(shè)其為真實(shí)軌道重新定軌生成協(xié)方差。而兩行軌道根數(shù)本身含有軌道誤差，這可能是關(guān)聯(lián)過(guò)程中關(guān)聯(lián)閾值較大的主要原因。同時(shí)兩行軌道根數(shù)是以真赤道平春分點(diǎn)(True Equator Mean Equinox,TEME)作為參考坐標(biāo)系，需要轉(zhuǎn)變?yōu)镴2000坐標(biāo)系分析。

圖4是20個(gè)衛(wèi)星中，相同衛(wèi)星不同歷元軌道根數(shù)進(jìn)行協(xié)方差關(guān)聯(lián)的關(guān)聯(lián)量k值的統(tǒng)計(jì)分布情況。可以看出20顆衛(wèi)星UCTs關(guān)聯(lián)量分布主要在60～80。分析發(fā)現(xiàn)當(dāng)UCTs關(guān)聯(lián)閾值限定為200時(shí)，20顆衛(wèi)星中相同衛(wèi)星不同歷元的軌道根數(shù)都可以成功關(guān)聯(lián)。關(guān)聯(lián)值大的原因作者分析如下：

1)首先文章中選用美國(guó)發(fā)布的兩行軌道根數(shù)，這種根數(shù)本身存在較大的軌道誤差。由于模擬中選取20顆衛(wèi)星的兩行軌道根數(shù)作為真實(shí)軌道，忽略了兩行軌道根數(shù)數(shù)據(jù)本身所含有的誤差，所以計(jì)算的關(guān)聯(lián)閾值較大。

2)在雷達(dá)觀測(cè)數(shù)據(jù)模擬過(guò)程中根據(jù)不同雷達(dá)測(cè)站的噪聲加入了測(cè)量誤差。

3)在精密定軌過(guò)程中由于弧段較短時(shí)不能進(jìn)行精密定軌確定協(xié)方差，因此在一些衛(wèi)星之間的關(guān)聯(lián)模擬過(guò)程中選用了較長(zhǎng)的弧段進(jìn)行定軌，弧段長(zhǎng)定軌協(xié)方差小，而馬氏距離計(jì)算時(shí)協(xié)方差的逆較大，關(guān)聯(lián)值較大。

假設(shè)選取精密的軌道根數(shù)作為真實(shí)軌道，短弧段情況下可以計(jì)算出軌道協(xié)方差，則關(guān)聯(lián)閾值會(huì)大幅下降。

圖5表示相同衛(wèi)星不同歷元間隔對(duì)應(yīng)的關(guān)聯(lián)量?？梢钥闯觯S著軌道時(shí)間間隔增加關(guān)聯(lián)量逐漸下降，相同軌道時(shí)間間隔，不同衛(wèi)星軌道所計(jì)算的關(guān)聯(lián)量也不相同。關(guān)于關(guān)聯(lián)閾值的選取需要針對(duì)不同衛(wèi)星軌道做進(jìn)一步的分析。

3 結(jié)束語(yǔ)

本文研究了基于協(xié)方差理論的UCTs動(dòng)態(tài)關(guān)聯(lián)算法，對(duì)比分析了基于雅克比轉(zhuǎn)換的線性協(xié)方差演化方法和基于UT轉(zhuǎn)換的UKF協(xié)方差演化方法。以400 km和1k500kkm軌道高度為例，初始軌道協(xié)方差相同的情況下，基于UT變換的協(xié)方差演化精度要高于基于雅克比轉(zhuǎn)換的線性協(xié)方差方法演化精度。應(yīng)用UKF方法軌道協(xié)方差演化一天的位置誤差分別是分米量級(jí)和厘米量級(jí)，應(yīng)用線性方法演化一天的位置誤差是米級(jí)。

把協(xié)方差演化方法應(yīng)用到UCTs動(dòng)態(tài)關(guān)聯(lián)中去，模擬印度一箭104星中的20顆衛(wèi)星的軌道關(guān)聯(lián)，當(dāng)軌道關(guān)聯(lián)量閾值選取為200時(shí)，20顆衛(wèi)星軌道可以全部關(guān)聯(lián)。

軌道協(xié)方差演化有很好的應(yīng)用前景，本文在空間目標(biāo)編目方面，為處理各非關(guān)聯(lián)軌道提出了一種動(dòng)態(tài)算法，其關(guān)聯(lián)量源于每個(gè)非關(guān)聯(lián)軌道的軌道誤差協(xié)方差矩陣。另外空間目標(biāo)碰撞評(píng)估的一個(gè)重要環(huán)節(jié)是計(jì)算兩個(gè)空間物體的碰撞概率。碰撞概率的計(jì)算使用了相對(duì)狀態(tài)概率密度函數(shù)的數(shù)值積分算法。每種算法需要為兩個(gè)對(duì)象輸入狀態(tài)向量的誤差信息。狀態(tài)向量的誤差信息通過(guò)協(xié)方差矩陣的方式來(lái)表示。本研究可以為將來(lái)的空間目標(biāo)編目和空間碎片預(yù)警提供方法參考。