(中國電子科技集團(tuán)公司第二十研究所,陜西西安710068)

0 引言

隨著艦載防空導(dǎo)彈技術(shù)的發(fā)展,現(xiàn)代反艦導(dǎo)彈在距離目標(biāo)艦艇30 km以內(nèi)時(shí),普遍采用機(jī)動(dòng)的方式降低艦載防空系統(tǒng)的跟蹤精度與命中概率[1]。反艦導(dǎo)彈在中段距離目標(biāo)20 km左右,一般采用蛇形機(jī)動(dòng)等方式規(guī)避艦空導(dǎo)彈,有的反艦導(dǎo)彈除了蛇形機(jī)動(dòng)外還采用降高或在極短時(shí)間內(nèi)將速度由亞音速迅速加速到2.5～3 Ma的超音速的方式來規(guī)避艦載防空導(dǎo)彈的打擊;在末段10 km以內(nèi),一般采用爬升俯沖或繼續(xù)降高飛行的方式規(guī)避防空武器的打擊,這對(duì)現(xiàn)代艦載防空武器與探測(cè)跟蹤系統(tǒng)提出了較高的要求。為了應(yīng)對(duì)反艦導(dǎo)彈的中端機(jī)動(dòng)等問題,需要對(duì)機(jī)動(dòng)目標(biāo)跟蹤算法改進(jìn)。典型的目標(biāo)跟蹤濾波器都是單模型的,因此只能對(duì)具有單一運(yùn)動(dòng)狀態(tài)的目標(biāo)有較好的跟蹤性能。但是實(shí)際的目標(biāo)運(yùn)動(dòng)經(jīng)常由多種運(yùn)動(dòng)狀態(tài)組合而成,因此單一模型的濾波器對(duì)機(jī)動(dòng)目標(biāo)的跟蹤性能都不是太好。目前比較成熟且效果較好的目標(biāo)跟蹤算法是交互多模型(IMM)算法[2]。IMM算法是一種基于混合系統(tǒng)狀態(tài)估計(jì)的次優(yōu)算法。在k時(shí)刻利用交互多模型算法時(shí),認(rèn)為每一個(gè)模型濾波器都有可能成為當(dāng)前有效的系統(tǒng)模型濾波器,每一個(gè)濾波器的初始值都是基于前一個(gè)時(shí)刻所有模型濾波結(jié)構(gòu)的合成。最后的濾波結(jié)果是每一個(gè)濾波器的綜合輸出。

但是,交互多模型算法的性能與模型集的選取與單模型的跟蹤性能有很大關(guān)系。本文以基于勻速運(yùn)動(dòng)、勻加速運(yùn)動(dòng)、勻速轉(zhuǎn)彎和“當(dāng)前”統(tǒng)計(jì)模型這四種跟蹤模型組成的模型集為基礎(chǔ),對(duì)交互多模型算法進(jìn)行了仿真研究。為了在高機(jī)動(dòng)目標(biāo)條件下獲得較高的跟蹤精度,將無味濾波算法(UKF)應(yīng)用到了IMM算法中,并通過仿真對(duì)UKF-IMM算法進(jìn)行了仿真驗(yàn)證。

1 基于無味濾波的交互多模型算法

1.1 交互多模型算法

IMM算法是一種基于軟切換的多模型算法,在該算法中存在多個(gè)目標(biāo)運(yùn)動(dòng)模型,每個(gè)運(yùn)動(dòng)模型對(duì)應(yīng)于不同的目標(biāo)機(jī)動(dòng)輸入水平。在計(jì)算出模型為正確的后驗(yàn)概率之后,就可以通過對(duì)各模型正確的狀態(tài)估計(jì)加權(quán)求和來給出最終的目標(biāo)運(yùn)動(dòng)狀態(tài)估計(jì)。這里加權(quán)因子為模型正確的后驗(yàn)概率。交互多模型的具體步驟[3]如下:

1)計(jì)算目標(biāo)在k時(shí)刻的預(yù)測(cè)位置

F為狀態(tài)轉(zhuǎn)移矩陣,C為權(quán)系數(shù)。

2)相互作用

式中,pij為模型開關(guān)概率,通常認(rèn)為其服從馬爾可夫過程,μi為模型轉(zhuǎn)換概率。

3)濾波

式中,S為濾波量測(cè)預(yù)測(cè)協(xié)方差矩陣,z,h分別為量測(cè)輸入和量測(cè)預(yù)測(cè)值。

4)模型概率更新

5)形成綜合輸出

其中,

1.2 模型集的選取

考慮到目標(biāo)常見的運(yùn)動(dòng)方式有勻速、勻加速和勻速轉(zhuǎn)彎三種,因此在確定模型集時(shí)首先加入以上三種模型。另外,為了在目標(biāo)運(yùn)動(dòng)狀態(tài)發(fā)生突變時(shí)濾波器可以快速跟上運(yùn)送狀態(tài)的突變,在模型集內(nèi)加入了“當(dāng)前統(tǒng)計(jì)模型”。具體的模型[4]為勻速運(yùn)動(dòng)模型、勻加速運(yùn)動(dòng)模型、勻速轉(zhuǎn)彎模型和“當(dāng)前”統(tǒng)計(jì)模型。

