IIR數(shù)字濾波器的設計流程圖DOC

1、第 頁共17頁目錄TOC o 1-5 h z HYPERLINK l bookmark0 o Current Document 目錄0 HYPERLINK l bookmark2 o Current Document 前言1 HYPERLINK l bookmark4 o Current Document 數(shù)字濾波器簡介1 HYPERLINK l bookmark6 o Current Document 使用數(shù)字濾波器的原因1 HYPERLINK l bookmark8 o Current Document 設計的原理和內容1工程概況2正文2 HYPERLINK l bookmark10 o

2、Current Document 3.1設計的目的和意義2 HYPERLINK l bookmark12 o Current Document 3.2目標和總體方案2 HYPERLINK l bookmark14 o Current Document 設計方法和內容3 HYPERLINK l bookmark16 o Current Document 硬件環(huán)境3 HYPERLINK l bookmark18 o Current Document 軟件環(huán)境3 HYPERLINK l bookmark20 o Current Document 3.6IIR數(shù)字濾波器設計思路3 HYPERLINK

3、l bookmark22 o Current Document IIR數(shù)字濾波器的設計流程圖3 HYPERLINK l bookmark26 o Current Document IIR數(shù)字濾波器設計思路4 HYPERLINK l bookmark40 o Current Document 3.9設計IIR數(shù)字濾波器的兩種方法4 HYPERLINK l bookmark42 o Current Document 3.10雙線性變換法的基本原理5 HYPERLINK l bookmark60 o Current Document 3.11用雙線性變換法設計IIR數(shù)字濾波器的步驟6 HYPERLI

4、NK l bookmark62 o Current Document 3.12程序源代碼和運行結果6低通濾波器63.12.3帶通濾波器103.12.4帶阻濾波器13 HYPERLINK l bookmark64 o Current Document 3.13結論15存在的問題15解決方案16 HYPERLINK l bookmark66 o Current Document 致謝1616參考文獻前言1.1數(shù)字濾波器簡介數(shù)字濾波器是一種用來過濾時間離散信號的數(shù)字系統(tǒng)，通過對抽樣數(shù)據進行數(shù)學處理來達到頻域濾波的目的。可以設計系統(tǒng)的頻率響應，讓它滿足一定的要求，從而對通過該系統(tǒng)的信號的某些特定的頻率

5、成分進行過濾，這就是濾波器的基本原理。如果系統(tǒng)是一個連續(xù)系統(tǒng)，則濾波器稱為模擬濾波器。如果系統(tǒng)是一個離散系統(tǒng)，則濾波器稱為數(shù)字濾波器。信號通過線性系統(tǒng)后，其輸出信號就是輸入信號和系統(tǒng)沖激響應的卷積。從頻域分析來看，信號通過線性系統(tǒng)后，輸出信號的頻譜將是輸入信號的頻譜與系統(tǒng)傳遞函數(shù)的乘積。除非為常數(shù)，否則輸出信號的頻譜將不同于輸入信號的頻譜，某些頻率成分較大的模，因此，中這些頻率成分將得到加強，而另外一些頻率成分的模很小甚至為零，中這部分頻率分量將被削弱或消失。因此，系統(tǒng)的作用相當于對輸入信號的頻譜進行加權。使用數(shù)字濾波器的原因數(shù)字濾波器具有比模擬濾波器更高的精度，甚至能夠實現(xiàn)后者在理論上也無法

6、達到的性能。數(shù)字濾波器相比模擬濾波器有更高的信噪比。數(shù)字濾波器還具有模擬濾波器不能比擬的可靠性。根據其沖擊響應函數(shù)的時域特性可將數(shù)字濾波器分為iir(有限長沖擊響應)和fir(無限長沖擊響應)。設計的原理和內容在windows環(huán)境下進行語言信號采集，通過iir數(shù)字濾潑器的設計，數(shù)字帶濾波器就是用軟件來實現(xiàn)上面的濾波過程，可以很好的克服模擬濾波器的缺點，數(shù)字帶濾波器的參數(shù)一旦確定，就不會發(fā)生變化。IIR型有較好的通帶與阻帶特性，所以，在一般的設計中選用IIR型。IIR型又可以分成Butterworth型濾波器，ChebyshevII型濾波器和橢圓型濾波器等。IIR數(shù)字濾波器的設計一般是利用目前已

