基于matlab的iir數(shù)字濾波器的設計與仿真畢業(yè)論文

1、<p><b>　　目錄</b></p><p><b>　　摘 要1</b></p><p>　　ABSTRACT2</p><p><b>　　前 言3</b></p><p>　　第1章 濾波器的基本概念4</p><p>　　1

2、.1 濾波器的定義4</p><p>　　1.2 濾波器的發(fā)展4</p><p>　　1.3 濾波器的分類5</p><p>　　第2章 濾波器的原理與設計7</p><p>　　2.1 模擬濾波器的原理與設計7</p><p>　　2.2 數(shù)字濾波器的原理與設計8</p><p>

3、　　第3章 IIR數(shù)字濾波器的設計10</p><p>　　3.1 IIR數(shù)字濾波器簡介10</p><p>　　3.2 IIR數(shù)字濾波器的主要技術指標10</p><p>　　3.3 IIR數(shù)字濾波器的設計過程12</p><p>　　3.3.1用脈沖相應不變法設計IIR數(shù)字濾波器13</p><p>

4、;　　3.3.2雙線性變換法設計IIR數(shù)字濾波器16</p><p>　　第4章 IIR數(shù)字濾波器的實現(xiàn)20</p><p>　　4.1 matlab簡介20</p><p>　　4.1.1 matlab的概況20</p><p>　　4.1.2 matlab的語言特點20</p><p>　　4.1.

5、3 matlab的功能22</p><p>　　4.2 利用模擬原型進行IIR數(shù)字濾波器的設計和仿真23</p><p>　　4.2.1巴特沃斯(Butterworth)和橢圓數(shù)字低通濾波器的設計24</p><p>　　4.2.2切比雪夫（Chebyshev）數(shù)字濾波器的設計30</p><p>　　4.3 IIR數(shù)字濾波器的直接

6、設計法34</p><p>　　第5章 IIR數(shù)字濾波器的綜合設計和仿真40</p><p>　　5.1數(shù)字低通濾波器設計實例與仿真40</p><p>　　5.2數(shù)字帶通濾波器設計實例與仿真42</p><p><b>　　結 論45</b></p><p><b>　　致

7、謝46</b></p><p><b>　　參考文獻47</b></p><p><b>　　畢業(yè)設計小結48</b></p><p><b>　　摘 要</b></p><p>　　在現(xiàn)代通信系統(tǒng)中，由于信號中經?；煊懈鞣N噪聲和干擾，所以很多信號分析都是基于

8、濾波器而進行的，而數(shù)字濾波器是通過數(shù)值運算實現(xiàn)濾波，具有處理精度高、穩(wěn)定、靈活、不存在阻抗匹配問題等優(yōu)點，可以實現(xiàn)模擬濾波器無法實現(xiàn)的特殊濾波功能。數(shù)字濾波器根據(jù)其沖激響應函數(shù)的時域特性，可分為兩種，即無限長沖激響應(IIR)數(shù)字濾波器和有限長沖激響應(FIR)數(shù)字濾波器。實現(xiàn)IIR濾波器的階次較低，所用的存儲單元較少，效率高，精度高，而且能夠保留一些模擬濾波器的優(yōu)良特性，因此應用很廣。Matlab軟件以矩陣運算為基礎，把計算、可視化及

9、程序設計有機融合到交互式工作環(huán)境中，并且為數(shù)字濾波的研究和應用提供了一個直觀、高效、便捷的利器。尤其是Matlab中的信號處理工具箱使各個領域的研究人員可以直觀方便地進行科學研究與工程應用。</p><p>　　本文根據(jù)模擬濾波器的設計原理，提出了IIR數(shù)字濾波器的設計方法，并在MATLAB環(huán)境下實現(xiàn)了IIR 數(shù)字濾波器的設計和仿真。其主要內容概括為：首先對濾波器的原理和設計進行了介紹；接著描述了IIR數(shù)字濾波器

10、的基本概念，其中包括系統(tǒng)的描述、系統(tǒng)的傳遞函數(shù)、系統(tǒng)的模型；接著簡單介紹MATLAB，并對數(shù)字濾波器在MATLAB環(huán)境下如何實現(xiàn)進行了介紹；重點描述了IIR數(shù)字濾波器的設計過程，最后對IIR濾波器進行仿真。</p><p>　　關鍵詞： MATLAB， IIR數(shù)字濾波器， 模擬濾波器</p><p><b>　　ABSTRACT</b></p>&

11、lt;p>　　In modern communication system, because the signal is often mixed with various noise and interference, so a lot of signal analysis are based on the filter and the digital filter, and through numerical arithmeti

12、c to realize filtering, has high processing precision, stable, flexible, there is no problem of impedance matching and other advantages, can realize analog filter can achieve special filtering function.Digital filter bas

13、ed on the impulse response function in the time domain characteristics</p><p>　　KEY WORDS: MATLAB, IIR digital filter, Dimulation filt</p><p><b>　　前 言</b></p><p>　　濾波

14、是信號處理的基礎，濾波運算是信號處理中的基本運算，濾波器的設計也就相應成為數(shù)字信號處理的最基本問題之一。濾波器可廣義的理解為一個信號選擇系統(tǒng)，它讓某些信號成分通過又阻止或衰減另一些成分。在更多的情況下，濾波器可理解為選頻系統(tǒng)，如低通、高通、帶通、帶阻。濾波器可分為三種：模擬濾波器、采樣濾波器和數(shù)字濾波器。模擬濾波器可以是由RLC構成的無源濾波器，也可以是加上運放的有源濾波器，是連續(xù)時間系統(tǒng)；采樣濾波器由電阻、電容、電荷轉移器件、運放等組

15、成，屬于離散時間系統(tǒng)，幅度連續(xù)；數(shù)字濾波器由加法器、乘法器、存儲延遲單元、時鐘脈沖發(fā)生器和邏輯單元等數(shù)字電路構成，精度高，穩(wěn)定性好，不存在阻抗匹配問題，可以時分復用。</p><p>　　數(shù)字濾波器是具有一定傳輸選擇特性的數(shù)字信號處理裝置，其輸入、輸出均為數(shù)字信號，實質上是一個由有限精度算法實現(xiàn)的線性時不變離散系統(tǒng)。它的基本工作原理是利用離散系統(tǒng)特性對系統(tǒng)輸入信號進行加工和變換，改變輸入序列的頻譜或信號波形，讓有

16、用頻率的信號分量通過，抑制無用的信號分量輸出。數(shù)字濾波器和模擬濾波器有著相同的濾波概念，數(shù)字濾波器根據(jù)其單位脈沖響應的性質可分為無限長單位脈沖響應濾波器IIR和有限長單位脈沖響應濾波器(FIR)兩種，按其頻率響應特性可分為低通、高通、帶通、帶阻等類型,與模擬濾波器相比，數(shù)字濾波器除了具有數(shù)字信號處理的固有優(yōu)點外，還有靈活性強、穩(wěn)定性好、濾波精度高等優(yōu)點。</p><p>　　第1章 濾波器的基本概念</p&

