隨著科學(xué)技術(shù)的發(fā)展，頻率估計(jì)在在液位測(cè)量器中的應(yīng)用。是自動(dòng)計(jì)量領(lǐng)域中重要的檢測(cè)與控制參數(shù)之一，液位測(cè)量技術(shù)經(jīng)過不斷的發(fā)展，已經(jīng)得到了很大的進(jìn)步。目前提高頻率估計(jì)的校正方法有很多，Rife插值法[3]利用兩個(gè)采樣點(diǎn)的比值來(lái)估計(jì)峰值的位置，但在噪聲背景中，容易造成插值方向錯(cuò)誤，引起較大的估計(jì)誤差。參考文獻(xiàn)[4]提出的FFT+FT譜連續(xù)細(xì)化法所需的計(jì)算量較大?；谏鲜龇治隹梢园l(fā)現(xiàn)，Rife比值法和三角形校正法在頻譜校正上存在著各自的問題。本文在分析兩者優(yōu)缺點(diǎn)的基礎(chǔ)上，結(jié)合兩者的優(yōu)勢(shì)，提出了一種魯棒性更強(qiáng)的滑動(dòng)頻率估計(jì)算法。該方法假設(shè)噪聲是加性高斯白噪聲，根據(jù)信噪比的變化情況，利用單次FFT后峰值譜線周圍的多條譜線信息的三角形法和Rife比值法估計(jì)信號(hào)頻率。

　　1 算法分析

　　1.1 Rife比值法

　　Rife比值法是利用主瓣峰頂附近兩條譜線的幅度比值進(jìn)行頻譜校正的，信號(hào)的實(shí)際頻率與估計(jì)頻率之間的相對(duì)偏差δ為（加矩形窗）：

　　RIFE為開發(fā)和設(shè)計(jì)Web應(yīng)用程序提供了另外一種途徑。它基于Java平臺(tái)構(gòu)建，但它提供了所有需要的工具和AIPs用統(tǒng)一的風(fēng)格，快速地實(shí)現(xiàn)所有 web程序各開發(fā)角色之間的任務(wù)聯(lián)系。RIFE被設(shè)計(jì)成在整個(gè)開發(fā)周期中完全分離了開發(fā)角色之間的任務(wù)，每個(gè)developer,dba和 designer只需關(guān)注它們自己的任務(wù)，整個(gè)開發(fā)工作隨時(shí)都可以很容易的與其它團(tuán)隊(duì)的工作相結(jié)合，使得這樣子開發(fā)出來(lái)的將是一個(gè)松耦合的健壯的系統(tǒng)。 但是在有噪聲的情況下，當(dāng)|?啄|比較小時(shí)，可能會(huì)出現(xiàn)FFT頻譜中其他譜線的幅度超過主瓣內(nèi)次大值的情況，從而造成頻率插值方向相反，引起較大的頻率估計(jì)誤差。當(dāng)N一定時(shí)，不同信噪比下存在相應(yīng)的臨界值？啄0。當(dāng)|δ|<δ0時(shí)，Rife比值法的頻率估計(jì)誤差比較大，特別是在|δ|接近于0時(shí)，頻率估計(jì)誤差將大大增加。

　　1.2 三角形法

　　三角形法（triangle method）是一種儲(chǔ)量計(jì)算法。它的實(shí)質(zhì)是把形狀不規(guī)則的礦體，人為地簡(jiǎn)化為許多便于計(jì)算體積的三棱柱狀體。即在儲(chǔ)量平面圖所圈定的礦體范圍內(nèi)，以直線連接各相鄰勘探工程，把礦體分為一系列緊密連接的三角形塊段。再依據(jù)三角形塊段頂點(diǎn)的勘探工程資料，分別計(jì)算各地段的礦產(chǎn)儲(chǔ)量。這種儲(chǔ)量計(jì)算法不僅不能反映礦體的真實(shí)特點(diǎn)，而且計(jì)算過程繁雜，實(shí)際工作中很少應(yīng)用。

　　三角形法主要是根據(jù)幾何原理進(jìn)行頻譜校正，即用直線分別連接主瓣內(nèi)峰值的左右譜線時(shí)， 可近似形成一個(gè)三角形，通過三角形的比例關(guān)系，可以得到譜線號(hào)修正量的式子am,即：

　　該方法校正原理簡(jiǎn)單，計(jì)算量小，較高。選取不同的n值，所獲得的測(cè)量也不同。隨著信噪比的增加，三角形法的測(cè)量有所增加，但提高幅度不大，受信噪比變化的影響較小。它不同于Rife插值法，整個(gè)頻段誤差的變化范圍較小，穩(wěn)定性好，但是頻率估計(jì)誤差比Rife插值法在|?啄|不接近0時(shí)大。

　　1.3 滑動(dòng)頻率估計(jì)算法

