基于去特征處理的信源數(shù)估計算法

吳 娜,司偉建,焦淑紅,吳 迪

(哈爾濱工程大學(xué)信息與通信工程學(xué)院,黑龍江哈爾濱150001)

基于去特征處理的信源數(shù)估計算法

吳 娜,司偉建,焦淑紅,吳 迪

(哈爾濱工程大學(xué)信息與通信工程學(xué)院,黑龍江哈爾濱150001)

在實際測向系統(tǒng)中,由于陣列中各陣元的相互干擾,陣列接收的噪聲為空域相關(guān)色噪聲,由于系統(tǒng)熱噪聲和散彈噪聲的存在,導(dǎo)致陣列接收的噪聲同時包含空域相關(guān)色噪聲與白噪聲。傳統(tǒng)針對空域相關(guān)色噪聲的信源數(shù)估計算法在白噪聲背景下的估計性能往往下降。針對這一問題,提出了一種基于去特征處理的信源數(shù)估計方法。該方法對數(shù)據(jù)協(xié)方差矩陣進(jìn)行去特征處理,得到一組與特征值相對應(yīng)的新數(shù)據(jù)協(xié)方差矩陣,利用新數(shù)據(jù)協(xié)方差矩陣在信號和噪聲子空間上的投影矩陣構(gòu)造判決函數(shù),進(jìn)而實現(xiàn)信源數(shù)估計。仿真實驗及對比數(shù)據(jù)表明該方法在包含有白噪聲成分的色噪聲背景下能夠有效進(jìn)行信源數(shù)估計,其突出優(yōu)勢表現(xiàn)在小快拍數(shù)、低信噪比條件下仍能準(zhǔn)確估計信源數(shù),而沒有明顯增加算法的運算量。

信源數(shù)估計;色噪聲;去特征處理;判決函數(shù)

0 引 言

空間譜估計技術(shù)是陣列信號處理領(lǐng)域的一個重要分支,一直以來因其優(yōu)越的空域參數(shù)特性而備受學(xué)者青睞。然而,當(dāng)信源數(shù)與實際信源數(shù)不符時,大多數(shù)空間譜估計算法的測向性能嚴(yán)重下降甚至失效。因此信源數(shù)的準(zhǔn)確估計是空間譜估計算法實現(xiàn)其高分辨性能的必要前提。

目前比較有代表性的信源數(shù)估計算法有信息論方法[12]、蓋氏圓(Gerschgorin’s disk estimation,GDE)方法[3]、平滑秩法[4]、矩陣分解法[5]及正則相關(guān)技術(shù)(canonical correlate technology,CCT)[67]等方法。基于信息論的方法包括Akaike信息論(Akaike’s information crilerion, AIC)準(zhǔn)則、最小描述長度(minimum description length, MDL)準(zhǔn)則等方法,這類方法的優(yōu)點在于計算量較小,但是低信噪比條件下存在欠估計問題、高信噪比條件下存在過估計問題,且在色噪聲背景下算法失效。為了將信息論準(zhǔn)則推廣應(yīng)用到色噪聲背景下,文獻(xiàn)[8]提出一種對角加載技術(shù)平滑色噪聲背景下的信源數(shù)估計方法,即MAIC法。文獻(xiàn)[9]提出一種基于對角加載的MDL準(zhǔn)則信源數(shù)估計方法,即MMDL法。這些修正的信息論方法可以在一定程度上有效克服這些缺點。GDE方法在色噪聲背景下能夠準(zhǔn)確估計信源數(shù),但是該方法在低信噪比小快拍數(shù)條件下效果較差,且該方法在求解過程中得到的數(shù)據(jù)不能應(yīng)用于后續(xù)處理而造成運算效率較低。CCT在低信噪比及色噪聲背景下可以有效估計信源數(shù)目,但該算法需要兩個相互獨立的空間分離陣列,降低了陣列有效孔徑,進(jìn)而限制了該算法的實際推廣應(yīng)用[10- 19]。

