劉曉娣,周新力,肖金光

(海軍航空工程學院電子信息工程系,山東 煙臺　264001)

?

基于空間平滑的單快拍波達方向估計算法

劉曉娣,周新力,肖金光

(海軍航空工程學院電子信息工程系,山東 煙臺264001)

0引言

波達方向(DOA)估計在雷達、無線通信、聲納等領域有著廣泛應用，經(jīng)典的子空間類算法(如MUSIC算法[1]、ESPRIT算法[2])能夠?qū)崿F(xiàn)DOA的超分辨估計，但要求接收數(shù)據(jù)的協(xié)方差矩陣的秩等于信源個數(shù)，且估計的準確性以采用大量快拍數(shù)據(jù)獲取高精度的信號子空間或噪聲子空間為前提。而在實際應用中，相干信號普遍存在，且往往只能獲得有限的快拍數(shù)據(jù)，尤其是在處理短時突發(fā)信號或運動目標信號時，只有單快拍可用，此時傳統(tǒng)的DOA估計算法往往失效。

近年來，采用單次快拍數(shù)據(jù)的相干信源DOA估計得到關注[3-12]，并取得了一定的研究成果。文獻[5]直接采用單次快拍數(shù)據(jù)構(gòu)造Toeplitz矩陣實現(xiàn)DOA估計(簡稱DTOP)，實時性好；文獻[6]在文獻[5]的基礎上，對單次快拍數(shù)據(jù)進行相關處理后再構(gòu)造Toeplitz矩陣，提高了低信噪比條件下的DOA估計性能；文獻[7]利用虛擬陣列接收數(shù)據(jù)構(gòu)造偽協(xié)方差矩陣，并利用酉變換實現(xiàn)DOA估計，降低了計算復雜度。但文獻[5-7]中的算法要求所有入射信號均具有零初始相位，這在實際應用中幾乎無法滿足[8]。文獻[10]在文獻[5]的基礎上，以陣列的中心陣元作為參考陣元，利用單次快拍數(shù)據(jù)構(gòu)造偽協(xié)方差矩陣實現(xiàn)DOA估計(簡稱ENDTOP)，克服了所有入射信號均需具有零初始相位的局限性，但該算法要處理M個信源，需要的接收陣元個數(shù)N≥2M-1，嚴重制約了其信號處理能力。文獻[13]提出一種利用單次快拍接收信號重排的偽協(xié)方差矩陣構(gòu)造模型，文獻[5-6,10]的偽協(xié)方差矩陣構(gòu)造方法可以統(tǒng)一在該模型下，該模型主要針對基于特征值分解的空間譜估計。

鑒于此，本文從實際應用出發(fā)，針對單次快拍情況下相干信源DOA估計的解相干和實時性問題，在經(jīng)典空間平滑算法[14]的基礎上，提出了基于空間平滑的單快拍DOA估計算法(SS-FBSS)。

1信號模型和空間平滑算法

1.1信號模型

假設一個標準均勻線陣，包含N個相同的全向陣元，陣元間距d≤λ/2(λ是入射波長)?，F(xiàn)有M個遠場窄帶信源(M

(1)

式(1)中，si(t)表示第i個窄帶遠場信號源，nk(t)表示第k陣元接收到的理想高斯白噪聲。則陣列接收的數(shù)據(jù)可以表示為：

X(t)=[x1(t),x2(t),…,xN(t)]T=AS(t)+N(t)

(2)

假設前r個信源相干，

si(t)=αis1(t)(i=1,2,…,r)

(3)

其中復常數(shù)αi表示si(t)相對于s1(t)的相干系數(shù)。

陣列的協(xié)方差矩陣表示為:

(4)

1.2空間平滑算法

空間平滑算法[15]是目前一種有效的解相干算法，它利用均勻線陣的平移不變性，將一個N元均勻線陣分成相互重疊的l個子陣，每個子陣的陣元數(shù)為m，即有N=l+m-1。采用前向空間平滑算法得到l個子陣，如圖1所示，對各個子陣的子協(xié)方差矩陣進行算術平均，得到前向空間平滑的協(xié)方差矩陣

(5)

前后向空間平滑(FBSS)的相關矩陣為：

(6)

式中J是反對角線元素為1的交換矩陣。

圖1　前向空間平滑算法原理Fig.1　The method of forward spatial smoothingalgorithm

空間平滑算法是通過求解各個子陣協(xié)方差矩陣的均值來消除信源之間的相干性的，為保證計算精度，往往需采用大量快拍數(shù)據(jù)。

2基于空間平滑的單快拍DOA估計算法

考慮單次快拍情況下的解相干問題，同時提高算法的實時性，本文提出基于空間平滑的單快拍DOA估計算法(SS-FBSS)。

2.1基于空間平滑的單快拍數(shù)據(jù)解相干預處理

以第1個子陣為參考子陣，第k個子陣的接收數(shù)據(jù)為:

(7)

直接采用t時刻的單次快拍數(shù)據(jù)，以l個子陣的接收數(shù)據(jù)構(gòu)造維數(shù)為m×l的矩陣

(8)

同理，采用后向空間平滑技術劃分l個子陣，第i個子陣的接收數(shù)據(jù)為

(9)

采用t時刻的單次快拍數(shù)據(jù)進行后向空間平滑，以l個子陣構(gòu)造維數(shù)為m×l的矩陣

(10)

式(10)中,J為反對角線元素為1的交換矩陣，下標為矩陣的維數(shù)。

利用式(8)、式(10)構(gòu)造偽協(xié)方差矩陣：

(11)

(12)

2.2DOA估計

考慮R為centro-Hermitian矩陣，利用酉變換可將其轉(zhuǎn)化成實矩陣。定義酉矩陣Q：

(13)

其中下標表示矩陣的維數(shù)，上述矩陣具有如下兩個性質(zhì)：

JQ*=Q

(14)

QQH=I

(15)

由式(13)定義的矩陣Q對偽協(xié)方差矩陣RG做酉變換

(16)

(17)

在此基礎上，采用MUSIC、ESPRIT算法可實現(xiàn)DOA估計。

1)MUSIC算法

