基于互質(zhì)陣虛擬陣列空間平滑的相干信號(hào)DOA估計(jì)方法

韋 娟, 嚴(yán)世安, 寧方立

(1. 西安電子科技大學(xué)通信工程學(xué)院, 陜西 西安 710071; 2. 西北工業(yè)大學(xué)機(jī)電學(xué)院, 陜西 西安 710072)

0 引 言

波達(dá)方向(direction of arrival,DOA)估計(jì)是一種陣列接收信號(hào)獲取信號(hào)入射方向信息的技術(shù),被廣泛應(yīng)用在雷達(dá)、聲納、無(wú)線通信等領(lǐng)域。目前多以非相干信號(hào)為研究對(duì)象,但在實(shí)際信號(hào)傳播過(guò)程中,由于多徑傳播和同頻干擾造成大量相干信號(hào)源存在,因此相干信號(hào)的DOA估計(jì)亦是一研究熱點(diǎn)。

均勻線陣相干信號(hào)DOA估計(jì)算法中,以空間平滑、Toeplitz矩陣重構(gòu)為基礎(chǔ)的解相干算法,均存在陣列孔徑損失問(wèn)題?？梢酝ㄟ^(guò)增加陣元數(shù)來(lái)提高DOA估計(jì)算法的角度自由度(degrees of freedom,DOF),但該方法會(huì)增加系統(tǒng)部署成本;亦可引入如互質(zhì)陣列、嵌套陣列等非均勻稀疏陣列構(gòu)建虛擬陣列以拓展陣列孔徑?；ベ|(zhì)陣列是由兩個(gè)陣元數(shù)分別為互質(zhì)數(shù)的均勻子陣相互穿插形成,結(jié)合了稀疏陣列的優(yōu)勢(shì)和質(zhì)數(shù)的性質(zhì),能夠突破奈奎斯特采樣定理的限制,陣列部署靈活,具有廣闊的應(yīng)用前景。

前向空間平滑(forward spatial smoothing,FOSS)算法是最早被提出來(lái)估計(jì)相干信號(hào)DOA的,但該方法孔徑損失嚴(yán)重。前后向空間平滑(forward and backward spatial smooth,FBSS)算法利用了后向子陣列,增加了平滑次數(shù),增大了陣列孔徑,擁有比FOSS算法更高的DOF,但對(duì)信號(hào)沒有進(jìn)行完全的解相干。Han等提出了一種類旋轉(zhuǎn)不變子空間信號(hào)參數(shù)估計(jì)(estimating signal parameters via rotational invariance techniques-like,ESPRIT-like)算法,將數(shù)據(jù)協(xié)方差矩陣的任意一行構(gòu)造成Toeplitz矩陣,能夠完全消除信號(hào)的相干性,利用ESPRIT-like算法進(jìn)行DOA估計(jì),但該算法估計(jì)精度低,亦存在陣列孔徑損失。Stoica等提出的最大似然(maximum likelihood,ML)估計(jì)算法以及Ottersten等提出的加權(quán)子空間擬合(weighted subspace fitting, WSF)算法均可以在不損失陣列孔徑的前提下對(duì)相干信號(hào)進(jìn)行DOA估計(jì),但這兩種算法包含多維非線性搜索,計(jì)算復(fù)雜度高,且初值設(shè)置對(duì)于DOA測(cè)向影響較大。Wang等提出的通過(guò)四階累積重構(gòu)Toeplitz矩陣進(jìn)行相干信號(hào)DOA估計(jì),DOA估計(jì)精度較高,但計(jì)算復(fù)雜度高。Zhang等提出了一種基于多重Toeplitz矩陣重構(gòu)的DOA估計(jì)方法,對(duì)相關(guān)矩陣的平方進(jìn)行加權(quán)求和,以形成滿秩等效數(shù)據(jù)協(xié)方差矩陣。該算法無(wú)需提前消噪,但低信噪比、低快拍條件下估計(jì)性能不佳。唐曉杰等提出了一種前后向多重Toeplitz矩陣重構(gòu)(forward and backward multi-Toeplitz,FBMT)算法,構(gòu)造了 Toeplitz 矩陣,通過(guò)誤差最小準(zhǔn)則構(gòu)造代價(jià)函數(shù),實(shí)現(xiàn)相干信號(hào)的 DOA 估計(jì)。該方法無(wú)需信源數(shù)先驗(yàn)信息,但陣列孔徑損失嚴(yán)重,低信噪比條件下估計(jì)性能較差。Peng等提出一種將協(xié)方差矩陣子矩陣對(duì)角線元素之和的共軛作為加權(quán)因子的自加權(quán)空間平滑(self-weighted spatial smoothing,SWSS)算法,信號(hào)估計(jì)精度較高,由于采用空間平滑預(yù)處理,依舊存在陣列孔徑損失。

