李朝海，李 明，張 偉

(電子科技大學電子工程學院 成都 611731)

.通信與信息工程.

一種MIMO雷達更窄主瓣波束形成方法

李朝海，李 明，張 偉

(電子科技大學電子工程學院 成都 611731)

介紹了一種獲得更窄主瓣虛擬收發(fā)波束的方法。該方法通過天線陣元發(fā)射正交編碼的波形，然后在接收端對回波信號進行處理。MIMO雷達由于靠稀疏陣列避免旁瓣的不利影響，所以可提高角分辨率和多普勒分辨率。但受到如空間、成本或者其他實際因素的限制，須將陣列布置成緊湊陣列，而緊湊陣列對低速目標檢測效果較差。通過提出的加窗方式，緊湊陣列利用虛擬收發(fā)波束形成可獲得更窄的主波束，提高對低速目標的檢測效果。通過仿真驗證，對歸一化頻率為0.062 5的低速目標，信干噪比提升系數(shù)(SINR IF)約提升7.3 dB。

緊湊陣列; MIMO雷達; 稀疏陣列; 加窗方式

MIMO雷達相比于傳統(tǒng)雷達有許多優(yōu)勢[1-3]，如具有更高的角分辨率和多普勒分辨率等。當MIMO雷達發(fā)射正交信號[4-7]，發(fā)射模式和接收模式均可進行陣列處理，在接收端通過匹配濾波進行信號提取，提取到的信號可用于提高雜波的空間分辨率。利用不同發(fā)射天線和不同接收天線的相位偏差可構成一個新的虛擬陣列導向矢量，合理設計天線位置，MIMO雷達利用少量的天線就可獲得一個非常長的嚴格采樣的導向矢量數(shù)組[8]。

MIMO雷達使用稀疏陣列可避免旁瓣帶來的不利影響[9]。但設置成稀疏陣列需要更大空間，制約了在機載雷達上的應用，且發(fā)射和接收不能使用相同天線，增加了成本。若使用緊湊陣列則可以很好地解決這些問題。本文利用新的算法使用緊湊陣列獲得比傳統(tǒng)波束形成方案更窄的主波束。

1 MIMO雷達信號模型

假設有TN個發(fā)射單元，間隔均為Td，有RN個接收單元，間隔均為Rd。采用單站式MIMO雷達，在每個發(fā)射單元，一個相關處理間隔發(fā)射脈沖的個數(shù)為L個。假設發(fā)射信號相互正交，系統(tǒng)內部沒有雜波干擾，雷達系統(tǒng)平臺的運動速度為v[10]。發(fā)射天線和接收天線有相同的方位角θ和俯仰角?。設為波長。在每個脈波重復間隔(PRI)rT，系統(tǒng)發(fā)射波束集合為：

式中，γ=dT/dR；fd=2vsinθcos??Tr/λ是移動平臺的歸一化多普勒頻率；fs=dRsinθcos?/λ是空間頻率；ξ(θ)是在θ處的反射系數(shù)。假設dR設置合理，保證?0.5≤fs≤0.5，從而避免空間頻率混疊。

本文假設發(fā)射導向矢量和接收導向矢量分別為：

定義多普勒矢量為：

對于第l個脈沖，接收單元的雜波回聲為：

式中，Yl∈CNR×1。接收端信號處理是為了分離不同信號的相對時間延遲特性，從而為后續(xù)處理提供信號向量。每個接收到的信號包含所有發(fā)射波形，設sH為接收端的濾波器組，有：

式中，X∈CNR×NT。使用向量形式保持第l個脈沖的過濾雜波數(shù)據(jù)：

式中，x∈CNRNT×1；“?”表示Kronecker積。距離環(huán)l被分為NC(NC>>NTNRL)個雜波片，那么第l個脈沖的雜波數(shù)據(jù)分離形式為：

然后將第L個脈沖的整個雜波回波保存為向量形式：

2 更窄主瓣波束形成方法

接收波束通過相位變換疊加處理來利用天線陣元的輸出實現(xiàn)延遲疊加波束形成，從而實現(xiàn)數(shù)字合成[12]。MIMO雷達的正交信號可在接收端通過匹配濾波進行提取，提高雜波的空間分辨率[13]。維數(shù)為NTNR的陣列導向矢量可以等效為發(fā)射導向矢量與接收導向矢量的卷積。假設NR=NT=N ，對于波束形成，只考慮等間距的線性陣列。已經證明在匹配濾波之后實際獲得的天線波束相比于對應的傳統(tǒng)相控雷達具有相同的天線增益和空間分辨率。同時數(shù)字濾波器的系數(shù)矢量按如下方式設置：

