龔柏春, 金 鑫,,*, 張仁勇, 陳修橋

(1. 南京航空航天大學(xué)航天學(xué)院, 江蘇 南京 210016; 2. 中國(guó)科學(xué)院空間應(yīng)用工程與 技術(shù)中心, 北京 100038; 3. 中國(guó)人民解放軍32032部隊(duì), 北京 100094)

0 引 言

近年來,航天技術(shù)快速發(fā)展,各國(guó)的航天活動(dòng)日益頻繁,近地軌道上衛(wèi)星數(shù)量急劇增長(zhǎng),擁擠的太空環(huán)境給衛(wèi)星的安全帶來了威脅[1-4]。根據(jù)空間安全威脅的來源和性質(zhì),可將其分為兩類:第一類源于自然環(huán)境因素,如小行星碎片撞擊;第二類則來自人為威懾和威脅,如衛(wèi)星武器[1-9]。因此,實(shí)現(xiàn)對(duì)空間非合作目標(biāo)的態(tài)勢(shì)感知與跟蹤，對(duì)保障我國(guó)現(xiàn)役衛(wèi)星的安全而言具有重要意義。

基于星載測(cè)量傳感器實(shí)現(xiàn)對(duì)非合作目標(biāo)的軌道機(jī)動(dòng)檢測(cè)已經(jīng)成為了天基態(tài)勢(shì)感知與跟蹤領(lǐng)域需要突破的關(guān)鍵技術(shù),已經(jīng)得到了國(guó)內(nèi)外的廣泛關(guān)注和研究[9-12]。根據(jù)技術(shù)實(shí)現(xiàn)途徑,現(xiàn)有的空間非合作目標(biāo)軌道機(jī)動(dòng)檢測(cè)方法主要可分為兩類,現(xiàn)綜述如下。

第一類研究側(cè)重于機(jī)動(dòng)特征量的選取[13-17],不對(duì)測(cè)量方法有過高或特定的要求,如Qin等[15]根據(jù)衛(wèi)星軌道半長(zhǎng)軸的長(zhǎng)期變化,提出了地球同步軌道、傾斜地球同步軌道和中地球軌道衛(wèi)星面內(nèi)軌道機(jī)動(dòng)的檢測(cè)方法，以及地球同步軌道衛(wèi)星面外軌道機(jī)動(dòng)的檢測(cè)方法。Huang等[16]使用半長(zhǎng)軸和偏心率作為檢測(cè)量來判斷目標(biāo)是否機(jī)動(dòng),但其方法僅適用于檢測(cè)在近地點(diǎn)或遠(yuǎn)地點(diǎn)沿速度方向的脈沖機(jī)動(dòng)。Roberts等[17]提出了使用公開的地球同步經(jīng)度數(shù)據(jù)庫進(jìn)行神經(jīng)網(wǎng)絡(luò)訓(xùn)練來預(yù)測(cè)衛(wèi)星狀態(tài),以星下點(diǎn)的經(jīng)緯度和衛(wèi)星高度作為特征量進(jìn)行檢測(cè),但該方法僅適用于地球同步軌道,具有較高的局限性。

第二類研究側(cè)重于對(duì)衛(wèi)星狀態(tài)參數(shù)測(cè)量及數(shù)據(jù)處理方法進(jìn)行探討[18-25],以此來檢測(cè)采用傳統(tǒng)測(cè)量方法因精度不夠而無法識(shí)別的軌道機(jī)動(dòng),或者直接通過特定測(cè)量/計(jì)算量的代價(jià)忽略機(jī)動(dòng)存在以實(shí)現(xiàn)穩(wěn)定跟蹤。王慶瑞等[23]將距離變化率設(shè)定為檢驗(yàn)參數(shù),提出了基于概率判決的軌道機(jī)動(dòng)檢驗(yàn)方法,基于Neyman-Pearson準(zhǔn)則，能夠生成機(jī)動(dòng)的判決門限，從而實(shí)現(xiàn)自適應(yīng)的檢驗(yàn)軌道機(jī)動(dòng),最終可以實(shí)現(xiàn)一般情況下正確率大于0.8的機(jī)動(dòng)檢測(cè)。Jia等[24]提出使用4顆相互協(xié)作的衛(wèi)星來解決卡爾曼濾波法在目標(biāo)機(jī)動(dòng)后檢測(cè)性能降低的問題,但該方法涉及到測(cè)量衛(wèi)星相互協(xié)作的問題,過程復(fù)雜,且需要考慮到某一個(gè)或多個(gè)測(cè)量衛(wèi)星無法觀測(cè)時(shí)對(duì)系統(tǒng)的調(diào)整問題。Liu等[25]提出了聯(lián)合濾波的方法,其利用主、輔助濾波器的配合消除了機(jī)動(dòng)對(duì)目標(biāo)跟蹤的影響?？梢?第二類方法計(jì)算量較大,給衛(wèi)星增加了負(fù)擔(dān),因此如何用普通的跟蹤濾波器對(duì)空間機(jī)動(dòng)目標(biāo)進(jìn)行跟蹤,是一個(gè)難點(diǎn)[26-28]。

