一種用于水下小尺度運(yùn)動(dòng)陣列的目標(biāo)波達(dá)方向估計(jì)方法

郭拓， 王英民， 張立琛

(西北工業(yè)大學(xué) 航海學(xué)院， 陜西 西安 710072)

一種用于水下小尺度運(yùn)動(dòng)陣列的目標(biāo)波達(dá)方向估計(jì)方法

郭拓， 王英民， 張立琛

(西北工業(yè)大學(xué) 航海學(xué)院， 陜西 西安 710072)

為解決水下小尺度運(yùn)動(dòng)陣列對(duì)目標(biāo)定位的孔徑和樣本不足問(wèn)題，采用合成孔徑方法實(shí)現(xiàn)孔徑的擴(kuò)展。通過(guò)分析相干信源情況下不同陣元數(shù)與樣本數(shù)之比時(shí)樣本協(xié)方差矩陣的譜分離特性，提出采用小快拍的主特征空間目標(biāo)波達(dá)方向估計(jì)方法。研究結(jié)果表明：仿真實(shí)驗(yàn)中，陣元數(shù)與樣本數(shù)之比為5時(shí)，該方法依然可以正確分辨4個(gè)目標(biāo)；水池試驗(yàn)中，陣元數(shù)與樣本數(shù)之比為4.8時(shí)，該方法同樣可以準(zhǔn)確分辨3個(gè)相鄰目標(biāo)。該方法不需要先驗(yàn)信息如信源個(gè)數(shù)等，解決了小尺度運(yùn)動(dòng)陣列在實(shí)際應(yīng)用時(shí)所暴露的孔徑和樣本方面的固有局限，可滿(mǎn)足實(shí)際應(yīng)用需求。

聲學(xué)； 小尺度運(yùn)動(dòng)陣； 主特征空間； 波達(dá)方向估計(jì)； 小快拍

Abstract: The underwater small-scale moving array is limited by its small aperture and sample size for direction of arrival (DOA) estimation of coherent target. The passive synthetic aperture technique is used to solve the problem of insufficient aperture. The spectral separation of sample covariance matrix is stu-died in the case of coherent target. A DOA estimation method for coherent target based on main feature space is proposed, which uses small snapshots. In the simulation, the proposed method can still distinguish the four targets correctly when the ratio of the number of sensors to the number of samples is 5; in the water tank experiment, the above ratio is 4.8, and it can identify three adjacent targets clearly. The proposed method can achieve better resolution for coherent targets in small samples, which meets the application needs of small-scale motions array of targeting, and requires no priori information of signal source number.

Key words: acoustics; small-scale moving array; main feature space; DOA estimation; small snapshot

0 引言

小尺度聲納基陣由于其體積小、成本低、基陣校準(zhǔn)和使用方便，廣泛應(yīng)用于魚(yú)雷、水雷、水下無(wú)人航行器等小尺度安裝平臺(tái)上。但該聲納平臺(tái)的運(yùn)動(dòng)特性導(dǎo)致在信號(hào)的局部平穩(wěn)過(guò)程中所獲得的樣本數(shù)較小，小尺度陣即陣列孔徑較小，因此基于小尺度運(yùn)動(dòng)陣列的目標(biāo)波達(dá)方向(DOA)估計(jì)成為一個(gè)關(guān)鍵問(wèn)題。

DOA估計(jì)作為陣列信號(hào)處理的一個(gè)重要分支，已經(jīng)廣泛用于雷達(dá)、聲納、通信、無(wú)源定位等多個(gè)領(lǐng)域[1-3]。陣列DOA估計(jì)經(jīng)過(guò)幾十年的發(fā)展，已經(jīng)形成了兩類(lèi)技術(shù)，即基于波束形成技術(shù)和基于子空間分解技術(shù)。其中：基于波束形成技術(shù)的典型方法有常規(guī)波束形成器(CBF)和最小方差無(wú)畸變響應(yīng)(MVDR)波束形成器；基于子空間分解技術(shù)的經(jīng)典方法[4]是多重信號(hào)分類(lèi)(MUSIC)算法。

