巴特沃斯濾波器課程設(shè)計

1、<p>　　《數(shù)字信號處理》課程設(shè)計報告</p><p>　　學(xué) 院 信息工程學(xué)院</p><p>　　專 業(yè) 通信工程專業(yè)</p><p><b>　　班 級 </b></p><p><b>　　學(xué) 號 </b></p><p>&l

2、t;b>　　學(xué)生姓名 </b></p><p><b>　　二 0一0年十二月</b></p><p>　　巴特沃斯高通數(shù)字濾波器設(shè)計</p><p>　　要求：設(shè)計一個數(shù)字高通濾波器，要求通帶截止頻率p=0.8rad，通帶衰減不大于3dB，阻帶截止頻率s=0.44dB，阻帶衰減不小于15dB。希望采用巴特沃斯型濾波器。<

3、;/p><p><b>　　一、課程設(shè)計目的：</b></p><p>　　1.了解巴特沃斯高通數(shù)字濾波器濾波器的概念及原理。</p><p>　　2、掌握巴特沃斯高通數(shù)字濾波器濾波器的設(shè)計方法。</p><p>　　3、綜合運用數(shù)字信號處理的理論知識對語音信號進(jìn)行時頻分析和濾波器設(shè)計，通過理論推導(dǎo)得出相應(yīng)結(jié)論，再利用 MA

4、TLAB 作為編程工具進(jìn)行計算機(jī)實現(xiàn)。</p><p>　　4、綜合運用專業(yè)及基礎(chǔ)知識，解決實際工程技術(shù)問題的能力。</p><p><b>　　二、設(shè)計原理：</b></p><p><b>　　1、數(shù)字濾波器</b></p><p>　　所謂數(shù)字濾波器，是指輸入輸出均為數(shù)字信號，通過數(shù)字運算處理

5、改變輸入信號所含頻率成分的相對比例，或者濾除某些頻率成分的數(shù)字器件或程序。正因為數(shù)字濾波通過數(shù)值運算實現(xiàn)濾波，所以數(shù)字濾波器處理精度高、穩(wěn)定、體積小、重量輕、靈活、不存在阻抗匹配問題，可以實現(xiàn)模擬濾波器無法實現(xiàn)的特殊濾波功能。</p><p>　　按照不同的分類方式，數(shù)字濾波器可以有很多種類型，但總起來可以分為兩大類：經(jīng)典濾波器和現(xiàn)代濾波器。經(jīng)典濾波器的特點是其輸入信號中有用的頻率成分和希望濾除的成分分別占有不同

6、的頻帶，通過一個合適的選頻濾波器濾除干擾，得到純凈信號，達(dá)到濾波目的。但是，如果信號和干擾的頻譜相互重疊，則經(jīng)典濾波器無法有效濾除干擾，最大限度恢復(fù)信號，這就需要現(xiàn)代濾波器?，F(xiàn)代濾波器是根據(jù)隨機(jī)信號的一些統(tǒng)計特性，在某種最佳準(zhǔn)則下，最大限度抑制干擾，同時最大限度恢復(fù)信號，達(dá)到最佳的濾波效果的目的。</p><p><b>　　2、模擬濾波器</b></p><p>　

7、　模擬濾波器的理論和設(shè)計方法已經(jīng)發(fā)展的相當(dāng)成熟，且有多種典型的模擬濾波器供我們選擇，如巴特沃斯濾波器、切比雪夫濾波器、橢圓濾波器、貝塞爾濾波器等。這些濾波器有嚴(yán)格的設(shè)計公式、現(xiàn)成的曲線和圖表供設(shè)計人員使用，而且所涉及的系統(tǒng)函數(shù)都滿足電路實現(xiàn)條件。這些典型的濾波器各有特點：巴特沃斯濾波器具有單調(diào)下降的幅頻特性；切比雪夫濾波器的幅頻特性在通帶或阻帶有等波紋特性，可以提高選擇性；貝塞爾濾波器通帶內(nèi)有較好的線性相位特性；橢圓濾波器的特性相對前三

8、者是最好的，但通帶和阻帶內(nèi)均呈現(xiàn)等波紋特性，相位特性的非線性也稍嚴(yán)重。設(shè)計時，根據(jù)具體要求選擇濾波器的類型。</p><p>　　現(xiàn)在，我們分別使用、H、H表示模擬濾波器的單位沖激響應(yīng)、頻率響應(yīng)函數(shù)、系統(tǒng)函數(shù)，三者的關(guān)系為：</p><p>　　可以用、、中任意一個描述模擬濾波器，也可以用線性常系數(shù)微分方程描述模擬濾波器。但在設(shè)計濾波器時，設(shè)計指標(biāo)一般由幅頻響應(yīng)函數(shù)給出，而模擬濾波器設(shè)計就

