鄭玉偉，張春良，岳 夏，劉 凱，朱厚耀

（1.南華大學(xué)機(jī)械工程學(xué)院，湖南衡陽 421000；2.廣州大學(xué)機(jī)械與電氣工程學(xué)院，廣東廣州 510006）

一種基于二次型約束的最優(yōu)對角加載量的計(jì)算方法*

鄭玉偉1，張春良2，岳 夏2，劉 凱1，朱厚耀2

（1.南華大學(xué)機(jī)械工程學(xué)院，湖南衡陽 421000；2.廣州大學(xué)機(jī)械與電氣工程學(xué)院，廣東廣州 510006）

對角線加載技術(shù)可以使自適應(yīng)波束形成器在失配情況下仍然能保持良好的性能，但是加載量的確定卻沒有一個(gè)嚴(yán)格的標(biāo)準(zhǔn)。采用一個(gè)二次型約束來討論對角加載問題。將加載水平同權(quán)值模的約束條件聯(lián)系起來，通過迭代的方法求出最優(yōu)加載水平。該方法不需要知曉信號的先驗(yàn)信息，其具有更好的靈活性與適用性。仿真結(jié)果表明：最優(yōu)加載在不同信噪比，較少快拍情況下均具有較好的穩(wěn)健性，且在模糊環(huán)境下更容易實(shí)現(xiàn)。

陣列信號處理；波束形成；DOA失配；對角加載；魯棒性

0 引言

波束形成技術(shù)廣泛應(yīng)用于聲納、雷達(dá)和通訊系統(tǒng)中，長久以來研究者們致力于優(yōu)秀算法的研究，其中MVDR、MUSIC和ESPRIT的方法不斷將波束形成技術(shù)推向更高的高度。在波束形成器的設(shè)計(jì)過程中，信號的向?qū)噶恳约翱沼蜃V矩陣均采用理想的模型，但是在現(xiàn)實(shí)使用中波束形成器是按照離散的步進(jìn)方式改變向?qū)噶窟M(jìn)行空域掃描，并且采用快拍采樣協(xié)方差的方法估計(jì)空域譜矩陣。這將會造成波束形成器的失配，除此之外還有陣列擾動，波達(dá)方向失配等失配方式，都將造成波束形成的性能損失。

眾多學(xué)者在增強(qiáng)波束形成器的魯棒性方面也做出了許多卓越的成就。COX等人提出了一種對角加載的方法[1]，使得波束形成器在失配時(shí)更加穩(wěn)健。Feldman等人提出了基于特征空間方法[2]，但是當(dāng)子空間維數(shù)不確定或者SNR較低時(shí)，該方法有較為嚴(yán)重的性能衰落。Gershman等人基于加寬波束零陷帶的思想，提出了獨(dú)立數(shù)據(jù)微分約束法[3]。Bell等人提出基于貝葉斯方法的自適應(yīng)算法，能在陣列接收信號與期望信號DOA的先驗(yàn)信息之間取得一種平衡。其中，對角加載方法是最為常見一種，但是加載量的確定一直是一個(gè)比較困難的問題。文獻(xiàn)[4]分析了對角加載對自適應(yīng)陣列SINR和INR性能的影響，研究加載量的門限選取問題，該方法是對預(yù)計(jì)的SINR和INR的水平有有足夠的信息，而選擇一個(gè)固定的適合的對角加載量。但是在現(xiàn)實(shí)中沒有足夠的先驗(yàn)信息來確定固定的加載量。文獻(xiàn)[5]研究波束畸變的根本原因，根據(jù)接收信號協(xié)方差矩陣的結(jié)構(gòu)特點(diǎn)，討論了小特征值λ與加載量的關(guān)系。但是在門限的選取和加載量的確定仍然沒有明確的說明。

1 基于二次型約束的對角加載

假設(shè)線性均勻陣列，陣元數(shù)為M，有K個(gè)不相干的信號源。

信號模型：

其 中 ， X(t)=[x0(t)，x1(t)，···，xM-1(t)]T，A=[a(θ1),a(θ2)，···,a(θk)]，a(θ)為陣列流矢量，N為背景噪聲。

陣列輸出：

其中，w=[w1,w2,···,wM]T，為陣列信號的加權(quán)矢量。

輸出功率譜：

其中R為接收信號的協(xié)方差矩陣R=E[x(t)xH(t)]，在現(xiàn)實(shí)實(shí)現(xiàn)中，通常采用快拍采樣估計(jì)其協(xié)方差矩陣

最小方差無畸變響應(yīng)（Mvdr）的最優(yōu)權(quán)矢量

所以針對失配問題，在陣列設(shè)計(jì)時(shí)，增加一個(gè)敏感度約束，即T≤T0。

則，二次型約束優(yōu)化問題可以表示為，