　　矩形窗Rife比值法和三角形法都只需要FFT的運(yùn)算量，但在頻率估計(jì)上都存在各自的問題。算法的具體過程如下：

　?。?）通過仿真實(shí)驗(yàn)，實(shí)時(shí)確定現(xiàn)場(chǎng)所在環(huán)境下的信噪比，記Rife比值法的偏差閾值為δ0，并查找P（k）的值P（m），其中P（k）=X2（k）。

　　（2）由于三角形法的頻率估計(jì)誤差穩(wěn)定，因此估計(jì)的偏差am可信度高，當(dāng)|am|<δ0，使用Rife插值法估計(jì)效果較差，此時(shí)使用三角形法進(jìn)行頻率估計(jì)，按照公式（4）可得到信號(hào)的實(shí)際頻率估計(jì)。

　?。?）反之為|am|≥δ0時(shí)，三角形法對(duì)實(shí)際頻率的估計(jì)誤差大于矩形窗Rife插值法，因此，使用Rife插值法進(jìn)行頻率估計(jì)，按照公式（2）可得到信號(hào)的實(shí)際頻率估計(jì)。

　　可見，Rife比值法與三角形法的有效結(jié)合構(gòu)成了滑動(dòng)頻率估計(jì)算法，根據(jù)次大值與值的幅度比值來(lái)預(yù)先設(shè)定偏差閾值δ0。根據(jù)|δ|的不同，即頻率所在位置的不同，采用不同的估計(jì)方法。

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

　　為了驗(yàn)證所提算法的有效性，作如下仿真實(shí)驗(yàn)。

　　這里采用鋸齒波LFMCW雷達(dá)測(cè)量液位，根據(jù)被測(cè)頻率與目標(biāo)距離之間的關(guān)系，采用3種方法對(duì)實(shí)際距離進(jìn)行測(cè)量。設(shè)定液面與天線的距離為0~20 m，信噪比為12 dB。對(duì)所選的幾個(gè)被測(cè)距離分別采用Rife比值法、三角形法（n=3）、滑動(dòng)頻率估計(jì)算法進(jìn)行估值（δ0=0.12），各被測(cè)點(diǎn)進(jìn)行了1 000次測(cè)量得到平均值R和均方根誤差σR如表1所示。

　　從表1可以看出，當(dāng)被測(cè)距離接近離散頻譜譜線所對(duì)應(yīng)的距離時(shí)（R=8.045），即|?啄|較小時(shí)，Rife插值法的距離測(cè)量均方根誤差大大增加，為0.010 8。而此時(shí)的三角形法距離測(cè)量的均方根誤差較小，為0.000 5，且整個(gè)變化區(qū)間內(nèi)，均方根誤差變化較小。但是當(dāng)被測(cè)距離遠(yuǎn)離離散頻譜譜線所對(duì)應(yīng)的距離時(shí)，三角形法的估計(jì)誤差比Rife插值法大（R=8.020）。而本文所提出的滑動(dòng)頻率估計(jì)算法結(jié)合兩者的優(yōu)勢(shì)，根據(jù)|δ|的不同選取不同的估計(jì)方法，在整個(gè)所觀測(cè)的測(cè)量范圍內(nèi)，距離估計(jì)的均方根誤差較小，測(cè)量較高。

　　下面來(lái)觀察一下距離估計(jì)誤差與數(shù)據(jù)輸入長(zhǎng)度N（N=16，32，64，128，256，512，1 024）之間的關(guān)系。圖1為N不同時(shí)三種估計(jì)方法的性能比較。設(shè)R=4.176 m,SNR=15 dB，其他條件同上。由仿真結(jié)果可以看出，當(dāng)N取值較小時(shí)， Rife算法的頻率估計(jì)較低，而三角形法的頻率估計(jì)較高，且在所觀察的變化區(qū)間內(nèi)，三角形法的頻率估計(jì)性能穩(wěn)定，而Rife算法在N較大的時(shí)候，頻率估計(jì)效果較好。在結(jié)合兩者優(yōu)點(diǎn)的基礎(chǔ)上，本文提出的滑動(dòng)頻率估計(jì)算法在所觀察的區(qū)間內(nèi)受數(shù)據(jù)長(zhǎng)度的變化影響較小，測(cè)距較高，而且運(yùn)算簡(jiǎn)單，計(jì)算量不大。

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

　　本文通過對(duì)Rife比值法和三角形法的研究，建立了一種隨信噪比變化的頻譜校正模型。該方法有效地融合了多條譜線信息的貢獻(xiàn)，根據(jù)頻率所在位置的不同選取不同的估計(jì)方法，避免了Rife算法的插值方向錯(cuò)誤問題。實(shí)驗(yàn)結(jié)果顯示，該方法的測(cè)量高，即使在噪聲背景中，微波液位測(cè)量的誤差也能滿足工程需要，并且頻譜校正受數(shù)據(jù)長(zhǎng)度的變化影響較小。