胡蜀徽,唐 佳,王志華,冉財章

(中國電子科技集團(tuán)公司第二十九研究所,四川成都610036)

0 引言

MIMO雷達(dá)各天線陣元發(fā)射相互正交的線性調(diào)頻信號,在空間疊加不能形成高增益波束,而是寬波束覆蓋[1],在接收端通過匹配濾波處理來恢復(fù)發(fā)射信號分量,并通過DBF來形成同時數(shù)字多波束,覆蓋發(fā)射波束區(qū)域。它具有抗截獲性能高、弱目標(biāo)檢測能力強(qiáng)、距離分辨率高和空間搜索效率高[2]等特點,其信號模型及處理方式類似于步進(jìn)頻率雷達(dá)。步進(jìn)頻率雷達(dá)在不增加單通道發(fā)射帶寬的前提下,通過脈沖間信號的合成來獲得距離上的高分辨。MIMO雷達(dá)在兼顧前者優(yōu)點的同時還能完成測速和相參積累功能。

步進(jìn)頻率雷達(dá)發(fā)射一串載頻跳變的線性調(diào)頻脈沖信號,然后再通過對脈沖回波進(jìn)行匹配濾波和IFFT處理獲得合成距離高分辨的效果。不同步進(jìn)頻率信號的參數(shù)選取會影響信號處理的方法,使得IFFT合成處理結(jié)果有一定冗余,從而造成目標(biāo)距離走動并產(chǎn)生虛假目標(biāo)。

MIMO雷達(dá)的處理方式是在數(shù)字下變頻和匹配濾波后將各個發(fā)射信號分量進(jìn)行相參合成,其信號處理過程包含了頻譜拼接和相位補(bǔ)償。合成后的一維距離像非常清晰,且副瓣很低,不會因為冗余造虛假目標(biāo),也不需要在距離上進(jìn)行裁剪拼接。

1 MIMO雷達(dá)的信號模型及距離維高分辨合成方法

1.1 MIMO雷達(dá)的信號模型

設(shè)發(fā)射脈沖串?dāng)?shù)目為N,天線個數(shù)為M,各天線發(fā)射頻率分量滿足正交關(guān)系[3],間隔為fΔ,子脈沖寬度為Tp,脈沖重復(fù)周期為Tr,調(diào)頻斜率μ=B/Tp,單個頻點帶寬為B,則接收到第i個重頻的線性調(diào)頻信號模型為

式中,rect(t/Tp)是寬度為Tp的矩形脈沖;f k是第k個信號分量的載頻,f k=f0+kfΔ,f0為中頻頻率,兩個不同頻率的信號分量之間呈正交關(guān)系[4]?？値挒锽s=B+(M-1)fp,信號的采樣率為

靜止目標(biāo)的雷達(dá)回波信號模型為

式中,τ=2R0/c,R0為該目標(biāo)離雷達(dá)的視線距離。

1.2 匹配濾波處理

經(jīng)過低通濾波器將全部頻譜分量提取出來,然后進(jìn)行數(shù)字下變頻,得到

不同頻率的匹配濾波信號模型為

將數(shù)字下變頻后的信號分別與不同頻率的匹配濾波信號進(jìn)行卷積,得到匹配濾波結(jié)果[5]。

式中,D=BTp為時寬帶寬積。

1.3 距離維高分辨合成處理

將不同頻率的匹配濾波的信號直接相參合成[6],得到的信號為

輸出信號的帶寬與匹配濾波器輸入信號的帶寬相等,為(M-1)fΔ+B。

利用等比數(shù)列求和公式：

得到合成后的信號包絡(luò)為

輸出信號包絡(luò)具有周期sinc函數(shù)的形式,輸出脈沖寬度近似為發(fā)射信號帶寬的倒數(shù),即T o=1/[(M-1)fΔ+B]。因此,合成后的最小可分辨距離為

如果取fΔ=B,則距離分辨率為ΔR=c/(2MB),相對于單個頻率的信號提高了M倍。

1.4 速度影響及消除方法

目標(biāo)相對于雷達(dá)的徑向速度是影響正交線性調(diào)頻信號合成效果的一個重要因素,會引起目標(biāo)峰值幅度降低,包絡(luò)走動,旁瓣幅度升高,虛假目標(biāo)出現(xiàn),從而導(dǎo)致距離維高分辨效果變差[7]。

設(shè)目標(biāo)的徑向速度為v,則雷達(dá)回波的信號模型為