本文提出一種基于去特征處理的信源數(shù)估計方法,該方法對數(shù)據(jù)協(xié)方差矩陣進(jìn)行特征分解,然后對其進(jìn)行去特征值處理,構(gòu)造新的數(shù)據(jù)協(xié)方差矩陣,將新的數(shù)據(jù)協(xié)方差矩陣分別投影到信號子空間和噪聲子空間,對投影得到的數(shù)據(jù)進(jìn)行統(tǒng)計分類,最后根據(jù)判據(jù)函數(shù)值的大小來判定信源個數(shù)。計算機仿真實驗表明本文算法可以實現(xiàn)空域相關(guān)色噪聲[10]及白噪聲并存背景下的信源數(shù)估計,并且在低信噪比、低快拍數(shù)條件下該算法的信源數(shù)估計準(zhǔn)確率仍然很高。

1 信號模型

本節(jié)以均勻線陣為例假設(shè)D個相互獨立的遠(yuǎn)場窄帶信號入射到空間M元均勻線陣,各信號源的到達(dá)方向角分別為(θ1,θ2,…,θD),則天線陣的接收數(shù)據(jù)的矢量形式為

式中,X(t)=[x1(t),x2(t),…,xM(t)]T為陣列的M×1維接收數(shù)據(jù)矢量;S(t)=[s1(t),s2(t),…,sD(t)]T為信號的D×1維矢量;N(t)=[n1(t),n2(t),…,nM(t)]T為M×1維噪聲矢量;A=[a(θ1),a(θ2),…,a(θD)]為陣列的M×D維導(dǎo)向矢量陣,其中a(θi)(i=1,2,…,D)為導(dǎo)向矢量,具體表達(dá)式如下:

式中,ω為信號源載波角頻率;d為相鄰陣元的間距;c為光速。

設(shè)信號與噪聲相互獨立,陣列輸出的數(shù)據(jù)協(xié)方差矩陣為

式中,RS=E[S(t)SH(t)]為信號協(xié)方差矩陣;RN=E[N(t)· NH(t)]為噪聲協(xié)方差矩陣。

對數(shù)據(jù)協(xié)方差矩陣RX進(jìn)行特征值分解,有

式中,Σ=diag(λ1,λ2,…,λM)為特征值組成的對角陣,且滿足λ1≥λ2≥…≥λM;U=[e1,e2,…,eM],ei(i=1,2,…,M)為特征值λi對應(yīng)的特征向量;ΣS=diag(λ1,λ2,…,λD)是由前D個較大特征值組成的對角陣;US=[e1,e2,…,eD]是由前D個較大特征值對應(yīng)的特征向量張成的信號子空間;ΣN= diag(λD+1,λD+2,…,λM)是由后M-D個較小的特征值組成的對角陣;UN=[eD+1,eD+2,…,eM]是由后M-D個較小特征值對應(yīng)的特征向量張成的噪聲子空間。

將陣列導(dǎo)向矢量陣按行分塊可以得到A的另一種形式:

式中

式中,T表示矩陣轉(zhuǎn)置。

則陣列輸出的數(shù)據(jù)協(xié)方差矩陣可以表示為

已知各信號源相互獨立,則RS為滿秩對角陣,因此有

2 算法原理

本節(jié)首先給出了去特征處理的定義。數(shù)據(jù)協(xié)方差矩陣的特征值可以分為對應(yīng)信號子空間的特征值和對應(yīng)噪聲子空間的特征值,根據(jù)這一劃分,將本文去特征處理得到的新數(shù)據(jù)協(xié)方差矩陣分為兩類進(jìn)行討論,對應(yīng)每一個新的數(shù)據(jù)協(xié)方差矩陣構(gòu)造了一個函數(shù),同理得到了兩組不同性質(zhì)的函數(shù),有這兩組函數(shù)構(gòu)造判據(jù)函數(shù),進(jìn)而對信源數(shù)進(jìn)行估計。本節(jié)最后總結(jié)給出了算法的具體步驟,對原理進(jìn)行分析,進(jìn)而將算法擴展到色噪聲情況下。

