(長(zhǎng)沙理工大學(xué) 電氣與信息工程學(xué)院，長(zhǎng)沙 410114)

1 引 言

信源數(shù)的判定是陣列信號(hào)處理領(lǐng)域的關(guān)鍵問題[1-2]。準(zhǔn)確的信源數(shù)信息是其他陣列信號(hào)參數(shù)進(jìn)行高分辨估計(jì)的前提，若信源數(shù)未知或估計(jì)不準(zhǔn)，勢(shì)必導(dǎo)致后續(xù)相關(guān)算法的性能大幅下降，因此準(zhǔn)確的信源數(shù)估計(jì)具有一定的實(shí)際意義。

許多學(xué)者根據(jù)不同的準(zhǔn)則提出了不同的估計(jì)方法，其中最常見的一種是基于信息論準(zhǔn)則的估計(jì)算法，包括Akaike信息論(Akaike′s Information Crilertion,AIC[3])和最小長(zhǎng)度描述準(zhǔn)則(Mininum Description Length,MDL[4])。這些方法都是通過組合特征值分解、最大似然函數(shù)和一些罰函數(shù)來進(jìn)行檢測(cè)，優(yōu)點(diǎn)是計(jì)算簡(jiǎn)單，不需認(rèn)為設(shè)置門限，但當(dāng)接收信號(hào)包含色噪聲時(shí)易失效。文獻(xiàn)[5]利用接收信號(hào)的延時(shí)相關(guān)信息重構(gòu)協(xié)方差矩陣并結(jié)合對(duì)角加載技術(shù)對(duì)信源數(shù)進(jìn)行估計(jì)，有效克服了上述缺點(diǎn)。另一種是通過比較蓋氏圓半徑大小來估計(jì)信源數(shù)，稱為蓋氏圓盤估計(jì)法(Gerschgorin′s Disk Estimation，GDE[6])，在色噪聲環(huán)境下能有效實(shí)現(xiàn)信源數(shù)估計(jì)，缺點(diǎn)是低信噪比時(shí)檢測(cè)效果較差。文獻(xiàn)[7]直接對(duì)寬帶信號(hào)陣列進(jìn)行快速傅里葉變換，得到了一種在色噪聲環(huán)境下估計(jì)寬帶相干信源數(shù)的方法。

由于現(xiàn)有信源數(shù)估計(jì)算法都是針對(duì)陣列信號(hào)，不能直接用于單通道接收模型，結(jié)合間隔抽樣和刀切法，提出了一種改進(jìn)算法。該算法首先通過間隔抽樣實(shí)現(xiàn)單通道接收信號(hào)空間的矢量化，接著采用刀切法對(duì)此矢量化空間重構(gòu)協(xié)方差矩陣以便高效利用現(xiàn)有數(shù)據(jù)和抑制噪聲的影響。實(shí)驗(yàn)結(jié)果表明，與常規(guī)算法相比，該算法在信噪比較低和采樣數(shù)較少的情況下具有較好的檢測(cè)性能。

2 信號(hào)模型描述

2.1 信號(hào)模型

對(duì)一個(gè)陣列接收信號(hào)，假設(shè)有p個(gè)通道同時(shí)接收信號(hào)，接收信號(hào)中包含q個(gè)源信號(hào)，那么信號(hào)接收模型可表示為[8-10]

x(t)=As(t)+n(t) 。

(1)

式中：x(t)∈p×l為觀測(cè)信號(hào)，A∈p×q為混合矩陣，s(t)∈q×l為源信號(hào)，n(t)∈p×l為噪聲。

由于噪聲與源信號(hào)相互獨(dú)立，且是均值為0的隨機(jī)過程，則觀測(cè)信號(hào)的協(xié)方差矩陣表示為

(2)

(3)

式中：λ1>λ2>…>λq>λq+1=…=λp是R按照降序排列出來的特征值，ui是λi對(duì)應(yīng)的特征向量。

2.2 常規(guī)的信源數(shù)估計(jì)算法

常規(guī)的信源數(shù)估計(jì)算法包括基于AIC、MDL和GDE準(zhǔn)則的信源數(shù)估計(jì)算法。AIC和MDL準(zhǔn)則由似然函數(shù)和罰函數(shù)兩部分組成：

AIC(k)=2N(M-k)lbf(k)+2k(2M-k) ，

(4)

MDL(k)=N(M-k)lbf(k)+0.5(2M-k)lbN。

(5)

