張狀和,韓 東,劉德亮

(陸軍工程大學(xué)石家莊校區(qū) 導(dǎo)彈工程系，石家莊 050003)

共形陣列天線[1]是一種和事先規(guī)定的物體外形共形，而且又不會(huì)帶來(lái)額外負(fù)擔(dān)的天線，這些物體例如飛機(jī)、導(dǎo)彈、高速列車、衛(wèi)星等，共形陣列天線可以滿足上述物體的空氣動(dòng)力學(xué)需求，而且還能有效減少荷載的大小和重量；另一方面，在天線集成時(shí)，諸如城市等通信環(huán)境中亦較少地被干擾。傳統(tǒng)的波達(dá)方向(Dimensional Direction of Arrival,DOA)估計(jì)研究，其仿真實(shí)驗(yàn)大多設(shè)置在線陣、平面陣環(huán)境下，因此研究共形陣DOA估計(jì)方法具有重要意義。目前，在共形陣高精度、高分辨率DOA估計(jì)方面，國(guó)內(nèi)外有很多有價(jià)值的研究成果。

文獻(xiàn)[2]通過(guò)將共形陣分解為較小的子陣列，使用插值技術(shù)，改進(jìn)了多重信號(hào)分類算法(Multiple Signal Classification,MUSIC)，理論上適用于任意形狀的共形陣；文獻(xiàn)[3]使用交叉偶極子建立了四元數(shù)表示的錐面共形陣列模型，通過(guò)秩損原理實(shí)現(xiàn)了高精度2-DDOA估計(jì)；文獻(xiàn)[4]充分利用柱面載體的單曲率特征，應(yīng)用旋轉(zhuǎn)不變參數(shù)估計(jì)技術(shù)(estimation of signal parameters via rotational invariance techniques,ESPRIT)，實(shí)現(xiàn)了盲極化2-DDOA估計(jì)；文獻(xiàn)[5-6]考慮到共形陣載體曲率對(duì)方向圖的影響，利用歐拉旋轉(zhuǎn)變換公式，建立了共形天線陣列導(dǎo)向矢量模型。文獻(xiàn)[7]則提出了一種基于張量技術(shù)的共形陣2-DDOA估計(jì)，該方法利用共形陣接收數(shù)據(jù)的多維信息構(gòu)建三階張量，通過(guò)協(xié)方差張量的高階奇異值分解來(lái)獲得期望信號(hào)子空間，進(jìn)而利用子空間類算法實(shí)現(xiàn)2-D DOA估計(jì)。上述研究均引入了經(jīng)典的子空間類算法來(lái)解決共形陣DOA估計(jì)問(wèn)題，但是這類算法在對(duì)樣本協(xié)方差矩陣進(jìn)行特征分解時(shí)，往往需要大量的快拍數(shù)據(jù)支撐，以保證信號(hào)子空間與噪聲子空間不發(fā)生混疊。但是在實(shí)際的工程應(yīng)用中，信號(hào)經(jīng)常處于不能長(zhǎng)時(shí)間穩(wěn)定或快速時(shí)變的環(huán)境中，采集大量的快拍數(shù)據(jù)會(huì)導(dǎo)致處理時(shí)間過(guò)長(zhǎng)、增大與真實(shí)樣本的誤差等問(wèn)題。相比于子空間類算法，稀疏類算法利用真實(shí)頻率的稀疏性，在實(shí)現(xiàn)高分辨率的同時(shí)，還降低了算法對(duì)快拍數(shù)的需求，2010年P(guān)etre Stoica等[8-10]等關(guān)注了信號(hào)真實(shí)頻率在頻域的稀疏性，提出了基于協(xié)方差稀疏迭代的譜估計(jì)方法(sparse iterative covariance-based estimation method,SPICE)。相比于其他稀疏類方法，SPICE計(jì)算負(fù)擔(dān)較小，且無(wú)需精細(xì)地選擇用戶參數(shù),甚至在只有一個(gè)快拍數(shù)據(jù)的條件下，也可以達(dá)成高精度、高分辨率DOA估計(jì)的目的。

柱面共形陣作為最常見(jiàn)的共形陣形式之一，在通信、雷達(dá)系統(tǒng)的天線集成設(shè)計(jì)等領(lǐng)域有著廣闊的應(yīng)用前景。本文在應(yīng)用歐拉旋轉(zhuǎn)變換建立柱面共形陣的窄帶快拍數(shù)據(jù)模型基礎(chǔ)上，應(yīng)用SPICE方法，將基于統(tǒng)計(jì)可靠的加權(quán)數(shù)據(jù)協(xié)方差擬合代價(jià)函數(shù)轉(zhuǎn)換為凸優(yōu)化問(wèn)題，構(gòu)造循環(huán)最小化器，具有全局收斂特性；最后迭代收斂，繪制出每個(gè)掃描點(diǎn)的三維能量譜圖，譜峰所對(duì)應(yīng)方位角和俯仰角即為信號(hào)源2-D DOA的估計(jì)值。

