孫華普+鄭旺+呂衛(wèi)華+張小飛

摘要： 通信領(lǐng)域中非圓（NC）信號的波達方向 （DOA）估計已經(jīng)得到廣泛的研究和應(yīng)用。 考慮到非圓信號的多重信號分類算法（MUSIC）計算復(fù)雜度較高， 提出線陣中基于降維NC-MUSIC的非圓信號DOA估計算法。 與傳統(tǒng)的需要二維搜索的2D-NC-MUSIC算法相比， 該算法只需要一維搜索， 降低了計算復(fù)雜度。 仿真結(jié)果表明， 該算法優(yōu)于傳統(tǒng)的MUSIC算法， 并且與2D-NC-MUSIC算法的估計性能非常接近。

關(guān)鍵詞： 線性陣列； 非圓信號； MUSIC； DOA估計

中圖分類號： TN911.7文獻標(biāo)識碼： A文章編號： 1673-5048（2016）04-0036-06

Abstract： Direction of arrival （DOA） estimation for noncircular （NC） signals has been widely studied and used in communications. Considering the high computational complexity of NC multiple signal classification （MUSIC） algorithm， reduceddimension NCMUSIC algorithm is proposed for DOA estimation of NC signals. Compared with conventional 2DNCMUSIC algorithm which needs twodimension search， the proposed algorithm only requires onedimension search， which reduces the computational complexity. Simulation results illustrate that the proposed algorithm outperforms the conventional MUSIC algorithm and has close estimation performance to the 2DNCMUSIC algorithm.

Key words： linear array； NC signal； MUSIC； DOA estimation

0引言

在陣列信號處理中， DOA估計是一個主要的研究方向， 已廣泛應(yīng)用于通信、 雷達、 聲納和醫(yī)學(xué)圖像等領(lǐng)域[1-4]。 傳統(tǒng)的DOA估計算法包括MUSIC算法[5] 、 ESPRIT算法[6-8] 、 Capon算法[9-10]及PM算法[11]等， 均已經(jīng)有了較為成熟的理論， 但估計性能卻越來越難以滿足要求。 為了提高估計性能， 開始利用信號的非圓特征估計信號的DOA[12-19]。 非圓信號是指在星座圖中振幅只有同相分量而正交分量為零的信號， 諸如二進制相移鍵控（Binary Phase Shift Keying， BPSK）和調(diào)幅（Amplitude Modulation， AM）信號等 [20]。 當(dāng)信源發(fā)出非圓信號時， 利用非圓特性， 可以有效地將接收的數(shù)據(jù)矩陣維數(shù)加倍， 從而在提高參數(shù)估計性能的同時還能估計出更多的信號源。 因此， 非圓信號的非圓特性已被廣泛用于提高DOA估計性能。 文獻[12]提出一種非圓MUSIC算法實現(xiàn)DOA估計（NC-MUSIC）。 文獻[13]提出非圓信號求根MUSIC算法（Root-NC-MUSIC）， 避免了譜峰搜索。 文獻[14]提出一種用ESPRIT實現(xiàn)非圓信號的定位算法， 該算法（NC-ESPRIT）無需譜峰搜索， 復(fù)雜度較低。 文獻[17]則對非圓信號DOA估計提出一種非圓PM算法， 該算法（NC-PM）比傳統(tǒng)的PM算法[12]有著更好的性能。

上述算法的計算復(fù)雜度均較高， 尤其是文獻[12]提出的2D-NC-MUSIC算法， 需要對DOA和非圓相位進行二維搜索， 計算復(fù)雜度較高， 難以應(yīng)用到實際中。

因此， 為了降低計算復(fù)雜度， 在文獻[21]面陣降維MUSIC算法的啟發(fā)下， 結(jié)合2D-NC-MUSIC算法二維搜索的特征， 提出了降維NC-MUSIC算法。 仿真結(jié)果證明了降維NC-MUSIC算法的有效性。

本文其余部分的結(jié)構(gòu)如下： 第2節(jié)介紹了信號的數(shù)據(jù)模型； 第3節(jié)提出了降維NC-MUSIC算法； 第4節(jié)對算法進行性能分析； 在第5節(jié)中， 仿真結(jié)果證明了算法的有效性； 第6節(jié)為本文的總結(jié)。