GDE準(zhǔn)則本質(zhì)是通過比較特征向量來估計(jì)信源數(shù)。該方法首先對(duì)采樣協(xié)方差矩陣Rj進(jìn)行酉變換：

RM=THRjT。

(6)

經(jīng)變換后可得

(7)

矩陣RM中半徑非零的蓋氏圓盤的個(gè)數(shù)即為信源數(shù)的估計(jì)值。在實(shí)際工程應(yīng)用中，采用如下函數(shù)來估計(jì)信源數(shù)：

(8)

式中：D(N)為與快拍數(shù)N有關(guān)的調(diào)整因子，取值在0～1之間。分別計(jì)算GDE(k)在1≤k≤p之間的值，當(dāng)GDE(k)的值首次小于零時(shí)停止計(jì)算，信源數(shù)的估計(jì)值取k-1。

在色噪聲環(huán)境下，基于AIC和MDL的信源數(shù)估計(jì)算法的檢測(cè)性能急劇下降，甚至失效。本文引入對(duì)角加載技術(shù)[11-12]，對(duì)噪聲特征值進(jìn)行某種修正，即

αi=λi+λD。

(9)

式中：λi(i=1,2,…,p)為采樣協(xié)方差矩陣原有的特征值，λD為加載量，αi為加載后的特征值。λD的選取是對(duì)角技術(shù)的關(guān)鍵。文獻(xiàn)[13]提出了一種λD選取的經(jīng)驗(yàn)值：

(10)

2.3 刀切法的基本原理

刀切法(Jackknife)的目的是降低估計(jì)的偏差，其基本思想是從樣本集中一次刪去若干個(gè)數(shù)據(jù)，對(duì)一系列剩余樣本計(jì)算估計(jì)值，在估計(jì)值中得到最優(yōu)估計(jì)[14]。假設(shè)樣本X=(x1,x2,…,xn)獨(dú)立同分布，θ=θ(x1,x2,…,xn)為未知參數(shù)θ的一個(gè)有偏估計(jì)，則θ的偏差為

(11)

(12)

θ的刀切法估計(jì)的偏差為[15]

(13)

結(jié)合式(11)和式(13)可知，估計(jì)偏差由o(n-1)變成o(n-2)，復(fù)雜度降低，準(zhǔn)確性得到提高。所以，刀切法估計(jì)可以近似為無偏估計(jì)。

3 采用刀切法的單通道信源數(shù)估計(jì)算法

在單通道接收信號(hào)時(shí)，只有一個(gè)傳感器可以利用，觀測(cè)數(shù)據(jù)只能獲得一個(gè)，假設(shè)有兩個(gè)源信號(hào)被單個(gè)傳感器接收，則式(1)可表示為

x(t)=a1(t)s1(t)+a2(t)s2(t)+n(t)。

(14)

式中：a1、a2為比例因子；s1(t)、s2(t)為源信號(hào)；n(t)為噪聲。

假設(shè)輸入信號(hào)x(t)滿足采樣條件，忽略噪聲的影響，對(duì)其進(jìn)行周期為Δ的采樣，則輸入信號(hào)x(t)的離散形式為

x(kΔ)=a1s1(kΔ)+a2s2(kΔ) 。

(15)

對(duì)離散信號(hào)x(kΔ)進(jìn)行周期為T的采樣，令t=kT+iΔ，i=1,2,…,M則新序列x(i)(kT)為

x(i)(kT)=x(kT+iΔ) 。

(16)

式中：x(kT+iΔ)=a1s1(kT+iΔ)+a2s2(kT+iΔ),i=1,2,…,M。

由此可知，如果序列x(kΔ)按周期T進(jìn)行抽樣，且新序列的采樣周期T是Δ的整數(shù)倍，即T/Δ=M(M∈)，則通過M倍抽樣，可將單通道信號(hào)分解為M通道信號(hào)，用公式表示為

(17)

將單通道接收信號(hào)升維后，可運(yùn)用現(xiàn)有的陣列信號(hào)處理算法對(duì)信源數(shù)進(jìn)行估計(jì)。為了充分利用接收到的數(shù)據(jù)，實(shí)現(xiàn)更加精準(zhǔn)的估計(jì)，采用刀切法對(duì)信源數(shù)目進(jìn)行估計(jì)。在刀切法中，對(duì)一個(gè)包含N個(gè)元素的集合X，X=[X(1),X(2),…,X(N)]從其中隨機(jī)選取M個(gè)元素組成一個(gè)新的集合XM:

(18)