1.3 UKF算法

無味濾波算法(UKF)是一種基于Unscented變換的非線性狀態(tài)估計(jì)算法,Unscented變換是一種近似計(jì)算,建立了非線性變換的隨機(jī)變量的統(tǒng)計(jì)特性的新方法,它建立的依據(jù)是:近似一個(gè)概率分布比近似一個(gè)任意的非線性變換(或函數(shù))要容易[4]。UKF算法相比EKF算法能精確估計(jì)任意非線性函數(shù)的二階泰勒近似解;估計(jì)誤差為三階及三階以上高階矩項(xiàng)截?cái)嗾`差。無味濾波算法具體步驟見文獻(xiàn)[5]。

2 仿真驗(yàn)證

2.1 仿真場景

仿真場景如下:雷達(dá)為三坐標(biāo)搜索雷達(dá),數(shù)據(jù)更新周期為2 s;距離測(cè)量精度為150 m;方位與俯仰角測(cè)量精度為0.3°。目標(biāo)軌跡的仿真采用典型反艦導(dǎo)彈對(duì)艦艇的攻擊過程彈道,具體目標(biāo)的運(yùn)動(dòng)仿真參數(shù)如表1所示,目標(biāo)的軌跡如圖1所示。

表1 目標(biāo)軌跡仿真參數(shù)

圖1 目標(biāo)運(yùn)動(dòng)軌跡曲線

2.2 仿真結(jié)果

為了比較采用交互多模型算法的效果,對(duì)基于擴(kuò)展卡爾曼濾波算法的交互多模型(EKFIMM)與基于無味濾波的交互多模型算法(UKFIMM)進(jìn)行比較。經(jīng)過100次蒙特卡羅仿真,得到EKF-IMM和UKF-IMM的仿真結(jié)果如表2所示。

表2 仿真結(jié)果對(duì)比

由表2中數(shù)據(jù)可以看出,采用UKF-IMM算法可以顯著地提高對(duì)高機(jī)動(dòng)目標(biāo)的跟蹤精度。采用UKF-IMM和EKF-IMM算法對(duì)雷達(dá)測(cè)量和目標(biāo)跟蹤在距離、方位、俯仰的誤差進(jìn)行仿真,仿真結(jié)果分別如圖2和圖3所示。由仿真結(jié)果可以看出,在目標(biāo)蛇形機(jī)動(dòng)段,采用EKF-IMM算法在方位上會(huì)出現(xiàn)跟蹤滯后,方位誤差曲線出現(xiàn)震蕩,而采用UKF-IMM算法穩(wěn)態(tài)性能較好。

圖2 UKF-IMM仿真結(jié)果

圖3 EKF-IMM仿真結(jié)果

3 結(jié)束語

針對(duì)實(shí)際運(yùn)動(dòng)中目標(biāo)運(yùn)動(dòng)的復(fù)雜多變,單模型的目標(biāo)跟蹤濾波器無法跟蹤的問題,提出一種基于UKF交互多模型算法。經(jīng)仿真驗(yàn)證,與EKF-IMM算法相比,UKF-IMM算法的跟蹤性能更優(yōu),可以顯著地提高對(duì)高機(jī)動(dòng)目標(biāo)的跟蹤精度。該算法的有效性在實(shí)際的陸地與海上雷達(dá)組網(wǎng)試驗(yàn)中也得到了驗(yàn)證。

[1]張峰,田康生,息木林.彈道導(dǎo)彈運(yùn)動(dòng)建模與跟蹤研究[J].彈箭與制導(dǎo)學(xué)報(bào),2012,32(3):53-58.

[2]張華睿,楊宏文,郁文賢.多目標(biāo)跟蹤時(shí)的雷達(dá)波形調(diào)度問題[J].雷達(dá)科學(xué)與技術(shù),2012,10(6):650-655.ZHANG Hua-rui,YANG Hong-wen,YU Wen-xian.Optimal Waveform Scheduling for Multi-Target Tracking[J].Radar Science and Technology,2012,10(6):650-655.(in Chinese)

[3]BLOM H A P,BAR-SHALOM Y.The Interacting Multiple Model Algorithm for Systems with Markovian Switching Coefficients[J].IEEE Trans on Automatic Control,1988,33(8):780-783.

[4]周猛飛,肖成松,逯振飛,等.基于UKF的滾動(dòng)時(shí)域估計(jì)方法研究[J].計(jì)算機(jī)與應(yīng)用化學(xué),2013,30(12):1406-1410.

[5]LI X R,JILKOV V P.Survey of Maneuvering Target Tracking Part I:Dynamic Models[J].IEEE Trans on Aerospace and Electronic Systems,2003,39(4):1333-1364.