IIR濾波器零相位數字濾波及其應用

摘要： 本文介紹了一種利用Delphi7實(shí)現數字濾波器零相位濾波的方法，與一般差分濾波方法相比，該方法不僅可以避免相移，而且還能改善起始部分的波形畸變，這一點(diǎn)在信號處理中具有重要的應用價(jià)值。
關(guān)鍵詞： 虛擬儀器信號處理零相位數字濾波

Zerophase Digital Filter and Application
Fan XinhaiAn GangZhang Chuanqing

Abstract： This paper introduces a method to realize zerophase filtering for digital filters using Delphi7. Comparing with general difference filter, the zerophase digital filter can not only avoid phaseshifting, but also improve the wave aberration of section start. This method has important practical value in signal processing.
Keywords： virtual instrument, signal processing, zerophase digital filtering.

0引言
在動(dòng)態(tài)測試信號處理過(guò)程中，濾波器是常用的測試儀器之一。它常被用于抗混濾波，以避免傅立葉變換時(shí)在頻域產(chǎn)生混疊，或從具有多種頻率成分的復雜信號中，將感興趣的頻率成分提取出來(lái)，而將不感興趣的頻率成分衰減掉。在傳統測試儀器中，濾波器的功能通常需要依靠硬件系統來(lái)實(shí)現。
隨著(zhù)數字信號處理技術(shù)的不斷完善，計算機硬件技術(shù)的日新月異以及軟件技術(shù)飛速發(fā)展，測試儀器系統的設計思想發(fā)生了重大改變。部分傳統的專(zhuān)用測試設備會(huì )逐步被以計算機和應用軟件為核心的虛擬儀器所代替[1]。虛擬儀器的出現標志著(zhù)“軟件即儀器(The soft is the instrument)”時(shí)代的到來(lái)。在計算機輔助測試系統(CAT，Computer Aided Test)中，以往模擬濾波器(AF，Analog Filter)的功能，可用數字濾波器(DF，Digital Filter)來(lái)替代。數字濾波器的實(shí)現不但比模擬濾波器容易的多，而且還能獲得較理想的濾波器性能。

2數字濾波器的時(shí)域描述與分類(lèi)
對于一個(gè)線(xiàn)性移不變(LSI，Linear Shift Invariant)離散時(shí)間系統，如圖1所示，可用如下差分方程表示：y(n)+∑Nk=1a(k)y(n-k)=∑Mr=0b(r)x(n-r)(1)式中a(k),b(r)為方程系數。

圖1LSI系統
如果a(k)，k=1,2,Λ,N不全為零，該系統為無(wú)限沖擊響應(IIR)系統。若a(k)均為零，并令b(0)=1，則y(n)=∑Nr=1b(r)x(n-r)+x(n)(2)該系統為有限沖擊響應(FIR)系統。
由此，數字濾波器在從實(shí)現方法上，有IIR濾波器和FIR濾波器之分。這兩類(lèi)濾波器無(wú)論在性能上，還是在設計方法上都有著(zhù)很大的區別。IIR數字濾波器與FIR濾波器相比，前者保留了模擬濾波器的優(yōu)點(diǎn)，幅頻特性較好，但存在相位失真。后者相頻特性較好，可實(shí)現線(xiàn)性相位，但在相同指標要求下要比前者的階數高的多。

3差分數字濾波及其存在的問(wèn)題
IIR數字濾波器的設計方法主要有雙線(xiàn)性Z變換法和沖激響應不變法；FIR數字濾波器的設計方法主要有窗函數法、頻率抽樣法和切比雪夫逼近法等。有關(guān)數字濾波器設計的詳細方法見(jiàn)文獻[2]，這里就不再贅述。
對于設計好的濾波器，確定出濾波器的系數，通過(guò)對差分方程式(1)或式(2)的疊代算法，可以實(shí)現IIR濾波器或FIR濾波器的數字濾波?？紤]到程序的通用性，這里就按式(1)編寫(xiě)一個(gè)差分濾波函數DiffFilter(b,a,x)，調用參數為三個(gè)數組，b用于存放濾波器的系數b(r)(r=0,1,…,M)；a用于存放濾波器的系數a(k)(k=0,1,…,N)，且a(0)=1；x用于存放輸入信號，返回后用于存放濾波后的輸出信號。如果將調用參數a均賦零值，則可實(shí)現FIR濾波器的數字濾波。
在此，用一實(shí)例進(jìn)行說(shuō)明。圖2為某型坦克發(fā)動(dòng)機轉速1000rpm時(shí)，在排氣口測得的噪聲信號(采樣點(diǎn)數1024，采樣頻率4KHz)。從信號的波形可以看出，它具有一定的周期性，主要是以發(fā)動(dòng)機周期性排氣產(chǎn)生的噪聲為主。

圖2某型坦克發(fā)動(dòng)機噪聲信號
采用窗函數(Hanning窗)法設計出某200階FIR低通濾波器，截止頻率為200Hz，分析頻率為2KHz，其幅頻曲線(xiàn)及相頻曲線(xiàn)如圖3所示。

圖3低通濾波器的幅頻特性和相頻特性
利用該FIR低通數字濾波器對圖2所示發(fā)動(dòng)機噪聲信號進(jìn)行差分濾波，濾波后的波形如圖4所示。

圖4差分濾波后的發(fā)動(dòng)機噪聲信號
對比圖3和圖4所示波形可以看出，差分濾波后，信號變得“光滑”了許多，原來(lái)的“毛刺”被濾掉了。但同時(shí)也可以看出差分濾波方法存在兩個(gè)問(wèn)題：一是濾波后的信號相對與原信號而言，發(fā)生了相移；二是濾波后的信號在起始部分，波形畸變較為嚴重。在信號處理中，如果對信號的相位有特殊的要求，相移問(wèn)題需要引起高度的注意。而起始部分的畸變是由于疊代過(guò)程中，沒(méi)有考慮濾波器的初始條件，剛開(kāi)始點(diǎn)數少，沒(méi)能用到濾波器全部系數的緣故。大多情況下，這種畸變可以接受，但當數據較短，而濾波器的階數又較高時(shí)，這種畸變會(huì )帶來(lái)較大的負面影響。
如果輸入信號為一單位脈沖信號，即
x(n)=1,n=1
0,n≠1
那么，濾波器的輸出為其脈沖響應。本文所用200階FIR低通數字濾波器的脈沖響應如圖5所示。

圖5濾波器的脈沖響應3零相位數字濾波算法及其實(shí)現