根據(jù)拉格朗日乘子算法解得

上式與參考文獻(xiàn)[7]中加載模型一致。同時(shí)可以看出，加載水平λ與T的選擇有關(guān)。在前面提到的文獻(xiàn)中，加權(quán)水平都是根據(jù)信號先驗(yàn)信息直接給出的，并沒有將其與加權(quán)的模的原約束值以及敏感度T聯(lián)系起來。

下面將通過二次型約束問題，求出最優(yōu)加載值與敏感度以及加權(quán)模值得關(guān)系。

在上式中，施加的一次約束條件，是約束期望信號無畸變的通過波束形成器，但是在現(xiàn)實(shí)實(shí)現(xiàn)中，接收信號是由可能的多個(gè)期望信號、干擾信號以及噪聲復(fù)合的信息，為了使全部的期望信號按照一定的準(zhǔn)則通過波束形成器，獲得最佳的波束形成，需要施加更多的約束條件。那么該問題的一次約束就轉(zhuǎn)化為了靈活度更大、適用性更廣的線性約束問題。下面將通過研究二次型線性約束問題，確定最佳的加載水平。

LCMV優(yōu)化問題：

一次約束條件問線性約束：

二次約束條件：

其中，C是一個(gè)N×M維的約束矩，g為約束值。例如，如果希望感興趣的K個(gè)信號無畸變地通過波束形成器，則此時(shí)的約束矩陣C=[a(θ1),a(θ2),···,a(θK)]，約束值g=[1,1,···,1]Tk×1。

首先研究T0所允許的取值范圍。

先解出等式約束問題，即wHw=T0。

線性約束的一個(gè)等效的實(shí)現(xiàn)結(jié)構(gòu)為廣義旁瓣對消器[7]。則把最優(yōu)解wopt表示為廣義旁瓣對消器的形式，

其中PC是wopt向約束子空間的投影PC=C[CTC]-1CH，P⊥C是wopt向約束空間的正交空間投影。

所以

故，有

因此，

即

所以設(shè)置的約束值不能小于gH(CHC)-1g?，F(xiàn)在解優(yōu)化問題，

由Lagrange乘子算法解出，

定義

并為了強(qiáng)調(diào)不等式約束關(guān)系，定義

有廣義旁瓣對消結(jié)構(gòu)，最優(yōu)權(quán)矢量可以表示如下：

其中，阻塞矩陣B滿足BHC=0。注意B并不是唯一的，構(gòu)造B的一個(gè)方法是P⊥C來確定，并且結(jié)果滿足BHB=I，

則

首先把式（22）寫為

其中，Rz=BHRB，是z的譜矩陣，z是阻塞矩陣的輸出，且pz=BHRwq，

當(dāng)λ≥0時(shí)，對角加載數(shù)據(jù)矩陣(R+λI)是正定的，在上式的中的導(dǎo)數(shù)值是負(fù)值。即λ＞0時(shí)，權(quán)矢量的模隨著λ增大而減小。

現(xiàn)在求滿足不等式的約束情況的λ：

或

接下來將使用數(shù)值方法求解對角加載量的與二次型約束之間的關(guān)系。

開始有λ0=0；令λ1=λ0+Δλ且有

算法總結(jié)：

（1）由二次型約束構(gòu)建線性約束問題；

（2）由約束條件確定二次約束的門限值T0；

（3）構(gòu)造阻塞矩陣B，使其滿足BHC=0, BHB=I，并且求出wq與a表達(dá)式；

（4）根據(jù)求得的公式和條件，使用迭代法求出最優(yōu)加載量λopt。

2 仿真分析

仿真1，采用16陣元均勻線陣，陣列間距為半波長，，背景噪聲為高斯白噪聲。DOA（波達(dá)方向）分別取-20，-10，0，20，30，60。期望的DOA為0，快拍次數(shù)為1 000。比較在T0條件下對角加載前后的無畸變約束下波束形成方向圖。這里取T0=0.15，因?yàn)闊o畸變約束下的白噪聲最大增益為均勻加權(quán)時(shí)的增益N，對于任意非均勻加權(quán)，都有

由圖1～4可以看出，在較低信噪比的情況下，陣列信號的擾動不是很明顯，加載前后的波束方向如沒有明顯變化，但是，隨著信噪比的提高，陣列信號失配逐漸增大，波束方向圖失真也更加嚴(yán)重，由圖4可以清楚的看到，加載前方向圖嚴(yán)重失真，幾乎無法獲得期望的DOA，而加載后的波束方向圖仍然具有良好的保形性。對失配環(huán)境具有較好的穩(wěn)健性。再者，比較不同信噪比條件下，加載后的波束方向圖，不管信噪比如何變化，加載后的波束方向圖都沒有太大的變化，所以在不知道信號先驗(yàn)信息的情況下，二次型約束最優(yōu)加載具有更好的適用性。

