基于MATLAB的數(shù)字濾波器設(shè)計畢業(yè)論文設(shè)計(原創(chuàng)版)[1]

1、編號 淮安信息職業(yè)技術(shù)學院畢業(yè)論文題 目基于MATLAB的數(shù)字濾波器設(shè)計學生姓名。學 號。系 部。專 業(yè)。班 級。指導教師。 顧問教師。 。摘 要本論文介紹了FIR數(shù)字濾波器的設(shè)計方法，即窗函數(shù)法。在此基礎(chǔ)上，用MATLAB實現(xiàn)IIR數(shù)字濾波器。介紹了IIR數(shù)字濾波器的傳統(tǒng)設(shè)計思想與步驟，及其計算機輔助設(shè)計方法。以一數(shù)字帶通濾波器為例,著重說明了基于MATLAB的三種濾波器的實現(xiàn)手段：模擬低通原型、合適模擬帶通及直接原型,為數(shù)字濾波器設(shè)計帶來全新的實現(xiàn)手段。關(guān)鍵詞：濾波 IIR濾波器 FIR濾波器 MATLAB目 錄第一章  前言 1.1 MATLAB 軟件簡介 (4)1.2數(shù)字濾波

2、器技術(shù)的發(fā)展狀況 (5)第二章  數(shù)字濾波器的基本概念2.1數(shù)字濾波器的概況 (6)2.2 FIR 數(shù)字濾波器的基本概念 (6)2.2.1 FIR 數(shù)字濾波器的窗函數(shù)設(shè)計法 (7)2.2.2 窗函數(shù)設(shè)計法的步驟 (7)2.3 MATLAB環(huán)境下的實例 (9)2.3.1高通濾波器的設(shè)計 (9)2.3.2低通濾波器的設(shè)計 (10)第三章  IIR數(shù)字濾波器的設(shè)計過程及方法3.1 IIR濾波器的基本特點 (13)3.2 IIR濾波器的設(shè)計思路與步驟 (14)3.3 IIR 濾波器的設(shè)計 (14)3.4 IIR濾波器設(shè)計方法MATLAB的實現(xiàn) (15)3.4.1 基于模擬低通原型的

3、MATLAB實現(xiàn) (15)3.4.2基于合適類型模擬濾波器的MATLAB實現(xiàn) (16)3.4.3 基于直接原型變換法的MATLAB實現(xiàn) (18)總  結(jié)(19)參考文獻(20)第一章  前言1.1. MATLAB簡介MATLAB (Matrix Laboratory)為美國Mathworks公司1983年首次推出的一套高性能的數(shù)值分析和計算軟件，其功能不斷擴充，版本不斷升級，1992年推出劃時代的4.0版，1993年推出了可以配合Microsoft Windous使用的微機版，95年4.2版，97年5.0版，99年5.3版，5.X版無論是界面還是內(nèi)容都有長足的進展，其幫助信

4、息采用超文本格式和PDF格式，可以方便的瀏覽。至2001年6月推出6.1版，2002年6月推出6.5版，繼而推出6.5.1版, 2004年7月MATLAB7和Simulink6.0被推出，目前的最新版本為7.1版。     MATLAB將矩陣運算、數(shù)值分析、圖形處理、編程技術(shù)結(jié)合在一起，為用戶提供了一個強有力的科學及工程問題的分析計算和程序設(shè)計工具，它還提供了專業(yè)水平的符號計算、文字處理、可視化建模仿真和實時控制等功能，是具有全部語言功能和特征的新一代軟件開發(fā)平臺。     MATLAB 已發(fā)展成為適合眾多學科，多種工作平臺、功能強大

5、的大型軟件。在歐美等國家的高校，MATLAB已成為線性代數(shù)、自動控制理論、數(shù)理統(tǒng)計、數(shù)字信號處理、時間序列分析、動態(tài)系統(tǒng)仿真等高級課程的基本教學工具。成為攻讀學位的本科、碩士、博士生必須掌握的基本技能。在設(shè)計研究單位和工業(yè)開發(fā)部門，MATLAB被廣泛的應用于研究和解決各種具體問題。在中國，MATLAB也已日益受到重視，短時間內(nèi)就將盛行起來，因為無論哪個學科或工程領(lǐng)域都可以從MATLAB中找到合適的功能。 MATLAB6.1 的新功能 （1） 全新的開發(fā)環(huán)境GUIDE(GUI Development Environment)：由早期的單一命令窗口發(fā)展為將一些常用的交互式工作界面高度的集中于操作桌