式中，κ滿足歸一化的要求；s0是來自期望方向θ0處的目標信號，有：

通過式(12)和式(13)，可以獲得具有最大信噪比(SNR)的濾波器輸出。而收發(fā)波束的指向可通過數(shù)字濾波器獲得。將式(12)改寫為：

式中，κ=κ1κ2，均為常數(shù)；wT為發(fā)射波束的數(shù)字濾波器的系數(shù)向量；wR為接收波束的權向量。利用式(12)和式(14)中的最佳濾波器，通過虛擬收發(fā)波束形成，可獲得虛擬天線增益：

上式為離散傅里葉變換(DFT)[14]的形式。本文設置dR=dT=d，d/λ=0.5，NR=NT=N ,可得：

從式(15)和式(16)知，G(θ,θ0)可以看做是兩個長度均為N的方波離散傅里葉變換(DFT)的積。圖1闡述了一維均勻線性陣列導向矢量的卷積性質。

MIMO雷達的虛擬收發(fā)波束形成可以看作是兩個方波的離散傅里葉變換的積。兩個方波的卷積是三角波，其長度為2N?1，如圖1所示。圖1中的振幅是具有相同的相位中心的陣元的數(shù)量。根據(jù)傅里葉理論，發(fā)射波束的主瓣寬度被擴寬了。為了更明確這一點，改寫MIMO雷達接收端經過匹配濾波后的信號的矩陣如下：

式中，γ=dT/dR，vT為發(fā)射導向矢量；vR為接收導向矢量，它們均已在式(3)和式(4)中進行了定義。如果考慮到噪聲，則匹配濾波后的輸出信號為：

可以發(fā)現(xiàn)式(19)中有很多元素相等，如圖2所示。

根據(jù)離散傅里葉變換的性質，相同長度的方波主瓣比三角波更窄。在直覺上可以通過權值反轉變換，將三角波變換為方波，如圖3所示。

與式(18)或式(19)對應的權值轉換矩陣為：

d=stack(D)。對應于虛擬收發(fā)波束的數(shù)字濾波器系數(shù)向量為：

如果使用窗函數(shù)，由于權向量wnew與期望信號s不能完美匹配，所以SNR有一定的損失。根據(jù)式(18)中的一維波束形成加窗的信噪比損失因子，定義虛擬收發(fā)波束加窗的信噪比損失因子為ηeff，它是加窗后的信噪比(SNRtap)與最優(yōu)輸出信噪比(SNRopt)之比值。為便于計算，ηeff通常定義為：

式中，s0是期望信號，其定義如式(13)?？紤]錐形的情況，wnew=s0*d，保證單位增益約束為wHw=1，vec(t)必須滿足：

圖4為信噪比損失因子ηeff與陣元數(shù)目變化關系。從圖4可以發(fā)現(xiàn)，隨著陣元數(shù)目的增加，ηeff的減小速率將會降低。對應于16陣元、32陣元、80陣元，其ηeff分別為?2.515、?3.193、?3.998 dB。更窄的主瓣會帶來SNR的損失，可以根據(jù)需要合理選擇。更窄的波束可以提高角分辨率和多普勒分辨率，抑制來自目標的主瓣雜波干擾。通常，接收端的剩余旁瓣會帶來足夠的干擾掩蓋期望的返回信號，在波束形成中是不允許存在過高的旁瓣。本文也可以使用切比雪夫窗函數(shù)。直接將切比雪夫窗函數(shù)式(21)寫為：

將式(25)對應的加窗方式命名為方式1。根據(jù)式(20)中提出的加權方式，該加權方式將會改變切比雪夫窗函數(shù)的形狀，旁瓣可能得不到有效地抑制?？梢酝ㄟ^改變切比雪夫窗函數(shù)來克服該問題。根據(jù)圖4描述的相位中心，將wch,2N?1擴展到N2個點，wch,2N?1是2N?1個點的切比雪夫窗函數(shù)的權向量。即由相同相位的匹配濾波器有相同的切比雪夫窗權系數(shù)為：

將式(26)對應的加窗方式命名為方式2，將式(26)代入式(25)，替換

可以得到虛擬天線增益為：

為了評估減小的主瓣雜波，虛擬收發(fā)波束形成被配置成串聯(lián)T-S模式。時間權向量首先在每個天線通道的每個匹配濾波器的輸出上進行計算，有：

