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

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

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

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

4、bp、lp2bs等函數(shù)實(shí)現(xiàn)了從模擬低通到高通的轉(zhuǎn)化，并對圖形做了對比分析關(guān)鍵詞：MATLAB；濾波器；巴特沃斯；系統(tǒng)函數(shù)目 錄1 課題描述12模擬濾波器的設(shè)計(jì)22.1模擬低通濾波器的設(shè)計(jì)指標(biāo)及逼近方法2 2.2 巴特沃斯低通濾波器的設(shè)計(jì)3 2.3 模擬濾波器的轉(zhuǎn)化原理6 2.4低通到高通的頻率轉(zhuǎn)換73設(shè)計(jì)過程73.1軟件介紹73.2設(shè)計(jì)內(nèi)容83.3設(shè)計(jì)步驟84程序運(yùn)行結(jié)果及分析10總 結(jié)11參考文獻(xiàn)121 課題描述濾波器（filter），是一種用來消除干擾雜訊的器件，將輸入或輸出經(jīng)過過濾而得到純凈的直流電。對特定頻率的頻點(diǎn)或該頻點(diǎn)以外的頻率進(jìn)行有效濾除的電路，就是濾波器，其功能就是得到某些特定

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

6、通帶或阻帶有波動，可以提高選擇性；貝塞爾濾波器通帶內(nèi)有較好的線性相位特性；橢圓濾波器的選擇性相對前三種是最好的。本課題正是在設(shè)計(jì)好低通濾波器的基礎(chǔ)上在逐步轉(zhuǎn)變?yōu)楦咄V波器的，由于在設(shè)計(jì)的過程中運(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)大功能，通過設(shè)置參通過設(shè)置參數(shù)和調(diào)用工具箱函數(shù)butter,cheby,cheby2等可以直接計(jì)算出模擬濾波器，并畫出其幅頻

7、特性，實(shí)現(xiàn)了從模擬低通濾波器到模擬高通的幅度特性對比演示，清晰地比較了這兩種濾波器的功能及用途。2模擬濾波器的設(shè)計(jì)2.1模擬低通濾波器的設(shè)計(jì)指標(biāo)及逼近方法在本課題中，分別用、表示模擬濾波器的單位沖擊響應(yīng)、系統(tǒng)函數(shù)、頻率響應(yīng)函數(shù)，三者的關(guān)系如下： （2.1.1） （2.1.2） 可以用、中任意一個(gè)描述模擬濾波器，也可以在線性廠系數(shù)微分方程描述模擬濾波器。但是設(shè)計(jì)模擬濾波器時(shí)，設(shè)計(jì)指標(biāo)一般有幅頻相應(yīng)函數(shù)給出，而模擬濾波器設(shè)計(jì)就是根據(jù)設(shè)計(jì)指標(biāo)，求系統(tǒng)函數(shù)。 工程實(shí)際中通常用所謂的損耗函數(shù)（也稱為衰減函數(shù)）來描述濾波器的幅頻響應(yīng)特性，對歸一化幅頻響應(yīng)函數(shù)，定義如下（其單位是分貝，用表示）： （2.1.

