一種穩(wěn)健的自適應(yīng)波束形成方法

1、收稿日期 :2004-04-12; 修回日期 :2004-08-05。作者簡(jiǎn)介 :楊莘元 (1944- , 男 , 教授 , 博士生導(dǎo)師 , 主要研究方向?yàn)橥ㄐ判盘?hào)處理 , 自適應(yīng)陣列信號(hào)處理 , 多目標(biāo)檢測(cè) , 跟蹤及抗干擾技術(shù)。一種穩(wěn)健的自適應(yīng)波束形成方法楊莘元 , 謝 紅 , 陳四根(哈爾濱工程大學(xué)信息與通信工程學(xué)院 , 黑龍江 哈爾濱 150001摘 要 :基于線性約束最小方差 (linearly constrained minimum variance , LC M V 波束形成原理 , 根據(jù)貝葉斯估計(jì) 理論 , 提出了一種基于加權(quán)的波束形成方法 。 該波束形成方法能夠很好地克服 L

2、C M V 波束形成對(duì)期望信號(hào) DOA誤差特別敏感的弱點(diǎn) 。 仿真試驗(yàn)表明 , 該方法性能穩(wěn)定 , 對(duì)期望信號(hào) DOA 誤差具有很好的穩(wěn)健性 。關(guān)鍵詞 :陣列信號(hào)處理 ; 自適應(yīng)波束形成 ; 線性約束最小方差 中圖分類(lèi)號(hào) :T N957. 52 文獻(xiàn)標(biāo)識(shí)碼 :AR obust adaptive beam Y 2, , N Si 2gen(College o f Information , Harbin Engineering Univer sity , Harbin 150001, China :linearly constrained m inimum variance (LC M V ad

3、aptive beam form ing , an adaptive beam by using Bayesian posterior probability weighted sum is proposed in terms of the Bayesian estimation theory. method is developed to im prove robustness to pointing error and overcome the shortcom ing of LC M V beam form ing s sensitivity to pointing error. The

4、 results of simulation indicate that the performance of the pro 2posed beam form ing method is robust to the pointing error.K ey w ords :array signal processing ; adaptive beam form ing ; linearly constrained m inimum variance1 引 言在陣列信號(hào)處理中 , 自適應(yīng)波束形成是一種用于估計(jì)未知信號(hào)的有效方法 。線性約束最小方差 (linearly con 2strained

5、minimum variance , LC M V 波束形成技術(shù)在理想的情 況下能夠提供沒(méi)有畸變的期望信號(hào)輸出 , 同時(shí)能夠最大程 度地抑制干擾和噪聲 。 但是 ,LC M V 波束形成過(guò)于依賴于期 望信號(hào)的來(lái)波方向 , 對(duì)來(lái)波方向的誤差特別敏感 1, 如果 LC M V 波束形成中所利用的觀察方向和實(shí)際期望信號(hào)來(lái)波 方向不完全匹配 ,LC M V 波束形成的性能就會(huì)急劇下降 , 而 且這種不匹配的情況經(jīng)常發(fā)生 , 如 DOA 估計(jì)的不精確 , 信 號(hào)源或者是天線陣列的移動(dòng)導(dǎo)致的來(lái)波角度的變化等等 。 研究人員提出了很多方法來(lái)克服這一對(duì)指向誤差的敏感問(wèn) 題 , 主要有兩類(lèi) 2, 一類(lèi)是通過(guò)圍繞

6、假設(shè)的 DOA 拓寬陣列 天線方向圖的主波束 , 另一類(lèi)是在實(shí)時(shí)的環(huán)境中 , 通過(guò)學(xué)習(xí) 和估計(jì)的方法從觀察的數(shù)據(jù)中獲得信號(hào)的 DOA , 以提高波 束形成的穩(wěn)健性 , 在低信噪比和信號(hào)來(lái)波波動(dòng)較快的情況 下 , 這種波束形成的性能變得很差 3。本文基于 LC M V 波束形成準(zhǔn)則 , 根據(jù)貝葉斯方法推導(dǎo) 出一種穩(wěn)健的波束形成方法 。 該方法對(duì)于來(lái)波方向的誤差 具有很好的穩(wěn)健性 。2 信號(hào)模型及波束形成原理211 信號(hào)模型設(shè)有 q 個(gè)信號(hào)入射到天線陣上 , 入射角度分別是 1, 2, , q , 天線陣由 M 個(gè)陣元組成 。假設(shè)天線陣和信號(hào)源 位于同一個(gè)平面中 , 所有的信號(hào)相對(duì)于天線陣為遠(yuǎn)場(chǎng)窄

7、帶 信號(hào) , 具有相同的中心頻率 。 令 s i (t 表示第 i 個(gè)信號(hào)的復(fù) 包絡(luò) , x i (t 表示為第 i 個(gè)陣元所接收信號(hào)的復(fù)包絡(luò) , 則X (t =x 1(t , x 2(t , , x M (t T為天線陣接收信號(hào)的復(fù)包絡(luò)矢量 ,T 為矩陣的轉(zhuǎn)置運(yùn)算 。在加性噪聲的情況下 , 陣 列接收矢量可以表示為X (t = qi =1a (is i(t +N (t (1式中 :a (i DOA 為 i 的信號(hào)導(dǎo)向矢量 ; N (t 陣 列的接收噪聲矢量 。 式 (1 的矩陣表示為X (t =A (S (t +N (t (2 =1, 2, , q T, A ( =a (1 , a (2 ,

8、 ,a (q T, S (t =s 1(t , s 2(t , , s q (t T, 假設(shè)信號(hào)干擾以及噪聲之間互不相關(guān) 。 則接收信號(hào)的相關(guān)矩陣為R X =EXX H =AR s A H +2n I(3式中 :R s =ESS H 信號(hào)相關(guān)矩陣 ,H 矩陣的共軛2005年 2月 系統(tǒng)工程與電子技術(shù)Feb. 2005第 27卷 第 2期Systems Engineering and Electronics V ol. 27 N o. 2文章編號(hào) :10012506X (2005 0220244203轉(zhuǎn)置運(yùn)算 , 2n, I 噪聲的方差和單位矩陣 。2. 2 波束形成原理 6波束形成的輸出表示為

9、 y (t =W H X (t , 當(dāng)期望信號(hào) 的 DOA 已知時(shí) , 對(duì)于 LC M V 波束形成就是要選擇一個(gè)權(quán) , 使 得陣列輸出功率最小 , 同時(shí)滿足信號(hào)方向的增益為常數(shù)的 約束條件 , 即為minWW H R X Wa (0 H W =1(4 式中 :0 期望信號(hào)的 DOA 。 LC M V 波束形成的權(quán)矢量為 W =R -1X a (0 /a (0 H R -1X a (0 (5 當(dāng) 0不能精確已知時(shí) , 一種方法可以用來(lái)降低這種波束形 成對(duì)指向誤差的敏感性 , 即通過(guò)強(qiáng)加約束條件使得陣列主 瓣變寬變平坦 , 這時(shí)權(quán)函數(shù)滿足的條件是minWW H R X WC H W =f式中 :

10、C M ×K 的矩陣 ,件 , f K ×1矢量 , ,權(quán)矢量為W R -X (C H R -1X C -1f (7 附加額外的約束條件可以提高波束形成對(duì)指向誤差的 穩(wěn)健性 , 但是這種做法同時(shí)也妨礙了噪聲的抑制能力 , 因?yàn)?它減少了自適應(yīng)的自由度 , 或者說(shuō)由于主瓣變寬 , 噪聲抑制 能力自然就降低 。另一方法是從已接收的數(shù)據(jù)中估計(jì)信號(hào)的 DOA , 0利 用估計(jì)值 0去計(jì)算波束形成的權(quán)矢量W =R -1X a (0 /a (0 H R -1X a (0 (8 這一技術(shù)要求對(duì) 0的估計(jì)具有很好的性能 , 估計(jì)誤差大顯 然會(huì)導(dǎo)致波束形成性能的急劇降低 。3 穩(wěn)健的波束形成

11、方法假設(shè) 0 是一個(gè)先驗(yàn)概率密度函數(shù) (PDF 為 p ( 的 隨機(jī)變量參數(shù) , p ( 反應(yīng)了信號(hào)源 DOA 的不確定性 , 一般 假設(shè)信號(hào) DOA 只在一個(gè)區(qū)間內(nèi)離散的取值 , 即 =01,02, , 0m , 這區(qū)間包含有 0, 但是并不包含干擾信號(hào)的 DOA , 一般是以信號(hào) DOA 估計(jì)值為中心 , 以陣列天線的主波 束寬度為取值區(qū)間來(lái)定義 , 取值的個(gè)數(shù)決定了計(jì)算的復(fù)雜度和角度分辨率 。 令 XL表示 L 個(gè)快拍的陣列接收信號(hào) 。期望信號(hào)的 M MSE 估計(jì)即是在給定 XL下的期望信號(hào)的條 件均值 4,5。y (t =E s 0(t |X L =E E s 0(t |X L , =

12、 mi =1p (i |X L E s 0(t |X L , i (9 式中 :E 最佳線性估計(jì)算子 , E s 0(t |X L , 在已知 XL, 的條件下對(duì) s 0(t 的估計(jì)值 , E 對(duì) 求數(shù)學(xué)期望 , p (i|X L 在給定觀察信號(hào)下的后驗(yàn)概率密度 , 根據(jù)貝葉斯公式可得p (i |X L =( ( mi =1p (i p (X L |i (10 當(dāng)信號(hào)為聯(lián)合高斯信號(hào)時(shí) , E s 0(t |XL, i 就是指向?yàn)?i的波束形成 , 所以 y (t = mi =1p (i |X L W H (i X (t 。 所以這種波束形成的權(quán)矢量為W = mi =1p (i |X L W H

13、 (i (11 這樣權(quán)矢量其實(shí)是一個(gè)后驗(yàn)概率的加權(quán)和 。在高斯假設(shè)的前提下 , p (XL|i 是一個(gè)高斯概率密度函數(shù) , 均值為 0, 協(xié)方差為 RX(i L iLl =M |X(i |-t l R -1X (i X (t l =-RX(i |-L exp Ll =1(-X H (t l R -1X (i X (t l (12 將陣列信號(hào)的相關(guān)矩陣寫(xiě)成如下形式R X =2s a (i H a (i +R IN根據(jù)矩陣行列式和矩陣求逆定理|R X |=|R IN |(1+2s a H (i R -1IN a (i (13 R -1X =R -1IN -2-1( H ( -11+2s a H (

14、i R -1IN a (i (14 所以p (i |X L =cp (i (1+2s (i Lexp2( 21+2s (i H ( -1-1(i2(15(i=a H (i R -1IN a (i , c 為與角度 i 無(wú)關(guān)的歸一化因 子 , 使得定義域內(nèi)的概率總和為 1。可以看到 , 當(dāng) S NR 很 高 , L 很大時(shí) , 這一后驗(yàn)概率密度將會(huì)在真實(shí)的 DOA 處等 于 1, 而其它值處為 0, 這時(shí)的波束形成就等價(jià)于一個(gè)空間 維納濾波器 。 如果當(dāng) S NR 很低 , L 較小時(shí) , 后驗(yàn)概率密度將 幾乎等于先驗(yàn)概率密度 , 這時(shí)波束形成表現(xiàn)為各個(gè)空間維 納濾波器的先驗(yàn)概率的加權(quán)平均 。因

15、為函數(shù)中 RIN未知 , 后驗(yàn)概率密度估計(jì)起來(lái)也十分 困難 。 下面介紹一種近似的后驗(yàn)概率密度估計(jì)方法 ?？紤]沒(méi)有干擾的情況 , 則 RIN=2n I , 所以2(1+2s (i =2n222n(1+M 2s /2n >(16 后驗(yàn)概率密度近似為p (i |X L cp (i exp L H ( -1-1(i2(17 此時(shí)指數(shù)項(xiàng)正好是理想的最小方差響應(yīng)不變 (M VDR 波束形成的輸出功率 , 并且此波束形成的指向?yàn)?i。在 的內(nèi)指數(shù)的作用是將這種輸出功率的大小關(guān)系進(jìn)行非線性的放大 。 由于 RIN未知 , 利用 LC M V 波束形成的輸出功率來(lái) 估計(jì)原后驗(yàn)概率密度函數(shù)中的指數(shù)項(xiàng) 第

16、27卷 第 2期 一種穩(wěn)健的自適應(yīng)波束形成方法 245 p (i |X L =cp (i exp L (a H (i R X a (i -1(18 這種估計(jì)是在預(yù)先得到的期望信號(hào)的 DOA 的較小范 圍內(nèi)進(jìn)行的 , 認(rèn)為在估計(jì)域內(nèi)不包含有干擾信號(hào) , 因此可以 利用此式估計(jì)后驗(yàn)概率密度 。4 性能對(duì)比分析為了進(jìn)一步說(shuō)明加權(quán)波束形成的性能 , 這里對(duì)典型的 LC M V 波束形成和本文所提出的后驗(yàn)概率加權(quán)波束形成的 性能進(jìn)行了仿真對(duì)比 。 仿真時(shí) , 天線陣陣元數(shù)為 10, 陣元間 隔為 015的載波波長(zhǎng) , 信號(hào)假設(shè)為遠(yuǎn)場(chǎng)窄帶信號(hào) , 快拍數(shù)為 1000。圖 1為兩種波束形成的陣列方向圖 ,

17、其中虛線表示的 是 LC M V 波束形成的陣列方向圖 , 實(shí)線表示的是加權(quán)波束 形成的方向圖 。 期望信號(hào)信噪比 S NR =0dB ,DOA 為 -; 有兩個(gè)同頻干擾 ,DOA 分別為 -40°,20=20dB 。 LC M V +1即是說(shuō)使用 LC M V 9°1中的虛線可以看 到 , , 但在實(shí)際的信號(hào) 方向上也有一個(gè) -30dB 左右的下陷 。 加權(quán)波束形成 -9°為 中心 , 對(duì) ±2. 5°的區(qū)間 , 以 015°的步長(zhǎng)進(jìn)行加權(quán)運(yùn)算 , 在實(shí)線 表示的方向圖中 , 在 -10°處陣列增益為 0dB , 同時(shí)在干

18、擾方 向上也有比 -40dB 還多的深零陷。 這說(shuō)明加權(quán)波束形成具 有比 LC M V 不可比擬的抗方向誤差的穩(wěn)健性。在圖 1試驗(yàn)相同的條件下 , 圖 2表示的是近似的后驗(yàn) 概率密度圖 , 可以看到 , 雖然由于預(yù)先得到的期望信號(hào)的 DOA 有 +1°誤差 , 但是在后驗(yàn)概率密度上仍然在實(shí)際的 DOA 上體現(xiàn)出高概率 ; 圖 3是陣列天線的部分輸出信號(hào) , 同 樣曲線 為加權(quán)波束形成 , 曲線 表示 LC M V 波束形成的 輸出信號(hào) , 對(duì)比可以看到前者明顯優(yōu)于后者 , 后者的輸出接 近噪聲狀態(tài) 。 圖 1 兩種波束形成的方向圖 圖 2 近似后驗(yàn)概率密度 其它條件都相同 , 實(shí)驗(yàn)還

19、對(duì)不同的角度誤差的情況下 兩者的輸出信噪比進(jìn)行了對(duì)比 , 誤差的范圍為 0°215°, 實(shí) 驗(yàn)步長(zhǎng)為 0105°, 信噪比對(duì)比曲線如圖 4。上方較為平坦的 是加權(quán)波束形成的輸出信噪比曲線 , 下方呈現(xiàn)明顯下降趨 勢(shì)的是 LC M V 波束形成的輸出信噪比曲線 , 再一次表明加 權(quán)波束形成對(duì)期望信號(hào)的 DOA 誤差具有很好的穩(wěn)定性 。 圖 4 兩種波束形成的 輸出信噪比5 結(jié) 論LC M V 波束形成是一種較為常見(jiàn)的陣列波束形成方法 ,但是該方法要求非常精準(zhǔn)地已知期望信號(hào)的 DOA , 對(duì)期望 信號(hào)的 DOA 的誤差十分敏感 , 本文推導(dǎo)了一種后驗(yàn)概率加 權(quán)的波束形

20、成方法 。假設(shè)已知一個(gè)期望信號(hào) DOA 粗略的 估計(jì)值 , 在一個(gè)離散的區(qū)間內(nèi) , 根據(jù)陣列接收到的數(shù)據(jù) , 近 似地算出區(qū)間內(nèi)各值的后驗(yàn)概率 , 然后加權(quán)綜合得到穩(wěn)健 的波束形成權(quán)矢量 。實(shí)驗(yàn)證明 , 這一方法能夠有效地克服 LC M V 波束形成的指向誤差敏感的問(wèn)題 。參考文獻(xiàn) :1C ox H. Res olving power and sensitivity to m ism atch of optimum array pro 2cess orJ.Acoustic S ociety Am erican , 1973, 54(9 :771-785. 2Bell KL , E phraim Y, Van T rees H L. R obust adaptive beam form ing us 2ing