為了解決上述問題,本文針對(duì)空間非合作目標(biāo)近程態(tài)勢(shì)感知任務(wù),在構(gòu)建一種新型機(jī)動(dòng)檢測(cè)特征量的基礎(chǔ)上,提出利用無源光學(xué)相機(jī)僅測(cè)量角度信息的機(jī)動(dòng)檢測(cè)和跟蹤方法。下文將在第1節(jié)和第2節(jié)中建立相對(duì)動(dòng)力學(xué)模型和傳感器測(cè)量模型,在第3節(jié)中設(shè)計(jì)基于新型特征量的機(jī)動(dòng)檢測(cè)算法,在第4節(jié)中進(jìn)行數(shù)值仿真實(shí)驗(yàn)與性能分析。

1 相對(duì)運(yùn)動(dòng)動(dòng)力學(xué)模型

如圖1所示,建立當(dāng)?shù)卮怪碑?dāng)?shù)厮?local vertical local horizontal, LVLH)參考坐標(biāo)系,原點(diǎn)在感知衛(wèi)星的質(zhì)心,X軸由地心指向衛(wèi)星質(zhì)心,Z軸垂直于軌道,Y軸滿足右手定則。

圖1 地心坐標(biāo)系(左)和LVLH系(右)Fig.1 Geocentric coordinate system (left) and LVLH frame (right)

在近圓軌道、兩個(gè)航天器之間的距離足夠小且不考慮額外攝動(dòng)的情況下,可以采用線性化的Clohessy-Wiltshire方程來建模航天器之間的相對(duì)運(yùn)動(dòng)[29-30]。ti時(shí)刻的相對(duì)軌道狀態(tài)xi可以用to時(shí)刻的狀態(tài)向量x0表示:

xi=Φix0

(1)

式中:Φi為狀態(tài)轉(zhuǎn)移矩陣,具體形式參見文獻(xiàn)[29]。

2 傳感器測(cè)量模型

本文假設(shè)感知衛(wèi)星自身的軌道和姿態(tài)已知,星上僅有無源光學(xué)相機(jī)可以用來對(duì)非合作目標(biāo)進(jìn)行相對(duì)測(cè)量。根據(jù)工程實(shí)際約束情況,相機(jī)偏離衛(wèi)星質(zhì)心安裝,偏心測(cè)量如圖2所示。因此,可以直接在LVLH參考坐標(biāo)系下建立俯仰角和偏航角的測(cè)量模型:

(2)

式中:[x,y,z]T是目標(biāo)ti時(shí)刻的相對(duì)位置;[dx,dy,dz]T為相機(jī)的偏置量;[n1,n2]T為測(cè)量誤差,設(shè)其為零均值高斯白噪聲。

圖2 相機(jī)偏置測(cè)量示意圖Fig.2 Illustration diagram of camera offset measurement

3 基于相對(duì)角動(dòng)量的機(jī)動(dòng)檢測(cè)算法

本文所研究的非合作目標(biāo)軌道機(jī)動(dòng)僅測(cè)角檢測(cè)方法的基本邏輯關(guān)系如圖3所示。首先以式(2)所示的視線角測(cè)量量結(jié)合式(1)所示的動(dòng)力學(xué)模型進(jìn)行濾波估計(jì),獲得非合作目標(biāo)的相對(duì)軌道參數(shù),然后基于相對(duì)軌道參數(shù)進(jìn)行軌道機(jī)動(dòng)待檢量的構(gòu)造,最后通過對(duì)待檢量進(jìn)行檢驗(yàn)判別,來確定機(jī)動(dòng)時(shí)間,并下達(dá)濾波系統(tǒng)重新初始化指令。

