鄧海鵬IIR有源濾波器的設計實現(xiàn)分析

1、IIR 有源濾波器的設計與實現(xiàn)摘要：數(shù)字濾波器是數(shù)字信號處理中的重要組成局部 .它主要分為有限脈沖響應， (FIR) 和無限沖激響應(IIR)兩大類。本文主要介紹了對有限脈沖響應(FIR)和無限沖激響應(IIR)數(shù)字濾波器的原理與設計思路、設計方法、設計步驟進展了詳細的設計。在 MATLAB 環(huán)境下，通過對一個模擬原型法中的完全設計函數(shù)對低通 Butterworth IIR 數(shù)字濾波器的實現(xiàn)和用Ha ing窗函數(shù)法對低通FIR數(shù)字濾波器進展實現(xiàn)。實例發(fā)現(xiàn)， IIR 數(shù)字濾波器和 FIR 數(shù)字濾波器各有優(yōu)點，但是 FIR 數(shù)字濾波器比 IIR 數(shù)字濾波器 有穩(wěn)定、誤差小、運算效率高等優(yōu)勢。關鍵詞

2、. MATLAB 數(shù)字濾波器 IIRMatlab IIR based on the filter whose designAbstract： figures filter is the number of important part of the signalprocessing. it is limited mainly impulse response (FIR) and infiniteresponse at provoking (IIR) two categories. this paper mainly introducesimpulse response to a limited

3、 (FIR) and infi nite response at provoking(IIR) figures filter the principles and design methods,desi gn andimplementation procedures. detailed information on matlab environment, for example by the method of analysis to IIR figures filter and figuresare based on the design approach. fo：undIIR figure

4、s filter and FIRfigures filter have advantages, but the FIR figures filter than IIR figures filter a stable and error, the operation efficiency, etc.Keywords ： MATLAB IIR the function目錄引言 31 IIR數(shù)字濾波器 41.1數(shù)字濾波器的設計方法 41.2 IIR數(shù)字濾波器的設計思路 41.3利用模擬濾波器設計IIR數(shù)字濾波器的步驟 51.4 IIR濾波器的典型設計法 52 IIR濾波器的MATLAB的實現(xiàn)92.1

5、濾波器的MATLAB的實現(xiàn)9小結 10參考文獻 11引言隨著信息技術的迅猛開展，數(shù)字信號處理已成為一個極其重要的學科和技術 領域，在通信、語音、圖像、自動控制、雷達、軍事、航空航天、醫(yī)療和家用電 器等眾多領域得到了廣泛的應用。在數(shù)字信號處理的根本方法中。通常會涉及到 變換、濾波、頻譜分析、調制解調和編碼解碼等處理。其中，數(shù)字濾波是是數(shù)學 信號處理的重要環(huán)節(jié)，它在數(shù)字信號處理中占據(jù)著重要的地位，它具有可靠性好、 精度高、愛活性大、體積小、重量輕等優(yōu)點。隨著數(shù)字技術的開展，數(shù)字濾波器 越來越受到人們的重視，廣泛地應用于各個領域。MATLAB是于1984年由美國Mathworks公司推出，該軟件具有

6、使用簡單、方便， 易編程，語言簡練， 數(shù)據(jù)庫可任意擴大， 采用全新數(shù)據(jù)類型和面向對象編程技術等特 點，有強大的數(shù)值分析、矩陣運算、圖形繪制、數(shù)據(jù)處理等功能，因此已被廣泛應用 在數(shù)學、科研和工程設計的各個領域。如今 ATLAB 已經成為信號與信號處理，特 別是數(shù)字信號處理，Digital Sig nal Processi ngQSP應用中分析與仿真設計的主要工具。 本文將在ATLAB環(huán)境下，對有限脈沖響應(FIR)數(shù)字濾波器和無限脈沖響應(IIR)數(shù)字 濾波器的設計原理和設計方法及步驟進展分析，并結合實例根據(jù)對FIR數(shù)字濾波器和IIR數(shù)字濾波器的實現(xiàn)來對它們進展分析比擬。1.IIR數(shù)字濾波器II

7、R濾波器有以下幾個特點 單位沖激響應h(n)是無限長的， 系統(tǒng)函數(shù)Hz在有限z平面0<|z|v %上極點存在， 構造上存在著輸出到輸入的反應，也就是構造上是遞歸型的。1.1. IIR數(shù)字濾波器的設計方法1模擬原型法，采用經典低通濾波器作為連續(xù)域上的設計模型，使其逼近性能指 標，而后將模擬濾波器系統(tǒng)函數(shù)變換為所需的數(shù)字濾波器系統(tǒng)函數(shù)，最終得到IIR數(shù)字濾波器2直接設計法，基于給定的濾波器參數(shù)直接在離散域上尋找適宜的數(shù)字濾波器，不同于基于模擬濾波器的方法，該方法并不局限于常規(guī)的濾波器類型， 如低通、高通帶 通和帶阻等，而可以設計任意的頻率響應.3參數(shù)建模法，尋找接近于所需要設計的濾波器通用模

