徐 哲，熊 剛

（1.中國人民解放軍91977 部隊，北京 100001；2.中國電子科技集團公司第三十研究所，四川 成都 610041）

0 引言

直接序列擴譜（Direct Sequence Spread Spectrum，DSSS）信號具有頻帶寬、抗干擾性能好、傳輸速率高和隱蔽性強等多種特點[1]，已經(jīng)被軍事通信與民用通信領(lǐng)域廣泛采用，如無人機數(shù)據(jù)鏈通信系統(tǒng)、全球定位系統(tǒng)（Global Positioning System，GPS）、美國第三代國防衛(wèi)星通信系統(tǒng)（Defense Satellite Communication System III，DSCS III）、跟蹤與數(shù)據(jù)中繼衛(wèi)星系統(tǒng)（Tracking and Data RelaySatellite System，TDRSS）、碼分多址（Code Division Multiple Access，CDMA）通信和無線局域網(wǎng)（Wireless Local Area Network，WLAN）等[2]，由此低信噪比下DSSS 信號的檢測技術(shù)成為了頻譜監(jiān)測、網(wǎng)電對抗等應(yīng)用中的重要研究課題。直擴通信信號的存在性檢測為后續(xù)參數(shù)估計和解擴的前提，但在實際環(huán)境中電磁環(huán)境日益復(fù)雜，且可提供的目標(biāo)信號先驗信息較少，可以視作近似盲檢測范疇。

由于DSSS 信號具有低功率譜和低工作信噪比的性質(zhì)，其經(jīng)常表現(xiàn)為淹沒在電磁空間噪聲里，導(dǎo)致很多傳統(tǒng)的檢測方法在運用于DSSS 信號時效果不理想。這個問題引起了眾多研究機構(gòu)和學(xué)者的關(guān)注，例如文獻[3]中描述的方法都是基于已知DSSS 信號的周期特征和一些先驗假設(shè)條件的，否則無法檢測；文獻[4]是根據(jù)短碼DSSS 信號的符號循環(huán)平穩(wěn)特進行檢測，不能用于長碼DSSS 信號。因此，現(xiàn)有研究的不足也對DSSS 信號的檢測算法分析及改進提出了較高要求[5]。

針對上述問題，本文提出了一種DSSS 信號的檢測新方法，該算法具有檢測率較高，且不需要太多的目標(biāo)先驗信息的特點。該方法基于對接收的DSSS 信號頻域空間協(xié)方差矩陣優(yōu)化計算思路，能夠在低信噪比條件下正確檢測出DSSS 信號，穩(wěn)健性強。此外，該方法通過對信號的交叉相位線性特征進行深入分析，實用性較好，能夠適用于長碼直擴、短碼DSSS 信號，可以實現(xiàn)成功的檢測。

1 信號模型分析

下面對本文研究的DSSS 信號目標(biāo)對象系統(tǒng)和算法的數(shù)學(xué)模型進行說明。DSSS 信號的收發(fā)原理框圖如圖1 所示。

在通道下變頻、濾波和采樣處理后，接收到的DSSS 信號可以表示為：

式中：{d(kTc)}表示原始的信息符號；k為整數(shù)，k=1,2,3,…,n；{c(kTc)} 表示直接擴譜的碼序列；Tc表示擴譜碼間隔；S為DSSS 信號能量；{n(kTc)}表示信道高斯噪聲。擴譜碼序列可用{c1,…,cN}，ci∈{±1} 表示；相移鍵控（Phase Shift Keying，PSK）調(diào)制后得到的符號可用di表示，di∈{±1}。

短碼DSSS 信號和長碼DSSS 信號的區(qū)別為：短碼是在擴譜周期內(nèi)僅調(diào)制一個符號碼元[6]；長碼是指對應(yīng)于各個擴譜周期可調(diào)制多符號碼元。[7]長碼DSSS 信號的數(shù)學(xué)表達式為：

式中：dk表示擴譜符號序列，是獨立同分布的；n(t)為高斯噪聲；Ts為DSSS 信號的擴譜符號周期；cu（u=0.1,…,KN-1）表示其中的第u個碼元符號；Tc為符號長度，Tc=Ts/N；N表示長碼序列相對于符號的等分長度。長碼擴譜信號持續(xù)長度設(shè)為T0，則T0=KTs=KNTc；K表示數(shù)量，K=1,2,…,n。

