基于視覺注意機制的認知雷達數(shù)據(jù)關(guān)聯(lián)算法

王樹亮，畢大平，阮懷林，周陽

國防科技大學(xué) 電子對抗學(xué)院，合肥 230037

雜波背景下的機動目標跟蹤是目標跟蹤領(lǐng)域研究的熱點和難點，其關(guān)鍵是解決目標機動運動狀態(tài)和量測數(shù)據(jù)來源的不確定性[1]。針對機動運動狀態(tài)的不確定性，典型的機動模型構(gòu)建方法[2-3]有Singer模型、Jerk模型、“當(dāng)前”統(tǒng)計(CS)模型、交互式多模型(IMM)和變結(jié)構(gòu)多模型(VSMM)等。而針對量測數(shù)據(jù)來源不確定性，比較典型的算法[4]有概率數(shù)據(jù)關(guān)聯(lián)(PDA)和聯(lián)合概率數(shù)據(jù)關(guān)聯(lián)(JPDA)算法，其中PDA算法容易實現(xiàn)，但局限于目標的關(guān)聯(lián)門不相交或沒有回波處于相交區(qū)域的環(huán)境；JPDA算法跟蹤效果比較理想，但當(dāng)目標數(shù)目和有效量測數(shù)增大時，算法的計算量呈指數(shù)增長甚至出現(xiàn)組合爆炸現(xiàn)象，工程上不易實現(xiàn)。

從工程應(yīng)用角度出發(fā)，將IMM和PDA相結(jié)合的IMM-PDA算法為當(dāng)前雜波背景下機動目標跟蹤普遍采用的方法[5]，然而IMM中各個子濾波器由于使用不同關(guān)聯(lián)波門可能導(dǎo)致目標失跟。為此，潘泉等學(xué)者提出了綜合IMM-PDA[6]算法，算法中各個子濾波器共享一個關(guān)聯(lián)門，即子濾波器的候選回波相同，從而提高了數(shù)據(jù)關(guān)聯(lián)的效率。在數(shù)據(jù)關(guān)聯(lián)算法中，關(guān)聯(lián)波門的設(shè)計是影響雷達計算資源的重要因素，若關(guān)聯(lián)波門設(shè)計過大，較多量測落入波門，計算消耗較大；若設(shè)計過小，當(dāng)目標發(fā)生機動變化時，容易造成目標失跟，因此有必要設(shè)計一種能夠隨目標機動動態(tài)調(diào)整的自適應(yīng)關(guān)聯(lián)波門。另外，IMM-PDA算法在應(yīng)用中還存在IMM模型集選取以及傳統(tǒng)PDA算法不適于多交叉目標跟蹤等問題，需要不斷進行研究改進。

近年來，結(jié)合認知科學(xué)和仿生學(xué)的相關(guān)思想，加拿大學(xué)者Haykin等首先明確提出了認知雷達的概念[7-8]，其本質(zhì)特征是融入人類“記憶、注意”等認知過程，通過對目標和環(huán)境特性的實時感知，自適應(yīng)調(diào)整雷達發(fā)射端波形結(jié)構(gòu)以及接收端數(shù)據(jù)處理模型[9-14]，從而增強了雷達智能化處理能力。人類“視覺注意”[15-19]的選擇性“過濾器”信息處理機制決定了人類能夠從眾多視覺信息中快速地選擇那些最重要、最有用、與當(dāng)前行為最相關(guān)的感興趣的視覺信息，而視覺分析中的“特征整合”理論又能夠幫助人類快速對目標進行識別。將人類視覺注意機制應(yīng)用到復(fù)雜背景下機動目標跟蹤的目的，就是通過設(shè)計自適應(yīng)的關(guān)聯(lián)波門使雷達保持對目標不失跟的前提下使雷達計算資源耗費較小，并且通過“特征整合”理論將不同關(guān)聯(lián)波門交叉區(qū)域公共量測進行分類，增強傳統(tǒng)PDA算法的環(huán)境適應(yīng)性。