式中：XM?X，x(t)∈X,t=1,2,…,M且M=rN,0.5≤r<1。重復(fù)上述過程k次，得到協(xié)方差矩陣Rj：

(19)

對(duì)得到的協(xié)方差矩陣Rj進(jìn)行特征值分解：

(20)

(21)

下面從理論上驗(yàn)證采用刀切法重構(gòu)的協(xié)方差矩陣Rj能有效抑制高斯噪聲的影響。

證明：當(dāng)整個(gè)樣本的信源數(shù)檢測(cè)概率p大于50%時(shí)，pj表示采用刀切法后信源數(shù)的檢測(cè)概率，則pj≥p。分兩種情況驗(yàn)證。

當(dāng)Z=2n時(shí)，即迭代次數(shù)為偶數(shù)時(shí)，采用刀切法后的檢測(cè)概率pj為

(22)

式中：p為正確檢測(cè)的概率，q為錯(cuò)誤檢測(cè)的概率。同理可得當(dāng)Z=2n+1時(shí)，pj≥p。所以不管迭代次數(shù)Z取何值，都有pj≥p，得證。

基于以上分析可知，本文提出的算法能夠適用于單通道模型，并能有效抑制高斯白噪聲和高斯色噪聲的影響。圖1給出了本文算法流程圖。

圖1 本文算法的流程框圖Fig.1 Flow chart of the proposed algorithm

利用刀切法提高信源數(shù)的檢測(cè)概率是以增加算法復(fù)雜度為代價(jià)的。設(shè)定陣元數(shù)M=8，快拍數(shù)T=1 000，迭代次數(shù)Z=20，采樣比例r=0.625。由于在多次采用刀切法和酉變換的過程中，計(jì)算相互獨(dú)立，本文算法可通過并行計(jì)算方式處理，此時(shí)常規(guī)算法的運(yùn)行時(shí)間是0.002 9 s，本文算法的運(yùn)行時(shí)間是0.003 2 s。常規(guī)算法GDE和本文算法的算法復(fù)雜度如表1所示。

表1 常規(guī)GDE算法與本文算法的復(fù)雜度分析Tab.1 Complexity analysis of GDE algorithm and improved algorithm

從表1中可知，雖然本文算法需要多次循環(huán)迭代處理，但都可以通過增加空間負(fù)荷來減少計(jì)算時(shí)間。對(duì)比常規(guī)算法和本文算法，除在取最大頻次過程中的時(shí)間復(fù)雜度有所不同外，其余基本一樣。

4 實(shí)驗(yàn)仿真與分析

為了最大程度模擬真實(shí)情況下的信源數(shù)檢測(cè)性能，分別在白噪聲和色噪聲環(huán)境下進(jìn)行仿真分析。由于常規(guī)算法的檢測(cè)性能會(huì)隨著快拍數(shù)的減少而下降，因此本文取快拍數(shù)T=2 000。

實(shí)驗(yàn)1：假設(shè)單通道接收信號(hào)中包含3個(gè)信號(hào)源，經(jīng)AD采樣后取快拍數(shù)T=2 000，采用間隔抽樣法構(gòu)建12個(gè)虛擬通道，使得單通道接收信號(hào)轉(zhuǎn)換成陣列信號(hào)的形式。設(shè)定信噪比變化范圍為-15～15 dB，以2 dB遞增，分別在白噪聲和色噪聲環(huán)境下進(jìn)行200次蒙特卡洛實(shí)驗(yàn)，基于AIC準(zhǔn)則和MDL準(zhǔn)則方法、GDE方法和本文算法的檢測(cè)概率如圖2所示。

(a)白噪聲環(huán)境下

(b)色噪聲環(huán)境下圖2 四種算法檢測(cè)性能與信噪比的關(guān)系Fig.2 The relationship between performance and SNR of four algorithms

從圖2中可以看出，當(dāng)信噪比低于-7 dB時(shí)，三種常規(guī)算法的性能急劇下降，而本文算法的檢測(cè)概率仍然較高。在色噪聲環(huán)境下時(shí)，基于信息論的AIC準(zhǔn)則和MDL準(zhǔn)則對(duì)信源數(shù)失去了準(zhǔn)確估計(jì)；當(dāng)信噪比為-7 dB時(shí)，GDE方法的檢測(cè)性能為0.7，而本文算法的檢測(cè)概率已收斂于1。因此當(dāng)通道噪聲為白或色噪聲時(shí)，相比其他三種常規(guī)算法，本文算法具有良好的檢測(cè)性能。