7、經很成熟的模擬濾波器的設計方法來進行設計，通常采用模擬濾波器原型有butterworth函數(shù)、chebyshev函數(shù)、bessel函數(shù)、橢圓濾波器函數(shù)等。IIR數(shù)字濾波器的設計步驟：按照一定規(guī)則把給定的濾波器技術指標轉換為模擬低通濾波器的技術指標；根據模擬濾波器技術指標設計為響應的模擬低通濾波器；很據脈沖響應不變法和雙線性不變法把模擬濾波器轉換為數(shù)字濾波器；*學院（4）如果要設計的濾波器是高通、帶通或帶阻濾波器，則首先把它們的技術指標轉化為模擬低通濾波器的技術指標，設計為數(shù)字低通濾波器，最后通過頻率轉換的方法來得到所要的濾波器。我這次做的信號與系統(tǒng)課程設計的任務是在MATLAB軟件學過的有關知

8、識設計簡單的IIR數(shù)字濾波器續(xù)系統(tǒng)，本次采用雙線性變換法設計IIR濾波器。在做課程設計中首先是對MATLAB軟件的了解和認識，掌握一些MATLAB軟件的基本常用函數(shù)的用法，對MATLAB軟件進行程序操作。同時利用MATLAB軟件也能對書本上的知識進行驗證，在MATLAB軟件下編寫函數(shù)程序，然后運行程序，與書本上的信號的求解進行對照分析和比較。對MATLAB軟件進行一定的了解和運用之后，開始做此次課正程設文計。數(shù)字濾波器是一個離散時間系統(tǒng)（按預定的算法，將輸入離散時間信號轉換為所要求的輸出離散時間信號的特定功能裝置）。應用數(shù)字濾波器處理模擬信號時，首先須對輸入模擬信號進行限帶、抽樣和模數(shù)轉換。數(shù)

9、字濾波器輸入信號的抽樣率應大于被處理信號帶寬的兩倍，其頻率響應具有以抽樣頻率為間隔的周期重復特性，且以折疊頻率即12抽樣頻率點呈鏡像對稱。為得到模擬信號，數(shù)字濾波器處理的輸出數(shù)字信號須經數(shù)模轉換、平滑。數(shù)字濾波器具有高精度、高可靠性、可程控改變特性或復用、便于集成等優(yōu)點。數(shù)字濾波器在語言信號處理、圖像信號處理、醫(yī)學生物信號處理以及其他應用領域都得到了廣泛應用。數(shù)字濾波器有低通、高通、帶通、帶阻和全通等類型。它可以是時不變的或時變的、因果的或非因果的、線性的或非線性的。應用最廣的是線性、時不變數(shù)字濾波器。3.1設計的目的和意義我們是通信工程的本科生，信號與系統(tǒng)是我們重要的必修課程。當代社會學要大

10、學培養(yǎng)出理論扎實，動手實踐能力強的大學生。所以，本次課程設計的目的就在于通過一次實踐性的活動加深對這門課程的理解，使我們在感性的認識上進一步升華為理性的認識。為后繼課程的學習打下堅實的基礎。增進對MATLAB的認識，加深對數(shù)字信號處理理論方面的理解，掌握數(shù)字信號處理中IIR濾波器的設計。馬克思主義唯物辯證法認為，實踐是連接客觀實在和人主觀意識的通道和橋梁。物質對意識的作用以及意識對物質的反作用都蘊含在實踐活動當中。也就是，實踐是檢驗真理的唯一標準。對這門課的學習狀況的好壞，用一次課程設計便可以檢驗出來。而這，就是本次我們進行設計的意義之所在。3.2目標和總體方案IIR被稱為無限長脈沖響應濾波器