圖3 算法邏輯框圖Fig.3 Algorithm logic block diagram

由于本文的重點(diǎn)在于機(jī)動(dòng)檢測(cè)算法設(shè)計(jì),所以對(duì)圖3中的僅測(cè)角定軌算法不過多論述,而是直接采用Gong等在文獻(xiàn)[31-32]中提出的方法。視線角與相對(duì)軌道參量的映射關(guān)系如下:

[rT,vT]T=f(az,el)

(3)

下文中的檢測(cè)算法設(shè)計(jì)與分析過程將直接采用r和v進(jìn)行。下面將重點(diǎn)先闡述基于相對(duì)角動(dòng)量的新型機(jī)動(dòng)待檢特征量的構(gòu)建及其敏感性分析,然后給出采用該特征量的軌道機(jī)動(dòng)自動(dòng)檢測(cè)算法。

3.1 機(jī)動(dòng)待檢特征量建模及其敏感性分析

傳統(tǒng)的軌道機(jī)動(dòng)檢測(cè)特征量存在敏感性差、適用性窄等問題,難以滿足空間對(duì)抗對(duì)態(tài)勢(shì)感知靈敏性、檢測(cè)快速性的要求。因此,本文提出一種基于相對(duì)角動(dòng)量概念的新型機(jī)動(dòng)檢測(cè)特征量。下面首先闡述該待檢特征量的定義,然后對(duì)其敏感性進(jìn)行分析、證明。

定義相對(duì)角動(dòng)量為

hrel=rrel×vrel

(4)

式中:

(5)

衛(wèi)星的機(jī)動(dòng)會(huì)直接影響衛(wèi)星的位置和速度矢量,而角動(dòng)量的計(jì)算與物體的位置和速度矢量都具有相關(guān)性,其通過將速度和位置矢量進(jìn)行叉乘的處理方式,不僅可以反映位置和速度矢量大小的變化,也可以直接通過角動(dòng)量數(shù)值的變化反映位置和速度矢量方向上的變化[33],因此相對(duì)角動(dòng)量對(duì)衛(wèi)星機(jī)動(dòng)更加敏感。此外,相對(duì)角動(dòng)量參數(shù)對(duì)衛(wèi)星所處軌道以及軌道上的位置沒有要求,所以也將更具普適性。

同時(shí),與地心慣性系下的軌道角動(dòng)量相比,由式(4)定義的軌道系下的相對(duì)角動(dòng)量對(duì)目標(biāo)機(jī)動(dòng)更具敏感性,現(xiàn)將原因分析如下。衛(wèi)星攜帶能量有限,所以衛(wèi)星機(jī)動(dòng)的沖量較小[34-35],這使其機(jī)動(dòng)后的短時(shí)間內(nèi)的速度變化量也很小,而衛(wèi)星在低地球軌道運(yùn)行時(shí)速度要達(dá)到7.5 km/s量級(jí)。因此,衛(wèi)星在機(jī)動(dòng)后短時(shí)間內(nèi)的慣性位置、速度改變量相較于其機(jī)動(dòng)前的位置大小、速度大小是極小的。而以地心為參考點(diǎn)的角動(dòng)量由相對(duì)于地心的位置和速度貢獻(xiàn),這使慣性系下的角動(dòng)量變化也很難在噪聲中被發(fā)現(xiàn)。然而,若取測(cè)量衛(wèi)星為參考點(diǎn),軌道系下的相對(duì)角動(dòng)量由相對(duì)位置與相對(duì)速度貢獻(xiàn),大大減小了變化量與原始基礎(chǔ)量之間的差距,因此可以有效減小相對(duì)角動(dòng)量的變化量被噪聲淹沒的風(fēng)險(xiǎn)。

下面對(duì)上述定性分析進(jìn)行數(shù)學(xué)證明。根據(jù)式(4)可以得到機(jī)動(dòng)前t1時(shí)刻和機(jī)動(dòng)完成后t2時(shí)刻的相對(duì)角動(dòng)量分別為

