代 振,吳吉偉,尹美方
(海軍士官學(xué)校,安徽 蚌埠,233000)
連續(xù)波(continuous wave,CW)信號作為一種比較簡單的發(fā)射信號波形,在雷達(dá)、聲吶以及聲自導(dǎo)魚雷中受到了廣泛應(yīng)用[1-3]。在高斯白噪聲背景下,檢測CW 信號的最佳檢測器為匹配濾波器(matched filter,MF)[4]。然而在實(shí)際檢測過程中,背景通常是有色的。高斯色噪聲背景下的最佳檢測器則為廣義匹配濾波器(generalized matched filter,GMF)。GMF 的思路是通過一定的變換,將色噪聲背景下的信號檢測問題轉(zhuǎn)化為白噪聲背景下的檢測問題,從而應(yīng)用MF 完成檢測[5]。
GMF 的核心在于對色噪聲進(jìn)行預(yù)白處理。常見的預(yù)白化方法主要包括Karhunen-Loève 變換法、Cholesky 分解法以及增強(qiáng)現(xiàn)實(shí)(augmented reality,AR)預(yù)白化法等[6-8]。但除AR 預(yù)白化法外,其他預(yù)白化法都要求色噪聲的協(xié)方差矩陣是先驗(yàn)已知的,這在實(shí)際過程中通常難以滿足。AR 預(yù)白化法雖然可以根據(jù)接收信號實(shí)時(shí)估計(jì)色噪聲的協(xié)方差矩陣,但也存在定階困難和計(jì)算量大等問題[9]。
為此,文中提出一種基于權(quán)向量的信號檢測方法,該方法在大數(shù)據(jù)記錄下(采樣點(diǎn)數(shù)足夠多時(shí))檢測CW 信號時(shí),無需進(jìn)行預(yù)白處理,并且檢測性能和GMF 相當(dāng)。
考慮如下檢測模型
式中:xn是接收信號;sn是待檢測信號;wn是零均值的平穩(wěn)高斯色噪聲,其協(xié)方差矩陣為Cw,是對稱正定toeplitz 矩陣。
易知X是N維高斯隨機(jī)矢量,在假設(shè)H0和H1下其概率分布為
對于上述信號檢測模型,根據(jù)經(jīng)典的統(tǒng)計(jì)信號檢測理論可知,其最佳檢測器為GMF,對應(yīng)的檢驗(yàn)統(tǒng)計(jì)量為[10]
由式(3)可以看出,GGMF是接收信號矢量X的線性組合,是一個(gè)線性檢測器。因此基于權(quán)向量構(gòu)造一個(gè)線性檢驗(yàn)統(tǒng)計(jì)量,即
式中,β=[β0,β1,···,βN-1]T為權(quán)向量。
由于Gβ也是X的線性組合,易知Gβ是一維高斯隨機(jī)變量,其概率分布為
由式(5)可知,在不同假設(shè)下Gβ的均值不等而方差相等。可以證明,以Gβ作為檢驗(yàn)統(tǒng)計(jì)量時(shí),其虛警概率Pf與檢測概率Pd總滿足[10]
由式(6)可以看出,在給定虛警概率時(shí),信號的檢測概率完全由檢測指數(shù)d2決定,其計(jì)算公式為
式中,E(G/Hj)與Var(G/Hj)分別是假設(shè)Hj(j=0,1)下檢驗(yàn)統(tǒng)計(jì)量G的均值與方差。
結(jié)合式(5)和式(7)可得G=Gβ=βTX時(shí)的檢測指數(shù)為
在雷達(dá)、聲吶以及聲自導(dǎo)魚雷中,通常采用NP 準(zhǔn)則作為最佳判決準(zhǔn)則,即給定虛警概率,使得檢測概率最大。由式(6)可知在虛警概率一定時(shí),檢測概率是檢測指數(shù)d2的單調(diào)遞增函數(shù),因此使得d2取最大值時(shí)的權(quán)向量就是N-P 準(zhǔn)則下的最優(yōu)權(quán)向量。下面推導(dǎo)最優(yōu)權(quán)向量的表達(dá)式。
在式(8)兩端同時(shí)對權(quán)向量 β求導(dǎo)可得
由于常系數(shù)對信號檢測沒有影響,為方便起見,可直接取最優(yōu)權(quán)向量為,將其代入式(4)可得最優(yōu)檢驗(yàn)統(tǒng)計(jì)量為
對比式(3)和式(11)可以發(fā)現(xiàn),根據(jù)最優(yōu)權(quán)向量構(gòu)造的最優(yōu)檢驗(yàn)統(tǒng)計(jì)量Gm和GGMF的表達(dá)式完全一樣,這與GMF 是高斯色噪聲背景下最佳線性檢測器的結(jié)論一致。
顯然,GMF 的檢測性能同樣可由式(12)確定的檢測指數(shù)決定。
當(dāng)待檢測信號sn為CW 信號時(shí),可寫為
式中:fc為歸一化后的信號頻率;A、θ分別為CW信號的幅值與相位。對于單次檢測而言,幅值A(chǔ)與相位 θ通常未知,但其值是確定的,因此可視為確知信號檢測。
如果是在高斯白噪聲背景下檢測CW 信號,其最佳檢測器為MF,檢驗(yàn)統(tǒng)計(jì)量為
對比式(14)和式(3)可以發(fā)現(xiàn),GMF 相比于MF 多出一個(gè)因子。由于Cw是色噪聲的協(xié)方差矩陣,所以的存在相當(dāng)于對色噪聲先進(jìn)行了預(yù)白處理,然后再進(jìn)行MF 檢測。這也是GMF 被稱為廣義匹配濾波的原因。
但是,對于色噪聲背景下的CW 信號檢測,如果直接應(yīng)用GMF,綜合式(11)和式(12)可知,無論是檢驗(yàn)統(tǒng)計(jì)量計(jì)算還是檢測性能分析,都不可避免地涉及到。然而,Cw一般不是先驗(yàn)已知的,通常需要進(jìn)行實(shí)時(shí)估計(jì)得到,這將大大增加檢測的計(jì)算量。為此考慮選擇合適的權(quán)向量,以避免預(yù)白處理,也就是避免計(jì)算。
顯然,一個(gè)最簡單的方法是繼續(xù)用MF 在色噪聲背景下檢測CW 信號,這意味著直接取權(quán)向量為
式(15)表達(dá)式雖然簡單,但色噪聲下的檢測性能并不是最優(yōu)的,需要進(jìn)一步分析。
由于MF 等價(jià)于權(quán)向量取為β=S,將其代入式(8)可得MF 的檢測指數(shù)為
定義MF 和GMF 的檢測指數(shù)之比為ηd,可知其表達(dá)式為
顯然 ηd<1,但如果 ηd越接近于1,表明MF 的檢測性能越接近于GMF。
下面證明,當(dāng)接收信號足夠長,即在大數(shù)據(jù)記錄下 ηd是等于1 的。
首先,對于CW 信號有
其次,根據(jù)文獻(xiàn)[10]可知,在大數(shù)據(jù)記錄下有
對于CW 信號而言,其周期譜圖只在fi=fc和fi=1-fc(實(shí)際上對應(yīng)的頻率是 -fc)處有非零值,而在其他頻率處的周期譜圖值均為零,并且有I(fc)=I(1-fc)。同時(shí)由于P(f)是偶函數(shù),故有P(fc)=P(1-fc)=P(-fc)。將二者代入式(19)可得
根據(jù)Parseval 定理,確知信號的時(shí)域能量等于其頻域能量,故有
將式(18)、(20)和(21)代入式(17)可得
表明在大數(shù)據(jù)記錄下,MF 和GMF 的檢測指數(shù)相等,即MF 和GMF 的檢測性能是一樣的。
假設(shè)高斯色噪聲背景wn可以用M階AR 模型來擬合,即
式中:ak為AR 模型系數(shù);un是均值為0、方差為1 的高斯白噪聲。
取采樣頻率為1 600 Hz,CW 信號歸一化頻率為0.15。采用4 階AR 模型來模擬高斯色噪聲,AR 模型參數(shù)取[-0.38,0.36,0.16,0.17]。采樣點(diǎn)數(shù)為N,仿真結(jié)果如下。
圖1 給出了 ηd與采樣點(diǎn)數(shù)N的關(guān)系。可以明顯看出,當(dāng)采樣點(diǎn)數(shù)N較小時(shí),ηd取值與1 相差較大;隨著N的增加,ηd取值也迅速增加;當(dāng)N>300時(shí),ηd取值較為平穩(wěn),接近于1。