1數(shù)據(jù)模型

一個有M個陣元的線性陣列如圖1所示。

降維NC-MUSIC算法的主要步驟如下：

（1） 通過式（4）擴展信號矩陣；

（2） 通過式（28）構(gòu)造擴展信號的協(xié)方差矩陣R^， 并對R^進行特征分解得到噪聲子空間U^N， 計算復(fù)雜度為O（4M2L+8M3）；

（3） 通過式（29）進行譜峰搜索， 估計信號的DOA， 計算復(fù)雜度為O（8M3-4M2K+8（M2+M）n）。

3性能分析

3.1復(fù)雜度分析

本文提出的降維NC-MUSIC算法的計算復(fù)雜度為O（4M2L+16M3-4M2K+8（M2+M）n）， 而2D-NC-MUSIC算法的計算復(fù)雜度為O（4M2L+16M3-4M2K+（4M2+2M）n2）， 其中n為搜索次數(shù)。 因此， 降維NC-MUSIC算法比2D-NC-MUSIC算法的計算復(fù)雜度低。

在K=3， n=6 000， M=8的情況下， 兩種算法的復(fù)雜度比較如圖2所示。 可以看出， 降維NC-MUSIC算法復(fù)雜度較低。

3.2主要優(yōu)點

本文提出的降維NC-MUSIC算法有如下優(yōu)點：

（1） 與傳統(tǒng)的MUSIC算法相比， 降維NC-MUSIC算法有較好的估計性能， 同時能夠估計更多的信號；

（2） 與2D-NC-MUSIC算法相比， 降維NC-MUSIC算法計算復(fù)雜度更低；

（3） 與2D-NC-MUSIC算法相比， 降維NC-MUSIC算法擁有相似的估計性能；

（4） 降維NC-MUSIC算法不需要估計非圓相位， 就能很好地估計DOA。

上述討論中， 假設(shè)非圓信號個數(shù)K已知。 假設(shè)在K未知的情況下， 也能通過信息論中的算法、 矩陣分解算法、 平滑排位算法或者蓋氏圓盤算法估計信號個數(shù)K[24-27]。

4仿真結(jié)果

通過仿真來評估降維NC-MUSIC算法的角度估計性能。

5結(jié)論

本文提出了降維NC-MUSIC算法估計非圓信號的DOA。 與傳統(tǒng)的MUSIC算法相比， 該算法擁有更好的估計性能， 且能夠估計更多的信號。 與2D-NC-MUSIC算法相比， 該算法只需要一維搜索， 在保持估計性能的同時， 有著較低的計算復(fù)雜度。 仿真結(jié)果證明了算法的有效性。

參考文獻：

[1] Krim H， Viberg M. Two Decades of Array Signal Processing Research： The Parametric Approach[J]. IEEE Signal Processing Magazine， 1996， 13（4）： 67-94.

[2] Applebaum S R. Adaptive Arrays[J].IEEE Transactions on Antennas and Propagation， 1976， 24（5）： 585-598.

[3] Li Jianfeng， Zhang Xiaofei. ClosedForm Blind 2DDOD and 2DDOA Estimation for MIMO Radar with Arbitrary Arrays[J]. Wireless Personal Communications， 2013， 69（1）： 175-186.

[4] Bekkerman I， Tabrikian J.Target Detection and Localization Using MIMO Radars and Sonars[J].IEEE Transactions on Signal Processing， 2006， 54（10）： 3873-3883.

[5] Schmidt R O. Multiple Emitter Location and Signal Parameter Estimation[J].IEEE Transactions on Antennas and Propagation， 1986， 34（3）： 276-280.

[6] Zheng Guimei， Chen Baixiao， Yang Minglei. Unitary ESPRIT Algorithm for Bistatic MIMO Radar [J]. Electronics Letters， 2012， 48（3）： 179-181.

[7] Bencheikh M L， Wang Y. Joint DODDOA Estimation Using Combined ESPRITMUSIC Approach in MIMO Radar [J]. Electronics Letters， 2010， 46（15）： 1081-1083.

[8] Zoltowski M D， Haardt M， Mathews C P. ClosedForm 2D Angle Estimation with Rectangular Arrays in Element Space or Beamspace via Unitary ESPRIT [J]. IEEE Transactions on Signal Processing， 1996， 44（1）： 326-328.