上述傳統(tǒng)DOA估計(jì)算法大都基于孔徑較大的陣列，其估計(jì)性能隨著陣列孔徑的減小會(huì)迅速下降，且都需要估計(jì)基陣的統(tǒng)計(jì)協(xié)方差矩陣。在實(shí)際數(shù)據(jù)處理中都以樣本協(xié)方差矩陣代替統(tǒng)計(jì)協(xié)方差矩陣，在小樣本情況下如此估計(jì)出的統(tǒng)計(jì)協(xié)方差矩陣誤差很大，會(huì)導(dǎo)致DOA估計(jì)算法失效。特別是對(duì)于相干目標(biāo)，經(jīng)典的高分辨方法如MUSIC隨著信號(hào)相干程度的增加，空間分別率逐漸下降，甚至完全失效。因此針對(duì)水下小尺度運(yùn)動(dòng)陣列對(duì)相干目標(biāo)DOA估計(jì)成為亟待解決的問(wèn)題，以實(shí)現(xiàn)在小孔徑情況下采用小樣本進(jìn)行方位估計(jì)。

針對(duì)小樣本數(shù)問(wèn)題，為了獲得與統(tǒng)計(jì)協(xié)方差矩陣更加相近的樣本協(xié)方差矩陣，通常情況下會(huì)對(duì)樣本協(xié)方差矩陣進(jìn)行修正，目前已經(jīng)發(fā)展了很多修正方法，如：1)對(duì)角加載方法[5]、收縮因子方法[6-8]，二者都是通過(guò)對(duì)樣本協(xié)方差矩陣的特征值進(jìn)行加權(quán)操作以改進(jìn)其性能；2)前后向平均[9]和空間光滑技術(shù)[10]，通過(guò)移動(dòng)窗口虛擬出更多的快拍，以改進(jìn)樣本協(xié)方差矩陣估計(jì)；3)偽隨機(jī)噪聲再采樣技術(shù)[11-13]，通過(guò)人為生成偽隨機(jī)噪聲信號(hào)加在接收到的信號(hào)上，以改進(jìn)樣本協(xié)方差矩陣估計(jì)；4)子空間泄露分析方法[14]，通過(guò)MUSIC算法的估計(jì)改進(jìn)所需的樣本協(xié)方差矩陣，以減少樣本協(xié)方差矩陣方位估計(jì)時(shí)信號(hào)子空間的泄漏。上述方法的思路都是設(shè)法改進(jìn)樣本協(xié)方差矩陣，以使其與大樣本時(shí)的樣本協(xié)方差矩陣更加相近，但沒(méi)有考慮樣本協(xié)方差矩陣與樣本數(shù)的關(guān)系，以及樣本數(shù)不同時(shí)協(xié)方差矩陣本身所具有的不同特性。

本文首先采用被動(dòng)合成孔徑技術(shù)進(jìn)行陣列的孔徑擴(kuò)展，解決陣列的孔徑不足問(wèn)題；對(duì)于小樣本問(wèn)題將打破上述通過(guò)改進(jìn)樣本協(xié)方差矩陣來(lái)改進(jìn)方位估計(jì)的思路，而是通過(guò)研究不同快拍時(shí)樣本協(xié)方差矩陣本身具有的一些特性如譜的分離特性，提出基于小樣本協(xié)方差矩陣的主特征空間DOA估計(jì)方法。

1 陣列接收信號(hào)模型

如圖1所示，水平線(xiàn)陣由N個(gè)陣元組成，假設(shè)有D個(gè)相互獨(dú)立的平面波入射，陣列第k次快拍X(k)為

X(k)=A(θ)S(k)+N(k),

(1)

