一種基于對角加載的自適應(yīng)波束形成算法
【專利摘要】本發(fā)明公開了一種基于對角加載的自適應(yīng)波束形成算法，涉及智能天線技術(shù)領(lǐng)域，首先對線陣接收陣元所采集到的采樣信號求其采樣協(xié)方差矩陣，做為樣本協(xié)方差矩陣的估計(jì)。然后利用對角加載技術(shù)對采樣協(xié)方差矩陣進(jìn)行重構(gòu)，使其滿足矩陣求逆引理公式，避免進(jìn)行矩陣求逆運(yùn)算。最后再結(jié)合最小均方誤差(MSE)準(zhǔn)則，得到方向權(quán)向量的最優(yōu)化解，運(yùn)用重構(gòu)的采樣協(xié)方差矩陣取代了迭代運(yùn)算，大大縮減了算法收斂時(shí)間。該算法不僅有效解決和優(yōu)化了自適應(yīng)數(shù)字波束形成算法的收斂時(shí)間問題，并且通過仿真實(shí)驗(yàn)驗(yàn)證了該算法在高低信噪比環(huán)境下性能都比較穩(wěn)定，同時(shí)也可以在一定程度上消除對模型誤差的敏感問題。
【專利說明】
一種基于對角加載的自適應(yīng)波束形成算法
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及智能天線技術(shù)領(lǐng)域，尤其涉及一種基于對角加載的自適應(yīng)波束形成算 法。
【背景技術(shù)】
[0002] 水下聲成像技術(shù)主要用于水下目標(biāo)探測與搜索，水底地貌繪制，海底沉船、救援、 黑盒子打撈等眾多軍事和民用領(lǐng)域。已經(jīng)在國防和民用領(lǐng)域取得廣泛的應(yīng)用。為了實(shí)現(xiàn)一 定距離下水下目標(biāo)的高清晰聲成像目的，就必須研究復(fù)雜水聲環(huán)境下的穩(wěn)定成像技術(shù)，并 且同時(shí)要滿足一定的成像幀率和較遠(yuǎn)的作用距離。這就要求在研究水聲成像的自適應(yīng)波束 形成算法時(shí)，不僅也需要考慮魯棒性問題，也要盡量提高系統(tǒng)的輸出信噪比，同時(shí)減小自適 應(yīng)波束形成算法的運(yùn)算量。
[0003] 自適應(yīng)波束形成技術(shù)通過調(diào)整權(quán)重向量來改變陣列的方向圖，使波束主瓣對準(zhǔn)期 望信號，旁瓣和零陷對準(zhǔn)干擾信號，從而提高輸出的信干噪比，以實(shí)現(xiàn)某準(zhǔn)則下的最佳接 收。
[0004] LMS自適應(yīng)波束形成算法是一種結(jié)構(gòu)簡單、算法復(fù)雜度低、易于實(shí)現(xiàn)和穩(wěn)定性高的 波束形成方法。但一直因?yàn)槠涫諗克俣容^慢，在工程應(yīng)用上受到一定程度的限制。為此，各 學(xué)者相繼以不同的調(diào)整策略:瞬時(shí)誤差、權(quán)矢量的前向預(yù)測以及平滑梯度矢量等，提出變步 長的LMS算法來平衡收斂速度和算法失調(diào)。雖然在平衡收斂速度和失調(diào)方面優(yōu)于經(jīng)典LMS算 法，但應(yīng)對突變能力較差。
[0005] 作為判斷最佳接收的準(zhǔn)則之一，均方誤差(MSE)性能量度由威德魯?shù)热颂岢觥２⒂?Wiener、H〇pf推導(dǎo)出最優(yōu)的維納解。經(jīng)典LMS算法正是在MSE準(zhǔn)則的基礎(chǔ)上，運(yùn)用最優(yōu)化方法 如:最速下降法、加速梯度算法等迭代運(yùn)算出最優(yōu)權(quán)向量的。
[0006] 基于以上所述，本文提出一種基于MSE準(zhǔn)則和對角加載技術(shù)的魯棒自適應(yīng)波束形 成算法。在該算法中，通過在采樣協(xié)方差矩陣對角線上人工注入白噪聲，即對角加載，重構(gòu) 采樣協(xié)方差矩陣。然后運(yùn)用矩陣求逆引理，避免矩陣求逆運(yùn)算和迭代運(yùn)算，而且將對角加載 系數(shù)轉(zhuǎn)化為LMS算法步長因子的函數(shù)。仿真結(jié)果表明此算法不僅能有效降低收斂時(shí)間，并在 高、低信噪比環(huán)境下能都表現(xiàn)出較好的性能，有較好的魯棒性。

