基于子陣的FDA-MIMO雷達去柵瓣方法

陳 浩， 葉澤浩， 吳華新， 謝 讜， 吳彩華

(空軍預警學院, 湖北武漢 430019)

0 引言

與相控陣相比，頻率分集陣列(Frequency Diverse Array, FDA)在陣元間引入了很小的載頻差(即頻率間隔)，使得FDA可將波束指向相同角度、不同的距離單元，即會形成距離-角度耦合波束[1]，然而，由于FDA雷達的天線方向圖在距離維上的周期性，導致在探測空間內(nèi)很容易形成多個柵瓣，主波束空間能量分散嚴重，干擾落在柵瓣的幾率大，極大地影響了其抗干擾的性能和目標探測能力。

當FDA采用對數(shù)[2]、平方或立方[3]、正弦[1]等非線性形式的頻率間隔以及整數(shù)頻偏策略[4]時，雖然可去除柵瓣，形成單峰值的“點”狀波束，但這種波束的距離維副瓣很高，抗距離依賴性干擾性能和目標距離維估計精度都很差，且對于頻率誤差很是敏感。為此，Xiong等提出可采用遺傳算法優(yōu)化頻率間隔的設置[5]，Wang提出了一種基于球函數(shù)的波束綜合方法[6]，Shao和Gao等提出了一種多載頻FDA雷達即各陣元發(fā)射頻率為步進多載頻信號的波束形成方法[7-8]，它們都能夠在目標位置處形成點狀單峰值波束，去除柵瓣，降低了副瓣，但形成的方向圖的距離維主瓣寬度較大，計算量大，并且均只適用于發(fā)射，也沒有得到其發(fā)射方向圖的閉合表達式，對于回波的接收也沒有展開深入研究。

考慮多輸入多輸出(Multiple-Input Multiple-Output, MIMO)雷達在目標探測、參數(shù)估計、空間分辨率等方面具有的許多優(yōu)勢[9]，本文結(jié)合FDA陣列與MIMO雷達，提出了一種基于子陣的頻率分集多輸入多輸出(Multiple-Input Multiple-Output with Frequency Diverse Array, FDA-MIMO)雷達去柵瓣的方法，推導出基于子陣的FDA-MIMO雷達天線方向圖閉合表達式，指出FDA-MIMO雷達采用子陣去除柵瓣和提高距離維分辨率的原理。仿真顯示，與采用線性頻率間隔的傳統(tǒng)方法相比，該方法不僅去除了柵瓣，也使得距離維主瓣寬度更窄，提高了距離維分辨率。

1 FDA-MIMO雷達信號模型

假設FDA-MIMO雷達是采用N個陣元的等距線陣，d為陣元間距，假設傳輸?shù)男盘柺沁h場窄帶信號，其與陣列法線方向夾角為θ，選定第一個陣元為參考陣元，從散射點到第一個陣元的距離為r，其陣列模型如圖1所示。

圖1 FDA-MIMO信號模型

第m個陣元通道的發(fā)射頻率為

fm=f0+(m-1)Δfm=1,2,…,N

(1)

式中,f0和Δf分別表示第一個陣元載波頻率和相鄰兩個陣元之間的頻率間隔。考慮為遠場窄帶信號，則第m個陣元與參考陣元的相位差為

(2)

由式(2)可以看出，F(xiàn)DA-MIMO雷達發(fā)射信號的相位差不僅與空間角有關(guān)，同時與傳播距離有關(guān)，因此，F(xiàn)DA-MIMO雷達天線方向圖具有距離-角度依賴性。注意到，式中第3項為陣元數(shù)的二次函數(shù)，當頻率間隔遠小于載頻，即NΔf?f0時，相比較前兩項帶來的相位差可忽略不計。

由于FDA-MIMO雷達不同陣元發(fā)射信號的包絡互相正交，則有

(3)

式中，φm(t)為第m個陣元發(fā)射信號的包絡，T為脈沖發(fā)射持續(xù)時間，( )*表示共軛計算，τ為任意的延遲時間。

第m個陣元接收的頻率為fn的信號可表示為

(4)

整個陣列接收的目標信號表示為

Y(θ,r)=ξ[y11,y12,…,y1N,y21,…，yNN]T=

ξb(θ)?a(θ,r,Δf)=ξv(θ,r,Δf)

(5)