2.1 去特征處理

定義1 矩陣B減去其特征值與單位陣的乘積稱為這個特征值對矩陣B進(jìn)行了去特征處理。

用數(shù)據(jù)協(xié)方差矩陣RX的M個特征值分別對其進(jìn)行去特征處理,可以得到M個新的數(shù)據(jù)協(xié)方差矩陣R'Xi:

當(dāng)噪聲為空間平穩(wěn)高斯白噪聲時,噪聲協(xié)方差矩陣為

式中,σ2為噪聲功率,將式(9)代入式(7),則數(shù)據(jù)協(xié)方差矩陣為

由式(7)可得下面分兩種情況進(jìn)行討論。

2.1.1 M≥i≥M-D

當(dāng)M≥i≥M-D時,即λi為噪聲子空間對應(yīng)的特征向量時,σ2=λi,則

可以證明,R'Xi與RX的特征向量相同,即U也為R'Xi的特征向量矩陣,即

求R'Xi在U中的投影,有

定義函數(shù):

式中,i為在構(gòu)成新的協(xié)方差矩陣時所減去的特征值對應(yīng)的編號。由以上的分析可知,當(dāng)(1≤k≤D)時,fi(k)的值較大,當(dāng)D+1≤k≤M時,fi(k)的值較小,接近于0。

2.1.2 1≤i≤D

當(dāng)1≤i≤D時,即λi為信號子空間對應(yīng)的特征值時,σ2?λi,根據(jù)第2.1.1節(jié)的推導(dǎo)過程,同理可得到

式中,uk(k=1,2,…,M)為單位矩陣對應(yīng)第k列向量。則

與式(13)類似,定義函數(shù):當(dāng)1≤k≤D時,ek為信號子空間對應(yīng)的特征向量,則(λi-σ2)I}為一個接近于0的值。因此,fi(k)的值為一個接近于0的數(shù)。

當(dāng)D+1≤k≤M時,ek為噪聲子空間對應(yīng)的特征向量,則

當(dāng)1≤i≤D時,fi(k)對應(yīng)信號子空間的特征值,與fM(k)的取值規(guī)律相反,則g(i)取值較大。

當(dāng)M≥i≥M-D時,fi(k)對應(yīng)信號子空間的特征值,與fM(k)的取值規(guī)律相同,則g(i)取值較小。

據(jù)此,可利用GDE判別準(zhǔn)則對g(i)取值大小進(jìn)行區(qū)分,進(jìn)而實現(xiàn)信源數(shù)的估計,只不過這里需要對GDE判別準(zhǔn)則稍作修正即可。即是一個絕對值較大(接近于λi-σ2)的負(fù)數(shù),取模后為一個較大的數(shù)值,即fi(k)為一個較大值。

由以上分析可知,當(dāng)1≤k≤D時,fi(k)的值較小;當(dāng)D+1≤k≤M時,fi(k)的值較大。

對比M≥i≥M-D和1≤i≤D兩種情況可知,在兩種情況下,f(k)的取值規(guī)律剛好相反,本文即根據(jù)這一區(qū)別來判斷信源數(shù)目。

2.2 由fi(k)構(gòu)造統(tǒng)計量——判決函數(shù)

設(shè)構(gòu)造R'Xi時減去最小特征值所對應(yīng)的函數(shù)為fM(k),定義判決函數(shù):

k=1

其中,α(L,SNR)是關(guān)于快拍數(shù)L及信噪比SNR有關(guān)的調(diào)整因子,通過計算機仿真統(tǒng)計可知,其取值范圍為0到1。值得注意的是,快拍數(shù)L對調(diào)整因子影響權(quán)重較大,當(dāng)快拍數(shù)逐漸趨于無窮大時,調(diào)整因子的數(shù)值逐漸趨于1。當(dāng)i在1～M中取值,所有使Δ(i)＞0的g(i)個數(shù)即為所求信源數(shù)目。