【發(fā)明內(nèi)容】

[0007] 為了改進(jìn)現(xiàn)有自適應(yīng)波束形成方法在收斂速度上偏慢的缺點(diǎn)，突破采樣頻率限 制，獲得更精準(zhǔn)的權(quán)重方向向量，本發(fā)明的第一目的是在LMS算法和MVDR的基礎(chǔ)上避免循環(huán) 迭代和矩陣求逆運(yùn)算，縮短算法收斂時(shí)間，使其能很好地在工程中得到應(yīng)用。該方法可以在 形成自適應(yīng)波束中通過重構(gòu)采樣協(xié)方差矩陣，并在此基礎(chǔ)上應(yīng)用矩陣求逆引理，從而避免 了求逆運(yùn)算和迭代運(yùn)算，具有較好的指向性能和干擾抑制能力。
[0008] 本發(fā)明的第二目的在于為減小各種誤差導(dǎo)致的副瓣電平升高、主瓣便宜、波束畸 變、SINR下降等問題，將對角加載技術(shù)引入自適應(yīng)波束形成中，并給出了確定加載系數(shù)的公 式。該方法實(shí)現(xiàn)簡單，有利于減少波束形成過程中的偏差，提高波束形成的準(zhǔn)確性和穩(wěn)健 性。
[0009]為實(shí)現(xiàn)上述目的，本發(fā)明提供一種基于對角加載的自適應(yīng)波束形成算法，所述方 法包含如下步驟：
[0010]步驟1:考慮平面空間的等距均勻線陣，設(shè)陣元數(shù)為M，陣元間距為d，其中(1 = λ/2(λ 為陣列接收單元接收信號的波長），假設(shè)有L個(gè)信源回波(M>L)，設(shè)波達(dá)方向?yàn)?^02,...，0L， 以陣列的第一個(gè)陣元作為基準(zhǔn)點(diǎn)，則在第k次快拍的采樣點(diǎn)m處的采樣值為：
[0012] 式中nm(k)表示第m個(gè)陣元上的噪聲，Sl(k)表示各信源回波在基準(zhǔn)點(diǎn)的基帶信號。 步驟2:各陣元在快拍k時(shí)刻接收到的信號分別為XKk)，X 2(k)，…，XM(k)，即:X(k) = [XKk)， X 2 ( k )，. . .，X μ ( k ) ] τ，此為陣列輸入矢量。得到協(xié)方差矩陣估計(jì)值為
。式中Κ表示陣列天線的快拍數(shù)，X(k)表示陣列天線上第k次快拍接 收到的信號(k=l，2,. . .，K)，上標(biāo)Η表示矩陣共輒轉(zhuǎn)置。
[0013] 步驟3:在時(shí)域中，陣列輸出為
[0014] y(t)=coTX(t) (2)
[0015] 參考信號d(t)與實(shí)際輸出信號的誤差為
[0016] ε ⑴=d(t)_y(t) =d(t)_ ω TX(t) (3)
[0017] 對(3)式求平方可得
[0018] e2(t) = d2(t)-2d(t)〇TX(t)+〇TX(t)XT(t)〇 (4)
[0019] 對上式兩邊取數(shù)學(xué)期望可得
[0020] (5；
[0021] 式中^表示對d(t)取數(shù)學(xué)期望，互相關(guān)矩陣Rxd為Rxd = E{d(k)XT(k)}，令 i/:(〇=S，:貝丨J
[0022] E{e2(t)} = S-2〇TRxd+oTRxxo (6)
[0023] 適當(dāng)選擇權(quán)重向量ω可使E{>2(t)}達(dá)到最小?？芍?6)是ω的二次函數(shù)，該函數(shù) 的極值是一個(gè)最小值，由式(6)對權(quán)重向量求梯度并令其為零，求出使E{> 2(t)}最小的ω 值，得到權(quán)重向量的最優(yōu)值滿足下式：
[0024] if)
[0025] 步驟4:在波束形成算法方面，LMS算法作為常步長LMS算法，其權(quán)向量的迭代公式 可表述為：
[0026] (〇{k +1) = o.){k) - //V(/;) C 8 )
[0027] 為了克服矩陣求逆等運(yùn)算，LMS算法采用最陡下降法求解式(8)，得到LMS算法的迭 代公式
[0028] ω (k+Ι) = ω (k)+yX(k)e*(k) (9)
[0029] 式中，μ為步長因子，可以控制自適應(yīng)的速率。通過分析可知，μ步長因子的取值范 圍滿足關(guān)系：
，可以證明當(dāng)?shù)螖?shù)無限增加時(shí)，權(quán)重向量的期望值可以收斂至 維納解。
[0030]步驟5:在增強(qiáng)自適應(yīng)波束形成器的魯棒性方面，對角加載技術(shù)被用于抑制方向圖 畸變。在本文所依據(jù)的信號模型基礎(chǔ)上，實(shí)際計(jì)算采樣協(xié)方差矩陣Rxx是由K次采樣信號得到 的估計(jì)值
[0034] 將對角加載技術(shù)運(yùn)用到LMS算法的權(quán)向量計(jì)算中，得到
[0035] 反,.=(?/ + 疋、） （12)
[0036]引理:令矩陣Ae CnXn的逆矩陣存在，并且X，y是兩個(gè)η X 1維向量，使得(A+xyH)可 逆，則
[0038] 將其推廣為矩陣之和求逆公式，即為：
[0039] (A+UBV) -1=A-〔A-bB (B+BVA-bB)-^VA-1
[0040] =A-〔A-bd+BVA-%)-hVA-1 (14)
[0041 ]因?yàn)椴蓸訁f(xié)方差見^Hermit ian矩陣，則由式（11)可以推導(dǎo)出
[0043] 式中，U為特征向量矩陣，Λ =diag( γι, γ2,...，Ym)，Yi為及#的特征值。
[0044] 根據(jù)上述矩陣求逆公式可以推導(dǎo)出
[0046]由式（10)可知及"是11個(gè)陣元的K次采樣數(shù)據(jù)相關(guān)矩陣的均值。將我"轉(zhuǎn)化為譜分解 形式后(式（15))應(yīng)用到上式的推導(dǎo)過程中，替代第k次采樣數(shù)據(jù)相關(guān)矩陣X(k)PXH(k)(P取 單位矩陣)得到
[0048] 然后將忒Π 7)代入結(jié)里Π
2)中，_，
[0050] 式中α表示加載系數(shù)，定義：
，其中0〈λ〈1。因此，對角加載系數(shù)的確定可以 由LMS算法中的步長因子μ和λ確定。
[0051] 最后，輸出自適應(yīng)波束為
[0052] y(k) = ω 〇ptTX(k) (19)
[0053] 本發(fā)明的有益效果是:通過在采樣協(xié)方差矩陣對角線上人工注入白噪聲，即對角 加載，重構(gòu)采樣協(xié)方差矩陣。然后運(yùn)用矩陣求逆引理，推導(dǎo)出一種魯棒的自適應(yīng)波束形成 算法。重構(gòu)后的采樣協(xié)方差矩陣滿足矩陣求逆引理的條件，推導(dǎo)出的權(quán)向量公式避免了矩 陣求逆和循環(huán)迭代，實(shí)現(xiàn)了快速收斂的目的。同時(shí)，重構(gòu)的采樣協(xié)方差矩陣引入了對角加載 因子，使得該算法的具有一定的魯棒性。通過實(shí)驗(yàn)的驗(yàn)證，該算法可以同時(shí)應(yīng)用在低信噪比 和高信噪比的各種復(fù)雜環(huán)境中，并且在收斂速度和輸出信噪比上相對傳統(tǒng)的MVDR和LMS算 法有很大的提升，使其可以應(yīng)用在水聲成像的復(fù)雜環(huán)境中，從而保證了成像的穩(wěn)定性和成 像的幀率。
【附圖說明】
[0054]圖1流程圖。
[0055] 圖2低信噪比（_3dB)環(huán)境下各自適應(yīng)波束指向圖。
[0056] 圖3高信噪比（30dB)環(huán)境下各自適應(yīng)波束指向圖 [0057]圖4高信噪比（40dB)環(huán)境下各自適應(yīng)波束指向圖
[0058]圖5低信噪比（_3dB)對角加載因子對波束形成的影響圖 [0059 ]圖6高信噪比(40dB)對角加載因子對波束形成的影響圖
【具體實(shí)施方式】
[0060] 以下將結(jié)合附圖，對本發(fā)明的優(yōu)選實(shí)施例進(jìn)行詳細(xì)的描述;應(yīng)當(dāng)理解，優(yōu)選實(shí)施例 僅為了說明本發(fā)明，而不是為了限制本發(fā)明的保護(hù)范圍。
[0061] 圖1是本發(fā)明算法的流程圖，如圖所示:本發(fā)明提供的一種基于對角加載的自適應(yīng) 波束形成方法，包括以下步驟：
[0062] S1:均勻線陣各陣元對信號進(jìn)行采樣X(k);
[0063] S2:對采樣信號求其采樣協(xié)方差矩陣
，做為樣本協(xié)方差矩陣 的估計(jì)；
[0064] S3:利用對角加載技術(shù)對采樣協(xié)方差矩陣進(jìn)行重構(gòu)&
[0065] S4:結(jié)合最小均方誤差(MSE)準(zhǔn)則，計(jì)算出方向權(quán)向量的最優(yōu)解= , 將S3中的及替代，并應(yīng)用矩陣求逆公式進(jìn)行推導(dǎo)，得到
[0066] S5:將得到的方向權(quán)值對采樣信號數(shù)據(jù)進(jìn)行加權(quán)求和，得到自適應(yīng)波束信號y(k) -〇 opt X (k) 〇
[0067] 具體實(shí)施步驟：
[0068] 步驟1:根據(jù)本算法依據(jù)的信號模型，第k次快拍的采樣點(diǎn)m的采樣值為：
[0070]各陣元在快拍k時(shí)刻接收到的信號分別為XKk)，X2(k)，…，XM(k)，即：X(k) = [Xi (k)，X2(k)，...，XM(k)]T;
[0071 ]步驟2:得到協(xié)方差矩陣估計(jì)值為
[0073]做為樣本協(xié)方差矩陣的估計(jì)；
[0074] 步驟3:將對角加載技術(shù)運(yùn)用到算法的權(quán)向量計(jì)算中，得到
[0075] 4, ={aI + RJ
[0076] ^是Μ個(gè)陣元的K次采樣數(shù)據(jù)相關(guān)矩陣的均值，α為對角加載系數(shù)。
[0077] 步驟4:將轉(zhuǎn)化為譜分解形式
，替代第k次采樣數(shù)據(jù)相關(guān) 矩陣X(k)PXH(k)(P取單位矩陣)并應(yīng)用到矩陣求逆公式的推導(dǎo)過程中，得到
[0079]將式（17)代入結(jié)果（12)中，則：
[0082] 步驟5:將得到的方向權(quán)值對采樣信號數(shù)據(jù)進(jìn)行加權(quán)求和，得到自適應(yīng)波束信號y (k) = ω 〇ptTX(k) 〇
[0083] 為了驗(yàn)證該算法的有效性，利用MATLAB仿真工具進(jìn)行算法仿真。仿真實(shí)驗(yàn)采用由 16個(gè)陣元組成的均勻線陣，陣元間隔為半個(gè)波長。假設(shè)期望信號和干擾的波達(dá)方向分別為 0°和40°，并且期望信號和干擾互不相干。噪聲均值為0,方差為1的加性高斯白噪聲。在仿真 實(shí)驗(yàn)中，將本文提出的算法DL-MSE與經(jīng)典的LMS算法、MVDR算法進(jìn)行對比分析。采樣數(shù)均為 500，LMS算法的迭代次數(shù)也為500 4 = 0 · 0005，λ = 〇 · 5。
[0084] 實(shí)驗(yàn)1:在該實(shí)驗(yàn)中，驗(yàn)證各種算法在低信噪比環(huán)境下的指向性能，取SNR = _3dB。 結(jié)果如圖2所示，在低信噪比的情況下，LMS算法性能嚴(yán)重失調(diào)，跟蹤效果變得很差。而本文 提出的算法DL-MSE和MVDR算法的效果比較接近，性能良好。在干擾方向上，DL-MSE算法稍差 于MVDR算法。
[0085]實(shí)驗(yàn)2:在該實(shí)驗(yàn)中，驗(yàn)證各種算法在高信噪比環(huán)境下的指向性能，取SNR分別為 30dB、40dB。結(jié)果如圖3、圖4所示，在高信噪比的情況下，MVDR算法性能嚴(yán)重失調(diào)，而本文提 出的算法DL-MSE算法和LMS算法的性能比較穩(wěn)定。在干擾抑制上，DL-MSE算法稍優(yōu)于LMS算 法，都能在干擾方向上形成零陷。隨著SNR的升高，對比圖3和圖4可以發(fā)現(xiàn)DL-MSE算法和LMS 算法都表現(xiàn)出較好的性能，能保持穩(wěn)定的跟蹤指向性能。對比Cox等提出的對角加載方法， 本文提出的方法解決了在較高信噪比（SNR)條件下，采用對角加載方法的自適應(yīng)波束形成 器有較嚴(yán)重的性能衰落。
[0086]實(shí)驗(yàn)3:在本文提出設(shè)計(jì)的自適應(yīng)波束形成算法中，作為對角加載的系數(shù)，α和LMS 算法步長因子之間滿足線性關(guān)系。當(dāng)系數(shù)λ取不同值時(shí)，加載系數(shù)隨之發(fā)生變化。圖5、圖6說 明了在低信噪比（_3dB)和高信噪比（40dB)兩種情況下λ取不同的值時(shí)，DL-MSE算法的性能 變化。在低信噪比時(shí)，加載系數(shù)
越大，算法對噪聲的抑制效果越好，并在干擾方向上 的形成較深的零陷。而在高信噪比的環(huán)境中，算法性能對加載系數(shù)的變化不敏感，但也同樣 能在干擾方向上形成較好的零陷。
[0087]在算法的收斂速度上，本文提出的算法由于避免了MVDR算法的矩陣求逆運(yùn)算和 LMS算法循環(huán)迭代更新權(quán)向量，所以在收斂速度上相對MVDR算法和LMS算法都有一定的優(yōu) 勢，表1列出了在采樣500次或者LMS迭代500次的情況下，自適應(yīng)波束收斂情況。從表1對比 可知，通過采樣后進(jìn)行協(xié)方差矩陣求逆運(yùn)算得到權(quán)向量的MVDR和通過循環(huán)迭代的LMS算法 在收斂速度上差距不大。而本文提出的DL-MSE算法在收斂速度上有很大的優(yōu)勢，說明本發(fā) 明的算法可以被應(yīng)用在實(shí)時(shí)性要求較高的場合中。
[0089] 表1各自適應(yīng)波束權(quán)向量形成時(shí)間對比
[0090] 最后說明的是，以上所述僅為本發(fā)明的較佳實(shí)例而已，并不用以限制本發(fā)明，凡在 本發(fā)明的精神和原則之內(nèi)所作的任何修改，等同替換和改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù) 范圍之內(nèi)。
【主權(quán)項(xiàng)】
1. 一種基于對角加載的自適應(yīng)波束形成算法，其特征在于:包括W下步驟： S1:均勻線陣各陣元對信號進(jìn)行采樣； S2:求其采樣協(xié)方差矩陣，作為樣本協(xié)方差矩陣的估計(jì)； S3:利用對角加載技術(shù)對采樣協(xié)方差矩陣進(jìn)行重構(gòu)； S4:結(jié)合最小均方誤差(MSE)準(zhǔn)則，計(jì)算出方向權(quán)向量的最優(yōu)解； S5:將得到的方向權(quán)值對采樣信號數(shù)據(jù)進(jìn)行加權(quán)求和，得到自適應(yīng)波束信號。2. 根據(jù)權(quán)利要求1所述的一種基于對角加載的自適應(yīng)波束形成算法，其 特征在于：所述步驟1中：各陣元對信號進(jìn)行采樣時(shí)，依據(jù)信號模型其中nm(k)表示第m個(gè)陣元上的噪聲，si(k) 表示各信源回波在基準(zhǔn)點(diǎn)的基帶信號，L表示信源個(gè)數(shù);得到第k次快拍的采樣值為Xi化），拉 化），···，Χμ化），即:X化）= [Xi化），X2化），...，Xm化）]τ，其中Μ表示陣元個(gè)數(shù)。3. 根據(jù)權(quán)利要求1所述的一種基于對角加載的自適應(yīng)波束形成算法，其特征在于:所述 步驟2中：對接收向量用二階統(tǒng)計(jì)量的平均值進(jìn)行表示，即信號自相關(guān)矩陣，作為協(xié)方差矩 陣估計(jì)值廷中Κ表示陣列天線的快拍數(shù)，X化)表示陣列天線上第k次 快拍接收到的信號化=1，2,...，Κ)，上標(biāo)Η表示矩陣共輛轉(zhuǎn)置。4. 根據(jù)權(quán)利要求1所述的一種基于對角加載的自適應(yīng)波束形成算法，其特征在于:所述 步驟3中：將對角加載技術(shù)運(yùn)用到算法的權(quán)向量計(jì)算中，得到式中4是Μ個(gè)陣元的Κ次采樣數(shù)據(jù)相關(guān)矩陣的均值，α為對角加載系數(shù)。5. 根據(jù)權(quán)利要求1所述的一種基于對角加載的自適應(yīng)波束形成算法，其特征在于:所述 步驟4中：得到方向權(quán)向量的最優(yōu)解。具體包括W下步驟： S4 1 :由于采樣協(xié)方差.《"為Hermi t ian矩陣，則可表示譜分解形式：妻中U為特征向量矩陣，Λ =diag(丫 1，丫 2, . . .，丫M)，丫 i為矣。的特 征值； S42:由矩陣求逆引理將[α?巧化)ΡχΗ化)]-1轉(zhuǎn)換關(guān)菱示. 543:將541步驟中的4替換542步驟中的乂（1〇口乂"（1〇"取單位矩陣）得到其中α表示加載系數(shù)； S44:定義對角加載系i廷中μ表示LMS算法中的步長因子，λ為常量，取值范圍為 0<入<1; S45:根據(jù)MSE最小均分誤差準(zhǔn)則，得到方向權(quán)向量最優(yōu)解6.根據(jù)權(quán)利要求1所述的一種基于對角加載的自適應(yīng)波束形成算法，其特征在于:所述 步驟5中：將得到的方向權(quán)值對采樣信號數(shù)據(jù)進(jìn)行加權(quán)求和，得到自適應(yīng)波束信號y化）= W opt X(k) 〇
【文檔編號】G01S7/539GK106093920SQ201610538089
【公開日】2016年11月9日
【申請日】2016年7月9日
【發(fā)明人】鄧正宏, 李學(xué)強(qiáng), 黃杰, 黃一杰, 付明月, 馬春苗
【申請人】西北工業(yè)大學(xué)