《電子技術(shù)應(yīng)用》
您所在的位置：首頁 > 測試測量 > 設(shè)計應(yīng)用 > 滑動頻率估計算法在微波測量液位中的應(yīng)用
滑動頻率估計算法在微波測量液位中的應(yīng)用
來源:電子技術(shù)應(yīng)用2011年第1期
朱 磊,， 董 亮,， 苗鳳娟,， 劉樹東
齊齊哈爾大學(xué) 通信與電子工程學(xué)院, 黑龍江 齊齊哈爾161006
摘要： 在微波測量液位的系統(tǒng)中,，正弦信號的頻率估計是實現(xiàn)液位測量的關(guān)鍵,。比較并分析了離散FFT變換的兩種正弦信號頻率估計方法：Rife比值法和三角形法,提出了一種滑動頻率估計方法,，通過實時確定比例因子,，得到了滑動頻譜校正算法,。實驗結(jié)果表明,，該方法頻率估計精度較高，估計性能優(yōu)于單一的Rife法,、三角形法,，適合于低信噪比條件下的實時信號處理。
中圖分類號: TN911.23
文獻標(biāo)識碼: A
文章編號： 0258-7998(2011)01-0074-03
Sliding frequency estimation algorithm in the application of microwave liquid level measurement
Zhu Lei, Dong Liang, Miao Fengjuan, Liu Shudong
Institute of Communications and Electronics Engineering, University of Qiqihar, Qiqihar 161006, China
Abstract： In the microwave liquid level measurement system,the frequency estimation of sine signal is the key to achieve liquid level measurement.The article compares sinusoidal frequency estimation methods of the discrete FFT transform,Rife ratio method and triangle method and proposed a sliding frequency estimation method.Through real-timely determining the scale factor,the sliding spectrum correction algorithm has been achieved.Experimental results show the new algorithm has high estimation accuracy, its estimation performance is better than single Rife method or triangle method,and it is suitable for real-time signal processing under low SNR.
Key words : liquid level measurement; sliding frequency estimation; data length; accuracy


   液位是自動計量領(lǐng)域中重要的檢測與控制參數(shù)之一,，液位測量技術(shù)經(jīng)過不斷的發(fā)展,，已經(jīng)得到了很大的進步。其中微波法測量液位是近幾年發(fā)展迅速而且日益成熟的一種方法,，它采用非接觸的測量方法,，通過測量微波傳播延遲時間的原理來測量液位。在微波液位測量系統(tǒng)中,，液位測量的關(guān)鍵在于正弦信號的頻率估計[1],。正弦信號的頻率估計有很多方法，其中最簡單的方法是FFT變換,，但是由于FFT的柵欄效應(yīng)[2],，使得正弦信號的頻譜發(fā)生泄漏，導(dǎo)致頻率估計精度不高,，無法滿足高精度的測量要求,。目前提高頻率估計精度的校正方法有很多，Rife插值法[3]利用兩個采樣點的比值來估計峰值的位置,，但在噪聲背景中,，容易造成插值方向錯誤，引起較大的估計誤差,。參考文獻[4]提出的FFT+FT譜連續(xù)細化法所需的計算量較大,。參考文獻[5-6]提出的三角形法是根據(jù)幾何原理進行頻譜校正的,，該方法運算簡單，受信噪比變化影響較小,。參考文獻[7]提出的相位差法精度較高,，但在進行相位估計時會遇到相位模糊的問題。
　　基于上述分析可以發(fā)現(xiàn),，Rife比值法和三角形校正法在頻譜校正上存在著各自的問題,。本文在分析兩者優(yōu)缺點的基礎(chǔ)上，結(jié)合兩者的優(yōu)勢,，提出了一種魯棒性更強的滑動頻率估計算法,。該方法假設(shè)噪聲是加性高斯白噪聲，根據(jù)信噪比的變化情況,，利用單次FFT后峰值譜線周圍的多條譜線信息的三角形法和Rife比值法估計信號頻率。
1 算法分析
1.1 Rife比值法
　　Rife比值法是利用主瓣峰頂附近兩條譜線的幅度比值進行頻譜校正的,，信號的實際頻率與估計頻率之間的相對偏差δ為(加矩形窗)：