王超，范雪冰，佟首峰,等.空間光到少模光纖的耦合效率及影響因素[J].光子學(xué)報(bào),2018,47(11)：1206001

實(shí)驗(yàn)2：假設(shè)一單通道接收信號(hào)中包含3個(gè)復(fù)指數(shù)信號(hào)，設(shè)定信噪比為-5 dB，采樣快拍數(shù)T=2 000，圖3表示了本文算法的檢測(cè)概率與嵌入維數(shù)的關(guān)系。從圖3中可以看出，本文算法的檢測(cè)概率隨嵌入維數(shù)的增加呈先上升后下降的趨勢(shì)。當(dāng)嵌入維數(shù)小于7時(shí)，檢測(cè)概率低于0.5，因?yàn)殚g隔抽樣后的協(xié)方差矩陣不能完全包括原單通道接收信號(hào)的全部信息；當(dāng)嵌入維數(shù)高于12時(shí)，檢測(cè)概率呈下降趨勢(shì)，因?yàn)榍度刖S數(shù)過多，各分量之間會(huì)相互干擾。

圖3 嵌入維數(shù)的選擇對(duì)算法性能的影響Fig.3 The influence of embedding dimension on algorithm performance

實(shí)驗(yàn)3：假設(shè)一單通道接收信號(hào)中包含3個(gè)復(fù)指數(shù)信號(hào)，在高斯白噪聲環(huán)境下設(shè)定信噪比為10 dB，采樣快拍數(shù)由100變換到1 000；在色噪聲環(huán)境下設(shè)定信噪比為-7 dB時(shí)，采樣快拍數(shù)由1 000以200遞增到3 000，每種實(shí)驗(yàn)進(jìn)行200次蒙特卡洛實(shí)驗(yàn)。圖4為四種算法在不同快拍數(shù)下的正確檢測(cè)概率。

(a)SNR=10 dB，白噪聲環(huán)境下

(b)SNR=-7 dB,色噪聲環(huán)境下圖4 四種算法估計(jì)性能與采樣點(diǎn)數(shù)的關(guān)系Fig.4 The relationship between the estimation performance of four algorithms and the number of sampling points

由圖4的仿真結(jié)果可知，當(dāng)通道噪聲為高斯白噪聲時(shí)，四種算法的檢測(cè)概率都隨快拍數(shù)的增大而提高。當(dāng)快拍數(shù)低于200時(shí)，AIC方法的檢測(cè)性能最好，本文算法略低于AIC算法但明顯高于MDL和GDE算法；當(dāng)快拍數(shù)高于300時(shí)，本文算法的檢測(cè)概率優(yōu)先收斂1。當(dāng)通道噪聲為高斯色噪聲時(shí)，本文算法的檢測(cè)概率明顯優(yōu)于其他三種算法，保持了良好的檢測(cè)性能。

實(shí)驗(yàn)4：為了獲得更加精確的估計(jì)結(jié)果，選擇200組實(shí)測(cè)數(shù)據(jù)分別在信號(hào)源為2個(gè)和3個(gè)的情況下采用四種算法對(duì)信源數(shù)進(jìn)行估計(jì)。表2和表3的數(shù)字表示四種算法估計(jì)不同信源數(shù)的次數(shù)。

表2 2個(gè)信源的估計(jì)結(jié)果Tab.2 Estimation results of two signals

表3 3個(gè)信源的估計(jì)結(jié)果Tab.3 Estimation results of three signals

從表2和表3可知，當(dāng)信源數(shù)個(gè)數(shù)為2時(shí)，AIC準(zhǔn)則、MDL準(zhǔn)則、GDE準(zhǔn)則和本文算法的估計(jì)正確概率分別為81%、85%、82.5%和92.5%；當(dāng)信源數(shù)個(gè)數(shù)為3時(shí)，四種算法的估計(jì)正確概率分別為65%、80%、58%和87%。由此可知，在不同信源數(shù)下，本文算法檢測(cè)的正確概率優(yōu)于其他三種算法。

5 結(jié)束語