式中，wt,l是對應于離散傅里葉變換系數(shù)的權向量；下標l表示第l個多普勒通道；下標t表示時域。多普勒處理器的輸出結果在空域中進行進一步處理以獲得檢驗統(tǒng)計量：

式中，下標s表示時域，有：

通過漸變模式得到的更窄的主瓣波束有很多優(yōu)勢，比如提升角分辨率和多普勒分辨率，具有更尖銳的機載雷達雜波凹口。這些優(yōu)勢都可以應用于機載雷達對地面低速運動目標的檢測。通常用單陣元、單脈沖的SINR值來衡量空時處理算法的性能。SINRin作為單陣元、單脈沖的SINR。通常情況下雜波和干擾比較強，所以SINRin相對是一個非常小的量。SINR的提高因子IFSINR被定義為：

SINRout(?)在串聯(lián)T-S模式下的配置為：

式中，st,l(fd)是期望目標的多普勒導向矢量：

Rl是第l個多普勒通道的協(xié)方差矩陣，有：式中，I是單位矩陣；R是整個協(xié)方差矩陣，有：

3 仿真結果

考慮一八陣元發(fā)射接收陣列，NT=NR= N=8，dT=dR=d=λ/2，主瓣方向θ0=0°，目標位于主瓣所在方向。載波頻率為f0=1.2 GHz ，脈沖重復頻率PRF為2 000 Hz，一個相關處理間隔發(fā)射脈沖數(shù)為M=64。雷達運動速度為v=125 m/s ，對于每一個天線陣元，每一個脈沖內的SNR=?5 dB，CNR=45 dB。式(25)對應方式1的波束形成方式，式(27)中對應方式2的波束形成方式，以及只考慮在接收波束形成中使用30 dB切比雪夫加窗(方式3)的相控陣雷達波束。對比3種不同的波束形式如圖5所示。

由圖5可知，由于方式1的切比雪夫窗函數(shù)形式的改變并沒能有效抑制旁瓣。而方式2的切比雪夫加窗過程由于采用虛擬的發(fā)射接收波束，旁瓣得到了有效的抑制。對于相控陣雷達系統(tǒng)，由于考慮到最大發(fā)射功率的要求，切比雪夫加窗只用在接收波束形成過程中，雖然旁瓣得到了抑制，但主瓣也相應的展寬了。相比于其他方法的信干噪比改善因子如圖6所示。為方便比較，在歸一化多普勒頻率處的最小的改善因子也被歸一化為0 dB。

圖6中所加切比雪夫窗為30 dB，本文提出的基于漸變處理的波束形成方法主要在小多普勒頻率范圍優(yōu)于沒有進行漸變處理的波束形成方法，這一點對于檢測低速運動目標非常有用。比如對于低速運動目標，在歸一化多普勒頻率為0.062 5時，新方法的信干噪比改善因子可以達到7.3 dB。沒有進行切比雪夫加窗時，新方法的信干噪比改善因子比沒有進行漸變處理時在旁瓣區(qū)域差，但可以通過簡單地修正切比雪夫窗函數(shù)來克服不足。在小多普勒頻率區(qū)域，用30 dB切比雪夫加窗形成的相控陣雷達接收波束效果比較差，特別是在歸一化多普勒頻率位于(?0.1,0.1)時。通過上面的分析可知，波束主瓣寬度越窄，對低速目標的檢測效果越好。

4 總 結

本文提出的MIMO雷達匹配濾波器輸出信號的幅度經過加窗處理之后，對低速運動物體的檢測效果有一定程度的提高。通過本文所提出的加窗方式可以獲得更窄的主瓣，顯著提高雜波抑制效果，改善角度分辨率，提升探測低速運動目標性能，因此可以提高對低速運動物體的檢測效果。對于其他陣元數(shù)目的MIMO雷達，采用該加窗方式之后檢測效果類似。方便將該方式應用于級聯(lián)時空處理。本文的方法也可以應用于其他波束域的STAP算法。這種新型加窗方式提高了對低速運動目標的檢測效果，完善了MIMO雷達的性能，對于現(xiàn)代機載雷達具有重要意義。

[1] BROOKNER E, RAYTHEON C. MIMO radar demystified and where it makes sense to use[C]//Phased Array Systems&Technology. Waltham, USA: IEEE, 2013.

[2] MOO P W, ZHEN Ding. Tracking performance of MIMO radar for accelerating targets[J]. Signal Processing, 2013, 61(21): 5205-5216.

