史 飛，?；郜?，王 濟(jì)，呂 鵬

(1.河北省電磁頻譜認(rèn)知與管控重點(diǎn)實(shí)驗(yàn)室，河北 石家莊 050081;2.中國(guó)電子科技集團(tuán)公司第五十四研究所，河北 石家莊 050081)

0 引言

頻譜觀測(cè)數(shù)據(jù)是人們感知電磁環(huán)境的重要依據(jù)之一，通過(guò)對(duì)頻譜觀測(cè)數(shù)據(jù)的分析處理，能夠掌握其中信號(hào)的整體分布情況和變化規(guī)律，進(jìn)而支撐認(rèn)知無(wú)線通信、無(wú)線電監(jiān)測(cè)、通信對(duì)抗等多個(gè)應(yīng)用領(lǐng)域。基于頻譜觀測(cè)數(shù)據(jù)的信號(hào)檢測(cè)通常包括功率譜估計(jì)、數(shù)據(jù)預(yù)處理、檢測(cè)算子設(shè)計(jì)等，以獲取信號(hào)分布參數(shù)，如信噪比、中心頻率、信號(hào)帶寬。復(fù)雜電磁環(huán)境中，信號(hào)檢測(cè)范圍內(nèi)各通信節(jié)點(diǎn)分布位置、發(fā)射功率及傳播路徑等各不相同，檢測(cè)節(jié)點(diǎn)接收到的各個(gè)信號(hào)強(qiáng)度存在較大差異；通信系統(tǒng)規(guī)格類型眾多，通信信號(hào)占用頻帶有寬有窄；有效的頻譜資源日趨緊張，帶內(nèi)信號(hào)疏密分布也差異較大。面對(duì)上述諸多問(wèn)題，相關(guān)學(xué)者針對(duì)各自應(yīng)用場(chǎng)景，設(shè)計(jì)了相應(yīng)的信號(hào)盲檢測(cè)方法。其中，能量檢測(cè)[1]工程應(yīng)用廣泛，但難以設(shè)定最優(yōu)能量門限以同時(shí)適應(yīng)高信噪比信號(hào)與低信噪比信號(hào)；形態(tài)學(xué)處理方法[2]能夠跟蹤噪底變化趨勢(shì)，但需要進(jìn)行形態(tài)尺度的調(diào)優(yōu)；子帶劃分法[3]需要事先確定子帶分布；匹配檢測(cè)[4]、循環(huán)平穩(wěn)分析[5]通常用于窄帶信號(hào)檢測(cè)分析。因此，本文提出了卡爾曼濾波聯(lián)合累積的數(shù)據(jù)平滑方法，提取了能夠良好反映信號(hào)分布的曲率特征，設(shè)計(jì)了高效靈活的正則匹配檢測(cè)算子，實(shí)驗(yàn)仿真分析證明該方法具有良好的低信噪比適應(yīng)性、檢測(cè)動(dòng)態(tài)范圍及準(zhǔn)確性。

1 盲信號(hào)檢測(cè)

本文提出的盲檢測(cè)方法，首先對(duì)頻譜數(shù)據(jù)進(jìn)行卡爾曼濾波處理，融合利用頻譜數(shù)據(jù)預(yù)測(cè)量與觀測(cè)量的統(tǒng)計(jì)信息降低頻譜估計(jì)方差，從而達(dá)到平滑效果以更好地反映頻譜內(nèi)信號(hào)的能量分布。隨后，將能量分布的變化特征轉(zhuǎn)化為曲率描述，減少噪聲非平穩(wěn)特性引入的頻譜緩變對(duì)檢測(cè)性能的影響。在此基礎(chǔ)上，對(duì)曲率特征進(jìn)行符號(hào)化表達(dá)，并采用高效、靈活的正則表達(dá)式進(jìn)行信號(hào)特征序列匹配，實(shí)現(xiàn)信號(hào)盲檢測(cè)。

1.1 卡爾曼濾波平滑

工程中常用的平滑濾波方法有滑動(dòng)平均[6]、FFT濾波[7]、形態(tài)濾波[8]等，但此類方法都沒(méi)有充分利用通信信號(hào)自身的先驗(yàn)信息?？紤]接收信號(hào)r(t)可表示為：

(1)

式中，si(t)為頻率相互不重疊的第i個(gè)發(fā)送信號(hào)，n(t)為高斯白噪聲，K為信號(hào)個(gè)數(shù)。因此，接收功率譜R(ω)為：

(2)

由貝葉斯濾波演化而來(lái)的卡爾曼濾波[9]，能夠有效融合觀測(cè)信息與預(yù)測(cè)信息，通過(guò)觀測(cè)統(tǒng)計(jì)不斷修正預(yù)測(cè)統(tǒng)計(jì)，在服從正態(tài)分布的不確定性干擾中發(fā)現(xiàn)最優(yōu)狀態(tài)，減小預(yù)測(cè)狀態(tài)估計(jì)方差，從而達(dá)到平滑效果。因此，可分別建立離散頻譜預(yù)測(cè)方程和觀測(cè)方程為：

Xk=FXk-1+Qk，

Yk=HXk+Rk,

(3)

(4)

(5)

以包含3個(gè)不同信噪比、不同帶寬QPSK信號(hào)的寬帶采樣數(shù)據(jù)為例，數(shù)據(jù)長(zhǎng)度32 768點(diǎn)，滑動(dòng)窗長(zhǎng)4 096點(diǎn)，窗重疊率50%，通過(guò)Pwelch方法估計(jì)的功率譜如圖1所示，經(jīng)過(guò)卡爾曼平滑后的功率譜如圖2所示?？梢钥闯?，頻譜中各信號(hào)分布趨勢(shì)保持一致的同時(shí)，抖動(dòng)等干擾特征明顯減小，達(dá)到了良好的平滑效果。