式中：X(k)=[x1(k),x2(k),…,xN(k)]T；A(θ)為陣列流形矩陣，A(θ)=[a(θ1),a(θ2),…,a(θD)]，θj為第j個(gè)平面波的水平入射角，a(θj)=[1,ej2πdsin (θj)/λ,…,ej2π(N-1)dsin (θj)/λ]為第j個(gè)入射平面波的導(dǎo)向向量,d為陣元間距，其值等于波長(zhǎng)λ的一半 ；S(k)=[s1(k),s2(k),…,sD(k)]T,sj(k)為第j個(gè)平面波信號(hào)的第k次快拍；N(k)=[n1(k),n2(k),…,nN(k)]為陣列接收到的加性高斯白噪聲，ni(k)為第i個(gè)陣元接收到的噪聲。

圖1 陣列接收信號(hào)模型Fig.1 Array signal model

當(dāng)信號(hào)被假設(shè)具有零均值時(shí)，陣列接收信號(hào)的協(xié)方差矩陣R=E[X(k)XH(k)]，即協(xié)方差矩陣需要通過(guò)求集總平均得到。而實(shí)際應(yīng)用中只能得到有限的采樣數(shù)據(jù)，如有M次觀(guān)測(cè)樣本X(1),X(2),…,X(M)，則用這M個(gè)樣本估計(jì)的采樣協(xié)方差矩陣為

(2)

2 基于主特征空間的DOA估計(jì)方法

2.1 樣本協(xié)方差矩陣特征的譜分離特性

令y=N/M，y表示陣元數(shù)與樣本數(shù)之比，文獻(xiàn)[15-17]采用隨機(jī)矩陣?yán)碚撗芯繕颖緟f(xié)方差矩陣譜的極限分布，研究表明樣本協(xié)方差矩陣具有譜分離特性。文獻(xiàn)[18]詳細(xì)研究了樣本協(xié)方差矩陣的譜分離特性，并給出樣本協(xié)方差矩陣的譜隨y變化的特點(diǎn)，研究了y為1.00(樣本數(shù)等于陣元數(shù))、y為0.01(樣本數(shù)是陣元數(shù)的100倍)以及在y由1.00逐漸變小至0.01的過(guò)程中的特征值分布表現(xiàn)形式，研究發(fā)現(xiàn)：特征值分布表現(xiàn)為一簇一簇的，每簇中有幾個(gè)大小相近的特征值，其中每個(gè)簇的中心漸近于統(tǒng)計(jì)協(xié)方差矩陣R的一個(gè)真實(shí)特征值。每個(gè)簇的大小取決于y，即估計(jì)R的樣本數(shù)，如果每個(gè)陣元所取的樣本數(shù)小，則每個(gè)簇將會(huì)很大，的許多特征值將匯集在一個(gè)簇中，無(wú)法分離；反之，如果每個(gè)陣元所取的樣本數(shù)足夠大，則每個(gè)簇只對(duì)應(yīng)R的一個(gè)真實(shí)特征值，這樣R的譜就得到很好的分離。

仿真實(shí)驗(yàn)采用30個(gè)陣元，4個(gè)相干目標(biāo)，噪聲為高斯白噪聲，可得如圖2所示的特征值隨不同y的變化曲線(xiàn)。由圖2可知：當(dāng)y很小時(shí)，即樣本數(shù)遠(yuǎn)遠(yuǎn)大于陣元數(shù)即y很小時(shí)，的特征值會(huì)在最大特征值序數(shù)處陡峭升高，即最大特征值遠(yuǎn)大于其他特征值；當(dāng)y較大時(shí)，特征值不再陡峭升高，緊鄰最大特征值處還有較大的特征值。圖2中特征值隨y變化的特征，是相干源情況下樣本協(xié)方差譜分離特性的具體體現(xiàn)。從子空間的角度來(lái)解釋上述相干源時(shí)的譜分離特性，當(dāng)樣本數(shù)遠(yuǎn)遠(yuǎn)大于陣元數(shù)即y很小時(shí)，樣本協(xié)方差矩陣得到較好的估計(jì)，相干信號(hào)的信息全部集中在一個(gè)最大特征值所對(duì)應(yīng)的子空間上；當(dāng)y較大時(shí)，樣本協(xié)方差矩陣沒(méi)有得到很好的估計(jì)，此時(shí)信號(hào)子空間發(fā)生了泄漏，使得除了最大特征值之外還出現(xiàn)一些連續(xù)的較大特征值，但是總體而言所有相干信號(hào)的主要信息仍然集中于最大特征值所對(duì)應(yīng)的子空間中。本文將這個(gè)最大特征值所對(duì)應(yīng)的子空間稱(chēng)為主特征空間，而不是MUSIC算法所假設(shè)的那樣，即信號(hào)子空間就是s個(gè)信源所對(duì)應(yīng)的s個(gè)較大特征值所構(gòu)成的空間。通過(guò)上述對(duì)相干源情況下不同樣本數(shù)時(shí)譜分離特性的仿真實(shí)驗(yàn)與分析，顯然信號(hào)信息主要集中在主特征空間而不是MUSIC算法所使用的s個(gè)較大特征值構(gòu)成的空間，從而很好地解釋了MUSIC算法對(duì)相干目標(biāo)DOA估計(jì)時(shí)失效的原因。