11、，結構中有反饋，所以理論上脈沖響應永遠不為零。實現(xiàn)相同的指標比IIR濾波器需要的階數(shù)較少。非線性相位，用于對相位不太敏感的場合，如語音。設計時先設計模擬濾波器，再由一定的轉換原則轉換為數(shù)字濾波器。本次設計的目標在于將IIR數(shù)字濾波器利用MATLAB的知識編寫相應程序實現(xiàn)高通、帶通或帶阻濾波的功能。于是特制訂了一個總體的方案。由于時間只有十天，故做了如下的計劃安排，將這項工程分為兩大部分：程序的設計和程序的調試。首先在程序的設計部分由分為幾個步驟：第一步：查閱有關IIR數(shù)字濾波器設計的資料，用三天的時間。第二步：對設計高通、帶通或帶阻濾波器時所用到的參數(shù)的作用進行更深層的了解，用第2頁共17頁一

12、到四天的時間。第三步：設計高通、帶通或帶阻濾波器的代碼，兩天的時間。其次，進行程序的調試。用一天。3.3設計方法和內容“工欲善其事，必先利其器”。有了總體方案后必須用一個事半功倍的設計方法來提高程序設計的效率。IIR數(shù)字濾波器具有無限寬的沖激響應，與模擬濾波器相匹配。所以IIR濾波器的設計可以采取在模擬濾波器設計的基礎上進一步變換的方法。IIR通常是非線性的，但是目前也有準線性相位設計方法得到IIR數(shù)字濾波器的系數(shù)，其結果是使得通帶內的相位波動維持在一個工程可接受的范圍內oIIR比FIR最大的優(yōu)點是達到同樣的矩形系數(shù)所需的階數(shù)少，往往5階的IIR濾波器就可以比擬數(shù)十上百階的FIR濾波器。設計方

13、面，我們將考慮映射為數(shù)字濾波器的方法。首先，重點研究由低通濾波器設計數(shù)字低統(tǒng)濾波器的方法，然后將這些設計轉換成更通用的頻率選擇濾波器的方法。硬件環(huán)境微型計算機：聯(lián)想臺式品牌機中央處理器：Pentuim4主頻：3.0GHz主存容量：512M硬盤容量：80G軟件環(huán)境WindowsXP操作系統(tǒng)MATBLE軟件程序3.6IIR數(shù)字濾波器設計思路IIR數(shù)字濾波器可用一個n階差分方程表示y(n)二工bx(n-r)+Eay(n-k)rk不難看出,數(shù)字濾波器與模擬濾波器的設計思路相仿,其設計實質也是尋找一組系數(shù)b,a，去逼近所要求的頻率響應，使其在性能上滿足預定的技術要求;不同的是模擬濾波器的設計是在S平面上

14、用數(shù)學逼近法去尋找近似的所需特性H(S),而數(shù)字濾波器則是在Z平面尋找合適的H(z)oIIR數(shù)字濾波器的單位響應是無限長的，而模擬濾波器一般都具有無限長的單位脈沖響應，因此與模擬濾波器相匹配。由于模擬濾波器的設計在理論上已十分成熟，因此數(shù)字濾波器設計的關鍵是將H(S)H(Z),即,利用復值映射將模擬濾波器離散化。已經證明，沖擊響應不變法和雙線性變換法能較好地擔當此任，則在此基礎上，數(shù)字濾波器的設計就可首先歸結為模擬濾波器的設計了。3.7IIR數(shù)字濾波器的設計流程圖*學院第 頁共17頁圖3-1數(shù)字濾波器的設計流程圖3.8IIR數(shù)字濾波器設計思路IIR數(shù)字濾波器可用一個n階差分方程表示y(n)=工

15、bx(n一r)+rk或用它的Z域系統(tǒng)函數(shù)：工bZ-rH(z)=于1一厶az-kk對照模擬濾波器的傳遞函數(shù)：bSm+bSM-i+bH(S)=mm-10-aSm+aSm-1+amm-10不難看出,數(shù)字濾波器與模擬濾波器的設計思路相仿,其設計實質也是尋找一組系數(shù)b,a,去逼近所要求的頻率響應，使其在性能上滿足預定的技術要求;不同的是模擬濾波器的設計是在S平面上用數(shù)學逼近法去尋找近似的所需特性H(S),而數(shù)字濾波器則是在Z平面尋找合適的H(z)。IIR數(shù)字濾波器的單位響應是無限長的，而模擬濾波器一般都具有無限長的單位脈沖響應,因此與模擬濾波器相匹配。由于模擬濾波器的設計在理論上已十分成熟,因此數(shù)字濾波

