IIR數(shù)字濾波器的設計

1、吉林建筑大學電氣與電子信息工程學院數(shù)字信號處理課程設計報告設計題目： IIR數(shù)字濾波器的設計 專業(yè)班級： 信科121 學生姓名： 學 號： 指導教師： 王蓉暉 王超 設計時間： 2015.1.122015.1.16 教師評語：成績 評閱教師 日期 目 錄 摘要.第1章 緒論.11.1設計的作用、目的.11.2設計任務及要求.11.3設計內(nèi)容.1第2章 IIR數(shù)字低通波器的原理.22.1數(shù)字低通濾波器的設計原理.22.2 IIR數(shù)字濾波器采用方法的原理.3第3章 IIR數(shù)字濾波器仿真分析.83.1設計步驟.83.2數(shù)字低通濾波器MATLAB編程及幅頻特性曲線.113.3優(yōu)缺點.14總結(jié).15參考

2、文獻.16摘 要 數(shù)字濾波器，又名“無限脈沖響應數(shù)字濾波器”，或“遞歸濾波器”。遞歸濾波器，也就是IIR數(shù)字濾波器，IIR數(shù)字濾波器在設計上可以借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現(xiàn)成的設計數(shù)據(jù)或圖表可查，其設計工作量比較小，對計算工具的要求不高。在設計一個IIR數(shù)字濾波器時，我們根據(jù)指標先寫出模擬濾波器的公式，然后通過一定的變換，將模擬濾波器的公式轉(zhuǎn)換成數(shù)字濾波器的公式。 一般認為具有無限的脈沖響應。 熟悉巴特沃斯模擬濾波器的頻率特性掌握用沖激響應不變法、雙線性變換法設計IIR濾波器的原理及方法熟悉相應的MATLAB編程。第1章 緒 論1.1設計的作用、目的 掌

3、握脈沖響應不變法和雙線性變換法設計IIR數(shù)字濾波器的具體方法和原理，熟悉雙線性變換法和脈沖響應不變法設計低通、帶通IIR數(shù)字濾波器的計算機編程。觀察雙線性變換法和脈沖響應不變法設計的數(shù)字濾波器的頻域特性，了解雙線性變換法和脈沖響應不變法的特點和區(qū)別。 課程設計是理論學習的延伸，是掌握所學知識的一種重要手段，對于貫徹理論聯(lián)系實際、提高學習質(zhì)量、塑造自身能力等于有特殊作用。本次課程設計一方面通過MATLAB仿真設計內(nèi)容，使我們加深對理論知識的理解，同時增強其邏輯思維能力，另一方面對課堂所學理論知識作一個總結(jié)和補充。1.2設計任務及要求1.2.1設計任務： 已知通帶截止頻率,通帶最大衰減,阻帶截止頻

4、率，阻帶最小衰減，T=1ms，按照以上技術(shù)要求，用脈沖響應不變法和雙線性變換法設計巴特沃斯數(shù)字低通濾波器，并觀察所設計數(shù)字濾波器的幅頻特性曲線。1.2.2通過課程設計各環(huán)節(jié)的實踐，應達到如下要求： （1）掌握雙線性變換法及脈沖響應不變法設計IIR數(shù)字濾波器以及窗函數(shù)法設計FIR數(shù)字濾波器的原理、具體方法及計算機編程； （2）觀察雙線性變換法、脈沖響應不變法及窗函數(shù)法設計的濾波器的頻域特性，了解各種方法的特點； （3）用MATLAB畫出三種方法設計數(shù)字濾波器的幅頻特性曲線，記錄帶寬和衰減量，檢查結(jié)果是否滿足要求。1.3設計內(nèi)容 已知通帶截止頻率,通帶最大衰減,阻帶截止頻率，阻帶最小衰減，T=1m