8、型，分為時域和頻域 建模兩類方法數(shù)字濾波器的設計思路IIR數(shù)字濾波器可用一個n階差分方程y(n) brx(n r) + v(h-k)或用它的Z域聚統(tǒng)函數(shù)對照槻擬濾波器的傳遞函數(shù)宀廠十乩5小十十花 十"-IS""十十心不難看出，數(shù)字濾波器與模擬濾波器的設計思路相仿，其設計實質也是尋找一組 系數(shù)b,a,去逼近所要求的頻率響應，使其在性能上滿足預定的技術要求，不同的是 模擬濾波器的設計是在S平面上用數(shù)學逼近發(fā)去尋找近似的所需特性H(S)而數(shù)字濾波器那么是在Z平面尋找適宜的H(Z),IIR數(shù)字濾波器的單位響應應是無限長的，而模擬濾波器一般都具有無限長的單位脈沖響應，因此與

9、模擬濾波器相匹配。由于模擬 濾波器的設計在理論上已經什么成熟，因此數(shù)字濾波器的設計的關鍵是將H(S)H(Z),即利用復值映射將模擬濾波器離散化。已經證明，沖擊響應不變法和雙線性變換法能 較好地擔當此任，那么在此根底上，數(shù)字濾波器的設計就可首先歸結為模擬濾波器的 設計了。數(shù)字濾波器設計步驟如下列圖所示，字潑技指模擬濾 注曙技 術抬標料應 模擬濾 數(shù)字 的濾 波耕離 濾波 迪曙f 熬化 f 器 沒計圖I數(shù)亍滝戲益詵計步腮1.3利用模擬濾波器設計IIR數(shù)字濾波器的步驟(1) 將給定的數(shù)字濾波器的性能指標，按某一變換(映射)規(guī)那么轉換成相應的模擬濾波器性能指標。(2) 如果要設計的不是數(shù)字濾波器，那么

