基于二次虛擬擴(kuò)展的高分辨率波達(dá)方向估計(jì)方法

趙 宇，李文興，毛曉軍

(哈爾濱工程大學(xué) 信息與通信工程學(xué)院， 哈爾濱 150001)

基于二次虛擬擴(kuò)展的高分辨率波達(dá)方向估計(jì)方法

趙 宇，李文興，毛曉軍

(哈爾濱工程大學(xué) 信息與通信工程學(xué)院， 哈爾濱 150001)

為進(jìn)一步提高天線陣波達(dá)方向估計(jì)的分辨率，在四階量多重信號(hào)分類方法的基礎(chǔ)上，提出一種高分辨率的波達(dá)方向估計(jì)方法. 利用陣列接收數(shù)據(jù)的四階矩量進(jìn)行虛擬陣列擴(kuò)展，再利用陣列接收數(shù)據(jù)的共軛進(jìn)行虛擬陣列擴(kuò)展，實(shí)現(xiàn)二次虛擬擴(kuò)展；將擴(kuò)展后的陣列導(dǎo)向矢量和協(xié)方差矩陣用于波達(dá)方向估計(jì)，與原陣列導(dǎo)向矢量和協(xié)方差矩陣相比，相當(dāng)于構(gòu)造了更多的虛擬陣元，并擴(kuò)展了陣列的孔徑. 仿真結(jié)果表明：與四階量多重信號(hào)分類等波達(dá)方向估計(jì)方法相比，所提出的方法在波達(dá)方向估計(jì)中成功概率更高，均方誤差更低，具有更高的分辨率. 所提方法通過二次虛擬擴(kuò)展，構(gòu)造了更多的虛擬陣元，有效地提高了天線陣波達(dá)方向估計(jì)的分辨率.

波達(dá)方向估計(jì)；多重信號(hào)分類；四階量法；虛擬陣列；分辨率

波達(dá)方向(direction-of-arrival, DOA)估計(jì)在雷達(dá)、聲吶、無線通信等眾多應(yīng)用領(lǐng)域有著重要作用[1-2]. Schmidt等[3]提出的多重信號(hào)分類(multiple signal classification, MUSIC)是最具有代表性的DOA估計(jì)方法，是一種基于信源協(xié)方差矩陣特征分解的子空間類方法，具有較高的分辨率.

隨著科技的快速發(fā)展，對(duì)DOA估計(jì)的分辨率也提出更高的要求. 虛擬天線陣是一種研究構(gòu)造虛擬陣元的方法和將實(shí)際天線陣轉(zhuǎn)換為虛擬天線陣的先進(jìn)技術(shù)，主要包括內(nèi)插變換法[4-5]和四階量法[6-8]. 將虛擬天線陣技術(shù)與MUSIC方法結(jié)合，可以增加天線陣的自由度，有效提高天線陣DOA估計(jì)的分辨率. 內(nèi)插變換法通過在實(shí)際陣元間插入虛擬陣元，提高了MUSIC方法的分辨率[5]，但是內(nèi)插變換法不利于陣列孔徑擴(kuò)展，對(duì)分辨率的提高有限. 基于四階量的MUSIC(fourth-order MUSIC, FOMUSIC)方法能夠有效地?cái)U(kuò)展陣列孔徑，并提高陣列的分辨率[6]，是一種高分辨率的DOA估計(jì)方法.

FOMUSIC方法能夠擴(kuò)展陣列孔徑，但是運(yùn)算量較大[9]. Akkar等[10]利用正交分解代替特征分解； Lie等[11]用二階和四階協(xié)方差矩陣進(jìn)行混合估計(jì)；Wang等[12]則去除了四階量矩陣中的部分冗余信息，都減小了FOMUSIC方法的運(yùn)算量. Liao等[13]利用基于行列式的空間譜，實(shí)現(xiàn)了未知互耦情況下的DOA估計(jì)；宋海巖等[14]用一組四階量矩陣進(jìn)行聯(lián)合對(duì)角化，能處理相干信號(hào)，且與單個(gè)四階量矩陣相比獲得了更高的分辨率；Li[15]等通過增強(qiáng)信號(hào)子空間，減小了內(nèi)插變換法的變換誤差，與常規(guī)內(nèi)插變換法相比分辨率更高，且運(yùn)算量較??；Shan等[16]提出了共軛擴(kuò)展多重信號(hào)分類(conjugate augmented MUSIC, CAM)方法，通過擴(kuò)展陣列導(dǎo)向矢量以及協(xié)方差矩陣所含信息，大幅提高了分辨率，是一種分辨率極高的虛擬天線陣DOA估計(jì)方法.