5、s，按照以上技術(shù)要求，用脈沖響應不變法和雙線性變換法設計巴特沃斯數(shù)字低通濾波器，并觀察所設計數(shù)字濾波器的幅頻特性曲線。第2章 IIR數(shù)字低通波器的原理2.1數(shù)字低通濾波器的設計原理 數(shù)字濾波器可以分為兩大類：一類是經(jīng)典濾波器，即一般的濾波器特點是輸入信號中有用的頻率成分和希望濾去的頻率成分各占不同的頻率帶，通過一個合適的選頻濾波器達到綠寶锝目的，這種濾波器主要有維納濾波器爾曼濾波器，自適應濾波器等。 數(shù)字濾波器是一個離散時間系統(tǒng)。應用數(shù)字濾波器處理模擬信號時，首先須對輸入模擬信號進行限帶、抽樣和模數(shù)轉(zhuǎn)換。數(shù)字濾波器輸入信號的抽樣率應大于被處理信號帶寬的兩倍。數(shù)字濾波器的頻率響應具有以抽樣頻率為

6、間隔的周期重復特性，且以折疊頻率（即二分之一抽樣頻率點）呈鏡像對稱。為得到模擬信號，數(shù)字濾波器處理的輸出數(shù)字信號須經(jīng)數(shù)模轉(zhuǎn)換、平滑。 數(shù)字濾波器具有高精度、高可靠性、可程控改變特性或復用、便于集成等優(yōu)點。數(shù)字濾波器在語聲信號處理、圖像信號處理、醫(yī)學生物信號處理以及其他應用領(lǐng)域（如通信、雷達、聲納、儀器儀表和地震勘探等）都得到了廣泛的應用。 數(shù)字濾波器有低通、高通、帶通、帶阻和全通等類型。它可以是時不變的或時變的、因果的或非因果的、線性的或非線性的。如果數(shù)字濾波器的內(nèi)部參數(shù)不隨時間而變化，則稱為時不變的，否則為時變的。如果數(shù)字濾波器在某一給定時刻的響應與在此時刻以后的激勵無關(guān)，則稱為因果的，否則

7、為非因果的。如果數(shù)字濾波器對單一或多個激勵信號的響應滿足線性條件，則稱為線性的，否則為非線性的。應用最廣的是線性、時不變數(shù)字濾波器。 數(shù)字濾波器可以按所處理信號的維數(shù)分為一維、二維或多維數(shù)字濾波器。一維數(shù)字濾波器處理的信號為單變量函數(shù)序列，例如時間函數(shù)的抽樣值。二維或多維數(shù)字濾波器處理的信號為兩個或多個變量函數(shù)序列。例如，二維圖像離散信號是平面坐標上的抽樣值。 與模擬濾波器相同，數(shù)字濾波器從功能上看可分為低通、高通、帶通和帶阻幾類。因為它們的單位脈沖響應是非因果且無限長的，所以實際上理想濾波器是不可能實現(xiàn)的。與模擬濾波器不同的是數(shù)字濾波器的傳輸函數(shù)都是以2為周期的，濾波器的低通頻帶處于2的整數(shù)

8、倍處，而高頻頻帶處于的奇數(shù)倍附近。理想濾波器就是一個讓輸入信號中的某些有用頻譜分量無任何變化的通過，同時又能完全抑制另外那些不需要的成分的具有某種選擇性的器件、網(wǎng)絡或以計算機硬件支持的計算機程序。根據(jù)對不同信號的處理可將濾波器分為模擬濾波器和數(shù)字濾波器兩種。模擬濾波器和數(shù)字濾波器的概念相同，只是信號的形式和實現(xiàn)濾波方法不同。數(shù)字濾波器是指輸入輸出都是數(shù)字信號的濾波器。濾波器的濾波原理局勢根據(jù)信號與噪聲占據(jù)不同的頻帶，將噪聲的頻率放在濾波器的阻帶中，而由于阻帶的響應為零，這樣就濾去了噪聲。一個理想濾波器將。是一個無法實現(xiàn)的非因果系統(tǒng)，我們只能用一個穩(wěn)定的因果系統(tǒng)函數(shù)去逼近工程需要的確定的性能要求