圖2 相干目標(biāo)的譜分離性Fig.2 Spectral separation of coherent target

綜上所述，無(wú)論快拍數(shù)大小，目標(biāo)為相干信號(hào)時(shí)其樣本協(xié)方差矩陣具有很好的譜分離特性，下面將基于譜分離特性提出相干信號(hào)的目標(biāo)DOA估計(jì)方法。

2.2 基于被動(dòng)合成孔徑的主特征空間DOA估計(jì)方法

被動(dòng)合成孔徑聲納有幾種經(jīng)典算法，本文采用較實(shí)用的擴(kuò)展拖曳陣測(cè)量(ETAM)方法[19]。

假設(shè)目標(biāo)是窄帶聲源，f0為聲源頻率，則目標(biāo)信號(hào)可表示為s(tg)=Aexp (j2πf0tg)，A為信號(hào)幅度，tg為第g個(gè)時(shí)刻點(diǎn)。tg時(shí)刻第i個(gè)陣元接收到θj的信號(hào)可表示為

(3)

式中：f為目標(biāo)頻率；c為水下聲速；εi,g為噪聲。

由于陣列與目標(biāo)之間的相對(duì)運(yùn)動(dòng)，接收到的目標(biāo)頻率f為

f=f0(1±vsinθj/c)，

(4)

式中：v為陣列與目標(biāo)之間的相對(duì)運(yùn)動(dòng)速度。

將(4)式代入(3)式，得

(5)

將(5)式化簡(jiǎn)，略去sinθ的2階項(xiàng)，整理得

(6)

如果陣列不移動(dòng)，則tg+τ時(shí)刻(τ為時(shí)間間隔)第i個(gè)陣元接收到的信號(hào)可表示為

(7)

如果選擇合適的間隔時(shí)間τ，使得vτ=qd，即運(yùn)動(dòng)的距離正好等于q個(gè)陣元間距d，則tg+τ時(shí)刻第i個(gè)陣元接收到的信號(hào)變換為

(8)

而tg時(shí)刻第i+q個(gè)陣元接收到的信號(hào)為

(9)

比較(8)式與(9)式發(fā)現(xiàn)，二式只相差一個(gè)相位補(bǔ)償項(xiàng)，說(shuō)明可以通過(guò)第i個(gè)陣元的運(yùn)動(dòng)來(lái)虛擬出tg時(shí)刻的第i+q個(gè)陣元。當(dāng)陣列勻速運(yùn)動(dòng)時(shí)，時(shí)刻tg的第N-q+p個(gè)陣元與時(shí)刻tg+τ的第p個(gè)陣元在空間中同一個(gè)位置，即不同的陣元在不同的時(shí)刻接收空間相同位置上的信號(hào)。N-q為兩個(gè)時(shí)刻重疊陣元數(shù)，利用N-q個(gè)重疊陣元計(jì)算兩個(gè)時(shí)刻的相位補(bǔ)償，q即為一次合成孔徑虛擬出的陣元數(shù)，N為陣列物理陣元個(gè)數(shù)，經(jīng)過(guò)u次合成則可以虛擬出q×u個(gè)陣元，經(jīng)過(guò)u次孔徑合成的總陣元數(shù)為N+q×u.