6、面； （2） 在命令窗口增加了錯誤跟蹤功能(error display message and abort function)； （3） 提供了許多新的數(shù)值處理功能(numerical Treatment Techniques)，更新了部分函數(shù)的功能和算法，增加了許多新函數(shù)； （4） 在圖形窗口增加了曲線擬合、數(shù)據(jù)統(tǒng)計等交互工具；(Curve Fitting, Data Analyzing ) （5） 引入了類與對象及函數(shù)句柄等概念;(Object，Handle) 用戶可以創(chuàng)建自己定義的類函數(shù)和函數(shù)句柄； （6） 改進了編輯/調(diào)試器的界面及功能; (Editor/Debugger)增加了行號和

7、書簽等功能； （7） 屬性編輯器功能更強(Property Debugger)，使用更方便； （8） 建立了一個與以前完全不同的圖形用戶界面（GUI）圖形窗口，(Graphical Window)使用更加方便靈活； （9） 增加了虛擬現(xiàn)實工具箱，(Virtual Reality Toolbox)使用標準的虛擬現(xiàn)實建模語言（VRML）技術(shù)，可以創(chuàng)建由MATLAB和simulink環(huán)境驅(qū)動的三維動畫場景； （10） 在應用程序接口方面增加了與Java的接口(Interface for Java)，并為二者的數(shù)據(jù)交換提供了相應的程序庫。1.2數(shù)字濾波器技術(shù)的發(fā)展狀況及簡介  &#

8、160; 數(shù)字濾波器是數(shù)字信號處理理論的一部分。數(shù)字信號處理主要是研究用數(shù)字或符號的序列來表示信號波形，并用數(shù)字的方式去處理這些序列，把它們改變成在某種意義上更為有希望的形式，以便估計信號的特征參量，或削弱信號中的多余分量和增強信號中的有用分量。具體來說，凡是用數(shù)字方式對信號進行濾波、變換、調(diào)制、解調(diào)、均衡、增強、壓縮、固定、識別、產(chǎn)生等加工處理，都可納入數(shù)字信號處理領(lǐng)域。 濾波器的種類很多，從功能上可以分為低通、高通、帶通和帶阻濾波器，上述每種濾波器又可以分為模擬濾波器和數(shù)字濾波器。如果濾波器的輸入輸出都是數(shù)字信號，則這樣的濾波器稱之為數(shù)字濾波器，它通常通過一定的運算關(guān)系改變輸入信號所含頻率

9、成分的相對比例或者濾除某些頻率成分來實現(xiàn)濾波。根據(jù)數(shù)字濾波器沖激響應的時域特性，可將數(shù)字濾波器分為兩種，即無限長沖激響應（IIR）濾波器和有限長沖激響應（FIR）濾波器。有數(shù)字信號處理的一般理論可知，IIR濾波器的特征是具有無限持續(xù)時間的沖激響應，而FIR濾波器使沖激響應只能持續(xù)一定的時間。第二章 數(shù)字濾波器的基本概念2.1數(shù)字濾波器的概況數(shù)字濾波技術(shù)是數(shù)字信號處理的一個重要組成部分，濾波器的設(shè)計是信號處理的核心問題之一。根據(jù)FIR濾波器的原理，提出了FIR濾波器的窗函數(shù)設(shè)計法，并對常用的幾種窗函數(shù)進行了比較。給出了在MATLAB環(huán)境下，用窗函數(shù)法設(shè)計FIR濾波器的過程和設(shè)計實例。仿真結(jié)果表明