基于此,對(duì)互質(zhì)陣虛擬陣列的連續(xù)部分進(jìn)行空域平滑,提出一種基于互質(zhì)陣列虛擬陣列空間平滑(coprime array virtual array spatial smoothing,CASS)的相干信號(hào)DOA估計(jì)方法,該算法不需要信源數(shù)先驗(yàn)信息,利用互質(zhì)陣虛擬陣列可以實(shí)現(xiàn)DOF拓展的特點(diǎn),搭載空間平滑算法進(jìn)行相干信號(hào)DOA估計(jì)。在低信噪比環(huán)境下對(duì)于信號(hào)的DOA估計(jì)精度以及魯棒性擁有較好表現(xiàn)。

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

1.1 互質(zhì)陣列接收信號(hào)模型

互質(zhì)陣列利用虛擬域信號(hào)處理方法,在虛擬域上形成一個(gè)增廣的虛擬陣列,對(duì)該虛擬陣列信號(hào)進(jìn)行與均勻線陣等價(jià)的統(tǒng)計(jì)信號(hào)處理,能夠有效提升DOA估計(jì)的DOF。

互質(zhì)陣列結(jié)構(gòu)如圖1所示,該陣列由兩個(gè)均勻線陣嵌套構(gòu)成。

圖1 互質(zhì)陣列Fig.1 Coprime array

其中,,為互質(zhì)的兩個(gè)數(shù),虛擬陣元間距=2?；ベ|(zhì)子陣1為2元均勻線陣,陣元間距為,子陣1陣元位置為{|0≤≤2-1}。互質(zhì)子陣2為元均勻線陣,陣元間距為,子陣2陣元位置為{|0≤≤-1}。

互質(zhì)陣列物理陣元真實(shí)位置為=[,,…,2+-2]。

互質(zhì)陣虛擬等效差聯(lián)合陣列如圖2所示。

圖2 互質(zhì)陣列虛擬陣列Fig.2 Virtual array of coprime array

虛擬陣元位置為={ ±(-)|0≤≤2-1, 0≤≤-1}。其中,虛擬陣列陣元位置-(+-1)與(+-1)之間無(wú)空洞,可以視為均勻線陣。

個(gè)信號(hào)入射到互質(zhì)陣列接收信號(hào)為

(1)

對(duì)于第個(gè)相干信號(hào)可表示為

()=(),=2,3,…,

式中:()為原始參照信號(hào);為相干系數(shù)。

1.2 空間前后向平滑處理

互質(zhì)陣列協(xié)方差矩陣為

(2)

對(duì)接收協(xié)方差矩陣進(jìn)行向量化處理:

(3)

式中:=vec (),為單位矩陣;=[()?(),…,()?()],?表示Kronecker內(nèi)積,(·)表示共軛。

通過(guò)在Toeplitz矩陣結(jié)構(gòu)中重新排列的元素可獲得虛擬陣列信號(hào)的等價(jià)協(xié)方差。由于信噪比和快拍數(shù)對(duì)協(xié)方差矩陣的影響,同一波程差對(duì)應(yīng)的協(xié)方差矩陣元素不同,故對(duì)同一波程差對(duì)應(yīng)的元素取平均值:

(4)

式中:(diff)表示同一波程差diff對(duì)應(yīng)的第個(gè)協(xié)方差矩陣元素;(diff)表示波程差diff相同的協(xié)方差元素總和。

在互質(zhì)虛擬陣列空洞位置插入天線,將聯(lián)合差陣列變?yōu)榫鶆蚓€陣后,假設(shè)空洞處天線無(wú)接收信號(hào)。根據(jù)波程差元素對(duì)協(xié)方差矩陣進(jìn)行擴(kuò)展,構(gòu)成Toeplitz矩陣∈(2-+1)×(2-+1)。

(5)

式中：1<,<2-+1。

盡管矩陣包含所有互質(zhì)陣列信號(hào)信息,但由于稀疏陣列特性,是低秩矩陣,把看作部分?jǐn)?shù)據(jù)信息缺失的協(xié)方差矩陣,利用低秩矩陣恢復(fù)理論,對(duì)矩陣中的零元素進(jìn)行近似填充。

創(chuàng)建優(yōu)化問(wèn)題:

(6)

式中:()是以復(fù)向量為第一列的Hermitian Toeplitz矩陣,∈(2-+1)×1;為協(xié)方差矩陣擬合誤差。

由于矩陣秩函數(shù)是非連續(xù)、非凸函數(shù),利用跡范數(shù)最小化對(duì)式(6)進(jìn)行凸松弛表示:

(7)

構(gòu)建映射矩陣,∈(2-+1)×(2-+1)。

(8)

由于,向量的首元素?zé)o法進(jìn)行共軛拓展,故將第一個(gè)元素初值先定義為實(shí)數(shù),經(jīng)CVX工具箱計(jì)算得到的整體優(yōu)化向量后,再對(duì)其初值進(jìn)行優(yōu)化。

經(jīng)凸優(yōu)化重構(gòu)后的矩陣()可看作由2-+1個(gè)陣元組成的均勻陣列的協(xié)方差矩陣。為提高算法估計(jì)精度,對(duì)重構(gòu)后的Toeplitz矩陣()進(jìn)行前后向空間平滑分成個(gè)子陣列。為避免虛擬陣列內(nèi)插零陣元造成協(xié)方差秩虧損,滑窗長(zhǎng)度須不小于連續(xù)差聯(lián)合子陣正向長(zhǎng)度(+-1)?；伴L(zhǎng)度設(shè)置為=+-1。

FOSS處理如圖3所示。

圖3 FOSSFig.3 FOSS

前向平滑時(shí),定義第個(gè)子陣列接收數(shù)據(jù)為

(9)

第個(gè)子陣接收數(shù)據(jù)協(xié)方差矩陣可表示為

(10)

FOSS后的協(xié)方差矩陣可通過(guò)各子陣的協(xié)方差矩陣求平均而得

(11)

(12)

式中:為置換矩陣,其反對(duì)角線元素均為1,其余位置全為0。

將FOSS矩陣及后向空間平滑矩陣取平均得到FBSS矩陣(∈(+-1)×(+-1));

(13)

由于除2以外的質(zhì)數(shù)均為奇數(shù),當(dāng)=2時(shí),為奇數(shù),則+-1為奇數(shù);當(dāng)>2時(shí),,均為奇數(shù),則+-1為奇數(shù)。即對(duì)于互質(zhì)陣列,子陣陣元數(shù)<時(shí),連續(xù)位置+-1均為奇數(shù)。則構(gòu)建的Toeplitz矩陣為奇數(shù)階。

Toeplitz矩陣可以看作陣元位置為{-,…,0,…,}均勻線性陣列的輸出協(xié)方差矩陣:

(14)

式中:=(-1)2。

對(duì)于相干信號(hào),又有

(15)

式中:

(16)

2 基于后向平滑的相干信號(hào)估計(jì)

2.1 重構(gòu)后向平滑Toeplitz矩陣

選取的第行:

(17)

式中:=[,1,…,,],且信號(hào)協(xié)方差矩陣中不會(huì)存在全零行,因此中不存在零元素。利用長(zhǎng)度為+1的滑窗截取的元素,得到行向量:

,=[(,-),(,1-),…,(,-)], 0≤≤

(18)

取第一個(gè)行向量,0:

(19)

式中:

利用向量重構(gòu)Toeplitz矩陣:

(20)

(21)

式中:+1,表示第條對(duì)角線元素為1,其余元素均為零的+1階方陣;=diag{}表示矩陣重構(gòu)后的虛擬信號(hào)協(xié)方差矩陣。

忽略噪聲影響,可以寫為

(22)

2.2 凸優(yōu)化去噪處理

對(duì)于第個(gè)信號(hào),必存在∈1×(+1)和其余-1個(gè)導(dǎo)向矢量張成的信號(hào)子空間正交:

(23)

式中：range{·}表示張成的子空間。

(24)

將式(24)代入式(22)得

(25)

根據(jù)正交性質(zhì),有

(26)

考慮到噪聲影響,式(26)可表示為

(27)

為減小誤差,利用每次平滑的子陣列構(gòu)造代價(jià)函數(shù),得到凸優(yōu)化模型:

(28)

(29)

構(gòu)建映射矩陣:

(30)

(31)

又

(32)

(33)

式(33)的根為

(34)

式中：(·)表示偽逆。

將式(34)代入式(28)得

(35)

將()()進(jìn)行特征值分解:

(36)

(37)

()=+1-max eig{()()}

(38)

式中:max eig{·}表示特征值分解后最大特征值。

構(gòu)建譜函數(shù)():

(39)

對(duì)()進(jìn)行譜搜索,對(duì)每個(gè)()()求其最大特征值,()譜峰對(duì)應(yīng)角度即為DOA估計(jì)方向。

算法步驟歸納如下:

根據(jù)式(2)計(jì)算陣列數(shù)據(jù)協(xié)方差矩陣=E{()()};

根據(jù)式(4)對(duì)虛擬陣列元素進(jìn)行整合,根據(jù)式(5)構(gòu)建矩陣;

利用CVX工具箱,根據(jù)式(7)重構(gòu);

對(duì)矩陣進(jìn)行后向平滑矩陣。用+1長(zhǎng)度的滑窗截取,根據(jù)式(17)得到;

根據(jù)式(30)、式(31)計(jì)算映射矩陣和();

根據(jù)式(39)計(jì)算空間譜函數(shù)(),通過(guò)譜峰搜索得到信號(hào)DOA。

2.3 算法DOF分析

由2+-1個(gè)陣元組成的互質(zhì)陣列經(jīng)Toeplitz重構(gòu)產(chǎn)生(),因平滑子陣陣元數(shù)目大于等于相干信號(hào)數(shù)時(shí),可以有效解相干。則對(duì)()進(jìn)行空間平滑時(shí),>。且由于對(duì)互質(zhì)陣虛擬陣列進(jìn)行了空洞填充處理,為保證信號(hào)信息的完整,滑窗長(zhǎng)度不小于虛擬連續(xù)陣列正向長(zhǎng)度(+-1),即≥+-1。則取滑窗長(zhǎng)度=+-1,經(jīng)前后向空間平滑預(yù)處理后,生成矩陣∈(2+1)×(2+1),其中=(+-2)2。

矩陣不存在秩虧損,對(duì)矩陣第行進(jìn)行空間平滑處理生成,此時(shí)滑窗長(zhǎng)度只需不小于信源數(shù)即可完成對(duì)的解相干?；伴L(zhǎng)度取最小值,即假設(shè)噪聲子空間僅由一條導(dǎo)向矢量構(gòu)成,滑窗長(zhǎng)度為+1。

綜上,2+-1個(gè)陣元組成的互質(zhì)陣列,其虛擬陣經(jīng)空間平滑后最多可以完成(+)2個(gè)相干角度估計(jì)。