經(jīng)過(guò)虛擬孔徑合成后，物理陣和虛擬陣列接收到的信號(hào)可表示為

Xe(k)=
[x1(k),x2(k),…,xN(k),xN+1(k)，…，xN+q×u(k)]T.

(10)

經(jīng)ETAM算法進(jìn)行孔徑擴(kuò)展后的實(shí)際使用樣本協(xié)方差矩陣可以用(11)式估計(jì)：

(11)

(12)

式中：λi為e的特征值；ei為該特征值所對(duì)應(yīng)的特征向量；E為以所有特征向量按列排所構(gòu)成的特征矩陣。

2.1節(jié)已將樣本協(xié)方差矩陣的最大特征值所對(duì)應(yīng)的空間定義為主特征空間，此處定義只包括最大特征值所對(duì)應(yīng)的空間，并非主特征空間只包含最大特征值所對(duì)應(yīng)的空間，也可以根據(jù)信源的特點(diǎn)擴(kuò)展主特征空間。定義e可以劃分為

(13)

式中：λm為主特征空間的特征值；em為與其對(duì)應(yīng)的特征向量；Em為主特征空間特征向量構(gòu)成的矩陣；Eo為特征空間E中去除主特征后的特征信息矩陣。

根據(jù)子空間理論，信號(hào)子空間與噪聲子空間正交，且導(dǎo)向向量張成的空間與信號(hào)子空間張成的空間是同一空間，即Ena(θ)=0，Esa(θ)=1，其中En為噪聲子空間，Es為信號(hào)子空間。

主特征空間中包含的主要是信號(hào)信息，即主特征子空間Em中主要包含的是Es空間的信息，當(dāng)樣本數(shù)較小時(shí)會(huì)發(fā)生子空間泄漏，有部分En的信息進(jìn)入Em，使Ema(θ)=1，則arccos(Ema(θ))=0，故定義其倒數(shù)為主特征空間方法的方位譜估計(jì)，即

(14)

下面將合成孔徑后，采用主特征空間方法進(jìn)行方位估計(jì)的方法簡(jiǎn)稱(chēng)為ETAM-MFS.

3 仿真實(shí)驗(yàn)

通過(guò)數(shù)值仿真模擬水下小尺度運(yùn)動(dòng)陣列目標(biāo)定位，首先采用被動(dòng)合成孔徑算法ETAM進(jìn)行孔徑擴(kuò)展，然后采用主特征空間DOA估計(jì)方法估計(jì)目標(biāo)方位。仿真采用一個(gè)半波長(zhǎng)均勻分布的6元線(xiàn)陣模擬小尺度陣列，陣列運(yùn)動(dòng)速度為0.07 m/s，入射的信號(hào)為4個(gè)遠(yuǎn)場(chǎng)窄帶信號(hào)，頻率為6.25 kHz，采樣頻率50 kHz，4個(gè)信號(hào)的入射角分別為-20°、-10°、0°、10°，信噪比皆為0 dB，噪聲為復(fù)高斯白噪聲，快拍數(shù)取6，ETAM算法經(jīng)過(guò)8次合成得到24個(gè)虛擬陣元，總線(xiàn)陣為30個(gè)陣元，陣元數(shù)與樣本數(shù)之比y為5.