式中m為幅值最大值處的離散頻率索引值,，fs為采樣頻率，N為采樣點數(shù),。該方法的優(yōu)點是插值公式簡單,，且不考慮噪聲影響的情況下，F(xiàn)FT主瓣內(nèi)次大值的幅度永遠大于其他旁瓣的幅度,，因此插值不會出現(xiàn)方向錯誤,。但是在有噪聲的情況下，當(dāng)|?啄|比較小時,，可能會出現(xiàn)FFT頻譜中其他譜線的幅度超過主瓣內(nèi)次大值的情況,，從而造成頻率插值方向相反，引起較大的頻率估計誤差,。當(dāng)N一定時,，不同信噪比下存在相應(yīng)的臨界值?啄0。當(dāng)|&delta;|<&delta;0時,，Rife比值法的頻率估計誤差比較大,特別是在|&delta;|接近于0時,，頻率估計誤差將大大增加。
1.2 三角形法
　　三角形法主要是根據(jù)幾何原理進行頻譜校正,，即用直線分別連接主瓣內(nèi)峰值的左右譜線時, 可近似形成一個三角形,，通過三角形的比例關(guān)系，可以得到譜線號修正量的式子am,即:

    該方法校正原理簡單,，計算量小,，精度較高。選取不同的n值,，所獲得的測量精度也不同,。隨著信噪比的增加,，三角形法的測量精度有所增加，但提高幅度不大,，受信噪比變化的影響較小,。它不同于Rife插值法，整個頻段誤差的變化范圍較小,，穩(wěn)定性好,，但是頻率估計誤差比Rife插值法在|?啄|不接近0時大。
1.3 滑動頻率估計算法
    矩形窗Rife比值法和三角形法都只需要一次FFT的運算量,，但在頻率估計精度上都存在各自的問題,。因此，本文提出了滑動頻率估計算法,，將兩種算法有效地結(jié)合起來,，在有噪聲的情況下，根據(jù)比例因子,，即頻率所在位置的不同（&delta;的不同）,，自適應(yīng)地采取相應(yīng)的算法。算法的具體過程如下：
    (1)通過仿真實驗,，實時確定現(xiàn)場所在環(huán)境下的信噪比,，記Rife比值法的偏差閾值為&delta;0，并查找P(k)的最大值P(m),，其中P(k)=X2(k),。
    (2)由于三角形法的頻率估計誤差穩(wěn)定，因此估計的偏差am可信度高,，當(dāng)|am|<&delta;0,，使用Rife插值法估計效果較差，此時使用三角形法進行頻率估計,，按照公式(4)可得到信號的實際頻率估計,。
    (3)反之為|am|&ge;&delta;0時，三角形法對實際頻率的估計誤差大于矩形窗Rife插值法,，因此,，使用Rife插值法進行頻率估計，按照公式(2)可得到信號的實際頻率估計,。
    可見,，Rife比值法與三角形法的有效結(jié)合構(gòu)成了滑動頻率估計算法，根據(jù)次大值與最大值的幅度比值來預(yù)先設(shè)定偏差閾值&delta;0,。根據(jù)|&delta;|的不同,，即頻率所在位置的不同，采用不同的估計方法,。
2 仿真實驗
    為了驗證所提算法的有效性,，作如下仿真實驗,。
    這里采用鋸齒波LFMCW雷達測量液位，根據(jù)被測頻率與目標(biāo)距離之間的關(guān)系,，采用3種方法對實際距離進行測量,。設(shè)定液面與天線的距離為0~20 m，信噪比為12 dB,。雷達的調(diào)制周期Ts=1 ms,，調(diào)制帶寬B=600 MHz，fs=256 kHz,，N=256,。對所選的幾個被測距離分別采用Rife比值法、三角形法(n=3),、滑動頻率估計算法進行估值（&delta;0=0.12）,，各被測點進行了1 000次測量得到平均值R和均方根誤差&sigma;R如表1所示。


