劉 尚,蔣金華,段海洋,杜飛飛

(1.江南機(jī)電設(shè)計研究所,貴陽 550009;2.西北工業(yè)大學(xué),西安 710072)

0 引 言

傳統(tǒng)的多目標(biāo)方位估計(Direction of Arrival,DOA)方法采用波束形成和最小方差無畸變響應(yīng)(Minimum Variance Distortionless Response,MVDR)方法以及一系列如MUSIC和ESPRIT等高分辨率估計方法。這些方法利用均勻線性陣列進(jìn)行多目標(biāo)方位估計,難以解決因空時欠采樣帶來的陣元間相位差模糊的問題,部分信號方位角估計效果不理想。

壓縮感知(Compressed Sensing,CS)理論是現(xiàn)代信息科學(xué)領(lǐng)域中一個全新的研究方向[1-2],對雷達(dá)信號處理、聲納成像[3]、模式識別、地質(zhì)勘探、信道估計等都有重要意義,它提供了一種新的信號采集和參數(shù)估計方法,能夠以遠(yuǎn)少于經(jīng)典奈奎斯特(Nyquist)采樣理論所需的測量數(shù)據(jù)精確恢復(fù)滿足稀疏性的信號或估計其相應(yīng)的參數(shù)[4]。信號在時域和空域的稀疏性特征滿足信號稀疏重建理論的基本條件,利用隨機(jī)采樣矩陣可以實現(xiàn)用較低采樣率的時域和空域數(shù)據(jù)恢復(fù)及估計信號的效果。相關(guān)學(xué)者從信號處理方法和陣元幾何設(shè)計方面開展了研究:文獻(xiàn)[5]提出基于時延、FFT和ESPRIT算法的高頻段空間2D到達(dá)方向無模糊估計方法,但該方法估計信號數(shù)目受陣元數(shù)目限制;文獻(xiàn)[6]基于窄帶信號設(shè)計非均勻陣列,且陣元布放在均勻劃分的以半波長為單位的網(wǎng)格點上,沒有解決陣元間距小于半波長的問題;文獻(xiàn)[7]提出利用CS理論對導(dǎo)向向量泰勒分解后,使用交替下降的最小二乘法實現(xiàn)信號重建,但此法利用混合范數(shù)的約束條件,重建過程中計算量較大,算法精度也受限。

為了解決上述空時欠采樣條件下多目標(biāo)DOA存在的問題,本文基于CS理論,采用空間網(wǎng)格角度劃分形式實現(xiàn)信號在空間域的稀疏性表示,在此稀疏域下將實現(xiàn)信號與其空間位置的一一對應(yīng),對陣元接收信號在測量矩陣上進(jìn)行壓縮投影采樣,通過改進(jìn)MVDR空間功率譜估計中的壓縮測量陣列接收信號的空間自相關(guān)矩陣和陣列流行向量的稀疏值,實現(xiàn)多目標(biāo)方位角估計,最后通過仿真驗證了本文方法的有效性。

1 多目標(biāo)信號方位估計模型

信號方位估計是根據(jù)信號來波估計信號方位角的信號處理技術(shù),是在空域濾波和空間譜估計基礎(chǔ)上發(fā)展起來的[8]?？沼驗V波及空間譜估計理論框架如圖1、圖2所示。

用ωi表示加權(quán)系數(shù),則相應(yīng)的濾波器輸出可表示為

(1)

圖3給出了K個遠(yuǎn)場信號入射到某空間的均勻線列陣模型,其中陣元數(shù)為M。

圖3 均勻線列陣模型

以均勻線列陣為例,陣列所接收信號的時延差為

(2)

兩個陣元間的相位差為

(3)

式中,d為兩個陣元間距;c為信號在空氣中的傳播速度;θ為遠(yuǎn)場信號的入射角度;λ為信號波長。

可以看出,只要知道信號的相位延遲,就可以得到信號的來波方向,進(jìn)而實現(xiàn)信號方位估計。

2 信號空域稀疏表示

均勻線列陣在某時刻t的接收信號的數(shù)學(xué)模型為

x(t)=As(t)+n(t)

(4)

式中,A為M×K維陣列流形矩陣,且

A=[a(θ1),a(θ2),a(θ3),…,a(θK)]

(5)

陣列流形向量為

(6)

式中,θi為第i個信號的入射方位角。

目標(biāo)源的輸入信號矢量為

(7)

因此陣元接收信號可以表示為

=[a(θ1),a(θ2),a(θ3),…,a(θk)]S(t)+N(t)

=AS(t)+N(t)

(8)

由式(8)可知,陣列流行矩陣的每一列a(θi)都對應(yīng)一個目標(biāo)信號源的空間位置,即矩陣A包含了目標(biāo)的方位信息。陣列流行矩陣的構(gòu)建也決定了信號在空間域的稀疏化表示方式。