16、器設計的關鍵是將H(S)H(Z),即,利用復值映射將模擬濾波器離散化。已經證明,沖擊響應不變法和雙線性變換法能較好地擔當此任,則在此基礎上,數(shù)字濾波器的設計就可首先歸結為模擬濾波器的設計了。3.9設計IIR數(shù)字濾波器的兩種方法IIR數(shù)字濾波器的的設計就是在給定了濾波器的技術指標后，確定濾波器的階數(shù)n和系數(shù)ai,bi。在滿足技術指標的條件下，濾波器的結束應盡可能低，因為濾波器的階數(shù)越低，實現(xiàn)濾波器的成本就越低。在設計IIR濾波器時，常用的方法是利用模擬濾波器來設計數(shù)字濾波器。廣泛采取這種方法的因素有：1，模擬濾波器設計技術已非常成熟；2，可得閉合形式的解；3，關于模擬濾波器設計有完整的設計公式和

17、圖表可以利用很查閱。為實現(xiàn)從模擬濾波器到數(shù)字濾波器的轉換，需要從系統(tǒng)的描述方法來考慮轉單位采樣響應h(n)換問題，無論是模擬濾波器還是數(shù)字濾波器，描述系統(tǒng)的基本方法都有四種。濾波器描述系統(tǒng)的方法因此，IIR濾波器的設計方法是首先將數(shù)字濾波器的技術指標轉化為對應模擬濾波器的技術指標，然后設計滿足技術指標的模擬濾波器Ha(s),然后將設計出的模擬濾波器Ha(s)轉換為滿足技術指標的數(shù)字濾波器H(z)。將Ha(s)轉換成H(z)的最終目的，是希望數(shù)字濾波器的頻率響應H(ej3)盡量接近模擬濾波器Ha(jQ)。將系統(tǒng)函數(shù)H(z)從s平面轉換到z平面的方法有很多種，但工程上常用的有兩種：一種是使數(shù)字濾波

18、器的h(n)近似于模擬濾波器的ha(t),可導出脈沖響應不變法；另一種使數(shù)字濾波器的差分方程近似于模擬濾波器的微分方程，由此可導出雙線性變換法。3.10雙線性變換法的基本原理脈沖響應不變法使得數(shù)字濾波器在時域上能夠較好的模仿模擬濾波器，但是由于從s平面到z平面的映射z二esT具有多值性，使得設計出來的數(shù)字濾波器不可避免的出現(xiàn)頻譜混迭現(xiàn)象。為了克服脈沖響應不變法可能產生的頻譜混跌效應的缺點，我們使用一種新的變換雙線性變換。雙線性變換法可認為是基于對微分方程的積分，利用對積分的數(shù)值逼近的道德。仿真濾波器的傳遞函數(shù)H(s)為dsk將展開為部份分式的形式，并假設無重復幾點，則H(s)仝-Assk=1p

19、k那么，對于上述函數(shù)所表達的數(shù)字信號處理系統(tǒng)來講，其仿真輸入x(t)和模擬輸出y(t)有如下關系心)-spy(t)=Ax(t)利用差分方程來代替導數(shù)，即同時令y(t)x(t)=n這樣，便可將上面的微分方程寫為對應的差分方程形式y(tǒng)(n)-y(n-1)專y(n)+y(n-1)=A為+x(n-1)兩邊分別取z變換，可得H(z)=竺=_AX2x-sT1+z-ip這樣，通過上述過程，就可得到雙線性變換中的基本關系，如下所示21-z-1s=xT1+z-i所謂的雙線性變換，僅是指變換公式中s與z的關系無論是分子部份還是分母部份都是線性的。3.11用雙線性變換法設計IIR數(shù)字濾波器的步驟MATLAB中設計II