10、還需要將步驟(1)中變換所得到相應的(高通，帶通，帶阻)模擬濾波器性能指標變換成模擬低通濾波器的性能指標這是因為模 擬低通濾波器才有圖形和表格可以利用.(3) 用所得到的模擬低通濾波器的性能指標，禾U用某種模擬濾波器逼近方法，設計 并查表求得此模擬低通濾波器的系統(tǒng)函數(shù)，以它作為設計數(shù)字濾波器的樣本。(4) 利用與步驟(1)，(2)中的同一變換規(guī)那么，將此作為“樣本的模擬原型低通濾波器的系統(tǒng)函數(shù)最終變換成所需的數(shù)字各型濾波器的系統(tǒng)函數(shù)H(z)步驟中的變換規(guī)那么就是從模擬濾波器數(shù)字化成為數(shù)字濾波器的方法，也就是要把s平面映射到z平面，使模擬系統(tǒng)函數(shù)Ha(s變換成所需的數(shù)字濾波器的系統(tǒng)函數(shù)H(z)+

11、這種由復變量s到復變量z之間的映射(變換)關系，必須滿足以下兩條根本要求，第一，H(z)的頻率 響應必須有可能模仿Ha(s)的頻率響應，即s平面的虛軸JQ必須映射到z平面的單位 圓ejw上，也就是頻率軸要對應。第二，果穩(wěn)定的Ha(s)應能映射成因果穩(wěn)定的H(z).也就是s平面的左半平面Res<o必須映射到z平面單位圓的內部z<11.4 IIR濾波器的典型設計法常見的模擬濾波器類型包括 Butterworth.巴特沃斯Chebyshev切比雪夫.I型、 ChebysheVI型和橢圓濾波器.基于這些模擬濾波器.ATLAB信號處理工具箱提供了 對應的IIR數(shù)字濾波器設計函數(shù)。(1)完全設

12、計函數(shù)MATLAB信號處理工具箱提供了幾個完全設計函數(shù).只要一條語句就可輕松地產生 任何階數(shù)、任何類型.低通、高通、帶通、帶阻.的IIR數(shù)字濾波器。butter函數(shù)。Butter函數(shù)用于設計Butterworth模擬和數(shù)字濾波器.該類型濾波器的特點是在通帶內具有最大平坦的幅度特性并且在通帶和阻帶內幅度都隨著頻率的增大而單調下降。函數(shù)的用法如下.b.a=butter(n,Rp,Wn).ba=butter(n,Rp,Wn,' ftype ')b.a=butter(n,Rp,Wn,' .a=bi)ttEr(n,Rp,Wn,' ftype ' , ' s

13、')z.p.k=butter(n,Rp,Wn,' ftype ' , ' s')A.B.C.D=butter(n,Rp,Wn,' ftype ' , ' s')函數(shù)輸出可以有濾波器的三種形式.傳遞函數(shù)形式.零極點增益形式和狀態(tài)方程形 式。a) 數(shù)字域b.a.=butter(n,Rp,Wn可設計出n階低通數(shù)字Butterworth濾波器.其截止頻率由 Wn 確定.通帶內的波紋由Rp確定.截比頻率是濾波器幅度下降到-Rp分貝處的頻率，且 是歸一化之后的X圍0, 1。通帶波紋Rp越小，可得到更寬的變換寬度。b,a，=butter

14、(n,Rp,Wn,' ftype')可設計高通、帶通和帶阻濾波器,Ftype=high時，可設計截止頻率為 Wn的高通濾波器，F(xiàn)type=band時，可設計帶通濾波器，Wn=W1 W2其帶通為 W1< 3 <W2。當 Wn=W1 W2時，ftype默認，butter函數(shù)可產生一個 2n階的數(shù)字帶通濾波器，其帶通為W1< 3 <W2。Ftype=stop時，可設計帶阻濾波器， 這時 Wn=W1 W2,且阻帶為 W1< 3 <W2。b) 模擬域b,a=butter(n,Rp,Wn, s')可設計出截止頻率為 Wn的n階低通模擬 Butte

15、rworth濾 波器。即H(G二型二盡I刃+執(zhí)2)卍】 + /+處f + 】) /f 何21 + £ + 鞏丹 +1)其截止頻率 Wn>0。注意這時Wn必須以弧度為單位。函數(shù)buttord可在給定濾波 器性能情況下，選擇模擬或數(shù)字 Butterworth濾波器最小的階數(shù)使用格式如下：n,Wn=buttord(Wp,Ws,Rp,Rs)數(shù)字域n,Wn=buttord(Wp,Ws,Rp,Rs； s')模擬域其中，Wp和Ws分別是通帶相阻帶的截止頻率，其值為0WWp(Ws)< 1 ,當其值為1時，表示0.5Fs Rp和Rs分別是通帶和阻帶內的波紋系數(shù)IS 2經奧設計出設計

16、就程圖Chebyl函數(shù)。Chebyl函數(shù)用于設計ChebyshevI型數(shù)字濾波器，該類型濾除波器幅頻特性在通帶內為等波紋，b,a= chebyl (n ,Rp,Wn)b,a= chebyl (n ,Rp,Wn,b,a= chebyl (n ,Rp,Wn,b,a= chebyl (n ,Rp,Wn,z,p,k= chebyl (n ,Rp,Wn,在阻帶單調下降。函數(shù)用法如下,'ftype ')'s')'ftype ' , ' s')'ftype ' , ' s')A, B, C, D= chebyl

17、(n,Rp,Wn, ' ftype ' , ' s')函數(shù)輸出可以有濾波器的三種形式，傳遞函數(shù)形式，零極點增益形式和狀態(tài)方程 形式。a)數(shù)字域b,a=cheby1(n,Rp,Wn可設計出n階低通數(shù)字ChebyshevI型濾波器，其截止頻率由Wn確定，通帶內的波紋由Rp確定，濾波器系統(tǒng)函數(shù)為,_侃】)十力(2)e ' +乙十何“十】)二"-1(-)截比頻率是濾波器幅度下降到-Rp分貝處的頻率，且是歸一化之后的 X圍0，10 通帶波紋Rp越小，可得到更寬的變換寬度。b,a= cheby1 (n,Rp,Wn： ftype')可設計 高通、帶通

