朱紅鵬,黃 勇,修建娟,關(guān) 鍵

(1.海軍航空工程學(xué)院電子信息工程系,山東煙臺(tái)264001; 2.海軍航空工程學(xué)院信息融合研究所,山東煙臺(tái)264001)

0 引言

雷達(dá)觀測(cè)下的弱小目標(biāo)檢測(cè)和跟蹤一直是一個(gè)難點(diǎn)和熱點(diǎn)問(wèn)題。在探測(cè)遠(yuǎn)距離目標(biāo)或在雜波密度比較大的情況下,目標(biāo)回波信號(hào)信噪比比較低,很難實(shí)現(xiàn)對(duì)目標(biāo)的準(zhǔn)確檢測(cè)和穩(wěn)定跟蹤。檢測(cè)前跟蹤方法(TBD)是實(shí)現(xiàn)雷達(dá)弱小目標(biāo)檢測(cè)和跟蹤的有效方法,基本思想是在作檢測(cè)決策前,通過(guò)在目標(biāo)軌跡上回波能量的積累來(lái)提高信噪比,實(shí)現(xiàn)弱小目標(biāo)檢測(cè)并回溯提取其運(yùn)動(dòng)軌跡。TBD是一種信號(hào)處理的思想,需要依賴(lài)具體的方法實(shí)現(xiàn)。

基于隨機(jī)有限集理論(RFS)的多目標(biāo)跟蹤方法能有效避免關(guān)聯(lián)計(jì)算,直接估計(jì)多個(gè)目標(biāo)的狀態(tài),在目標(biāo)數(shù)目時(shí)變的情況下仍然適用,是一種新興的、有前景的多目標(biāo)跟蹤方法。2003年,Mahler基于RFS理論提出了概率假設(shè)密度濾波器(PHD)[1]。此后,Vo給出了PHD濾波器的兩種收斂實(shí)現(xiàn),即序貫蒙特卡洛PHD(SMC-PHD)[2]和高斯混合PHD(GM-PHD)[3]。西安電子科技大學(xué)李翠蕓等針對(duì)多目標(biāo)跟蹤問(wèn)題提出了一種GMPHD前向后向平滑算法[4],利用量測(cè)數(shù)據(jù)對(duì)濾波值進(jìn)行平滑,提升目標(biāo)跟蹤精度。Nandakumaran等提出了SMC-PHD平滑濾波器[5],通過(guò)對(duì)前向?yàn)V波后的粒子進(jìn)行一定步數(shù)的后向平滑,在犧牲算法運(yùn)行效率的基礎(chǔ)上提升了目標(biāo)數(shù)目和目標(biāo)狀態(tài)的估計(jì)精度。

將PHD技術(shù)運(yùn)用于TBD領(lǐng)域還處在起步階段,有很多急需解決的問(wèn)題。Punithakumar最早將SMC-PHD方法運(yùn)用到紅外圖像的多目標(biāo)TBD問(wèn)題中[6],實(shí)現(xiàn)了未知數(shù)目的弱小目標(biāo)的檢測(cè)。在此基礎(chǔ)上,國(guó)防科技大學(xué)林再平在TBD問(wèn)題中引入SMC-PHD平滑濾波器,提升了紅外圖像中目標(biāo)估計(jì)的精度[7]。針對(duì)SMC-PHD-TBD算法估計(jì)精度不高的問(wèn)題,海軍航空工程學(xué)院關(guān)鍵等又提出了一種RBPF-PHD-TBD的方法,將目標(biāo)狀態(tài)空間降維分解,分別用線(xiàn)性和非線(xiàn)性濾波器進(jìn)行跟蹤,提升了算法的估計(jì)性能[8]?；诟倪M(jìn)的SMCPHD-TBD算法和標(biāo)準(zhǔn)SMC-PHD-TBD算法的優(yōu)點(diǎn)和粒子濾波的優(yōu)點(diǎn)一樣,對(duì)付非線(xiàn)性、非高斯情況比較理想,其缺點(diǎn)也和粒子濾波一樣存在粒子退化的問(wèn)題。另外,SMC實(shí)現(xiàn)的方法還存在峰值提取困難的問(wèn)題。高斯混合(GM)是一種使用起來(lái)比較方便的方法,峰值提取在這種方法的模式下實(shí)現(xiàn)起來(lái)非常容易,而且工程實(shí)現(xiàn)相對(duì)于粒子濾波方法要簡(jiǎn)單得多,計(jì)算量也比較小[3]。

