朱 磊，董 亮，孫振龍，劉樹東

（齊齊哈爾大學(xué) a.通信與電子工程學(xué)院；b.招生辦，黑龍江 齊齊哈爾 161006）

0 引言

正弦信號(hào)的頻率估計(jì)是信號(hào)處理領(lǐng)域中的一個(gè)非常重要的問題，而FFT算法的提出可快速有效地實(shí)現(xiàn)正弦信號(hào)的頻率估計(jì)問題，但是由于FFT的頻譜泄露和柵欄效應(yīng)，使得直接由FFT來估計(jì)信號(hào)的頻率效果不是很理想，因此必須進(jìn)行頻譜校正。1975年，John C.Burges采用插值法對加矩形窗的離散頻譜進(jìn)行校正[1]。1983年，Thomas Grandke提出了加漢明窗的內(nèi)插校正法[2]。國內(nèi)一些學(xué)者也提出了許多離散頻譜的校正方法：如相似最大似然頻率估計(jì)法[3]、比值法[4-8]、相位差法[9-10]、能量重心法[11]、FFT+DFT細(xì)化法等等。這些方法實(shí)現(xiàn)過程復(fù)雜，計(jì)算量大，不易于滿足DSP的實(shí)時(shí)信號(hào)處理，特別在有噪聲時(shí)，估計(jì)性能明顯下降。

針對上述情況，筆者提出了一種離散頻譜梯形窗幅度比值校正法。通過分析梯形窗函數(shù)的信號(hào)模型了解到，通過選取適當(dāng)?shù)膬蛇厓A斜度，可以獲得較好的幅頻特性。因此，本文在分析梯形窗時(shí)域、頻域特性的基礎(chǔ)上，利用離散頻譜中最高譜線和次高譜線兩條譜線的幅度比值關(guān)系，推導(dǎo)了加梯形窗的幅度比值頻譜校正公式，采用仿真數(shù)據(jù)對該算法的頻率校正精度進(jìn)行了分析，并與加矩形窗時(shí)的頻譜校正精度進(jìn)行了比較，誤差分析和實(shí)驗(yàn)結(jié)果表明該算法在信噪比較低時(shí)，估計(jì)精度高，性能明顯優(yōu)于加矩形窗的情況。

1 梯形窗函數(shù)模型

梯形窗函數(shù)的時(shí)域波形如圖1所示。具體函數(shù)形式為

相應(yīng)的幅值譜為圖2所示，具體的表達(dá)式為

由圖可見，梯形窗的主瓣寬度比較寬，絕對值最小的兩個(gè)零點(diǎn)之間的距離為8π/τ+τ1，但是旁瓣電平比較低。表1為幾個(gè)常見的窗函數(shù)幅頻特性對比表，由表1可以知道，梯形窗函數(shù)是一種性能較好的窗函數(shù)，其幅頻特性的主瓣寬度雖然比矩形窗函數(shù)幅頻特性的主瓣寬度稍寬，但是其旁瓣衰減速度大于矩形窗的衰減速度，旁瓣電平低。而且，通過選取適當(dāng)?shù)摩?可以獲得無負(fù)旁瓣，可以抵制變換中的高頻干擾和泄漏，特別適合于有噪聲的場合。

表1 幾種窗函數(shù)的幅頻特性

2 離散頻譜梯形窗幅度比值校正法

下面來推導(dǎo)加梯形窗的離散頻譜校正公式。設(shè)單頻正弦信號(hào)x(t)的形式為

式中：A,α,φ分別為信號(hào)的頻率、幅度、初相位?，F(xiàn)在來考察這個(gè)信號(hào)在[0，T]上加梯形窗之后的頻譜

進(jìn)一步整理可簡化為

這里的W(ω)就是式（2）的梯形窗函數(shù)的窗譜。在這里忽略負(fù)頻率泄露的影響，只考慮ω＞0時(shí)的情況。這樣有

如果信號(hào)的真實(shí)頻率并不是正好落在離散頻譜的整數(shù)根上，如圖2所示。這時(shí)利用FFT變換所得到的頻率估計(jì)值就會(huì)有誤差，就需要進(jìn)行頻譜校正。設(shè)真實(shí)頻率與FFT得到的最高幅值ab（譜線序號(hào)為k）之間的偏差為δ，幅值為X（k），次高譜線cd在k+1處，幅值為X（k+1），那么次高峰與真實(shí)頻率之間的偏差為1-δ。則

可以解出

式中：δ∈[-0.5,0.5]。

由式（8）可見，梯形窗的頻譜校正是利用了次高峰的峰值與最高峰的峰值之間的關(guān)系得到的。如果次高峰在最高峰的左邊，那么校正頻率為α?=(k-δ)Δf；如果次高峰在最高峰的右邊，那么校正頻率為α?=(k+δ)Δf。綜上所述，信號(hào)的實(shí)際頻率估計(jì)