將空間域劃分為{θ1,θ2,…,θNs},并假設(shè)每一個可能的位置θi,i=1,2,…,Ns都對應(yīng)一個潛在的目標(biāo)源信號si,如圖4所示,其中“·”表示空間中實際存在的信號,○表示實際不存在的信號,矩陣A的每一列對應(yīng)一個潛在目標(biāo)信號源的方位信息。為了體現(xiàn)目標(biāo)信號源的稀疏性,潛在的信號源個數(shù)比實際要多得多,即Ns?K,這樣就構(gòu)造了一個Ns×N維的稀疏信號s=[s1,s2,…,sNs]Τ,N表示快拍數(shù)。在目標(biāo)源信號s中只有實際存在目標(biāo)的K個位置有非零元素,而Ns-K個位置均為零。

圖4 信號空間稀疏化過程

對于給定的空間域稀疏化方式{θ1,θ2,…,θNs},式(5)中的陣列流行矩陣A得以確定,基于CS理論的DOA估計問題就轉(zhuǎn)化為利用已知陣列接收數(shù)據(jù)x和陣列流行矩陣A來重構(gòu)空間稀疏信號s=[s1,s2,…,sNs]T,該重構(gòu)過程將一個凸優(yōu)化問題轉(zhuǎn)化為簡單的線性規(guī)劃問題來求解[9],其中前K個最大的重構(gòu)分量就是空間上實際存在的信號源的重構(gòu)信號。根據(jù){θ1,θ2,…,θNs}與s=[s1,s2,…,sNs]Τ的一一對應(yīng)關(guān)系就能得到信號的DOA估計。

3 CS-MVDR方法

MVDR[10]方法使一些自由度在期望觀測方向形成波束,并利用剩余的自由度在干擾方向形成零陷,使非期望干擾最小,同時使增益在觀察方向上保持為常數(shù)1,在減小非期望干擾的同時構(gòu)建一個約束最優(yōu)化問題:

(9)

利用拉格朗日乘數(shù)法求解式(9),令目標(biāo)函數(shù)

(10)

對目標(biāo)函數(shù)中的w求導(dǎo),令導(dǎo)函數(shù)為零,得到最優(yōu)的權(quán)矢量為

w=λR-1a(θ0)

(11)

再利用wHa(θ0)=1可以求得常數(shù)

(12)

將式(12)代入式(11),得到最優(yōu)的權(quán)矢量表達(dá)式為

(13)

陣列輸出功率用空間譜描述為

(14)

根據(jù)壓縮感知信號空域稀疏表示理論,陣列流行矩陣A相當(dāng)于壓縮感知理論中的稀疏矩陣Ψ,因此為了保證能夠精確重構(gòu)信號源,并得到方位角的準(zhǔn)確估計,A必須滿足一定的稀疏重構(gòu)條件[11]。對于Ψ,當(dāng)陣元數(shù)M小于稀疏系數(shù)矢量的長度Ns時,A相當(dāng)于一個過完備基,即A是空域稀疏基。

設(shè)Φ為m×M,m

由于隨機(jī)稀疏測量矩陣滿足RIP準(zhǔn)則[12],因此利用壓縮測量數(shù)據(jù)y和感知矩陣Θ即可大概率重構(gòu)出陣列輸入信號s=[s1,s2,…,sNs]Τ[13-14]。DOA估計問題就是求解Ns×N維的稀疏信號s,其中的非零元素及其在向量中的位置信息分別表征目標(biāo)的幅度信息和角度信息,方位角估計就是在保證s在信號源真實方向相關(guān)的角度形成很明顯的峰值,而在其他區(qū)域是近似為零的低電平,尖峰數(shù)量代表信號源數(shù)量。

因此,根據(jù)式(14),基于壓縮感知改進(jìn)的CS-MVDR陣列輸出功率為

式中,b(θ)=Φa(θ)為壓縮測量陣列流行向量;Ry為壓縮測量陣列接收信號的空間自相關(guān)矩陣;a(θ)為傳統(tǒng)陣列流行向量。

4 仿真及分析

對非均勻線列陣進(jìn)行空時域欠采樣,利用Matlab仿真比較MVDR和CS-MVDR在低信噪比、低陣元個數(shù)情況下的多目標(biāo)方位角估計精度。

(1)仿真1:時域欠采樣下的空間譜圖

對傳統(tǒng)陣列進(jìn)行時域欠采樣、空域正常采樣。假設(shè)信號源個數(shù)為3,信號方位角分別為-60°、40°、50°,陣元間距選擇半波長d=λ/2,采樣頻率選擇小于3個信號源中心頻率中最大頻率的2倍(fs<2*fmax(1、2、3)),傳統(tǒng)線列陣的陣元個數(shù)分別為8、16、32,壓縮感知陣列陣元數(shù)為6、8,利用32陣列隨機(jī)抽取得到。3個目標(biāo)信號方位角的角度掃描范圍為[-90°,90°],即在角度網(wǎng)格劃分為N1=181。