9、是根據(jù)設(shè)計指標(biāo)，求系統(tǒng)函數(shù)。</p><p>　　工程上通常用所謂的損耗函數(shù)（也稱為衰減函數(shù)）來描述濾波器的幅頻響應(yīng)特性，對歸一化幅頻響應(yīng)函數(shù)，定義如下（其單位是分貝，用dB表示）：</p><p><b>　　dB</b></p><p>　　應(yīng)當(dāng)注意，損耗函數(shù)和幅頻特性函數(shù)只是濾波器幅頻響應(yīng)特性的兩種描述方法。損耗函數(shù)的優(yōu)點是對幅頻響應(yīng)的取

10、值非線性壓縮，放大了小的幅度，從而可以同時觀察通帶和阻帶的頻響特性的變化情況。</p><p>　　模擬低通濾波器的設(shè)計指標(biāo)參數(shù)有、、和。其中和分別稱為通帶邊界頻率和阻帶截止頻率，稱為通帶最大衰減（即通帶中允許的最大值），稱為阻帶最小衰減（即阻帶上允許的最小值），和的單位是dB。</p><p>　　對于單調(diào)下降的幅度特性，和可表示為：</p><p>　　濾波器的

11、技術(shù)指標(biāo)給定后，需要設(shè)計一個系統(tǒng)函數(shù)，希望其幅度平方函數(shù)滿足給定的指標(biāo)。一般的濾波器單位沖激響應(yīng)為實函數(shù)，因此</p><p>　　如果能由、、和求出，那么就可以求出，由此可求出所需要的。必須是因果穩(wěn)定的，因此極點必須落在s平面的左半平面，相應(yīng)的的極點必然落在右半平面。這就是由求所需要的的具體原則。</p><p><b>　　3、巴特沃斯濾波器</b></p&

12、gt;<p>　　巴特沃斯濾波器是電子濾波器的一種。巴特沃斯濾波器的特點是通頻帶的頻率響應(yīng)曲線最平滑。這種濾波器最先由英國工程師斯替芬·巴特沃斯（Stephen Butterworth）在1930年發(fā)表在英國《無線電工程》期刊的一篇論文中提出的。</p><p>　　巴特沃斯濾波器的特點是通頻帶內(nèi)的頻率響應(yīng)曲線最大限度平坦，沒有起伏，而在阻頻帶則逐漸下降為零。 在振幅的對數(shù)對角頻率的波特圖

13、上,從某一邊界角頻率開始,振幅隨著角頻率的增加而逐步減少,趨向負(fù)無窮大。 </p><p>　　一階巴特沃斯濾波器的衰減率為每倍頻6分貝，每十倍頻20分貝。二階巴特沃斯濾波器的衰減率為每倍頻12分貝、 三階巴特沃斯濾波器的衰減率為每倍頻18分貝、如此類推。巴特沃斯濾波器的振幅對角頻率單調(diào)下降，并且也是唯一的無論階數(shù)，振幅對角頻率曲線都保持同樣的形狀的濾波器。只不過濾波器階數(shù)越高，在阻頻帶振幅衰減速度越快。其他濾波

14、器高階的振幅對角頻率圖和低級數(shù)的振幅對角頻率有不同的形狀。</p><p>　　4、設(shè)計IIR數(shù)字濾波器的頻率變換法</p><p>　　（1）、歸一化模擬低通原型到數(shù)字高通濾波器的變換</p><p>　　設(shè)歸一化模擬低通原型濾波器的系統(tǒng)函數(shù)是，p為模擬域內(nèi)的拉普拉斯變量。在模擬域內(nèi)從低通到高通變換，即以代替p,有</p><p>　　進(jìn)行

15、反歸一化后，即以代入上式，得到反歸一化后的高通濾波器的傳輸函數(shù)</p><p><b>　　；</b></p><p>　　再進(jìn)雙線性變換，得數(shù)字高通濾波器的系統(tǒng)函數(shù)為</p><p><b>　??；</b></p><p>　　直接由歸一化低通原型變換為數(shù)字高通濾波器時，變換關(guān)系為</p&g

16、t;<p><b>　　，，</b></p><p>　　其中，為模擬濾波器3dB截止頻率；為數(shù)字濾波器3dB截止頻率。</p><p>　?。?）歸一化模擬低通原型到數(shù)字帶通濾波器的頻率變換</p><p>　　直接尋求從模擬低通到數(shù)字帶通之間的映射關(guān)系：</p><p><b>　　，<

