張朝星，趙興剛，韓曉軍

(1.解放軍32137部隊(duì)，河北 張家口 075000；2.解放軍66136部隊(duì)，北京 100042；3.解放軍31675部隊(duì)，河北 張家口 075000)

0 引 言

時(shí)頻分析[1]通過(guò)將一維時(shí)域信號(hào)轉(zhuǎn)換為二維時(shí)頻域圖像，從而提取有效的圖像特征來(lái)判斷有無(wú)目標(biāo)，是實(shí)現(xiàn)海雜波背景下目標(biāo)檢測(cè)[2]的一種有效方法。短時(shí)傅里葉變換(STFT)作為一種常用的時(shí)頻分布[3]，其基本思想是用一個(gè)確定寬度的窗函數(shù)沿時(shí)間軸逐段對(duì)信號(hào)進(jìn)行傅里葉分析，得到信號(hào)的譜圖。STFT具有許多優(yōu)點(diǎn)，然而也有其缺點(diǎn)和不足，根據(jù)不確定性原理[4]，其時(shí)間分辨率與頻率分辨率的乘積恒定，不能同時(shí)提高，并且窗的寬度不能根據(jù)不同頻率實(shí)現(xiàn)自適應(yīng)調(diào)整，時(shí)頻分辨力較低，導(dǎo)致使用STFT方法對(duì)海雜波背景下的弱小目標(biāo)進(jìn)行檢測(cè)時(shí)效果較差。

曲線擬合選擇一條合適的曲線來(lái)最佳地?cái)M合觀測(cè)數(shù)據(jù)，在歐式空間中，常見(jiàn)的曲線類(lèi)型有二維直線、拋物線、橢圓等。在功率譜的信息幾何理論[5]中，所有功率譜的集合可以看做一個(gè)黎曼流形[6]，則對(duì)回波數(shù)據(jù)進(jìn)行STFT處理后所得到的不同時(shí)間點(diǎn)的功率譜就對(duì)應(yīng)了流形上的若干點(diǎn)。類(lèi)比于歐式空間中的曲線擬合，本文提出了一種使用流形上的距離和測(cè)地線[7]對(duì)不同時(shí)間點(diǎn)的功率譜進(jìn)行曲線擬合以提高其時(shí)頻分辨力的方法。由于STFT是一種線性變換，我們選擇歐式空間中的直線在流形上的推廣(即測(cè)地線)對(duì)不同點(diǎn)的功率譜進(jìn)行擬合。在具體的實(shí)現(xiàn)過(guò)程中，歐式空間使得給定數(shù)據(jù)與所選擇曲線對(duì)應(yīng)點(diǎn)的歐式距離總的平方和最小(最小二乘法)[8]。轉(zhuǎn)換到流形上以后，實(shí)現(xiàn)方法不變，但歐式距離要替換為兩功率譜間的距離。我們選擇M-K運(yùn)輸問(wèn)題[9]中的運(yùn)輸距離作為兩功率譜之間的差異度量，得到了相應(yīng)的黎曼度規(guī)，對(duì)STFT后的譜圖實(shí)現(xiàn)了測(cè)地線擬合。擬合后的譜圖對(duì)于單一點(diǎn)的功率譜進(jìn)行了重構(gòu)，使其更加平滑和準(zhǔn)確。其次擬合所用測(cè)地線對(duì)不同時(shí)間點(diǎn)的功率譜實(shí)現(xiàn)了時(shí)域上的平滑，時(shí)間分辨率也得到提高，進(jìn)而整體提高了譜圖的時(shí)頻分辨力，STFT的檢測(cè)性能得到改善。

1 基于信息幾何的功率譜曲線擬合

1.1 運(yùn)輸距離及測(cè)地線

在上節(jié)中已提到，具體實(shí)現(xiàn)流形上的擬合時(shí)，要將最小二乘法中的歐式距離替換為2個(gè)不同時(shí)間點(diǎn)功率譜之間的距離，且該距離除了滿(mǎn)足對(duì)稱(chēng)性和三角不等式等距離的一般性質(zhì)外，還必須要滿(mǎn)足弱連續(xù)性[10]，以使得功率譜能夠收斂，即隨著2個(gè)功率譜的“靠近”或“遠(yuǎn)離”，該距離能相應(yīng)地減小或增大。因此，我們選用了M-K運(yùn)輸問(wèn)題中的運(yùn)輸距離來(lái)度量?jī)晒β首V之間的差異，對(duì)于2個(gè)歸一化的密度函數(shù)f0和f1，其運(yùn)輸距離定義如下[9]：

(1)

式中：X=[-π,π]；ψ(x)為一個(gè)f0和f1之間的保測(cè)度映射，即：|detψ(x)|f1(ψ(x))=f0(x)。

