IIR數字濾波器的設計與實現

張海蓉

（貴州師范學院，貴州 貴陽 550018）

IIR數字濾波器的設計與實現

張海蓉

（貴州師范學院，貴州 貴陽 550018）

數字濾波是數字信號處理的關鍵內容，可分為IIR和FIR兩大類。對于IIR數字濾波器的設計，需要借助模擬原型濾波器，再將模擬濾波器轉化為數字濾波器，文中采用的設計方法是脈沖響應不變法、雙線性變換法。根據IIR濾波器的特點，在MATLAB壞境下分別用雙線性變換法和脈沖響應不變法設計IIR低通數字濾波器，并比較兩種方法的優缺點。

IIR濾波器；MATLAB；數字濾波器；雙線性；脈沖響應

引言

數字濾波在諸多方面都涉及到，譬如通信、圖像編碼、語音編碼、雷達等領域中都有廣泛的應用。現今，它的應用取得了令人矚目的進展和成就。當然，相信它的未來也前途無量。

在文中用到的MATLAB是一種應用于科學計算領域的數學軟件，它主要包括數值計算和符號計算功能、繪圖功能編程功能以及應用工具箱的擴張功能。MATLAB推出的工具箱使各個領域的研究人員可以直觀方便地進行科學研究、工程應用，其中的信號處理、圖像處理、小波等工具箱為數字濾波研究的蓬勃發展提供了有力的工具。數字濾波器與模擬濾波器比較，它具有精度高、穩定、體積小、重量輕、靈活、不要求阻抗匹配和能實現模擬濾波器無法進行的特殊濾波等優勢[1]。本文主要介紹無限沖激響應數字濾波器（IIR）的設計原理、方法、步驟以及在MATLAB中的實現。

1 IIR數字濾波器的設計與基本結構

1.1 數字濾波器的定義及結構原理

數字濾波器利用離散時間系統對數字信號中時域離散、幅度量化的信號進行濾波就是數字濾波器的構成，具有完成信號濾波處理的功能，數字濾波器具備精度高、可靠性好、便于集成等優點[2]。數字濾波器可分為有限沖激響應（FIR）數字濾波器和無限沖激響應（IIR）數字濾波器[3]。IIR濾波器運算結構通常由延時、乘以系數和相加等基本運算組成，可以組合成直接型、正準型、級聯型、并聯型四種結構形式，都具有反饋回路。由于運算中的舍入處理，使誤差不斷累積，有時會產生微弱的寄生振蕩。

1.2 IIR濾波器的設計

濾波器就是依照所給定的指標來確定系數和階數。得到模擬濾波器指標是數字濾波器指標的轉換的，實現的主要方法有有兩種，一種是雙線性變換法，另一種是脈沖響應不變法，如圖3所示。

圖1 IIR數字濾波器的設計過程

1.3 IIR濾波器實現的基本結構

（1）IIR濾波器的直接型結構，如圖2所示:

圖2 直接型

優點：簡單、直觀，所用的延時器數量少，可以節省存儲單元或者寄存器。

缺點：作用于所用極點的重要因素就是改變某一個系數{ak}，同樣作用于所用的零點的關鍵因素就是改變某一個系數{bk}，更甚者有限字長的效應更易出現對外界事物反應很快的現象，針對這種結構，也輕易出現產生較大誤差和狀況的不穩固。大體上不汲取直接型結構，而是采取級聯型、并聯型等其他型的出現在三階以上的數字濾波器上。

（2）IIR濾波器的級聯型結構，如圖3所示:

圖3 級聯結構

優點：第一是所有系數變更僅作用子系統的零點和極點，易于獲取濾波器的零點和極點，對頻率特性也容易調改。第二是需要的存儲單元不用很多，在硬件實現過程中，用二階進行時分復用比直接型低對有限字長的影響。

缺點：當用二進制表示零、極點時，須采用有限位字長，對于各種不同實現方案，所帶來的誤差是有差別的，故需要優化各二階級聯的次序和極點，而且還不能直接調整零點。

（3）IIR濾波器的并聯型結構，如圖4所示：

圖4 并聯結構

優點：運算的速度快、每個基本節的偏差彼此不受影響，所以不僅擁有較小的誤差，還可調整極點的所在位置，其次傳輸準確的零點。

缺點：由于在子系統中的零點不是整個體系的零點，所以不能和級聯型一樣直接調動濾波器零點的所在位置。

2 IIR濾波器的實現與仿真

2.1 脈沖響應不變法設計

（1）基本原理

脈沖響應不變法是一種將模擬濾波器轉換成數字濾波器的基本方法。它通過對模擬濾波器的單位沖擊響應h（t）等間隔抽樣來獲得對應數字濾波器的單位脈沖響應h[k]，即

h[k]=h（t）|t=kT=h（kT）

（1）

在上述式子中，T作為抽樣間隔。若是知曉濾波器中模擬的系統函數H（s），那就可以轉換成數字的系統函數H（z），轉換步驟如5所示。

即先對H（s）進行Laplace逆變換獲得h（t）之后，再對h（t）等間隔抽樣得到h[k]，最后計算h[k]的z變換得到H（z）。

圖5 利用脈沖響應不變法將H（s）變換為H（z）

（2）設計方法

脈沖響應不變法常用的數字指標有{Ωp}、Ap、{Ωs}、As，設計的步驟如下：

①模擬濾波器的技術指標的被轉換是建立在數字濾波器的技術指標的基礎上。利用數字和模擬兩者頻率的間關系

（2）

然后將{Ωk}數字濾波器的頻率指標轉換為{ωk}模擬濾波器的頻率指標。

