張立民,鐘兆根, ,彭 耿

(1.海軍航空工程學(xué)院 電子信息工程系,山東 煙臺(tái) 264001;2.海軍裝備研究院,北京102249)

1 引 言

為了獲取詳細(xì)的敵方衛(wèi)星通信鏈路情報(bào)信息,進(jìn)而實(shí)施有效的干擾,其前提是對(duì)接收到的混合信號(hào)進(jìn)行有效檢測(cè)。針對(duì)現(xiàn)有復(fù)雜多變的衛(wèi)星信號(hào)環(huán)境特點(diǎn),如何實(shí)現(xiàn)非合作低信噪比條件下衛(wèi)星信號(hào)的快速盲檢測(cè),是現(xiàn)階段急需解決的問題。

目前,國內(nèi)外學(xué)者就這一問題已經(jīng)提出了許多行之有效的方法[1-2],但這些方法大都需要利用信號(hào)的某些先驗(yàn)知識(shí),不是嚴(yán)格意義上的盲檢測(cè)。針對(duì)像能量法、功率譜法等傳統(tǒng)盲信號(hào)檢測(cè)方法不能適應(yīng)低信噪比等缺陷,文獻(xiàn)[3-5]提出了一些針對(duì)性較強(qiáng)的方法,它們?cè)谝欢ǔ潭壬夏芴岣咝盘?hào)的檢測(cè)性能,但也存在著適應(yīng)信號(hào)類型單一、計(jì)算復(fù)雜和不能適應(yīng)多信號(hào)環(huán)境等缺陷。為了更好地解決這一問題,文獻(xiàn)[6]在前人研究的基礎(chǔ)上,提出了一種快速盲信號(hào)檢測(cè)方法,雖然該方法能滿足現(xiàn)實(shí)設(shè)備的需求,但是要應(yīng)用在衛(wèi)星等高速運(yùn)行的環(huán)境下,該方法計(jì)算量還是相對(duì)較高。為此,本文改進(jìn)了該方法的實(shí)現(xiàn)步驟,并從理論上給出了算法的檢測(cè)性能分析,進(jìn)一步驗(yàn)證了算法的可行性,并減少了計(jì)算量。

由于通常采用的無線電偵察接收機(jī)都是寬開的,多個(gè)信號(hào)可能同時(shí)進(jìn)入到接收機(jī)。本文在對(duì)上述方法進(jìn)行充分分析的基礎(chǔ)上,根據(jù)觀測(cè)矩陣奇異值與自協(xié)方差矩陣特征值之間的關(guān)系,提出了一種基于奇異值分解的多衛(wèi)星信號(hào)快速盲檢測(cè)方法,仿真結(jié)果驗(yàn)證了該方法對(duì)低信噪比、多信號(hào)環(huán)境的適應(yīng)性。同時(shí),該方法的計(jì)算量相對(duì)特征值方法大大減少,更適合應(yīng)用在星載設(shè)備上。

2 檢測(cè)算法的提出

本節(jié)在建立的信號(hào)模型基礎(chǔ)上,給出一種基于奇異值分解的盲檢測(cè)算法。

2.1 接收信號(hào)模型

針對(duì)現(xiàn)有的復(fù)雜電磁環(huán)境,衛(wèi)星偵察接收機(jī)的離散觀測(cè)信號(hào)可表示為

式中,I≥0為信號(hào)個(gè)數(shù);n(k)是零均值、方差為σ2的加性高斯白噪聲;si(k)為第i個(gè)數(shù)字調(diào)制衛(wèi)星信號(hào),且認(rèn)為各信號(hào)間、信號(hào)與噪聲間都是不相關(guān)的。

假設(shè)離散采樣之后的觀測(cè)數(shù)據(jù)x(k)含有N個(gè)采樣點(diǎn),即 k=1,2,3,…,N。為描述方便,首先對(duì)數(shù)據(jù)進(jìn)行零均值化處理,然后可利用重疊分段的方法,把長度為N的觀測(cè)數(shù)據(jù)分成N-M+1個(gè)相互重疊、長度為M、重疊長度為M-1的信號(hào)片段,為不失一般性,可令N-M+1≥M,將分段得到的數(shù)據(jù)寫成如下M×(N-M+1)維的Hankel矩陣形式:

