基于水聲信號的瞬時頻率估計算法

張 靖

(大連測控技術(shù)研究所,遼寧 大連 116013)

0 序言

在許多實際的信號處理應(yīng)用中,估計一個非平穩(wěn)時變過程的瞬時頻率(Instantaneous Frequency,IF)意義重大。一些信息探測系統(tǒng)只要系統(tǒng)與目標之間有相對運動,就會產(chǎn)生“多普勒效應(yīng)”,頻率會隨相對運動發(fā)生變化,傳播媒質(zhì)的擾動也會使頻率變化,雷達、聲吶、移動通信、醫(yī)療設(shè)備和天文觀測等都存在這一問題。

目前針對瞬時頻率的估計方法的核心為概論統(tǒng)計,統(tǒng)計中所用的統(tǒng)計量,主要是樣本平均值或樣本函數(shù)算數(shù)平均值,它們具有最小方差性、無偏性等相當好的特性。然而樣本平均值的抗干擾性很差。為了獲得穩(wěn)健統(tǒng)計解,發(fā)展了一種既保留所有數(shù)據(jù),又削弱污染源對統(tǒng)計結(jié)果影響的統(tǒng)計方法,即 M-估計[1],其能獲得穩(wěn)健性的瞬時頻率估計值。

目前對瞬時頻率的估計方法有很多,對于水聲信號受海洋聲場環(huán)境非常復(fù)雜、測量測備工作狀態(tài)不穩(wěn)定、信噪比不足等因素的影響,很難獲得理想的瞬時頻率特征[2-5]。對于隨時間呈線性變化的非平穩(wěn)信號來說,其穩(wěn)健的瞬時頻率估計是必要的。研究分析水聲信號的瞬時頻率估計方法對實船噪聲特征頻率提取,輻射噪聲“多普勒效應(yīng)”分析及噪聲源分析意義重大。本文提出一種綜合相位法和M估計方法獲得水聲信號的瞬時頻率,并利用計算機仿真分析研究該方法的效果。

1 相位差分方法

對于理想的信號,其離散水聲信號可以表示成r(n)=Aexp[jφ(n)]。

根據(jù)離散信息瞬時頻率定義,離散時間域?qū)?yīng)為差分運算,即可得到離散時間瞬時頻率的估計。有3種相位差分算子[6],它們分別是前向有限差分、后向有限差分和中心有限差分,與這3種相位差分算子相對應(yīng)的瞬時頻率估計見式(1)～式(3)。

當信噪比很高時,利用中心差分法直接估計某特征頻率的瞬時頻率,得到的頻移特征曲線較平滑,誤差較小,精度高。

由于水聲信號容易受到干擾,單純地采用時刻Ti-1和Ti+1的相位估計Ti時刻的瞬時頻率,由于相位離散度大引起的瞬時頻率估計誤差也很大,很難準確估計瞬時頻率。由于通過較長的平均時間可以有效地抑制噪聲的影響,因為當噪聲服從高斯分布時,噪聲方差將隨著平均時間的增加而線性減小。因此在低信噪比情況下,增加有效觀測時間是提高同步性能的方法之一,從而采用長時間多點相位中心差分法。

根據(jù)Tayker公式可得:

針對離散相位信息,分別取 Δ t=i,i+1;t=n0,則可得

因此可近似得到 1組 2f1(n0),2f2(n0),…,2fi(n0)。

2 基于M-估計方法分析瞬時頻率

對于M-估計方法,核心是在穩(wěn)健性和精確性的2個定量化泛函空間中,找出1個鞍點,即敏感性最差、精確性最好的解[6]。實際是取函數(shù) ψ(x),解方程

令wi=,則上式剛好寫成加權(quán)平均的形式t=∑xiwi,即形成自動加權(quán)的權(quán)函數(shù)平均形式。

通過求解方程可知,當滿足tn=tn-1時,tn-1為最優(yōu)解。