2.3 本文算法描述

總結(jié)本文算法步驟如下:

步驟1 對數(shù)據(jù)協(xié)方差矩陣RX進(jìn)行特征值分解,求得對應(yīng)特征值從大到小排列的特征空間U;

步驟2 根據(jù)式(7)構(gòu)造M個新的數(shù)據(jù)協(xié)方差矩陣R'Xi;

步驟3 分別求R'Xi在特征空間U上的投影矩陣;

步驟4 根據(jù)式(13)得到M個函數(shù)fi(k)(i=1,2,…,M);

步驟5 根據(jù)式(19)求得g(i),(i=1,2,…,M-1);

步驟6 設(shè)置調(diào)整因子α(L,SNR),根據(jù)式(20)進(jìn)行信源數(shù)估計。

2.4 算法原理分析

<1)，且各件產(chǎn)品是否為不合格品相互獨立．

構(gòu)造R'Xi時減去的λi為對應(yīng)噪聲子空間的特征值時,相當(dāng)于對原數(shù)據(jù)協(xié)方差矩陣進(jìn)行了去噪聲處理,此時得到的新數(shù)據(jù)協(xié)方差矩陣中信號的貢獻(xiàn)變大,也就是說新協(xié)方差矩陣在信號子空間上的投影值較大,而在噪聲子空間上的投影值較小;而構(gòu)造R'Xi時減去的λi為對應(yīng)信號子空間的特征值時,相當(dāng)于對原數(shù)據(jù)協(xié)方差矩陣進(jìn)行了去信號處理,此時得到的新數(shù)據(jù)協(xié)方差矩陣中信號的貢獻(xiàn)相對變小,新協(xié)方差矩陣在信號子空間上的投影值較小,而在噪聲子空間上的投影值較大。由以上分析可知,根據(jù)新協(xié)方差矩陣在信號子空間和噪聲子空間上的投影值的不同,可將經(jīng)過去特征處理得到的新協(xié)方差矩陣分為對應(yīng)于信號和噪聲的兩類,而對應(yīng)于信號的新協(xié)方差矩陣的個數(shù)即為信源數(shù)。本文算法就是利用這一思想增大了協(xié)方差矩陣在特征向量上的投影值的差異,構(gòu)造統(tǒng)計量fi(k)及判據(jù)函數(shù)g(i),從而判斷出信源數(shù)。

2.5 算法復(fù)雜度分析

本文算法的計算量主要集中在協(xié)方差矩陣的構(gòu)造、特征值分解及投影運算上。由M×1維的向量構(gòu)造協(xié)方差矩陣的計算量為C1=M2,構(gòu)造的協(xié)方差矩陣的維數(shù)為M× M,對該M×M維矩陣進(jìn)行特征分解,采用文獻(xiàn)[20]的快速子空間分解技術(shù)可將計算量降低至C2=M(L+2)2,求R'Xi在特征空間U上的投影矩陣的計算量為C31=M3,算法共有M次該投影運算,則求投影矩陣運算的計算量為C3= MC31=M4,則本文算法總計算量為C=C1+C2+C3= M2+M(L+2)2+M4,而經(jīng)典AIC算法和MDL算法僅需要進(jìn)行協(xié)方差矩陣的構(gòu)造和協(xié)方差矩陣特征分解,本文算法的投影矩陣運算引入了一定的計算量,導(dǎo)致最終算法的復(fù)雜度增加。但是復(fù)雜度的增加換來的是算法在兩種噪聲背景下及白噪聲和空域相關(guān)色噪聲并存背景下的估計性能的提升。

表1將本文算法與信息論類算法(AIC算法和MDL算法)的算法復(fù)雜度進(jìn)行了對比。

表1 算法復(fù)雜度對比

3 仿真實驗及結(jié)果分析