在一般情況下，ψ(x)可通過(guò)下式計(jì)算得到[9]：

F0(θ)=F1(ψ(θ))

(2)

(3)

當(dāng)f1是f0的一個(gè)小的增量時(shí)，可由運(yùn)輸距離得到一個(gè)黎曼度規(guī)，并稱(chēng)之為運(yùn)輸度規(guī)[11]：

(4)

(5)

定義了流形上的度規(guī)，就可得到f0和f1之間由梯度流決定的測(cè)地線[12]：

Fτ((1-τ)θ+τψ(θ))=F0(θ)

(6)

式中：ψ(θ)可以通過(guò)式(2)計(jì)算；Fτ為fτ的積累函數(shù)。

1.2 流形上的曲線擬合

有了2個(gè)功率譜之間的距離和測(cè)地線，我們就可以在流形上對(duì)功率譜進(jìn)行曲線擬合。假設(shè)給定一個(gè)功率譜序列：G={gτi(θ):θ∈[-π,π],i=0,1,…,n}，其中τi是一個(gè)指代時(shí)間的歸一化序列，且τ0=0,τn=1。這些功率譜通?？梢酝ㄟ^(guò)對(duì)觀測(cè)數(shù)據(jù)進(jìn)行短時(shí)傅里葉變換(STFT)得到，τi(i=0,1,…,n)可看做STFT中對(duì)應(yīng)時(shí)間窗的中點(diǎn)。令d代表兩功率譜之間的距離，則曲線擬合的過(guò)程就是尋找一條測(cè)地線fτ,τ∈[0,1]，使得下式取得最小值：

(7)

利用運(yùn)輸距離及其相應(yīng)的測(cè)地線來(lái)實(shí)現(xiàn)該過(guò)程，則整個(gè)測(cè)地線擬合的目標(biāo)就是要確定一條測(cè)地線fτ,τ∈[0,1]，使得下式取得最小值：

(8)

任意一條測(cè)地線都可以完全由兩“點(diǎn)”確定，在這里就是指f0和f1。另外，它也可以根據(jù)式(2)和式(6)中的Ψ確定。f0,f1和Ψ是根據(jù)STFT處理后得到的功率譜G來(lái)確定的，根據(jù)式(3),有：

(9)

根據(jù)式(6)，由ψ(θ)可得f0和f1之間測(cè)地線fτ的具體表達(dá)式為：

F0(θ)=F1(ψ(θ))=Fτi((1-τi)θ+τiψ(θ))

(10)

(11)

將式(11)代入式(9)，則目標(biāo)函數(shù)可寫(xiě)為：

(12)

(13)

(14)

所以Jg(f0,f1)可以通過(guò)以下的求和式進(jìn)行近似：

(15)

所以，功率譜的測(cè)地線擬合就轉(zhuǎn)化為了以式(15)為目標(biāo)函數(shù)的最優(yōu)化問(wèn)題，其線性約束為：

(16)

圖1 功率譜測(cè)地線擬合示意圖

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

本節(jié)采用加拿大McMaster大學(xué)提供的IPIX雷達(dá)#320號(hào)數(shù)據(jù)[14]對(duì)本文所提方法的檢測(cè)性能進(jìn)行驗(yàn)證，該數(shù)據(jù)文件包含14個(gè)距離單元，每個(gè)單元由131 072個(gè)采樣樣本構(gòu)成，其中第7個(gè)單元為主目標(biāo)單元，第6、8、9單元為次目標(biāo)單元，其余為純雜波單元。

首先將第1個(gè)距離單元的采樣數(shù)據(jù)作為雜波背景，通過(guò)STFT分析海雜波本身的譜特性。取觀測(cè)數(shù)據(jù)樣本長(zhǎng)度為256，窗函數(shù)寬度為32，時(shí)間窗重疊區(qū)域?qū)挾葹?6，脈沖重復(fù)頻率fr=1 000 Hz，可得第1個(gè)距離單元的海雜波二維譜圖，如圖2所示。

圖2 第1個(gè)距離單元的海雜波二維譜圖

從圖2中可以看出，海雜波譜是非對(duì)稱(chēng)的，由于海浪的移動(dòng)，使得譜的主瓣偏離零頻，雜波能量主要集中在[-200 Hz,0 Hz]這個(gè)區(qū)間內(nèi)，下面將仿真目標(biāo)信號(hào)加到第1個(gè)單元的純海雜波數(shù)據(jù)中，設(shè)雷達(dá)在該距離單元觀測(cè)的復(fù)包絡(luò)信號(hào)為：

x(n)=s(n)+v(n),n=1,…,N

(17)

式中：s(n)=aej(2πfdn/fr)，a為信號(hào)幅度，fd為目標(biāo)多普勒頻率。