[3] RANKIN G A, TIRKEL A Z, BUI L Q, et al. Radar imaging: Conventional and MIMO[C]//Communications and Electronic(ICCE). Hue, USA: IEEE, 2012.

[4] HE Qian, LEHMANN N H, BLUM R S, et al. MIMO radar moving target detection in homogeneous clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(3): 1290-1301.

[5] 牛朝陽, 張劍云. MIMO雷達正交波形集設計—線性調頻-位編碼混合波形[J]. 計算機工程與應用, 2012, 48(13): 133-137. NIU Zhao-yang, ZHANG Jian-yun. Orthogonal waveform set of hybird LFM and phase code for MIMO radar[J]. Computer Engineering and Application, 2012, 48(13): 133-137.

[6] RABIDEAU D J. Non-adaptive multiple-input, multipleoutput radar techniques for reducing clutter[J]. IET Radar Sonar Navig, 2009, 3(4): 304-313.

[7] 杜曉林, 王旭, 蘇濤, 等. MIMO雷達正交相位編碼信號的代數(shù)設計方法[J]. 西安電子科技大學學報(自然科學版), 2014, 41(2): 64-70. DU Xiao-lin, WANG Xu, SU Tao, et al. Orthogonal phase coding signal design for MIMO radar via an algebraic method[J]. Journal of Xidian University(Natural Science Edition), 2014, 41(2): 64-70.

[8] CHEN Chun-yang, VAIDYANATHAN P P. Beamforming issues in modern MIMO radars with doppler[C]//The 40th Asilomar Coference on Signals, Systems and Computers. Pacific Grove, CA, USA: IEEE, 2006.

[9] FORSYTHE K W, BLISS D W. MIMO radar signal processing[M]. New York, USA: Wiley, 2009.

[10] DENG Hai, BRAHAM H. A virtual antenna beamforming (VAB) approach for radar systems by using orthogonal coding waveforms[J]. IEEE Transactions on Antennas and Propagation, 2009, 57(2): 425-435.

[11] WANG Guo-hua, LU Yi-long. Clutter rank of mimo radar with a speical class of waveforms[C]//Waveform Diversity and Design Conference. Kissimmee, FL, USA: IEEE, 2009.

[12] MONDAL B, HEATH R W. Performance analysis of quantized beamforming MIMO systems[J]. IEEE Trans Signal Process, 2006, 54(12): 4753-4766.

[13] CHEN Chun-yang, VAIDYANATHAN P P. MIMO radar space-time adaptive processing using prolate spheroidal wave functions[J]. IEEE Transactions on Signal Processing, 2008, 56(2): 623-635.

[14] 周劍雄, 石志廣, 吳京. 離散時間信號頻域分析中的概念辨析[J]. 電氣電子教學學報, 2014, 36(1): 49-51. ZHOU Jian-xiong, SHI Zhi-guang, WU Jing. Conceptual discrimination in frequency domain analysis of discrete time signals[J]. Journal of EEE, 2014, 36(1): 49-51.

編 輯 稅 紅

An Approach to Achieving Narrower Mainlobe for MIMO Radar

LI Chao-hai, LI Ming, and ZHANG Wei
(School of Electronic Engineering, University of Electronic Science and Technology of China Chengdu 611731)

An innovative approach is introduced to form narrower mainlobe of virtual transmitting-receiving beamforming by transmitting orthogonal coding waveforms from the antenna elements and digitally processing their echoes at the receiver. The benefits such as increasing angle resolution and increasing Doppler resolution of MIMO radar mainly depend on the sparse arrays without the adverse effects of sidelobes. However, due to the actual limitation, such as room limitation for airborne, cost or other factors, we have to set the array in dense arrays configuration which is not good at detecting slow moving target. With the proposed amplitude taper approach, a narrower main beam can be achieved with dense arrays configuration after virtual transmitting-receiving beamforming and therefore the performance of detecting slow moving target is improved. Simulation result shows that for the normalized Doppler frequency 0.062 5, i.e., slow moving target, the signal to interference and noise improvement factor (SINR IF) is increased about 7.3 dB.

dense arrays; MIMO radar; sparse arrays; tapering

TN911.3

A

10.3969/j.issn.1001-0548.2015.03.001

2014 ? 04 ? 25；

2015 ? 01 ? 06

國家自然科學基金(61371184, 61201280, 61401062, 61301262)

李朝海(1972? )，男，高級工程師，主要從事相控陣雷達與自適應信號處理方面的研究.