變換后的導向向量為:

(18)

利用式(17)、(18)求解MUSIC算法的譜估計表達式為[1]:

(19)

2)ESPRIT算法

根據(jù)文獻[15]可知，對于復值的旋轉(zhuǎn)不變性有K2US=K1USΨ，則實值的旋轉(zhuǎn)不變性有

(20)

(21)

(22)

采用總體最小二乘法(TLS)求解式(20)可得Ψ′，對Ψ′進行特征值分解，有

(23)

2.3算法步驟總結(jié)

1)采用單次快拍數(shù)據(jù)根據(jù)式(8)、(10)、(11)構(gòu)造centro-Hermitian偽協(xié)方差矩陣RG；

2)利用式(13)定義的酉矩陣，將偽協(xié)方差矩陣RG變換到實值空間T，如式(16)；

4)采用MUSIC或ESPRIT算法求解各個信號的入射角。

2.4算法性能分析

1)計算復雜度與前后向空間平滑[15]算法相比，本文提出的SS-FBSS算法直接采用單次快拍數(shù)據(jù)的空間平滑子陣構(gòu)造偽協(xié)方差矩陣RG，而無需計算各個子陣的協(xié)方差矩陣，減小了運算量；并且該算法構(gòu)造的偽協(xié)方差矩陣RG是centro-Hermitian矩陣，進行酉變換處理后，算法中的奇異值分解由原來的復數(shù)域變?yōu)閷崝?shù)域，計算復雜度大大降低。

2)分辨能力在SS-FBSS算法中，構(gòu)造的偽協(xié)方差矩陣RG要實現(xiàn)“解相關”需要滿足m≥M且2l≥M，而陣元數(shù)N=l+m-1。由此可知，當有M個信號源、每個子陣的陣元數(shù)為M+1時，要分辨這些信源，SS-FBSS算法需要陣元個數(shù)N≥1.5M，其分辨能力與FBSS算法一致，比ENDTOP算法高。但該算法也是通過降低自由度來實現(xiàn)解相干的。

3)估計精度SS-FBSS算法利用前后向空間平滑子陣構(gòu)造偽協(xié)方差矩陣RG，克服了所有入射信號均具有零初始相位的局限性，且充分利用了單次快拍數(shù)據(jù)，估計精度高于DTOP算法、ENDTOP算法。

3數(shù)值仿真

仿真1假設一個均勻線陣，陣元數(shù)為5，陣元間距為λ/2。兩個相干實值信號，θ1=-20°，θ2=5°，信噪比SNR為10dB，分別采用SS-FBSS算法和DTOP[6]算法進行解相干處理，而后采用MUSIC算法實現(xiàn)DOA估計。圖2(a)、(b)分別給出了在α1=α2=1和α1=1、α2=exp(jπ/3)條件下，兩種算法處理后的MUSIC譜圖。由圖2可以看出，當在α1=α2=1時，SS-FBSS算法與DTOP算法皆能實現(xiàn)解相干，而在α1=1、α2=exp(jπ/3)情況下，DTOP算法已經(jīng)失效，而本文提出的SS-FBSS算法能準確地估計出波達角。這說明SS-FBSS算法克服了所有入射信號均具有零初始相位的局限性，無需對入射信號的初始相位作特別假定。