20、R數(shù)字濾波器的具體步驟如下：把給出的數(shù)字濾波器的性能指標轉換為模擬低通濾波器的性能指標；根據轉換后的性能指標，通過濾波器結束選擇函數(shù)，來確定濾波器的最小階數(shù)n和固有頻率wn；由最小階數(shù)n得到低通濾波器原型；由固有頻率wn把模擬低通濾波器轉換為模擬低通、高通、帶通或帶阻濾波器；運用雙線性變換法把模擬濾波器轉換成數(shù)字濾波器。3.12程序源代碼和運行結果3.12.1低通濾波器clearwp=100*2*pi;%通帶截止頻率ws=150*2*pi;%阻帶截止頻率rp=0.5;%通帶衰減rs=30;%阻帶衰減fs=2000;%采樣頻率n,wc=cheb2ord(wp,ws,rp,rs,s)%計算階數(shù)，與

21、截止頻率z,p,k=cheb2ap(n,rs);%建立切比雪夫2型數(shù)字濾波器%零極點轉換到空間狀態(tài)表達式a,b,c,d=zp2ss(z,p,k);%零極點轉換到空間狀態(tài)表達式at1,bt1,ct1,dt1=lp2lp(a,b,c,d,wc);%低通轉換到高通at2,bt2,ct2,dt2=bilinear(at1,bt1,ct1,dt1,fs);%雙線性變換num,den=ss2tf(at2,bt2,ct2,dt2)%空間狀態(tài)表達式轉換到傳遞函數(shù)%繪制幅頻、相頻圖(頻率響應特性圖)*學院-4000第 頁共17頁figure(1);freqz(num,den,128,fs);gridon;tit

