模擬濾波器設(shè)計(jì)演示的軟件實(shí)現(xiàn)

1、數(shù)字信號(hào)課程設(shè)計(jì)課程名稱數(shù)字信號(hào)處理課程設(shè)計(jì)實(shí)驗(yàn)項(xiàng)目 模擬濾波器設(shè)計(jì)演示的軟件實(shí)現(xiàn)實(shí)驗(yàn)儀器計(jì)算機(jī)學(xué)院/系別通信工程班級(jí)/學(xué)號(hào)學(xué)生姓名實(shí)驗(yàn)日期成 績(jī)指導(dǎo)教師題目四模擬濾波器設(shè)計(jì)演示的軟件實(shí)現(xiàn)一.設(shè)計(jì)目的熟悉和鞏固模擬濾波器的設(shè)計(jì)方法和原理；熟練掌握MATLAB：具軟件在工程設(shè)計(jì)中的使用； 熟練掌握模擬濾波器的設(shè)計(jì)及由 ALPF到各型模擬濾波器的幅頻特性對(duì)比。二.設(shè)計(jì)內(nèi)容 動(dòng)態(tài)演示由ALPF到模擬各型濾波器幅度特性曲線；分析并說(shuō)明由沖擊響應(yīng)不變法設(shè)計(jì) BSF和HPF加保護(hù)濾波器的必要性。三、設(shè)計(jì)要求 理論分析 覬ALPF勺設(shè)計(jì)及由ALP囹J模擬各型濾波器的設(shè)計(jì)過(guò)程； 動(dòng)態(tài)演示由ALP囹J模擬各型濾

2、波器幅度特性曲線；通帶、過(guò)渡帶和阻帶分別用不同的顏色；分析并說(shuō)明由沖擊響應(yīng)不變法設(shè)計(jì)BSFF口HPFm保護(hù)濾波器的必要性。四、實(shí)驗(yàn)儀器計(jì)算機(jī)1臺(tái)，安裝MATLAB件五、實(shí)驗(yàn)步驟 設(shè)計(jì)過(guò)程詳見教材相關(guān)內(nèi)容； 使用巴特沃斯濾波器，其階數(shù) N應(yīng)該根據(jù)實(shí)際參數(shù)計(jì)算（計(jì)算公式和方法如 教材所述），為方便作圖，這里指定階數(shù)為N=5,并假定通帶截止頻率p=1rad, 阻帶截止頻s=2 rad ,;分別用不同顏色曲線繪制通帶、過(guò)渡帶和阻帶。要求根據(jù)變換關(guān)系動(dòng)態(tài)演示低通濾波器和目標(biāo)濾波器的幅度特性。簡(jiǎn)要說(shuō)明采用沖擊響應(yīng)不變法對(duì) AHP標(biāo)口 ABS嚶字化時(shí)保護(hù)濾波器的作用。MATLAB是矩陣實(shí)驗(yàn)室”(MATrix

3、 LABoratoy )的縮寫，是一種科學(xué)計(jì)算 軟件，主要適用于矩陣運(yùn)算及控制和信息處理領(lǐng)域的分析設(shè)計(jì)，它使用方便，輸入簡(jiǎn)捷，運(yùn)算高效，內(nèi)容豐富。本課題在設(shè)計(jì)模擬濾波器的過(guò)程中，首先對(duì)設(shè)計(jì)的原理及方法做了非常詳細(xì) 的分析，特別是在設(shè)計(jì)高通濾波器的過(guò)程中，利用現(xiàn)有低通濾波器的系統(tǒng)函數(shù)， 在經(jīng)過(guò)簡(jiǎn)單的轉(zhuǎn)化之后，可以直接設(shè)計(jì)出這種濾波器，但在實(shí)際手工計(jì)算中特別 麻煩，所以本實(shí)驗(yàn)全是利用 MATLAB的中設(shè)計(jì)濾波器專用工具箱，在對(duì)低通的 技術(shù)指標(biāo)確定之后，通過(guò)調(diào)用buttord、butter等函數(shù)，實(shí)現(xiàn)了對(duì)巴特沃斯模擬低 通濾波器的設(shè)計(jì)，通過(guò)調(diào)用lp2hp、lp2bp、lp2bs等函數(shù)實(shí)現(xiàn)了從模擬低通