②在模擬和數字濾波器的衰減響應的坐標變換中，縱坐標衰減幅度沒有影響，無需變換，所以只需變換橫坐標的頻率軸，轉換后模擬濾波器的阻帶衰減As和通帶衰減Ap與數字濾波器相同。

③依照模擬濾波器的技術要求，設計其H（s）即系統函數，然后利用H（s）求解數字濾波器的系統函數H（z），并對所設計的數字濾波器的技術指標進行檢驗是否滿足要求[4]。

利用脈沖響應不變法及BW型低通濾波器設計一個數字低通濾波器滿足下列要求的，Ωp=0.2πrad，Ωs=0.6πrad，Ap=2dB,As=25dB

其轉換關系如下：

[numa,dena]=butter（N,wc,'s'）;%確定BW AF

[numd,dend]=impinvar（numa,dena,Fs）;%確定DF

其運行結果是Ap=1.0032dB,As=32.2549 dB。顯然，阻帶衰減不滿足設計要求，其原因是巴特沃斯型模擬濾波器不是帶限的，脈沖響應不變法所設計的數字濾波器存在頻譜混疊。其增益圖如圖6所示。

圖6 數字濾波器的增益響應

圖6是所設計的低通濾波器的增益響應，阻帶衰減不滿足所求要求時，其原因是原濾波器不是帶限的，被脈沖響應不變法所設計的數字濾波器呈現出頻譜混疊現象。又由于T的選擇對最終的設計結果沒有影響，所以不能通過減小T來減少頻譜混疊，有效的方法是待設計的阻帶最小衰減As。

2.2 雙線性不變法設計

（1）變換原理

雙線性變換法的基本原理是：把s域Ω值從0到±∞的分布壓縮至0到±π之間，這樣變換到z域就不會出現周期重復，就可以單值映射。

也就是說，不是直接從s域到z域，而是先將非帶限的H（s）映射為帶限的H（s’）,在模擬濾波器H（s）變換為數字濾波器H（z）的時候，再將s’映射到z域通過脈沖響應不變法的方法，即H（s）—H（s’）—H（z）。

（2）設計步驟

①將數字和模擬濾波器的頻率指標相互轉換，即{Ωk}轉為{ωk}。

②設計通帶截頻{ωp}通帶衰減Ap、阻帶衰減As、阻帶截頻{ωs}的模擬濾波器其中變換公式：

（3）

然而通帶最大衰減Ap和阻帶最小衰減As不變。

③利用雙線性不變法將模擬和數字濾波器系統函數進行轉換，即把H（s）通過相應的公式轉為H（z）。公式為：

（4）

轉換方法程序如下：

[N,Wc]=cheb1ord（Omgp,Omgs,Ap,As,'s'）%計算階數N和截止頻率ωc，[b,a]=cheby1（N,Ap,Wc,'s'）;%求模擬濾波器系數b和a，再由[bd,ad]=bilinear（b,a,0.5）;%用雙線性變換法求數字濾波器系數bd,ad例：用雙線性變換法設計IIR切比雪夫低通數字濾波器，wp=0.2πrad，ws=0.6πrad，Ap=3dB,As=40dB。

運行結果：

Omgp=0.3249

Omgs=1.3764

N=3

wo=0.3249

b=0 0 0 0.0086

a=1.0000 0.1941 0.0980 0.0086

圖7 雙線性變換法設計低通濾波器的分析

優點：由于z平面和s平面間的單一數值映射的關系，雙線性變換法不但能滿足對調換關系的兩點要求之外，而且還戰勝多值映射的關系，還可以消除頻率帶來的混疊，避免頻率失真現象[5]。

缺點:頻域的變換是非線性的，在高頻處有較大的失真，適用于低頻和中頻。

3 結論

采用脈沖響應不變法從時域上來設計IIR濾波器，美中不足的一定是由于頻率響應致使周期延拓而造成不可避免的混疊失真，因此這種方法只適合用來設計帶限的濾波器，如衰減特性很好和帶通濾波器的低通濾波器，高于頻率的折疊部分衰減越大，就會使失真情況減小，混疊效應越小。若從頻域出發的話，采取雙線性變換法。它是從s平面和z平面是一對一的關系，但克服了混疊的現象，成為這種方法的最大優點。

[1]王靖，李永全.數字橢圓濾波器的Matlab設計與實現[J].現代電子技術，2007（6）：4-5+8.

[2]袁煒.數字濾波器算法研究及其在DSP系統中的實現[D].武漢理工大學，2007

[3]辛修芳.數字濾波器在車輛電容稱重裝置中的應用[J].佳木斯大學學報：自然科學版，2012（3）：376-378.

[4]王淑嬌.數字濾波器的設計與參數選擇[J].科技資訊，2011（35）：32.

[5]封向東，李興源.遞歸數字濾波器的設計和實現[J].成都科技大學學報，1995（6）：84-92.

[責任編輯：袁向芬]

Design and implementation of IIR digital filter

ZHANG Hai-rong

（Guizhou Education University, Guiyang, Guizhou, 550018）

The digital filter are the important content digital signal processing,which may be divided into IIR and the FIR.To design the IIR digital filter , we need the help of analog prototype filter, and then transform analog filter into digital filter.In the paper we use the design of the pulse response invariable method and the bilinear method.Based on the IIR filter characteristic, we use the bilinearity method of transformation and respuesta invariable under the MATLAB bad boundary to design IIR,and compare the advantages and disadvantages of both.

Bigital filter IIR； MATLAB； Bigital filter； Bilinearity; Impulse response

2016-11-29

張海蓉（1983-），女，陜西涇陽人，碩士，貴州師范學院講師，研究方向：信號與信息處理、智能控制系統、嵌入式系統。

TN911.72

A

1674-7798（2016）12-0008-04