圖3為CBF算法給出的目標(biāo)DOA估計(jì)，圖4為MVDR算法給出的目標(biāo)DOA估計(jì)，圖5為MUSIC算法給出的目標(biāo)DOA估計(jì)，圖6為ETAM-MFS方法給出的目標(biāo)DOA估計(jì)。由圖3～圖6可以發(fā)現(xiàn)，在30個(gè)陣元、6快拍情形下，CBF方法、MVDR方法與MUSIC方法都不能正確估計(jì)出目標(biāo)方位，而ETAM-MFS方法可以正確分辨出目標(biāo)方位。由此可見(jiàn)，本文提出的方法可以在小快拍情況下很好地分辨相干目標(biāo)，且不需要像MUSIC算法那樣預(yù)先擁有目標(biāo)個(gè)數(shù)的先驗(yàn)信息，這是相干源時(shí)樣本協(xié)方差矩陣的譜分離特性帶來(lái)的好處，即信號(hào)信息全部集中于主特征空間，只有很少的信號(hào)泄漏于其他空間。

圖3 CBF方法DOA估計(jì)Fig.3 DOA estimation of CBF

圖4 MVDR方法DOA估計(jì)Fig.4 DOA estimation of MVDR

圖5 MUSIC方法DOA估計(jì)Fig.5 DOA estimation of MUSIC

圖6 ETMA-MFS方法DOA估計(jì)Fig.6 DOA estimation of ETMA-MFS

4 水池試驗(yàn)

為了驗(yàn)證所提算法在小尺度運(yùn)動(dòng)陣列小樣本的情況下能夠分辨相干目標(biāo)，在消聲水池(長(zhǎng)20 m, 寬8 m，深7 m，6面消聲)進(jìn)行了相關(guān)試驗(yàn)。試驗(yàn)采用由6個(gè)水聽(tīng)器組成的均勻線(xiàn)陣，設(shè)計(jì)頻率為6.25 kHz，水平布放于水下2.5 m深處，陣元間距為半波長(zhǎng)0.12 m，陣列勻速運(yùn)動(dòng)速度為0.07 m/s，使用電機(jī)勻速拉動(dòng)；3個(gè)發(fā)射源發(fā)射頻率均為6.25 kHz的連續(xù)波信號(hào)，與接收陣在同一深度，3個(gè)源按遠(yuǎn)場(chǎng)條件布放，與陣列形成的角度分別為24°、33°、44° ；采樣頻率為50 kHz. 首先采用被動(dòng)合成孔徑算法ETAM合成6次，共得到24個(gè)陣元，快拍數(shù)取5，陣元數(shù)與樣本數(shù)之比y為4.8. 試驗(yàn)系統(tǒng)如圖7所示，試驗(yàn)現(xiàn)場(chǎng)如圖8所示。

圖7 水池試驗(yàn)系統(tǒng)Fig.7 Water tank experimental system

圖8 水池試驗(yàn)現(xiàn)場(chǎng)Fig.8 Water tank experimental site

圖9～圖12為采用被動(dòng)合成孔徑ETAM算法得到24個(gè)虛擬陣元后，分別使用CBF方法、MVDR方法、MUSIC方法與MFS方法得到的目標(biāo)DOA估計(jì)，圖中目標(biāo)的真實(shí)方位用黑點(diǎn)標(biāo)出。由圖9可知，CBF方法不能清晰分辨3個(gè)目標(biāo)，且具有明顯的方位誤差，其中物理陣所對(duì)應(yīng)的是實(shí)際物理6個(gè)陣元沒(méi)有進(jìn)行合成孔徑得到的目標(biāo)DOA估計(jì)。由圖10可知，MVDR方法只能估計(jì)出3個(gè)目標(biāo)中的2個(gè)。由圖11可知，MUSIC方法也只能估計(jì)出3個(gè)目標(biāo)中的2個(gè)。從圖12可以發(fā)現(xiàn)，ETAM-MFS算法可以很好地分辨3個(gè)目標(biāo)，且沒(méi)有方位誤差。與其他3種方法相比，首先，本文所提出的方法不需要做協(xié)方差矩陣求逆運(yùn)算，減少了計(jì)算量；其次，不需要先驗(yàn)信息如信源個(gè)數(shù)；最后，最重要的是，本文所提出的方法在小樣本情形下表現(xiàn)優(yōu)異，能夠準(zhǔn)確分辨目標(biāo)。