令信噪比分別取0 dB和4 dB,目標(biāo)多普勒頻率分別取-100 Hz(主雜波譜區(qū))和100 Hz(非主雜波譜區(qū)),分別給出相應(yīng)的STFT變換后的譜圖，如圖3所示。

圖3 海雜波加仿真信號(hào)后STFT三維譜圖

從圖3(a)和圖3(c)中可以看出，當(dāng)目標(biāo)頻率處于非主雜波區(qū)(fd=100 Hz)時(shí),可以非常輕易地從STFT變換后的時(shí)頻圖像中判斷目標(biāo)的存在。但當(dāng)目標(biāo)頻率處于主雜波區(qū)(fd=-100 Hz)時(shí)，如圖3(b)和圖3(d)所示，信噪比較高時(shí)，還可以從譜圖中對(duì)目標(biāo)進(jìn)行辨別；但在低信噪比(SNR=0 dB)時(shí)，目標(biāo)被雜波淹沒(méi)，受譜圖分辨力的限制，很難判斷目標(biāo)的存在。

然后就可以用Matlab中的最優(yōu)化工具箱對(duì)該問(wèn)題進(jìn)行求解。對(duì)應(yīng)圖2中的4種情況，擬合得到的4條測(cè)地線的2個(gè)端點(diǎn)如圖4所示。

圖4 擬合得到的測(cè)地線端點(diǎn)圖

圖5 擬合處理后的譜圖

從圖4中單個(gè)時(shí)間點(diǎn)的功率譜已經(jīng)可以初步判斷各個(gè)情況下目標(biāo)的存在，但當(dāng)目標(biāo)處于主雜波區(qū)時(shí)，還不夠明顯。根據(jù)式(11)，有了2個(gè)端點(diǎn)就可以確定一條測(cè)地線，進(jìn)而就得到了測(cè)地線上各個(gè)時(shí)刻的功率譜，如圖5所示。

從圖5中可以看出，當(dāng)目標(biāo)處于主雜波區(qū)(fd=-100 Hz)并且信噪比較低時(shí)(SNR=0 dB)，與圖3(d)相比，圖5(d)中可以明顯地看到一條凸出的“山脊”，其所在頻率位置與目標(biāo)多普勒頻率相近，進(jìn)而就可以較容易地判斷目標(biāo)的存在。

以上結(jié)果都是將仿真目標(biāo)信號(hào)加到純雜波數(shù)據(jù)上進(jìn)行處理的，下面直接對(duì)實(shí)際有目標(biāo)的第7個(gè)距離單元數(shù)據(jù)進(jìn)行STFT變換，并進(jìn)一步進(jìn)行測(cè)地線擬合，給出處理結(jié)果如圖6所示。

圖6 第7個(gè)單元(主目標(biāo)單元)STFT+測(cè)地線擬合處理后的譜圖

IPIX實(shí)測(cè)數(shù)據(jù)是由雷達(dá)對(duì)漂浮于海面直徑約1 m的金屬球持續(xù)照射得到的，該金屬球會(huì)隨著海浪的波動(dòng)進(jìn)行漂移，但其運(yùn)動(dòng)速度與海浪運(yùn)動(dòng)速度相比相對(duì)較小，因而體現(xiàn)在功率譜上即其能量主要集中在零頻附近，且其運(yùn)動(dòng)是無(wú)規(guī)律的，目標(biāo)RCS和速度都會(huì)隨著時(shí)間推移而發(fā)生變化，因而其功率譜也會(huì)隨時(shí)間變化。從圖5中可以非常明顯地看到，在零頻附近有一條凸出的脊線，而且這條脊線并不像加仿真目標(biāo)信號(hào)時(shí)形狀比較規(guī)則，而是在時(shí)域上有一定起伏和變化，這與實(shí)際目標(biāo)狀況是相符的，因而就可以判斷目標(biāo)的存在。

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

本文根據(jù)功率譜的信息幾何理論，提出利用流形上的距離和測(cè)地線對(duì)STFT變換后的不同時(shí)間點(diǎn)的功率譜實(shí)現(xiàn)曲線擬合，以提高其時(shí)頻分辨力；并通過(guò)推導(dǎo),將該擬合過(guò)程轉(zhuǎn)換為可用數(shù)值方法解決的具有線性約束的最優(yōu)化問(wèn)題，使其容易求解；最后在仿真實(shí)驗(yàn)中，分別利用IPIX實(shí)測(cè)純海雜波數(shù)據(jù)加仿真目標(biāo)信號(hào)和本身帶目標(biāo)的海雜波數(shù)據(jù)2種方式對(duì)所提方法的性能進(jìn)行了測(cè)試，驗(yàn)證了本文方法的有效性。