基于以上文獻研究成果，本文重點結(jié)合人類“視覺注意”機制，以綜合IMM-PDA算法為基礎(chǔ)，研究雜波背景下雷達接收端數(shù)據(jù)處理模型的改進，① 采取假設(shè)目標最大機動水平已知的CS模型和勻速運動(CV)模型作為IMM的模型集，關(guān)聯(lián)波門中心和大小隨目標機動變化動態(tài)調(diào)整，從而使雷達計算資源消耗和跟蹤成功率得到兼顧；② 在傳統(tǒng)利用目標位置特征基礎(chǔ)上，進一步提取目標運動特征，在計算量增加不大的情況下，增強了傳統(tǒng)PDA算法的環(huán)境適應(yīng)性。

1 雷達跟蹤模型

目標運動的離散狀態(tài)方程表示為

X(m)(k)=F(m)(k|k-1)X(m)(k-1)+

(1)

pmn(k-1)?p{Mn(k)|Mm(k-1)}

(2)

目標的量測方程為

Y(k)=h(X(k))+V(k)

(3)

(4)

2 自適應(yīng)關(guān)聯(lián)門

注意的選擇性過濾器模型認為信息的加工需要經(jīng)歷“感覺記憶、過濾器、探測器和記憶”4個階段，如圖1所示[20]

過濾器能夠根據(jù)刺激的物理特性來識別注意到的信息，同時僅讓注意到的信息進入下一個階段的探測器，從而節(jié)約人類感知資源。雜波背景下的目標跟蹤實質(zhì)上也是一個特殊的信息加工過程，雷達接收機對掃描區(qū)域內(nèi)獲取的大量包含目標和雜波的數(shù)據(jù)進行處理，根據(jù)目標的運動特性設(shè)置一個類似“過濾器”的關(guān)聯(lián)波門，對落入關(guān)聯(lián)門內(nèi)的量測按照一定的數(shù)據(jù)關(guān)聯(lián)算法(類似探測器的加工過程)對目標狀態(tài)進行估計，最后將結(jié)果進行存儲(類似于人類“記憶”)，以形成目標的軌跡。

圖1 注意過濾器模型[20]Fig.1 Filtering model of attention[20]

2.1 IMM模型集的選取

IMM所包含的模型集，可以隨目標機動進行切換，然而其在模型集的選擇上存在這樣一個困境，即為了提高濾波精度需要盡可能多地增加模型，但是太多的模型除了增加計算量之外，反而降低了估計器的性能[21]。于是很多學(xué)者嘗試將適于強機動目標的速度自適應(yīng)濾波模型和CV模型作為IMM的模型集來取得跟蹤性能的折衷[12,21],本文正是基于這種思想，將速度自適應(yīng)的CS模型和CV模型作為模型集，由此帶來的好處是CS模型可以為了確保能及時地追蹤戰(zhàn)場目標的“機動”變化，而取較大的加速度極限值，所帶來的弱機動目標跟蹤精度不高的問題通過CV模型的交互競爭予以彌補。

2.2 基于傳統(tǒng)IMM-PDA的關(guān)聯(lián)波門

關(guān)聯(lián)波門當(dāng)前k時刻波門中心和大小，分別由運動模型的一步預(yù)測中心zk|k-1和新息協(xié)方差陣Sk來決定，落入波門的有效量測集合可表示為[22]

(5)

式中：zk為量測；γ為門限。傳統(tǒng)IMM-PDA算法中各模型使用各自的預(yù)測中心和新息協(xié)方差陣所形成的關(guān)聯(lián)波門。

1) 基于CV模型的關(guān)聯(lián)波門

CV模型結(jié)構(gòu)簡單，當(dāng)目標始終保持勻速運動，無機動時，可以有效跟蹤目標，且計算資源消耗較小。然而，當(dāng)目標發(fā)生機動時，這種簡單處理可能會導(dǎo)致關(guān)聯(lián)門內(nèi)無任何量測，從而引起目標失跟。圖2為CV模型跟蹤示意圖(目標在k-2時刻至k-1時刻勻速直線運動，k-1時刻至k時刻做轉(zhuǎn)彎機動)。

圖2 CV模型跟蹤示意圖Fig.2 Diagram of tracking with CV model

2) 基于CS模型的關(guān)聯(lián)波門