17、gt;<p><b>　　濾波器的定義</b></p><p>　　濾波器是一種用來消除干擾雜訊的器件，凡是可以使信號中特定的頻率成分通過，而極大地衰減或抑制其他頻率成分的裝置或系統(tǒng)都稱之為濾波器，相當于頻率“篩子”。 濾波器的功能就是允許某一部分頻率的信號順利的通過，而另外一部分頻率的信號則受到較大的抑制，它實質上是一個選頻電路。</p><p>　　

18、濾波器中，把信號能夠通過的頻率范圍，稱為通頻帶或通帶；反之，信號受到很大衰減或完全被抑制的頻率范圍稱為阻帶；通帶和阻帶之間的分界頻率稱為截止頻率；理想濾波器在通帶內的電壓增益為常數(shù)，在阻帶內的電壓增益為零；實際濾波器的通帶和阻帶之間存在一定頻率范圍的過渡帶。</p><p><b>　　濾波器的發(fā)展</b></p><p>　　在近代電信設備和各類控制系統(tǒng)中，濾波器應

19、用極為廣泛；在所有的電子部件中，使用最多，技術最為復雜的要算濾波器了。濾波器的優(yōu)劣直接決定產品的優(yōu)劣，所以，對濾波器的研究和生產歷來為各國所重視。 </p><p>　　1917年美國和德國科學家分別發(fā)明了LC濾波器，美國第一個多路復用系統(tǒng)的出現(xiàn)，20世紀50年代無源濾波器日趨成熟，60年代起由于計算機技術、集成工藝和材料工業(yè)的發(fā)展，濾波器發(fā)展上了一個新臺階，并且朝著低功耗、高精度、小體積、多功能、穩(wěn)定可靠和價廉

20、的方向努力，70年代以后的主攻方向為小體積、多功能、高精度、穩(wěn)定可靠導致RC有源濾波器、數(shù)字濾波器、開關電容濾波器和電荷轉移器等各種濾波器的飛速發(fā)展，到70年代后期，上述幾種濾波器的單片集成已被研制出來并得到應用，80年代，致力于各類新型濾波器的研究，努力提高性能并逐漸擴大應用范圍，90年代至現(xiàn)在主要致力于把各類濾波器應用于各類產品的開發(fā)和研制。當然，對濾波器本身的研究仍在不斷進行。 </p><p>　　我國于

21、50年代后廣泛使用濾波器，當時主要用于話路濾波和報路濾波。經過半個世紀的發(fā)展，我國濾波器在研制、生產和應用等方面已踏入國際發(fā)展步伐，但由于缺少專門研制機構，集成工藝和材料工業(yè)跟不上來，使得我國許多新型濾波器的研制應用與國際發(fā)展有一段距離。</p><p><b>　　濾波器的分類</b></p><p>　　安縣好處理的類型分，濾波器分為模擬濾波器和數(shù)字濾波器。模擬濾

22、波器由電阻，電容，電感等原器件構成；而數(shù)字濾波器則可由數(shù)字電路實現(xiàn)或軟件實現(xiàn)。按脈沖響應來分類，數(shù)字濾波器分為IIR和FIR，即無限沖激響應濾波器和有限沖激響應濾波器；其中IIR網(wǎng)絡中有反饋回路，F(xiàn)IR網(wǎng)絡中沒有反饋回路。</p><p><b>　　從小的方面分：</b></p><p>　　a.按所處理的信號分為模擬濾波器和數(shù)字濾波器兩種。</p>

23、<p>　　b.按所通過信號的頻段分為低通、高通、帶通和帶阻濾波器四種。</p><p>　　低通濾波器：它允許信號中的低頻或直流分量通過，抑制高頻分量或干擾和噪聲。 </p><p>　　高通濾波器：它允許信號中的高頻分量通過，抑制低頻或直流分量。</p><p>　　帶通濾波器：它允許一定頻段的信號通過，抑制低于或高于該頻段的信號、干擾和噪聲。 <

24、;/p><p>　　帶阻濾波器：它抑制一定頻段內的信號，允許該頻段以外的信號通過。</p><p>　　圖1.1 由通過信號的頻段劃分的四種濾波器</p><p>　　c.按所采用的元器件分為無源和有源濾波器兩種: </p><p>　　無源濾波器：僅由無源元件(R、L 和C)組成的濾波器，它是利用電容和電感元件的電抗隨頻率的變化而變化的原理構成