為了進(jìn)一步提高天線陣DOA估計(jì)的分辨率，本文提出一種基于二次虛擬擴(kuò)展的多重信號(hào)分類(twice virtual expansion MUSIC, TVEM)方法. 該方法先求出陣列接收數(shù)據(jù)的四階矩量，以及相應(yīng)擴(kuò)展的陣列導(dǎo)向矢量，再求出陣列接收數(shù)據(jù)及其共軛，以及相應(yīng)擴(kuò)展的陣列導(dǎo)向矢量，實(shí)現(xiàn)二次虛擬擴(kuò)展，構(gòu)造了更多的虛擬陣元，進(jìn)一步提高了天線陣DOA估計(jì)的分辨率.

1 天線陣模型及MUSIC方法

(1)

式中：M為信源數(shù)，(·)T表示轉(zhuǎn)置，A=[a1(θ),a2(θ), …,aM(θ)]為陣列導(dǎo)向矢量矩陣，S(t)為信號(hào)復(fù)包絡(luò)，N(t)為天線陣的噪聲.

陣列接收數(shù)據(jù)的協(xié)方差矩陣為

協(xié)方差矩陣Rxx是Hermitian矩陣，特征分解后為

Rxx=UΣUH,

式中Σ=diag[λ1,…,λN]是Rxx的特征值按降序排列組成的對(duì)角陣，U為特征向量矩陣.

將Rxx的特征向量根據(jù)特征值的大小分為信號(hào)和噪聲兩個(gè)部分：

式中US為大特征值所對(duì)應(yīng)的信號(hào)子空間，UN為小特征值所對(duì)應(yīng)的噪聲子空間.

根據(jù)協(xié)方差矩陣信號(hào)子空間與噪聲子空間的正交性，可以由噪聲子空間估計(jì)出信號(hào)的波達(dá)方向，在實(shí)際計(jì)算中，Rxx通常用有限次快拍數(shù)K的采樣協(xié)方差矩陣代替：

MUSIC算法的空間譜為

式中a(θ)為陣列導(dǎo)向矢量. 再通過對(duì)PMUSIC譜峰搜索即可估計(jì)出信號(hào)的波達(dá)方向.MUSIC方法是一種分辨率較高的DOA估計(jì)方法，但是當(dāng)信源數(shù)超過實(shí)際陣元數(shù)時(shí)，MUSIC方法將會(huì)失效，此外，MUSIC方法的分辨率也受到實(shí)際陣列孔徑的限制.

2 FOMUSIC方法

用四階量替代二階量進(jìn)行DOA估計(jì)，可以保留更豐富的波達(dá)方向信息，獲得更好的DOA估計(jì)性能[7]. 設(shè)x(t)為零均值平穩(wěn)實(shí)隨機(jī)過程，則其四階累積量矩陣為

(2)

式中(·)*表示共軛. 式(2)可以表示為[7]C4x=E{[X(t)?X*(t)][X(t)?X*(t)]H}- E{X(t)?X*(t)}E{[X(t)?X*(t)]H}-E{X(t)XH(t)}?E{[X(t)XH(t)]*}.

(3)

式中?表示克羅內(nèi)克積，第一項(xiàng)為四階矩量協(xié)方差矩陣：

D4x=E{[X(t)?X*(t)][X(t)?X*(t)]H}.

式中：X(t)為N×K維矩陣，X(t)?X*(t)為N2×K維矩陣，D4x為N2×N2維矩陣. 由式(3)可以得到陣列接收數(shù)據(jù)的四階協(xié)方差矩陣.

根據(jù)克羅內(nèi)克積的性質(zhì)[17]，將式(1)代入X(t)?X(t)*：

X(t)?X*(t)=[AS(t)+N(t)]?

[AS(t)+N(t)]*=[AS(t)]?

[AS(t)]*+N(t)?N*(t)=

(A?A*)[S(t)?S*(t)]+

N(t)?N*(t).

(4)

