一種基于非一致響應圓陣的DOA估計改進算法?

周興建，黃興莉，邵懷宗

（1.中國西南電子技術研究所，成都610036；2.電子科技大學通信與信息工程學院，成都611731）

一種基于非一致響應圓陣的DOA估計改進算法?

周興建1，??，黃興莉2，邵懷宗2

（1.中國西南電子技術研究所，成都610036；2.電子科技大學通信與信息工程學院，成都611731）

由于傳統(tǒng)波達方向（DOA）估計算法對接收天線陣列的陣元一致性要求高，即幅度和相位響應要求一致，在工程應用中受到較大限制。針對該問題，結合均勻圓陣的特點，通過引入陣元幅度和相位不一致因子對傳統(tǒng)MUSIC算法的模型加以改進，提出了一種改進的MUSIC算法。該算法降低了陣元的一致性要求，可通過迭代求解的方法來完成對多個信源來波方向進行估計。仿真和分析結果表明，提出的算法可有效地進行多個目標的DOA估計。

均勻圓陣；DOA估計；MUSIC算法；迭代處理

1 引言

基于MUSIC算法［1－2］的信號測向與多信號檢測在無線電頻譜監(jiān)測和管理中得到了廣泛的應用。過去30年中，在同一空域中利用基于天線陣列MUSIC算法來對多個目標信號的來波方向進行估計的算法得到了國內(nèi)外學者的廣泛研究。雖然計算機仿真和一些實驗系統(tǒng)的測試結果均表明：基于特征分解的MUSIC算法波達方向估計與傳統(tǒng)的波達方向

估計技術相比具有更好的性能和分辨率，但在實際中仍然存在很多因素限制其工程應用，一個主要原因是系統(tǒng)的校正［3－5］，由于組成陣列需要多個傳感器才能完成，多個傳感器的及每個傳感器到達中心處理器之間的通路一致性都需校正，而基于MUSIC算法的基本假設是均勻圓陣或均勻線陣，即組成陣列的傳感器位置是精確的。但現(xiàn)實情況是傳感器位置有誤差且需校正處理，其中對通路校正和傳感器位置校正的研究［2－3，6］較多，而對于傳感器一致性校正的研究很少有報道，就問題本身而言，傳感器要精確校正在工程中有較大的困難。而通常的MUSIC算法中，組成陣列的天線元的一致性的好壞，即其幅度響應和相位響應是否一致，是基于本征分解的傳統(tǒng)MUSIC類算法能否正常工作的關鍵。關于這方面的工作，雖然A.Paulraj等［2］人對未知增益和相位的陣列進行研究，但是這些研究中做了一個基本的假設，即在整個陣列中不同陣元具有不同的增益，而對每個獨立陣元來說，其增益是全向一致的，該假設限制了其應用環(huán)境。

本文結合工程實際，對常規(guī)MUSIC算法的模型進行了改進，得到一種新的迭代算法。該算法可實現(xiàn)在傳感器陣列各個陣元之間幅度和相位非一致情況下的DOA（Direction of Arrival）估計。

2 算法模型

假設信源個數(shù)為N，天線數(shù)目為M，第n個信源發(fā)射的信號為sn（t），那么第m個天線接收到的信號為

其中，amn、φmn分別表示第m個天線對第n個信源的增益和相位；τmn是第m個天線相對第n個信源的時間延遲，滿足以下關系：

其中，dmn為第m個天線和第n個信號源之間的距離；xm、ym分別為第m個天線在坐標系上的橫坐標和縱坐標；θn是第n個信號源的入射方位角，也是需要估計的參數(shù)。如果以圓點表示天線，天線對信源的增益以及方向通過封閉曲線來表示，那么不同的amn、φmn可以將圓陣分為以下3種情況。

（1）理想的陣列

如圖1（a）所示，天線是各向一致的且天線陣列流行是全向性的，即各天線的增益和相位響應完全一致。這種理想情況下的接收天線中的參數(shù)amn、φmn可分別用常數(shù)來表示，即

（2）陣元增益非一致陣列

即各個天線陣元的增益是各向一致的，但不同的陣元之間的增益和相位是有差異的，如圖1（b）所示。此時的參數(shù)可表示為

（3）不一致陣列

如圖1（c）所示，組成陣列的每個天線陣元的響應隨著角度的不同其增益是一致的，自然其天線的陣列流行不是全向性的。

此時，對同一個天線而言，不同信源的增益和相位也是不相等的，即