CS模型是一種參數(shù)自適應(yīng)機動模型[23]，基本思想是在每一種具體的戰(zhàn)術(shù)場合，人們所關(guān)心的僅是機動加速度的“當(dāng)前”概率密度，即目標機動的當(dāng)前可能性，并認為其概率密度近似服從修正的瑞利分布。

“當(dāng)前”加速度a(k)為正時，有

(6)

“當(dāng)前”加速度a(k)為負時，有

(7)

其中：α為機動頻率常數(shù)；Q0為一對稱矩陣，詳見文獻[23]實際應(yīng)用時為了確保對機動目標的有效跟蹤，往往取較大的加速度極限值，所帶來的問題是對弱機動目標的跟蹤精度反而下降，且關(guān)聯(lián)波門體積維持在一個相對較大的水平。圖3為CS模型跟蹤示意圖。

通過對CV模型和CS模型的跟蹤機制分析，可以看出基于傳統(tǒng)IMM-PDA算法框架，可能存在的問題有：① 在目標發(fā)生機動時，由于CV模型預(yù)測波門內(nèi)無有效量測，導(dǎo)致跟蹤失敗；② 在目標處于勻速運動或弱機動時，由于CS模型取較大的關(guān)聯(lián)波門可能使落入波門內(nèi)的量測過多，一方面容易造成關(guān)聯(lián)跟蹤精度下降，另一方面對雷達計算資源消耗較大。

圖3 CS模型跟蹤示意圖Fig.3 Diagram of tracking with CS model

2.3 基于綜合IMM-PDA的自適應(yīng)關(guān)聯(lián)門設(shè)計

綜合IMM-PDA算法[6]，首先，對各個子濾波器k-1時刻的狀態(tài)進行輸入交互；然后，對各個子濾波器執(zhí)行一步預(yù)測，獲取各子濾波器的量測預(yù)測值和新息方差；最后，根據(jù)模型預(yù)測概率進行加權(quán)。采取CS模型和CV模型作為模型集的綜合量測預(yù)測值(即關(guān)聯(lián)門中心):

(8)

(9)

Vk=cnzγnz/2|Sk|

(10)

式中：nz為目標量測向量維數(shù)；cnz為與量測向量維數(shù)有關(guān)的一個常量。由式(8)～式(10)容易看出，關(guān)聯(lián)門的中心位置和大小取決于兩個模型的實時預(yù)測概率、預(yù)測位置及新息協(xié)方差矩陣。若當(dāng)前目標處于高度機動狀態(tài)，則在綜合IMM-PDA的模型交互中，CS模型占主導(dǎo)地位，模型交互后的綜合新息方差較大，即波門擴大；反之當(dāng)目標處于弱機動或勻速運動時，其波門相應(yīng)縮小，由此實現(xiàn)了關(guān)聯(lián)門的自適應(yīng)，圖4為自適應(yīng)關(guān)聯(lián)波門原理示意圖(目標在k-1時刻至k+1時刻轉(zhuǎn)彎機動，k+1時刻至k+2時刻做勻速直線運動)。

圖4 自適應(yīng)關(guān)聯(lián)波門原理示意圖Fig.4 Diagram of adaptive association gate principle

3 基于“特征整合”的優(yōu)化PDA

Treisman和Gormican[19]在視覺分析過程中提出了特征整合理論。圖5給出了客體被知覺的過程圖。

在Treisman的理論中，客體圖像加工過程的第一步是前注意階段，此時客體(如一個紅色的滾動的圓球)被分解為獨立的特征，如顏色、形狀、運動方式等，然后在集中注意階段，各獨立的特征經(jīng)過特征整合成為一個完整的客體。

本文基于Treisman的特征整合理論，綜合考慮目標位置特征和運動特征，得到一種優(yōu)化的PDA(OPDA)算法。算法步驟為

步驟1位置特征

圖5 Treisman理論的特征整合Fig.5 Treisman’s theory of feature integration

(11)

若目標第k時刻的距離、方位新息方差為Sk,則第g個有效回波量測隸屬于目標i的隸屬度為

(12)

步驟2運動特征

(13)

(14)

圖6 目標運動方向制約圖Fig.6 Direction restriction of target motion

4 仿真實驗