22、le(幅頻、相頻圖)；%繪制脈沖響應特性圖figure(2);impz(num,den,128,fs);grid;title(脈沖響應特性圖)%濾波檢驗figure(3);t=0:0.0005:0.2;x=sin(2*pi*50*t)+sin(2*pi*200*t);y=filter(num,den,x);plot(t,x,:,t,y,-);grid;legend(XSignal,YSignal);title(濾波檢驗)運行結果：n=6500-50-100幅頻、相頻圖1002003004005006007008009001000Frequency(Hz)Bd(edutingawc=879.25

23、59num=0.0287-0.10850.2038-0.24470.2038-0.10850.0287den=1.0000-4.44998.4145-8.61765.0302-1.58410.210301002003004005006007008009001000Frequency(Hz)1000-100-200-300圖3-2*學院第 頁共17頁0.120.10.080.060.040.020-0.02-0.04脈沖響應特性圖0.010.050.060.020.030.04nT(seconds)圖3-3濾波檢驗21.510.50-0.5-1-1.5_XSignalYSignal-IInli1

24、-20I0.020.040.060.080.10.120.140.160.180.2圖3-43.12.2高通濾波器wp=100*2*pi;%通帶截止頻率ws=150*2*pi;%阻帶截止頻率rp=0.5;%通帶衰減rs=30;%阻帶衰減fs=2000;%采樣頻率n,wc=cheb2ord(wp,ws,rp,rs,s);%計算階數(shù)，與截止頻率z,p,k=cheb2ap(n,rs);%建立切比雪夫2型數(shù)字濾波器a,b,c,d=zp2ss(z,p,k);%零極點轉換到空間狀態(tài)表達式at1,bt1,ct1,dt1=lp2hp(a,b,c,d,wc)%低通轉換到高通at2,bt2,ct2,dt2=bil

25、inear(at1,bt1,ct1,dt1,fs)%雙線性變換num,den=ss2tf(at2,bt2,ct2,dt2);%空間狀態(tài)表達式轉換到傳遞函數(shù)%繪制幅頻、相頻圖(頻率響應特性圖)figure(1);freqz(num,den,128,fs);gridon;title(幅頻、相頻圖)%繪制脈沖響應特性圖figure(2);impz(num,den,128,fs);grid;title(脈沖響應特性圖)%濾波檢驗figure(3);t=0:0.0005:0.1;x=sin(2*pi*50*t)+sin(2*pi*200*t);y=filter(num,den,x);plot(t,x,:

26、,t,y,-);legend(XSignal,YSignal);title（濾波檢驗）運行結果：n=6wc=879.2559num=0.5365-3.06887.4574-9.8501den=1.0000-4.52158.8272-9.47317.4574-3.06880.53655.8745-1.99140.2878幅頻、相頻圖圖3-5脈沖響應特性圖圖3-6圖3-73.12.3帶通濾波器clear;%下通帶截止頻率wpl=80HZ,上通帶截止頻率wp2=170HZ，阻帶下限頻率wsl=100HZ，阻帶上限頻率ws2=150HZ；解得wp=85HZ,ws=139Hzwp=85*2*pi;%通帶

27、截止頻率wpws=139*2*pi;%阻帶截止頻率wsrp=0.5;%通帶衰減rprs=30;%阻帶衰減rsfs=2000;%采樣頻率n,wc二cheb2ord(wp,ws,rp,rs,s)%計算階數(shù)n,與截止頻率wnz,p,k=cheb2ap(n,rs);%建立切比雪夫2型數(shù)字濾波器a,b,c,d=zp2ss(z,p,k);%零極點模型轉換到空間狀態(tài)模型at1,bt1,ct1,dt1=lp2bp(a,b,c,d,wc,50*2*pi);%低通轉換到高通at2,bt2,ct2,dt2=bilinear(at1,bt1,ct1,dt1,fs);%用雙線性變換法實現(xiàn)模擬濾波器到數(shù)字濾波器的轉換nu

28、m,den=ss2tf(at2,bt2,ct2,dt2);%空間狀態(tài)模型轉換到傳遞函數(shù)模型%繪制幅頻、相頻圖(頻率響應特性圖)figure(1);freqz(num,den,128,fs);gridon;title(幅頻、相頻圖)%繪制脈沖響應特性圖figure(2);impz(num,den,128,fs);grid;title(脈沖響應特性圖)%濾波檢驗figure(3);t=0:0.0005:0.3;x=sin(2*pi*50*t)+sin(2*pi*150*t);y=filter(num,den,x);plot(t,x,:,t,y,-);legend(XSignal,YSignal);

29、title(濾波檢驗)運行結果：n=5wc=849.6593num=0.0100-0.07200.2250-0.37870.32070.0000-0.32070.3787-0.22500.0720-0.0100den=1.0000-8.773735.4059-86.4271141.2110-161.2926130.4129-73.714527.8888-6.38240.6718Frequency(Hz)Frequency(Hz)圖3-8脈沖響應特性圖圖3-9濾波檢驗21.510.50-0.5-1-1.5LIXSignalYSignal-200.050.10.150.20.250.30.35圖3

30、-103.12.4帶阻濾波器clear;%下通帶截止頻率wpl=80HZ,上通帶截止頻率wp2=170HZ，阻帶下限頻率wsl=100HZ，阻帶上限頻率ws2=150HZ；解得wp=89HZ,ws=225Hzwp=89*2*pi;%通帶截止頻率wpws=225*2*pi;%阻帶截止頻率wsrp=0.5;%通帶衰減rprs=30;%阻帶衰減rsfs=2000;%采樣頻率n,wc二cheb2ord(wp,ws,rp,rs,s);%計算階數(shù)n,與截止頻率wn%建立切比雪夫2型數(shù)字濾波器%計算階數(shù)n,與截止頻率wnz,p,k=cheb2ap(n,rs);%建立切比雪夫2型數(shù)字濾波器%計算階數(shù)n,與截止頻率wna,b,c,d=zp2ss(z,p,k);%零極點模型轉換到空間狀態(tài)模型at1,bt1,ct1,dt1=lp2bs(a,b,c,d,wc,100*2*pi);%低通轉換到高通at2,bt2,ct2,dt2=bilinear(at1,bt1,ct1,dt1,fs);%用雙線性變換法實現(xiàn)模擬濾波器到數(shù)字濾波器的轉