1 柱面共形陣列數(shù)據(jù)模型

假設(shè)有K個(gè)互不相關(guān)的遠(yuǎn)場(chǎng)窄帶信源入射到M×N元柱面共形陣上，陣元數(shù)為D=N×M，其結(jié)構(gòu)如圖1所示。坐標(biāo)原點(diǎn)的位置即為參考陣元位置，定義n為由下至上的圓環(huán)序號(hào)，m為每一圓環(huán)上逆時(shí)針?lè)较虻年囋蛱?hào)，R表示圓柱半徑，H表示圓柱的高，d為上下兩陣元間的距離；歐拉旋轉(zhuǎn)變換在文獻(xiàn)[5-6]中有詳細(xì)推導(dǎo)，在此共形陣建模的基礎(chǔ)上，以全局坐標(biāo)系為參考，給出如下柱面共形陣列導(dǎo)向矢量：

(1)

u=sin(θ)cos(φ)X+sin(θ)sin(φ)Y+cos(θ)Z

(2)

|gi||pl|cos(θigk)=gi·ql=giθkθ+giφkφ

(3)

其中：Bi,i=1,2,…。D表示陣元位置坐標(biāo)矢量；u表示信源單位矢量坐標(biāo)；(θ，φ)分別表示信源相對(duì)全局直角坐標(biāo)的方位角與俯仰角；ri,i=1,2,…。D表示陣列中第i個(gè)陣元對(duì)入射信號(hào)源的響應(yīng)；gi表示陣元方向圖，kθ,kφ分別表示信源來(lái)波極化狀態(tài)。

圖1 柱面共形陣列結(jié)構(gòu)

接收信號(hào)模型可以表示為

X(t)=A(θ,φ)S(t)+E(t)

(4)

各向量分別定義：

X(t)=[x1(t),x2(t),…,xD(t)]T

(5)

S(t)=[s1(t),s2(t),…,sK(t)]T

(6)

E(n)=[e1(t),e2(t),…,eD(t)]T

(7)

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

(8)

其中:X表示快拍數(shù)據(jù)模型，S表示陣列入射信號(hào)，E表示加性高斯白噪聲，A(θ,φ)表示柱面共形陣列流型矩陣。

2 柱面共形陣2-D DOA算法

通常情況下，實(shí)際感興趣的信源數(shù)是未知的，因此需要把可能存在的目標(biāo)區(qū)域內(nèi)的波達(dá)角度均分為G個(gè)掃描點(diǎn)，每個(gè)點(diǎn)對(duì)應(yīng)潛在的信源，對(duì)每個(gè)信源進(jìn)行能量估計(jì)，得到能量譜，從而得到信源波達(dá)角度的估計(jì)值。實(shí)際上，G的取值往往遠(yuǎn)大于區(qū)域內(nèi)信源的個(gè)數(shù)，導(dǎo)致S中大部分元素的數(shù)值是零，與此同時(shí)，其中的第k個(gè)非零列也對(duì)應(yīng)了位于(θk,φk)處的信源，因此稀疏類算可以適用于DOA估計(jì)。將式(4)整理為稀疏源接收疊加形式：

(9)

假設(shè)：

(10)

(11)

(12)

AH=[a1,…,aKI]=[a1,…,aK+D]

(13)

(14)

(15)

式中||·||表示矩陣的Frobenius范數(shù)，R-1/2為Hermitian矩陣R-1的正定平方根。式(15)可以改寫為如下形式[9]：

tr[(I-xxHR-1)(R-xxH)]=

(16)

其中

(17)

可以看出式(16)、式(17)與下式等價(jià)：

(18)

(19)

式(18)和式(19)是等價(jià)的，式(19)是一個(gè)加權(quán)l(xiāng)1范數(shù)約束問(wèn)題，同樣也是一個(gè)凸優(yōu)化問(wèn)題，構(gòu)造一個(gè)等價(jià)約束最優(yōu)化問(wèn)題來(lái)求解式(19)：

(20)

利用拉格朗日乘子法求解式(20)，構(gòu)造代價(jià)函數(shù)：

J(Q,L)=tr(QHP-1Q)+tr(LH(I-QHA))

(21)

其中，L為拉格朗日乘數(shù)矩陣，求式(21)的梯度可得聯(lián)立方程：

(22)

(23)

聯(lián)立式(22)、式(23)可得：

L=(AHPHA)-1=R-1

(24)

將式(24)代入式(23)可得：

Q=PAR-1

(25)

將式(35)代入式(20)可得：