在擴譜信號接收過程中，設(shè)濾波器的頻譜響應(yīng)為G(f)，T表示擴頻周期，γ(f)表示信號的頻域變換，f表示頻率，當(dāng)T取值較大時，DSSS 信號的平均功率表達式為：

其方差表達式為：

式中：R^(τ)表示信號的解析式，根據(jù)中心極限準(zhǔn)則，R^(τ)服從高斯分布。

由于信號的高斯均值為0。當(dāng)τ值較小時，可得到：

然后，對信號的檢測模型開展分析[8]。其二元假設(shè)檢驗的數(shù)學(xué)表達式為：

式中：H0和H1分別表示目標(biāo)信號不存在和目標(biāo)信號存在兩種情況；N表示接收的目標(biāo)信號樣本數(shù)量。設(shè)w[n]表示在時間和空間上不相關(guān)的、白色的、零均值高斯噪聲，其方差為σ2。s[n]表示待檢測的目標(biāo)信號，w[n]表示噪聲。

對于H1假設(shè)，如圖2 所示，DSSS 信號檢測可采用多通道結(jié)構(gòu)[9]，使用M個處理單元組成的空間分布陣列來進行分析。其中，F(xiàn)FT（Fast Fourier Transform）表示快速傅里葉變換。

圖2 本文基于相位線性分析的DSSS 信號多通道檢測模型

從上圖中可看出，本文中設(shè)計的DSSS 信號檢測接收模型表達式為：

式中：m表示處理單元的序號；Dm表示第m個通道和參考通道之間的時間偏移量；x[n]表示基帶源信號；c[n]表示長度為L的長碼擴譜碼。對于長碼DSSS 信號，L通常應(yīng)大于擴頻系數(shù)因子。式（7）的頻域變換計算式為：

式中：k表示頻率譜線的索引值；K表示頻域譜線的總數(shù)量（K≤N）；?表示卷積。在頻域中M個通道的輸出向量可表示為：

進一步得出，第k個頻域譜線的空間協(xié)方差矩陣可定義為：

式中：(·)H表示矩陣的共軛轉(zhuǎn)置；A[k]代表流形矩陣；I[k]表示單位矩陣。設(shè)S[k]表示接收的單一寬帶DSSS 信號矩陣，則信號的協(xié)方差矩陣為：

2 本文檢測算法分析

若存在DSSS 信號時，交叉通道協(xié)方差項的相位對應(yīng)于信號帶寬是呈線性的[10]。圖3 和圖4 分別給出了交叉協(xié)方差項的幅度響應(yīng)和相位響應(yīng)。

圖3 本文算法的交叉協(xié)方差項的幅度響應(yīng)

圖4 本文算法的交叉協(xié)方差項的相位響應(yīng)

從圖4 中可以看出，DSSS 信號的空間協(xié)方差相位具有明顯的線性特征。

本文的檢測思路主要是通過先構(gòu)建圖2 所示的多通道處理模型；然后對空間協(xié)方差矩陣進行計算，提取出空間分布通道的相位線性特征。在最小均方意義下，DSSS 信號的線性斜率與帶寬內(nèi)的相位響應(yīng)匹配，并對相位響應(yīng)函數(shù)實現(xiàn)最佳擬合，從而提取出線性度特征。在DSSS 信號存在的情況下，線性擬合的曲線與信號相位是基本重合的，而在無信號存在，僅有噪聲的情況下，擬合偏差量較大。

DSSS 信號是寬帶信號，在低信噪比環(huán)境中常常被背景噪聲掩蓋；因此，本文對DSSS 信號的檢測采用了寬帶信號分解的思路，信號被分解成多個窄帶通道，然后處理窄帶內(nèi)的各個信號，再將各頻帶的處理數(shù)據(jù)進行組合，得到最終結(jié)果。

（1）將DSSS 信號平均劃分為不重疊的多個通道子帶，其窄帶帶寬對應(yīng)的快拍數(shù)為B；

（2）計算快照B時間內(nèi)信號的離散傅里葉變換（Discrete Fourier Transform，DFT），得到如下空間估計協(xié)方差矩陣：