式中，ξ為目標信號的幅度信息，上標T表示轉(zhuǎn)置，?為Kronecker積，v(θ,r,Δf)為陣列導向矢量，a(θ,r,Δf)為發(fā)射導向矢量，b(θ)為接收導向矢量，具體表達式如下[10]：

a(θ,r,Δf)=aθ(θ)⊙ar(r,Δf)=

(6)

(7)

假如有一目標位于(θs,rs)，為了讓FDA-MIMO雷達波束指向目標所在位置(θs,rs)，需要在FDA-MIMO雷達接收時進行波束加權(quán)，接收陣元通道加權(quán)矩陣為

(8)

式中，

(9)

(10)

此時，接收到的第m個陣元發(fā)射的第n個頻率的信號可以表示為

(11)

不考慮幅度信息，最后形成的方向圖可以表示為

Pr(t;θ,r)=wr(θs,rs)Hy(t;θ,r)=

(12)

式中，y(t;θ,r)為到達目標位置處的信號矢量。

由式(12)取極值可得FDA-MIMO雷達雙程方向圖的峰值點，有下式：

(13)

式中，k,p為任意整數(shù)。

(14)

由上式可知，當d=λ0/2時，F(xiàn)DA-MIMO雷達波束在角度上與頻率間隔Δf無關(guān)，沒有周期性，即在角度維沒有柵瓣，但在距離上具有周期性，即存在柵瓣，且柵瓣出現(xiàn)的距離周期為c0/2Δf，頻率間隔越大，距離周期越短，在探測范圍內(nèi)形成的柵瓣越多，特別是當干擾落在柵瓣時，就會使得輸出SINR急劇下降[11]。

根據(jù)柵瓣的距離周期性可以通過限制頻率間隔的選取范圍，使得FDA-MIMO雷達在探測范圍內(nèi)沒有柵瓣[12]，有

Δf≤c0/2Rmax

(15)

式中，Rmax為FDA-MIMO雷達探測的最遠距離。由上式可知，雖然可以解決距離維柵瓣的問題，但雷達探測越遠，頻率間隔選取的范圍越小，而距離維主瓣寬度為c0/NΔf[13]，頻率間隔越小，距離維主瓣寬度越大，距離維分辨率越低。

為解決這個矛盾，本文提出基于子陣的FDA-MIMO雷達去柵瓣的方法，在去除FDA-MIMO雷達柵瓣的前提下，提高了距離維分辨率。

2 基于子陣的FDA-MIMO雷達去柵瓣原理

如圖2所示，假設FDA-MIMO雷達陣列為N個陣元組成的線陣，陣元間距為半個波長，整個陣列分成k個子陣，前N1個陣元構(gòu)成子陣1，采用頻率間隔Δf1；依次類推，中間Ni個陣元構(gòu)成子陣i，采用的頻率間隔為Δfi，最后Nk個陣元構(gòu)成子陣k，采用的頻率間隔為Δfi，且N=N1+…+Ni+…+Nk。

圖2 FDA-MIMO雷達子陣模型

陣列帶寬為

(16)

假設目標位于(θs,rs)，此時，不考慮時間和幅度因素，則整個陣列的雙程方向圖為

|P(θ,r)|=

(17)

從上式很明顯可以看出，最后合成的方向圖是k個子陣方向圖的疊加。可通過使所有子陣方向圖的柵瓣位置盡量不重合來去除FDA-MIMO雷達大頻率間隔產(chǎn)生的柵瓣，同時增大頻率間隔來提高距離維分辨率。

以任意子陣i形成的方向圖為例進行分析，有其方向圖

|P(θ,r)|=

(18)

式(18)取極大值時，對應子陣i形成的方向圖柵瓣，有

(19)

式中，k,p為任意非零整數(shù)。參照式(13)分析可知，子陣i方向圖的柵瓣只能位于目標所在角度，距離維方向圖柵瓣位置為

(20)

同理可知，任意子陣j的距離維柵瓣位置為

(21)

式中，q為任意非零整數(shù)。要使這兩個子陣形成的柵瓣不能在同一個位置處疊加，有

(22)

化簡可得

(23)

(24)

此時，雖然子陣的柵瓣會有部分重合，但重合的位置位于FDA-MIMO雷達探測范圍外，即在探測范圍內(nèi)去除了柵瓣，不僅如此，由于采用的子陣頻率間隔較大，所以合成的方向圖距離維主瓣寬度更窄，距離維分辨率更高。