圖9 水池試驗(yàn)合成孔徑前后DOA估計(jì)對(duì)比Fig.9 DOA estimations before and after synthetic aperture in water tank experiment

圖10 水池試驗(yàn)MVDR方法DOA估計(jì)Fig.10 DOA estimation of MVDR in water tank experiment

圖11 水池試驗(yàn)MUSIC算法DOA估計(jì)Fig.11 DOA estimation of MUSIC in water tank experiment

圖12 水池試驗(yàn)ETAM-MFS方法DOA估計(jì)Fig.12 DOA estimation of ETAM-MFS in water tank experiment

5 結(jié)論

為解決水下小尺度運(yùn)動(dòng)陣列的目標(biāo)DOA估計(jì)中的孔徑和樣本不足問(wèn)題，本文首先采用被動(dòng)合成孔徑算法ETAM彌補(bǔ)小尺度陣的孔徑不足，然后在對(duì)樣本協(xié)方差矩陣譜分離特性分析的基礎(chǔ)上提出了小快拍主特征空間DOA估計(jì)方法。在分析相干目標(biāo)樣本協(xié)方差矩陣的譜分離特性時(shí)，還得出了MUSIC算法對(duì)相干目標(biāo)分辨失效的原因。通過(guò)仿真實(shí)驗(yàn)和水池試驗(yàn)，將本文所提出的方法與傳統(tǒng)DOA估計(jì)方法(CBF算法、MVDR算法以及MUSIC算法)進(jìn)行了對(duì)比，結(jié)果發(fā)現(xiàn)：當(dāng)快拍數(shù)很小時(shí)，試驗(yàn)中陣元數(shù)與快拍之比接近5時(shí)，該方法依然可以準(zhǔn)確地分辨多個(gè)目標(biāo)，而其他3種方法或是出現(xiàn)多個(gè)偽目標(biāo)或是有目標(biāo)不能估計(jì)出。

本文所提出的方法實(shí)現(xiàn)了水下小尺度運(yùn)動(dòng)陣列的目標(biāo)方位估計(jì)，主要體現(xiàn)在：1)通過(guò)合成孔徑克服了小尺度陣列孔徑的不足；2)能夠?qū)崿F(xiàn)小快拍相干目標(biāo)的DOA估計(jì)，可以很好地彌補(bǔ)MUSIC算法的缺陷，且不需要先驗(yàn)信息如信源個(gè)數(shù)。該方法適合于水下小尺度運(yùn)動(dòng)陣列目標(biāo)定位的應(yīng)用需求，期望能在無(wú)人水下航行器、魚(yú)雷等應(yīng)用中有很好的應(yīng)用前景。

References)

[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] 吳振, 戴繼生, 朱湘臨, 等. 互耦效應(yīng)下一種基于實(shí)值稀疏表示的波達(dá)方向估計(jì)算法[J]. 兵工學(xué)報(bào), 2015, 36(2):294-298. WU Zhen,DAI Ji-sheng, ZHU Xiang-lin, et al. A real-valued sparse representation method for DOA estimation with unknown mutual coupling[J]. Acta Armamentarii, 2015, 36(2):294-298.(in Chinese)

[3] 鄒吉武, 孫大軍. 線(xiàn)陣雙基地聲納波束零點(diǎn)形成MUSIC算法[J]. 兵工學(xué)報(bào), 2010, 31(3):364-368. ZOU Ji-wu, SUN Da-jun. MUSIC algorithm of beam null forming on linear array of bi-static sonar[J]. Acta Armamentarii, 2010, 31(3):364-368.(in Chinese)

[4] 王永良, 陳輝, 彭應(yīng)寧, 等. 空間譜估計(jì)理論與算法[M]. 北京:清華大學(xué)出版社, 2004:2-5. WANG Yong-liang, CHEN Hui, PENG Ying-ning, et al. Theory and algorithm of spatial spectrum estimation[M]. Beijing: Tsinghua University Press, 2004:2-5.(in Chinese)