本文研究了雷達(dá)觀測(cè)下基于GM-PHD的TBD問(wèn)題,前期研究發(fā)現(xiàn)基于GM-PHD的TBD算法在信噪比較低的情況下存在目標(biāo)數(shù)目估計(jì)不準(zhǔn)確、狀態(tài)估計(jì)精度下降的問(wèn)題。鑒于此,本文提出了一種基于高斯混合概率假設(shè)密度平滑的檢測(cè)前跟蹤算法(SGM-PHD-TBD),首先給出了SGMPHD-TBD算法的目標(biāo)運(yùn)動(dòng)模型和雷達(dá)傳感器觀測(cè)模型,在此基礎(chǔ)上,本文闡述了SGM-PHD-TBD算法的具體實(shí)施步驟。仿真分析表明,該算法在低信噪比條件下能有效提升目標(biāo)數(shù)目和目標(biāo)位置的估計(jì)精度。

1 目標(biāo)運(yùn)動(dòng)模型和觀測(cè)模型

1.1 目標(biāo)運(yùn)動(dòng)模型

基于高斯混合概率假設(shè)密度濾波算法適用于高斯線(xiàn)性的目標(biāo)運(yùn)動(dòng)模型。假設(shè)k時(shí)刻目標(biāo)的數(shù)目為N k,第p個(gè)目標(biāo)的運(yùn)動(dòng)狀態(tài)為X p(k),則雷達(dá)弱小目標(biāo)運(yùn)動(dòng)模型為

1.2 雷達(dá)觀測(cè)模型

文獻(xiàn)[9]建立的雷達(dá)觀測(cè)模型區(qū)別于傳統(tǒng)的紅外觀測(cè)模型,將回波信號(hào)表示為距離-方位多普勒上的三維觀測(cè)。假設(shè)每一掃描周期得到的觀測(cè)數(shù)據(jù)大小為Nr×Nd×Nb個(gè)分辨單元,其中Nr,Nd和Nb分別表示距離單元、多普勒單元和方位單元數(shù),則k時(shí)刻的觀測(cè)值為

式中,觀測(cè)噪聲v k=vI,k+ivQ,k,vI,k和vQ,k均服從獨(dú)立的零均值高斯分布,方差均為σ2。分辨單元(i,j,l)處的觀測(cè)值為

第p個(gè)目標(biāo)在第(i,j,l)個(gè)單元內(nèi)的點(diǎn)擴(kuò)散函數(shù)為

式中:Rc,Dc和Bc表示與距離、多普勒和方位單元尺寸有關(guān)的常數(shù),分別根據(jù)帶寬、積累時(shí)間和波束寬度決定;Lr,Ld和Lb分別表示3個(gè)觀測(cè)維度上的損耗系數(shù);r k,d k和b k表示k時(shí)刻目標(biāo)所處的距離、多普勒和方位單元。

由文獻(xiàn)[10]可知,針對(duì)TBD的“標(biāo)準(zhǔn)”多目標(biāo)觀測(cè)模型下的推得的雷達(dá)觀測(cè)似然函數(shù)為

2 SGM-PHD-TBD算法

在線(xiàn)性高斯分布下,當(dāng)先驗(yàn)分布是混合加權(quán)的高斯和形式時(shí),那么更新的后驗(yàn)分布也可以表述成混合高斯的形式。后驗(yàn)密度的估計(jì)是通過(guò)高斯混合的均值、權(quán)值和協(xié)方差進(jìn)行遞推得到的。因?yàn)楦咚狗至康膫€(gè)數(shù)會(huì)隨著時(shí)間無(wú)限增長(zhǎng),所以需要在算法中對(duì)描述后驗(yàn)概率密度函數(shù)的高斯分量進(jìn)行處理:通過(guò)剪枝操作,去掉權(quán)值低的高斯分量;通過(guò)合并操作,將分布非常接近的高斯分布合并為一個(gè)高斯分量[11]。算法主要步驟如下。

1)初始化

在k=0時(shí)刻,定義初始PHD為

第i個(gè)高斯分量具有狀態(tài)權(quán)值均值和方差

2)預(yù)測(cè)