9、。數(shù)字濾波器根據(jù)其沖激響應函數(shù)的時域特性，可分為兩種，即無限長沖激響應(IIR)濾波器和有限長沖激響應(FIR)濾波器。IIR濾波器的特征是，具有無限持續(xù)時間沖激響應。這種濾波器一般需要用遞歸模型來實現(xiàn)，因而有時也稱之為遞歸濾波器。FIR濾波器的沖激響應只能延續(xù)一定時間，在工程實際中可以采用遞歸的方式實現(xiàn)，也可以采用非遞歸的方式實現(xiàn)。數(shù)字濾波器的設計方法有多種，如雙線性變換法、窗函數(shù)設計法、插值逼近法和Chebyshev逼近法等等。隨著MATLAB軟件尤其是MATLAB的信號處理工作箱的不斷完善，不僅數(shù)字濾波器的計算機輔助設計有了可能，而且還可以使設計達到最優(yōu)化。（1）確定指標在設計一個濾波器

10、之前，必須首先根據(jù)工程實際的需要確定濾波器的技術(shù)指標。在很多實際應用中，數(shù)字濾波器常常被用來實現(xiàn)選頻操作。因此，指標的形式一般在頻域中給出幅度和相位響應。幅度指標主要以兩種方式給出。第一種是絕對指標。它提供對幅度響應函數(shù)的要求，一般應用于FIR濾波器的設計。第二種指標是相對指標。它以分貝值的形式給出要求。在工程實際中，這種指標最受歡迎。對于相位響應指標形式，通常希望系統(tǒng)在通頻帶中具有線性相位。運用線性相位響應指標進行濾波器設計具有如下優(yōu)點：只包含實數(shù)算法，不涉及復數(shù)運算；不存在延遲失真，只有固定數(shù)量的延遲；長度為N的濾波器(階數(shù)為N-1)，計算量為N/2數(shù)量級。 (2)逼近確定了技術(shù)指標后，就

11、可以建立一個目標的數(shù)字濾波器模型。通常采用理想的數(shù)字濾波器模型。之后，利用數(shù)字濾波器的設計方法，設計出一個實際濾波器模型來逼近給定的目標。(3)性能分析和計算機仿真上兩步的結(jié)果是得到以差分或系統(tǒng)函數(shù)或沖激響應描述的濾波器。根據(jù)這個描述就可以分析其頻率特性和相位特性，以驗證設計結(jié)果是否滿足指標要求；或者利用計算機仿真實現(xiàn)設計的濾波器，再分析濾波結(jié)果來判斷。（4）濾波器的性能指標我們在進行濾波器設計時，需要確定其性能指標。一般來說，濾波器的性能要求往往以頻率響應的幅度特性的允許誤差來表征。以低通濾波器特性為例，頻率響應有通帶、過渡帶及阻帶三個范圍。在通帶內(nèi)： 在阻帶中： 其中為通帶截止頻率，Ap為

12、通帶誤差，為阻帶誤差。與模擬濾波器類似，數(shù)字濾波器按頻率特性劃分為低通、高通、帶通、帶阻、全通等類型，由于數(shù)字濾波器的頻率響應是周期性的，周期為。 2.2 IIR數(shù)字濾波器采用方法的原理利用模擬濾波器成熟的理論及其設計方法來設計IIR數(shù)字低通濾波器是常用的方法。2.2.1設計過程按照數(shù)字濾波器技術(shù)指標要求設計一個過度模擬低通濾波器，在按照一定的轉(zhuǎn)換關(guān)系將轉(zhuǎn)換成數(shù)字低通濾波器的系統(tǒng)函數(shù)。由此可見，設計的關(guān)鍵問題就是找到這種轉(zhuǎn)換關(guān)系，將s平面上的轉(zhuǎn)換成z平面上的。為了保證轉(zhuǎn)換后的穩(wěn)定且滿足技術(shù)指標要求，對轉(zhuǎn)換關(guān)系提出兩點要求： （1）因果穩(wěn)定的模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，仍是因果穩(wěn)定的。

