李梓正,曹司磊,王 瑤

(海軍航空大學(xué),山東 煙臺(tái) 264000)

0 引言

自適應(yīng)波束形成技術(shù)作為陣列信號(hào)處理的一重要分支,旨在通過(guò)一定的優(yōu)化準(zhǔn)則,計(jì)算陣列最優(yōu)加權(quán)矢量,達(dá)到有效抑制干擾、噪聲來(lái)獲得高處理增益并提高輸出信干噪比(signal to interference and noise ratio,SINR)的效果,在通信系統(tǒng)、航天、醫(yī)療、雷達(dá)等領(lǐng)域中均有廣泛應(yīng)用[1-2]。典型的波束形成器,如標(biāo)準(zhǔn)Capon波束形成器,在理想情況下可達(dá)到最優(yōu)性能[3],然而,實(shí)際工作場(chǎng)景中由于存在目標(biāo)來(lái)波方向誤差、陣列位置誤差等,導(dǎo)致波束形成器輸出性能下降嚴(yán)重。

目前,諸多學(xué)者在穩(wěn)健自適應(yīng)波束形成領(lǐng)域進(jìn)行了深入研究,提出的方法主要涵蓋:加載類算法[4-9]、不確定集算法[10-14]和協(xié)方差矩陣重構(gòu)類算法[15-20]等。

文獻(xiàn)[4]提出了對(duì)樣本協(xié)方差矩陣的對(duì)角元素加載一個(gè)常數(shù)以提高魯棒性的對(duì)角加載(loading sample matrix inversion,LSMI)算法,該算法普適性較強(qiáng),但對(duì)最佳加載值的選取并未有確定的方法,加載值過(guò)大易降低干擾抑制能力,加載值過(guò)小對(duì)算法的改善影響較為有限。文獻(xiàn)[8]提出了一種自適應(yīng)可變對(duì)角加載技術(shù)(automatic variable loading technology,AVL-RAB)當(dāng)輸入信噪比(signal to noise ratio,SNR)較低時(shí)以非均勻加載矩陣為主導(dǎo),當(dāng)輸入SNR較高時(shí)以可變對(duì)角加載矩陣為主導(dǎo),提升了高SNR下輸出性能,但算法抗干擾能力有所減弱。

文獻(xiàn)[12]提出了目標(biāo)導(dǎo)向矢量雙層估計(jì)的穩(wěn)健波束形成算法,首先將目標(biāo)信號(hào)大致方位區(qū)間上的導(dǎo)向矢量構(gòu)建一個(gè)線性空間,并將其用線性空間基向量組合表示,然后構(gòu)建不確定集約束校正目標(biāo)信號(hào)導(dǎo)向矢量,提高了波束形成器導(dǎo)向矢量失配誤差下的穩(wěn)健性,但在低SNR環(huán)境下,輸出性能下降嚴(yán)重。文獻(xiàn)[16]研究了穩(wěn)健Capon波束形成(robust capon beamforming,RCB)算法,建立導(dǎo)向矢量不確定集約束,求解不確定集系數(shù),當(dāng)導(dǎo)向矢量誤差較小時(shí),算法魯棒性較好,但當(dāng)導(dǎo)向矢量失配時(shí)性能下降嚴(yán)重。文獻(xiàn)[21]通過(guò)對(duì)目標(biāo)信號(hào)導(dǎo)向矢量進(jìn)行線性約束,提升了波束形成器的穩(wěn)健性,并添加了阻塞矩陣預(yù)選環(huán)節(jié)解決了線性約束帶來(lái)的自由度損失問(wèn)題。