再將式(4)代入D4x[7]：

D4x=E{[X(t)?X*(t)][X(t)?X*(t)]H}=

(A?A*)E{[S(t)?S*(t)]·

[S(t)?S*(t)]H}(A?A*)H+

E{[N(t)?N*(t)][N(t)?N*(t)]H}=

(A?A*)D4S(A?A*)H+D4N.

(5)

式中：A?A*為擴(kuò)展后的陣列導(dǎo)向矢量矩陣，D4S和D4N分別為信號(hào)和噪聲的四階矩量矩陣.

與推導(dǎo)四階矩量矩陣的式(5)類似，可由式(3)推導(dǎo)出四階累積量矩陣為

C4x=(A?A*)C4S(A?A*)H+C4N,

(6)

式中C4S和C4N分別為信號(hào)和噪聲的四階累積量矩陣. 由式(5)和式(6)可知，與二階協(xié)方差矩陣相比，四階協(xié)方差矩陣的維度得到了擴(kuò)展，所含波達(dá)方向信息也更豐富，實(shí)際陣列通過克羅內(nèi)克積運(yùn)算被轉(zhuǎn)換為虛擬陣列，陣列孔徑得到了擴(kuò)展.

由式(3)得到C4x，再對(duì)其特征分解得到噪聲子空間U4N，根據(jù)噪聲子空間與信號(hào)子空間的正交性，可知，F(xiàn)OMUSIC方法的空間譜為

FOMUSIC方法將實(shí)際陣列轉(zhuǎn)換為虛擬陣列，擴(kuò)展了陣列孔徑，能夠在信源數(shù)超過實(shí)際陣元數(shù)的情況下進(jìn)行DOA估計(jì)，與MUSIC方法相比，分辨率也更高.

3 TVEM方法

為了進(jìn)一步提高天線陣DOA估計(jì)的分辨率，本文在FOMUSIC方法的基礎(chǔ)上，提出了TVEM方法，該方法通過二次虛擬擴(kuò)展，構(gòu)造了更多的虛擬陣元，能夠獲得比FOMUSIC方法更高的分辨率.

先求出陣列接收數(shù)據(jù)的四階矩量，進(jìn)行第一次虛擬擴(kuò)展，陣列接收數(shù)據(jù)X(t)的四階矩量為

Xu1(t)=X(t)?X*(t).

(7)

根據(jù)克羅內(nèi)克積的性質(zhì)[17]，將式(1)代入式(7)得

Xu1(t)=X(t)?X*(t)=[AS(t)+N(t)]?

[AS(t)+N(t)]*=[AS(t)]?[AS(t)]*+

N(t)?N*(t)=(A?A*)[S(t)?S*(t)]+

加入死亡話題的親子談話這個(gè)中介變量后，采用線性回歸分析來考察父母來生信念和死亡話題的親子談話與兒童來生信念的關(guān)系。依次將父母的來生信念(BA量表得分)、死亡話題親子談話(PDCA量表得分)納入回歸方程，進(jìn)行回歸。采用Sobel檢驗(yàn)方法對(duì)心理理論的中介效應(yīng)進(jìn)行檢驗(yàn)，Z=-2.36，p<0.05。當(dāng)將親子談話納入回歸模型時(shí)，父母來生信念的預(yù)測(cè)作用不再顯著，表明親子談話起到顯著、完全的中介作用(具體指標(biāo)見表7)。

N(t)?N*(t),

式中Xu1為N2×K維矩陣，A?A*為擴(kuò)展后的陣列導(dǎo)向矢量矩陣. 擴(kuò)展后的陣列導(dǎo)向矢量為b(θ)=a(θ)?a*(θ)=[1,e,…,e]T?[1,e,…,e]H=

[1,e,…， e,e,…,e]T.

(8)

由式(8)可知，b(θ)中有大量冗余，只有2N-1個(gè)非重復(fù)元素，即第一次虛擬陣列擴(kuò)展的陣元位置，如圖1所示.

圖1 第一次虛擬擴(kuò)展的虛擬陣元位置

Fig.1 Position diagram of virtual array elements of the first extension

Fig.2 Position diagram of virtual array elements of the second extension

求出二次虛擬擴(kuò)展后的擴(kuò)展協(xié)方差矩陣：

(9)