(6)

當(dāng)機(jī)動(dòng)為脈沖機(jī)動(dòng)時(shí),由于機(jī)動(dòng)時(shí)間很短,即tm≈0, 其中tm為機(jī)動(dòng)時(shí)間。

因此,

(7)

則

hrel2=rrel2×vrel2=rrel1×(vrel1+Δvrel)=

rrel1×vrel1+rrel1×Δvrel

(8)

記

Δhrel=rrel1×Δvrel

(9)

與式(6)～式(9)的過程同理,可以得到地心慣性系下的角動(dòng)量為

(10)

由于常見的單個(gè)脈沖所能提供的加速度大小約為0.1～1 m/s2[34-35],而衛(wèi)星相對(duì)于地心的運(yùn)行速度至少達(dá)到7.5 km/s量級(jí),所以有

(11)

從而有

he1?Δhe

(12)

式中:h=norm(h)。

因此,以地心慣性系下的角動(dòng)量為檢驗(yàn)標(biāo)準(zhǔn),將很難在短時(shí)間內(nèi)從噪聲中辨別出是否機(jī)動(dòng)。

相反地,

(13)

從而有

(14)

由此可以看出,相對(duì)角動(dòng)量作為檢驗(yàn)標(biāo)準(zhǔn),將具有更優(yōu)良的特性。

當(dāng)考慮相對(duì)軌道估計(jì)誤差時(shí),式(8)變?yōu)?/p>

hrel2=(rrel2+nr)×(vrel2+nv)=

(rrel1+nr)×(vrel1+Δvrel+nv)=

rrel1×vrel1+rrel1×Δvrel+rrel1×nv+

nr×vrel1+nr×Δvrel+nr×nv

(15)

式中:nr和nv分別為相對(duì)位置和速度誤差。

(16)

結(jié)合式(15)和式(16)可得

‖rrel1×Δvrel+nr×Δvrel‖?

‖rrel1×nv+nr×vrel1+nr×nv‖

(17)

式(17)表明了由軌道機(jī)動(dòng)造成的相對(duì)角動(dòng)量的變化將顯著地區(qū)別于由相對(duì)軌道估計(jì)誤差引起的變化。因此,只要相對(duì)軌道估計(jì)誤差滿足式(17)所示的條件,相對(duì)角動(dòng)量對(duì)目標(biāo)機(jī)動(dòng)是非常敏感的。

3.2 軌道機(jī)動(dòng)點(diǎn)自動(dòng)判別算法

若目標(biāo)衛(wèi)星所受到的合力并非指向感知衛(wèi)星,則其相對(duì)角動(dòng)量會(huì)發(fā)生變化,變化規(guī)律為

(18)

式中:τ為目標(biāo)衛(wèi)星B受到的相對(duì)于感知衛(wèi)星A質(zhì)心的單位質(zhì)量力矩。

在目標(biāo)衛(wèi)星B沒有施加機(jī)動(dòng)時(shí),τ僅為地心引力以及相關(guān)攝動(dòng)力提供,因而τ的變化很平緩,即hrel的變化是連續(xù)的。換言之,如果hrel出現(xiàn)了突變的情況,那么一定是目標(biāo)衛(wèi)星B進(jìn)行了機(jī)動(dòng)(或被碰撞)。因此,檢測(cè)hrel的突變就能判別機(jī)動(dòng)。

本文設(shè)計(jì)了基于二次滑窗方差比值的自主尋找機(jī)動(dòng)點(diǎn)方法,算法流程如圖4所示。

圖4 機(jī)動(dòng)點(diǎn)檢測(cè)算法框圖Fig.4 Framework diagram of the maneuver detection algorithm

在每完成一次滑窗內(nèi)的計(jì)算后,滑窗向后移動(dòng)一個(gè)格,重新計(jì)算并記錄結(jié)果。將通過對(duì)相對(duì)角動(dòng)量進(jìn)行兩次方差計(jì)算后的數(shù)據(jù)記作Dh。

