基于matlab的iir數(shù)字濾波器的設(shè)計(jì)與仿真畢業(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章 濾波器的原理與設(shè)計(jì)7</p><p>　　2.1 模擬濾波器的原理與設(shè)計(jì)7</p><p>　　2.2 數(shù)字濾波器的原理與設(shè)計(jì)8</p><p>

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

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

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

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

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

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

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

10、的基本概念，其中包括系統(tǒng)的描述、系統(tǒng)的傳遞函數(shù)、系統(tǒng)的模型；接著簡單介紹MATLAB，并對數(shù)字濾波器在MATLAB環(huán)境下如何實(shí)現(xiàn)進(jìn)行了介紹；重點(diǎn)描述了IIR數(shù)字濾波器的設(shè)計(jì)過程，最后對IIR濾波器進(jìn)行仿真。</p><p>　　關(guān)鍵詞： 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、是信號處理的基礎(chǔ)，濾波運(yùn)算是信號處理中的基本運(yùn)算，濾波器的設(shè)計(jì)也就相應(yīng)成為數(shù)字信號處理的最基本問題之一。濾波器可廣義的理解為一個信號選擇系統(tǒng)，它讓某些信號成分通過又阻止或衰減另一些成分。在更多的情況下，濾波器可理解為選頻系統(tǒng)，如低通、高通、帶通、帶阻。濾波器可分為三種：模擬濾波器、采樣濾波器和數(shù)字濾波器。模擬濾波器可以是由RLC構(gòu)成的無源濾波器，也可以是加上運(yùn)放的有源濾波器，是連續(xù)時(shí)間系統(tǒng)；采樣濾波器由電阻、電容、電荷轉(zhuǎn)移器件、運(yùn)放等組

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

32、例或者濾除某些頻率成分的器件。數(shù)字濾波器和模擬濾波器相比，因?yàn)樾盘柕男问胶蛯?shí)現(xiàn)濾波的方法不同，數(shù)字濾波器具有比模擬濾波器精度高、穩(wěn)定、體積小、重量輕、靈活、不要求阻抗匹配等優(yōu)點(diǎn)。</p><p>　　數(shù)字濾波器從功能上分類:可分為低通濾波器、高通濾波器、帶通濾波器、帶阻濾波器。從濾波器的網(wǎng)絡(luò)結(jié)構(gòu)或者從單位脈沖響應(yīng)分類:可分為IIR濾波器(即無限長單位沖激響應(yīng)濾波器)和FIR濾波器(即有限長單位沖激響應(yīng)濾波器)。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è)計(jì)數(shù)字濾波器的思路與模擬濾波器相仿,設(shè)計(jì)實(shí)質(zhì)都是尋找一組系數(shù){b,a},去逼近所求的頻率響應(yīng),以便在性能上滿足預(yù)定的技術(shù)要求；不同的是數(shù)字濾波器是在Z平面尋找合適的H(Z), 模擬濾波器則是在S平面用數(shù)學(xué)逼近法尋找近似的所需特性H(S)。由于模擬濾波器的設(shè)計(jì)在理論上已十分成熟,因此數(shù)字濾波器設(shè)計(jì)的關(guān)鍵是將H(S)變換至H(Z),即利用復(fù)值映射將模擬濾波器離散

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

73、;p>　?。?）運(yùn)算符豐富。由于MATLAB是用C語言編寫的，MATLAB提供了和C語言幾乎一樣多的運(yùn)算符，靈活使用MATLAB的運(yùn)算符將使程序變得極為簡短。</p><p>　　（3）MATLAB既具有結(jié)構(gòu)化的控制語句（如for循環(huán)，while循環(huán)，break語句和if語句），又有面向?qū)ο缶幊痰奶匦浴?lt;/p><p>　?。?）程序限制不嚴(yán)格，程序設(shè)計(jì)自由度大。例如，在MATLAB

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

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

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

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

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

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

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

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

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

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

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

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

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