17、/b></p><p><b>　　，</b></p><p>　　其中，為模擬濾波器3dB截止頻率；分別是數(shù)字帶通濾波器的中心頻率、通帶上邊頻、通帶下邊頻和帶通的上阻帶下邊頻。</p><p>　　（3）歸一化模擬低通原型到數(shù)字帶阻濾波器的變換</p><p>　　直接尋求從模擬低通到數(shù)字帶阻之間的映射關(guān)系：&l

18、t;/p><p>　　其中，為模擬濾波器3dB截止頻率；分別為數(shù)字帶阻濾波器的中心頻率、阻帶上邊頻、阻帶下邊頻和帶阻的上通帶下邊頻。</p><p>　?。?）數(shù)字低通濾波器到數(shù)字低通的變換</p><p>　　由截頻為的數(shù)字低通濾波器系統(tǒng)函數(shù)得到截頻為的數(shù)字低通濾波器系統(tǒng)函數(shù)，變換關(guān)系為</p><p>　?。?）數(shù)字低通濾波器到數(shù)字高通變換&

19、lt;/p><p>　　由截頻的數(shù)字低通濾波器系統(tǒng)函數(shù)得到截頻為的數(shù)字高通濾波器系統(tǒng)函數(shù)，變換關(guān)系為</p><p>　?。?）數(shù)字低通濾波器到數(shù)字帶通變換</p><p>　　由截頻為的數(shù)字低通濾波器系統(tǒng)函數(shù)得到中心頻率為，帶通上、下截頻分別為的數(shù)字帶通濾波器系統(tǒng)函數(shù)，變換關(guān)系為</p><p>　　5、MATLAB模擬頻率變換函數(shù)</p

20、><p>　　MATLAB信號處理工具箱函數(shù)lp2lp、lp2bp、lp2hp、lp2bs是模擬頻率變換函數(shù)。</p><p>　?。?）低通模擬濾波器到低通模擬濾波器變換函數(shù)lp2lp</p><p><b>　　調(diào)用格式說明如下。</b></p><p>　　[bs,as]=lp2lp(b,a,Wl)</p>

21、<p>　　該函數(shù)將截止頻率為1rad/s（歸一化截止頻率）的模擬低通原型濾波器變換成截止頻率為Wl的低通濾波器。</p><p>　?。?）低通模擬濾波器到帶通模擬濾波器的變換函數(shù)lp2bp</p><p><b>　　調(diào)用格式如下。</b></p><p>　　[bs,as]=lp2bp(b,a,W0,Bw)</p>

22、;<p>　　該函數(shù)將截止頻率為1rad/s（歸一化截止頻率）的模擬低通原型濾波器變換成具有指定帶寬和中心頻率為的帶通濾波器。如果要求濾波器的低端截止頻率為，高端截止頻率為，則按下式可計算帶寬和中心頻率為的帶通濾波器。</p><p><b>　　Bw=W2-W1</b></p><p>　　W0=sqrt(W1*W2)</p><p

23、>　?。?）低通模擬濾波器到高通模擬濾波器變換函數(shù)lp2hp</p><p><b>　　調(diào)用格式如下。</b></p><p>　　[bs,as]=lp2hp(b,a,Wh)</p><p>　　該函數(shù)將截止頻率為1rad/s（歸一化截止頻率）的模擬低通原型濾波器變換成截止頻率為的高通濾波器。</p><p>　

24、?。?）低通模擬濾波器到帶阻模擬濾波器變換函數(shù)lp2bs</p><p><b>　　調(diào)用格式如下。</b></p><p>　　[bs,as]=lp2bs(b,a,Wh)</p><p>　　該函數(shù)將截止頻率為1rad/s（歸一化截止頻率）的模擬低通原型濾波器變換成具有指定帶寬和中心頻率為的帶阻濾波器。</p><p>

25、;　　三、設(shè)計步驟與流程框圖：</p><p><b>　　四、設(shè)計過程：</b></p><p><b>　　1、設(shè)計要求：</b></p><p>　　設(shè)計一個數(shù)字高通濾波器，要求通帶截止頻率p=0.8rad，通帶衰減不大于3dB，阻帶截止頻率s=0.44dB，阻帶衰減不小于15dB。希望采用巴特沃斯型濾波器。<

26、/p><p><b>　　2、設(shè)計過程：</b></p><p>　　（1）確定數(shù)字高通的技術(shù)指標(biāo)：</p><p><b>　　， </b></p><p><b>　　rad, </b></p><p>　　將高通數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成高通模擬濾波器