該襯板需要專業(yè)工廠定做，建議在對(duì)抗沖擊性能要求不大，對(duì)耐磨性能要求嚴(yán)格，而且形狀較簡(jiǎn)單的部位鋪設(shè)該襯板。但該襯板存在如下問題：① 表面不光滑，初始狀態(tài)摩擦阻力較大；② 襯板直接接口不規(guī)則，沉頭螺栓部位阻力大；③ 由于硬度太高，加工性能較困難。

通過計(jì)算方差的方式將相對(duì)角動(dòng)量突變轉(zhuǎn)化為數(shù)值大小上的突變。為了使數(shù)值突變更加明顯,本文采用二次滑窗的方式來處理相對(duì)角動(dòng)量。設(shè)置一個(gè)寬度為2N+1、步長(zhǎng)為1的數(shù)據(jù)滑窗(包括待檢測(cè)點(diǎn)ti時(shí)刻,前后再各取N個(gè)數(shù)據(jù)),N的選取與定軌精度相關(guān)。滑窗內(nèi)的方差為

(19)

雖然可以從Dj序列的突變來檢測(cè)機(jī)動(dòng),但由于目標(biāo)星的相對(duì)位置和機(jī)動(dòng)方式不同,Dj的范圍差異很大,無法設(shè)置固定的閥值來判斷是否機(jī)動(dòng)。為了解決這個(gè)問題,首先在濾波平穩(wěn)后,計(jì)算該點(diǎn)到檢測(cè)點(diǎn)Dhj的平均值:

(20)

式中:Dhi表示檢測(cè)點(diǎn)時(shí)刻tj的Dh。

然后,為了消除各種情況間的差異,計(jì)算各個(gè)時(shí)刻Dhj與Dhmean(j)的比值:

Ni=Dhj/Dhmean(j)

(21)

同時(shí),為了縮小未機(jī)動(dòng)時(shí)由濾波造成的浮動(dòng),并凸顯機(jī)動(dòng)時(shí)刻N(yùn)值的變化,可以將Dj的首位D0替換為相對(duì)于其的無窮大,如10 000D0。為了進(jìn)一步擴(kuò)大機(jī)動(dòng)后的N值與未機(jī)動(dòng)時(shí)的N值的差異,可以將所有的N值乘以同一倍數(shù),如200，則有

Ni,X/Y/Z≥Nf

(22)

式中:下標(biāo)X/Y/Z表示Ni由LVLH參考坐標(biāo)系下X、Y或Z方向計(jì)算得到,Nf為由實(shí)驗(yàn)得到的閥值。

若滿足式(22),則說明滑窗內(nèi)包含機(jī)動(dòng)時(shí)刻數(shù)據(jù),認(rèn)為t=i即為機(jī)動(dòng)時(shí)刻。

4 數(shù)值仿真驗(yàn)證

假設(shè)目標(biāo)星和感知衛(wèi)星都在近圓軌道上飛行,二者之間的初始相對(duì)軌道如表1所示。測(cè)角相機(jī)的安裝矢量為[1,1,3]Tm,安裝誤差為1%,測(cè)角誤差為零均值高斯白噪聲,其均方差將在仿真算例中具體說明,感知衛(wèi)星的姿態(tài)測(cè)量包含0.001 rad/axis的常值偏差和均方差為0.001 rad的白噪聲。參考軌道的積分步長(zhǎng)為1 s。

表1 目標(biāo)衛(wèi)星相對(duì)于測(cè)量衛(wèi)星的初始相對(duì)軌道參數(shù)

為了驗(yàn)證該方法的有效性,首先設(shè)置目標(biāo)衛(wèi)星在t=500 s時(shí)刻分別沿LVLH參考坐標(biāo)系下X軸、Y軸和Z軸添加ΔV=0.1 m/s的脈沖。之后,分別在t=500 s時(shí)刻和t=1 000 s時(shí)刻設(shè)置不同的脈沖機(jī)動(dòng)。

下面針對(duì)相對(duì)角動(dòng)量對(duì)軌道機(jī)動(dòng)的敏感性、敏感器精度對(duì)檢測(cè)的影響、機(jī)動(dòng)檢測(cè)相對(duì)誤差和雙脈沖機(jī)動(dòng)下的檢測(cè)結(jié)果與持續(xù)跟蹤觀測(cè)效果進(jìn)行仿真實(shí)驗(yàn)驗(yàn)證和性能分析。