tr(xHQHP-1Qx)=tr(xHR-1x)

(26)

可以驗(yàn)證式(19)與式(20)等價(jià)，而式(25)即為式(20)的解，上述推導(dǎo)也表明了對(duì)Q，P的循環(huán)最小化可使目標(biāo)函數(shù)式(20)快速收斂。

算法總結(jié)如下：

初始化：

迭代：

3 仿真實(shí)驗(yàn)與分析

為驗(yàn)證本文所提方法的有效性，對(duì)比2D-MUSIC算法展開實(shí)驗(yàn)。為了不失一般性，設(shè)置如圖1所示的接收陣列，M=5、N=10，相鄰圓環(huán)的間距為d=λ/2；噪聲為加性高斯白噪聲；信號(hào)源個(gè)數(shù)K=2，且相互獨(dú)立；來(lái)波信號(hào)角度設(shè)置為θ1=30°，φ1=40°，θ2=35°，φ2=45°；信源極化狀態(tài)kθ=7,kφ=5；假設(shè)在DOA估計(jì)之前已對(duì)感興趣的區(qū)域有一定認(rèn)知，所以提前設(shè)定方位角掃描范圍為「25°，40°?，俯仰角掃描范圍為「35°，50°?，掃描間隔設(shè)置為0.5° ；根據(jù)以往文獻(xiàn)經(jīng)驗(yàn)，迭代次數(shù)預(yù)設(shè)為10。單元方向圖設(shè)置參考文獻(xiàn)[5]：

當(dāng)0≤θ≤π/2時(shí)

(cosφ-jsinφ)

(27)

cosθ(sinφ-jcosφ)

(28)

其他，gθ(θ,φ)=gφ(θ,φ)=0

J0,J2分別是零階和二階第一類貝塞爾函數(shù)。

仿真實(shí)驗(yàn)1：

設(shè)置信噪比(Signal to Noise Ratio,SNR)為15 dB，快拍數(shù)N1=2,N2=200。

仿真結(jié)果表明(圖2)：當(dāng)快拍數(shù)足夠大時(shí)，兩種算法都能成功地把兩個(gè)遠(yuǎn)場(chǎng)信源分辨出來(lái)；而當(dāng)快拍數(shù)減小到2時(shí)，2D-MUSIC算法已經(jīng)失效，真實(shí)信號(hào)源的譜峰混疊在偽峰之中導(dǎo)致無(wú)法分辨出真實(shí)入射角度，但所提算法仍然可以分辨出兩個(gè)銳利的譜峰，且旁瓣水平較低，方位角和俯仰角的估計(jì)偏差也很小。

仿真實(shí)驗(yàn)2：

定義實(shí)驗(yàn)成功須滿足下式：

(29)

定義角度均方根誤差為

(30)

圖3(b)(d)設(shè)置快拍數(shù)N2=200，信噪比從-15 dB變化到15 dB；圖3(a)(c)設(shè)置信噪比為10 dB,快拍數(shù)從10變化到200，其他實(shí)驗(yàn)條件不變。100次蒙特卡洛實(shí)驗(yàn)觀察分辨成功概率和均方根誤差隨快拍數(shù)的變化曲線。

從圖3的實(shí)驗(yàn)結(jié)果可以看出：隨著快拍數(shù)及信噪比的增加，兩種算法的估計(jì)精度逐漸提高，均方根誤差則逐漸減小，說(shuō)明算法的性能在快拍數(shù)據(jù)充足，信噪比較高的情況下變現(xiàn)良好。但是當(dāng)快拍數(shù)與信噪比較低時(shí)，本文提出的算法估計(jì)誤差明顯小于2D-MUSIC算法，成功辨別出兩個(gè)信源的概率也高于2D-MUSIC算法。

圖2 能量譜估計(jì)仿真結(jié)果 (a)(c)2D-MUSIC;(b)(d)所提算法

圖3 均方根誤差及成功概率實(shí)驗(yàn)結(jié)果 (a)(c)快拍數(shù);(b)(d)信噪比

4 結(jié)論

將協(xié)方差稀疏迭代估計(jì)方法應(yīng)用到柱面共形陣2-D DOA估計(jì)上，實(shí)現(xiàn)了短快拍條件下的2-D DOA估計(jì)。所提算法充分利用了協(xié)方差矩陣的稀疏性，利用循環(huán)最小化方法解決全局最優(yōu)化問(wèn)題。相比于傳統(tǒng)子空間類共形陣2-D DOA估計(jì)算法，本文的方法不需要額外選擇用戶參數(shù)，在短快拍、低信噪比的條件下，成功地分辨分布較近的遠(yuǎn)場(chǎng)源，具有較高的估計(jì)精度，為解決短快拍條件下共形陣2-DDOA估計(jì)問(wèn)題提供了重要思路。