[9] 張保鋒. 幾種DOA估計算法的性能分析[J]. 現(xiàn)代電子技術(shù)， 2003（9）： 35-37.

[10] Stoica P， Hndel P， Sderstrm T. Study of Capon Method for Array Signal Processing[J]. Circuits， Systems and Signal Processing， 1995， 14（6）： 749-770.

[11] Marcos S， Marsal A， Benidir M. The Propagator Method for Source Bearing Estimation[J].Signal Processing， 1995， 42（2）： 121-138.

[12] Abeida H， Delmas J P. MUSICLike Estimation of Direction of Arrival for Noncircular Sources[J].IEEE Transactions on Signal Processing， 2006， 54（7）： 2678-2690.

[13] Chargé P， Wang Y D， Saillard J. A Noncircular Sources Direction Finding Method Using Polynomial Rooting[J].Signal Processing， 2001， 81（8）： 1765-1770.

[14] Zoubir A， Chargé P， Wang Yide. Non Circular Sources Localization with ESPRIT[C]∥Proceedings of European Conference on Wireless Technology， Munich， Germany， 2003.

[15] Steinwandt J， Roemer F， Haardt M. Performance Analysis of ESPRITType Algorithms for NonCircular Sources[C]∥International Conference on Acoustics， Speech， and Signal Processing， Vancouver， BC， 2013.

[16] Haardt M， Rmer F. Enhancements of Unitary ESPRIT for NonCircular Sources[C]∥Proceedings of 29th IEEE International Conference on Acoustics， Speech， and Signal Processing， Montreal， Quebec， Canada， 2004.

[17] 孫心宇， 周建江. 非圓傳播算子DOA估計算法[J]. 數(shù)據(jù)采集與處理， 2013， 28（3）： 313-318.

[18] Zhang Xiaofei， Cao Renzheng， Zhou Ming. NoncircularPARAFAC for 2DDOA Estimation of Noncircular Signals in Arbitrarily Spaced Acoustic VectorSensor Array Subjected to Unknown Locations[J]. EURASIP Journal on Advances in Signal Processing， 2013（1）： 1-10.

[19] Liu J， Huang Z T， Zhou Y Y. Azimuth and Elevation Estimation for Noncircular Signals[J]. Electronic Letters， 2007， 43（20）： 1117-1119.

[20] Delmas J P， Abeida H. CramerRao Bounds of DOA Estimates for BPSK and QPSK Modulated Signals[J]. IEEE Transactions on Signal Processing， 2006， 54（1）： 117-126.

[21] Zhang Xiaofei， Xu Lingyun， Xu Lei， et al. Direction of Departure （DOD） and Direction of Arrival （DOA） Estimation in MIMO Radar with ReducedDimension MUSIC[J]. IEEE Communications Letters， 2010， 14（12）： 1161-1163.

[22] 張小飛，汪飛，陳華偉. 陣列信號處理的理論與應(yīng)用[M].北京： 國防工業(yè)出版社， 2013.

[23] Zhong Manli， Fan Zheyi. DirectionofArrival Estimation for Noncircular Signals[C]∥Proceedings of the International Conference on Computer Networks & Communication Engineering， 2013， 30： 634-637.

[24] Wax M， Kailath T. Detection of the Signals by Information Theoretic Criteria[J]. IEEE Transactions on Acoustics， Speech， and Signal Processing， 1985， 33（2）： 387-392.

[25] Di Angzhao. Multiple Sources Location——A Matrix Decomposition Approach[J]. IEEE Transactions on Acoustics， Speech， and Signal Processing， 1985， 35（4）： 1086-1091.

[26] Huang Lei， Long Teng， Mao Erke， et al. MMSEBased MDL Method for Robust Estimation of Number of Sources without Eigendecomposition[J].IEEE Transactions on Signal Processing， 2009， 57（10）： 4135-4142.

[27] Wu H T， Yang J F， Chen F K. Source Number Estimator Using Gerschgorin Disks[C]∥Proceedings of ICASSP， Adelaide， Australia， 1994.

[9] Marcos S， Marsal A， Benidir M. The Propagator Method for Source Bearing Estimation [J].Signal Processing， 1995， 42（2）： 121-138.