為了驗證本文提出算法的性能,利用Matlab進(jìn)行了計算機仿真實驗。實驗分別驗證了本文算法在白噪聲背景下及空域相關(guān)色噪聲背景下的測向性能,為了驗證本文算法使用于白噪聲和空域相關(guān)色噪聲成分同時存在的情況,還對算法在這兩種背景條件下的測向性能進(jìn)行了對比,仿真條件如下:3個相互獨立的遠(yuǎn)場窄帶信號入射到8陣元均勻線陣,陣元間距為半波長,3個信源的方向角為[15°30°40°]。在白噪聲情況下,將本文算法與AIC準(zhǔn)則、MDL準(zhǔn)則、m AIC準(zhǔn)則、m MDL準(zhǔn)則及GDE準(zhǔn)則進(jìn)行對比。由于在空域相關(guān)色噪聲背景下AIC準(zhǔn)則和MDL準(zhǔn)則失效,故在色噪聲背景下僅對本文算法、mAIC準(zhǔn)則、m MDL準(zhǔn)則及GDE準(zhǔn)則進(jìn)行對比。每個實驗均進(jìn)行200次Monte Carlo仿真。算法的評價標(biāo)準(zhǔn)為成功檢測概率,即200次Monte Carlo仿真實驗中能夠成功檢測出信源數(shù)的次數(shù)的百分比。

實驗1 白噪聲背景下成功檢測概率隨信噪比的變化

設(shè)置快拍數(shù)為200,信噪比從-14 dB變化到10 dB,在白噪聲背景下進(jìn)行仿真,得到正確檢測概率隨信噪比的變化。

白噪聲背景下成功檢測概率隨信噪比的變化仿真結(jié)果如圖1所示?？梢钥闯?6種算法的成功檢測概率均隨SNR取值的變大而變大。其中本文算法、AIC準(zhǔn)則及MDL準(zhǔn)則在信噪比大于-4dB時成功檢測概率為100%,而mAIC準(zhǔn)則、m MDL準(zhǔn)則和GDE準(zhǔn)則在信噪比分別取-2 dB、0 dB和6 dB時成功檢測概率達(dá)到100%?？梢?本文算法在低信噪比下估計信源數(shù)的成功概率與信息論準(zhǔn)則算法相當(dāng)而明顯優(yōu)于mAIC準(zhǔn)則、m MDL準(zhǔn)則及GDE準(zhǔn)則。

圖1 白噪聲下正確信源數(shù)估計概率隨SNR變化

實驗2 空域相關(guān)色噪聲背景下成功檢測概率隨信噪比的變化

設(shè)置快拍數(shù)為200,信噪比以2 d B為步進(jìn)值從-10 dB變化到10 dB。色噪聲采用文獻(xiàn)[10]給出的空域相關(guān)色噪聲模型。令該色噪聲表達(dá)式中空間相關(guān)系數(shù)ρ=0.5。將本文算法與mAIC準(zhǔn)則、mMDL準(zhǔn)則及GDE準(zhǔn)則進(jìn)行對比,得到空域相關(guān)色噪聲背景下成功檢測概率隨信噪比的變化。

空域相關(guān)色噪聲背景下成功檢測概率隨信噪比的變化如圖2所示?？梢钥闯?隨著信噪比的增加,以上4種算法的成功檢測概率均變大。其中本文算法的成功檢測概率要高于m AIC準(zhǔn)則、m MDL準(zhǔn)則及GDE準(zhǔn)則對應(yīng)的成功檢測概率。本文算法在信噪比大于-3 dB時成功檢測概率達(dá)到100%,而m AIC準(zhǔn)則、m MDL準(zhǔn)則的信噪比均需達(dá)到0 dB、GDE準(zhǔn)則的信噪比需達(dá)到4 dB時,成功檢測概率才能達(dá)到100%?？梢?在色噪聲背景下本文算法的性能要優(yōu)于其他3種方法。

圖2 空域相關(guān)色噪聲下正確信源數(shù)估計概率隨SNR變化

實驗3 白噪聲背景下成功檢測概率隨快拍數(shù)的變化