式中,τi=2(R0-iv Tr)/c,R0為該目標(biāo)離雷達(dá)視線的距離,i為重頻序號,fr為射頻頻率,fd為多普勒頻率。

匹配濾波后信號的頻域表達(dá)式為

提取出對應(yīng)的相位項,可以得到

式中,φ1項為常數(shù)項,不會造成影響;φ2項為形成距離高分辨項的相位項;φ3項為各子脈沖間的多普勒項,在進(jìn)行距離合成時,會降低目標(biāo)的峰值,提高旁瓣的幅度,引起包絡(luò)的走動,當(dāng)運動參數(shù)較大時,將直接影響距離高分辨合成的效果;φ4項則會導(dǎo)致合成過后主瓣展寬。

因此,必須盡可能精確補(bǔ)償速度,才能保證合成的效果。步進(jìn)頻率測速方法有頻域互相關(guān)法、最小脈組誤差法、最小脈組相位差分法和正負(fù)頻率步進(jìn)信號組合法等[8],但大多計算量過大,測速精度不高。而MIMO雷達(dá)不需要精確測速,只需對每個正交信號的脈沖串進(jìn)行MTD處理,就可以補(bǔ)償速度帶來的相位項,然后再合成距離高分辨。

經(jīng)過數(shù)字下變頻和匹配濾波的結(jié)果為

當(dāng)積累時間不長的情況下,脈沖串之間的包絡(luò)走動可以忽略,則可以得到

對每一個頻率的脈沖信號進(jìn)行MTD處理,重排后對距離單元作FFT變換,可得

式中,Np為脈沖串個數(shù),多普勒頻率為fdk=2(fr+kfΔ)v/c。當(dāng)合成后的相對帶寬為窄帶時,可以認(rèn)為fdk相等。

當(dāng)fdk=fr/Npl時,獲得最大輸出,為

速度引起的相位項被補(bǔ)償,然后再進(jìn)行相參合成,就可以實現(xiàn)距離維的高分辨。

MIMO雷達(dá)距離高分辨合成的信號處理流程如圖1所示。

圖1 距離高分辨信號處理流程圖

2 仿真分析

根據(jù)上述理論分析,對基于4發(fā)4收的MIMO雷達(dá)進(jìn)行仿真。仿真參數(shù)如下：設(shè)正交線性調(diào)頻信號的起點頻率f0=10 MHz,射頻頻率fr=10 GHz,采樣頻率fs=100 MHz,調(diào)頻帶寬B=3 MHz,頻率間隔fΔ=3 MHz,脈沖寬度Tp=10μs,脈沖重復(fù)周期為100μs,脈沖個數(shù)為128。自相關(guān)結(jié)果如圖2所示,不同頻率信號之間的互相關(guān)結(jié)果如圖3所示。

從圖2可以得出,自相關(guān)的旁瓣電平為-13.7 dB。從圖3可以看出,相鄰兩個信號存在一定互相關(guān)分量[9],互相關(guān)峰值幅度大小如表1所示。圖3(a)、(d)及(f)中峰值約為-21 dB,其原因在于sinc函數(shù)旁瓣的影響[10]。

圖2 自相關(guān)結(jié)果

圖3 不同頻率信號間互相關(guān)結(jié)果

表1 信號互相關(guān)峰值幅度

設(shè)一目標(biāo)距離為16 km,回波信號為所有發(fā)射信號分量之和。對接收信號進(jìn)行匹配濾波的結(jié)果如圖4所示。從圖4可以看出,自相關(guān)旁瓣電平小于-13 dB,而互相關(guān)電平峰值小于-20 d B,互相關(guān)峰值的影響低于自相關(guān)峰值旁瓣,通過匹配濾波就可以分離出不同頻率的線性調(diào)頻信號。

設(shè)一目標(biāo)距離為16 km,另一目標(biāo)距離為16.035 km,合成處理前后的對比結(jié)果如圖5所示。根據(jù)仿真參數(shù),單個頻率的線性調(diào)頻信號的距離分辨率為50 m,合成后的距離分辨率為12.5 m,兩個目標(biāo)相差35m時,單個線性調(diào)頻信號的距離分辨率無法分辨,而合成后兩目標(biāo)凹點的幅度為-16.94 dB,低于旁瓣峰值-13.7 dB,能夠清晰地分辨出兩個目標(biāo)。

圖4 不同頻率的匹配濾波輸出

圖5 合成前后匹配濾波器輸出