4、到高 通的轉(zhuǎn)化，并對(duì)圖形做了對(duì)比分析關(guān)鍵詞：MATLAB ;濾波器；巴特沃斯；系統(tǒng)函數(shù)1課題描述5.2模擬濾波器的設(shè)計(jì)6.2.1 模擬低通濾波器的設(shè)計(jì)指標(biāo)及逼近方法 62.2 巴特沃斯低通濾波器的設(shè)計(jì)7.2.3 模擬濾波器的轉(zhuǎn)化原理 102.4 低通到高通的頻率轉(zhuǎn)換1.13設(shè)計(jì)過(guò)程1.13.1 軟件介紹1.13.2 設(shè)計(jì)內(nèi)容123.3 設(shè)計(jì)步驟124程序運(yùn)行結(jié)果及分析1.4總 結(jié)1.5參考文獻(xiàn)1.61課題描述濾波器（filter）,是一種用來(lái)消除干擾雜訊的器件，將輸入或輸出經(jīng)過(guò)過(guò)濾 而得到純凈的直流電。對(duì)特定頻率的頻點(diǎn)或該頻點(diǎn)以外的頻率進(jìn)行有效濾除的電 路，就是濾波器，其功能就是得到某些特定頻

5、率或消除某些特定頻率。濾波器分為模擬濾波器和數(shù)字濾波器，盡管數(shù)字濾波器有很多的優(yōu)勢(shì)，但設(shè)計(jì)數(shù)字濾波器 濾波器的基礎(chǔ)很多都是應(yīng)用模擬濾波器而設(shè)計(jì)的。模擬濾波器按幅度特性可分為低通， 高通，帶通和帶阻濾波器，但我們?cè)O(shè)計(jì) 濾波器時(shí)，總是先設(shè)計(jì)低通濾波器，再通過(guò)頻率變換將低通濾波器轉(zhuǎn)波換成希望 的濾器。模擬濾波器的理論和設(shè)計(jì)方法已發(fā)展的相當(dāng)成熟。且有若干典型的模擬濾波器供我們選擇，如巴特沃斯（Butterworth）濾波器，切比雪夫（Chebyshev港波 器，橢圓（Ellipse）濾波器，貝塞爾（Bessel）g波器等，這幾種濾波器各有特點(diǎn)，巴 特沃斯濾波器具有單調(diào)下降的幅頻特性；切比雪夫?yàn)V波器的幅

6、頻特性在通帶或阻 帶有波動(dòng)，可以提高選擇性；貝塞爾濾波器通帶內(nèi)有較好的線性相位特性；橢圓濾波器的選擇性相對(duì)前三種是最好的。本課題正是在設(shè)計(jì)好低通濾波器的基礎(chǔ)上在逐步轉(zhuǎn)變?yōu)楦咄V波器的，由于在設(shè)計(jì)的過(guò)程中運(yùn)算相當(dāng)大，所以在設(shè)計(jì)中使用了MATLAB軟件，MATLAB是MATrix LABoratory的縮寫，早期主要用于現(xiàn)代控制中復(fù)雜的矩陣、向量的各種 運(yùn)算。MATLAB以矩陣作為基本編程單元，它提供了各種矩陣的運(yùn)算與操作， 并有較強(qiáng)的繪圖功能。本課題利用了 MATLAB在設(shè)計(jì)濾波器上的強(qiáng)大功能，通 過(guò)設(shè)置參通過(guò)設(shè)置參數(shù)和調(diào)用工具箱函數(shù)butter,cheby,cheby方可以直接計(jì)算出模擬濾波器