則信號(hào)檢測(cè)問題可轉(zhuǎn)化為如下二元假設(shè)檢驗(yàn):

式中,P為檢測(cè)統(tǒng)計(jì)量,T1為檢測(cè)門限。

2.2 檢測(cè)量的設(shè)計(jì)

針對(duì)上述接收信號(hào)模型,在奇異值分解定理的基礎(chǔ)上,對(duì)觀測(cè)矩陣做如下分解:

其中,U和V為酉矩陣,VH為 V的共軛轉(zhuǎn)置,Δ為準(zhǔn)對(duì)角矩陣:

式中 ,D=diag(α1,α2,α3,…,αM),對(duì)角元素為 X 的非零奇異值 ,且 α1≥α2≥α3≥…≥αM>0。

下面給出觀測(cè)矩陣奇異值與自協(xié)方差矩陣特征值的關(guān)系,首先定義自協(xié)方差矩陣為

式中,符號(hào) E{·}表示期望。將式(4)代入式(6)可求得RX的特征值分解為

式中,□=ΔΔH也為對(duì)角矩陣,其對(duì)角元素為 RX的特征值,于是可知觀測(cè)矩陣的奇異值與自協(xié)方差矩陣的特征值是平方關(guān)系,RX的特征值為α2i,1≤i≤M。同時(shí),由式(1)和式(6)可知:

式中,RS為信號(hào)的自協(xié)方差矩陣,I為M階單位矩陣,且RS為半正定的Hermitian矩陣,假設(shè)其秩為D(D

式中,Λ為D階對(duì)角矩陣,其主對(duì)角線元素為RS的非零特征值(設(shè)為 λi(i=1,2,3,…,D),且按降序排列)。將該式代入式(8)可得RX的特征值為

根據(jù)陣列信號(hào)處理的相關(guān)定義可知,D個(gè)較大特征值對(duì)應(yīng)的特征向量張成空間稱為信號(hào)子空間,M-D個(gè)較小特征值對(duì)應(yīng)的特征向量張成空間稱為噪聲子空間,且信號(hào)子空間與噪聲子空間正交。

在前期的研究工作中,文獻(xiàn)[7]通過估計(jì)信號(hào)子空間的維數(shù)大小來判斷信號(hào)的有無,但是在低信噪比條件下,這種維數(shù)估計(jì)的方法可靠性不強(qiáng),一旦估計(jì)錯(cuò)誤整個(gè)算法就可能會(huì)失效。針對(duì)這一問題,本文采用文獻(xiàn)[6]設(shè)計(jì)的類似檢測(cè)統(tǒng)計(jì)量,通過理論分析并推導(dǎo)檢測(cè)算法的性能,使得其檢測(cè)性能不受空間維數(shù)估計(jì)的影響,即

式中,E為待估計(jì)的信號(hào)子空間維數(shù)。

2.3 門限等參數(shù)的確定

對(duì)于上面建立的檢測(cè)統(tǒng)計(jì)量,首先需要估計(jì)出信號(hào)子空間的維數(shù)E,然后通過推導(dǎo)得出檢測(cè)量的統(tǒng)計(jì)特性,進(jìn)而確定檢測(cè)門限。

信號(hào)子空間的估計(jì)方法在陣列信號(hào)處理中已比較成熟,如基于信息論準(zhǔn)則的方法和蓋氏圓準(zhǔn)則的方法等。本文采用文獻(xiàn)[8]給出的MDL方法來估計(jì)子空間維數(shù)E,不同的是,文中不再用特征值判決函數(shù)進(jìn)行計(jì)算,而是采用奇異值進(jìn)行處理,這種替換在文獻(xiàn)[9]中已得到了證明:

已被Anderson等[10]證明,對(duì)實(shí)自協(xié)方差矩陣有

式中,N(0,1)表示標(biāo)準(zhǔn)正態(tài)分布。對(duì)上式進(jìn)行變換,化為檢測(cè)量P的表達(dá)式,可得

于是在H0成立時(shí),μ0=E(P H0),通過給定的虛警概率Pfa就可以得到檢測(cè)門限gT,從而得出絕對(duì)檢測(cè)門限 T1=μ0+gT·σ1。

針對(duì)噪聲的隨機(jī)性而導(dǎo)致的信號(hào)漏檢問題,本文通過做重復(fù)實(shí)驗(yàn)來有效抑制信號(hào)漏檢,給出一種平均符號(hào)檢測(cè)方法:首先對(duì)上述的檢測(cè)量做100次重復(fù)實(shí)驗(yàn),分別為 Pi、T1i(i=1,2,3,…,100),則這100個(gè)值大于絕對(duì)門限的百分比為

式中,sgn為符號(hào)函數(shù),最后給出置信概率門限 T2(本文選取為80%),當(dāng) TRatio大于T2時(shí),則以大概率相信有信號(hào)存在,否則就沒有信號(hào)。

2.4 算法描述

下面給出盲信號(hào)檢測(cè)算法的實(shí)現(xiàn)步驟:

(1)首先對(duì)觀測(cè)信號(hào)進(jìn)行零均值化處理,按照式(2)給出的方式進(jìn)行分段處理,得出觀測(cè)矩陣 X;

(2)對(duì)觀測(cè)矩陣進(jìn)行奇異值分解,求出非零的奇異值 αi(α1≥α2≥α3≥…≥αM>0);

(3)采用式(12)～(13)的判決方法估計(jì)E;

(4)由E結(jié)合式(16)～(17)計(jì)算檢測(cè)量的均值和方差,進(jìn)而根據(jù) Pfa確定T1;

(5)根據(jù)式(11)計(jì)算檢測(cè)量P;

(6)重復(fù)步驟2～5做100次仿真實(shí)驗(yàn),根據(jù)式(18)計(jì)算平均符號(hào)概率TRatio;

(7)若TRatio大于檢測(cè)置信概率 T2則相信有信號(hào)存在,否則就沒有信號(hào)。

3 檢測(cè)性能和計(jì)算量分析

本節(jié)首先從理論上分析算法的檢測(cè)性能,根據(jù)檢測(cè)量與信噪比的關(guān)系得出本算法所能適應(yīng)的最低信噪比,然后對(duì)比算法的計(jì)算量。

3.1 檢測(cè)性能分析

由式(16)可知:

根據(jù)紐曼皮爾遜準(zhǔn)則,在給定虛警概率Pfa,可求出檢測(cè)門限T1,從而可得出檢測(cè)概率如下:

綜合式(19)對(duì)上式化簡可得

下面給出檢測(cè)量與信噪比的關(guān)系,進(jìn)而得出本算法所能適應(yīng)的最低信噪比。由文獻(xiàn)[11]可知,信噪比與觀測(cè)矩陣的奇異值滿足如下關(guān)系:

將式(11)代入上式整理得

若單次檢測(cè)到信號(hào),則P大于T1,即

從上面的式子可以看出,在給定虛警概率的條件下,單次檢測(cè)所能適應(yīng)的最低信噪比可表述為

令M=10、N從1000到9000,單次檢測(cè)所能適應(yīng)的最低信噪比如圖1所示。

圖1 不同條件下單次檢測(cè)所能適應(yīng)的最低信噪比Fig.1 The minimum SNR that a single detection can adapt

由該圖可知:仿真值與理論值之間雖存在一定的差異,但趨勢(shì)是基本吻合的,它們之間存在的差異主要是由奇異值的估計(jì)誤差等因素造成的。

3.2 計(jì)算量分析

由前面的分析可見,奇異值分解方法與文獻(xiàn)[6]給出的方法都可分為4步來計(jì)算:第一是生成觀測(cè)矩陣或協(xié)方差矩陣;第二是估計(jì)奇異值和特征值;第三是子空間維數(shù)估計(jì);第四是檢測(cè)比較分析。下面主要討論比較兩種方法的計(jì)算量,并假設(shè)加法、乘法和比較運(yùn)算的計(jì)算量相當(dāng),其結(jié)果如表1所示。

表1 計(jì)算量分析Table 1 Computational complexity analysis

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

本節(jié)首先就單個(gè)和多個(gè)信號(hào)時(shí)算法的檢測(cè)性能進(jìn)行仿真驗(yàn)證,然后將本文方法與典型盲檢測(cè)算法進(jìn)行性能對(duì)比,同時(shí)給出計(jì)算復(fù)雜度對(duì)比結(jié)果。

(1)仿真實(shí)驗(yàn)1:驗(yàn)證算法對(duì)單個(gè)信號(hào)的檢測(cè)性能

本次仿真實(shí)驗(yàn)分別驗(yàn)證算法對(duì)MPSK、MQAM信號(hào)的檢測(cè)性能,其參數(shù)設(shè)置為:N=5000,M=50,Pfa分別為 10-2、10-3、10-4,T2為 80%,采樣頻率200MHz,信號(hào)載頻70 MHz,碼元速率10Mb/s,擴(kuò)頻碼長1023,進(jìn)行10000次Monte Carlo仿真,其結(jié)果如圖2所示。

圖2 不同虛警概率時(shí)單個(gè)信號(hào)的盲檢測(cè)仿真結(jié)果Fig.2 Simulation results for different probability of false alarm

由圖 2可知:在 Pfa小于 10-2、SNR=-11 dB時(shí),MPSK、MQAM常用衛(wèi)星信號(hào)的檢測(cè)概率均達(dá)到90%以上,能較好滿足信號(hào)偵察的需求。

(2)仿真實(shí)驗(yàn)2:驗(yàn)證算法對(duì)多個(gè)信號(hào)的檢測(cè)性能

在Pfa為10-2、T2為80%、其他仿真參數(shù)同實(shí)驗(yàn)1,兩個(gè)信號(hào)混合的功率比為0 dB,頻譜重疊度分別為50%、100%。SNR為信號(hào)能量與噪聲能量比值的對(duì)數(shù),Monte Carlo仿真10000次,其結(jié)果如圖3所示。

圖3 信號(hào)個(gè)數(shù)、頻譜重疊度對(duì)算法性能的影響Fig.3 The influence of signal number and spectrum overlap

由圖3可知:信號(hào)個(gè)數(shù)、頻譜重疊度雖對(duì)算法性能有一定影響,但影響不大,究其原因是檢測(cè)量 P在多個(gè)信號(hào)情況下能將每個(gè)信號(hào)能量累加起來進(jìn)行考慮,即更充分地將信號(hào)能量集中起來利用。

(3)仿真實(shí)驗(yàn)3:與典型算法的性能對(duì)比

文獻(xiàn)[12]提出了一種較典型的DSSS信號(hào)盲檢測(cè)算法,選取其作為對(duì)比算法,仿真參數(shù)同實(shí)驗(yàn)1,Monte Carlo仿真10000次,其結(jié)果如圖4所示。

圖4 檢測(cè)性能對(duì)比曲線Fig.4 The comparison curve of detection performance

由圖4可知:相對(duì)于文獻(xiàn)[12](虛警概率為1%、SNR=-10 dB時(shí),檢測(cè)概率可達(dá)50%以上),本文算法的性能有明顯提高,且適應(yīng)的信號(hào)類型不僅局限于DSSS信號(hào),檢測(cè)性能的提高在于檢測(cè)量 P能將信號(hào)的能量有效集中起來利用。同時(shí),將本文的仿真結(jié)果與文獻(xiàn)[6]的仿真結(jié)果比較可知,兩者的性能差不多,但本文是在對(duì)檢測(cè)量進(jìn)行理論分析的基礎(chǔ)上得出的結(jié)論,其仿真結(jié)果更具備可靠性。

針對(duì)表1的計(jì)算量分析結(jié)果,令數(shù)據(jù)長度N從2048變化到20480時(shí),可得圖5所示的計(jì)算復(fù)雜度對(duì)比曲線。由圖可知:相對(duì)于特征值分解的盲檢測(cè)算法,本文的計(jì)算量一般為其1/10,具有相對(duì)較小的計(jì)算復(fù)雜度,更適合應(yīng)用在實(shí)時(shí)性要求較高的星載設(shè)備上。

圖5 與文獻(xiàn)算法的計(jì)算復(fù)雜度對(duì)比Fig.5 The comparison of computational complexity

5 結(jié) 論

針對(duì)日益復(fù)雜時(shí)變的多衛(wèi)星信號(hào)偵察環(huán)境,本文在分析、研究當(dāng)前通信信號(hào)盲檢測(cè)方法的基礎(chǔ)上,根據(jù)觀測(cè)信號(hào)奇異值與協(xié)方差矩陣特征值之間的關(guān)系,提出了基于奇異值分解的多衛(wèi)星信號(hào)快速盲檢測(cè)算法,并從理論上推導(dǎo)了檢測(cè)算法的檢測(cè)性能,給出了檢測(cè)量與信噪比的關(guān)系。從計(jì)算復(fù)雜度的對(duì)比分析可知,相比于文獻(xiàn)[6]的特征值方法,本文方法確實(shí)能在很大程度上減少計(jì)算量,并達(dá)到快速實(shí)時(shí)的效果。仿真結(jié)果表明,算法的性能有明顯提高,且適應(yīng)的信號(hào)類型不僅局限于直擴(kuò)信號(hào),相比于功率譜和能量檢測(cè)等傳統(tǒng)方法,該算法能將多個(gè)信號(hào)的能量集中起來進(jìn)行考慮,因此在多信號(hào)檢測(cè)時(shí)能適應(yīng)更低的信噪比。這為星載偵察設(shè)備的研制和空間信息對(duì)抗技術(shù)的研究提供了重要的技術(shù)支持,具有重要的意義。

[1]Barkat M.Signal detection and estimation[M].2nd ed.Norwood,MA,USA:Artech House,2005.

[2]HillD A,Bodie J B.Carrier detection of PSK signals[J].IEEE Transactions on Communications,2001,49(3):487-496.

[3]Cheng X,Vanhaverbeke F,Guan Y L.Blind combining for weighted energy detectuin of UWB signals[J].Electronics Letters,2011,47(1):55-57.

[4]Tsung Hui Chang,ChienWei Hsin,Wing Kin Ma,et al.A linear farctional semidefinite relaxation approach to maximumlikelihood detection of higher-order QAM OSTBC in unknown channels[J].IEEE Transactions on Signal Processing,2010,58(4):2315-2326.

[5]BarbieriA.Blind Per-State detection of DPSK over correlated fading channels[J].IEEE Transactions on Vehicular Technology,2010,59(5):2320-2327.

[6]彭耿,黃知濤,陸鳳波,等.基于特征值分析的通信信號(hào)快速盲檢測(cè)方法[J].電子學(xué)報(bào),2009,37(9):1950-1958.PENG Geng,HUANG Zhi-tao,LU Feng-bo,et al.A fast blind detection algorithm of communication signals based on eigenvalue analysis[J].Acta Electronica Sinica,2009,37(9):1950-1958.(in Chinese)

[7]沙志超,彭耿,徐神,等.基于信號(hào)子空間維數(shù)估計(jì)的單通道盲信號(hào)檢測(cè)[J].信號(hào)處理,2009,25(8A):94-97.SHA Zhi-chao,PENG Geng,XU Shen,et al.Single channel blind signal detection based on signal subspace dimension estimation[J].Signal Processing,2009,25(8A):94-97.(in Chinese)

[8]Jensen J,Heusdens R.Improved subspace-based singlechannel speech enhancement using generalized super-Gaussian priors[J].IEEE Transactions on Audio,Speech,and Language Processing,2007,15(3):862-872.

[9]劉玲,曾孝平,曾浩.基于奇異值的信源數(shù)估計(jì)方法[J].計(jì)算機(jī)工程,2009,35(12):34-37.LIU Ling,ZENG Xiao-ping,ZENG Hao.Signal source number estimation method based on singular value[J].Computer Engineering,2009,35(12):34-37.(in Chinese)

[10]Anderson T W.An introduction to multivariate statistical analysis[M].New York:Wiley,1958.

[11]詹亞鋒,曹志剛,馬正新.無線數(shù)字通信的盲信噪比估計(jì)[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,43(7):957-960.ZHAN Ya-feng,CAO Zhi-gang,MA Zheng-xin.Blind SNR estimations inwireless digital communications[J].Journalof Tsinghua University(Science and Technology),2003,43(7):957-960.(in Chinese)

[12]BurelG,Bouder C,Berder O.Detection of direct sequence spreadspectrum transmissionswithout prior knowledge[C]//Proceedings of 2001 IEEE Global Telecommunications Conference.San Antonio,USA:IEEE,2001:236-239.