樣本協(xié)方差矩陣中含有期望信號(hào)分量,會(huì)極大影響算法的性能,剔除期望信號(hào)分量時(shí)易引入其他誤差,影響后續(xù)的計(jì)算過(guò)程;不確定集算法約束的大小難以確定,在各種誤差存在的情況下性能波動(dòng)較大。針對(duì)上述問(wèn)題,本文中提出基于干擾加噪聲協(xié)方差矩陣優(yōu)化重構(gòu)的穩(wěn)健波束形成算法。首先,通過(guò)Capon功率譜得到目標(biāo)導(dǎo)向矢量的預(yù)估值,重構(gòu)干擾加噪聲協(xié)方差矩陣;然后,建立干擾導(dǎo)向矢量不確定集約束,干擾加噪聲矩陣的優(yōu)化估計(jì)值;繼而,根據(jù)子空間擴(kuò)展知識(shí),構(gòu)造導(dǎo)向矢量誤差不確定集約束模型,循環(huán)迭代更新目標(biāo)導(dǎo)向矢量和干擾加噪聲協(xié)方差矩陣,最終得到二者穩(wěn)態(tài)估計(jì)值并計(jì)算加權(quán)矢量。

1 陣列信號(hào)模型

考慮陣元數(shù)為M的理想一維均勻線陣,假設(shè)空間中有K個(gè)入射信號(hào),目標(biāo)信號(hào)從θ0方向入射,K-1個(gè)干擾信號(hào)從θi(i=1,2,…,K-1)方向入射,假設(shè)入射信號(hào)與噪聲之間相互統(tǒng)計(jì)獨(dú)立,則陣列在t時(shí)刻的接收信號(hào)可表示為

(1)

a(θi)=[1,ej(2πd/λ)sinθi,…,ej(M-1)(2πd/λ)sinθi]H

(2)

其中:a(θi)和si(t)分別為第i個(gè)射信號(hào)的導(dǎo)向矢量和復(fù)包絡(luò),n(t)為噪聲信號(hào),λ為信號(hào)波長(zhǎng),d為陣元間距。

陣列接收信號(hào)的協(xié)方差矩陣可表示為

R=E[X(t)XH(t)]=Rs+Ri,n

(3)

(4)

(5)

實(shí)際工作中,矩陣R一般是未知的,故采用樣本協(xié)方差矩陣代替如式(6)所示:

(6)

式中:P為快拍數(shù)。

2 干擾加噪聲協(xié)方差矩陣優(yōu)化重構(gòu)

2.1 干擾加噪聲協(xié)方差矩陣預(yù)估計(jì)

假設(shè)目標(biāo)大致出現(xiàn)的區(qū)域?yàn)棣?在該區(qū)域內(nèi)對(duì)Capon功率譜求積分,則有

(7)

將矩陣Q特征值分解,以其最顯著特征值對(duì)應(yīng)的特征向量作為目標(biāo)導(dǎo)向矢量的估計(jì)值。

(8)

式中:Umax是Q特征值分解后最顯著特征值對(duì)應(yīng)特征向量。

從而可將期望信號(hào)的自相關(guān)矩陣估計(jì)為

(9)

(10)

式中:σd表示對(duì)角加載量。

2.2 干擾加噪聲協(xié)方差矩陣優(yōu)化重構(gòu)

(11)

(12)

通過(guò)拉格朗日乘數(shù)法對(duì)式(12)求解,則干擾導(dǎo)向矢量可表示為

(13)

式中:參數(shù)μ的值由如下計(jì)算式給出

(14)

得出干擾導(dǎo)向矢量的優(yōu)化估計(jì)值后,繼而求解干擾信號(hào)功率。

(15)

設(shè)a(θl)與a(θk)為入射方向不同的信號(hào)導(dǎo)向矢量,在信號(hào)入射方位稀疏的條件下,不同角度的信號(hào)導(dǎo)向矢量相互正交[25],則當(dāng)l≠k時(shí),有aH(θl)a(θk)=0。在式(15)的左右兩端分別乘aH(θl)和a(θl)得

(16)

進(jìn)而可得出干擾信號(hào)功率表達(dá)式:

(17)