其中，Am、φm是圖中虛圓形表示的天線的基本增益和相位；αmn、βmn分別是增益和相位的衰減值；當αmn＝1、βmn＝0時表示天線的陣列流形是全向性，不存在增益和相位衰減。從式（8）、（9）中可以看出，相對于情況二，此時的某個天線對信號源的增益和相位在不同方向都出現(xiàn)了不同的衰減，正是這種衰減使得天線的陣列流形不是全向性的，其DOA的估計也與全向性條件下的估計不一樣。

圖1 五元陣的3種情況比較Fig.1 Comparison of five-element array under three conditions

3 改進的MUSIC算法

對于窄帶信號而言，對觀察的數(shù)據(jù)進行傅里葉變換后可得

其中，Sn是第n個信源的傅里葉變換，Vm是噪聲的傅里葉變換。假設快拍數(shù)為K，則

式中，k＝0，1，…，K－1。

其矩陣可以表示為X（k）＝Γ*A*S（k）＋V（k），k＝0，1…K－1（12）其中，Γ＝diag｛a1e－jω0φ1，a2e－jω0φ2，…aMe－jω0φM｝表示天線基本的增益和相位（圖中虛線圓形所示），Amn＝αmne－jω0βmne－jω0τmn表示非全向性陣列流形的方位矩陣。根據(jù)式（2）和式（3）可知

因此，需要估計的參數(shù)是θn，即信號的波達方向。但從上式中可以看到，接收信號有兩個未知矩陣，一是Γ矩陣，另一個是A矩陣，A矩陣不僅包含了信號的波達方向，也包括了陣元對信源的增益和相位的衰減信息。根據(jù)文獻［2］，采用以下方法進行波達方向的估計：

（1）DOA的估計，采用基本的MUSIC算法進行估計；

（2）相位和增益的估計，將波達方向固定，最小化代價函數(shù)：

其中，δ＝［Γ11，Γ22，…，ΓMM］T，Q（n）＝diag｛a（θn）｝，a（θn）是利用MUSIC算法估計出來的波達方向。利用文獻的準則簡化代價函數(shù)，可得相位和增益估計的方法為

改進的MUSIC算法的實現(xiàn)可分為如下幾個步驟。

Step 1：參數(shù)初始化，迭代次數(shù)j＝0；

Step 2：計算信號的自相關矩陣~Rx＝

Step 3：對相關矩陣~Rx進行特征值分解，并根據(jù)特征值的大小關系，選擇相應的特征向量構造信號子空間＾U；

Step 4：估計信號的DOA：

Step 5：估計相位和增益：

Step 6：計算代價函數(shù)Jjc＝δHZjδ，若代價函數(shù)小于閾值，則停止計算，否則令j＝j＋1，并回到Step 4。

4 仿真與分析

為了驗證改進算法的性能，本節(jié)對其進行了相應的仿真。仿真采用均勻圓陣，其中信源個數(shù)為2，天線數(shù)目為5。其他參數(shù)包括：信號源入射角θ1、θ2；Γ矩陣，包括天線基本的增益和相位矩陣，分別用T-gain和T-phase來表示；天線增益和相位的衰減因子α、β。

仿真1：θ1＝－15°，θ2＝5°，T-gain＝［10.58 0.82151.0770.7970］；T-phase＝［0－3.8

14.63 4.4－3.4］°；α＝1，β＝0。不同迭代次數(shù)的仿真結果如圖2所示。其中圖2（a）是傳統(tǒng)MU

SIC算法，其估計出來的DOA存在很大的誤差，而圖2（b）、（c）、（d）是采用改進后的迭代算法，可以看到，隨著迭代次數(shù)的增加，DOA估計越來越準確。

圖2 仿真1的仿真結果圖Fig.2 Comparison of DOA precision under case 1

仿真2：θ1＝－15°，θ2＝5°，T-gain＝［10.58 0.82151.0770.7970］；T-phase＝［0－3.8 14.634.4－3.4］；α是0.9～1之間且服從均勻分布的隨機數(shù)，β＝0。仿真結果如圖3所示。從圖3中可以看出，在僅存在小幅度增益衰減的情況下，經(jīng)過多次迭加之后的DOA估計越來越準確。

圖3 仿真2的仿真結果圖Fig.3 Comparison of DOA precision under case 2

仿真3：θ1＝－15°，θ2＝5°，T-gain＝［10.58 0.82151.0770.7970］；T-phase＝［0－3.8

14.63 4.4－3.4］；α是0.9～1之間且服從均勻分布的隨機數(shù)，β是0～7°之間且服從均勻分布的隨機相位。仿真結果如圖4所示。從圖4中可以看出，在存在小幅度相位和增益衰減時，增加迭代次數(shù)可使準確地估計出DOA。

圖4 仿真3的仿真結果圖Fig.4 Comparison of DOA precision under case 3