13、 （2）數(shù)字濾波器的頻率響應模仿模擬濾波器的頻響特性，s平面的虛軸映射為z平面的單位圓，響應的頻率之間呈線性關(guān)系。 一個無限脈沖響應系統(tǒng)的傳遞函數(shù)具有可選擇極點和可選擇零點（而有線脈沖響應系統(tǒng)只具有可選擇零點以及位于原點的固定極點）。模擬系統(tǒng)同樣具有極點和零點，甚至經(jīng)常只有極點（多項式濾波器或全幾點濾波器）。由此得出結(jié)論是，通過把PN圖從s域映像變換到z域中，可以用一個無限脈沖（IIR）系統(tǒng)對模擬濾波器進行數(shù)字仿真。 目前，IIR數(shù)字濾波器設計最通用的方法是借助于模擬濾波器的設計方法。模擬濾波器設計已經(jīng)有了一套相當成熟的方法，它不但有完整的設計公式，而且還有較為完整的圖表供查詢，因此，充分利用

14、這些已有的資源將會給數(shù)字濾波器的設計帶來很大方便。 將高通、帶通或者帶阻數(shù)字濾波器的技術(shù)指標先轉(zhuǎn)化為低通濾波器的技術(shù)指標，然后按上述步驟(2)設計出模擬低通濾波器H(s)，再由沖擊響應不變法或雙線性變換將H(s)轉(zhuǎn)換為所需的H(z)。 s - z 映射的方法有：沖激響應不變法、階躍響應不變法、雙線性變換法等。實際中有三種廣泛應用的濾波器原型，既巴特沃斯低通濾波器、切比雪夫低通濾波器、橢圓低通濾波器。2.2.2巴特沃斯低通濾波器這種濾波器的特征是其通帶和阻帶都有平坦的幅度響應。N階低通濾波器的平方幅度響應的表達式為：通帶： 使信號通過的頻帶阻帶：抑制噪聲通過的頻帶過渡帶：通帶到阻帶間過渡的頻率范

15、圍c ：截止頻率。 過渡帶為零理想濾波器 阻帶|H(j )|=0 通帶內(nèi)幅度|H(j)|=cons. H(j)的相位是線性的其中N為濾波器的階數(shù)，c是邊緣頻率(單位為rad/s)。取N=1,2,10，可得到平方幅頻響應。N增加，通帶和阻帶的近似性越好，過渡帶越陡。通帶內(nèi)，分母/c<1, ( /c)2N<1，A(2)1。過渡帶和阻帶，/c>1， ( /c)2N>1, 增加， A(2)快速減小。=c, ,幅度衰減，相當于3db衰減點。振幅平方函數(shù)的極點 可見，巴特沃斯濾波器的振幅平方函數(shù)有2N個極點，它們均勻?qū)ΨQ地分布在|S|=c的圓周上。考慮到系統(tǒng)的穩(wěn)定性，知DF的系統(tǒng)函

16、數(shù)是由S平面左半部分的極點(SP3，SP4，SP5)組成的，它們分別為：系統(tǒng)函數(shù)為：歸一化的三階BF： 如果要還原的話，則有2.2.3脈沖響應不變法設計IIR數(shù)字濾波器 利用模擬濾波器成熟的理論及其設計方法來設計IRR數(shù)字低通濾波器是通常的方法。 從濾波器的脈沖響應出發(fā)，使數(shù)字濾波器的單位脈沖響應模仿模擬濾波器的沖擊響應，即h(n)是ha(t)的采樣值。設T為采樣周期，變換過程®¾¾¾®¾¾®¾利用模擬濾波器來設計數(shù)字濾波器，也就是使數(shù)字濾波器能模仿模擬濾波器的特性，這種模仿可以從不同的角度出發(fā)。 脈沖響