FDA-MIMO雷達距離維方向圖主要是由陣元個數(shù)與頻率間隔之積即帶寬Bp確定的，帶寬越大，距離維主瓣寬度越窄，距離維副瓣越低[13]。在FDA-MIMO雷達總陣元數(shù)一定時，若劃分的子陣數(shù)越多，為滿足子陣需采用互質(zhì)大頻率間隔的條件，會導致采取的頻率間隔變多、變大，整個陣列帶寬變大，合成的距離維副瓣電平會進一步降低。需要指出的是，由于子陣距離維副瓣電平本身較高，陣列合成的方向圖副瓣是子陣方向圖柵瓣與柵瓣或副瓣的疊加，所以整個陣列合成方向圖的副瓣會隨著子陣個數(shù)的增加降低到一定程度后，不會再降低。

3 仿真結(jié)果與分析

參數(shù)設置：FDA-MIMO雷達共有120個陣元，分成4個陣元數(shù)為30的子陣，假設探測最遠距離為200 km，目標位于(0°,100 km)。子陣頻率間隔依次為7，11，13，17 kHz；根據(jù)限定頻率間隔選取范圍去除柵瓣的方法(以下簡稱常規(guī)方法)，選取頻率間隔Δf=c/2Rmax=750 Hz。兩種方法的仿真圖分別如圖3和圖4所示。

(a) FDA-MIMO雷達方向圖

(b) 0°處切面(距離維)方向圖圖3 常規(guī)方法仿真圖

(a) FDA-MIMO雷達方向圖

(b) 0°處切面(距離維)方向圖圖4 基于子陣法仿真圖

從圖3和圖4可以發(fā)現(xiàn)，本文提出的基于子陣去柵瓣的方法與常規(guī)方法均可以在探測范圍內(nèi)去除柵瓣，形成指向目標位置處的點狀單峰值波束，而本文方法與常規(guī)方法相比，距離維主瓣寬度更窄，分辨率更高，但可以發(fā)現(xiàn)本文方法形成的波束副瓣在-12 dB左右，明顯比常規(guī)方法形成的副瓣高。

改變FDA-MIMO雷達陣元為120個，分成6個子陣，子陣頻率間隔改為7，11，13，17，19，23 kHz，其他參數(shù)不變，仿真圖如圖5所示。不改變參數(shù)設置，采用切比雪夫加窗(-30 dB)處理，兩種方法0°處切面(距離維)方向圖如圖6所示。

(a) FDA-MIMO雷達方向圖

(b) 0°處切面(距離維)方向圖圖5 基于子陣的去柵瓣方法仿真圖

(a) 常規(guī)方法

(b) 基于子陣法圖6 兩種方法的0°處切面(距離維)方向圖

對比圖4(b)和圖5(b)可以發(fā)現(xiàn)，本文提出的基于子陣的去柵瓣方法，副瓣會隨著子陣個數(shù)的增加逐漸降低，且降低為-15 dB，并且隨著子陣頻率間隔的增大，距離維主瓣寬度會減小，距離維分辨率會進一步提高。對比圖6(a)與圖3(b)以及圖5(b)與圖6(b)可以發(fā)現(xiàn)，在采用加窗處理后，本文方法的距離維副瓣電平幾乎沒有任何變化，而常規(guī)方法的距離維副瓣明顯下降。這主要是本文提出的基于子陣的去柵瓣方法形成的距離維副瓣，實質(zhì)上是子陣方向圖柵瓣與柵瓣或副瓣的疊加，采用加窗處理對柵瓣電平是沒有影響的，所以疊加后距離維副瓣電平幾乎沒有影響。

4 結(jié)束語

本文在推導出基于子陣的FDA-MIMO雷達雙程方向圖閉合表達式的基礎上，提出了一種基于子陣的FDA-MIMO雷達去柵瓣的方法，不僅去除了柵瓣，也提高了距離維主瓣分辨率。該方法提供了一種利用副瓣疊加去柵瓣的思路，不僅適用于FDA-MIMO雷達，也同樣適用于FDA雷達，然而通過仿真可以發(fā)現(xiàn)此方法形成的方向圖距離維副瓣仍然較高，采用加窗處理效果不明顯，需進一步研究降低副瓣電平的方法。