費太勇,譚賢四,曲智國,鄒立平,王 紅

(1.空軍預(yù)警學(xué)院,湖北武漢430019;

2.中國人民解放軍93253部隊,遼寧大連116019)

0 引言

多輸入多輸出(MIMO)雷達是利用多個發(fā)射天線同步地發(fā)射正交的波形,同時使用多個接收天線接收回波信號,并集中處理的一種新型雷達體制[1]。其中,統(tǒng)計MIMO雷達(以下簡稱MIMO雷達)具有良好的空間分集增益、結(jié)構(gòu)增益和極化分集增益[2]。眾多學(xué)者對MIMO雷達通過空間分集增益來提高雷達系統(tǒng)的檢測性能進行了深入的研究[3-14]。然而這些研究包含了一個重要的前提條件,即MIMO雷達探測目標時信道間的距離(發(fā)射天線到目標再到接收天線的距離)近似相等。事實上,MIMO雷達探測目標時并不一定處處滿足上述條件,尤其是當MIMO雷達波長較長時為獲得空間分集增益MIMO雷達天線間的間距較遠,此時MIMO雷達在探測某些位置的目標時信道間的距離差異較大而不能再認為近似相等。本文將對這種信道間的距離差異對MIMO雷達檢測性能的影響進行研究。

1 MIMO雷達信號模型

假設(shè)MIMO雷達有M個發(fā)射天線和N個接收天線,其中第m個發(fā)射天線記為T m,第n個接收天線記為Rn,目標記為o′。根據(jù)雙基地雷達方程[15], MIMO雷達接收天線R n與發(fā)射天線T m形成的信道nm接收到的目標回波復(fù)信號可表示為

式中:PT為MIMO雷達總的發(fā)射脈沖功率;GTGR分別為發(fā)射天線和接收天線的增益;λ為信號波長;r o′Tm為o′到T m的距離,r o′Rn為o′到R n的距離;ζmn為目標在nm信道中的RCS所對應(yīng)的電壓散射系數(shù)[4],假設(shè)目標雷達截面積(Radar Cross-Section,RCS)為σ0且為Swerling 1型,則ζnm服從均值為0、方差為σ0的復(fù)高斯分布[16],假設(shè)MIMO雷達收發(fā)天線的位置滿足文獻[3]中空間分集的條件,則各信道的電壓散射系數(shù)是獨立同分布的;s m(t)為T m發(fā)射的信號,‖s m(t)‖2=1;τTm,τRn為o′到T m,R n的距離所對應(yīng)的信號時延;rT,rR分別為o′到發(fā)射天線幾何中心和接收天線幾何中心的距離;Δr Tm=r o′Tm-rT,Δr Rn=r o′Rn-rR。

式(1)可進一步寫成

2 MIMO雷達的似然比檢測

2.1 MIMO雷達信號處理流程

MIM O雷達信號處理流程如圖1所示。接收天線R n接收的信號可表示為

式中,n n(t)為高斯白噪聲。由于發(fā)射信號的正交性,N個接收天線接收的信號通過匹配濾波器組,可將M個發(fā)射信號所對應(yīng)的目標回波進行分離,從而得到N×M個不相關(guān)的目標回波,即對目標形成N×M個獨立的觀測通道。由于目標回波在每個通道的時延可能是不同的,所以需要對目標進行空間配準和包絡(luò)對齊[7,17],然后再進行似然比檢測。

圖1 MIMO雷達信號處理流程圖

2.2 MIMO雷達的似然比函數(shù)

當某檢測單元存在目標時,nm信道在該檢測單元的采樣值可近似為

式中:ET=PTTs為MIMO雷達發(fā)射脈沖的能量;Ts為發(fā)射脈沖寬度;n nm為高斯噪聲,方差為σn。故相應(yīng)的信噪比可表示為

MIMO雷達的N×M個通道在目標檢測單元的采樣值的求模輸出為

由于統(tǒng)計MIMO雷達各信道的電壓散射系數(shù)是獨立同分布的,故無信號(H0)和有信號(H1)時x的概率密度函數(shù)為

相應(yīng)的似然比函數(shù)為

對式(8)兩邊取對數(shù),得對數(shù)似然函數(shù)比為