在不考慮衍生目標(biāo)的情況下,假設(shè)k-1時(shí)刻目標(biāo)的后驗(yàn)概率為以下高斯混合形式:

若vS,k|k-1(x),γk(x)分別表示存活目標(biāo)和新生目標(biāo)的強(qiáng)度函數(shù)的預(yù)測(cè),那么k時(shí)刻的目標(biāo)預(yù)測(cè)狀態(tài)的強(qiáng)度函數(shù)也是高斯混合的形式:

式中,pS,k表示k時(shí)刻目標(biāo)存活的概率,若存活目

3)更新

假設(shè)k時(shí)刻的先驗(yàn)概率密度是高斯混合形式,那么k時(shí)刻的后驗(yàn)概率密度也是高斯混合的。更新后的后驗(yàn)概率密度v k(x)是由未檢測(cè)到的項(xiàng)

(1-pD,k)v k|k-1和|Z k|檢測(cè)到的項(xiàng)vD,k(·;z)組成的,檢測(cè)的任一項(xiàng)z∈Z k如下:

式中,pD,k表示k時(shí)刻目標(biāo)檢測(cè)的概率。若k時(shí)刻雜波的強(qiáng)度函數(shù)為κk,量測(cè)值為z k,量測(cè)噪聲協(xié)方差為R k,那么高斯分量的權(quán)值、均值和協(xié)方差的計(jì)算公式為

4)平滑

定義在高斯混合假設(shè)條件下,D k|l(x)為

則v k|l(x)為高斯混合形式,表示為

本文采用一階平滑的算法實(shí)現(xiàn)雷達(dá)弱小目標(biāo)的檢測(cè)前跟蹤,即l=k+1,平滑公式為

5)合并修剪

在濾波過(guò)程中,高斯分量的個(gè)數(shù)會(huì)隨著時(shí)間無(wú)限增長(zhǎng),因此需要對(duì)描述后驗(yàn)概率密度函數(shù)的高斯分量進(jìn)行處理:通過(guò)剪枝操作去掉權(quán)值低的分量;通過(guò)合并將分布非常接近的高斯分量直接合并成一個(gè)高斯分量。設(shè)合并門(mén)限為U,若第i個(gè)分量和第j個(gè)分量的均值滿(mǎn)足:

就可以把兩者合并成一個(gè)分量。合并以后需要進(jìn)行修剪操作,舍棄權(quán)值低于截?cái)嚅撝礣的高斯分量,權(quán)值更新后重復(fù)步驟2)～步驟5)。

3 仿真與分析

3.1 仿真環(huán)境

仿真設(shè)置目標(biāo)的距離單元均勻分布在[0,10 km],多普勒單元隨機(jī)分布在[-10 m/s, 10 m/s],為了簡(jiǎn)化計(jì)算,方位單元只取一個(gè)波位。設(shè)置雷達(dá)的檢測(cè)概率pD,k=pD=0.98,雜波密度λc=12.5×10-4m-2。信噪比計(jì)算式為SNR=10 lg(P/2σ2),可依據(jù)目標(biāo)功率和信噪比推算背景噪聲方差σ2。假設(shè)此區(qū)域內(nèi)存在4個(gè)目標(biāo),目標(biāo)的幅度是固定的且均為20,目標(biāo)的存活概率PS,k= PS=0.99,共仿真40幀數(shù)據(jù),采樣周期T=1 s。目標(biāo)運(yùn)動(dòng)近似為線(xiàn)性運(yùn)動(dòng),運(yùn)動(dòng)方程同式(1)。目標(biāo)運(yùn)動(dòng)噪聲的功率譜密度q1=0.001,目標(biāo)初始狀態(tài)矩陣及存活時(shí)間如表1所示。

表1 目標(biāo)運(yùn)動(dòng)狀態(tài)表

本文分別以目標(biāo)數(shù)目估計(jì)標(biāo)準(zhǔn)差和位置最優(yōu)子模型分配距離(OSPA)作為目標(biāo)數(shù)目和目標(biāo)狀態(tài)估計(jì)精度的評(píng)價(jià)標(biāo)準(zhǔn)[12],選取的位置估計(jì)OSPA距離參數(shù)為c=5,p=2。在目標(biāo)檢測(cè)和跟蹤問(wèn)題中,信噪比低于10 dB的目標(biāo)被稱(chēng)為弱小目標(biāo),仿真建立在3種不同的環(huán)境下,信噪比分別為9 dB, 8 d B,7 d B。因?yàn)橐淮畏抡娴碾S機(jī)性較大,所以采用100次蒙特卡洛實(shí)驗(yàn)求統(tǒng)計(jì)均值的方法來(lái)驗(yàn)證算法的性能。