設(shè)置信噪比為5 dB,快拍數(shù)以10為步進(jìn)值從10變化到230。設(shè)置空間噪聲為白噪聲,對本文算法、AIC準(zhǔn)則、MDL準(zhǔn)則、改進(jìn)mAIC準(zhǔn)則、改進(jìn)m MDL準(zhǔn)則及GDE方法進(jìn)行仿真,得到正確檢測概率隨快拍數(shù)的變化。

白噪聲背景下成功檢測概率隨快拍數(shù)的變化仿真結(jié)果如圖3所示,可以看出,6種方法的成功檢測概率均隨快拍數(shù)的增多而變大,本文算法在相同快拍數(shù)的情況下,較其他算法的成功檢測概率大,但其檢測性能較AIC準(zhǔn)則和MDL準(zhǔn)則差別不大。本文算法、AIC準(zhǔn)則及MDL準(zhǔn)則當(dāng)快拍數(shù)增加到30時,成功檢測概率均達(dá)到100%,而m AIC準(zhǔn)則、mMDL準(zhǔn)則及GDE方法的檢測性能則要明顯差于以上3種方法?？梢姳疚乃惴ㄔ诎自肼暠尘跋碌臋z測性能并沒有下降,甚至在較低快拍數(shù)的情況下的性能要略優(yōu)于AIC準(zhǔn)則和MDL準(zhǔn)則。結(jié)合實驗1的仿真結(jié)果可知,本文算法適用于白噪聲背景條件下,且估計性能與傳統(tǒng)的信息論方法相當(dāng)。

圖3 白噪聲下正確信源數(shù)估計概率隨快拍數(shù)變化

實驗4 空域相關(guān)色噪聲背景下成功檢測概率隨快拍數(shù)的變化

與實驗2的空域相關(guān)色噪聲背景相同,取信噪比為10 dB,快拍數(shù)從1依次變化到240。對4種方法的檢測成功概率進(jìn)行仿真。

空域相關(guān)色噪聲背景下成功檢測概率隨快拍數(shù)的變化仿真結(jié)果如圖4所示?？梢钥闯?4種算法的成功檢測概率均隨著快拍數(shù)的增加而變大。在低快拍數(shù)的情況下,本文算法的成功檢測概率明顯高于其他兩種方法,在快拍數(shù)大于40時成功檢測概率達(dá)到100%,而其他3種算法則需要快拍數(shù)達(dá)到80、100甚至更大時,成功檢測概率才能達(dá)到100%?？梢?在空域相關(guān)色噪聲背景下,本文算法的性能優(yōu)于其他3種方法。結(jié)合實驗2的仿真結(jié)果可知,本文算法適用于色噪聲環(huán)境,且在一定條件下,性能優(yōu)于其他3種算法。

圖4 空域相關(guān)色噪聲下正確信源數(shù)估計概率隨快拍數(shù)變化

實驗5 在白噪聲和空域相關(guān)色噪聲兩種噪聲背景下的檢測性能對比

為了驗證本文算法可以適用于包含白噪聲成分的色噪聲環(huán)境,本文通過計算機仿真對白噪聲、空域相關(guān)色噪聲及包含有白噪聲成分的空域相關(guān)色噪聲3種噪聲背景下的檢測性能進(jìn)行了對比仿真,對比情況如圖5和圖6所示。

本文算法在3種噪聲背景下的成功檢測概率隨SNR變化的對比情況如圖5所示?？梢钥闯?在白噪聲背景下,本文算法的成功檢測概率達(dá)到100%時的信噪比為-4 dB,色噪聲情況下則為-2 dB,二者僅相差2 dB,而在包含有白噪聲成分的色噪聲背景情況下,本文算法的成功檢測概率介于單純白噪聲和單純空域相關(guān)色噪聲背景之間,接近于白噪聲背景。

圖5 正確信源數(shù)估計概率隨SNR的變化

成功檢測概率在3種噪聲背景下隨快拍數(shù)變化對比情況如圖6所示?？梢钥闯?在快拍數(shù)大于50時,3種噪聲背景下的成功檢測概率均達(dá)到100%。