2.3 檢驗統(tǒng)計量的分布

當檢測單元無信號時:

根據(jù)文獻[18]的結(jié)論,獨立同分布的指數(shù)分布隨機變量的加權(quán)和近似服從伽馬分布,即

式中,伽馬分布X～Γ(α,β)的概率密度函數(shù)為

由式(10)、(13)、(17)得

3 信道距離差異對檢測性能的影響

式(2)、(5)表明信道間的距離差異會導(dǎo)致信道間目標回波信噪比的起伏,這種信噪比的起伏可用向量表示為

當MIMO雷達信道間無距離差異時,c為全1向量,此時無信噪比的起伏。距離差異越大,信噪比的起伏越大,為此定義信噪比的起伏度為

由式(20)的推導(dǎo)過程可知,檢驗統(tǒng)計量TMIMO中的參數(shù)α0,β0,α1,β1都是c的函數(shù),在參考信噪比ρ一定的條件下,TMIMO的檢測性能完全由c確定。為了定量研究c對TMIMO檢測性能的影響,可采用文獻[3,19]中的偏離系數(shù)d來對TMIMO的檢測性能進行度量,即

d越大表示檢驗統(tǒng)計量TMIMO越容易區(qū)分H0假設(shè)和H1假設(shè),即檢測性能越好。

因為伽馬分布Γ(α,β)的均值為E(X)=αβ,方差為D(X)=αβ2,由式(13)～(15)和式(17)～(19),式(23)可變?yōu)?/p>

當c為全1向量,Δ=0,式(24)變?yōu)?/p>

式(26)與文獻[3]中的式(41)是一樣的,這進一步驗證了似然比函數(shù)及其統(tǒng)計特性推導(dǎo)的正確性。

4 仿真實例與數(shù)據(jù)分析

在xoy平面內(nèi),假設(shè)MIMO雷達的天線布置如圖2所示,天線直角坐標分別為R1(-20 km,0),T1(-10 km,0),R2(0,0),T1(10 km,0),R3(20 km, 0),收發(fā)天線幾何中心均為坐標原點o。假設(shè)目標以50 km的半徑繞原點飛行,則目標的極坐標為o′(50 km,θ)。目標方位角θ與信噪比起伏度Δ的關(guān)系如圖3所示。

圖2 MIMO雷達天線位置布置和目標探測示意圖

圖3 目標方位角與信噪比起伏度的關(guān)系

由于MIMO雷達天線是線性布置,當目標位于線陣的法線方向時,信道間的距離差異最小,此時信噪比起伏度最小;當目標偏離法線的角度越大,信道間的距離差異越大,此時信噪比的起伏度越大。

當參考信噪比分別為ρ=3,ρ=5,ρ=10時,信噪比起伏度與檢驗統(tǒng)計量TMIMO的偏離系數(shù)d的關(guān)系曲線如圖4所示。從圖4可以看出,隨著信號起伏度的增加,偏離系數(shù)d都有不同程度的下降,即檢測性能有不同程度的下降。但參考信雜比ρ不同偏離系數(shù)d下降的幅度是不同的。下面進一步考察偏離系數(shù)d與參考信噪比的關(guān)系。

圖4 信噪比起伏度與偏離系數(shù)的關(guān)系

當目標方位角分別為θ=π/2,θ=5π/18,θ=π/9時,由圖3可知其對應(yīng)信噪比起伏度分別為Δ=0.08,Δ=0.42,Δ=0.84,此時偏離系數(shù)d隨參考信雜比ρ的變化如圖5所示。當參考信噪比(ρ＜2)較小時,信噪比起伏度對偏離系數(shù)d的影響不大,而且此時信噪比起伏度越大偏離系數(shù)d反而越大;但隨著參考信噪比的增加,信噪比起伏度對偏離系數(shù)d的影響逐漸增大,例如當ρ=10 dB時,Δ=0.84較Δ=0時的偏離系數(shù)從d≈10下降到d≈6,下降幅度超過1/3。

5 結(jié)束語

圖5 參考信噪比與偏離系數(shù)的關(guān)系