　從表1可以看出,，當(dāng)被測距離接近離散頻譜最大譜線所對應(yīng)的距離時（R=8.045）,，即|?啄|較小時，Rife插值法的距離測量均方根誤差大大增加,，為0.010 8。而此時的三角形法距離測量的均方根誤差較小,，為0.000 5,，且整個變化區(qū)間內(nèi)，均方根誤差變化較小,。但是當(dāng)被測距離遠離離散頻譜最大譜線所對應(yīng)的距離時,，三角形法的估計誤差比Rife插值法大（R=8.020）。而本文所提出的滑動頻率估計算法結(jié)合兩者的優(yōu)勢,，根據(jù)|&delta;|的不同選取不同的估計方法,，在整個所觀測的測量范圍內(nèi)，距離估計的均方根誤差較小,，測量精度較高,。
    下面來觀察一下距離估計誤差與數(shù)據(jù)輸入長度N(N=16，32,，64,，128，256,，512,，1 024)之間的關(guān)系。圖1為N不同時三種估計方法的性能比較,。設(shè)R=4.176 m,SNR=15 dB,，其他條件同上,。由仿真結(jié)果可以看出，當(dāng)N取值較小時,， Rife算法的頻率估計精度較低,，而三角形法的頻率估計精度較高，且在所觀察的變化區(qū)間內(nèi),，三角形法的頻率估計性能穩(wěn)定,，而Rife算法在N較大的時候，頻率估計效果較好,。在結(jié)合兩者優(yōu)點的基礎(chǔ)上,，本文提出的滑動頻率估計算法在所觀察的區(qū)間內(nèi)受數(shù)據(jù)長度的變化影響較小，測距精度較高,，而且運算簡單,，計算量不大。

    再來觀察一下不同估計方法隨信噪比的變化情況,。此時R=14.565 m,，其他條件同上。圖2是三種估計方法的距離估計的均方根誤差與信噪比的變化情況,。從圖中可以看出,，Rife插值法在整個信噪比區(qū)間均方根誤差變化較大，受信噪比的影響較大,，當(dāng)信噪比較低時,，Rife插值法引起了更大的估計誤差。而三角形法在整個信噪比區(qū)間變化較小,，即使在信噪比較低時,，其均方根誤差也不大，但在信噪比較高時,，估計誤差大于Rife插值算法,。而本文所采用的滑動頻率估計算法，克服了上述兩種算法的不足,，在所觀測的整個信噪比變化區(qū)間內(nèi),，距離測量的均方根誤差都較小，測量精度高,。

    本文通過對Rife比值法和三角形法的研究,，建立了一種隨信噪比變化的頻譜校正模型。根據(jù)信噪比的變化情況,，利用正弦信號采樣序列的FFT頻譜圖上譜峰附近的多條譜線信息,，結(jié)合Rife比值法和三角形法各自的優(yōu)點，提出了一種滑動頻率估計方法,。該方法有效地融合了多條譜線信息的貢獻,，根據(jù)頻率所在位置的不同選取不同的估計方法,，避免了Rife算法的插值方向錯誤問題。實驗結(jié)果顯示,，該方法的測量精度高,，即使在噪聲背景中，微波液位測量的誤差也能滿足工程需要,，并且頻譜校正精度受數(shù)據(jù)長度的變化影響較小,。
參考文獻
[1] 陳先忠,柳瑾.FMCW微波液位測量的快速頻率估計算法[J].傳感技術(shù)學(xué)報，2005,18(4):901-905.
[2] 朱磊,董亮,劉樹東.基于Quinn 算法與改進的Rife 算法的正弦信號頻率估計[J].大慶石油學(xué)院學(xué)報,2010,34(1): 98-101.
[3] RIFE D C, VINCENT G A. Use of the discrete Fourier  transform in the measurement of frequencies and levels of  tones[J]. Bell. Sys. Tech. J., 1970,49(2):197-228.
[4] 劉進明,應(yīng)懷樵. FFT譜連續(xù)細化分析的傅里葉變換法[J].振動工程學(xué)報,1995,8(2):162-166.
[5] 曹延偉,張昆帆,江志紅,等.一種穩(wěn)健的離散頻譜校正方法[J].電子與信息學(xué)報,2005,27(9):1353-1356.
[6] 張昆帆,王蘭云,趙擁軍.基于窗函數(shù)的離散頻譜校正法[J].現(xiàn)代雷達,2007,29(9):59-62.
[7] 齊國清,賈欣樂.基于DFT相位的正弦波頻率和初相的高精度估計方法[J].電子學(xué)報,2001,29(9):1164-1167.

此內(nèi)容為AET網(wǎng)站原創(chuàng),，未經(jīng)授權(quán)禁止轉(zhuǎn)載,。