[5] Carlson B D. Covariance matrix estimation errors and diagonal loading in adaptive arrays[J]. IEEE Transactions on Aerospace Electronic Systems, 1988, 24(4):397-401.

[6] Chen Y, Wiesel A, Eldar Y C, et al. Shrinkage algorithms for MMSE covariance estimation[J]. IEEE Transactions on Signal Processing, 2010, 58(10):5016-5029.

[7] Tong J, Schreier P J, Guo Q, et al. Shrinkage of covariance matrices for linear signal estimation using cross-validation[J]. IEEE Transactions on Signal Processing, 2016, 64(11):2965-2975.

[8] Auguin N, Morales-Jimenez D, Mckay M, et al. Robust shrinkage M-estimators of large covariance matrices[C]∥Proceedings of 2016 IEEE Statistical Signal Processing Workshop. Palma de Mallorca, Baleares, Spain:IEEE, 2016:1-4.

[9] Pillai S U, Kwon B H. Forward/backward spatial smoothing techniques for coherent signal identification[J]. IEEE Transactions on Acoustics Speech & Signal Processing, 1989, 37(1):8-15.

[10] Pal P, Vaidyanathan P P. A novel array structure for directions-of-arrival estimation with increased degrees of freedom[C]∥Proceedings of 2010 IEEE International Conference on Acoustics Speech and Signal Processing. Dallas, TX, US:IEEE, 2010:2606-2609.

[11] Gershman A B, Bohme J F. Improved DOA estimation via pseudo-random resampling of spatial spectrum[J]. IEEE Signal Processing Letters, 1997, 4(2):54-57.

[12] Vasylyshyn V. Removing the outliers in root-MUSIC via pseudo-noise resampling and conventional beamformer[J]. Signal Processing, 2013, 93(12):3423-3429.

[13] Qian C, Huang L, So H C. Improved unitary root-MUSIC for DOA estimation based on pseudo-noise resampling[J]. IEEE Signal Processing Letters, 2014, 21(2):140-144.

[14] Shaghaghi M, Vorobyov S A. Subspace leakage analysis and improved DOA estimation with small sample size[J]. IEEE Tran-sactions on Signal Processing, 2015, 63(12):3251-3265.

[15] Bai Z D, Silverstein J W. No eigenvalues outside the support of the limiting spectral distribution of large-dimensional sample covariance matrices[J]. Annals of Probability, 1998, 26(1):316-345.

[16] Bai Z D, Silverstein J W. Exact separation of eigenvalues of large dimensional sample covariance matrices[J]. Annals of Probability, 1999, 27(3):1536-1555.

[17] 李華, 白志東, 肖玉山. 大維隨機(jī)矩陣的漸進(jìn)特征[J]. 東北師大學(xué)報(bào):自然科學(xué)版, 2014, 46(4):1-8. LI Hua, BAI Zhi-dong, XIAO Yu-shan. The asymptotic properties of the large dimension random matrix[J]. Journal of Northeast Normal University:Natural Science Edition, 2014, 46(4): 1-8.(in Chinese)

[18] Mestre X, Lagunas M. Modified subspace algorithms for DoA estimation with large arrays[J]. IEEE Transactions on Signal Processing, 2008, 56(2): 598-614.

[19] Stergiopoulos S, Sullivan E J. Extended towed array processing by an overlap correlator[J]. Journal of the Acoustical Society of America, 1989, 86(1): 158-171.

NovelTargetDirection-of-arrivalEstimationMethodforUnderwaterSmall-scaleMovingArray

GUO Tuo, WANG Ying-min, ZHANG Li-chen

(School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an 710072, Shaanxi, China)

TB566

A

1000-1093(2017)09-1779-07

10.3969/j.issn.1000-1093.2017.09.015

2017-01-17

國(guó)家自然科學(xué)基金項(xiàng)目(61401362)

郭拓(1986—)，男，博士研究生。E-mail:guotuonwpu@126.com

王英民(1963—)，男，教授，博士生導(dǎo)師。E-mail:ywang@nwpu.edu.cn