進行M次Monte Carlo仿真實驗，算法性能評價指標選?。孩?距離、跟蹤速度均方根誤差(RMSE)；② 目標跟蹤失敗率，若第L(L≤M)次跟蹤時某采樣時刻距離估計偏離真實位置大于或等于3σr(其中σr為雷達傳感器的距離測量誤差標準差)，就認為第L次跟蹤失敗；③ 計算耗時。

1) 目標運動軌跡

目標1初始位置為(1.25×104,1.5×104) m，初始速度為(-100，-50) m/s，目標在前17 s做勻速直線運動，從第18～43 s做角速度ω=4.77 (°)/s的勻角速度左轉(zhuǎn)彎運動，第44～53 s繼續(xù)做勻速直線運動，第54～80 s做角速度ω=4.77 (°)/s的勻角速度右轉(zhuǎn)彎運動，從第81～100 s再做勻速直線運動。

目標2初始位置為(1.0×104,1.5×104) m，初始速度為(100，-50) m/s，目標在前17 s做勻速直線運動，從第18～43 s做角速度ω=4.77 (°)/s的勻角速度右轉(zhuǎn)彎運動，第44～53 s繼續(xù)做勻速直線運動，第54～80 s做角速度ω=4.77 (°)/s的勻角速度左轉(zhuǎn)彎運動，從第81～100 s再做勻速直線運動。

2) 模型參數(shù)

假設(shè)CS模型最大加速度為30 m/s2，機動頻率常數(shù)選為1/60。IMM算法中兩模型的初始概率為1/2，模型轉(zhuǎn)移概率矩陣為

(15)

3) 關(guān)聯(lián)波門參數(shù)

假設(shè)目標的檢測概率為1，門概率為0.99，選取橢圓波門，關(guān)聯(lián)門限取4個Sigma點。

4) 量測參數(shù)

假設(shè)距離、速度和方位的測量精度分別為50 m，5 m/s和0.1°，雜波密度為ρ(每單位體積內(nèi)的虛假量測數(shù)，本文中單位取每平方千米)。

4.1 自適應(yīng)關(guān)聯(lián)波門算法驗證

比較傳統(tǒng)IMM-PDA和綜合IMM-PDA算法性能。圖7給出了綜合IMM-PDA算法中CS模型和CV模型的實時概率，圖8為橢圓關(guān)聯(lián)波門隨時間變化的面積曲線。

表1給出了兩種算法在不同雜波密度背景下，500次跟蹤中的失敗跟蹤次數(shù)，圖9給出了3種雜波密度背景下(雜波密度1、2、3分別對應(yīng)ρ=0.01,ρ=0.10和ρ=0.50)的計算時間直方圖對比。

由表1可知，當(dāng)取雜波密度ρ=0.10時，傳統(tǒng)IMM-PDA算法的目標失跟率達到了72.2%，已嚴重失效，從圖9計算時間耗費來看，傳統(tǒng)IMM-PDA算法中由于兩個模型分別在各自關(guān)聯(lián)波門內(nèi)進行量測關(guān)聯(lián)，計算量較大。本文基于綜合IMM-PDA框架下波門自適應(yīng)算法，能夠根據(jù)目標機動自適應(yīng)改變預(yù)測中心和新息協(xié)方差(如圖7所示，兩模型概率隨目標機動調(diào)整，圖8為其關(guān)聯(lián)波門面積變化曲線)，使目標保持了較好跟蹤效果(失跟率為1.8%)，且由于兩模型共用同一個波門內(nèi)量測，計算耗時明顯較低。

圖7 綜合IMM-PDA模型的實時概率Fig.7 Real time probability of integrated IMM-PDA

圖8 橢圓關(guān)聯(lián)波門的面積曲線Fig.8 Curve of area of association gate of ellipse

表1 目標跟蹤失敗次數(shù)Table 1 Number of target tracking loss

圖9 兩種算法的計算耗時直方圖Fig.9 Histogram of computing time of two algorithms

4.2 OPDA算法驗證

4.2.1 OPDA算法和傳統(tǒng)PDA算法性能比較