4.1 相對(duì)角動(dòng)量對(duì)機(jī)動(dòng)反應(yīng)的仿真結(jié)果

為了驗(yàn)證本文所提出的相對(duì)角動(dòng)量對(duì)目標(biāo)衛(wèi)星機(jī)動(dòng)的敏感性,設(shè)置相機(jī)的測(cè)角均方差為3×10-4rad,并以兩種不同相對(duì)軌跡為例進(jìn)行驗(yàn)證。圖5、圖6分別為V-Bar靜止軌道添加Z方向脈沖、跳躍軌道添加X方向脈沖后相對(duì)角動(dòng)量的變化。脈沖添加點(diǎn)均在t=500 s時(shí)刻,脈沖大小為0.1 m/s。

圖5 V-bar靜止軌道添加Z方向脈沖后的相對(duì)角動(dòng)量變化Fig.5 Relative angular momentum change on V-bar stationary orbit with the pulse added on Z-axis

圖6 跳躍軌道添加X方向脈沖后的相對(duì)角動(dòng)量的變化Fig.6 Relative angular momentum change on hop orbit with the pulse on X-axis

在圖5、圖6中,前期為濾波收斂過程造成的數(shù)值波動(dòng),在濾波收斂之后到脈沖機(jī)動(dòng)之前,各個(gè)軸上的相對(duì)角動(dòng)量基本穩(wěn)定。但是非常明顯地，在機(jī)動(dòng)脈沖施加點(diǎn)t=500 s之后,圖5中X軸、Y軸和圖6中Y軸的相對(duì)角動(dòng)量發(fā)生了突變。這就表明了相對(duì)角動(dòng)量對(duì)軌道機(jī)動(dòng)具有優(yōu)良的敏感能力。

4.2 測(cè)角精度對(duì)脈沖機(jī)動(dòng)檢測(cè)的影響結(jié)果

該仿真用于驗(yàn)證僅測(cè)角導(dǎo)航情況下的機(jī)動(dòng)檢測(cè)有效性。首先,在測(cè)角均方差為3×10-4rad的情況下得到合適的閥值Nf=Nf_-4;接下來,在取閥值為Nf_-4的情況下,進(jìn)行敏感器測(cè)角均方差分別為3×10-5rad、3×10-4rad、3×10-3rad,脈沖機(jī)動(dòng)大小分別為ΔV=0.095～0.100 m/s和ΔV=0.950～1.000 m/s,初始相對(duì)軌道隨機(jī)、機(jī)動(dòng)方向隨機(jī)的300次打靶仿真,統(tǒng)計(jì)結(jié)果如表2所示。由表2可知,在閥值均取Nf_-4的情況下,測(cè)角均方差為3×10-4rad的機(jī)動(dòng)檢測(cè)效果最好,測(cè)角均方差為3×10-5rad時(shí)效果稍微減弱,測(cè)角均方差為3×10-3rad時(shí)小機(jī)動(dòng)的檢測(cè)效果明顯下降,這是因?yàn)樵谑褂貌煌鹊南鄼C(jī)時(shí),需要對(duì)應(yīng)地調(diào)整Nf以達(dá)到最優(yōu),但本文為了說明所提方法具有一定的魯棒性,采用了統(tǒng)一的Nf_-4進(jìn)行不同精度相機(jī)下的仿真。

從表2中可以看出,采用針對(duì)3×10-4rad測(cè)角精度設(shè)置的最優(yōu)經(jīng)驗(yàn)閥值時(shí),1 m/s量級(jí)的機(jī)動(dòng)檢測(cè)誤差小于60 s的概率大于95%。這相比于文獻(xiàn)[23]提出的基于距離變化率的概率判別方法在相近條件下的80%成功率在相近測(cè)量精度水平下得到了提高;與文獻(xiàn)[36]所提方法的92%左右的成功率相比,本文方法也有一定程度的提高,因此本文所提出的檢測(cè)方法在相關(guān)的應(yīng)用場(chǎng)合具有較強(qiáng)的優(yōu)勢(shì)。

表2 機(jī)動(dòng)檢測(cè)成功率