通過以上5個實驗結(jié)果可以看出,本文算法提高了色噪聲背景下的估計性能,在適用于空間色噪聲的同時,在白噪聲背景下的檢測性能也沒有明顯的下降。因此,本文算法可用于包含有白噪聲成分的空域相關(guān)色噪聲環(huán)境下進(jìn)行信源數(shù)估計。

圖6 正確信源數(shù)估計概率隨快拍數(shù)的變化

4 結(jié) 論

實際測向系統(tǒng)中的背景噪聲均比較復(fù)雜,很難在理論上及工程上將實際環(huán)境中的噪聲單純的分為白噪聲或色噪聲,而往往是兩者并存的情況。本文針對這一現(xiàn)象提出一種基于去特征處理的信源數(shù)估計方法,仿真實驗結(jié)果驗證了該算法不僅提高了空域相關(guān)色噪聲背景下的估計性能,同樣在白噪聲背景下也具有較好的性能。通過對白噪聲、空域相關(guān)色噪聲及包含有白噪聲成分的空域相關(guān)色噪聲背景下算法的性能進(jìn)行仿真對比,驗證了本文算法適用于包含有白噪聲成分的空域相關(guān)色噪聲環(huán)境,并且在低信噪比和小快拍數(shù)情況下且具有較好的估計性能。但是本文算法的計算復(fù)雜度要略高于傳統(tǒng)算法,并且在利用判決函數(shù)進(jìn)行信源數(shù)估計時需要人工設(shè)置閾值參數(shù),因此如何適當(dāng)?shù)慕档陀嬎銖?fù)雜度及如何通過自適應(yīng)算法來設(shè)置閾值參數(shù)將是下一步研究的重點。

[1]Max M,Tailath T.Detection of signals by information theoretic criteria[J].IEEE Trans.on ASSP,1985,33(2):387-392.

[2]Ye Z F,Xiang L,Xu X.Improvement of source number estimation based on information theoretic criteria[J].Chinese Journal of Radio Science,2007,22(4):593-598.(葉中付,向利,徐旭.基于信息論準(zhǔn)則的信源個數(shù)估計算法改進(jìn)[J].電波科學(xué)學(xué)報, 2007,22(4):593 598.)

[3]Dong S M,Liang G L.A modified Gerschgorin disks estimation method for source number estimation[J].Journal of Harbin Engineering University,2013(4):440 444.(董姝敏,梁國龍.改進(jìn)的蓋爾圓源數(shù)目估計方法[J].哈爾濱工程大學(xué)學(xué)報,2013 (4):440-444.)

[4]Cozzens J H,Sousa M J.Source enumeration in a correlate signed environment[J].IEEE Trans.on Signal Processing,1994,42 (2):304-317.

[5]Huang L,Wu S,Li X.Reduced-rank MDL method for source enumeration in high-resolution array processing[J].IEEE Trans. on Signal Processing,2007,55(12):5658-5667.

[6]Chen W,Reilly J P.Detection of the number of signals in noise with banded covariance matrices[J].IEEE Trans.on Signal Processing,1992,42(5):377-380.

[7]Wong K M,Wu Q,Stoica P.Generalized correlation decomposition applied to array processing in unknown noise environments[M].2nd ed.Haykin:Prentice Hall,1995.

[8]Zhang J,Liao G S,Wang J.Performance improvement of source number detection using diagonal loading[J].Chinese Journal of Electronics,2004,32(12):2094-2097.(張杰,廖桂生,王玨.對角加載對信號源數(shù)檢測性能的改善[J].電子學(xué)報,2004,32 (12):2094-2097.)

[9]Xie J L,Si X C.Determining the number of sources based on diagonal loading to the covariance matrix[J].Systems Engineering and Electronics,2008,29(5):596-600.(謝紀(jì)嶺,司錫才.基于協(xié)方差矩陣對角加載的信源數(shù)估計方法[J].系統(tǒng)工程與電子技術(shù),2008,29(5):596 600.)