文章研究了MIMO雷達信道間距離差異對檢測性能的影響,研究表明MIMO雷達信道間距離差異會導(dǎo)致信道間信噪比的起伏,進而影響MIMO雷達系統(tǒng)的檢測性能。在信道信噪比較低時這種影響較小,但隨著信噪比的提高,這種影響較大,會導(dǎo)致檢測性能嚴重下降。故MIMO雷達在進行天線布置時,既需要考慮增大天線間距來獲得空間分集增益,又需要考慮增大天線間距所帶來的目標探測時信道間的距離差異以及對檢測性能的影響。研究結(jié)果對統(tǒng)計MIMO雷達的系統(tǒng)配置和實際應(yīng)用具有參考價值。

[1]陳浩文,黎湘,莊釗文,等.多發(fā)多收雷達系統(tǒng)分析與應(yīng)用[M].北京:科學(xué)出版社,2015:2-5.

[2]陳浩文,黎湘,莊釗文.一種新興的雷達體制:MIMO雷達[J].電子學(xué)報,2012,40(6):1190-1198.

[3]FISHLER E,HAIMOVICH A,BLUM R,et al. Spatial Diversity in Radars:Models and Detection Performance[J].IEEE Trans on Signal Processing, 2006,54(3):823-838.

[4]湯俊,伍勇,彭應(yīng)寧,等.MIMO雷達對空域Rician起伏目標檢測性能研究[J].中國科學(xué)F輯:信息科學(xué), 2009,39(8):866-874.

[5]何茜.MIMO雷達檢測與估計理論研究[D].成都:電子科技大學(xué),2010.

[6]廖羽宇.統(tǒng)計MIMO雷達檢測理論研究[D].成都:電子科技大學(xué),2012.

[7]CHELLAPPA R,THEODORIDIS S.Academic Press Library in Signal Processing,Volume 2:Communications and Radar Signal Processing[M].Oxford:Academic Press,2014:667-713.

[8]LIU Jun,LI Hongbin,HIMED B.Persymmetric Adaptive Target Detection with Distributed MIMO Radar[J].IEEE Trans on Aerospace and Electronic Systems,2015,51(1):372-381.

[9]SORIA F R C,GARCIA J S,BARBOZA F M M. Improved Detection of SM-SMux Signals for MIMO Channels[J].IEEE Latin America Transactions, 2015,13(1):43-47.

[10]LIU Su,ZHANG Gong,ZHANG Jiandong,et al. Transmit and Receive Gain Optimization for Distributed MIMO Radar[J].Wireless Personal Communications,2015,85(4):1969-1986.

[11]谷文堃,王黨衛(wèi),馬曉巖,等.分布式OFDM-MIMO雷達非相參積累目標檢測方法[J].系統(tǒng)工程與電子技術(shù),2015,37(10):2266-2271.

[12]谷文堃,王黨衛(wèi),鄭岱堃,等.分布式OFDM-MIMO雷達高速微弱目標多幀檢測[J].華中科技大學(xué)學(xué)報(自然科學(xué)版),2016,44(2):65-69.

[13]王懷軍,許紅波,陸珉,等.MIMO雷達技術(shù)及其應(yīng)用[J].雷達科學(xué)與技術(shù),2009,7(4):245-249.

[14]馬紅星,陳思佳.水面雜波背景的統(tǒng)計MIMO雷達檢測方法研究[J].雷達科學(xué)與技術(shù),2015,13(5): 479-484.

[15]CHERNYAK V S.Fundamentals of Multisite Radar Systems:Multistatic Radar and Multiradar Systems [M].Boca Raton,FL:CRC Press,2011:25-36.

[16]RICHARDS M A.Fundamentals of Radar Signal Processing[M].New York:McGraw-Hill Education,2005:41-83.

[17]谷文堃,王黨衛(wèi),彭尚,等.分布式OFDM-MIMO雷達包絡(luò)對齊方法[J].華中科技大學(xué)學(xué)報(自然科學(xué)版),2014,42(10):19-24.

[18]BOX G E P.Some Theorems on Quadratic Forms Applied in the Study of Analysis of Variance Problems,I.Effect of Inequality of Variance in the One-Way Classification[J].Annals of Mathematical Statistics,1954,25(2):290-302.

[19]SCHARF L L.Statistical Signal Processing:Detection,Estimation and Time Series Analysis[M].2nd ed.New York:Pearson Education,2002.