3.2 仿真結(jié)果及分析

實(shí)驗(yàn)結(jié)果如圖1～圖3所示。

圖1 SNR=9 dB兩種算法處理效果對(duì)比

圖2 SNR=8 dB兩種算法處理效果對(duì)比

對(duì)比兩種算法處理效果,主要從目標(biāo)數(shù)目估計(jì)準(zhǔn)確度、狀態(tài)估計(jì)精度和算法的運(yùn)行效率三個(gè)方面進(jìn)行分析。

圖3 SNR=7 dB兩種算法處理效果對(duì)比

在目標(biāo)數(shù)目估計(jì)方面。由圖1(a)、圖2(a)和圖3(a)可以看出,兩種算法在信噪比較高時(shí)對(duì)目標(biāo)數(shù)目估計(jì)的準(zhǔn)確度均較高;隨著信噪比的降低,兩種算法的準(zhǔn)確度均逐步下降,且GM-PHD-TBD算法對(duì)目標(biāo)數(shù)目估計(jì)的偏差越來(lái)越大,而SGMPHD-TBD算法保持了對(duì)目標(biāo)數(shù)目良好的估計(jì)性能。圖1(b)、圖2(b)和圖3(b)表示的是兩種算法100次蒙特卡洛實(shí)驗(yàn)?zāi)繕?biāo)數(shù)目估計(jì)的標(biāo)準(zhǔn)差的均值以及40幀數(shù)據(jù)目標(biāo)數(shù)目估計(jì)標(biāo)準(zhǔn)差的均值,統(tǒng)計(jì)值越小代表目標(biāo)數(shù)目估計(jì)越準(zhǔn)確。統(tǒng)計(jì)圖1(b)、圖2(b)和圖3(b)中不同信噪比條件下兩種算法對(duì)目標(biāo)數(shù)目估計(jì)的標(biāo)準(zhǔn)差均值,如表2所示。

表2 不同SNR條件下兩種算法目標(biāo)數(shù)目估計(jì)標(biāo)準(zhǔn)差均值

從表2可以看出,SGM-PHD-TBD算法在SNR=9 dB時(shí),相比GM-PHD-TBD算法對(duì)目標(biāo)數(shù)目估計(jì)的精度可以提升19.4%,在SNR＜9 dB時(shí)可以提升30%以上,說(shuō)明SGM-PHD-TBD算法對(duì)弱小目標(biāo)數(shù)目估計(jì)更有優(yōu)勢(shì)。

在跟蹤精度方面。圖1(c)、圖2(c)和圖3(c)表示的是兩種算法100次蒙特卡洛實(shí)驗(yàn)OSPA均值以及40幀數(shù)據(jù)OSPA的均值,統(tǒng)計(jì)值越小表明目標(biāo)的狀態(tài)估計(jì)越精確。由圖可知,當(dāng)SNR＜9 d B時(shí),SGM-PHD-TBD算法對(duì)目標(biāo)狀態(tài)估計(jì)明顯優(yōu)于GM-PHD-TBD算法。統(tǒng)計(jì)圖1(c)、圖2(c)和圖3(c)不同信噪比條件下兩種算法的OSPA誤差均值,如表3所示。

表3 不同SNR條件下兩種算法OSPA均值

從表3可以看出,SGM-PHD-TBD算法相比GM-PHD-TBD算法在目標(biāo)狀態(tài)的估計(jì)精度方面提升了10%以上,且信噪比越低提升的幅度越大,體現(xiàn)出改進(jìn)算法在低信噪比條件下良好的目標(biāo)狀態(tài)估計(jì)能力。

在運(yùn)算效率方面,由于在原方法上引入了平滑濾波器,所以算法的復(fù)雜度增加了。統(tǒng)計(jì)不同SNR條件下兩種算法100次蒙特卡洛實(shí)驗(yàn)運(yùn)算時(shí)間均值,如表4所示。

表4 不同SNR條件下兩種算法運(yùn)算時(shí)間