式（9）就是離散頻譜梯形窗幅度比值法的校正公式。

通過上述分析，可以得到新算法的實(shí)現(xiàn)步驟如下：

1）令 X（m）=FFT（x（n），其中 x（n）是 x（t）的采樣序列。

2）查找abs（X（m））的最大值X（k），次大值X（k-1）或X（k+1）。

4）按照公式（9）估計(jì)信號(hào)的頻率。

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

根據(jù)以上分析進(jìn)行仿真實(shí)驗(yàn)。設(shè)信號(hào)形式如式（3）所示。信號(hào)的幅值為A=1，φ=0，采樣點(diǎn)數(shù)N=512，采樣頻率fs=512 Hz。加矩形窗后的信號(hào)的頻率校正誤差和本文所提的梯形窗幅度比值法的頻率校正誤差的對比曲線如圖3所示。可以看出，無噪聲時(shí)，矩形窗的校正精度要比梯形窗的校正精度高，但是當(dāng)信號(hào)質(zhì)量明顯下降時(shí)，即噪聲比較大時(shí)，對本文所提的離散頻譜梯形窗的校正精度影響較小，而對加矩形窗的校正精度影響較大。這說明梯形窗的抑制噪聲的能力較強(qiáng)。

圖4給出了兩種方法的頻率估計(jì)的均方根誤差RMSE與SNR的變化關(guān)系，此時(shí)，α=102.125 Hz。從實(shí)驗(yàn)結(jié)果可見，無論是加矩形窗還是加梯形窗，這兩種算法的校正精度都會(huì)隨著信噪比的升高而逐漸趨近于真實(shí)頻率。這是因?yàn)檫@兩種算法都是從窗函數(shù)的數(shù)學(xué)定義出發(fā)來推導(dǎo)頻率的校正公式，因此校正頻率的精度較高。但同時(shí)也能看出，在信噪比較低時(shí)，梯形窗的校正效果要高于矩形窗；而信噪比高時(shí)，情況正好相反。這是因?yàn)樘菪未暗呐园觌娖奖容^低，可以有效地抑制干擾，但主瓣寬度比矩形窗寬，因此在信噪比高時(shí)，其頻率分辨率不如加矩形窗時(shí)的頻率分辨率高。

4 結(jié)論

筆者從梯形窗的數(shù)學(xué)定義出發(fā)，給出了一種加梯形窗的離散頻譜幅度比值校正法，利用最高譜線和次高譜線兩條譜線的幅度比值關(guān)系，推導(dǎo)了加梯形窗的離散頻譜校正公式，對加矩形窗和加梯形窗時(shí)的單頻正弦信號(hào)的頻率校正誤差進(jìn)行了比較。仿真結(jié)果表明，該頻譜校正算法的校正精度在低信噪比時(shí)，明顯優(yōu)于加矩形窗的校正精度，而且運(yùn)算簡單，適合于低信噪比條件下的實(shí)時(shí)信號(hào)處理等場合。

[1]BURGESS T.On digital spectrum analysis of peri odic signals[J].J.A coust Soc.Am，1975，58（3）：556-567.

[2]GRANDKE T.Interpolation algorithms for discrete fourier transforms of weighted signals[J].IEEE Trans.Instrumentation and Measurement，1983，32（2）：350-355.

[3]黃柏圣，陳晨，許家棟.一種新的復(fù)正弦信號(hào)二維頻率估計(jì)算法[J].計(jì)算機(jī)仿真，2009，26（3）：198-201.

[4]張強(qiáng)，張頻，張明童.加三角窗的頻譜校正[J].振動(dòng)與沖擊，2009，28（2）：96-98.

[5]齊國清，賈欣樂.插值FFT估計(jì)正弦信號(hào)的精度分析[J].電子學(xué)報(bào)，2004，32（4）：625-628.

[6]齊國清.幾種基于FFT的頻率估計(jì)方法精度分析[J].振動(dòng)工程學(xué)報(bào)，2006，19（1）：86-91.

[7]謝明，丁康.頻譜分析的校正方法[J].振動(dòng)工程學(xué)報(bào)，1994，7（2）：172-179.

[8]丁康，謝明.離散頻譜三點(diǎn)卷積幅值校正法的誤差分析[J].振動(dòng)工程學(xué)報(bào)，1996，9（1）：92-98.

[9]齊國清，賈欣樂.基于DFT相位的正弦波頻率和初相的高精度估計(jì)方法[J].電子學(xué)報(bào)，2001，29（9）：1164-1167.

[10]黃云志，徐科軍.基于相位差法的頻譜校正方法的研究[J].振動(dòng)與沖擊，2005，24（2）：77-79.

[11]丁康，江利旗.離散頻譜的能量重心校正法[J].振動(dòng)工程學(xué)報(bào)，2001，14（3）：354-358.