10、，設(shè)計的FIR濾波器的各項性能指標均達到了指定要求，設(shè)計過程簡便易行。該方法為快速、高效地設(shè)計FIR濾波器提供了一個可靠而有效的途徑。 隨著信息時代的到來，數(shù)字信號處理已經(jīng)成為當今一門極其重要的學科和技術(shù)，并且在通信、語音、圖像、自動控制等眾多領(lǐng)域得到了廣泛的應用。在數(shù)字信號處理中，數(shù)字濾波器占有極其重要的地位，它具有精度高、可靠性好、靈活性大等特點?，F(xiàn)代數(shù)字濾波器可以用軟件或硬件兩種方式來實現(xiàn)。軟件方式實現(xiàn)的優(yōu)點是可以通過濾波器參數(shù)的改變?nèi)フ{(diào)整濾波器的性能。MATLAB是一種面向科學和工程計算的語言，它集數(shù)值分析、矩陣運算、信號處理和圖形顯示于一體，具有編程效率高、調(diào)試手段豐富、擴充能力強等

11、特點。MATLAB的信號處理工具箱具有強大的函數(shù)功能，它不僅可以用來設(shè)計數(shù)字濾波器，還可以使設(shè)計達到最憂化，是數(shù)字濾波器設(shè)計的強有力工具。 2.2 FIR數(shù)字濾波器的基本概念 R濾波器的數(shù)學表達式為：式中，N是FIR濾波器的抽頭數(shù)，x(n)表示第n時刻的輸入樣本；h(i)是FIR濾波器的第i級抽頭系數(shù)。普通的直接型FIR濾波器結(jié)構(gòu)如圖1所示。FIR濾波器實質(zhì)上是一個分節(jié)的延遲線，把每一節(jié)的輸出加權(quán)累加，便得到濾波器的輸出。對于FIR濾波器，幅度上只需滿足以下兩個條件之一，就能構(gòu)成線性相位FIR濾波器。h(n)=h(N-1-n) (2)h(n)=-h(N-1-n) (3)(沒有式一？)式（2）稱

12、為第一類線性相位的幅度條件（偶對稱），式（3）稱為第二類線性相位的幅度條件（奇對稱）。2.2.1 FIR濾波器的設(shè)計 根據(jù)沖激響應的時域特性，數(shù)字濾波器可分為無限長沖激響應濾波器（IIR）和有限長沖激響應濾波器（FIR）。FIR的突出優(yōu)點是：系統(tǒng)總是穩(wěn)定的、易于實現(xiàn)線性相位、允許設(shè)計多通帶（或多阻帶）濾波器，但與IIR相比，在滿足同樣阻帶衰減的情況下需要的階數(shù)較高。濾波器的階數(shù)越高，占用的運算時間越多，因此在滿足指標要求的情況下應盡量減少濾波器的階數(shù)。FIR濾波器的基本結(jié)構(gòu)可以理解為一個分節(jié)的延時線，把每一節(jié)的輸出加權(quán)累加，可得到濾波器的輸出。FIR濾波器的沖激響應h(n)是有限長的，數(shù)學上M

13、階FIR濾波器可以表示為：FIR濾波器的設(shè)計問題實質(zhì)上是確定能滿足所要求的轉(zhuǎn)移序列或脈沖響應的常數(shù)的問題，設(shè)計方法主要有窗函數(shù)法、頻率采樣法和等波紋最佳逼近法等。2.2.2 窗函數(shù)設(shè)計法的步驟窗函數(shù)設(shè)計法是一種通過截短和計權(quán)的方法使無限長非因果序列成為有限長脈沖響應序列的設(shè)計方法。通常在設(shè)計濾波器之前，應該先根據(jù)具體的工程應用確定濾波器的技術(shù)指標。在大多數(shù)實際應用中，數(shù)字濾波器常常被用來實現(xiàn)選頻操作，所以指標的形式一般為在頻域中以分貝值給出的相對幅度響應和相位響應。用窗函數(shù)法設(shè)計FIR濾波器的步驟如下：（1）根據(jù)過渡帶寬及阻帶衰減要求，選擇窗函數(shù)的類型并估計窗口長度N（或階數(shù)M=N-1）。窗函

14、數(shù)類型可根據(jù)最小阻帶衰減AS獨立選擇，因為窗口長度N對最小阻帶衰減AS沒有影響。在確定窗函數(shù)類型以后，可根據(jù)過渡帶寬小于給定指標確定所擬用的窗函數(shù)的窗口長度N。設(shè)待求濾波器的過渡帶寬為，它與窗口長度N近似成反比。窗函數(shù)類型確定后，其計算公式也確定了，不過這些公式是近似的，得出的窗口長度還要在計算中逐步修正。原則是在保證阻帶衰減滿足要求的情況下，盡量選擇較小的N。在N和窗函數(shù)類型確定后，即可調(diào)用MATLAB中的窗函數(shù)求出窗函數(shù)wd(n)。（2）根據(jù)待求濾波器的理想頻率響應求出理想單位脈沖響應hd(n)。如果給出待求濾波器的頻率響應為Hd(ej)，則理想的單位脈沖響應可以用下面的傅里葉反變換式求出

15、：在一般情況下，hd(n)是不能用封閉公式表示的，需要采用數(shù)值方法表示。從=0到=2采樣N點，采用離散傅里葉反變換(IDFT)即可求出。（3）計算濾波器的單位脈沖響應h(n)。它是理想單位脈沖響應和窗函數(shù)的乘積，即h(n)=hd(n)·wd(n)，在MATLAB中用點乘命令表示為h=hd·wd。（4）驗算技術(shù)指標是否滿足要求。為了計算數(shù)字濾波器在頻域中的特性，可調(diào)用freqz子程序，如果不滿足要求，可根據(jù)具體情況，調(diào)整窗函數(shù)類型或長度，直到滿足要求為止。使用窗函數(shù)法設(shè)計時要滿足以下兩個條件：（1）窗譜主瓣盡可能地窄，以獲得較陡的過渡帶；（2）盡量減少窗譜的最大旁瓣的相對幅度

16、，也就是使能量盡量集中于主瓣，減小峰肩和紋波，進而增加阻帶的衰減。根據(jù)工程經(jīng)驗，給定的濾波器指標參數(shù)一般為通帶截止頻率p、阻帶截止頻率s、實際通帶波動Rp和最小阻帶衰減As。窗函數(shù)設(shè)計的經(jīng)驗公式為：在實際工程中常用的窗函數(shù)有五種，即矩形窗、三角窗、漢寧窗、海明窗和凱澤窗。這些窗函數(shù)在MATLAB中分別用boxcar、triang、hanning、hamming、kaiser實現(xiàn)，它們之間的性能比較如表1所示。表1 5種窗函數(shù)性能比較窗類型旁瓣峰值主瓣峰值最小阻帶衰減矩形窗13dB4/M21dB三角窗25dB8/M25dB漢寧窗31dB8/M44dB海明窗41dB8/M53dB凱澤窗57dB12

17、/M74dB2 .3 MATLAB環(huán)境下的設(shè)計實例2.3.1 高通濾波器的設(shè)計用窗函數(shù)設(shè)計高通濾波器，性能指標如下：通帶截止頻率s=0.2，阻帶截止頻率p=0.3，實際通帶波動Rp=0.25dB,最小阻帶衰減As=70dB。分析：從表1可以看出凱澤窗能提供74dB的最小阻帶衰減，所以選用凱澤窗進行設(shè)計，程序主要部分如下：As=70;s=0.2*;p=0.3*tr_width=p-s;%計算過渡帶寬M=ceil(As-7.95)*2*/(14.36*tr_width)+1)+1;按凱澤窗計算濾波器長度disp(濾波器的長度為,num2str(M);beta=0.1102*(As-8.7); %計

18、算凱澤窗的值n=0:1:M-1;disp(線性相位斜率為,num2str(beta); w_kai=(kaiser(M,beta); %求凱澤窗函數(shù) c=(s+p)/2; hd=ideal_lp(,M)-ideal_lp(c,M); %求理想脈沖響應 h=hd*w_kai; %設(shè)計的脈沖響應為理想脈沖響應與窗函數(shù)乘積 db,mag,pha,grd,=freqz_m(h,1); delta_=2*/1000; Rp=-(min(db(p/delta_+1:1:501); disp(實際通帶波動為,num2str(Rp);%以下為作圖程序 As=-round(max(db(1:1:s/delta_

19、+1); disp(最小阻帶衰減為,num2str(As);subplot(1,1,1);subplot(2,2,1);stem(n,hd);title(理想脈沖響應);axis(0 M-1 -0.4 0.8);ylabel(hd(n);subplot(2,2,2);stem(n,w_kai);title(凱澤窗); axis(0 M-1 0 1.1);ylabel(wd(n);subplot(2,2,3);stem(n,h);title(實際脈沖響應);axis(0 M-1 -0.4 0.8);xlabel(n);ylabel(h(n);subplot(2,2,4);plot(/,db);t

20、itle(幅度響應/dB);axis(0 1 -100 10);grid;xlabel(以為單位的頻率);ylabel(分貝數(shù)/dB);程序運行結(jié)果如圖1所示。實際通帶波動為0.04369,最小阻帶衰減為70，濾波器長度為89，線性相位斜率為6.7553,符合設(shè)計要求。（下圖是否能通過仿真完成？）2.3.2 低通濾波器的設(shè)計用窗函數(shù)設(shè)計低通濾波器，性能指標如下：通帶截止頻率p=0.1，阻帶截止頻率s=0.25，實際通帶波動Rp=0.10dB,最小阻帶衰減As=40dB。分析：從表1可以看出，漢寧窗、海明窗和凱澤窗能提供大于40dB的最小阻帶衰減。但漢寧窗的旁瓣峰值較小，而主瓣寬度和海明窗一樣。

21、可以使濾波器的階數(shù)較少，所以選用漢寧窗進行設(shè)計，程序主要部分如下：p=0.10*;s=0.25*;tr_width=s-p; %計算過渡帶寬M=ceil(6.6*/tr_width)+1; %按漢寧窗計算濾波器長度disp(濾波器的長度為,num2str(M);n=0:M-1;c=(s+p)/2; %截止頻率取為兩邊緣頻率的平均值hd=ideal_lp(c,M); %求理想脈沖響應w_han=(hanning(M); %求漢寧窗函數(shù)h=hd*w_han; %設(shè)計的脈沖響應為理想脈沖響應與窗函數(shù)乘積db,mag,pha,grd,=freqz_m(h,1);%以下為作圖語句delta_=2*/10

22、00;Rp=-(min(db(1:1: p/delta_+1);disp(實際通帶波動為,num2str(Rp); %以下為作圖程序As=-round(max(db(s/delta_+1:1:501);disp(最小阻帶衰減為,num2str(As);subplot(221)stem(n,hd);title(理想沖擊響應),axis(0 M-1 -0.1 0.3);ylabel(hd(n);subplot(222)stem(n,w_han);title(漢寧窗),axis(0 M-1 0 1.1);ylabel(wd(n);subplot(223)stem(n,h);title(實際沖擊響應)

23、, axis(0 M-1 -0.1 0.3);xlabel(n);ylabel(h(n);subplot(224);plot(/,db);title(幅度響應（db);axis(0 1 -100 10),grid;xlabel(以為單位的頻率);ylabel(分貝數(shù));仿真結(jié)果如圖2所示。實際通帶波動為0.076565,最小阻帶衰減為44，濾波器長度為67，符合設(shè)計要求。與其他高級語言的程序設(shè)計相比，MATLAB環(huán)境下可以更方便、快捷地設(shè)計出具有嚴格線性相位的FIR濾波器，節(jié)省大量的編程時間，提高編程效率，且參數(shù)的修改也十分方便，還可以進一步進行優(yōu)化設(shè)計。相信隨著版本的不斷提高，MATLAB在

24、數(shù)字濾波器技術(shù)中必將發(fā)揮更大的作用。同時，用MATLAB計算有關(guān)數(shù)字濾波器的設(shè)計參數(shù)，如H(z)、h(n)等，對于數(shù)字濾波器的硬件實現(xiàn)也提供了一條簡單而準確的途徑和依據(jù)。第三章 IIR數(shù)字濾波器設(shè)計過程及方法 3.1 IIR數(shù)字濾波器的基本特點 1.IIR數(shù)字濾波器的系統(tǒng)函數(shù)可以寫成封閉函數(shù)的形式。2.IIR數(shù)字濾波器采用遞歸型結(jié)構(gòu)，即結(jié)構(gòu)上帶有反饋環(huán)路。IIR濾波器運算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運算組成，可以組合成直接型、正準型、級聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式，都具有反饋回路。由于運算中的舍入處理，使誤差不斷累積，有時會產(chǎn)生微弱的寄生振蕩。 3.IIR數(shù)字濾波器在設(shè)計上可以借助成熟的模

25、擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現(xiàn)成的設(shè)計數(shù)據(jù)或圖表可查，其設(shè)計工作量比較小，對計算工具的要求不高。在設(shè)計一個IIR數(shù)字濾波器時，我們根據(jù)指標先寫出模擬濾波器的公式，然后通過一定的變換，將模擬濾波器的公式轉(zhuǎn)換成數(shù)字濾波器的公式。 4.IIR數(shù)字濾波器的相位特性不好控制，對相位要求較高時，需加相位校準網(wǎng)絡(luò)。 在MATLAB下設(shè)計IIR濾波器可使用Butterworth函數(shù)設(shè)計出巴特沃斯濾波器，使用Cheby1函數(shù)設(shè)計出契比雪夫I型濾波器，使用Cheby2設(shè)計出契比雪夫II型濾波器，使用ellipord函數(shù)設(shè)計出橢圓濾波器。下面主要介紹前兩個函數(shù)的使用。 與FIR濾波器的設(shè)計

26、不同，IIR濾波器設(shè)計時的階數(shù)不是由設(shè)計者指定，而是根據(jù)設(shè)計者輸入的各個濾波器參數(shù)（截止頻率、通帶濾紋、阻帶衰減等），由軟件設(shè)計出滿足這些參數(shù)的最低濾波器階數(shù)。在MATLAB下設(shè)計不同類型IIR濾波器均有與之對應的函數(shù)用于階數(shù)的選擇。 IIR單位響應為無限脈沖序列FIR單位響應為有限的;iir幅頻特性精度很高，不是線性相位的，可以應用于對相位信息不敏感的音頻信號上； fir幅頻特性精度較之于iir低，但是線性相位，就是不同頻率分量的信號經(jīng)過fir濾波器后他們的時間差不變。這是很好的性質(zhì)。 另外有限的單位響應也有利于對數(shù)字信號的處理，便于編程，用于計算的時延也小，這對實時的信號處理很重要。 數(shù)字

27、濾波器是具有一定傳輸選擇特性的數(shù)字信號處理裝置,其輸入、輸出均為數(shù)字信號,實質(zhì)上是一個由有限精度算法實現(xiàn)的線性時不變離散系統(tǒng)。它的基本工作原理是利用離散系統(tǒng)特性對系統(tǒng)輸入信號進行加工和變換,改變輸入序列的頻譜或信號波形,讓有用頻率的信號分量通過,抑制無用的信號分量輸出。數(shù)字濾波器和模擬濾波器有著相同的濾波概念,根據(jù)其頻率響應特性可分為低通、高通、帶通、帶阻等類型,與模擬濾波器相比,數(shù)字濾波器除了具有數(shù)字信號處理的固有優(yōu)點外,還有濾波精度高(與系統(tǒng)字長有關(guān))、穩(wěn)定性好(僅運行在0與l兩個電平狀態(tài))、靈活性強等優(yōu)點。數(shù)字濾波器按單位脈沖響應的性質(zhì)可分為無限長單位脈沖響應濾波器IIR和有限長單位脈沖

28、響應濾波器(FIR)兩種。本文介紹(IIR)數(shù)字濾波器的設(shè)計與分析。數(shù)字濾波器技術(shù)指標指標參數(shù)變換模擬濾波器技術(shù)指標相應模擬濾波器設(shè)計模擬濾波器離散化數(shù)字濾波器3.2 IIR數(shù)字濾波器設(shè)計思路與步驟   IIR 數(shù)字濾波器可用一個n階差分方程y(n)=brx(n-r)+aky(n-k),或用它的Z域系統(tǒng)函數(shù): 對照模擬濾波器的傳遞函數(shù):    不難看出,數(shù)字濾波器與模擬濾波器的設(shè)計思路相仿,其設(shè)計實質(zhì)也是尋找一組系數(shù)b,a,去逼近所要求的頻率響應,使其在性能上滿足預定的技術(shù)要求;不同的是模擬濾波器的設(shè)計是在S平面上用數(shù)學逼近法去尋找近

29、似的所需特性H(S),而數(shù)字濾波器則是在Z平面尋找合適的H(z)。IIR數(shù)字濾波器的單位響應是無限長的,而模擬濾波器一般都具有無限長的單位脈沖響應,因此與模擬濾波器相匹配。由于模擬濾波器的設(shè)計在理論上已十分成熟,因此數(shù)字濾波器設(shè)計的關(guān)鍵是將H(S)H(Z),即,利用復值映射將模擬濾波器離散化。已經(jīng)證明,沖擊響應不變法和雙線性變換法能較好地擔當此任,則在此基礎(chǔ)上,數(shù)字濾波器的設(shè)計就可首先歸結(jié)為模擬濾波器的設(shè)計了。   數(shù)字濾波器的設(shè)計步驟如圖3.1所示。  圖3.1 數(shù)字濾波器設(shè)計步驟3.3 IIR數(shù)字濾波器設(shè)計方法   IIR數(shù)字濾波器的設(shè)計方

30、法有多種,可歸納為下述兩種。1 .傳統(tǒng)設(shè)計方法   根據(jù)前述設(shè)計思路,首先設(shè)計一個模擬原型濾波器(截止頻率為1rad/s的低通濾波器),然后在模擬域(S平面)進行頻率變換,將模擬原形濾波器轉(zhuǎn)換成所需類型(指定截止頻率的低通、高通、帶通、帶阻)的模擬濾波器,再將其數(shù)字離散化,從S平面映射至Z平面,得到所需技術(shù)指標的數(shù)字濾波器。   上述過程中,也可先將模擬原型離散化,得到數(shù)字原型濾波器,繼而在數(shù)字域(Z平面)進行頻率變換,得到所需類型的數(shù)字濾波器。   模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換可在時域進行也可在頻域?qū)崿F(xiàn),時域轉(zhuǎn)換的關(guān)鍵是要使數(shù)字濾波

31、器與模擬濾波器時域響應的采樣值相等,以保持其瞬態(tài)特性不變,常用的是沖擊響應不變法。頻域變換法必須使得數(shù)字濾波器在-范圍內(nèi)的幅頻特性與模擬濾波器在-/T/T 范圍內(nèi)的幅頻特性一致,即保證S平面與z平面上幅頻特性的一一單值對應關(guān)系,常用的是雙線性變換法。2.計算機輔助設(shè)計方法   傳統(tǒng)設(shè)計方法思路清晰,步驟詳盡,可參閱公式、手冊循章而行。但由于計算繁瑣,手工計算大多只能用來進行簡單低階選頻濾波器(如LP,HP,BP及BS等)的設(shè)計。計算機輔助設(shè)計方法是集電路理論、網(wǎng)絡(luò)圖論、數(shù)值分析、矩陣運算、元件建模、優(yōu)化技術(shù)、高級計算機語言等多交叉學科于一身的新領(lǐng)域,它把計算機的快速、高精度

32、、大存儲容量、嚴格的邏輯判斷和優(yōu)良的數(shù)據(jù)處理能力與人的思維創(chuàng)造能力充分結(jié)合起來,極大地簡化了數(shù)字濾波器的設(shè)計過程。在優(yōu)秀科技應用軟件MATLAB的信號處理工具箱中,提供了一整套模擬、數(shù)字濾波器的設(shè)計命令和運算函數(shù),方便準確,簡單易行,使得設(shè)計人員除了可按上述傳統(tǒng)設(shè)計步驟快速地進行較復雜高階選頻濾波器的計算、分析外,還可通過原型變換法直接進行各種典型數(shù)字濾波器設(shè)計,即應用MATLAB設(shè)計工具從模擬原型直接變換成滿足原定頻域指標要求的數(shù)字濾波器。3 .4各種設(shè)計方法的MATLAB實現(xiàn)3.4.1 基于模擬低通原型的MATLAB實現(xiàn)通過模擬低通原型濾波器進行數(shù)字帶通設(shè)計程序：fp=480 %模擬低通通

33、帶上限頻率fs=520 %模擬低通阻帶下限頻率wp=2*pi*fp %模擬低通通帶上限角頻率ws=2*pi*fs %模擬低通阻帶下限角頻率rp=3 %通帶波動rs=20 %阻帶衰減%巴特沃斯模擬低通原型濾波器設(shè)計n,wn=buttord(wp,WS,rp ,rs,'s')z,p,k=buttap(n) %模擬低通原型零、極點系數(shù)b1,a1=zp2tf(z,p,k) %零、極點系數(shù)轉(zhuǎn)換為傳遞函數(shù)%巴特沃斯模擬低通原型濾波器頻率響應hl,w1=freqs(b1,a1)mag1=abs(h1)模擬低通原型濾波器幅頻特性曲線subplot(221);semilogx(w1,mag1)f

34、w=40 %模擬帶通濾波器帶寬頻率bw=2*pi*fw %模擬帶通濾波器帶寬角頻率由模擬低通原型變換為模擬帶通濾波器b2,a2=lp2bp(b1,a1,wn,bw) %模擬帶通濾波函數(shù)系數(shù)巴特沃斯模擬帶通濾波器頻率響應h2,w2= freqs(b2,a2)mag2 =abs(1l2)模擬帶通濾波器幅頻特性曲線(db)subplot(222);plot(20*log10(mag2)沖擊響應不變法進行離散化設(shè)計fo=2000 %采樣頻率bz,az=impinvar(b2,a2,2000) %數(shù)字帶通濾波函數(shù)系數(shù)巴特沃斯型數(shù)字帶通濾波器頻率響應hz,w= freqz(bz,az)magz=abs(h

35、z)phz=unwrap(angle(hz)subplot(223);plot(magz) %數(shù)字帶通濾波器幅頻特性曲線subplot(224);flot(plot) %數(shù)字帶通濾波器相頻特性曲線 3.4.2 基于合適類型模擬濾波器的MATLAB實現(xiàn)通過合適類型模擬濾波器進行數(shù)字帶通設(shè)計程序fp= 480,520;fs=450,550 %模擬通帶、阻帶頻率wp=480,520*pi*2 %模擬通帶角頻率ws=450,550*pi*2 %模擬阻帶角頻率rp=3;rs=20 %通帶波動、阻帶衰減巴特沃斯型模擬帶通濾波器設(shè)計n,wn=buttord (wp,ws,rp,rs,'s&

36、#39;)b,a=butter(n,wn,'s') %模擬帶通濾波函數(shù)系數(shù)巴特沃斯型模擬帶通濾波器頻率響應ha,w= freqs(b,a)ma=abs(ha);pha=unwrap(angle(ha)subplot(421);plot(w/(2*pi),ma) %模擬幅頻曲線subplot(423);plot(w/(2 pi),pha) %模擬相頻曲線沖擊響應不變法進行離散化設(shè)計fo=5000 %采樣頻率bn,an=impinvar(b,a,5000) %數(shù)字帶通濾波函數(shù)系數(shù)巴特沃斯型數(shù)字帶通濾波器頻率響應hz,w=freqz(bn,an)mz=abs(hz);phz=unwrap(angle(hz)subplot(422);plot(w,mz) %數(shù)字濾波器幅頻曲線subplot(424);plot(w,phz) %數(shù)字濾波器相頻曲線hi=impz(bn,an) %數(shù)字濾波器沖擊響應subplot(425),plot(hi) %沖擊響應曲線n=0:300;t=n/foxl=2*square(2*pi*500*t) %500Hz方波信號subplot(426);plot(x1) %500Hz方波波形yi=conv(hi,x1) %時域卷積輸出subplot(427);plot(yi) %卷積輸