17、應不變法是從濾波器的脈沖響應出發(fā)，使數(shù)字濾波器的單位脈沖響應序列h(n)模仿模擬濾波器的沖激響應ha(t)，即將ha(t)進行等間隔采樣，使h(n)正好等于ha(t)的采樣值，滿足h(n)=ha(nT)式中,T是采樣周期。 如果令Ha(s)是h其中T為采樣間隔。如果以Ha(s)及H(z)分別表示ha(t)的拉氏變換及h(n)的Z變換，則a(t)的拉普拉斯變換，H(z)為h(n)的Z變換，利用采樣序列的Z變換與模擬信號的拉普拉斯變換的關(guān)系 脈沖響應不變法設計中，模擬頻率與數(shù)字頻率之間的轉(zhuǎn)換關(guān)系是線性的()。同時，它可以保持脈沖響應不變，。因此，這一方法往往用于低通時域數(shù)字濾波器設計及相應的模擬系

18、統(tǒng)數(shù)字仿真設計。a.頻率混疊現(xiàn)象 數(shù)字濾波器的頻率響應是模擬濾波器頻率響應的周期延拓。只有當模擬濾波器的頻率響應是限帶的，且?guī)抻谡郫B頻率以內(nèi)時，才能使數(shù)字濾波器的頻率響應在折疊頻率以內(nèi)重現(xiàn)模擬濾波器的頻率響應，而不產(chǎn)生混疊失真。即|w|< 但是，任何一個實際的模擬濾波器頻率響應都不是嚴格限帶的（非理想）， 變換后就會產(chǎn)生周期延拓分量的頻譜交疊，即產(chǎn)生頻率響應的混疊失真。這時數(shù)字濾波器的頻響就不同于原模擬濾波器的頻響，而是有一定的失真。當模擬濾波器的頻率響應在折疊頻率以上處衰減越大、越快時，變換后頻率響應混疊失真就越小。這時采用脈沖響應不變法設計的數(shù)字濾波器才能得到良好的效果。 如圖2-

19、1所示這種頻譜混疊現(xiàn)象會使設計出的數(shù)字濾波器在附近的頻率響應特性程度不同的偏離模擬濾波器在附近的頻率特性，嚴重時使數(shù)字濾波器不滿足給定指標。為此，希望設計的濾波器是帶限濾波器。高通和帶阻濾波器不適合這種方法設計。綜上所述，脈沖響應不變法的優(yōu)點是頻率坐標變換是線性的，如果不考慮頻率混疊現(xiàn)象，用這種方法設計的數(shù)字濾波器會很好地重現(xiàn)原模擬濾波器的頻率特性。另一個優(yōu)點是數(shù)字濾波器的單位脈沖響應完全模仿模擬濾波器的單位脈沖響應，時域特性逼近好。其缺點是會產(chǎn)生頻譜混疊現(xiàn)象，使數(shù)字濾波器的頻響偏離模擬濾波器的頻響。脈沖響應不變法適合低通、帶通濾波器的設計，不適合高通、帶阻濾波器的設計。圖2-1 脈沖響應不變

20、法的頻率混疊現(xiàn)象2.2.4采用雙線性變換法設計IIR數(shù)字濾波器 脈沖響應不變法的缺點主要是會產(chǎn)生頻譜混疊現(xiàn)象，使數(shù)字濾波器偏離模擬濾波器的頻響特性。為了克服這一缺點，采用非線性頻率壓縮方法，將整個模擬頻率軸壓縮。=2/T*tan(w/2)T是采樣區(qū)間。轉(zhuǎn)到z平面得先確定數(shù)字低通濾波器的技術(shù)指標：帶通邊界頻率帶通最大衰減阻帶截止頻率阻帶最小衰減將邊界頻率轉(zhuǎn)換S平面的虛軸單值地映射于z平面的單位圓上，s平面的左半平面完全映射到z平面的單位圓內(nèi)。雙線性變換不存在頻率混疊問題。  在MATLAB中，可用函數(shù)bilinear實現(xiàn)從模擬濾波器到數(shù)字濾波器的雙線性變換映射。  雙線性變換