假設(shè)噪聲扇區(qū)為Θn,則可由下式計(jì)算剩余噪聲[26]

(18)

式中:N為噪聲扇區(qū)中的采樣點(diǎn)數(shù)量,a(θn)為角度θn對(duì)應(yīng)的標(biāo)稱導(dǎo)向矢量。則噪聲功率可估計(jì)為

(19)

(20)

(21)

3 目標(biāo)導(dǎo)向矢量估計(jì)

(22)

基于此,為最大化陣列輸出功率,并減小干擾及噪聲對(duì)目標(biāo)導(dǎo)向矢量的影響,建立約束關(guān)系如下:

(23)

(24)

4 算法求解及步驟

設(shè)參數(shù)β為收斂系數(shù),參數(shù)itmax為最大迭代次數(shù),設(shè)

(25)

則收斂條件可表示為

(26)

式中:t為已完成的迭代次數(shù),當(dāng)t≤itmax時(shí),以式(26)第1、2項(xiàng)不等式約束為終止判決條件;而當(dāng)t>itmax時(shí),迭代終止。

(27)

式中:ω為加權(quán)矢量。

分析式(26)易知,β和itmax的大小會(huì)影響迭代結(jié)果的精度和迭代的次數(shù);β過(guò)小或itmax過(guò)大會(huì)導(dǎo)致迭代次數(shù)增多,算法迭代時(shí)間過(guò)長(zhǎng);β過(guò)大或itmax過(guò)小會(huì)導(dǎo)致迭代精度下降。本文中算法流程如圖1所示。

圖1 算法流程圖

步驟4由式(26)判斷是否滿足迭代終止條件,不滿足則重復(fù)上述步驟。

5 仿真實(shí)驗(yàn)

考慮陣元數(shù)為8的理想一維均勻線陣,陣元間距設(shè)置為0.5λ。本節(jié)中,假設(shè)目標(biāo)信號(hào)實(shí)際來(lái)波方向?yàn)?,干擾信號(hào)來(lái)波方向?yàn)?10°和20°,干噪比均設(shè)為30 dB,噪聲均為高斯白噪聲,目標(biāo)信號(hào)與干擾信號(hào)間統(tǒng)計(jì)無(wú)關(guān)。將本文中算法與LSMI算法、RCB算法、AVL-RAB算法分析對(duì)比,每次仿真均進(jìn)行100次獨(dú)立的蒙特卡洛實(shí)驗(yàn)。

5.1 抗目標(biāo)來(lái)波方向誤差算法仿真及分析

本小節(jié)中,設(shè)期望信號(hào)來(lái)波方向誤差服從[-3,3]區(qū)間內(nèi)的隨機(jī)分布,則實(shí)際信號(hào)來(lái)波方向服從[2,8]區(qū)間內(nèi)的隨機(jī)分布,快拍數(shù)設(shè)為200,輸入SNR取[-30 dB,40 dB];分析各算法輸出性能隨快拍數(shù)變化趨勢(shì)時(shí),輸入SNR設(shè)為1 dB。將各算法100次蒙特卡洛實(shí)驗(yàn)的輸出SINR做統(tǒng)計(jì)平均,仿真結(jié)果如圖2所示。

圖2 目標(biāo)來(lái)波方向誤差下仿真對(duì)比

分析圖2(a)可知,當(dāng)目標(biāo)來(lái)波方向誤差存在時(shí),LSMI算法在10 dB輸出性能達(dá)到最優(yōu),隨著SNR提高,算法性能下降嚴(yán)重;RCB算法用于求取權(quán)值的樣本協(xié)方差矩陣仍包含期望信號(hào)在高SNR輸入下,性能下降明顯,曲線斜率變小;AVL-RAB算法采用了可變加載值加載矩陣表現(xiàn)整體優(yōu)于RCB算法,但高SNR輸入下性能仍有不足;本文中算法輸出性能總體較好,尤其在高SNR輸入下輸出性能仍穩(wěn)定提升。