3 實(shí)驗(yàn)分析

本節(jié)通過(guò)仿真實(shí)驗(yàn)將CASS算法與文獻(xiàn)[15]提出的FBSS算法、文獻(xiàn)[18-20]提出的ML算法、文獻(xiàn)[28]提出的FBMT算法、文獻(xiàn)[29]提出的SWSS算法進(jìn)行比較,對(duì)不同條件下的均方根誤差(root mean square error, RMSE)進(jìn)行分析。RMSE可以通過(guò)500次蒙特卡羅實(shí)驗(yàn)得到:

(40)

式中:為信源數(shù)。

設(shè)信號(hào)為零均值高斯信號(hào),噪聲為零均值高斯白噪聲。由于FBSS算法、ML算法、SWSS算法需要信號(hào)源數(shù)先驗(yàn)信息,而FBMT算法和CASS算法無(wú)需信號(hào)源數(shù)。故仿真實(shí)驗(yàn)中,FBSS算法、ML算法、SWSS算法的信號(hào)源數(shù)是已知的,FBMT算法和CASS算法信源數(shù)先驗(yàn)信息均未知。

為方便DOF的比較,CASS算法構(gòu)建互質(zhì)陣列子陣陣元數(shù)分別為2=6,=5,總陣元數(shù)為10。作為不損失陣列孔徑算法對(duì)照,ML算法采取10元均勻線陣。由于CASS算法、ML算法可以實(shí)現(xiàn)相干信號(hào)的最大測(cè)向數(shù)為9。FBSS算法、FBMT算法、SWSS算法構(gòu)建的均勻陣列陣元設(shè)置為17,即完成9個(gè)相干信號(hào)DOA估計(jì)的FMBT算法所需要的最小陣元數(shù)。

3.1 全相干信號(hào)的估計(jì)性能分析

3.1.1 不同信噪比條件下的算法性能

假設(shè)9個(gè)入射信號(hào)均相干,入射角度分別為[-66°,-42°,-30°,-15°, 0°, 10°, 25°, 42°, 56°]。相干系數(shù)分別為=[1, 02ejπ3, ejπ4, 04ejπ6, ejπ5, 06ej2π5, ejπ7, 08ej2π7, ej3π7],正則化參數(shù)=1.5×10,快拍數(shù)為500,掃描間隔為0.1°。全相干信號(hào)中信噪比對(duì)算法估計(jì)精度的影響如圖4所示。分析圖4可知, FMBT算法由于已經(jīng)達(dá)到其最大信號(hào)測(cè)量角度極限,高信噪比條件下測(cè)向性能與FBSS算法接近。ML算法由于測(cè)向信號(hào)較多,難以保證經(jīng)多次角度代入處理收斂到的局部最小值是真實(shí)信號(hào)值。SWSS算法由于利用信號(hào)協(xié)方差信息,在低質(zhì)量信號(hào)中提取出有效信息能力較強(qiáng),估計(jì)精度較高。在信噪比為-12～-4 dB條件下,CASS算法由于對(duì)原互質(zhì)陣列進(jìn)行同波程差數(shù)據(jù)整合以及前后向空間平滑,增強(qiáng)了信號(hào)強(qiáng)度,估計(jì)精度高。CASS算法由于填補(bǔ)虛擬陣列空洞,存在信息缺失,故在8 dB信噪比條件下估計(jì)性能略差。

圖4 全相干信號(hào)信噪比對(duì)不同算法估計(jì)精度的影響Fig.4 Influence of signal to noise ratio of fully coherent signal on the estimation accuracy of different algorithms