表2為綜合IMM-PDA和綜合IMM-OPDA算法在不同雜波密度背景下，500次跟蹤中的失敗跟蹤次數(shù)，圖10給出了相應(yīng)的計算耗時直方圖對比(雜波密度1、2、3分別對應(yīng)ρ=0.01,ρ=0.10和ρ=0.50)。圖11和圖12給出了兩種算法在雜波密度為1時，目標不失跟情況下的距離估計RMSE對比曲線和速度估計RMSE對比曲線。

傳統(tǒng)PDA算法在跟蹤交叉機動目標時，由于關(guān)聯(lián)波門重疊，貝葉斯全概率公式應(yīng)用前提被破壞，由表2看出，在雜波密度為ρ=0.50時，其目標失跟率高達18.4%。OPDA算法提取目標運動特征，并基于特征整合理論對重疊區(qū)域公共量測進行分類，由表2看出，其跟蹤失敗率得到明顯改善，而且由圖11和圖12距離、速度估計RMSE對比可以看出，在穩(wěn)定不失跟情況下，其跟蹤精度也明顯優(yōu)于PDA算法。

表2 目標跟蹤失敗次數(shù)Table 2 Number of target tracking loss

圖10 不同雜波密度計算耗時直方圖Fig.10 Histogram of computing time of different clutter density

圖11 兩種算法距離估計RMSE曲線Fig.11 Curves of range estimation RMSE of two algorithms

圖12 兩種算法速度估計RMSE曲線Fig.12 Curves of velocity estimation RMSE of two algorithms

由圖10算法在不同雜波背景下的計算耗時直方圖對比，可以看出綜合IMM-OPDA算法相比綜合IMM-PDA算法的計算量增加不大。可以認為本文提出的綜合IMM-OPDA算法是一種在計算代價和穩(wěn)定有效跟蹤之間的折衷選擇。

4.2.2 采樣周期T(或數(shù)據(jù)率)對算法性能影響

圖13為綜合IMM-OPDA算法在不同雜波密度背景下，隨采樣周期T的跟蹤失敗次數(shù)對比曲線。圖14給出了不同采樣周期(采樣周期1，2，3分別對應(yīng)T=0.5、1.0和1.5 s)下的跟蹤計算耗時直方圖對比。

圖13 跟蹤失敗次數(shù)隨采樣周期的變化曲線Fig.13 Curves of number of tracking loss vs sampling period

圖14 不同采樣周期對應(yīng)計算耗時直方圖Fig.14 Histogram of computing time of different sampling period

OPDA算法應(yīng)用的關(guān)鍵是基于“目標運動方向在短時間內(nèi)變化應(yīng)該不會太大”這一思想，如果采樣周期較大(如采樣周期3時)，在雜波密度較小(ρ=0.01和ρ=0.10)時，其跟蹤失敗率仍能保持在10%以內(nèi)，但當(dāng)雜波密度較大(ρ=0.50)時，其跟蹤失敗率達到32.4%。然而如果采樣周期過小(如采樣周期1時)，如圖14所示，其計算代價也較高，因此在實際工程應(yīng)用時要綜合考慮精度、時間等要求進行折衷選擇。

5 結(jié) 論

當(dāng)前，雷達正在向智能化方向發(fā)展，本文提出的基于視覺注意機制的雷達數(shù)據(jù)處理算法正是對這一趨勢的階段性探索。主要研究結(jié)論有：

1) 基于注意選擇性“過濾器”信息處理機制，通過CS模型與CV模型實時交互的綜合IMM-PDA算法，能夠根據(jù)目標機動的變化自適應(yīng)調(diào)整關(guān)聯(lián)波門，不僅提高了跟蹤成功率，而且使雷達計算耗時明顯降低。

2) 基于注意“特征整合”理論，提取目標位置特征和運動特征，將關(guān)聯(lián)門交叉區(qū)域內(nèi)的量測進行分類，優(yōu)化了傳統(tǒng)PDA算法，在計算資源消耗沒有太大增加的基礎(chǔ)上，使算法的環(huán)境適應(yīng)性得到明顯增強，可以認為其是一種在計算代價和穩(wěn)定有效跟蹤間的折衷選擇。

3) 綜合IMM-OPDA算法在工程應(yīng)用時要根據(jù)需要盡可能選擇合適的采樣周期，以平衡跟蹤精度和雷達計算耗時。