21、是一種非線性變換，這種非線性引起的幅頻特性畸變可通過預畸變得到校正。第3章 IIR數(shù)字濾波器仿真分析3.1設計步驟3.1.1脈沖響應不變法設計步驟 （1）頻率歸一化問題 沖響應不變法是實現(xiàn)模擬濾波器數(shù)字化的一種直觀而常用的方法。它特別適合 于對濾波器的時域特性有一定要求的場合。 脈沖響應不變法的核心原理是通過對連續(xù)函數(shù)ha(t)等間隔采樣得到離散序列ha(nT)。令h(n)= ha(nT) ，T為采樣間隔。它是一種時域上的轉(zhuǎn)換方法。 一個模擬濾波器的傳遞函數(shù)可以用有理分式表達式表示為：通過反拉普拉斯變換我們就可以得到它的沖激相應：脈沖響應不變法就是要保證脈沖響應不變，即：上式的沖激相應序列作變

22、換，就可以得到數(shù)字濾波器的傳遞函數(shù)： 一般來說，的分母多項式階次總是大于分子多項式的階次。假定 的沒有多重極點，則式（5-7）就可分解為：式中，均為復數(shù)，是的極點。其拉氏變換為 一般來說，的分母多項式階次總是大于分子多項式的階次。數(shù)字濾波器的傳遞函數(shù)H(z)經(jīng)過合并簡化，成為一般形式的有理分式傳遞函數(shù)在討論采樣序列z變換與模擬信號拉氏變換之間關(guān)系的有關(guān)章節(jié)中，我們已經(jīng)知道每一個s平面上寬度為/T的水平條帶將重迭映射到z平面上。因此脈沖響應不變法將s平面映射到z平面，不是一個簡單的一一對應的關(guān)系。對于高采樣頻率(T小)的情況，數(shù)字濾波器在頻域可能有極高的增益。為此我們采用利用模擬濾波器來設計II

23、R低通數(shù)字濾波器是常用的方法，稱之為模擬一數(shù)字轉(zhuǎn)換法。（2）利用在MATLAB設計IIR數(shù)字濾波器可分以下幾步來實現(xiàn) 按一定規(guī)則將數(shù)字濾波器的技術(shù)指標轉(zhuǎn)換為模擬低通濾波器的技術(shù)指標；根據(jù)轉(zhuǎn)換后的技術(shù)指標使用濾波器階數(shù)函數(shù)，確定濾波器的最小階數(shù)N；利用最小階數(shù)N產(chǎn)生模擬低通濾波原型；利用截止頻率把模擬低通濾波器原型轉(zhuǎn)換成模擬帶通原型；利用沖激響應不變法或雙線性不變法把模擬濾波器轉(zhuǎn)換成數(shù)字濾波器。（3）相關(guān)參數(shù) ,T=1ms N=63.1.2雙線性變換法設計步驟 根據(jù)以上IIR數(shù)字濾波器設計方法，采用雙線性變換法基于MATLAB設計一個巴特沃斯數(shù)字低通濾波器，帶截止頻率,通帶最大衰減,阻帶截止頻率

24、，阻帶最小衰減，T=1ms根據(jù)任務,確定性能指標:在設計低通濾波器之前,首先根據(jù)工程實際的需要確定濾波器的設計指標: 通帶截止頻率 阻帶截止頻率 通帶最大衰減 阻帶最小衰減 =1.2566 =1.8849根據(jù)公式=2/T*tan(w/2) 進行預修正，將數(shù)字濾波器指標轉(zhuǎn)換成模擬濾波器設計指標，對雙線性變換法一般T=1ms=(2/T)*tan(/2)=(2/T)*tan(/2)=1dB=25dB;N=63.2數(shù)字低通濾波器MATLAB編程及幅頻特性曲線3.2.1MATLAB語言編程（1）脈沖響應不變法：Fs=1000; %采樣頻率為周期倒數(shù)Wp=400*pi;Ws=600*pi; %設置歸一化通