依據(jù)估計準則的統(tǒng)計背景,M-估計可以分為3類[7]:一是基于概率密度(或誤差分布)的M估計,如最小二乘估計等;二是與概率密度無關(guān)的M估計(經(jīng)驗方法),權(quán)函數(shù)通過數(shù)值模擬獲得,而無須顧及是否知道誤差的概率密度,如丹麥法等;三是基于隨機模型驗后估計的穩(wěn)健估計,如李德仁選權(quán)迭代法。在參數(shù)回歸估計中,不同的權(quán)函數(shù)形式產(chǎn)生了不同的估計方法。為了實現(xiàn)參數(shù)估計的穩(wěn)健性,一般遵循大誤差以小權(quán),小誤差以大權(quán),等價權(quán)公式可以通過數(shù)值模擬經(jīng)驗確定的原則。等價權(quán)是殘差的函數(shù),參數(shù)的迭代求解中殘差的改變引起等價權(quán)的改變。本文采用的權(quán)函數(shù)[8]為:

3 仿真分析

本文通過仿真模擬不同信噪比下水下目標勻速直線運動產(chǎn)生的單頻線譜輻射噪聲,分析其由于多普勒效應(yīng)產(chǎn)生的瞬時頻率特征。

仿真數(shù)據(jù):信噪比分別為6 dB和9 dB;信號頻率100 Hz;采樣時間60 s;目標速度6 m/s;目標最近距離30 m;其建立的仿真信號如圖1所示。

對仿真建立的信號分別采用不同長度的權(quán)函數(shù)進行分析,其結(jié)果如圖2和圖3所示。

通過上圖可以看出,對于同1組噪聲信號,權(quán)函數(shù)的長度越長,其瞬時頻率估計值越接近真值;對比不同信噪比的噪聲信號,采用相同權(quán)函數(shù)長度分析情況下,信噪比越高,瞬時頻率估計值越理想。同時應(yīng)該指出,該方法權(quán)函數(shù)長度越長,要求的噪聲數(shù)據(jù)量越大,因此計算量大耗時較長,一般適合對信號的后續(xù)分析,不適合對信號的實時跟蹤。

4 結(jié) 語

本文通過研究建立了一種穩(wěn)健的瞬時頻率估計方法,通過計算機仿真分析,研究了不同信噪比、不同權(quán)函數(shù)長度對瞬時頻率估計的影響。綜合分析表明,該方法能有效的估計運動目標速度較低情況下產(chǎn)生的瞬時頻率。

圖3 不同權(quán)函數(shù)長度對應(yīng)IF估計值與真值對比(SNR=6)Fig.3 The IF estimate contrast with the true value in different weight function length(SNR=6)

[1]王安幫,等.M-估計的權(quán)函數(shù)選取及其抗差性比較[J].北京測繪,2008,10(1):8-10.WANG An-bang,et al.Selection of M-estimated weight functions and robust comparision[J].Beijing Surveying and Mapping.2008,10(1):8-10.

[2]李紅娟,孫超.水聲通信自適應(yīng)多普勒頻移補償仿真研究[J].計算機仿真,2007,24(2):302-304.LI Hong-juan,SUN Chao.Adaptive doppler compensation forunderwateracousticcommunication[J].Computer Simulation,2007,24(2):302-304.

[3]李君橋,等.基于多普勒頻移的水下運動目標測距技術(shù)[J].聲學技術(shù),2008,27(5):758-760.LI Jun-qiao,et al.Passive ranging of underwater motional targets based on Doppler frequency shift[J].Technical Acoustics,2008,27(5):758-760.

[4]莊曉明,等.基于多普勒頻移的水下超聲探測目標回波信號處理方法[J].艦船電子對抗,2008,31(6):88-91.

[5]DJURIC P M,KAY S M.Parameter estimation of chirp signals[J].IEEE Transactions on Acoustics,Speech and Signal Pro2 cessing,1990,38(12):2118-2126.

[6]孫曉奇.非平穩(wěn)信號瞬時頻率估計[D].碩士論文,2009.

[7]王文軍.極大似然估計的穩(wěn)健改進[J].武漢工業(yè)學院學報,1990,(4):35-41.

[8]黃曉紅,等.改進的相位展開算法及其在瞬時頻率估計中的應(yīng)用[J].電子學報,2009,(10):2266-2272.