27、的設(shè)計指標(biāo)：令T=2 s,預(yù)畸變校正得到模擬邊界頻率：</p><p>　　模擬低通濾波器的技術(shù)指標(biāo)計算如下：對通帶邊界頻率歸一化，即</p><p><b>　　==1，</b></p><p>　　將=1和-代入，求出歸一化低通濾波器的租帶截止頻率</p><p>　　=/=3.7190, </p>

28、<p>　　設(shè)計歸一化模擬濾波器G（p）</p><p>　　=5.5463， =/=3.719</p><p>　　N=lg/lg=1.3043, 取N=2</p><p>　　利用頻率變換公式將G（p）轉(zhuǎn)換成模擬高通</p><p>　　用雙線性變換法將模擬高通轉(zhuǎn)換成數(shù)字高通</p><p>　　利用MA

29、TLAB設(shè)計巴特沃斯高通濾波器</p><p><b>　　源程序：</b></p><p><b>　　頻率變換函數(shù)編程：</b></p><p>　　wpz=0.8;wsz=0.44;rp=3;rs=15;</p><p>　　[N,wc]=buttord(wpz,wsz,rp,rs);<

30、/p><p>　　[Bz,Az]=butter(N,wc,'high');</p><p>　　w=linspace(0,pi,1024);</p><p>　　h=freqz(Bz,Az,w);</p><p>　　plot(w/pi,20*log10(abs(h)))</p><p>　　axis([0

31、,1,-100,0]);grid;</p><p>　　xlabel('歸一化頻率');</p><p>　　ylabel('幅度/dB');</p><p>　　數(shù)字濾波器耗損函數(shù)如圖所示:</p><p><b>　　N=2</b></p><p>　　Bz=[

32、0.1326 -0.2653 0.1326]</p><p>　　Az=[1.0000 0.7394 0.2699]</p><p><b>　　五、結(jié)果分析：</b></p><p>　　由最終所得到的運行結(jié)果以及該濾波器的幅度響應(yīng)函數(shù)圖像可以看出，我們所設(shè)計的濾波器完全符合最初的要求，是一個巴特沃斯高通濾波器，通帶截止頻率p=0.

33、8rad，通帶衰減不大于3dB，阻帶截止頻率s=0.44dB，阻帶衰減不小于15dB。</p><p><b>　　六、總結(jié)：</b></p><p>　　確定好所需要設(shè)計的濾波器用途后，首先對其進(jìn)行性能需求分析，明確該系統(tǒng)應(yīng)該達(dá)到的各種性能指標(biāo)，其次，擬定多種濾波器類型，對這些方案采用Matlab進(jìn)行仿真嗎，進(jìn)行綜合分析和比較，選擇出最佳的濾波器類型作為本設(shè)計方案。

34、</p><p>　　系統(tǒng)思維 必須有系統(tǒng)的設(shè)計思維，把每一個細(xì)節(jié)都放到整個系統(tǒng)中考慮，考慮整個系統(tǒng)設(shè)計的可行性、穩(wěn)定性和功能的實現(xiàn)，這樣才不會局限在細(xì)節(jié)上，才能快速的完成性能優(yōu)越的硬件設(shè)計。</p><p>　　把握細(xì)節(jié) 系統(tǒng)也是由細(xì)節(jié)構(gòu)成的，在把握整個系統(tǒng)思維的基礎(chǔ)上把握每一個細(xì)節(jié)，因為每一個細(xì)節(jié)都有可能決定整個系統(tǒng)的性能。勇于嘗試，系統(tǒng)即使經(jīng)過了非常嚴(yán)格的論證也仍然有可能存在問

35、題，或許面對問題一時沒有很好的解決方案，但是有一個或許可行的大膽的想法，不要猶豫，試一試，勇于嘗試往往能夠找到更好的解決方法。</p><p>　　保持自信 無論遇到什么困難，我都相信自己一定能夠找到解決的方法，有的時候只要稍微再用一點力，結(jié)果就會不一樣。</p><p><b>　　七、參考資料：</b></p><p>　　【1】丁玉美，

36、高西全， 《數(shù)字信號處理》 西安電子科技大學(xué)出版社，2001</p><p>　　【2】胡廣數(shù)， 《數(shù)字信號處理實踐教程》 華中科技大學(xué)出版社，2002</p><p>　　【3】劉順蘭，吳杰， 數(shù)字信號處理教程，北京科學(xué)出版社，2003</p><p>　　【4】樓順天，李博菡 《基于MATLAB的系統(tǒng)分析與設(shè)計——信號處理》 西安電子科技大學(xué)出版社，2