對(duì)Ru進(jìn)行特征分解：

式中UuN為噪聲子空間.

二次虛擬擴(kuò)展后的虛擬陣列導(dǎo)向矢量為

(10)

由式(10)和圖1、圖2可知，通過二次虛擬擴(kuò)展，實(shí)際陣列被轉(zhuǎn)換為虛擬陣列，陣列的導(dǎo)向矢量得到了進(jìn)一步的擴(kuò)展，構(gòu)造了更多的虛擬陣元.TVEM方法的空間譜為

TVEM方法先求出陣列接收數(shù)據(jù)的四階矩量，以及相應(yīng)擴(kuò)展的陣列導(dǎo)向矢量，再求出陣列接收數(shù)據(jù)及其共軛，以及相應(yīng)擴(kuò)展的陣列導(dǎo)向矢量，通過二次虛擬擴(kuò)展，將實(shí)際陣列轉(zhuǎn)換為虛擬陣列.TVEM方法分辨率提高的原因可以從數(shù)學(xué)和物理意義兩方面進(jìn)行分析: 從數(shù)學(xué)上看，如式(6)、(9)所示，TVEM方法中協(xié)方差矩陣得到擴(kuò)展，其中包含了更多的波達(dá)方向信息和冗余，從而提高了分辨率；從物理意義上看，如圖1、2所示，TVEM方法構(gòu)造了更多的虛擬陣元，擴(kuò)展了陣列孔徑，因此獲得了比FOMUSIC方法更高的分辨率.

4 運(yùn)算量分析

將常規(guī)MUSIC方法、FOMUSIC方法、CAM方法[16]以及本文提出的TVEM方法的計(jì)算復(fù)雜度進(jìn)行比較，如表1所示.

表1 4種方法計(jì)算復(fù)雜度對(duì)比

由表1可知，在未去冗余的情況下，F(xiàn)OMUSIC法、CAM法和TVEM法這3種基于四階量的虛擬天線陣DOA估計(jì)方法的計(jì)算復(fù)雜度遠(yuǎn)高于MUSIC法，這是因?yàn)樗碾A量矩陣中包含大量冗余，其中CAM法計(jì)算復(fù)雜度比TVEM法高，TVEM法計(jì)算復(fù)雜度比FOMUSIC法高. 在天線陣為均勻直線陣的情況下，能夠去除四階量矩陣中的冗余[16]，此時(shí)TVEM法的計(jì)算復(fù)雜度與CAM方法十分接近，而TVEM法的計(jì)算復(fù)雜度明顯高于FOMUSIC法，這是因?yàn)門VEM法通過二次虛擬擴(kuò)展，構(gòu)造了比FOMUSIC法更多的虛擬陣元，從而獲得了更高的分辨率，其代價(jià)是運(yùn)算量增大.

5 仿真分析

考慮遠(yuǎn)場(chǎng)窄帶情況下的來波信號(hào)，信號(hào)之間互不相關(guān)，天線陣為四元均勻直線陣，陣元間距為半波長，快拍數(shù)為200. 下面通過三組仿真實(shí)驗(yàn)，比較TVEM方法、FOMUSIC方法和CAM方法[16]的分辨率，三組仿真結(jié)果均通過200次重復(fù)獨(dú)立實(shí)驗(yàn)得到. 分辨率的評(píng)價(jià)指標(biāo)采用均方誤差和成功概率. 成功概率的計(jì)算方法為：獨(dú)立重復(fù)實(shí)驗(yàn)中，波達(dá)方向估計(jì)的誤差低于門限值的次數(shù)與總次數(shù)之比，仿真中門限值設(shè)為0.2°.

5.1 在信源數(shù)少于實(shí)際陣元數(shù)的情況下，比較3種方法的分辨率

仿真1 信源數(shù)為3，來波方向?yàn)?45°、0°和40°，其中天線陣法線方向設(shè)為0°，左側(cè)為-90°到0，右側(cè)為0°到90°.

圖3為3種方法的均方誤差(root mean square error, RMSE)隨輸入信噪比RSN的變化曲線，其中RSN表示天線陣輸入的期望信號(hào)功率與噪聲功率之比；圖4為3種方法的成功概率隨輸入信噪比的變化曲線.

圖3 輸入信噪比變化時(shí)3種方法的均方誤差曲線