4.3 單脈沖機(jī)動(dòng)下的機(jī)動(dòng)檢測(cè)誤差仿真結(jié)果

本文在t=700 s添加常見的大小范圍為0.1～1 m/s的脈沖機(jī)動(dòng)，對(duì)各種情況進(jìn)行了驗(yàn)證。設(shè)置測(cè)角均方差為3×10-4rad,對(duì)每種情況進(jìn)行了重復(fù)20次的計(jì)算,并將ΔV分別為0.1 m/s和1 m/s時(shí)的平均誤差在表3中給出,其中不包含漏檢的情況,包含所有機(jī)動(dòng)被檢測(cè)到的情況。由表3可知,本文所提出的軌道機(jī)動(dòng)檢測(cè)方法的相對(duì)檢測(cè)誤差平均值在60 s以內(nèi)。

4.4 雙脈沖機(jī)動(dòng)下的機(jī)動(dòng)檢測(cè)與跟蹤結(jié)果

為了驗(yàn)證本文方法對(duì)目標(biāo)的多次機(jī)動(dòng)跟蹤能力,這里以橢圓繞飛為例進(jìn)行仿真驗(yàn)證。設(shè)置測(cè)角均方差為3×10-4rad,假定目標(biāo)機(jī)動(dòng)兩次,第一次在第500 s施加Z方向0.1 m/s的機(jī)動(dòng)脈沖,第二次在1 000 s施加X方向1 m/s的機(jī)動(dòng)脈沖。

仿真結(jié)果如圖7和圖8所示,算法給出的第一次機(jī)動(dòng)時(shí)間發(fā)生在第459 s,誤差為41 s,第二次機(jī)動(dòng)時(shí)間發(fā)生在第951 s,誤差為49 s?？梢娫趦煞N情況下,目標(biāo)星的機(jī)動(dòng)都能被檢測(cè)出來,且確定的機(jī)動(dòng)時(shí)間誤差較小。在檢測(cè)到目標(biāo)機(jī)動(dòng)并重啟相對(duì)軌道估計(jì)濾波器之后,估計(jì)誤差均趨于穩(wěn)定,實(shí)現(xiàn)了持續(xù)跟蹤定軌。相反地,濾波器無重啟情況下的估計(jì)誤差是發(fā)散的。

表3 機(jī)動(dòng)檢測(cè)結(jié)果

圖7 相對(duì)位置估計(jì)誤差對(duì)比曲線Fig.7 Comparison curve of relative position estimation error

圖8 相對(duì)速度估計(jì)誤差對(duì)比曲線Fig.8 Comparison curve of relative velocity estimation error

5 結(jié)束語

本文針對(duì)近圓軌道非合作目標(biāo)近程感知任務(wù),提出了一種基于新型特征量的目標(biāo)軌道機(jī)動(dòng)僅測(cè)角檢測(cè)及跟蹤方法,進(jìn)行了新型特征量的敏感性分析與量化證明,建立了整個(gè)算法流程框架,并進(jìn)行了數(shù)值仿真驗(yàn)證與性能分析。主要結(jié)論如下:

基于本文所提算法能夠在近程感知任務(wù)中實(shí)現(xiàn)對(duì)目標(biāo)軌道機(jī)動(dòng)的快速檢測(cè)。當(dāng)采用測(cè)角精度不低于3×10-4rad的無源相機(jī)進(jìn)行測(cè)量時(shí), 0.1 m/s～1 m/s量級(jí)范圍內(nèi)的小脈沖機(jī)動(dòng)可被有效檢測(cè),機(jī)動(dòng)時(shí)間檢測(cè)誤差在60 s以內(nèi)的概率超過95%。判定目標(biāo)機(jī)動(dòng)之后,通過重新初始化僅測(cè)角相對(duì)導(dǎo)航濾波器能夠?qū)崿F(xiàn)對(duì)目標(biāo)的持續(xù)、穩(wěn)定跟蹤。本文所提算法比現(xiàn)有算法的檢測(cè)成功率更高,且算法更為簡(jiǎn)便、適用性更寬。

對(duì)于無先驗(yàn)信息情況下的橢圓軌道非合作目標(biāo)機(jī)動(dòng)遠(yuǎn)程感知檢測(cè),將是下一步的研究?jī)?nèi)容。