仿真實驗4：θ1＝－15°，θ2＝－10°，T-gain＝［1 0.580.82151.0770.7970］；T-phase＝［0－3.814.634.4－3.4］；α＝1，β＝0。仿真結果如圖5所示。從圖5可以看出，當兩個信源之間的入射角相差較小的時候，傳統(tǒng)MUSIC算法的DOA估計將存在非常大的誤差，幾乎不能工作，迭代算法可有效完成DOA的估計。

圖5 仿真4的仿真結果圖Fig.5 Comparison of DOA precision under case 4

5 總結

針對傳統(tǒng)MUSIC算法在DOA估計中對接收天線陣列的陣元一致性要求高，與實際工程應用存在差距的問題，文中結合均勻圓陣的特點，通過引入陣元幅度和相位不一致因子，提出了一種改進的迭代求解新算法。該算法可降低天線陣元的一致性要求，與實際應用能緊密結合，可實現(xiàn)對多個信源的來波方向進行估計。文中對提出的算法進行了不同情況下的仿真，表明該算法在天線陣列增益及相位不一致的情況下，可有效地完成對多個目標DOA的有效估計。

［1］He Jin，Swamy N S，Ahmad M O.Efficient application of music algorithm under the coexistence of far-field and nearfield sources［J］.IEEE Transactions on Signal Processing，2012，60（4）：2066－2070.

［2］Paulraj A，Kailath T.Direction of arrival estimation by eigenstructure methods with unknown sensor gain and phase［C］／／Proceedings of 1985 IEEE ICASSP.Tampa：IEEE，1985：640－643.

［3］Friedlander B，Weiss A J.Direction finding in the presence of mutual coupling［J］.IEEE Transactions on Antennas and Propagation，1991，39（3）：273－284.

［4］Nemirovsky S，Doron M A.Sensitivity of MUSIC and Root-MUSIC to gain calibration errors of 2D arbitrary array configuration［C］／／Proceedings of Sensor Array and Multichannel Signal Processing Workshop.IEEE，2004：594－598.

［5］Ye Zhongfu，Liu Chao.On the resiliency of MUSIC direction finding against antenna sensor coupling［J］.IEEE Transactions on Antennas and Propagation，2008，56（2）：371－380.

［6］Friedlander B.A sensitivity analysis of the MUSIC algorithm［J］.IEEE Transactions on Acoustics Speech，and Signal Processing，1990，38（10）：1740－1751.

ZHOU Xing-jian was born in Guang′an，Sichuan Province，in 1970.He is now a senior engineer.His research concerns communication and data-link technology.

Email：zxjzjy＠163.com

黃興莉（1986—），女，四川內(nèi)江人，碩士研究生，主要從事寬帶無線通信信號處理方面的研究工作；

HUANG Xing-li was born in Neijiang，Sichuan Province，in 1986.She is now a graduate student.Her research concerns wireless communication and signal processing.

邵懷宗（1969—），男，四川巴中人，電子科技大學副教授、碩士生導師，主要從事通信信號處理、陣列信號處理和DSP技術方面的科研和教學工作。

SHAO Huai-zong was born in Bazhong，Sichuan Province，in 1969.He is now an associate professor.His research concerns communication signal processing and DSP technology.

An Improved DOA Estimation Algorithm Based on Nonuniform Response Antenna Circular Arrays

ZHOU Xing-jian1，HUANG Xing-li2，SHAO Huai-zong2
（1.Southwest China Institute of Electronic Technology，Chengdu 610036，China；2.Institute of Communication and Information Engineering，University of Electronic Science and Technology of China，Chengdu 611731，China）

The traditional direction of arrival（DOA）estimation algorithm has high conformity with array amplitude and phase element of receiving arrays.Consequently it is limited seriously in engineering application.In this paper，an improved model of MUSIC algorithm is proposed based on the variances of uniform circular array element′s amplitude and phase，and it reduces the request of the conformity of the array element.The improved algorithm can be utilized to estimate the direction of the receiving waves from multiple information sources through iteration procedure.The simulation results indicate that the proposed algorithm can estimate the DOA from different targets effectively.

uniform circular array；DOA estimation；MUSIC algorithm；iteration calculation

date：2013－03－02；Revised date：2013－07－15

??通訊作者：zxjzjy＠163.comCorresponding author：zxjzjy＠163.com

TN97

A

1001－893X（2013）09－1164－05

周興建（1970—），男，四川廣安人，高級工程師，主要研究方向為通信及數(shù)據(jù)鏈技術；

10.3969／j.issn.1001－893x.2013.09.009

2013－03－02；

2013－07－15