87、計(jì)巴特沃斯數(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值越大，通帶和阻帶的近似特性就越好。在截止頻率處，幅度平方響應(yīng)為=0處的1/2，相當(dāng)于幅度響應(yīng)的3dB衰減點(diǎn)。</p><p><b>　　其系統(tǒng)函數(shù)為 </b></p><p><b>　?。?.2）</b></p><p>　　式中，為歸一化常數(shù)，一般；為s平面左半平面的極點(diǎn)。</p>&

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

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

91、lt;b>　　將去歸一化最終得到</b></p><p>　　通過計(jì)算可以總結(jié)出過程太麻煩，而且容易出錯，結(jié)果不直觀。</p><p>　　下面用M程序設(shè)計(jì)來實(shí)現(xiàn)。</p><p>　　用脈沖響應(yīng)不變法設(shè)計(jì)的巴特沃斯數(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; %設(shè)計(jì)指標(biāo)</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è)計(jì)模擬低通原型的零極點(diǎn)增益參數(shù)</p><p>　　[bp,ap]=zp2tf(z,p,k); %將零極點(diǎn)增益轉(zhuǎn)換成分子分母參數(shù)</p><p>　　[bs,as]=lp2lp(bp,ap,Wn*pi*Fs);%將低通原型轉(zhuǎn)換為模擬低通</p><p>　　[bz,az]=impinvar(bs,as,Fs);

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

96、s(H))); %繪制幅頻響應(yīng)</p><p>　　grid on; %加坐標(biāo)網(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>　　運(yùn)行后的波形如下：</b></p><p>　　圖4.2 典型濾波器在M

98、atlab上運(yùn)行波形</p><p><b>　　運(yùn)行結(jié)果：</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>　　可以看出：通過編程，結(jié)果非常直觀，題目中要求的通帶邊界頻fp=2.1kHZ，通帶最大衰減Rp=0.5dB；阻帶邊界頻率fs=8kHZ，阻帶最小衰減Rs=30dB,采樣頻率為Fs=20kHZ，等設(shè)計(jì)指標(biāo)都達(dá)到，濾波器設(shè)計(jì)符合要求。</p><p>　　B 橢圓低通濾波器</p><p>　　橢圓低通濾

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

101、小衰減等參數(shù)確定。</p><p>　　雙線性變換法設(shè)計(jì)步驟如下：</p><p>　　（1）首先寫出該濾波器的系統(tǒng)函數(shù)</p><p>　?。?）利用雙線性變換法轉(zhuǎn)換，數(shù)字濾波器的系統(tǒng)函數(shù)為</p><p>　　用雙線性變換法設(shè)計(jì)橢圓數(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); %按頻率轉(zhuǎn)換公式進(jìn)行轉(zhuǎn)換 </p><p>　　[N,Wn]=ellipord(Wp,Ws,Rp,Rs,'s'); %計(jì)算模擬濾波器的最小階數(shù) </p><p>　　[z,p,k]=ellipap(N,Rp,Rs)

104、;%設(shè)計(jì)模擬原型濾波器 </p><p>　　[Bap,Aap]=zp2tf(z,p,k); %零點(diǎn)極點(diǎn)增益形式轉(zhuǎn)換為傳遞函數(shù)形式 </p><p>　　[b,a]=lp2lp(Bap,Aap,Wn); %低通轉(zhuǎn)換為低通濾波器的頻率轉(zhuǎn)化 </p><p>　　[bz,az]=bilinear(b,a,fs); %運(yùn)用雙線性變換法得到數(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 頻率響應(yīng)'); </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è)計(jì)的橢圓數(shù)字低通濾波器運(yùn)行波形</p><p><b>　　運(yùn)行結(jié)果：</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è)計(jì)要求，設(shè)計(jì)是成功的。同時(shí)了解到，橢圓數(shù)字濾波器的降斜度比較大，通帶和阻帶內(nèi)均為等波紋，同樣的性能指標(biāo)，橢圓濾波器可以更低的階數(shù)來實(shí)現(xiàn)。</p><

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

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

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

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

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

114、<p><b>　　程序如下：</b></p><p>　　a=10000;b=40000;AP=3;AS=35; %技術(shù)指標(biāo)</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); %求零、極點(diǎn)，增益</p><p>　　syms rad; %定義基本符號對象</p><p>　　hs1=K/(i*rad-P(1))/(i*rad-P(2)