圖4 輸入信噪比變化時(shí)3種方法的成功概率曲線

Fig.4 Probability of target resolution versusRSNofthethreemethods

從圖3可以看出:當(dāng)RSN低于3dB時(shí)，TVEM法的均方誤差高于CAM法，當(dāng)RSN高于3dB時(shí)，TVEM法的均方誤差低于CAM法；總體上，TVEM法的均方誤差略低于CAM方法，且二者的均方誤差都明顯低于FOMUSIC法. 從圖4可以看出:當(dāng)RSN低于3dB時(shí)，TVEM法的成功概率低于CAM法，當(dāng)RSN高于3dB時(shí)，TVEM法的成功概率高于CAM法；總體上，TVEM法的成功概率略高于CAM法，且二者的成功概率都明顯高于FOMUSIC法. 由此可知，在信源數(shù)少于實(shí)際陣元數(shù)的情況下，TVEM法的分辨率略高于CAM法，明顯高于FOMUSIC法.

5.2 在波達(dá)角度較為接近的情況下，比較3種方法的分辨率

仿真2 信源數(shù)為2，來波方向?yàn)?3°和3°.

圖5為3種方法的均方誤差隨輸入信噪比的變化曲線. 圖6為3種方法的成功概率隨輸入信噪比的變化曲線.

從圖5可以看出,在波達(dá)角度較為接近的情況下，隨著RSN不斷提高，3種DOA估計(jì)方法的均方誤差都不斷下降，3種方法均方誤差曲線的變化過程都可以大致分為3個(gè)階段. 第一階段，RSN相對(duì)較低時(shí)，均方誤差較大且下降較為緩慢；第二階段，隨著RSN的提高，均方誤差快速下降；第三階段，RSN提高到一定值后，均方誤差相對(duì)較小且平穩(wěn)下降. 這是因?yàn)楫?dāng)RSN較低時(shí)，DOA估計(jì)方法無法分辨波達(dá)角度較為接近的信源，因而均方誤差很大且下降較為緩慢；隨著RSN不斷提高，DOA估計(jì)方法可以在一定程度上分辨波達(dá)角度較為接近的信源，此時(shí)均方誤差快速下降；當(dāng)RSN達(dá)到一定值后，DOA估計(jì)方法可以相對(duì)有效地分辨波達(dá)角度較為接近的信源，則均方誤差相對(duì)較小且平穩(wěn)下降. 對(duì)比圖5中3種方法的均方誤差曲線可以知道，TVEM法的均方誤差曲線在RSN超過9dB后進(jìn)入第三階段，CAM法的均方誤差曲線在RSN超過11dB后進(jìn)入第三階段，F(xiàn)OMUSIC法的均方誤差曲線在RSN超過17dB后進(jìn)入第三階段. 此外，在第二、三階段，TVEM法的均方誤差略低于CAM法，明顯低于FOMUSIC法. 由此可知，TVEM法的均方誤差略低于CAM法，明顯低于FOMUSIC法.

圖6 輸入信噪比變化時(shí)3種方法的成功概率曲線

Fig.6 Probability of target resolution versusRSNofthethreemethods

從圖6可以看出:當(dāng)RSN低于5dB時(shí)，TVEM法的成功概率略低于CAM法，當(dāng)RSN高于5dB時(shí)，TVEM法的成功概率高于CAM法；總體上，TVEM法的成功概率略高于CAM法，且二者的成功概率都明顯高于FOMUSIC法. 由圖5和圖6可知，在波達(dá)角度較為接近的情況下，TVEM法的分辨率略高于CAM法， 明顯高于FO-MUSIC法.

5.3 在信源數(shù)超過實(shí)際陣元數(shù)的情況下，比較3種方法的分辨率

仿真3 信源數(shù)為5，來波方向?yàn)?50°、-25°、0°、15°和45°.

圖7為3種方法的均方誤差隨輸入信噪比的變化曲線，圖8為3種方法的成功概率隨輸入信噪比的變化曲線.

圖7 輸入信噪比變化時(shí)3種方法的均方誤差曲線

圖8 輸入信噪比變化時(shí)3種方法的成功概率曲線

Fig.8 Probability of target resolution versusRSNofthethreemethods