圖2　SS-FBSS算法與DTOP算法的MUSIC譜圖Fig.2　MUSIC spectrum of SS-FBSS and DTOP algorithms

仿真23個完全相干信源入射到均勻線陣，陣元數(shù)為15，陣元間隔為λ/2，入射角分別為θ1=-20°、θ2=5°、θ3=43°，信噪比SNR從0dB到20dB以0.5dB的間隔變化，分別采用SS-FBSS算法、ENDTOP[11]算法和FBSS算法[15]進行解相干處理，而后采用ESPRIT算法實現(xiàn)DOA估計。每個SNR做200次蒙特卡羅實驗，其中FBSS算法采用150次快拍數(shù)據(jù)，仿真結(jié)果如圖3所示，圖3(a)為均方根誤差隨SNR變化曲線，圖3(b)為成功概率隨SNR變化曲線(成功概率定義為均方根誤差小于0.2°的概率)。從圖3中可以看出，本文提出的SS-FBSS算法的均方根誤差、成功概率與ENDTOP算法基本一致，略優(yōu)于ENDTOP算法。在低信噪比條件下，該算法的DOA估計誤差較大，特別是信噪比在0～6dB時，隨著信噪比的增大，均方根誤差迅速減小，成功概率由不足10%迅速增大至65%。 而在信噪比大于6dB時，該算法的估計性能接近于采用150次快拍數(shù)據(jù)的FBSS算法。

圖3　對完全相干信源的估計性能隨SNR的變化曲線Fig.3　Estimation performance with SNR for coherent signals

仿真33個信源入射到一個均勻線陣，前兩個為相干信源，α1=1、α2=exp(jπ/6)，第3個信源與第1個不相關，其他條件與仿真2相同，仿真結(jié)果如圖4所示。從圖4可以看出，該算法對部分相干信源的估計性能與對完全相干信源的估計性能基本一致，略優(yōu)于ENDTOP算法，在低信噪比情況下估計誤差較大，且隨著信噪比的增大估計誤差迅速減小。但在信噪比大于6dB時，該算法的估計性能接近于采用150次快拍數(shù)據(jù)的FBSS算法。通過比較圖3與圖4還可以發(fā)現(xiàn)，與完全相干信源相比，F(xiàn)BSS算法對部分相干且相干系數(shù)為復數(shù)的信源DOA估計性能明顯降低，在低信噪比情況下，差異尤為明顯。

圖4　對部分相干信源的估計性能隨SNR的變化曲線Fig.4　Estimation performance with SNR for partial coherent signals

4結(jié)論

本文提出了基于空間平滑的單快拍DOA估計算法。該算法直接采用單次快拍數(shù)據(jù)的前后向空間平滑子陣，構(gòu)造centro-Hermitian偽協(xié)方差矩陣，實現(xiàn)解相干，而無需對入射信號的初始相位作特別假定，也無需計算各個子陣的協(xié)方差矩陣，減小了計算量；利用酉變換將復數(shù)域的奇異值分解實數(shù)化，降低了計算復雜度；最后結(jié)合MUSIC、ESPRIT等子空間類算法，實現(xiàn)了相干信源的DOA估計。數(shù)值仿真表明：該算法無需對入射信號的初始相位作特別假定；對于完全相干信源或部分相干信源，估計性能優(yōu)于ENDTOP算法，且在信噪比大于6dB的情況下，接近于采用150次快拍數(shù)據(jù)的FBSS算法，能夠有效實現(xiàn)單快拍DOA的快速估計，特別適用于對實時性要求較高的DOA估計場合。但該算法是通過降低自由度來實現(xiàn)解相干的，且在低信噪比情況下，估計誤差較大，針對這些問題還需進一步研究。

參考文獻:

[1]SchmidtRO.Multipleemitterlocationandsignalparameterestimation[J].IEEETransactionsonAntennasandPropagation, 1986,34(3): 276-280.

[2]RoyRKailathT.ESPRIT-estimationofsignalparametersviarotationalinvariancetechniques[J].IEEETransactionsonAcoustics,SpeechandSignalProcessing, 1989, 37(7): 984-995.

[3]AlcalaBE,RecioFR.Performanceanalysisofmetaheuristicapproachesforsingle-snapshotDOAestimation[J].IEEEAntennasandWirelessPropagationLetters, 2013, 12(1): 166-169.

[4]梁浩, 崔琛, 伍波, 等. 雙基地MIMO雷達相干目標角度估計算法[J]. 系統(tǒng)工程與電子技術, 2014, 36(6): 1068-1074.

[5]謝鑫, 李國林, 劉華文. 采用單次快拍數(shù)據(jù)實現(xiàn)相干信號DOA估計[J]. 電子與信息學報, 2010, 32(3): 604-608.