18、和帶阻濾波器，F(xiàn)type=high時，可設計截止頻率為 Wn的高通濾波器， Ftype=band時，可設計帶通濾波器，Wn=W1 W2，其帶通為 W1< 3 <W2。當 Wn=W1 W2時，ftype缺省，cheby1函數(shù)可產生一個2n階的數(shù)字帶通濾波器，其帶通為 W1< 3 <W2 0Ftype=stop時，可設計帶阻濾波器，這時 Wn=W1 W2,且阻帶為 W1< 3 <W2。b)模擬域b,a= cheby1 (n,Rp,Wri, s')可設計出截止頻率為 Wn的n階低通模擬 Chebychev I型濾波器。即H(s)二二加1)弋十風2),八+上

19、十加“十1)其截止頻率 Wn>00注意這時 Wn必須以弧度為單位。cheby2函數(shù)。Cheby2函 數(shù)用于設計Chebyshev型數(shù)字濾波器，該類型濾波器的頻率響應特性在通帶內單調 下降，在阻帶為等波紋。Chebychev U型濾波器的下降斜度比I型小，但阻帶內波動 小。函數(shù)用法如下,b,a= cheby2 (n ,Rp,Wn) b,a= cheby2 (n ,Rp,Wn, b,a= cheby2 (n ,Rp,Wn, b,a= cheby2 (n ,Rp,Wn, z,p,k= cheby2 (n ,Rp,Wn,'ftype ')'s')'ftyp

20、e ' , ' s')'ftype ' , ' s')A，B，C， D= cheby2 (n,Rp,Wn,' ftype ' , ' s')函數(shù)輸出可以有濾波器的三種形式，傳遞函數(shù)形式，零極點增益形式和狀態(tài)方程 形式。ellip函數(shù)。Ellip函數(shù)用于設計橢圓數(shù)字濾波器，相對于 Butterworth濾波器和Chybeshev 濾波器而言，橢圓濾波器具有更陡峭的滾降特性，即在給定的階和波紋要 求下，橢圓濾波器能夠獲得最窄的過度寬帶，但在通帶和阻帶內都為等波紋形式。函數(shù)用法如下，b,a= ellip(n,Rp

21、,Rs,Wn)b,a= ellip (n,Rp,Rs,Wn,' ftype ')b,a= ellip(n,Rp,Rs,Wn, ' s')b,a= ellip(n,Rp,Rs,Wn, z,p,k= ellip (n,Rp,Rs,Wn,ftype ', 's')'ftype ', 's')ftype ', 's')A， B， C， D= ellip (n,Rp,Rs,Wn,函數(shù)輸出可以有濾波器的三種形式，傳遞函數(shù)形式，零極點增益形式和狀態(tài)方程 形式。其中， Wn 為截止頻率，歸一化之后的

22、值， X 圍0， 1。 Rpz 指定通帶波紋， Rs指定阻帶波紋。函數(shù)b,a= ellip (n,Rp,Rs,Wn' ftype'河設計高通、帶通和帶阻濾 波器， Ftype=high 時，可設計截止頻率為 Wn 的高通濾波器， Ftype=band 時，可設 計帶通濾波器，Wn=W1 W2的通帶為 W1< 3 <W2。當 Wn=W1 W2時，ftype缺省， ellip函數(shù)可產生一個-2n階的數(shù)字帶通濾波器，其帶通為W1< 3 <W2。 Ftype=stop時 ， 可 設 計 帶 阻 濾 波 器 ， 這 時 Wn=W1 W2, 且 阻 帶 為 W1&l

23、t; 3 <W2 。 b,a=ellip(n,Rp,Rs,Wn： s')和b,a=ellip(n,Rp,Rs,Wn： ftype' ,' s')用來設計低通、 高通、帶通和帶阻橢圓濾波器。這時 Wn 必須以弧度為單位。2.1IIR 濾波器的 MATLAB 實現(xiàn)利用ATLAB進展IIR數(shù)字低通濾波器設計，首先是在給定濾波器性能情況下，確 定通帶和阻帶的截止頻率，Butterworth濾波器buttord函數(shù)計算出滿足指定性能的模擬Butterworth濾波器的階和截止頻率。有了上 述的前提和根底，我們可利用ATLAB語言中Butter函數(shù)設計出模擬域Butt

24、erworth濾波 器。然后再利用i piny盯函數(shù)，用脈沖響應不變法實現(xiàn)模擬到數(shù)字濾波器的轉換。設 計一個8階低通Butterworth濾波器，要求截止頻率300Hz給定采樣頻1000Hz.fs = 1000;n = 10;Wn = 300*2/fs;b,a， = butter(n,Wn,' low' ) = 256;freqz(b,a, ,fs)程序運行結果如下列圖所示，品一三一亍p匸200 =00Fruncy (FL)-10uDo20O50i(5気S旦爲荒200300Ffwqu&ncy (Hz)小結、數(shù)字濾波器以他優(yōu)越的性能，在數(shù)字信號處理領域中占有很重要的地位，是現(xiàn)代電子技術中必須掌握的設計技術。本文對IIR數(shù)字濾波器和FIR數(shù)字濾波器的的概 念、特點、設計思路、設計方法以及設計步驟進展了詳細的設計。并在MATLAB的環(huán)境下，通過一個模擬原型法中的完全設計函數(shù)對低通ButterworthllR數(shù)字濾波器的實現(xiàn)和用Ha ing窗函數(shù)法對低通FIR數(shù)字濾波器進展實現(xiàn)。可以發(fā)現(xiàn)，設計結果的 各項性能指標均到達指定要求，驗