全相干信號(hào)-12 dB條件下各算法的歸一化功率譜如圖5所示。分析圖5可知,在低信噪比條件下,FBSS算法由于沒有進(jìn)行完全解相干,估計(jì)精度較低;ML算法在對(duì)多個(gè)相干信號(hào)進(jìn)行DOA估計(jì)時(shí),需要選定合適的初值,DOA估計(jì)精度受初值設(shè)置影響。FBMT算法沒有對(duì)信號(hào)信息進(jìn)行預(yù)處理,噪聲影響較大,不能準(zhǔn)確估計(jì)DOA;SWSS算法在較大噪聲干擾下,無(wú)法通過(guò)加權(quán)子空間算法估計(jì)準(zhǔn)確的DOA。CASS算法經(jīng)過(guò)預(yù)處理,去噪相對(duì)徹底,算法精度較高。

圖5 全相干信號(hào)-12 dB條件下的歸一化功率譜Fig.5 Normalized power spectrum of fully coherent signal at -12 dB

3.1.2 不同快拍數(shù)條件下的算法性能

設(shè)置信噪比為0 dB,其余仿真環(huán)境同第3.1.1節(jié)。全相干信號(hào)中快拍數(shù)對(duì)算法估計(jì)精度的影響如圖6所示。

圖6 全相干信號(hào)快拍數(shù)對(duì)不同算法估計(jì)精度的影響Fig.6 Influence of snapshots number of fully coherent signal on the estimation accuracy of different algorithms

分析圖6可知,FBSS算法估計(jì)性能受快拍數(shù)影響較大,且隨著快拍數(shù)增加,估計(jì)精度提升明顯。在快拍數(shù)為20時(shí),ML算法、FBMT算法、SWSS算法以及CASS算法估計(jì)精度遠(yuǎn)高于FBSS算法。且CASS算法在快拍數(shù)為20時(shí)依舊保持良好的算法精度,估計(jì)性能最佳。

3.2 部分相干信號(hào)的估計(jì)性能分析

3.2.1 不同信噪比條件下的算法性能

假設(shè)入射角度為[-66°, -42°, -30°, -15°, 0°, 10°, 25°, 42°, 56°]的9個(gè)信號(hào),前6個(gè)信號(hào)為相干信號(hào),后3個(gè)信號(hào)是與其他信號(hào)完全不相干的信號(hào)。正則化參數(shù)τ=1.5×10,快拍數(shù)為500,掃描間隔為0.1°。部分相干信號(hào)中信噪比對(duì)算法估計(jì)精度的影響如圖7所示。

圖7 部分相干信號(hào)信噪比對(duì)不同算法估計(jì)精度的影響Fig.7 Influence of signal to noise ratio of partially coherent signal on the estimation accuracy of different algorithms

分析圖7可知,對(duì)于部分相干信號(hào),FBSS算法在低信噪比條件下估計(jì)精度依舊較低。ML算法、FBMT算法與SWSS算法估計(jì)精度均低于CASS算法。CASS算法在低信噪比條件下估計(jì)角度均方誤差較小。

部分相干信號(hào)-12 dB條件下各算法的歸一化功率譜如圖8所示。

圖8 部分相干信號(hào)-12 dB條件下的歸一化功率譜Fig.8 Normalized power spectrum of partially coherent signal at -12 dB

分析圖8可知,在部分相干信號(hào)中,FBSS波峰尖銳,完成所有DOA估計(jì),但與真實(shí)角度偏差較大。 ML算法經(jīng)多維非線性搜索,收斂值不能保證為全局最佳。FBMT算法在已達(dá)到其相干信源數(shù)估計(jì)極限的情況下,功率譜函數(shù)波峰較平;SWSS算法經(jīng)自加權(quán)處理,在低信噪比條件下能夠有效去噪,波峰較尖銳。CASS算法由于凸優(yōu)化去噪相對(duì)徹底,其混合信號(hào)估計(jì)性能與全相干信號(hào)估計(jì)性能相對(duì)穩(wěn)定。

3.2.2 不同快拍數(shù)條件下的算法性能

設(shè)置信噪比為0 dB,其余仿真環(huán)境同第3.2.1節(jié)。部分相干信號(hào)中快拍數(shù)對(duì)算法估計(jì)精度的影響如圖9所示。