25、帶和阻帶截止頻率Ap=1;As=25; %設置通帶最大和最小衰減N,Wc=buttord(Wp,Ws,Ap,As,'s'); %調(diào)用butter函數(shù) B,A=butter(N,Wc,'s'); %設計巴特沃斯濾波器W=(0:pi:1000); %指定一段頻率值D,C=impinvar(B,A,Fs); %調(diào)用脈沖響應不變法W=(0:0.001:pi)Hz=freqz(D,C,W); %返回頻率響應 plot(W,abs(Hz)/abs(Hz(1); %繪出幅頻特性曲線grid;xlabel('頻率');ylabel('幅值dB')

26、;（2）雙線性變換法 wp=0.4*pi; %通帶截止頻率 ws=0.6*pi; %阻帶截止頻率 Ap=1; %通帶最大衰減 As=25; %阻帶最小衰減 T=0.001;fs=1/T; Wp=2/T*tan(wp/2);Ws=2/T*tan(ws/2); N,Wn=buttord(Wp,Ws,Ap,As,'s'); %巴特沃斯濾波器 b,a=butter(N,Wn,'s'); numd,dend=bilinear(b,a,fs); %雙線性變換法 w=linspace(0,pi,500); h=freqz(numd,dend,w); plot(w,20*log

27、10(abs(h); xlabel('頻率'); ylabel('幅度/dB'); grid; title('雙線性變換法');3.2.2幅頻特性曲線（1） 脈沖響應不變法圖3-1Ap=1,As=25,Fs=1000,Wc=1.369268611959261e+03,Wp=1.256637061435917e+03Ws=1.884955592153876e+03A=1,7.885303666046043e+03,3.108900695287958e+07,8.000392885207492e+10,1.475920968172557e+14,2.

28、020932255451204e+17,2.053890440162035e+20,1.496408320076508e+23,7.116050029613046e+25,1.691990325921431e+28B=0,0,0,0,0,0,0,0,0,1.691990325921431e+28C=1,-1.919940662665034,2.532351531105067,-2.205257262678870,1.386786124116544,-0.630894421427259,0.204467001672835,-0.044995465406303,0.006049838435517,

29、-3.762323452047202e-04D=-7.767084753140807e-13,1.677227578906053e-04,0.015295223041337,0.099548037576103,0.144380467187304,0.061131885122289,0.007470107777323,1.973357622301677e-04,3.656905171283142e-07,0(2)雙線性變換法圖3-2Ap=1,As=25, wp=1.256637061435917,ws=1.884955592153876a=1,6.584973310037336e+03,2.16

30、8093674695203e+07,4.525576614274616e+10,6.297650306260731e+13,5.555907639115123e+16,2.450764030490121e+19 b=0,0,0,0,0,0,2.450764030490121e+193.3優(yōu)缺點 用脈沖響應不變法設計IIR數(shù)字低通濾波器的缺點是會產(chǎn)生頻譜混疊現(xiàn)象，使數(shù)字濾波器的頻響偏離模擬濾波器的頻響特性。產(chǎn)生的原因是模擬低通濾波器不是帶限于折疊頻率/T，在離散化（采樣）后產(chǎn)生了頻譜混疊，再通過映射關(guān)系使數(shù)字濾波器在=附近形成頻譜混疊。為了克服這個缺點，可以采用非線性頻率壓縮方法，將整個模擬頻率

31、軸壓縮到/T之間,轉(zhuǎn)換到Z平面上。 用脈沖響應不變法設計IIR數(shù)字定濾波器的優(yōu)點是頻率變換關(guān)系是線性的，既=T，如果不出現(xiàn)頻譜混疊現(xiàn)象，用這種方法設計的數(shù)字濾波器會很好地重現(xiàn)原模擬濾波器的頻響特性。另外一個優(yōu)點是數(shù)字濾波器的單位脈沖響應完全模仿模擬濾波器的單位沖激響應波形，時域特性逼近好。但是，有限階的模擬濾波器不可能是理想帶限的，所以，脈沖響應不變法的最大缺點是會產(chǎn)生不同程度的頻譜混疊失真，其適用于低通、帶通濾波器的設計，不適用于高通、帶阻濾波器的設計。 用雙線性變換法設計IIR數(shù)字低通濾波器的缺點是之間的非線性關(guān)系，是數(shù)字濾波器頻響曲線不能保真的模仿模擬濾波器的頻響曲線形狀。這種非線性影響

32、的實質(zhì)問題是：如果刻度是均勻的，則其映像的刻度是不均勻的，而是隨增加越來越密。因此，如果模擬濾波器的頻響具有片段常數(shù)特性，則主要是數(shù)字濾波器頻響特性曲線的轉(zhuǎn)折點頻率值與模擬濾波器特性轉(zhuǎn)折點的頻率之成非線性關(guān)系。當然，對于不是片段常數(shù)的相應特性仍有非線性失真。因此，雙線性變換法適合片斷常數(shù)特性的濾波器的設計。 用雙線性變換法設計IIR數(shù)字低通濾波器的優(yōu)點是可用簡單的公式代替Ha(S)直接轉(zhuǎn)換成H（z），這是改變換法的優(yōu)點。但當階數(shù)稍高時，將H（z）整理成需要的形式，也不是一件簡單的工作?？偨Y(jié)通過這次數(shù)字信號處理的設計，讓我感受到“溫故而知新”是的重要性，我們總是以為知識學過了，就不用再花費時間去

33、復習，其實，這是錯誤的。我感覺我們在學習知識的時候不但要聽懂，學會，還要把它轉(zhuǎn)化成自己知識寶庫里的一部分，并能夠舉一反三。同時，我還深刻的感覺到，在我們學好本專業(yè)課程的同時，我們還要多去圖書館或者上網(wǎng)查閱一些相關(guān)的知識，而且我們還要自學一些軟件，這樣我們就在無形之中提高了自身的競爭力。比如這次使用的MATLAB 軟件，我們曾在上課的時候?qū)W過，但是由于很久都沒有去再接觸它，我們已經(jīng)有點淡忘，通過本次設計，我們不但對之前的操作進行了復習，而且通過查閱相關(guān)書籍又學到了很多，發(fā)現(xiàn) MATLAB 具有相當強大的矩陣運算和操作功能，其程序編寫也比其他語言要簡單。通過這次設計，使我發(fā)現(xiàn)了自己身上真的有好多不

34、足，希望在以后的不斷學習中能夠得到提高，我們要“學以致用”，并能夠熟練地通過書籍和網(wǎng)絡查找知識，而且最重要的是我們以后要多花費一些時間在自學知識與軟件的操作上。本次課程設計可以說是獲益匪淺。通過查閱了很多資料，了解了許多匯編程序的思想，擴展了自己的視野，不再僅僅局限于書本中幾條簡短的程序，而且更重要的是明白態(tài)度：仔細謹慎，精益求精。 在實際操作過程中犯的一些錯誤還會有意外的收獲，感覺設計很有意思。在具體操作中對這學期所學的技術(shù)的理論知識得到鞏固，平時在做題目過程中，總更多的是被動的通過題目來熟悉再理解，而這次是通過實際操作,用Matlab畫圖，發(fā)現(xiàn)錯誤和問題，然后在理解的基礎上再熟記于心。這是一種飛躍的進步。達到設計的基本目的，也發(fā)現(xiàn)自己的不足之出，在以后的上機中應更加注意，同時體會到科