分析圖2(b)可知,本文中算法與AVL-RAB算法總體表現(xiàn)較好,尤其是在小塊拍數(shù)下,快拍數(shù)大于50次后,基本達(dá)到收斂,隨快拍數(shù)增大二者輸出SINR提高速度變緩,且逐步接近理論最優(yōu)值。

5.2 抗陣元位置誤差算法仿真及分析

設(shè)每個(gè)陣元實(shí)際位置相對(duì)理想位置偏移誤差在[-0.05λ,0.05λ]區(qū)間內(nèi)隨機(jī)分布,期望信號(hào)來(lái)波方向誤差設(shè)置為3,輸入信噪比取[-30 dB,40 dB],快拍數(shù)為200;仿真各算法輸出性能隨快拍數(shù)變化趨勢(shì)時(shí),輸入SNR設(shè)為0 dB。將各算法100次蒙特卡洛實(shí)驗(yàn)的輸出SINR做統(tǒng)計(jì)平均,仿真結(jié)果如圖3所示。

圖3 陣列位置誤差下仿真對(duì)比

分析圖3(a)可知,當(dāng)陣列位置誤差存在時(shí),LSMI算法性能隨輸入SNR提高下降嚴(yán)重,說(shuō)明其抗陣列位置誤差干擾性能較差,AVL-RAB算法輸出性能總體優(yōu)于RCB算法,但小塊拍數(shù)下性能有所下降,是因正交矩陣特征向量計(jì)算不準(zhǔn)確導(dǎo)致,本文中算法對(duì)干擾加噪聲協(xié)方差矩陣進(jìn)行了優(yōu)化重構(gòu),表現(xiàn)優(yōu)于其他3種算法,具備對(duì)陣列位置誤差的穩(wěn)健性。

分析圖3(b)可知,本文中算法在小快拍數(shù)下輸出性能較其他算法有明顯優(yōu)勢(shì),且收斂速度較快,當(dāng)快拍數(shù)大于70次后接近收斂,并逐漸接近理論最優(yōu)值。

5.3 抗陣列互耦誤差算法仿真及分析

設(shè)互耦矩陣為C,對(duì)理想一維均勻線陣,矩陣C為Toeplitz形式,本小節(jié)中,僅考慮相鄰陣元間的互耦效應(yīng),互耦系數(shù)設(shè)為0.5exp-jπ/2,其余實(shí)驗(yàn)條件設(shè)置如上,無(wú)特殊說(shuō)明,將100次蒙特卡洛實(shí)驗(yàn)的輸出做統(tǒng)計(jì)平均,仿真結(jié)果如圖4所示。

圖4 陣列互耦誤差下仿真對(duì)比

分析圖4(b)可知,本文中算法輸出SINR隨快拍數(shù)穩(wěn)健提升,且在較小快拍數(shù)條件下較其他算法仍有優(yōu)勢(shì),在快拍數(shù)大于70次后基本收斂。

6 結(jié)論

本文中提出了基于干擾加噪聲協(xié)方差矩陣優(yōu)化重構(gòu)的穩(wěn)健波束形成算法。算法的基本思想是基于目標(biāo)導(dǎo)向矢量和干擾加噪聲協(xié)方差矩陣的聯(lián)系,稀疏重構(gòu)協(xié)方差矩陣;采用循環(huán)迭代的方法求解建立的目標(biāo)導(dǎo)向矢量誤差優(yōu)化模型,得出最優(yōu)加權(quán)矢量。仿真實(shí)驗(yàn)表明:該算法抗文中所述的3種失配誤差能力較強(qiáng),穩(wěn)健性較好,尤其在快拍數(shù)較小或SNR較高的條件下輸出性能較之其他算法明顯提升。但本文中算法綜合復(fù)雜度較高,降低算法的復(fù)雜度是后續(xù)研究中的重點(diǎn)內(nèi)容。