25、的。這類濾波器的優(yōu)點是：電路比較簡單，不需要直流電源供電，可靠性高；缺點是：通帶內的信號有能量損耗，負載效應比較明顯，使用電感元件時容易引起電磁感應，當電感L較大時濾波器的體積和重量都比較大，在低頻域不適用。 </p><p>　　有源濾波器：由無源元件(一般用R和C)和有源器件(如集成運算放大器）組成。這類濾波器的優(yōu)點是：通帶內的信號不僅沒有能量損耗，而且還可以放大，負載效應不明顯，多級相聯(lián)時相互影響很小，利用

26、級聯(lián)的簡單方法很容易構成高階濾波器，并且濾波器的體積小、重量輕、不需要磁屏蔽(由于不使用電感元件）；缺點是：通帶范圍受有源器件(如集成運算放大器）的帶寬限制，需要直流電源供電，可靠性不如無源濾波器高，在高壓、高頻、大功率的場合不適用。</p><p>　　第2章 濾波器的原理與設計</p><p>　　濾波器是一種用來過濾時間離散信號的數(shù)字系統(tǒng)，通過對抽樣數(shù)據(jù)進行數(shù)學處理來達到頻域濾波的目

27、的。可以設計系統(tǒng)的頻率響應，讓它滿足一定的要求，從而對通過該系統(tǒng)的信號的某些特定的頻率成分進行過濾，這就是濾波器的基本原理。如果系統(tǒng)是一個連續(xù)系統(tǒng)，則濾波器稱為模擬濾波器。如果系統(tǒng)是一個離散系統(tǒng)，則濾波器稱為數(shù)字濾波器。</p><p>　　2.1 模擬濾波器的原理與設計</p><p>　　模擬濾波器可以分為無源和有源濾波器，濾波器是一種能使有用信號順利通過而同時對無用頻率信號進行抑制（

28、或衰減）的電子裝置。工程上常用它來做信號處理、數(shù)據(jù)傳送和抑制干擾等。以往主要采用無源元件R、L和C組成模擬濾波器，六十年代以來，集成運放獲得了迅速地發(fā)展，由它和R、C組成的有源濾波器，具有不用電感、體積小、重量輕等優(yōu)點。此外，由于集成運放的開環(huán)電壓增益和輸入阻抗均很高，而輸出阻抗又很低，而且，由其構成的有源濾波器還具且一定的電壓放大和緩沖作用。因此，基于放大器和R、C構成的有源濾波器應用日益廣泛。模擬濾波器的設計，就是用模擬系統(tǒng)的系統(tǒng)函

29、數(shù)去逼近所要求的理想特性。標準的模擬低通濾波器的設計公式有巴特沃思和切比雪夫等，它們都是根據(jù)幅度平方函數(shù)來確定的。</p><p>　　為逼近理想低通濾波器，其模擬理想低通濾波器的幅度特性可用幅度平方函數(shù)表示，即</p><p><b>　?。?.1）</b></p><p>　　式中，為所設計的模擬濾波器的系統(tǒng)函數(shù)，它是s的有理函數(shù)；是其穩(wěn)態(tài)

30、響應，即濾波器頻率特性為濾波器的穩(wěn)態(tài)振幅特性。</p><p>　　由幅度平方函數(shù)確定的方法是：</p><p>　?。╝） 在中，令，得到；</p><p>　　(b) 將的有理式進行分解，得到零點和極點。如果系統(tǒng)函數(shù)是最小相位函數(shù)，則s平面左半平面的零點和極點都屬于，而任何在虛軸上的極點和零點都是偶次的，其中一半屬于</p><p>　　

31、(c) 根據(jù)具體情況，比較和的幅度特性，確定增益常數(shù)，這樣就完全確定了。</p><p>　　在模擬濾波器中，低通濾波器的設計是最基本的，高通，帶通和帶阻濾波器等可以用頻帶變換方法由低通濾波器轉變得到。</p><p>　　2.2 數(shù)字濾波器的原理與設計</p><p>　　所謂數(shù)字濾波器，是指輸入、輸出均為數(shù)字信號，通過一定運算關系改變輸入信號所含頻率成分的相對比

32、例或者濾除某些頻率成分的器件。數(shù)字濾波器和模擬濾波器相比，因為信號的形式和實現(xiàn)濾波的方法不同，數(shù)字濾波器具有比模擬濾波器精度高、穩(wěn)定、體積小、重量輕、靈活、不要求阻抗匹配等優(yōu)點。</p><p>　　數(shù)字濾波器從功能上分類:可分為低通濾波器、高通濾波器、帶通濾波器、帶阻濾波器。從濾波器的網(wǎng)絡結構或者從單位脈沖響應分類:可分為IIR濾波器(即無限長單位沖激響應濾波器)和FIR濾波器(即有限長單位沖激響應濾波器)。I

33、IR數(shù)字濾波器可用一個n階差分方程表示：</p><p><b>　?。?.2）</b></p><p>　　或用它的Z域系統(tǒng)函數(shù):</p><p><b>　?。?.3）</b></p><p>　　對照模擬濾波器的傳遞函數(shù):</p><p><b>　?。?.4

34、）</b></p><p>　　可以看出,設計數(shù)字濾波器的思路與模擬濾波器相仿,設計實質都是尋找一組系數(shù){b,a},去逼近所求的頻率響應,以便在性能上滿足預定的技術要求；不同的是數(shù)字濾波器是在Z平面尋找合適的H(Z), 模擬濾波器則是在S平面用數(shù)學逼近法尋找近似的所需特性H(S)。由于模擬濾波器的設計在理論上已十分成熟,因此數(shù)字濾波器設計的關鍵是將H(S)變換至H(Z),即利用復值映射將模擬濾波器離散

35、化。</p><p>　　已經證明,脈沖響應不變法和雙線性變換法能較好地擔當此任。IIR數(shù)字濾波器的單位響應是無限長的,而模擬濾波器一般都具有無限長的單位脈沖響應,故與模擬濾波器相類似。</p><p>　　第3章 IIR數(shù)字濾波器的設計</p><p>　　3.1 IIR數(shù)字濾波器簡介</p><p>　　數(shù)字濾波器在數(shù)字信號處理的各種應

36、用中發(fā)揮著十分重要的作用。它是通過對采樣數(shù)據(jù)信號進行數(shù)學運算處理來達到濾波的目的。其中無限沖擊響應數(shù)字濾波器也稱IIR是采用對離散采樣數(shù)據(jù)作差分方程運算來進行濾波。IIR數(shù)字濾波器的優(yōu)點在于可以利用模擬濾波器設計的結果，然后采用雙線性變換法或沖激響應不變法將模擬濾波器轉換成數(shù)字濾波器，而模擬濾波器的設計方便簡單、有大量的圖表可查。IIR數(shù)字濾波器具有以下特點：</p><p>　?。?） IIR數(shù)字濾波器的系統(tǒng)函

37、數(shù)可以寫成封閉函數(shù)的形式。  </p><p>　?。?） IIR數(shù)字濾波器采用遞歸型結構，即結構上帶有反饋環(huán)路。IIR濾波器運算結構通常由延時、乘以系數(shù)和相加等基本運算組成，可以組合成直接型、正準型、級聯(lián)型、并聯(lián)型四種結構形式，都具有反饋回路。由于運算中的舍入處理，使誤差不斷累積，有時會產生微弱的寄生振蕩。  </p><p>　?。?） IIR數(shù)字濾波器在設計上可以

38、借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現(xiàn)成的設計數(shù)據(jù)或圖表可查，其設計工作量比較小，對計算工具的要求不高。在設計一個IIR數(shù)字濾波器時，我們根據(jù)指標先寫出模擬濾波器的公式，然后通過一定的變換，將模擬濾波器的公式轉換成數(shù)字濾波器的公式。  </p><p>　?。?） IIR數(shù)字濾波器的相位特性不好控制，對相位要求較高時，需加相位校準網(wǎng)絡。</p><p&g

39、t;　　3.2 IIR數(shù)字濾波器的主要技術指標</p><p>　　在設計IIR數(shù)字濾波器之前，首先了解數(shù)字濾波器的技術指標，才能將其轉化成滿足要求的模擬濾波器技術指標，再將模擬濾波器轉化成數(shù)字濾波器。數(shù)字濾波器的主要技術指標為：</p><p>　　(1) 特征頻率參數(shù)</p><p>　　濾波器的頻率參數(shù)主要有：</p><p>　?、?/p>

40、通帶截止頻率：為通帶與過渡帶的邊界點，在該點信號增益到規(guī)定的下限。②阻帶截頻：為阻帶于過渡帶的邊界點，在該點信號衰減到規(guī)定的下限。③轉折頻率：為 信號功率衰減到(約為3)時的頻率，但在多數(shù)情況下也常以作為通帶或阻帶截頻。④當電路沒有損耗時，固有頻率，就是其諧振頻率。</p><p><b>　　(2) 增益和衰減</b></p><p>　　濾波器在通帶內的增益并不是

41、常數(shù)。①對低通濾波器來說，通帶增益一般指時的增益；對于高通而言，通帶增益指時的增益；對于帶通而言，通帶增益則指中心頻率處的增益。②對帶阻濾波器而言，應給出阻帶衰減，衰減定義為增益的倒數(shù)。③通帶增益變化指通帶內各點增益的最大變化量，如果以為單位，則指增益值的變化量。</p><p>　　(3) 阻尼系數(shù)與品質因數(shù)</p><p>　　阻尼系數(shù)是表征濾波器對角頻率為信號的阻尼作用，同時是濾波器

42、中表示能量衰減的一項重要指標，也是與傳遞函數(shù)的極點實部大小相關的一項系數(shù)。阻尼系數(shù)的倒數(shù)1/稱為品質因數(shù)Q，是評價帶通和帶阻濾波器頻率選擇特性的一個重要指標，它可以表示為。其中為帶通或帶阻的帶寬，為中心頻率，在多數(shù)情況下中心頻率等于固有頻率。</p><p><b>　　(4) 靈敏度</b></p><p>　　濾波器的濾波電路是由很多元件構成的，每個元件的參數(shù)值變

43、化都會影響它的性能。當濾波器的某一個性能指標對另一個性能指標變化的靈敏度記為，定義為：</p><p><b>　　(3.1)</b></p><p>　　靈敏度是濾波電路設計中的一個重要參數(shù)，但該靈敏度與測量儀器或電路系統(tǒng)靈敏度概念不同，該靈敏度越小，標志著電路容錯能力越大，穩(wěn)定性也就越高。</p><p><b>　　(5) 群延

44、時函數(shù)</b></p><p>　　當濾波器幅頻特性滿足設計要求時，為保證輸出信號失真不超過允許范圍，對其相頻特性也相應提出了一定要求。在濾波器設計中，常用群延時函數(shù)評價信號經濾波后相位失真程度，越接近常數(shù)，信號相位失真越小。</p><p>　　3.3 IIR數(shù)字濾波器的設計過程</p><p>　　IIR數(shù)字濾波器的設計一般有兩種方法：一個是借助模

45、擬濾波器的設計方法進行。其設計步驟是，先設計模擬濾波器，再按照某種方法轉換成數(shù)字濾波器。這種方法比較容易一些，因為模擬濾波器的設計方法已經非常成熟，不僅有完整的設計公式，還有完善的圖表供查閱；另外一種直接在頻率或者時域內進行，由于需要解聯(lián)立方程，設計時需要計算機做輔助設計。其設計步驟是：先設計過渡模擬濾波器得到系統(tǒng)函數(shù)，然后將按某種方法轉換成數(shù)字濾波器的系統(tǒng)函數(shù)。另外，還有一些典型的優(yōu)良濾波器類型可供我們使用。</p>&

46、lt;p>　　為了保證轉換后的穩(wěn)定且滿足技術指標要求，對轉換關系提出兩點要求:</p><p>　　（1）因果穩(wěn)定的模擬濾波器轉換成數(shù)字濾波器，仍是因果穩(wěn)定的。</p><p>　?。?）數(shù)字濾波器的頻率相應模仿模擬濾波器的頻響特性，s平面的虛軸映射為z平面的單位圓，相應的頻率之間呈線性關系。</p><p>　　利用模擬濾波器成熟的理論設計IIR數(shù)字濾波器的

47、過程是：</p><p>　　(1)確定數(shù)字低通濾波器的技術指標：通帶邊界頻率、通帶最大衰減、阻帶截止頻率、阻帶最小衰減。</p><p>　　(2)將數(shù)字低通濾波器的技術指標轉換成相應的模擬低通濾波器的技術指標。</p><p>　　(3)按照模擬低通濾波器的技術指標設計過渡模擬低通濾波器。</p><p>　　(4)用所選的轉換方法，將模

48、擬濾波器轉換成數(shù)字低通濾波器系統(tǒng)函數(shù)。</p><p>　　IIR數(shù)字濾波器的設計流程圖如下：</p><p>　　圖3.1 IIR數(shù)字濾波器的設計步驟流程圖</p><p>　　成熟的數(shù)字濾波器設計方法主要有脈沖響應不變法和雙線性變換法。</p><p>　　3.3.1用脈沖相應不變法設計IIR數(shù)字濾波器</p><p

49、><b>　　一 設計原理</b></p><p>　　利用模擬濾波器來設計數(shù)字濾波器，也就是使數(shù)字濾波器能模仿模擬濾波器的特性，這種模仿可以從不同的角度出發(fā)。脈沖響應不變法是從濾波器的脈沖響應出發(fā)，使數(shù)字濾波器的單位脈沖響應序列h(n)模仿模擬濾波器的沖激響應ha(t)，即將ha(t)進行等間隔采樣，使h(n)正好等于ha(t)的采樣值，滿足h(n)=ha(nT)，式中T是采樣周期

50、。</p><p>　　如果令Ha(s)是ha(t)的拉普拉斯變換，H(z)為h(n)的Z變換，利用采樣序列的Z變換與模擬信號的拉普拉斯變換的關系得</p><p><b>　　(3.2)</b></p><p>　　則可看出，脈沖響應不變法將模擬濾波器的S平面變換成數(shù)字濾波器的Z平面，這個從s到z的變換z=esT是從S平面變換到Z平面的標準變

51、換關系式。</p><p>　　圖3.2 脈沖響應不變法的映射關系</p><p>　　由（3.2）式，數(shù)字濾波器的頻率響應和模擬濾波器的頻率響應間的關系為</p><p><b>　　(3.3)</b></p><p>　　這就是說，數(shù)字濾波器的頻率響應是模擬濾波器頻率響應的周期延拓。正如采樣定理所討論的，只有當模擬濾

52、波器的頻率響應是限帶的，且?guī)抻谡郫B頻率以內時，即</p><p><b>　　(3.4)</b></p><p>　　才能使數(shù)字濾波器的頻率響應在折疊頻率以內重現(xiàn)模擬濾波器的頻率響應，而不產生混疊失真，即</p><p><b>　　(3.5)</b></p><p>　　但是，任何一個實際的模

53、擬濾波器頻率響應都不是嚴格限帶的，變換后就會產生周期延拓分量的頻譜交疊，即產生頻率響應的混疊失真。這時數(shù)字濾波器的頻響就不同于原模擬濾波器的頻響，而帶有一定的失真。當模擬濾波器的頻率響應在折疊頻率以上處衰減越大、越快時，變換后頻率響應混疊失真就越小。這時，采用脈沖響應不變法設計的數(shù)字濾波器才能得到良好的效果。</p><p>　　圖3.3脈沖響應不變法中的頻響混疊現(xiàn)象</p><p>　　

54、對某一模擬濾波器的單位沖激響應ha(t)進行采樣，采樣頻率為fs，若使fs增加，即令采樣時間間隔（T=1/fs）減小，則系統(tǒng)頻率響應各周期延拓分量之間相距更遠，因而可減小頻率響應的混疊效應。</p><p>　　二 脈沖響應不變法的優(yōu)點</p><p>　　從以上討論可以看出，脈沖響應不變法使得數(shù)字濾波器的單位脈沖響應完全模仿模擬濾波器的單位沖激響應，也就是時域逼近良好，而且模擬頻率Ω和

55、數(shù)字頻率ω之間呈線性關系ω=ΩT。因而，一個線性相位的模擬濾波器（例如貝塞爾濾波器）通過脈沖響應不變法得到的仍然是一個線性相位的數(shù)字濾波器。</p><p>　　脈沖響應不變法的最大缺點是有頻率響應的混疊效應。所以，脈沖響應不變法只適用于限帶的模擬濾波器(例如，衰減特性很好的低通或帶通濾波器)，而且高頻衰減越快，混疊效應越小。至于高通和帶阻濾波器，由于它們在高頻部分不衰減，因此將完全混淆在低頻響應中。如果要對高通

56、和帶阻濾波器采用脈沖響應不變法，就必須先對高通和帶阻濾波器加一保護濾波器，濾掉高于折疊頻率以上的頻率，然后再使用脈沖響應不變法轉換為數(shù)字濾波器。當然這樣會進一步增加設計復雜性和濾波器的階數(shù)。</p><p>　　3.3.2雙線性變換法設計IIR數(shù)字濾波器</p><p><b>　　一 設計原理</b></p><p>　　脈沖響應不變法的主

57、要缺點是產生頻率響應的混疊失真。這是因為從S平面到Ｚ平面是多值的映射關系所造成的。為了克服這一缺點，可以采用非線性頻率壓縮方法，將整個頻率軸上的頻率范圍壓縮到-π/T～π/T之間，再用z=esT轉換到Z平面上。也就是說，第一步先將整個S平面壓縮映射到S1平面的-π/T～π/T一條橫帶里；第二步再通過標準變換關系z=es1T將此橫帶變換到整個Z平面上去。這樣就使S平面與Z平面建立了一一對應的單值關系，消除了多值變換性，也就消除了頻譜混疊現(xiàn)

58、象，映射關系如圖3.4所示。</p><p>　　圖3.4雙線性變換的映射關系</p><p>　　為了將S平面的整個虛軸jΩ壓縮到S1平面jΩ1軸上的-π/T到π/T段上，可以通過以下的正切變換實現(xiàn)</p><p><b>　　(3.6)</b></p><p>　　式中,T仍是采樣間隔。</p><

59、;p>　　當Ω1由-π/T經過0變化到π/T時，Ω由-∞經過0變化到+∞，也即映射了整個jΩ軸。將式（2-6）寫成</p><p><b>　　(3.7)</b></p><p>　　將此關系解析延拓到整個S平面和S1平面，令jΩ=s，jΩ1=s1，則得</p><p><b>　　(3.8)</b></p&g

60、t;<p>　　再將S1平面通過以下標準變換關系映射到Z平面 </p><p><b>　　(3.9)</b></p><p>　　從而得到S平面和Z平面的單值映射關系為：</p><p><b>　　(3.10)</b></p><p><b>　　(3.11)</

61、b></p><p>　　式（3.10）與式（3.10）是S平面與Z平面之間的單值映射關系，這種變換都是兩個線性函數(shù)之比，因此稱為雙線性變換</p><p>　　式（3.6）與式（3.10的雙線性變換符合映射變換應滿足的兩點要求。</p><p>　　首先,把z=ejω，可得</p><p><b>　?。?.12）</

62、b></p><p>　　即S平面的虛軸映射到Z平面的單位圓。</p><p>　　其次，將s=σ+jΩ代入式（2-12），得</p><p><b>　　（3.13）</b></p><p><b>　　因此</b></p><p><b>　?。?.14）

63、</b></p><p>　　由此看出，當σ<0時，|z|<1；當σ>0時，|z|>1。也就是說，S平面的左半平面映射到Z平面的單位圓內，S平面的右半平面映射到Z平面的單位圓外，S平面的虛軸映射到Z平面的單位圓上。因此，穩(wěn)定的模擬濾波器經雙線性變換后所得的數(shù)字濾波器也一定是穩(wěn)定的。</p><p>　　二 雙線性變換法優(yōu)缺點</p>&l

64、t;p>　　雙線性變換法與脈沖響應不變法相比，其主要的優(yōu)點是避免了頻率響應的混疊現(xiàn)象。這是因為S平面與Z平面是單值的一一對應關系。S平面整個jΩ軸單值地對應于Z平面單位圓一周，即頻率軸是單值變換關系。這個關系如式（3.10）所示，重寫如下：</p><p><b>　?。?.15）</b></p><p>　　上式表明，S平面上Ω與Z平面的ω成非線性的正切關系，

65、如圖3.5所示。</p><p>　　由圖3.5看出，在零頻率附近，模擬角頻率Ω與數(shù)字頻率ω之間的變換關系接近于線性關系；但當Ω進一步增加時，ω增長得越來越慢，最后當Ω→∞時，ω終止在折疊頻率ω=π處，因而雙線性變換就不會出現(xiàn)由于高頻部分超過折疊頻率而混淆到低頻部分去的現(xiàn)象，從而消除了頻率混疊現(xiàn)象。</p><p>　　圖3.5雙線性變換法的頻率變換關系</p><p&

66、gt;　　但是雙線性變換的這個特點是靠頻率的嚴重非線性關系而得到的，如式（3.12）及圖3.4所示。由于這種頻率之間的非線性變換關系，就產生了新的問題。首先，一個線性相位的模擬濾波器經雙線性變換后得到非線性相位的數(shù)字濾波器，不再保持原有的線性相位了；其次，這種非線性關系要求模擬濾波器的幅頻響應必須是分段常數(shù)型的，即某一頻率段的幅頻響應近似等于某一常數(shù)（這正是一般典型的低通、高通、帶通、帶阻型濾波器的響應特性），不然變換所產生的數(shù)字濾波器

67、幅頻響應相對于原模擬濾波器的幅頻響應會有畸變，如圖3.6所示:</p><p>　　圖3.6雙線性變換法幅度和相位特性的非線性映射</p><p>　　對于分段常數(shù)的濾波器，雙線性變換后，仍得到幅頻特性為分段常數(shù)的濾波器，但是各個分段邊緣的臨界頻率點產生了畸變，這種頻率的畸變，可以通過頻率的預畸變來加以校正。也就是將臨界模擬頻率事先加以畸變，然后經變換后正好映射到所需要的數(shù)字頻率上。<

68、;/p><p>　　第4章 IIR數(shù)字濾波器的實現(xiàn)</p><p>　　利用MATLAB設計濾波器,可以按照設計要求非常方便地調整設計參數(shù),極大地減輕了設計的工作量,有利于濾波器設計的最優(yōu)化。MATLAB因其強大的數(shù)據(jù)處理功能被廣泛應用于工程計算,其豐富的工具箱為工程計算提供了便利,利用MATLAB信號處理工具箱可以快速有效地設計各種數(shù)字濾波器,設計簡單方便。</p><p

69、>　　4.1 MATLAB簡介</p><p>　　4.1.1 MATLAB的概況</p><p>　　MATLAB是矩陣實驗室（Matrix Laboratory）之意。除具備卓越的數(shù)值計算能力外，它還提供了專業(yè)水平的符號計算，文字處理，可視化建模仿真和實時控制等功能。</p><p>　　MATLAB的基本數(shù)據(jù)單位是矩陣，它的指令表達式與數(shù)學,工程中常用的

70、形式十分相似,故用MATLAB來解算問題要比用C,FORTRAN等語言完成相同的事情簡捷得多.</p><p>　　MATLAB擁有數(shù)百個內部函數(shù)的主包和三十幾種工具包(Toolbox).工具包又可以分為功能性工具包和學科工具包.功能工具包用來擴充MATLAB的符號計算,可視化建模仿真,文字處理及實時控制等功能.學科工具包是專業(yè)性比較強的工具包,控制工具包,信號處理工具包,通信工具包等都屬于此類.</p&g

71、t;<p>　　開放性使MATLAB廣受用戶歡迎.除內部函數(shù)外,所有MATLAB主包文件和各種工具包都是可讀可修改的文件,用戶通過對源程序的修改或加入自己編寫程序構造新的專用工具包.</p><p>　　4.1.2 MATLAB的語言特點</p><p>　　MATLAB最突出的特點就是簡潔。MATLAB用更直觀的，符合人們思維習慣的代碼，代替了C和FORTRAN語言的冗長代

72、碼。MATLAB給用戶帶來的是最直觀，最簡潔的程序開發(fā)環(huán)境。以下簡單介紹一下MATLAB的主要特點：</p><p>　?。?）語言簡潔緊湊，使用方便靈活，庫函數(shù)極其豐富。MATLAB程序書寫形式自由，利用起豐富的庫函數(shù)避開繁雜的子程序編程任務，壓縮了一切不必要的編程工作。由于庫函數(shù)都由本領域的專家編寫，用戶不必擔心函數(shù)的可靠性?？梢哉f，用MATLAB進行科技開發(fā)是站在專家的肩膀上。</p><

73、;p>　?。?）運算符豐富。由于MATLAB是用C語言編寫的，MATLAB提供了和C語言幾乎一樣多的運算符，靈活使用MATLAB的運算符將使程序變得極為簡短。</p><p>　?。?）MATLAB既具有結構化的控制語句（如for循環(huán)，while循環(huán)，break語句和if語句），又有面向對象編程的特性。</p><p>　?。?）程序限制不嚴格，程序設計自由度大。例如，在MATLAB

74、里，用戶無需對矩陣預定義就可使用。</p><p>　?。?）程序的可移植性很好，基本上不做修改就可以在各種型號的計算機和操作系統(tǒng)上運行。</p><p>　?。?）MATLAB的圖形功能強大。在FORTRAN和C語言里，繪圖都很不容易，但在MATLAB里，數(shù)據(jù)的可視化非常簡單。MATLAB還具有較強的編輯圖形界面的能力。</p><p>　　（7）MATLAB的缺

75、點是，它和其他高級程序相比，程序的執(zhí)行速度較慢。由于MATLAB的程序不用編譯等預處理，也不生成可執(zhí)行文件，程序為解釋執(zhí)行，所以速度較慢。</p><p>　?。?）功能強大的工具箱是MATLAB的另一特色。MATLAB包含兩個部分：核心部分和各種可選的工具箱。核心部分中有數(shù)百個核心內部函數(shù)。其工具箱又分為兩類：功能性工具箱和學科性工具箱。功能性工具箱主要用來擴充其符號計算功能，圖示建模仿真功能，文字處理功能以及

76、與硬件實時交互功能。功能性工具箱用于多種學科。而學科性工具箱是專業(yè)性比較強的，如control,toolbox,signl proceessing toolbox,commumnication toolbox等。這些工具箱都是由該領域內學術水平很高的專家編寫的，所以用戶無需編寫自己學科范圍內的基礎程序，而直接進行高、精、尖的研究。</p><p>　　（9）源程序的開放性。開放性也許是MATLAB最受人們歡迎的特

77、點。除內部函數(shù)以外，所有MATLAB的核心文件和工具箱文件都是可讀可改的源文件，用戶可通過對源文件的修改以及加入自己的文件構成新的工具箱。</p><p>　　4.1.3 MATLAB的功能</p><p>　　MATLAB包含的內容非常豐富，功能強大，可以概括為以下幾個方面：</p><p>　?。?）可以在多種操作系統(tǒng)下運行，如DOS、Windows 95/98

78、/2000/2000/NT、Compaq Alpha、LinuxSun Solaris等。</p><p>　?。?）有超過500種的數(shù)學、統(tǒng)計、科學及工程方面的函數(shù)，使用簡單快捷，并且有很強的用戶自定義函數(shù)的能力。</p><p>　?。?）有強大的圖形繪制和可視化功能，可以進行視覺數(shù)據(jù)處理和分析，進行圖形、圖像的顯示及編輯，能夠繪制二維、三維圖形，使用戶可以制作高質量的圖形，從而寫出圖

79、文并茂的文章。</p><p>　　（4）有和用其他高級語言（如C，C++，F(xiàn)ORTRAN，JAVA）編寫的外部程序相接口的能力，也可把MATLAB程序轉換成上述高級語言的子程序。</p><p>　?。?）有從外部文件及外部硬件設備讀入數(shù)據(jù)的能力。</p><p>　?。?）有豐富的網(wǎng)絡資源，從相關的Web網(wǎng)站可以直接獲得全套的MATLAB聯(lián)機幫助文件和說明書的電

80、子文檔，還可以獲得各類技術支持與幫助。</p><p>　?。?）有豐富的工具箱〔toolbox〕。各個領域的專家學者將眾多學科領域中常用的算法編寫為一個個子程序，即m文件，這些m文件包含在一個個工具箱中。其工具箱可以分為兩大類，即功能性工具箱和科學性工具箱。功能性工具箱主要用來擴充MATLAB的符號計算、圖形可視化、建模仿真、文字處理等功能以及與硬件實時交互的功能。學科性工具箱是按學科領域來分類的，如信號處理、

81、控制、通信、神經網(wǎng)絡圖像處理、系統(tǒng)辨識、魯棒控制、模糊邏輯、小波等工具箱。</p><p>　　MATLAB中的信號處理工具箱內容豐富，使用簡便。在數(shù)字信號處理中常用的算法，如FFT，卷積，相關，濾波器設計，參數(shù)模型等，幾乎都只用一條語句即可調用。數(shù)字信號處理常用的函數(shù)有波形的產生、濾波器的分析和設計、傅里葉變換、Z變換等。</p><p>　　MATLAB是一個功能十分強大的系統(tǒng)，是集數(shù)

82、值計算、圖形管理、程序開發(fā)為一體的環(huán)境。除此之外，MATLAB還具有很強的功能擴展能力，與它的主系統(tǒng)一起，可以配備各種各樣的工具箱，以完成一些特定的任務。用戶可以根據(jù)自己的工作任務，開發(fā)自己的工具箱。在國際學術界，MATLAB已經被確認為準確、可靠的科學計算標準軟件。在許多國際一流學術刊物上，（尤其是信息科學刊物），都可以看到MATLAB的應用。 在設計研究單和工業(yè)部門，MATLAB被認作進行高效研究、開發(fā)的首選軟件工具。如美國Nati

83、onal Instruments公司信號測量、分析軟件LabVIEW，Cadence公司信號和通信分析設計軟件SPW等，或者直接建筑在MATLAB之上，或者以MATLAB為主要支撐。 </p><p>　　4.2 利用模擬原型進行IIR數(shù)字濾波器的設計和仿真</p><p>　　用MATLAB進行模擬原型的數(shù)字濾波器的設計，一般步驟如下:</p><p>　?。?）

84、按一定規(guī)則將給出的數(shù)字濾波器的技術指標轉換成模擬低通濾波器的技術指標；</p><p>　　（2）根據(jù)轉換后的技術指標使用濾波器階數(shù)選擇函數(shù)，確定最小階數(shù)N和固有頻率Wn，根據(jù)選用的模擬低通濾波器的類型可分別用: buttord,cheblord，cheb2ord，ellipord等函數(shù)；</p><p>　?。?）運用最小階數(shù)N產生模擬濾波器原型，模擬低通濾波器的創(chuàng)建函數(shù)有:buttap

85、,cheblap,cheb2ap,ellipap,besselap等；</p><p>　?。?）運用固有頻率Wn把模擬低通濾波器原型轉換成模擬低通、高通、帶通、帶阻濾波器，可分別用函數(shù)lp2lp,lp2hp,lp2bp,lp2bs；</p><p>　?。?）運用沖激響應不變法或雙線性變換法把模擬濾波器轉換成數(shù)字濾波器，分別用函數(shù)impinva和bilinear來實現(xiàn)。低通Chebysh

86、evl型數(shù)字濾波器的設計:設計中需要限定其通帶上限臨界頻率Wp,阻帶臨界濾波頻率Ws，在通帶內的最大衰減Rp，阻帶內的最小衰減Rs。</p><p>　　其步驟如圖4.1所示。</p><p>　　圖4.1 數(shù)字濾波器設計步驟</p><p>　　4.2.1巴特沃斯(Butterworth)和橢圓數(shù)字低通濾波器的設計</p><p>　　設

87、計巴特沃斯數(shù)字低通濾波器和橢圓數(shù)字低通濾波器，要求通帶邊界頻fp=2.1kHZ，通帶最大衰減Rp=0.5dB；阻帶邊界頻率fs=8kHZ，阻帶最小衰減Rs=30dB,采樣頻率為Fs=20kHZ。</p><p>　　A 巴特沃思低通濾波器的幅度平方函數(shù)為</p><p><b>　?。?.1）</b></p><p>　　式中，N為正整數(shù)，稱為

88、濾波器的階數(shù)。N值越大，通帶和阻帶的近似特性就越好。在截止頻率處，幅度平方響應為=0處的1/2，相當于幅度響應的3dB衰減點。</p><p><b>　　其系統(tǒng)函數(shù)為 </b></p><p><b>　?。?.2）</b></p><p>　　式中，為歸一化常數(shù)，一般；為s平面左半平面的極點。</p>&

89、lt;p>　　低通巴特沃斯濾波器設計步驟如下：</p><p><b>　　確定階數(shù)N。</b></p><p><b>　　取N=4</b></p><p><b>　　求極點</b></p><p><b>　　，，，，</b></p>

90、;<p>　　歸一化低通原型系統(tǒng)函數(shù)為</p><p>　　由N=4直接查表得到：</p><p><b>　　極點： </b></p><p>　　歸一化低通濾波器系統(tǒng)函數(shù)為</p><p>　　式中，0.0000,0.0999,0.1914,0.0252</p><p>&

91、lt;b>　　將去歸一化最終得到</b></p><p>　　通過計算可以總結出過程太麻煩，而且容易出錯，結果不直觀。</p><p>　　下面用M程序設計來實現(xiàn)。</p><p>　　用脈沖響應不變法設計的巴特沃斯數(shù)字低通濾波器的M程序如下：</p><p><b>　　fp=2100;</b><

92、;/p><p><b>　　fs=8000;</b></p><p><b>　　Fs=20000;</b></p><p><b>　　Rp=0.5;</b></p><p><b>　　Rs=30;</b></p><p>　　T=

93、1/Fs; %設計指標</p><p>　　W1p=fp/Fs*2;W1s=fs/Fs*2;%求歸一化頻率</p><p>　　[N,Wn]=buttord(W1p,W1s,Rp,Rs,'s');</p><p>　　%確定butterworth的最小介數(shù)N和頻率參數(shù)Wn</p><p>　　[z

94、,p,k]=buttap(N); %設計模擬低通原型的零極點增益參數(shù)</p><p>　　[bp,ap]=zp2tf(z,p,k); %將零極點增益轉換成分子分母參數(shù)</p><p>　　[bs,as]=lp2lp(bp,ap,Wn*pi*Fs);%將低通原型轉換為模擬低通</p><p>　　[bz,az]=impinvar(bs,as,Fs);

95、 %用脈沖響應不變法進行模數(shù)變換</p><p>　　sys=tf(bz,az,T); %給出傳輸函數(shù)H(Z)</p><p>　　[H,W]=freqz(bz,az,512,Fs); %生成頻率響應參數(shù)</p><p>　　subplot(2,1,1);</p><p>　　plot(W,20*log10(ab

96、s(H))); %繪制幅頻響應</p><p>　　grid on; %加坐標網(wǎng)格</p><p>　　xlabel('頻率/Hz'); </p><p>　　ylabel('振幅/dB');</p><p>　　subplot(2,1,2); </p&g

97、t;<p>　　plot(W,abs(H)); grid on; </p><p>　　xlabel('頻率/Hz'); </p><p>　　ylabel('振幅/H');</p><p><b>　　運行后的波形如下：</b></p><p>　　圖4.2 典型濾波器在M

98、atlab上運行波形</p><p><b>　　運行結果：</b></p><p><b>　　N =4</b></p><p>　　bz = 0.0000 0.0999 0.1914 0.0252</p><p>　　az= 1.0000 -1.4336 1.09

99、84 -0.4115 0.0627</p><p>　　可以看出：通過編程，結果非常直觀，題目中要求的通帶邊界頻fp=2.1kHZ，通帶最大衰減Rp=0.5dB；阻帶邊界頻率fs=8kHZ，阻帶最小衰減Rs=30dB,采樣頻率為Fs=20kHZ，等設計指標都達到，濾波器設計符合要求。</p><p>　　B 橢圓低通濾波器</p><p>　　橢圓低通濾

100、波器采用有限零點設計來逼近理想的低通濾波器的特性。因為零點的位置與橢圓函數(shù)的許多特性有關，所以稱之為橢圓低通濾波器。它的幅度平方函數(shù)為：</p><p><b>　　(4.3)</b></p><p>　　式中的是雅可比橢圓函數(shù)，是與通帶衰減有關的參數(shù)，階數(shù)N等于通帶和阻帶內最大點和最小點的總和。階數(shù)和系統(tǒng)函數(shù)由通帶的截止頻率，通帶內最大衰減，阻帶截止頻率及阻帶內的最

101、小衰減等參數(shù)確定。</p><p>　　雙線性變換法設計步驟如下：</p><p>　?。?）首先寫出該濾波器的系統(tǒng)函數(shù)</p><p>　?。?）利用雙線性變換法轉換，數(shù)字濾波器的系統(tǒng)函數(shù)為</p><p>　　用雙線性變換法設計橢圓數(shù)字低通濾波器的M程序如下：</p><p><b>　　fs=20000

102、;</b></p><p>　　wp=2*pi*2100/fs; </p><p>　　ws=2*pi*8000/fs; </p><p><b>　　Rp=0.5; </b></p><p><b>　　Rs=30; </b></p><p><b>

103、　　Ts=1/fs; </b></p><p>　　Wp=2/Ts*tan(wp/2);Ws=2/Ts*tan(ws/2); %按頻率轉換公式進行轉換 </p><p>　　[N,Wn]=ellipord(Wp,Ws,Rp,Rs,'s'); %計算模擬濾波器的最小階數(shù) </p><p>　　[z,p,k]=ellipap(N,Rp,Rs)

104、;%設計模擬原型濾波器 </p><p>　　[Bap,Aap]=zp2tf(z,p,k); %零點極點增益形式轉換為傳遞函數(shù)形式 </p><p>　　[b,a]=lp2lp(Bap,Aap,Wn); %低通轉換為低通濾波器的頻率轉化 </p><p>　　[bz,az]=bilinear(b,a,fs); %運用雙線性變換法得到數(shù)字濾波器傳遞函數(shù) </p&

105、gt;<p>　　[H,f]=freqz(bz,az,512,fs); </p><p>　　subplot(2,1,1);</p><p>　　plot(f,20*log10(abs(H))); </p><p>　　title('N=2 頻率響應'); </p><p><b>　　grid on;

106、 </b></p><p>　　xlabel('頻率/Hz'); </p><p>　　ylabel('振幅/dB'); </p><p>　　subplot(2,1,2); </p><p>　　plot(f,abs(H)); grid on; </p><p>　　x

107、label('頻率/Hz'); </p><p>　　ylabel('振幅/H');</p><p>　　圖4.3 雙線性變換法設計的橢圓數(shù)字低通濾波器運行波形</p><p><b>　　運行結果：</b></p><p><b>　　N=2</b></p&g

108、t;<p>　　bz= 0.1213 0.1662 0.1213</p><p>　　az= 1.0000 -0.9889 0.4218</p><p>　　可以看出，濾波器的參數(shù)均符合設計要求，設計是成功的。同時了解到，橢圓數(shù)字濾波器的降斜度比較大，通帶和阻帶內均為等波紋，同樣的性能指標，橢圓濾波器可以更低的階數(shù)來實現(xiàn)。</p><

109、p>　　4.2.2切比雪夫（Chebyshev）數(shù)字濾波器的設計</p><p>　　切比雪夫Ⅰ型濾波器在通帶內幅度特性是等波紋的，在阻帶內是單調的。切比雪夫Ⅱ濾波器在通帶內是單調的，在阻帶內是等波紋的。切比雪夫Ⅰ型低通濾波器的幅度平方函數(shù)為：</p><p><b>　　(4.4)</b></p><p>　　上式結果是一個正實數(shù)，式中

110、的是一個小于1的正數(shù)，它與通帶的波紋有關，越大，波紋越大，為通帶的截止頻率，時N階切比雪夫多項式，定義為：</p><p><b>　　(4.5)</b></p><p>　　切比雪夫濾波器的傳遞函數(shù)參數(shù)確定：，，N。當N大于或等于1時，切比雪夫多項式遞推公式為：</p><p>　　切比雪夫濾波器的階數(shù)N等于通帶內等幅波動的次數(shù)，即最大值和最

111、小值的總數(shù)，可由下式確定</p><p>　　， (4.6)</p><p>　　其中均為與衰減有關的參數(shù)，Ap為通帶內的最大衰減，As為阻帶內的容許衰減，它們的單位為。</p><p>　　因為濾波器的幅度平方函數(shù)為：</p><p><b>　　(4.7)</b&g

112、t;</p><p>　　進而化簡得到：，又因為，由式(4-7)得到：</p><p>　　故濾波器的階數(shù)N為： (4.8)</p><p>　　若要求阻帶截止頻率上的衰減越大，或過渡帶內幅度特性越陡，則所需的階數(shù)N越高。從式中可知切比雪夫濾波器的幅度平方函數(shù)也只有極點沒有零點，且只需求出左半平面的極點，所以一旦求出，，

113、N，就能確定。</p><p>　　利用MATLAB設計切比雪夫濾波器的實例。</p><p>　　A、現(xiàn)在以MATLAB的函數(shù)cheblap(語法：［z,p,k］=cheblap(n，rp)，其中n為階數(shù)，rp為通帶的幅度差，z,p,k分別代表零點，極點，增益)為工具設計切比雪夫Ⅰ濾波器。設計技術指標為：通帶截止頻率，通帶的最大衰減，阻帶截止頻率，阻帶的最小衰減。</p>

114、<p><b>　　程序如下：</b></p><p>　　a=10000;b=40000;AP=3;AS=35; %技術指標</p><p>　　T1=sqrt(10^(0.1*AP)-1);</p><p>　　T2=sqrt(10^(0.1*AS)-1);</p><p>　　N=ceil(a

115、cosh(T2/T1)/acosh(b/a)); %求階數(shù)</p><p>　　[Z,P,K]=cheb1ap(N,AP); %求零、極點，增益</p><p>　　syms rad; %定義基本符號對象</p><p>　　hs1=K/(i*rad-P(1))/(i*rad-P(2)