式中：B表示快拍數(shù)量，b=1,…,B；P表示每個快拍對應(yīng)的頻率譜線的索引數(shù)，p=1,…,P；Yb,p代表第b個快拍，第p個頻率分量對應(yīng)通道的輸出；進一步地，表示M×M維矩陣，且其協(xié)方差矩陣有P個。

（3）對于第p根頻率譜線，的交叉項相位值為：

式中：矩陣的行索引用i表示，i=1,…,M-1；列索引用j表示，j=2,…,M，i＜j。從而φp(i,j)可由P個頻率分量和相位響應(yīng)得出。信號帶寬對應(yīng)的交叉信道協(xié)方差項相位響應(yīng)如圖4 所示。當(dāng)檢測出DSSS 信號時，相位響應(yīng)曲線對于信號帶寬內(nèi)的部分都呈現(xiàn)線性特性。

（4）通過使用線性模型，可將全部頻率譜線對應(yīng)的相位響應(yīng)函數(shù)實現(xiàn)曲線擬合，得到線性模型的數(shù)學(xué)表達式為：

式中：F表示全部的頻率譜線集合變量；f表示頻率譜線索引；α(i,j)表示初始點；β(i,j)表示擬合曲線斜率。從信號中心頻率的兩側(cè)對稱地選擇信號帶寬內(nèi)的F頻點（F＜P）。式（15）中表示的模型也可以用第i個、第j個交叉項協(xié)方差矩陣表示為：

且有θ(i,j)表示未知參數(shù)向量，θ(i,j)=[α(i,j),β(i,j)]T；因此，各交叉項的擬合曲線參數(shù)可得出的數(shù)學(xué)估計表達式為：

基于通道相位線性分析的DSSS 信號檢測方法的主要處理思路步驟為：

（1）計算所有交叉信道協(xié)方差項的相位響應(yīng)和最佳擬合曲線之間的線性度偏差。為定量提取特征，可采用標(biāo)準(zhǔn)化的總均方根誤差，計算方式為：

（2）采用式（18）中的總誤差作為檢驗統(tǒng)計量，即各通道的信號相位響應(yīng)與擬合曲線之間的線性度偏差量。

（3）根據(jù)均方根值的分布選取合理的檢測閾值γthreshold，判決是否存在直擴DSSS 信號，如下所示：

γthreshold也與標(biāo)準(zhǔn)正態(tài)分布表和虛警率（False Alarm Rate，F(xiàn)AR）有關(guān)。TotalRMSE的值是對基于相位空間協(xié)方差矩陣的度量，并不依賴于噪聲功率。

3 仿真分析

本節(jié)對本文算法進行MATLAB 仿真試驗驗證。

首先設(shè)置仿真參數(shù)：接收的長碼/短碼DSSS信號的采樣率為200 MHz，載波頻率為60 MHz，擴頻碼元速率為10.24 Mb/s，樣本點數(shù)取為8 192，擴頻碼為M 序列。Monte Carlo 仿真次數(shù)設(shè)為1 000 次，信道噪聲為高斯白噪聲。

其次，對本文設(shè)計的基于通道相位線性分析的DSSS 信號檢測算法開展仿真，并選用過去的一些方法思路包括功率譜檢測算法、基于倒譜的檢測法和本文新方法進行比較，檢測性能曲線如圖5所示。

圖5 本文新算法與過去算法的性能對比曲線

分析圖5 可得出，本文的檢測算法可實現(xiàn)對DSSS 信號的正確檢測，且在低信噪比條件下，新方法的檢測效果相對優(yōu)于過去的功率譜法和基于倒譜的檢測方法，性能有所提高。

4 結(jié)語

在低信噪比、復(fù)雜電磁環(huán)境中的DSSS 信號檢測問題一直是非協(xié)作信號處理研究中的難點問題，需要對過去的一些算法加以改進，并探索新的方法。本文提出了一種基于通道相位線性分析的DSSS 信號檢測算法。仿真試驗驗證了新方法的性能，比起過去的功率譜算法、倒譜算法等具有更高的檢測概率，相對更優(yōu)，可為DSSS 信號的監(jiān)測提供一種更有效的手段，并對網(wǎng)電對抗技術(shù)領(lǐng)域的研究提供參考。