[10]Cedervall SP.Detection tests for array processing in unknown correlated noise fields[J].IEEE Trans.on Signal Processing, 1997,49(5):2351 2362.

[11]Xiao M L,Wei P,Tai H M.Estimation of the number of sources based on hypothesis testing[J].Journal of Communications and Networks,2012,14(5):481 999.

[13]Candan C,Koc S.Beamspace approach for detection of the number of coherent sources[C]∥Proc.of the IEEE Radar Conference,2012.

[14]Cheng W,Lee S,Zhang Z S,et al.Independent component analysis based source number estimation and its comparison for mechanical systems[J].Journal of Sound and Vibration, 2012,331(23):5153-5167.

[15]Ye H X,Hu X P.A new mechanical vibration source number estimation method based on EMD-SVD-DE[C]∥Proc.of the International Conference on Information Science,Automation and Material System,2011:1379-1382.

[16]Wang P,Zhang G J,Zhang W D,et al.Source number estimation algorithm using Gerschgorin disks for acoustic vector sensor array[J].Journal of Computational Information Systems, 2012,8(7):2921-2929.

[17]Xu L,Zeng C,Liao G S,et al.DOA and source number estimation method for strong and weak signals based on eigen beamforming[J].Journal of Electronics and Information Technology,2011,33(2):321-325.

[18]Lu Z H,Zoubir A M.Flexible detection criterion for source enumeration in array processing[J].IEEE Trans.on Signal Processing,2013,61(6):1303-1314.

[19]Lu Z H,Zoubir A M.Generalized Bayesian information criterion for source enumeration in array processing[J].IEEE Trans.on Signal Processing,2013,61(6):1470-1480.

[20]Xu G,Kailath T.Fast subspace decomposition[J].IEEE Trans.on Signal Processing,1994,42(3):539-

New source number estimation method based on feature eliminated process

WU Na,SI Wei-jian,JIAO Shu-hong,WU Di
(School of Information and Communication Engineering,Harbin Engineering University,Harbin 150001,China)

In actual direction-finding systems,the impinging noise is a mixture of white noise and colored noise.The estimation performance of the traditional algorithm used for colored noise tends to decline when it is used under the background of white noise.Accordingly,a new source number estimation method based on the feature eliminated process is proposed in this study.The method performs the feature eliminated process on data covariance firstly.A group of new covariance matrixes according to the eigen-values are obtained.Using the projection matrixes of signal subspace and noise subspace to construct the criterion function,then get the number of signals.Simulation results show that the method proposed can effectively estimate the number of signals under the white noise and colored noise background.The algorithm can estimate the number of sources with high accuracy rate under the condition of small snapshot and low SNR without obvious computation increase.

source number estimation;colored noise;feature eliminated process;decision function

TN 911.7

A

10.3969/j.issn.1001-506X.2015.03.06

吳 娜(1986-),女,博士研究生,主要研究方向為寬帶系統(tǒng)的信號檢測、處理與識別。

E-mail:wunahmc@163.com司偉建(1971-),男,研究員,博士,主要研究方向為寬帶系統(tǒng)的信號檢測、處理與識別。

E-mail:swj0418@263.com焦淑紅(1966-),女,教授,博士,主要研究方向為無源定位寬帶信號處理。E-mail:jiaoshuhong@hrbeu.edu.cn吳 迪(1987-),女,博士研究生,主要研究方向為寬帶系統(tǒng)的信號檢測、處理與識別。

E-mail:wudi7212@126.com

網(wǎng)址:www.sys-ele.com

1001-506X(2015)03-0509-06

2013 12 02;

2014 06 10;網(wǎng)絡(luò)優(yōu)先出版日期:2014 09 26。

網(wǎng)絡(luò)優(yōu)先出版地址:http:∥w ww.cnki.net/kcms/detail/11.2422.TN.20140926.1557.016.html

中央高校基本科研業(yè)務(wù)費專項資金項目(HEUCF140803)資助課題