圖1 輸入原始功率譜

圖2 平滑后功率譜

1.2 曲率特征分析與提取

觀測(cè)的頻譜數(shù)據(jù)可以認(rèn)為是在不同頻率上的順序估計(jì)值排列而成的一組有序的數(shù)字序列[10]。此時(shí)，可通過(guò)綜合應(yīng)用多個(gè)專業(yè)領(lǐng)域的數(shù)據(jù)分析方法，挖掘數(shù)據(jù)序列中隱含的各種信息，典型的一些順序序列特征提取方法有鏈碼特征[11]、樣條方法[12]、尺度空間特征表示[13]等。對(duì)于頻譜觀測(cè)數(shù)據(jù)，為了完成信號(hào)檢測(cè)，需要從中發(fā)現(xiàn)、提取能夠反映信號(hào)功率分布變化的特征，如能夠反映信號(hào)起始的功率上升特征、信號(hào)結(jié)束的功率下降特征等，本文設(shè)計(jì)了基于最小二乘圓擬合的曲率特征提取方法，用于區(qū)分信號(hào)邊界。

(6)

曲線f的曲率κ定義為：

κ=|f"|/(1+f′2)3/2。

(7)

x2+y2+ax+by+c=0,

(8)

則基于{pi}的擬合誤差E為：

(9)

(10)

(11)

可解得：

(12)

此時(shí)，可以將曲率超過(guò)門限的點(diǎn)集標(biāo)識(shí)為起始部分并用符號(hào)∩表示，將低于門限的點(diǎn)集標(biāo)識(shí)為終止部分并用符號(hào)∪表示，對(duì)于曲率接近零的部分認(rèn)為噪底或信號(hào)平穩(wěn)部分并用符號(hào)≈表示，從而得到符號(hào)化特征序列ζ。

1.3 基于正則匹配的特征檢測(cè)

正則表達(dá)式(Regular Expression)[14]定義了字符串搜索匹配的一種方法，可以用來(lái)檢測(cè)一個(gè)字符串中是否包含一個(gè)特定子字符串。與字符串精確匹配方法相比，基于正則表達(dá)式匹配規(guī)則定義更靈活便捷。

在得到的符號(hào)化特征序列ζ中，可定義一個(gè)信號(hào)檢測(cè)準(zhǔn)則表達(dá)式為ψ：

ψ=′∩{I,}·*?∪{J,} ′,

(13)

式中，I、J均為正整數(shù)，定義了信號(hào)起始、終止部分的最小持續(xù)長(zhǎng)度。圖3顯示了基于該準(zhǔn)則的頻譜檢測(cè)結(jié)果，能夠完成不同信噪比、帶寬及間距分布信號(hào)的準(zhǔn)確檢測(cè)。更為復(fù)雜的檢測(cè)條件可基于正則表達(dá)式規(guī)則進(jìn)行靈活定義。

圖3 匹配檢測(cè)結(jié)果

2 仿真實(shí)驗(yàn)及性能對(duì)比

對(duì)在高斯色噪聲信道條件下的寬帶頻譜進(jìn)行仿真處理，信號(hào)在60 MHz帶寬范圍內(nèi)隨機(jī)分布于50個(gè)頻點(diǎn)，調(diào)制速率9.6～500 ksps，調(diào)制方式包括BPSK、QPSK、MSK，信噪比范圍-2～15 dB，信噪比動(dòng)態(tài)范圍包括0 dB、10 dB，信號(hào)間隔范圍2.5～250 kHz，改變各信號(hào)的頻率位置、調(diào)制速率、調(diào)制樣式、信噪比等參數(shù)，在不同條件下進(jìn)行1 000次蒙特卡洛檢測(cè)仿真，當(dāng)各信號(hào)檢測(cè)中心頻率偏差不超過(guò)5%帶寬時(shí)，判定檢測(cè)結(jié)果正確。采用傳統(tǒng)的能量檢測(cè)法與本文提出的檢測(cè)方法對(duì)比，結(jié)果如圖4所示?？梢园l(fā)現(xiàn)，本文方法在大于2 dB信噪比時(shí)，信號(hào)檢測(cè)準(zhǔn)確率即可達(dá)到100%，顯著優(yōu)于能量檢測(cè)法，且具有良好的信號(hào)檢測(cè)動(dòng)態(tài)范圍，而能量檢測(cè)法在低信噪比時(shí)受信號(hào)動(dòng)態(tài)變化影響，檢測(cè)性能進(jìn)一步惡化。

圖4 檢測(cè)準(zhǔn)確率曲線

3 結(jié)束語(yǔ)

本文設(shè)計(jì)了一種基于卡爾曼平滑的正則匹配頻譜搜索方法，能夠充分利用信號(hào)功率譜形狀的先驗(yàn)信息通過(guò)卡爾曼濾波實(shí)現(xiàn)對(duì)頻譜觀測(cè)數(shù)據(jù)的有效平滑。在此基礎(chǔ)上，提出了針對(duì)頻譜累積分布變化的曲率特征表達(dá)方法，并完成符號(hào)化轉(zhuǎn)化處理，最后基于靈活的正則表達(dá)式設(shè)計(jì)實(shí)現(xiàn)信號(hào)的準(zhǔn)確檢測(cè)。該方法能夠準(zhǔn)確反映頻譜觀測(cè)數(shù)據(jù)中隱含的信號(hào)分布趨勢(shì)規(guī)律，克服噪聲干擾的影響，相比傳統(tǒng)能量門限檢測(cè)方法，具有更優(yōu)的搜索動(dòng)態(tài)范圍、檢測(cè)準(zhǔn)確率和自適應(yīng)性。