王 楠 孫進平 王文光
(北京航空航天大學(xué)電子信息工程學(xué)院,北京 100191)
近幾年來,多輸入多輸出(Multiple Input Multiple Output, MIMO)雷達因為其獨特的工作體制受到了世界各國學(xué)者們的廣泛研究。不同于傳統(tǒng)的相控陣雷達,MIMO雷達配備多個發(fā)射天線和接收天線,各發(fā)射天線和接收天線相互交叉組成發(fā)射-接收天線組。大量的理論和實踐研究表明,相比于傳統(tǒng)體制的雷達,MIMO雷達能夠提高目標檢測性能[1],改善角度分辨率[2],提供更多的空間自由度和更好的空間覆蓋性[3]等。MIMO雷達根據(jù)天線空間構(gòu)型可以分為兩大類型:第一類為集中式MIMO雷達,其發(fā)射天線和接收天線距離很近,可以看作是在同一個地點。第二類是分布式MIMO雷達,其發(fā)射天線和接收天線距離很遠,所有天線都不在同一個位置。本文主要討論分布式MIMO雷達的擴展目標檢測。
分布式MIMO雷達由于發(fā)射接收天線分布在不同的空間位置,因此具有良好的空間多樣性,可以提高目標的探測性能[4]。文獻[5-7]在均勻雜波環(huán)境下,設(shè)計了基于采樣協(xié)方差矩陣(SCM)的檢測器,對待測單元的目標進行檢測。在SCM檢測器的基礎(chǔ)上,Chong等人設(shè)計了魯棒采樣協(xié)方差矩陣(RSCM)檢測器[8],RSCM檢測器采用復(fù)合高斯分布模擬雜波環(huán)境。文獻[9]采用廣義對稱協(xié)方差模型,用以減少分布式MIMO雷達目標檢測中的訓(xùn)練數(shù)據(jù)量。所有這些算法有一個共同的特點,都是假設(shè)分布式MIMO雷達處在均勻干擾信號的環(huán)境下的。然而在現(xiàn)實情況中,由于分布式MIMO雷達中每一個發(fā)射-接收天線組的地理位置不同,所以他們之間的環(huán)境特性是相互獨立的,如果都采用同樣的干擾信號協(xié)方差矩陣顯然是不合理的。
此外在高分辨率雷達系統(tǒng)應(yīng)用中,擴展目標的檢測一直是一個關(guān)鍵問題。文獻[10]提出了一種改進的基于有序統(tǒng)計的廣義似然比檢驗(OS-GLRT)的穩(wěn)健算法,用于稀疏分布的擴展目標檢測。文獻[11]考慮了高斯白噪聲下擴展目標的檢測,作者假設(shè)干擾信號協(xié)方差矩陣未知,利用先驗信息和訓(xùn)練數(shù)據(jù)估計干擾信號的協(xié)方差矩陣,并采用GLRT和兩步GLRT檢測器對擴展目標進行檢測。而文獻[12]在處理擴展目標信號的時候,將其看作多級子目標信號進行處理。雖然有關(guān)擴展目標檢測的文獻很多,但是在非均勻干擾信號環(huán)境下,有關(guān)分布式MIMO雷達擴展目標檢測的研究卻很少。
本文首先根據(jù)雷達系統(tǒng)和工作環(huán)境的實際特性,建立先驗信息矩陣,然后令各個發(fā)射-接收天線組的干擾協(xié)方差矩陣服從以先驗信息矩陣為基礎(chǔ)的逆Wishart分布。最后根據(jù)傳統(tǒng)的Wald檢測器準則,設(shè)計了一種基于知識輔助的Wald檢測器(KA-Wald),以應(yīng)用于分布式MIMO雷達擴展目標檢測中。仿真結(jié)果表明,在小樣本的情況下,本文設(shè)計的KA-Wald檢測器在檢測性能上優(yōu)于傳統(tǒng)Wald檢測器,而與基于知識輔助的廣義似然比檢驗(KA-GLRT)檢測器[13]相比,檢測性能相近,但是計算效率更高。
文中采用的符號以及含義主要有:(·)T表示轉(zhuǎn)置運算。(·)H表示轉(zhuǎn)置共軛運算。etr{·}表示exp{tr(·)},tr(·)表示矩陣的跡運算,exp{·}表示指數(shù)運算。CN表示復(fù)高斯分布。CW-1表示逆Wishart分布。
由于本文討論的是分布式MIMO雷達的運動擴展目標檢測,雷達結(jié)構(gòu)如圖1所示,所以首先做出如下假設(shè):
圖1 分布式MIMO雷達結(jié)構(gòu)Fig.1 Distributed MIMO radar configuration
(1)分布式MIMO雷達中存在M個發(fā)射天線和N個接收天線;
(2)發(fā)射天線發(fā)射的信號波形相互正交;
(3)每個發(fā)射天線發(fā)射L個脈沖信號;
(4)擴展目標的散射點數(shù)為H,并且H≥1。
假設(shè)經(jīng)匹配濾波后第n個接收天線接收到第m個發(fā)射天線的信號為
Xmn=[xmn,1,xmn,2,…,xmn,H]
(1)
與傳統(tǒng)雷達目標檢測類似,將目標檢測問題轉(zhuǎn)化為如下二元假設(shè)問題:
m=1,2,…,Mn=1,2,…,Nk=1,2,…,K
(2)
其中,αmn=[αmn,1,…,αmn,H]T包含擴展目標各個散射點信號的幅度信息。Nmn=[nmn,1,…,nmn,H]表示待測單元的干擾信號(包含雜波和噪聲),其列向量假定是滿足獨立同分布的零均值高斯復(fù)矢量,即
nmn,h~CN(0,Rmn)
(3)
(4)
pmn=[1,e-j2πfmn,…,e-j2π(L-1)fmn]T
(5)
其中,fmn表示運動目標相對于m發(fā)射-n接收天線組的歸一化多普勒頻率。
對于分布式MIMO雷達來說,如果需要估計出所有發(fā)射-接收天線組的干擾協(xié)方差矩陣,那么需要的訓(xùn)練數(shù)據(jù)量就是傳統(tǒng)檢測算法所需訓(xùn)練數(shù)據(jù)量的M×N倍,這對于一個大型的分布式MIMO雷達系統(tǒng)(發(fā)射天線和接收天線很多)來說顯然是一個很大數(shù)據(jù)量。此外在非均勻干擾信號的環(huán)境下,可用訓(xùn)練數(shù)據(jù)的數(shù)量會大幅度減少,最終導(dǎo)致目標檢測性能的下降[14]。因此,為了減少小樣本帶來的檢測性能損失,引入知識輔助模型,建立關(guān)于干擾信號的先驗信息矩陣是很有必要的。
在選擇先驗信息分布的時候,主要考慮兩個方面的要求。第一,先驗分布要能夠反映發(fā)射-接收天線組干擾協(xié)方差矩陣Rmn的先驗信息;第二,先驗分布不能太復(fù)雜,為了能夠得到后驗分布的解析表達式,通常選取具有共軛特性的先驗分布。
Maio等人在文獻[15]中對非均勻干擾環(huán)境下協(xié)方差矩陣先驗分布建模的研究結(jié)果表明,逆Wishart分布能夠反映實際的環(huán)境信息。因此本文采用逆Wishart分布作為每個發(fā)射-接收天線組干擾協(xié)方差矩陣Rmn的先驗分布。
Rmn~CW-1(μ,(μ-L)Σ)
(6)
其中,Σ是一個L×L的半正定矩陣,包含了Rmn的先驗信息。在實際檢測中,結(jié)合雷達工作參數(shù)、先驗的地面散射系數(shù)、內(nèi)部雜波運動等信息可事先計算出Σ。μ>L,表示逆Wishart分布自由度。增加μ,可以增加干擾協(xié)方差矩陣中包含的先驗信息,從而提高檢測性能(但這種檢測性能的提升不會隨著μ的增加而無限增大)。根據(jù)式(6)可知協(xié)方差矩陣Rmn的概率密度函數(shù)為
(7)
其中
(8)
這里Γ(·)表示伽馬函數(shù)。根據(jù)逆Wishart的統(tǒng)計特性可知
E[Rmn]=Σ
(9)
為方便推導(dǎo)KA-Wald算法的檢測表達式,采用如下參數(shù)表示:
(2)θsmn,h是一個L2維的實列向量,它由協(xié)方差矩陣Rmn的對角線元素和對角線以外元素的實部和虛部組成(重復(fù)的不包含)。
與傳統(tǒng)Wald檢測準則類似,擴展目標的Wald檢測器可采用如下判決式
(10)
(11)
其中
[J-1(θmn,h)]θrmn,h,θrmn,h=
(Jθrmn,h,θrmn,h(θmn,h)-Jθrmn,h,θsmn,h(θmn,h)
(12)
根據(jù)實際的檢測問題可知
Jθrmn,h,θsmn,h(θmn,h)=02,L2
(13)
(14)
所以,式(12)可以簡化為
(15)
(16)
推導(dǎo)KA-Wald檢測器的表達式主要分為兩步:
根據(jù)式(2)所示的二元假設(shè)檢驗?zāi)P图跋嚓P(guān)假定,可知在H1假設(shè)條件下的聯(lián)合概率密度函數(shù)為
f(xmn,h,ymn,h|Rmn,αmn,h)=
(17)
C1=(xmn,h-αmn,hpmn)(xmn,h-αmn,hpmn)H+S
(18)
(19)
根據(jù)式(7)和(17)可以推導(dǎo)出H1假設(shè)條件下αmn,h的似然函數(shù)。
f(xmn,h,ymn,h|αmn,h)=
|C1+(μ-L)Σ|-(K+μ+1)
(20)
f(Rmn|xmn,h,ymn,h,αmn,h)=
(21)
對式(21)取對數(shù),并關(guān)于Rmn求導(dǎo),令求導(dǎo)結(jié)果等于零,可以得到H1假設(shè)條件下Rmn的最大后驗估計。
(22)
(23)
將式(22)帶入式(21)中就可以得到αmn,h的最大似然估計。
(24)
進而得出
(25)
將式(23)和式(25)帶入式(16),可得KA-Wald檢測器判決式為
(26)
為了進一步簡化KA-Wald檢測器的表達式,令
(27)
其中,e=(S+(μ-L)Σ)-1/ 2pmn,z=(S+(μ-L)Σ)-1/ 2xmn,h。根據(jù)矩陣求逆定理,可得
(28)
從而得到KA-Wald檢測器的判決式為
(29)
Yongchan Gao等人[13]最先提出分布式MIMO雷達的擴展目標檢測問題。通過建立先驗信息矩陣,模擬不同發(fā)射-接收天線組的干擾信號特性,最后采用基于知識輔助的廣義似然比檢驗(KA-GLRT)算法檢測擴展目標。KA-GLRT的檢測準則為
(30)
KA-GLRT檢測器判決式為
(31)
對比文獻[13]提出的KA-GLRT檢測器,本文所提KA-Wald檢測器的主要優(yōu)點在于:
(1)KA-GLRT檢測器的數(shù)學(xué)推導(dǎo)比較復(fù)雜,因為KA-GLRT算法采用的是傳統(tǒng)廣義似然比檢驗的檢測準則,所以在推導(dǎo)檢測器表達式的時候,需要推導(dǎo)H1條件和H0條件下關(guān)于αmn,h的似然函數(shù)和協(xié)方差矩陣Rmn。而KA-Wald檢測器只需要推導(dǎo)H1條件下關(guān)于αmn,h的似然函數(shù)和協(xié)方差矩陣Rmn,即KA-Wald算法需要估計的未知參數(shù)更少,推導(dǎo)過程更加簡單。
(2)對比式(29)和式(31)左邊可以看出,KA-Wald檢測器中檢測統(tǒng)計量的計算比KA-GLRT檢測器更為簡潔,實際應(yīng)用中所需運算量更小,計算效率更高。
本章主要通過10/Pfa次獨立仿真實驗來比較:(1)傳統(tǒng)Wald檢測器與KA-Wald檢測器的檢測性能;(2)KA-GLRT檢測器和KA-Wald檢測器的檢測性能和計算效率。
假設(shè)分布式MIMO雷達的發(fā)射天線數(shù)為M=2,接收天線數(shù)為N=2。發(fā)射天線發(fā)射脈沖數(shù)L=20。參數(shù)自由度μ=24。虛警率為Pfa=10-2。脈沖重復(fù)頻率為PRF=500 Hz,信號載頻為fc=1 GHz。目標的運動速度為108 km/h。目標所有的散射點信號幅度相同。先驗信息矩陣Σ為指數(shù)形矩陣,矩陣的第i行j列元素等于ρ|i-j|,實驗選取單元延遲相關(guān)系數(shù)ρ=0.9。
圖2是KA-Wald檢測器和傳統(tǒng)Wald檢測器在不同訓(xùn)練樣本下的檢測性能曲線。從圖2中可以看出,在小樣本的情況下,KA-Wald檢測器的檢測性能要優(yōu)于傳統(tǒng)Wald檢測器。這是因為KA-Wald檢測器引入了先驗信息矩陣,有效地補償了非均勻干擾信號環(huán)境下有限訓(xùn)練樣本數(shù)量帶來的性能損失問題,使得檢測器可以利用較少的訓(xùn)練數(shù)據(jù)獲取較好的檢測性能。
圖2 檢測率隨信噪比變化曲線(H=2)Fig.2 Probability of detection versus SNR (H=2)
圖3是KA-GRLT檢測器和KA-Wald檢測器在不同訓(xùn)練樣本數(shù)據(jù)下的檢測性能,圖4則表示兩種算法在不同訓(xùn)練樣本數(shù)據(jù)下的計算效率。
圖3 目標檢測概率隨SNR變化曲線Fig.3 Probability of detection versus SNR
從圖3中可以看出,(1)KA-Wald檢測器和KA-GLRT檢測器的目標檢測率都隨著信噪比的增加而提高,且上升趨勢相同;(2)KA-Wald檢測器和KA-GLRT檢測器在信噪比相同的情況下,都是訓(xùn)練樣本數(shù)據(jù)越多,目標檢測概率越高;(3)KA-Wald檢測器和KA-GLRT檢測器都可以在少量訓(xùn)練樣本數(shù)據(jù)情況下,獲得良好的目標檢測性能;(4)在同樣的檢測環(huán)境下(信噪比、樣本數(shù)據(jù)、擴展目標散射點數(shù)),KA-Wald算法的目標檢測率與KA-GLRT算法幾乎保持一致,這說明本文設(shè)計的KA-Wald檢測器在檢測性能上與KA-GLRT檢測器相近。
圖4 目標計算效率統(tǒng)計圖Fig.4 Calculation efficiency