本文針對(duì)現(xiàn)有信源數(shù)估計(jì)存在的問題，提出了一種改進(jìn)算法。該算法首先采用間隔抽樣法有效地解決了現(xiàn)有算法無法直接用于單通道信源數(shù)估計(jì)的問題，然后利用刀切法重構(gòu)多個(gè)協(xié)方差矩陣，并多次進(jìn)行酉變換，充分利用有限信息，避免了因舍去相關(guān)的信號(hào)分量造成誤差。仿真實(shí)驗(yàn)表明，與原有算法相比，在低信噪比和低快拍數(shù)下，本文提出的改進(jìn)算法在白噪聲和色噪聲環(huán)境中的檢測(cè)性能得到了較大改善。本文算法的不足在于只研究了高斯噪聲下的單通道信源數(shù)估計(jì)問題，因此還需要對(duì)非高斯噪聲的情況進(jìn)一步探索。

：

[1] PAZOS S,HURTADO M,MURAVCHIK C.On sparse methods for array signal processing in the presence of interference[J].IEEE Antennas & Wireless Propagation Letters,2015,14(5):1165-1168.

[2] LIU W.Antenna array signal processing for quaternion-valued wireless communication systems[C]//Proceedings of 2014 IEEE Benjamin Franklin Symposium on Microwave and Antenna Sub-systems for Radar,Telecommunications,and Biomedical Applications(BenMAS).Philadelphia,PA,USA:IEEE,2014:1-3.

[3] CHAKRABARTI A,GHOSH J K.AIC,BIC and recent advances in model selection[J].Philosophy of Statistics,2011(2):583-605.

[4] NANNEN V.A short introduction to model selection,kolmogorov complexity and minimum description length(MDL)[J].Computer Science,2010,46(2):446-464.

[5] 程晨,韓惠蓮,司偉建.基于延時(shí)處理的對(duì)角加載信源數(shù)估計(jì)方法[J].云南大學(xué)學(xué)報(bào)(自然科學(xué)版),2014(6):823-829.

CHENG Chen，HAN Huilian，SI Weijian.The sources number estimation using delay correlation processing based on diagonal Loading[J].Journal of Yunnan University(Natural Sciences Edition),2014(6):823-829.(in Chinese)

[6] LU Z,GAO M,JIANG H.Source number estimation for minimum redundancy arrays with Gerschgorin disk estimator[C]//Proceedings of 2012 IEEE 11th International Conference on Signal Processing(ICSP).Beijing:IEEE,2012:311-314.

[7] 周圍,王新青,閆霄翔,等.空間平滑和蓋氏半徑變換的寬帶相干信源數(shù)估計(jì)[J].重慶郵電大學(xué)學(xué)報(bào)(自然科學(xué)版),2014(3):346-351.

ZHOU Wei，WANG Xinqing,YAN Xiaoxiang，et al.Estimation of wideband coherent source number based on spatial smoothing and transformed Gerschgorin radii[J].Journal of Chongqing University of Posts and Telecommunications(Natural Science Edition),2014(3):346-351.(in Chinese)

[8] CHENG W,LEE S,ZHANG Z,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.

[9] SHOJAFAR M,ABOLFAZLI S,MOSTAFAEI H,et al.Improving channel assignment in multi-radiowireless mesh networks with learning automata[J].Wireless Personal Communications,2015，82(1):61-80.

[10] SHOJAFAR M,CORDESCHI N,BACCARCLLI E.Energy-efficient adaptive resource management for real-time vehicular cloud services[J].IEEE Transactions on Cloud Computing,2016,2(3):1-14.

[11] 伍警,蔣鴻宇,李兵,等.非均勻天線組陣中引入對(duì)角加載的合成算法[J].電訊技術(shù),2017,57(6):710-715.

WU Jing，JIANG Hongyu，LI Bing,et al.A combining algorithm using diagonal loading for non-uniform antenna arraying[J].Telecommunication Engineering,2017,57(6):710-715.(in Chinese)

[12] CHENG Q,PAL P,TSUJI M,et al.An MDL algorithm for detecting more sources than sensors using outer-products of array output[J].IEEE Transactions on Signal Processing,2014,62(24):6438-6453.

[13] 謝紀(jì)嶺，司錫才.基于協(xié)方差矩陣對(duì)角加載的信源數(shù)估計(jì)方法[J].系統(tǒng)工程與電子技術(shù)，2008,30(1)：46-49.

XIE Jiling,SI Xicai.Determining the number of sources based on diagonal loading to the covariance matrix[J].Systems Engineering & Electronics,2008,30(1)：46-49.(in Chinese)

[14] HAN K,NEHORAI A.Improved source number detection and direction estimation with nested arrays and ULAs using jackknifing[J].IEEE Transactions on Signal Processing,2013,61(23):6118-6128.

[15] HANSEN B E,RACINE J S.Jackknife model averaging[J].Journal of Econometrics,2012,167(1):38-46.