8、3）模擬低通濾波器的設(shè)計(jì)指標(biāo)參數(shù)有、和。其中和分別表示通帶邊界頻率和阻帶截止頻率，稱為通帶最大衰減，稱為最帶最小衰減，和的單位為。 模擬濾波器的技術(shù)指標(biāo)給定后，需要設(shè)計(jì)一個(gè)系統(tǒng)函數(shù)，希望其幅度平方函數(shù)滿足給定的指標(biāo)。一般濾波器的單位脈沖響應(yīng)為實(shí)函數(shù)，因此 （2.1.4）如果能由、和求出，那么就可以求出，由此可求出所需要的。必須是因果平穩(wěn)的，因此極點(diǎn)必須落在平面的左半平面，相應(yīng)的的極點(diǎn)必然落在右半平面。這就是由求所需求的的具體原則，即模擬低通濾波器的逼近方法。因此幅度平方函數(shù)在模擬濾波器的設(shè)計(jì)中起著很重要的作用。2.2 巴特沃斯低通濾波器的設(shè)計(jì) 巴特沃斯低通濾波器的幅度平方函數(shù)用下式表示： (2

9、.2.1)式中，N為濾波器的階數(shù)。當(dāng)，；時(shí)，是3截止頻率。以替換，將幅度平方函數(shù)寫成的函數(shù)： (2.2.2)復(fù)變量，此式表明幅度平方函數(shù)有2N個(gè)極點(diǎn)，極點(diǎn)用下式表示： (2.2.3)式中，。個(gè)極點(diǎn)等間隔分布在半徑為的圓上（該圓成為巴特沃斯園），間隔是。 為形成因果穩(wěn)定的濾波器，個(gè)極點(diǎn)中只取平面左半平面的個(gè)極點(diǎn)構(gòu)成，而右半平面的個(gè)極點(diǎn)構(gòu)成。的表達(dá)式為： (2.2.4)設(shè)N=3，極點(diǎn)有6個(gè)，它們分別為： (2.2.5) 取平面左半平面的極點(diǎn)組成系統(tǒng)函數(shù)，即： (2.2.6) 由于不同的技術(shù)指標(biāo)對應(yīng)的邊界頻率和濾波器幅頻特性不同，為使設(shè)計(jì)公式和圖表統(tǒng)一，將頻率歸一化。巴特沃斯濾波器采用對截止頻率歸一

10、化，歸一化以后的系統(tǒng)函數(shù)為： (2.2.7)式中，令，稱為歸一化頻率；令，p稱為歸一化復(fù)變量，這樣巴特沃斯濾波器的歸一化低通原型系統(tǒng)函數(shù)為： (2.2.8)式中，為歸一化極點(diǎn)，用下式表示： (2.2.9) 將極點(diǎn)展開可得到的的分母的階多項(xiàng)式，用下式表示： (2.2.10) 歸一化原型系統(tǒng)函數(shù)的系數(shù)，以及極點(diǎn)，可由下表得到:巴特沃斯歸一化低通濾波器分母多項(xiàng)式的因式分解 所以只要求出巴特沃斯濾波器的階數(shù)N和截止頻率，就可以求出濾波器的系統(tǒng)函數(shù)。所以巴特沃斯濾波器的設(shè)計(jì)實(shí)質(zhì)上就是根據(jù)設(shè)計(jì)指標(biāo)求階數(shù)和截止頻率的過程。下面介紹確定的方法。確定技術(shù)指標(biāo)：、和。將代入幅度平方函數(shù)中得： (2.2.11) 同

11、理可得： (2.2.12) 由以上兩式可得： (2.2.13) 令: ， (2.2.14) 則由下式表示： (2.2.15) 用上式求出的可能有小數(shù)部分，應(yīng)該取大于或等于的最小整數(shù)。關(guān)于截止頻率，如果技術(shù)指標(biāo)中沒有給出，可以有式求出得: (2.2.16) (2.2.17) 根據(jù)階數(shù)直接查上表可得，再將去歸一化，將代入得到實(shí)際的濾波器系統(tǒng)函數(shù): (2.2.18) 這里的為截止頻率。2.3 模擬濾波器的轉(zhuǎn)化原理 首先介紹各種頻率變化公式。從域到域映射的可逆變換記為。低通系統(tǒng)函數(shù)與之間的轉(zhuǎn)換關(guān)系為 (2.3.1) (2.3.2) 2.4 低通到高通的頻率轉(zhuǎn)換從低通到高通濾波器的映射關(guān)系為： (2.3

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

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

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

15、s=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=

16、2;rp=2,as=30; %設(shè)置濾波器的參數(shù)n,wc=buttord(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 onxlabel('頻率(Hz)');ylabel('幅度')綜合繪圖程序：wp=1,ws=2;rp=2,as=30; %設(shè)置濾

17、波器的參數(shù)n,wc=buttord(wp,ws,rp,as,'s'); %計(jì)算濾波器階數(shù)n和3dB截止頻率wcb,a=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/1000plot(n11,abs(hk(1:n1),'y-');hol

18、d onn12=(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 onbh,ah=lp2hp(b,a,2); %低通到

19、高通的轉(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=(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),'

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