7、，并畫出其幅頻特性，實(shí)現(xiàn)了從模擬低通濾波器到模擬高通的幅度特 性對(duì)比演示，清晰地比較了這兩種濾波器的功能及用途。2模擬濾波器的設(shè)計(jì)2.1 模擬低通濾波器的設(shè)計(jì)指標(biāo)及逼近方法在本課題中，分別用兒、Ha(S)、Ha(j)表示模擬濾波器的單位沖擊響應(yīng)、 系統(tǒng)函數(shù)、頻率響應(yīng)函數(shù)，三者的關(guān)系如下:Ha(s) LTha(t)ha(t)estdt(2.1.1 )Ha(j ) FTha(t)ha(t)e j tdt(2.1.2)可以用ha、Ha(s)、Ha(j)中任意一個(gè)描述模擬濾波器，也可以在線性廠系 數(shù)微分方程描述模擬濾波器。但是設(shè)計(jì)模擬濾波器時(shí)，設(shè)計(jì)指標(biāo)一般有幅頻相應(yīng)函數(shù)Ha(j )給出，而模擬濾波器設(shè)

8、計(jì)就是根據(jù)設(shè)計(jì)指標(biāo)，求系統(tǒng)函數(shù)Ha(s) 0工程實(shí)際中通常用所謂的損耗函數(shù)(也稱為衰減函數(shù))A()來(lái)描述濾波器的幅頻響應(yīng)特性，對(duì)歸一化幅頻響應(yīng)函數(shù)，A()定義如下(其單位是分貝，用dB表小)：A( ) 20lgHa(j )10lgHa(j )2dB(2.1.3)模擬低通濾波器的設(shè)計(jì)指標(biāo)參數(shù)有p、 p、，和s。其中p和s分別表示通帶邊界頻率和阻帶截止頻率，p稱為通帶最大衰減，s稱為最帶最小衰減，p和s的單位為dB。模擬濾波器的技術(shù)指標(biāo)給定后，需要設(shè)計(jì)一個(gè)系統(tǒng)函數(shù)Ha ,希望其幅度平方函數(shù)滿足給定的指標(biāo)。一般濾波器的單位脈沖響應(yīng)為實(shí)函數(shù)，因此Ha(j )2Ha(s)Ha( s)|s j Ha(j

9、)Ha(j )(2.1.4)如果能由p、 p、 s和s求出Ha(j )2，那么就可以求出Ha(s)Ha(s平面的此可求出所需要的Ha(s)o Ha(s)必須是因果平穩(wěn)的，因此極點(diǎn)必須落在左半平面，相應(yīng)的Ha( S)的極點(diǎn)必然落在右半平面。這就是由Ha(S)Ha( S)求所 需求的Ha(s)的具體原則，即模擬低通濾波器的逼近方法。因此幅度平方函數(shù)在 模擬濾波器的設(shè)計(jì)中起著很重要的作用。2.2 巴特沃斯低通濾波器的設(shè)計(jì)巴特沃斯低通濾波器的幅度平方函數(shù)Ha(j )2用下式表小:Ha(j )2121N(2.2.1)式中，N為濾波器的階數(shù)。當(dāng),Ha(j )1.c時(shí),Ha(jc是3dB截止頻率。以S替換,

10、將幅度平方函數(shù)Ha(j )2寫成S的函數(shù):Ha(S)Ha( S)1(2.2.2)c復(fù)變量sj ，此式表明幅度平方函數(shù)有2N個(gè)極點(diǎn)，極點(diǎn)Sk用下式表示:Sk( 1)2N(j c).J 2k 1 cej *(2.2.3)式中，k 0，1,2，,2N 1。2N個(gè)極點(diǎn)等間隔分布在半徑為c的圓上(該圓成rad為巴特沃斯因)，間隔是N為形成因果穩(wěn)定的濾波器，2N個(gè)極點(diǎn)中只取S平面左半平面的N個(gè)極點(diǎn)構(gòu)成Ha(s),而右半平面的N個(gè)極點(diǎn)構(gòu)成Ha( S)。Ha(s)的表達(dá)式為:NHa(s)(2.2.4)S Skk 0設(shè)N=3,極點(diǎn)有6個(gè)，它們分別為:SoS3,2ej3 sce, * S1,1e 1sce7 s4

11、c,S2c,S5,2 ce %,1ce3(2.2.5)取s平面左半平面的極點(diǎn)S1、S2組成系統(tǒng)函數(shù)Ha，即:(2.2.13)Ha(S)(2.2.6)為使設(shè)計(jì)公式3c和圖表統(tǒng)一，將頻率歸一化巴特沃斯濾波器采用對(duì)3dB截止頻率c歸一化,歸一化以后的系統(tǒng)函數(shù)為:Ga(二)c1N 1SSk(2.2.7)c ,人稱為歸一化頻率；令一化復(fù)變量，這樣巴特沃斯濾波器的歸一化低通原型系統(tǒng)函數(shù)為:、1Ga( P)n 1P Pkk 0(2.2.8)隊(duì)為歸一化極點(diǎn)，用下式表示：1 2k 1j (2 TN)Pk e 2 2N k 0,1, ,N 1(2.2.9)將極點(diǎn)展開可得到的Ga(p)的分母p的N階多項(xiàng)式，用下式表

12、示:1Ga( P)-NNHN 2PbN 1 PbN 2 Pbo(2.2.10)歸一化原型系統(tǒng)函數(shù)Ga(P) Ha(P)的系數(shù)bk,以及極點(diǎn)Pk,可由下表得到：巴特沃斯歸一化低通濾波器分母多項(xiàng)式的因式分解埠因式階敬eta -螞修汨式山電國(guó)心）月年、12(+1. 4H2f+l3( + +P(+D4c#*+ Ol 7654A+1B47fi/+B56* 0. M 7b/+1)5+ L 4142計(jì)1 ”川+ 1.931+1J7（產(chǎn)口.4450/+1N/+L組7匕+1”尸十，的1十N口十】）8C十加值/十萬(wàn)+】，1.62+1 )(+1, 9fiS+l9副)即+力3+*口/+ *5MHW+L 融所以只要求出

13、巴特沃斯濾波器的階數(shù)N和3dB截止頻率c ,就可以求出濾波器的系統(tǒng)函數(shù)Ha（s） 0所以巴特沃斯濾波器的設(shè)計(jì)實(shí)質(zhì)上就是根據(jù)設(shè)計(jì)指標(biāo)求階數(shù)N和3dB截止頻率c的過(guò)程。下面介紹確定N的方法。確定技術(shù)指標(biāo)：p、s。c代入幅度平方函數(shù)中得:2NP1010(2.2.11)同理可得:2Ns1010(2.2.12)由以上兩式可得:1010 1，1010 1令:則N由下式表示:s sppksps1010 1P1010 1(2.2.14)(2.2.15)kspsp用上式求出的N可能有小數(shù)部分，應(yīng)該取大于或等于 N的最小整數(shù)。關(guān)于3dB截止頻率c,如果技術(shù)指標(biāo)中沒有給出，可以有式求出得p 100.1 P 1 前S

14、 100.1 S 1 而S根據(jù)階數(shù)N直接查上表可得Ha(p),再將Ha(p)去歸一化，將(2.2.16)(2.2.17)sc代入Ha(p)得到實(shí)際的濾波器系統(tǒng)函數(shù)Ha(s)小忤)| sp c(2.2.18)這里的c為3dB截止頻率。2.3模擬濾波器的轉(zhuǎn)化原理首先介紹各種頻率變化公式。從 p域到s域映射的可逆變換記為p F S低通系統(tǒng)函數(shù)Q(p)與Hd s之間的轉(zhuǎn)換關(guān)系為Hd s Q(p)|p FS(2.3.1)Q(p) Hd s |sF(2.3.2)2.4低通到高通的頻率轉(zhuǎn)換從低通到高通濾波器的映射關(guān)系為:(2.3.3)在虛軸上該映射關(guān)系簡(jiǎn)化為如下頻率變換公式:P Ph(2.3.4)其中，Ph

15、為希望設(shè)計(jì)的高通濾波器HHPS的通帶邊界頻率。所以只要將式 帶如式，就可以將通帶邊界頻率 P的低通濾波器的系統(tǒng)函數(shù)Q(P)轉(zhuǎn)換成通帶邊界頻率為Ph的高通濾波器系統(tǒng)函數(shù):Hhp S Ha(p)|(2.3.5)3設(shè)計(jì)過(guò)程3.1 軟件介紹3.1.1 MATLAB 特點(diǎn)MATLAB是一種科學(xué)計(jì)算軟件，適用于工程應(yīng)用個(gè)領(lǐng)域的分析設(shè)計(jì)與復(fù)雜 計(jì)算，它使用方便，輸入簡(jiǎn)捷，運(yùn)算高效且內(nèi)容豐富，很容易由用戶自行擴(kuò)展。 因此，當(dāng)前已成為美國(guó)和其他發(fā)達(dá)國(guó)家大學(xué)教學(xué)和科學(xué)研究中最常用且必不可少 的工具。MATLAB以矩陣作為基本編程單元，它提供了各種矩陣的運(yùn)算與操作，并 有較強(qiáng)的繪圖功能。MATLAB集科學(xué)計(jì)算、圖像

16、處理、聲音處理于一身，是一 個(gè)高度的集成系統(tǒng)，有良好的用戶界面，并有良好的幫助功能。MATLAB不僅流行于控制界，在機(jī)械工程、生物工程、語(yǔ)音處理、圖像處理、信號(hào)分析、計(jì)算 機(jī)技術(shù)等各行各業(yè)中都有極廣泛的應(yīng)用。MATLAB語(yǔ)言的特點(diǎn)1.編程效率高 2.用戶使用方便3.擴(kuò)充能力強(qiáng) 4.語(yǔ)句簡(jiǎn)單，內(nèi)涵豐富5.高效方便的矩陣 和數(shù)組運(yùn)算6.方便的繪圖功能數(shù)字信號(hào)處理的主要對(duì)象是數(shù)字信號(hào)，且是采用運(yùn)算的方法達(dá)到處理目的，因 此，其實(shí)現(xiàn)方法不同于模擬信號(hào)的實(shí)現(xiàn)方法，基本上可以分為兩種實(shí)現(xiàn)方法，即軟 件實(shí)現(xiàn)方法和硬件實(shí)現(xiàn)的方法。而硬件實(shí)現(xiàn)指的是選用合適的DSP芯片,配有適合芯片語(yǔ)言及任務(wù)要求的軟件，實(shí)現(xiàn)某種

17、信號(hào)處理功能的一種方法。數(shù)字信號(hào)處 理的特點(diǎn)1.靈活性2.高精度和高穩(wěn)定性3.便于大規(guī)模集成4.對(duì)數(shù)字信號(hào)可以存 儲(chǔ)、運(yùn)算；系統(tǒng)可以獲得高性能指標(biāo)3.2 設(shè)計(jì)內(nèi)容模擬濾波器設(shè)計(jì)演示一一從模擬低通濾波器到模擬高通的幅度特性對(duì)比演 示。使用巴特沃斯濾波器，其階數(shù)N應(yīng)該根據(jù)實(shí)際參數(shù)計(jì)算(計(jì)算公式和方法如 教材所述)，為方便作圖，這里指定階數(shù)為N=5,并假定通帶截止頻率Qp=1rad 阻帶截止頻率Q s=2rad,;3.3 設(shè)計(jì)步驟根據(jù)要求首先設(shè)計(jì)模擬低通濾波器，由于題中沒給通帶的最大衰減和阻帶的 最小衰減，所以在本題中先設(shè)出這兩個(gè)技術(shù)指標(biāo)，但最后計(jì)算N仍然為題目所要的5,進(jìn)過(guò)幾次試驗(yàn)，得出合適的值分

18、別為ap=2dB aS=30dB利用巴特沃斯設(shè)計(jì)的模擬濾波器的程序如下：低通濾波器的程序：wp=1,ws=2;rp=2, as=30;%設(shè)置濾波器的參數(shù)n,wc=buttord(wp,ws,rp,as,s) %計(jì)算濾波器階數(shù) n 和 2dB 截止頻率 wcb,a=butter(n,wc,s);%計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式系數(shù)f=0:1/100:0.5;wk=2*pi*f;hk=freqs(b,a,wk);plot(f,abs(hk);grid onxlabel(頻率(Hz);ylabel(幅度)低通到高通的轉(zhuǎn)化程序：wp=1,ws=2;rp=2,as=30;%設(shè)置濾波器的參數(shù)n,wc=b

19、uttord(wp,ws,rp,as,s)%計(jì)算濾波器階數(shù) n 和 3dB 截止頻率 wcb,a=butter(n,wc,s);%計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式系數(shù)bh,ah=lp2hp(b,a,2);%低通到高通的轉(zhuǎn)化 f=0:1/100:0.5;wk=2*pi*f; hk=freqs(bh,ah,wk); plot(f,abs(hk);grid on xlabel(頻率(Hz);ylabel(幅度) 綜合繪圖程序： wp=1,ws=2;rp=2,as=30;%設(shè)置濾波器的參數(shù)n,wc=buttord(wp,ws,rp,as,s); %計(jì)算濾波器階數(shù) n 和 3dB 截止頻率 wc b,a

20、=butter(n,wc,s);%計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式系數(shù)f=0:1/1000:1;wk=2*pi*f; hk=freqs(b,a,wk) n1=ceil(wp/(2*pi)/0.001); n2=ceil(ws-wp)/(2*pi)/0.001); n3=1/0.001-n1-n2; subplot(2,1,1) n11=1/1000:1/1000:n1/1000 plot(n11,abs(hk(1:n1),y-); hold on n12=(n1+1)/1000:1/1000:(n1+n2)/1000; plot(n12,abs(hk(n1+1):(n1+n2),g-); ho

21、ld on n13=(n1+n2+1)/1000:1/1000:(n1+n2+n3)/1000; plot(n13,abs(hk(n1+n2+1):(n1+n2+n3),b-);xlabel(頻率(Hz);ylabel(幅度);axis(0,0.6,0,1.2) title(低通濾波器);grid onbh,ah=lp2hp(b,a,2);%低通到高通的轉(zhuǎn)化f=0:1/1000:1;wk=2*pi*f;hk=freqs(bh,ah,wk);subplot(2,1,2)n11=1/1000:1/1000:n1/1000;plot(n11,abs(hk(1:n1),y-);hold onn12=(

22、n1+1)/1000:1/1000:(n1+n2)/1000;plot(n12,abs(hk(n1+1):(n1+n2),g-);hold onn13=(n1+n2+1)/1000:1/1000:(n1+n2+n3)/1000;plot(n13,abs(hk(n1+n2+1):(n1+n2+n3),b-);xlabel(頻率(Hz);ylabel(幅度);axis(0,0.6,0,1.2)title(高通濾波器);grid on4程序運(yùn)行結(jié)果及分析由以上程序可得低通、高通濾波器的幅度圖如下，運(yùn)行結(jié)果如下圖所示:斯1率什卬圖中我們可以清晰地看出各種濾波器的幅頻特性且用不同顏色表示出各種濾波 器的通帶和阻帶，在低通濾波器中，可以看出通帶的截止頻率的阻帶的截止頻率， 同時(shí)通過(guò)簡(jiǎn)單的轉(zhuǎn)化得出了高通濾波器的幅頻特性， 結(jié)果與要求一致，效果理想。 總 結(jié)就這門課程的設(shè)計(jì)需要另外的一門課程做輔助，現(xiàn)在的輔助工具是 MATLAB ,在設(shè)計(jì)的過(guò)程中讓我更好的學(xué)好，和對(duì)它也有了比較深的體會(huì)。通 過(guò)這次課程設(shè)計(jì)，我掌