[6]梁浩, 李小波, 王磊. 采用單次快拍數(shù)據(jù)實現(xiàn)信源DOA估計[J]. 數(shù)據(jù)采集與處理, 2013, 28(1): 58-63.

[7]閆金山, 彭秀艷, 王咸鵬. 基于酉變換的虛擬陣列DOA估計算法[J]. 哈爾濱工業(yè)大學學報, 2012, 44(4): 136-140.

[8]DelmasJP.Commentson‘ConjugateESPRIT(C-SPRIT)’ [J].IEEETransactionsonAntennasandPropagation, 2007, 55(2): 511.

[9]Jung-TaeKim,Sung-HoonMoon,DongSeogHan,etal.FastDOAestimationalgorithmusingpseudocovariancematrix[J].IEEETransactionsonAntennasandPropagation, 2005, 53(4): 1346-1351.

[10]王凌, 李國林, 謝鑫. 互耦效應下用單快拍數(shù)據(jù)實現(xiàn)相干信源完全解相干和解耦合[J]. 電子與信息學報, 2012, 34(10): 2532-2536.

[11]王凌, 李國林, 孟晶, 等. 用一次快拍數(shù)據(jù)實現(xiàn)二維完全解相干和解互耦[J]. 系統(tǒng)工程與電子技術, 2012, 34(11): 2208-2213.

[12]王凌, 李國林, 李靜, 等. 一種新的單次快拍二維ESPRIT算法[J]. 北京理工大學學報, 2013, 33(1): 99-104.

[13]蔣柏峰, 呂曉德, 向茂生. 一種基于陣列接收信號重排的單快拍DOA估計方法[J]. 電子與信息學報, 2014, 36(6): 1334-1339.

[14]PillaiSU,KwonBH.Forward/backwardspatialsmoothingtechniquesforcoherentsignalidentification[J].IEEETransactionsonAcoustics,SpeechandSignalProcessing, 1989, 37(1): 8-15.

[15]王永良, 陳輝, 彭英寧, 等. 空間譜估計理論與算法[M]. 北京: 清華大學出版社, 2004: 199.

摘要:針對單次快拍情況下相干信源波達方向(DOA)估計的解相干和實時性問題，提出了基于空間平滑的單快拍DOA估計算法。該算法直接采用單次快拍數(shù)據(jù)的前后向空間平滑子陣，構(gòu)造centro-Hermitian偽協(xié)方差矩陣，實現(xiàn)解相干，而無需對入射信號的初始相位作特別假定，也無需計算各個子陣的協(xié)方差矩陣，減小了計算量；利用酉變換將復數(shù)域的奇異值分解實數(shù)化，降低了計算復雜度；結(jié)合MUSIC、ESPRIT等子空間類算法，實現(xiàn)了相干信源的DOA估計。數(shù)值仿真表明該算法無需對入射信號的初始相位作特別假定，且對于相干信源能夠有效實現(xiàn)單快拍DOA的快速估計，適用于對實時性要求較高的場合。

關鍵詞:波達方向估計；相干信源；空間平滑；單次快拍；酉變換

Single Snapshot DOA Estimation Algorithm Based on Spatial SmoothingLIU Xiaodi, ZHOU Xinli, XIAO Jinguang

(Department of Electronic Information Engineering, Naval Aeronautics and Astronautics University, Yantai 264001, China)

Abstract:A direction-of-arrival (DOA) estimation algorithm based on spatial smoothing was proposed for the decorrelation and instant problems among DOA estimation of coherent signals when there was only one valid snapshot. The forward/backward spatial smoothing subarrays of single snapshot were directly used to construct the pseudocovariance matrix to accomplish decorrelation, which does not require that signals’ natural phases must be zero, and that the covariance of every subarray must be computed. The pseudocovariance, which was a centro-Hermitian matrix, was changed into real one by unitary transformation to reduce the computational complexity. Combined with the subspace kind algorithms such as MUSIC and ESPRIT, the DOA ccould be estimated. Simulation results showed that the algorithm did not require special assumption on signal’s natural phases, and could estimate the single snapshot DOAs of coherent signals effectively, which was suitable for real-time demand higher occasion.

Key words:direction-of-arrival estimation; coherent signals; spatial smoothing; single snapshot; unitary transformation

中圖分類號:TN911.23

文獻標志碼:A

文章編號：1008-1194(2015)06-0066-05

作者簡介:劉曉娣(1982—)，女，山東昌邑人，工程師，博士研究生，研究方向：無線電波傳播、信號處理。E-mail：xdl_1982@163.com。

基金項目:國家自然科學基金項目資助(61179016)

*收稿日期:2015-04-12