從表4可以看出,SGM-PHD-TBD算法相比GM-PHD-TBD算法運(yùn)行時(shí)間增加了50%左右,在對(duì)運(yùn)行時(shí)間要求不高的情況下,采用SGMPHD-TBD算法可以更好地實(shí)現(xiàn)對(duì)雷達(dá)弱小目標(biāo)的檢測(cè)和跟蹤。

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

本文針對(duì)GM-PHD算法在用于TBD過(guò)程中存在信噪比降低時(shí)目標(biāo)數(shù)目估計(jì)不準(zhǔn)、目標(biāo)狀態(tài)估計(jì)精度較低的問(wèn)題,提出了SGM-PHD-TBD算法,在給出該算法標(biāo)準(zhǔn)運(yùn)動(dòng)模型和雷達(dá)傳感器觀測(cè)模型基礎(chǔ)上,闡述了該算法的具體實(shí)施步驟。仿真表明,在不考慮增加時(shí)間成本的前提下,SGMPHD-TBD算法相比GM-PHD-TBD算法在低信噪比條件下對(duì)目標(biāo)數(shù)目估計(jì)的準(zhǔn)確度可以提升30%以上,對(duì)目標(biāo)的狀態(tài)估計(jì)的精度可以提升10%以上,能夠更好地實(shí)現(xiàn)對(duì)雷達(dá)弱小目標(biāo)的準(zhǔn)確檢測(cè)和穩(wěn)定跟蹤。

[1]MAHLER R.Multitarget Bayes Filtering via First-Order Multitarget Moments[J].IEEE Trans on Aerospace and Electronics Systems,2003,39(3):1152-1178.

[2]VO B N,MA W K.A Closed-Form Solution for the Probability Hypothesis Density Filter[C]∥8th International Conference on Information Fusion,Philadelphia,PA:IEEE,2005:856-863.

[3]VO B N,MA W K.The Gaussian Mixture Probability Hypothesis Density Filter[J].IEEE Trans on Signal Processing,2006,54(11):4091-4104.

[4]李翠蕓,江舟,李斌,等.未知雜波環(huán)境的GM-PHD平滑濾波器[J].西安電子科技大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015,42(5):98-104.

[5]NANDAKUMARAN N,PUNITHAKUMAR K, KIRUBARAJAN T.Improved Multi-Target Tracking Using Probability Hypothesis Density Smoothing[C]∥SPIE Conference on Signal and Data Processing of Small Targets,San Diego,CA:SPIE,2007: 1-8.

[6]PUNITHAKUMAR K,KIRUBARAJAN T,SINHA A.A Sequential Monte Carlo Probability Hypothesis Density Algorithm for Multitarget Track-Before-Detect[C]∥SPIE Conference on Signal and Data Processing of Small Targets,Bellingham,Washington: SPIE,2005:587-594.

[7]林再平,周一宇,安瑋,等.基于概率假設(shè)密度濾波平滑器的檢測(cè)前跟蹤算法[J].光學(xué)學(xué)報(bào),2012,32 (10):1-8.

[8]柳超,關(guān)鍵,黃勇,等.基于PHD的多目標(biāo)檢測(cè)前跟蹤改進(jìn)方法[J].雷達(dá)科學(xué)與技術(shù),2016,14(1):1-6. LIU Chao,GUAN Jian,HUANG Yong,et al.An Improved Multitarget Track-Before Detect Algorithm Based on Probability Hypothesis Density Filter[J]. Radar Science and Technology,2016,14(1):1-6.(in Chinese)

[9]樊孔帥,郭云飛,劉鑫,等.弱目標(biāo)檢測(cè)前跟蹤快速實(shí)現(xiàn)算法及性能分析[J].杭州電子科技大學(xué)學(xué)報(bào), 2013,33(3):61-64.

[10]童慧思,張顥,孟華平,等.PHD濾波器在多目標(biāo)檢測(cè)前跟蹤中的應(yīng)用[J].電子學(xué)報(bào),2011,39(9): 2046-2051.

[11]胡瑋靜,陳秀宏.多目標(biāo)跟蹤中一種改進(jìn)的高斯混合PHD濾波算法[J].計(jì)算機(jī)工程與應(yīng)用,2016,52 (2):244-255.

[12]吳鑫輝,黃高明,高俊,等.未知噪聲統(tǒng)計(jì)下多模型概率假設(shè)密度粒子濾波算法[J].控制與決策,2014, 29(3):475-480.