從圖7可以看出，TVEM法的均方誤差低于CAM法，CAM法的均方誤差低于FOMUSIC法. 從圖8可以看出，TVEM法的成功概率高于CAM法，CAM法的成功概率高于FOMUSIC法. 由此可知，在信源數(shù)超過實(shí)際陣元數(shù)的情況下，TVEM法的分辨率略高于CAM法，明顯高于FOMUSIC法.

由上述仿真結(jié)果可知，TVEM方法的分辨率略高于CAM方法，且TVEM方法的分辨率明顯高于FOMUSIC方法. 因此，TVEM方法是進(jìn)一步提高天線陣DOA估計(jì)分辨率的一種有效方法.

6 結(jié) 論

本文提出的基于二次虛擬擴(kuò)展的高分辨率波達(dá)方向估計(jì)方法，通過陣列接收數(shù)據(jù)的四階矩量和共軛，擴(kuò)展了陣列的導(dǎo)向矢量和協(xié)方差矩陣，構(gòu)造出更多的虛擬陣元，并擴(kuò)展了陣列的孔徑. 將實(shí)際天線陣轉(zhuǎn)換為虛擬天線陣，實(shí)現(xiàn)了二次虛擬擴(kuò)展，進(jìn)一步提高了天線陣DOA估計(jì)的分辨率. 仿真結(jié)果表明，本文所提出的方法與四階量多重信號(hào)分類等DOA估計(jì)方法相比，具有更低的均方誤差和更高的成功概率，是一種提高天線陣DOA估計(jì)分辨率的有效方法.

[1] WANG Yi, YANG Minglei, CHEN Baixiao, et al. Improved DOA estimation based on real-valued array covariance using sparse Bayesian learning [J]. Signal Processing, 2016, 129(C): 183-189. DOI: 10.1016/j.sigpro.2016.06.002.