圖9 部分相干信號(hào)快拍數(shù)對(duì)算法估計(jì)精度的影響Fig.9 Influence of snapshots number of partially coherent signal on the estimation accuracy of algorithm

分析圖9可知,FBSS算法估計(jì)精度隨快拍數(shù)增加而提高。ML算法、FBMT算法、SWSS算法、CASS算法隨快拍數(shù)增加估計(jì)精度變化曲線較為平穩(wěn)。CASS算法估計(jì)性能最佳,且在低快拍條件下性能穩(wěn)定。

3.3 運(yùn)算時(shí)間比較

實(shí)驗(yàn)仿真環(huán)境為Matlab R2020a平臺(tái),Inter i7-9570H處理器,16 G內(nèi)存。仿真條件同第3.1.2節(jié)。各算法運(yùn)算時(shí)間比較如圖10所示。

圖10 不同算法運(yùn)算時(shí)間比較Fig.10 Comparison of operation time of different algorithms

分析圖10可知, FBSS算法只經(jīng)過(guò)空間平滑處理,運(yùn)算速度最快;ML算法由于需進(jìn)行多維搜索,且DOA估計(jì)數(shù)目愈多,運(yùn)算時(shí)間愈長(zhǎng);FBMT算法由于進(jìn)行多次前后向空間平滑處理及特征分解處理,運(yùn)算時(shí)間也相對(duì)較長(zhǎng)。SWSS算法由于需要對(duì)信號(hào)協(xié)方差進(jìn)行加權(quán)累加,隨著精度要求的提高,迭代次數(shù)變多,計(jì)算復(fù)雜度增高;CASS算法采用互質(zhì)陣列作為接收陣列,陣元數(shù)較少,且相較FBMT算法省去多次前向平滑處理以及映射陣列的特征值分解,計(jì)算復(fù)雜度較低,運(yùn)算時(shí)間相對(duì)較短。

4 結(jié) 論

針對(duì)現(xiàn)有基于完全相干信號(hào)的DOA估計(jì)算法陣列孔徑損失嚴(yán)重,且在低信噪比環(huán)境下DOA計(jì)算精度較低的問(wèn)題,提出了一種基于互質(zhì)陣虛擬陣列空間平滑的DOA估計(jì)算法。該算法將互質(zhì)陣列進(jìn)行空洞填充,協(xié)方差矩陣按照波程差大小重構(gòu)Toeplitz矩陣,然后再通過(guò)后向空間平滑算法二次構(gòu)建Toeplitz矩陣,利用構(gòu)建的與信號(hào)子空間最大信號(hào)特征值相關(guān)的代價(jià)函數(shù)來(lái)搜索譜函數(shù),得到DOA。該算法優(yōu)點(diǎn)如下:

(1) CASS算法在比FBSS、FBMT、SWSS算法少陣元的條件下完成了比以上算法更準(zhǔn)確的相干信號(hào)的DOA估計(jì)。同時(shí),CASS算法估計(jì)精度又高于同所需陣元數(shù)相同的ML算法。CASS算法系統(tǒng)架置成本低、部署靈活,且計(jì)算復(fù)雜度低、性價(jià)比高;

(2) CASS算法在低快拍數(shù)、低信噪比環(huán)境下對(duì)完全相干、部分相干信號(hào)DOA估計(jì)精度更高,具有更好的魯棒性;

(3) CASS算法估計(jì)相干信號(hào)DOA時(shí),不需要信源數(shù)先驗(yàn)信息,更符合實(shí)際應(yīng)用環(huán)境。

由于互質(zhì)陣列虛擬陣列法在信號(hào)分離時(shí)利用空洞填充,存在一定信號(hào)損失,且互質(zhì)陣列虛擬陣列算法是在虛擬域內(nèi)進(jìn)行的二階等價(jià)虛擬陣列信號(hào)處理,故該算法在較高信噪比環(huán)境下DOA估計(jì)精度遜色于FBMT算法、SWSS算法。高信噪比條件下CASS算法的DOA估計(jì)精度仍存在提升空間。