圖1 SNR=-10的波束方向圖

圖2 SNR=0的波束方向圖

圖3 SNR=10的波束方向圖

圖4 SNR=20的波束方向圖

圖5 SNR=-10情況下的波束方向圖

圖6 SNR=10情況下的波束方向圖

仿真2，取快拍次數(shù)N=40，其他仿真條件與仿真以相同，比較SNR=-10，10情況下，加載前后的波束方向圖。

如圖5所示，在快拍次數(shù)較少時(shí)，即使在較低的信噪比的情況下，加載前已不能較好的實(shí)現(xiàn)波束形成，而加載后的情況有明顯改善；如圖六示，并與圖3比較，在快拍次數(shù)較多的情況下，加載前后的波束方向圖相近，只是加載后的保形性更好，但是，在快拍較少的情況下，加載前的波束形成已經(jīng)失效，而加載后的波束方向圖依然保持較好的性能。

3 結(jié)論

通過以上的仿真分析可知，在估計(jì)協(xié)方差矩陣的采樣快拍次數(shù)有限，且信號的先驗(yàn)信息未知的情況下，通過二次型約束確定的最優(yōu)加載量，能夠有效的提高波束形成的穩(wěn)健性，并且由于其不需要知曉信號先驗(yàn)信息，所以該方法具有更廣泛的適用性和靈活性。

［1］H.COX Resolving power and sensitivity to mismatch of optimum array processors［J］.J.Acoust.Soc.Am.，vol.ASSP-54，pp.771-785，September 1973.

［2］Feldman D D，Griffith L J.A project approach for robust adaptive beam-forming［J］.IEEE Trans Signal Process? ing，1994（42）：867-876.

［3］ Gershman A B， Nickel U，Bohme J F.A adaptive beam-forming algorithms with robustness against jammer motion［J］.IEEE Trans Signal Processing， 1997（45）：1878-1885.

［4］陳曉初.自適應(yīng)陣對角加載研究［J］.電子學(xué)報(bào)，1998，26（4）：29-35.

［5］張小飛，徐大專.自適應(yīng)對角線加載波束形成算法［J］.中國空間科學(xué)技術(shù)，2007，4（2）：66-71.

［6］Harry L.Van Trees.Optimum Array Processing Part IV of Detection， Estimation，and Modulation Theory［M］. ISBN 978-7-302-14760-2.

［7］張小飛，汪飛，徐大專.陣列信號處理的理論與應(yīng)用［M］.北京：國防工業(yè)出版社，2010.

［8］謝斌斌.穩(wěn)健波束形成算法研究［D］.成都：電子科技大學(xué)，2012.

［9］蔣留兵，羅良桂，車?yán)?頻率不變約束的對角加載穩(wěn)健寬帶波束形成算法［J］.現(xiàn)代雷達(dá)，2012，12（34）：41-45.

［10］史英春，鐘子發(fā)，鄒翔，等.基于范數(shù)優(yōu)化的對角加載穩(wěn)健自適應(yīng)波束形成［J］.電路與系統(tǒng)學(xué)報(bào)，2013，18（1）：218-224.

A Method of Calculating the Optimal Diagonal Loading Based on a Quadratic Constraint

ZHENG Yu-wei1，ZHANG Chun-liang2，YUE Xia2，LIU Kai1，ZHU Hou-yao2
（1.School of Mechanical Engineering，University of South China，Hengyang421000，China；2.School of Mechanical and Electrical Engineering，Guangzhou University，Guangzhou510006，China）

Diagonal loading technique can make the adaptive beam-forming keep good performance in the case of mismatch，but there is not a strict standard in the determination of the loading quantity.This paper uses a quadratic constraint to discuss diagonal loading problems.The loading level will associate with the constraint conditions of weight module，and then calculate the optimal loading quantity by iteration.This method does not need to know the specific priori information of signal，so it is more flexible and applicative.Simulation results show that optimal load has better robustness in the case of both different SNR and less snap，and it is easier to implement in a fuzzy environment.

array signal processing；beam-forming；DOA mismatch；diagonal loading；robustness

TN957.51

A

1009－9492(2014)01－0021－05

10.3969/j.issn.1009-9492.2014.01.006

鄭玉偉，男，1989年生，安徽六安人，碩士研究生。研究領(lǐng)域：故障診斷。

(編輯：向 飛)

張春良，男，1964年生，教授，博士生導(dǎo)師。研究領(lǐng)域：設(shè)備狀態(tài)監(jiān)測與故障診斷，制造過程自動化，微制造技術(shù)，激光加工技術(shù)，數(shù)控技術(shù)等。

*國家自然科學(xué)基金（編號：51275099/E050302）；廣東省自然科學(xué)基金項(xiàng)目（編號：S2012010009505）

2013－07－12