不同速度條件下的合成處理結(jié)果如圖6所示,隨著速度的增加,旁瓣峰值升高,目標(biāo)峰值降低。當(dāng)速度從200 m/s增加到800 m/s時,旁瓣峰值從-12.38 dB上升到-4.93 dB,目標(biāo)峰值從-1.22 dB下降到-5.88 dB,目標(biāo)已經(jīng)不能正確檢測,而且還會出現(xiàn)虛假目標(biāo),距離維上目標(biāo)峰值所在位置也出現(xiàn)了走動。

圖6 速度對合成效果的影響

對不同速度條件下的線性調(diào)頻脈沖串進(jìn)行MTD處理,補(bǔ)償速度引起的相位項,然后再進(jìn)行相參合成處理,距離高分辨的合成結(jié)果如圖7所示?？梢钥闯?速度引起的相位項的影響已經(jīng)消除,不同速度條件下旁瓣和目標(biāo)峰值與速度為0時一致,能夠正確檢測出兩個目標(biāo)。

3 結(jié)束語

首先給出了MIMO雷達(dá)的信號模型和處理方法,分析了OFD-LFM信號的自相關(guān)函數(shù)和互相關(guān)函數(shù),其特點是利用頻率之間的正交特性,通過匹配濾波分離出不同頻率的信號分量,然后相參合成實現(xiàn)距離維的高分辨。同時還分析了速度對合成效果的影響,以及通過MTD處理補(bǔ)償速度引起的相位項消除速度影響的方法。該合成方法計算量較小、易實現(xiàn),相比步進(jìn)頻率雷達(dá)需在距離上去冗余和測得速度后進(jìn)行速度補(bǔ)償?shù)奶幚矸绞礁邇?yōu)勢,有利于工程應(yīng)用。

圖7 MTD后的合成效果

[1]FISHLER E,HAIMOVICH A,BLUM R,et al.MIMO Radar：an Idea Whose Time Has Come[C]∥Pro-ceedings of the IEEE International Conference on Radar,Philadelphia：IEEE Signal Processing Society,2004：71-78.

[2]何子述,韓春林,劉波.MIMO雷達(dá)概念及其技術(shù)特點分析[J].電子學(xué)報,2005,33(12A)：2441-2445.

[3]李宏偉.MIMO雷達(dá)波形設(shè)計方法綜述[J].現(xiàn)代雷達(dá),2013,35(6)：12-14.

[4]陳正輝,嚴(yán)濟(jì)鴻,何子述.MIMO雷達(dá)OFDM-LFM波形設(shè)計與實現(xiàn)[J].雷達(dá)科學(xué)與技術(shù),2013,11(1)：77-81,86.CHEN Zheng-hui,YAN Ji-hong,HE Zi-shu.Design for MIMO Radar and Implementation of OFDM-LFM Waveform[J].Radar Science and Technology,2013,11(1)：77-81,86.(in Chinese)

[5]陳正中,李小波,梁浩,等.正交頻分LFM信號的MIMO雷達(dá)的匹配濾波技術(shù)[J].雷達(dá)科學(xué)與技術(shù),2013,11(2)：197-202.CHEN Zheng-zhong,LI Xiao-bo,LIANG Hao,et al.Matched Filtering Technology for MIMO Radar with OFD-LFM Signal[J].Radar Science and Technology,2013,11(2)：197-202.(in Chinese)

[6]裴英,肖文書.相干MIMO雷達(dá)關(guān)鍵技術(shù)研究[J].現(xiàn)代雷達(dá),2013,35(3)：9-12.

[7]丁海林,李亞超,高昭昭,等.線性調(diào)頻步進(jìn)信號的三種合成方法的對比與分析[J].火控雷達(dá)技術(shù),2007,36(4)：10-16.

[8]曹宇飛,屈曉光,黃培康.一種快速有效的步進(jìn)頻率測速方法[J].現(xiàn)代防御技術(shù),2009,37(1)：80-82.

[9]劉波,何子述,王海江.MIMO雷達(dá)中的DFCW設(shè)計及性能分析[J].電子科技大學(xué)學(xué)報,2010,39(5)：688-691.

[10]施祥同,王虎,陳建軍,等.OFDM雷達(dá)信號的寬帶模糊函數(shù)性能分析[J].雷達(dá)科學(xué)與技術(shù),2010,8(6)：554-558.SHI Xiang-tong,WANG Hu,CHEN Jian-jun,et al.Wideband Ambiguity Function of OFDM Radar Signal[J].Radar Science and Technology,2010,8(6)：554-558.(in Chinese)