[2] 閆鋒剛, 張薇, 金銘. 求根MUSIC初值設(shè)置和更新算法 [J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2015, 47(3): 88-92. DOI: 10.11918/j.issn.0367-6234.2015.03.015.

YAN Fenggang, ZHANG Wei, JIN Ming. A new method for setting and updating theinitiation of root-MUSIC [J]. Journal of Harbin Institute of Technology, 2015, 47(3): 88-92. DOI: 10.11918/j.issn.0367-6234.2015.03.015.

[3] SCHMIDT R O. Multiple emitter location and signal parameter estimation [J]. IEEE Transactions on Antennas and Propagation, 1986, 34(3): 276-280. DOI: 10.1109/TAP.1986.1143830.

[4] BELLONI F, RICHTER A, KOIVUNEN V. DOA estimation via manifold separation for arbitrary array structures [J]. IEEE Transactions on Signal Processing, 2007, 55(10): 4800-4810. DOI: dx.doi.org/10.1109/TSP.2007.896115

[5] WANG Yongliang, CHEN Hui, WAN Shanhu. An effective DOA method via virtual array transformation [J]. Science in China Series E, 2001, 44(1): 75-82. DOI: 10.1007/BF02916727.

[6] LI Xiukun, LI Tingting, GU Xinyu, et al. Array gain of fourth-order cumulants beamforming under typical probability denisity background [J]. Chinese Journal of Acoustics, 2015, 34(1): 15-26. DOI: 10.15949/j.cnki.0217-9776.2015.01.002.

[7] 丁齊, 魏平, 肖先賜. 基于四階累積量的DOA估計(jì)方法及其分析 [J]. 電子學(xué)報(bào), 1999, 27(3): 25-28. DOI: 10.3321/j.issn:0372-2112.1999.03.008.

DING Qi, WEI Ping, XIAO Xianci. Estimation and analysis of DOA based on fourth-order cumulant [J]. ACTA Electronica SINICA, 1999, 27(3): 25-28. DOI: 10.3321/j.issn:0372-2112.1999.03.008.

[8] ZHANG Xin, LIU Xiaoming, YU Haixia, et al. Improved MUSIC algorithm for DOA estimation of coherent signals via toeplitz and fourth-order-cumulants [J]. International Journal of Control and Automation, 2015, 8(10): 261-272. DOI: 10.14257/ijca.2015.8.10.25.

[9] 陳建, 王樹勛. 基于高階累積量虛擬陣列擴(kuò)展的DOA估計(jì) [J]. 電子與信息學(xué)報(bào), 2007, 29(5): 1041-1044.DOI: 10.3969/j.issn.1671-5896.2006.04.001

CHEN Jian, WANG Shuxun. DOA estimation of virtual array extension based on fourth-order cumulant [J]. Journal of Electronics and Information Technology, 2007, 29(5): 1041-1044. DOI: 10.3969/j.issn.1671-5896.2006.04.001.

[10]AKKAR S, HARABI F, GHARSALLAH A. MUSIC like algorithms for fast direction of arrival estimation [C]// 1stInternational Conference on Advanced Technologies for Signal and Image Processing. Sousse: IEEE, 2014: 550-554. DOI: 10.1109/ATSIP.2014.6834675.

[11]LIE S, LEYMAN A R, CHEW Y H. Fourth-order and weighted mixed order direction-of-arrival estimators [J]. IEEE Signal Processing Letters, 2006, 13(11): 691-694. DOI: 10.1109/LSP.2006.879456.

[12]WANG Qing, CHEN Hua, ZHAO Guohuang,et al. An improved direction finding algorithm based on toeplitz approximation [J]. Sensors, 2013, 13(1): 746-757. DOI: 10.3390/s130100746.

[13]LIAO Bin, CHAN S C. A cumulant-based method for direction finding in uniform linear arrays with mutual coupling [J]. IEEE Antennas and Wireless Propagation Letters, 2014, 13: 1717-1720. DOI: 10.1109/LAWP.2014.2352939.

[14] 宋海巖, 樸勝春. 基于高階累積量矩陣組正交聯(lián)合對(duì)角化的高分辨率方位估計(jì)方法 [J]. 電子與信息學(xué)報(bào), 2010, 32(4): 967-972. DOI: 10.3724/SP.J.1146.2008.01176.

SONG Haiyan, PIAO Shengchun. DOA estimation method based on orthogonal joint diagonalization of high-order cumulant [J]. Journal of Electronics and Information Technology, 2010, 32(4): 967-972. DOI: 10.3724/SP.J.1146.2008.01176.

[15]LI Wenxing, MAO Xiaojun, YU Wenhua, et al. An effective technique for enhancing direction finding performance of virtual arrays [J]. International Journal of Antenna and Propagation, 2014(3):1-7. DOI: 10.1155/2014/728463.

[16]SHAN Zhilong, YUM T S P. A conjugate augmented approach to direction-of-arrival estimation [J]. IEEE Transactions on Signal Processing, 2005, 53(11): 4104-4109. DOI: 10.1109/TSP.2005.857012.

[17]BREWER J W. Kronecker products and matrix calculus in system theory [J]. IEEE Transaction on Circuits and Systems, 1978, 25(9): 772-781. DOI: 10.1109/TCS.1978.1084534.

(編輯 王小唯， 苗秀芝)

High resolution approach of direction-of-arrival estimation based on twice virtual extension

ZHAO Yu, LI Wenxing, MAO Xiaojun

(College of Information and Communications Engineering, Harbin Engineering University, Harbin 150001, China)

In order to further improve the resolution of direction-of-arrival (DOA) estimation of antenna array, a high resolution approach of DOA estimation is proposed based on the fourth-order multiple signal classification (FO-MUSIC) approach. First, the antenna array is extended through the fourth-order moment of the

data. Then, the antenna array is extended through the conjugate value of the received data. When the expanded steering vector and the expanded covariance matrix are used for DOA estimation instead of the original steering vector and the original covariance matrix, virtual array elements are formed and the aperture of array is extended. It is shown by simulation results that, compared to the FO-MUSIC approach, the proposed method has a higher probability of target resolution, lower root mean square error and higher resolution. Hence, more virtual array elements are formed by the proposed approach based on twice virtual extension, and the resolution of DOA estimation of antenna array is effectively improved.

direction-of-arrival estimation; multiple signal classification; fourth-order statistics; virtual array; resolution

10.11918/j.issn.0367-6234.201612032

2016-12-08

國防預(yù)研項(xiàng)目(4010403020102， 4010103020103)

趙 宇(1987—)，男，博士研究生； 李文興(1960—)，男，教授，博士生導(dǎo)師

李文興，liwenxing@hrbeu.edu.cn

TN911.7

A

0367-6234(2017)05-0122-06