圖1 N-ηd曲線Fig.1 Curve of N-ηd
約束虛警概率為Pf=10-3,取采樣點(diǎn)數(shù)N分別為50、100、500,統(tǒng)一進(jìn)行20 000 次蒙特卡洛實(shí)驗(yàn),MF 與GMF 的檢測性能對比見圖2。圖中橫坐標(biāo)為信噪比,記為SSNR,且SSNR=10log(A2/σ2w)。

圖2 MF 與GMF 檢測性能對比Fig.2 Comparison of detection performance for MF and GMF
從圖2 首先可以看出,當(dāng)采樣點(diǎn)數(shù)比較少時(shí)(N=50),GMF 的檢測性能優(yōu)于MF,前者比后者高出1 dB 左右,這與GMF 是高斯色噪聲背景下的最優(yōu)檢測器一致。同時(shí)還可以看出,隨著采樣點(diǎn)數(shù)的增加,MF 的檢測性能逐漸逼近GMF,當(dāng)采樣點(diǎn)數(shù)較大時(shí)(N=500),MF 與GMF 的檢測性能曲線幾近完全重合,這與文中的結(jié)論一致。
針對高斯色噪聲背景下的信號檢測問題,提出了一種基于權(quán)向量的檢驗(yàn)統(tǒng)計(jì)量構(gòu)造方法,推導(dǎo)了N-P 準(zhǔn)則下的最優(yōu)權(quán)向量,證明根據(jù)最優(yōu)權(quán)向量確定的檢測器與GMF 一致。將權(quán)向量應(yīng)用到CW 信號檢測中,證明了大數(shù)據(jù)記錄下MF 和GMF的檢測性能是等價(jià)的,仿真結(jié)果證明所提方法在檢測CW 信號時(shí)可以避免預(yù)白化處理,明顯提高檢測效率。但是,檢測調(diào)頻信號時(shí),仍然要進(jìn)行預(yù)白化處理,下一步可針對色噪聲背景下調(diào)頻信號的檢測,研究相應(yīng)權(quán)向量的選擇。