圖5、圖6分別為信噪比0 dB、10 dB下利用壓縮感知的6陣元和8陣元進(jìn)行波達(dá)方向估計結(jié)果,采用隨機(jī)稀疏測量矩陣,分別與8、16、32陣元下的MVDR估計結(jié)果作比較。

圖5 時域欠采樣下SNR=0 dB時的空間譜圖

可以看出,8陣元MVDR估計在0 dB和10 dB下幾乎都不能分辨臨近的40°和50°目標(biāo),而6、8陣元CS-MVDR和16、32陣元MVDR均能正確估計出3個方位角度數(shù)。CS-MVDR方法在保證正確估計方位角精度的前提下減少了陣元數(shù)目,大大降低了計算量。相較于6陣元CS-MVDR,8陣元旁瓣起伏較小,方位估計穩(wěn)健性更強。

(2)仿真2:空時欠采樣下的空間譜圖

對傳統(tǒng)陣列進(jìn)行時域、空域欠采樣。信號源個數(shù)為3,信號方位角分別為-60°、40°、50°,陣元間距選擇為半波長的1.2倍即d=1.2×λ/2,采樣頻率選擇小于3個信號源中心頻率中最大頻率的2倍。信噪比為0 dB和10 dB時的仿真結(jié)果分別如圖7、圖8所示。

圖7 空時欠采樣下SNR=0 dB時的空間譜圖

圖8 空時欠采樣下SNR=10 dB時的空間譜圖

可以看出,8、16、32陣元MVDR估計在0 dB和10 dB下都不能正確估計3個目標(biāo)方位角,而10 dB下基于壓縮感知的6、8陣元陣列能正確估計出3個方位角度數(shù),但是在0 dB下6陣元CS-MVDR法旁瓣起伏較大,甚至超過主瓣高度,8陣元的旁瓣起伏相對較小。

(3)仿真3:CS-MVDR法同時估計多個不同方位角

依然對3個目標(biāo)信號源進(jìn)行空時欠采樣,信號方位角為45°、55°、65°。陣元間距選擇為半波長的1.2倍即d=1.2×λ/2。信噪比從-20 dB步進(jìn)到10 dB,步進(jìn)長度為2 dB,快拍數(shù)為300,蒙特卡洛仿真次數(shù)為500,統(tǒng)計分析不同信噪比下多目標(biāo)方位角估計的檢測概率(當(dāng)估計均值的絕對誤差小于1°時,該次檢測有效),結(jié)果如圖9～11所示。

圖9 信噪比與檢測概率關(guān)系曲線

圖10 信噪比與信號方位估計均值關(guān)系曲線

圖11 信噪比與信號方位估計均方差關(guān)系曲線

可以看出,在低信噪比條件下,由于噪聲對空時欠采樣CS-MVDR算法的影響,檢測概率相對較低,估計的信號方位角和真實值誤差較大。不過,當(dāng)SNR≥12 dB時,檢測概率都能達(dá)到100%;當(dāng)SNR≥-14 dB時,估值接近真實值,且估計的均方誤差為零。在空時欠采樣下,當(dāng)同時估計多個信號方位角時,在誤差收斂前不同角度間誤差相差較大,這與目標(biāo)入射角度來向和陣列流型矩陣A中對應(yīng)的目標(biāo)信號源方位信息有關(guān),將在后續(xù)研究中進(jìn)行深入分析。

5 結(jié)束語

針對傳統(tǒng)MVDR方法在空域上采集測量數(shù)據(jù)量有限及對噪聲背景相對敏感等問題,本文提出一種優(yōu)化改進(jìn)的CS-MVDR多目標(biāo)估計方法,通過將空間角度網(wǎng)格劃分形式實現(xiàn)信號在空時域的稀疏性表示,利用稀疏重構(gòu)空間稀疏向量方法實現(xiàn)空間上多信號入射角的方位估計。仿真結(jié)果表明:在低信噪比條件下該方法對多目標(biāo)方位估計精度更高,偏差更小,而且在不損失估計精度的前提下,可以在一定程度上減少雷達(dá)陣面的陣元數(shù),降低復(fù)雜系統(tǒng)的運算量。

壓縮感知理論為工程實踐中信號稀疏化采樣提供了一種新思路,把壓縮感知理論運用到多目標(biāo)方位估計中,可以解決部分工程應(yīng)用中空時欠采樣下的雷達(dá)目標(biāo)回波探測和識別問題,不但降低了數(shù)據(jù)采用率、傳輸量、存儲量,還減少了雷達(dá)陣面的總陣元數(shù)和T/R組件工程設(shè)計及結(jié)構(gòu)減重等問題,同時提高了參數(shù)估計的準(zhǔn)確性,克服了傳統(tǒng)DOA估計算法的不足。