帶通濾波器實(shí)現(xiàn)

課題介紹所謂數(shù)字濾波器，是指輸入輸出均為數(shù)字信號(hào)，通過(guò)數(shù)值運(yùn)算處理改變輸入信號(hào)所含頻率成分的相對(duì)比例，或者濾除某些頻率成分的數(shù)字器件或程序。數(shù)字濾波器處理精度高、穩(wěn)定、體積小、重量輕、靈活、不存在阻抗匹配問(wèn)題。橢圓濾波器具有較好的性能，本課結(jié)合MATLAB函數(shù)設(shè)計(jì)模擬橢圓濾波器。設(shè)計(jì)原理2.1常用濾波器比較典型的模擬濾波器有巴特沃斯(Butterworth)濾波器、切比雪夫(Chebyshev)濾波器和橢圓(Ellipse)濾波器等。其中，巴特沃斯濾波器又叫最平坦響應(yīng)濾波器，顧名思義，它的響應(yīng)最為平坦，通帶內(nèi)沒(méi)有波紋，其頻率響應(yīng)在通帶和阻帶中都是單調(diào)的，且在靠近零頻處最平坦，而在趨向阻帶時(shí)衰減單調(diào)增大，缺點(diǎn)是從通帶到阻帶的過(guò)渡帶寬，對(duì)于帶外干擾信號(hào)的衰減作用弱。切比雪夫?yàn)V波器又分為切比雪夫I型濾波器和切比雪夫II型濾波器。切比雪夫I型濾波器在整個(gè)通帶內(nèi)紋波最小，在阻帶內(nèi)隨頻率單調(diào)遞增；切比雪夫II型濾波器在通帶內(nèi)隨頻率光滑且單調(diào)遞增，零頻處最為平坦，在整個(gè)阻帶內(nèi)的紋波最小，它們的過(guò)渡帶較巴特沃斯濾波器陡峭。巴特沃斯濾波器和切比雪夫?yàn)V波器的傳輸函數(shù)都是一個(gè)常數(shù)除以一個(gè)多項(xiàng)式，為全極點(diǎn)網(wǎng)絡(luò)，所有的零點(diǎn)在無(wú)窮處，僅在無(wú)限大阻帶處衰減為無(wú)限大，而橢圓函數(shù)濾波器在有限頻率上既有零點(diǎn)又有極點(diǎn)。極零點(diǎn)在通帶內(nèi)產(chǎn)生等紋波，即它在整個(gè)通帶和阻帶上都具有最小的等紋波，這一點(diǎn)區(qū)別于在通帶和阻帶都平坦的巴特沃斯濾波器，以及通帶平坦、阻帶等波紋或是阻帶平坦、通帶等波紋的切比雪夫?yàn)V波器。同時(shí)，阻帶內(nèi)的有限傳輸零點(diǎn)減少了過(guò)渡區(qū)，可獲得極為陡峭的衰減曲線。

由圖可見(jiàn)，對(duì)于給定的階數(shù)和波紋要求，切比雪夫?yàn)V波器的選擇性優(yōu)于巴特沃斯濾波器，而橢圓函數(shù)濾波器優(yōu)于切比雪夫?yàn)V波器，它的過(guò)渡帶更窄，帶外抑制更加陡峭。橢圓濾波器的設(shè)計(jì)是基于橢圓有理函數(shù)，根據(jù)設(shè)計(jì)指標(biāo)由橢圓逼近的歸一化過(guò)程確定所需橢圓濾波器的階數(shù)及系統(tǒng)函數(shù)。Matlab是一套集數(shù)值計(jì)算、符號(hào)運(yùn)算及圖形處理等強(qiáng)大功能于一體的科學(xué)計(jì)算語(yǔ)言。作為強(qiáng)大的科學(xué)計(jì)算平臺(tái)，它幾乎能夠滿足所有的計(jì)算需求。Matlab的信號(hào)處理工具箱提供了設(shè)計(jì)橢圓濾波器的函數(shù)：ellipord函數(shù)和ellip函數(shù)。通過(guò)編程可以很容易由濾波器的技術(shù)指標(biāo)得到所需濾波器的階數(shù)，實(shí)現(xiàn)各種類型的橢圓濾波器，大大簡(jiǎn)化了橢圓濾波器的設(shè)計(jì)2.2橢圓濾波器的介紹橢圓濾波器(Ellipticfilter)又稱考爾濾波器(Cauerfilter),是在通帶和阻帶等波紋的一種濾波器。橢圓濾波器相比其他類型的濾波器,在階數(shù)相同的條件下有著最小的通帶和阻帶波動(dòng)。它在通帶和阻帶的波動(dòng)相同，這一點(diǎn)區(qū)別于在通帶和阻帶都平坦的巴特沃斯濾波器,以及通帶平坦、阻帶等波紋或是阻帶平坦、通帶等波紋的切比雪夫?yàn)V波器。幅值響應(yīng)在通帶和阻帶內(nèi)都是等波紋的，對(duì)于給定的階數(shù)和給定的波紋要求，橢圓濾波器能獲得較其它濾波器為窄的過(guò)渡帶寬，就這點(diǎn)而言，橢圓濾波器是最優(yōu)的，其振幅平方函數(shù)為H(jH(j⑵2a11+￡2R2(Q/Q)Np3-1)其中RN(x)是雅可比(Jacobi)橢圓函數(shù)，￡為與通帶衰減有關(guān)的參數(shù)。有上述分析可知橢圓濾波器有如下特點(diǎn)：橢圓帶通濾波器是一種零、極點(diǎn)型濾波器，它在有限頻率范圍內(nèi)存在傳輸零點(diǎn)和極點(diǎn)。橢圓帶通濾波器的通帶和阻帶都具有等波紋特性，因此通帶、阻帶逼近特性良好。對(duì)于同樣的性能要求，它比前兩種濾波器所需用的階數(shù)都低，而且它的過(guò)渡帶比較窄。2.3設(shè)計(jì)橢圓模擬帶通濾波器的步驟確定模擬濾波器的性能指標(biāo)：W,W,Rs,Rp。ps由性能指標(biāo)計(jì)算出濾波器階次N。通過(guò)歸一化及去歸一化求出的模擬濾波器Ha(s)。2.4模擬濾波器的MATLAB實(shí)現(xiàn)1.Matlab的信號(hào)處理軟件提供了設(shè)計(jì)橢圓濾波器的函數(shù)：ellipord函數(shù)和ellip函數(shù)。.Ellipord函數(shù)的功能是求濾波器的最小階數(shù)，其調(diào)用格式為：[n,Wn]二ellipord(Wp,Ws,Rp,Rs,'s') (4-1)其中各參量分別為：n-橢圓濾波器最小階數(shù)；Wp-橢圓濾波器通帶頻率；Ws-橢圓濾波器阻帶頻率；Rp-通帶波紋(dB)；Rs-阻帶最小衰減(dB)；本次設(shè)計(jì)中，用下面程序可確定濾波器階次：clearRp=1;Rs=50;Wp=[20003000]*2*pi;Ws=[17003300]*2*pi;[N,Wn]二ellipord(Wp,Ws,Rp,Rs,'s');運(yùn)行可得：N=6Wn二1.0e+004*1.2566 1.8850Ellip函數(shù)的功能是用來(lái)設(shè)計(jì)橢圓濾波器，其調(diào)用格式：[b,a]二ellip(n,Rp,Rs,Wp,'s') (4-2)[b,a]二ellip(n,Rp,Rs,Wp,'ftype') (4—3)其中：'ftype'=‘high'為高通濾波器，'ftype'='low'為低通濾波器，'ftype'='stop'為帶阻濾波器,s表示模擬濾波器的設(shè)計(jì)。返回長(zhǎng)度為n+1的濾波器系數(shù)行向量b和a,進(jìn)而求得傳遞函數(shù)H(z)：“B(z)b(l)+b⑵z-1+…b(n+1)z-nH= = zA(z) 1+a⑵z—1+…a(n+1)z—n3設(shè)計(jì)內(nèi)容3,1帶通濾波器實(shí)現(xiàn)的程序wp=2*pi*[20003000]； %性能指標(biāo)ws=2*pi*[17003300];Rp=1;Rs=50;[N,Wn]=ellipord(wp,ws,Rp,Rs,'s'); %估計(jì)濾波器的最小階數(shù)N[ba]=ellip(N,Rp,Rs,Wn,'s'); %利用ellip函數(shù)設(shè)計(jì)濾波器w=linspace(1,5000,1000)*2*piH=freqs(b,a,w);plot(w/(2*pi),20*log10(abs(H)); %畫(huà)出幅頻特性圖xlabel('Frequency(Hz)');ylabel('Magnitude(dB)');title(['ellipordAnalogbandpassfilterorderN=',num2str(N)]);gridon;運(yùn)行程序的輸出為：運(yùn)行可得：N=6Wn二1.0e+004*1.2566 1.8850

-100-120IlliI XI ]?' II II X, -100-120IlliI XI ]?' II II X, II I-DFig-1、FileEditViewInsertToolsDeslctopWindow旦elpI□呂日◎|氐|鷲竪紗⑥型銘彳忌|iSl匡i|■口0 500 1000 15002000 250030003500 40004500 5000Frequency(Hz)wp=2*pi*[20003000];ws=2*pi*[17003300];Rp=1;Rs=50;[N,Wn]=ellipord(wp,ws,Rp,Rs,'s');[ba]=ellip(N,Rp,Rs,Wn,'s');w=linspace(1,5000,1000)*2*pi;H=freqs(b,a,w);plot(w/(2*pi),20*log10(angle(H))/2*pi);xlabel('頻率')；ylabel('相位')；title([‘相位響應(yīng)N=',num2str(N)]);gridon;圖形如下：

3,2設(shè)計(jì)結(jié)果分析橢圓濾波器能得到較其它濾波器更窄的過(guò)渡帶寬，可以獲得對(duì)理想濾波器幅頻響應(yīng)的最好逼近，是一種性價(jià)比很高的濾波器。利用Matlab語(yǔ)言，其信號(hào)處理工具箱提供了豐富的設(shè)計(jì)方法，可以使得繁瑣的程序設(shè)計(jì)簡(jiǎn)化成函數(shù)的調(diào)用，只要以正確的指標(biāo)參數(shù)調(diào)用函數(shù)，就可以正確快捷地得到設(shè)計(jì)結(jié)果從而較方便地設(shè)計(jì)出橢圓濾波器。橢圓濾波器可以用較少的階數(shù)獲得很高的選擇特性，在設(shè)計(jì)過(guò)程中可以對(duì)比濾波器的特性，隨時(shí)更改橢圓濾波器通帶截止角頻率Wp,阻帶起始角頻率Ws,通帶波紋Rp,阻帶最小衰減Rs等參數(shù)，觀察濾波器的濾波效果。通過(guò)在設(shè)計(jì)中計(jì)算出的階次，考慮到實(shí)際應(yīng)用中計(jì)算機(jī)計(jì)算時(shí)間的限制，濾波器階次不可能過(guò)高,左右浮動(dòng)試驗(yàn)多次后得出6階為最佳階數(shù)。通過(guò)圖5-3,5-4可知,設(shè)計(jì)出的濾波器具有較理想的帶通性能,對(duì)高頻和低頻信號(hào)能有明顯的截止作用,同時(shí)也不會(huì)對(duì)中頻信號(hào)造成損失,很好的滿足了設(shè)計(jì)要求。4,心得體會(huì)通過(guò)一個(gè)星期的數(shù)字信號(hào)處理課程設(shè)計(jì),我對(duì)教材中所學(xué)知識(shí)有了更深的理解和認(rèn)識(shí)，教材中的基本定理和原理對(duì)我的設(shè)計(jì)起到了很好的指導(dǎo)作用，同時(shí)學(xué)習(xí)應(yīng)用了數(shù)字信號(hào)處理軟件Matlab，感嘆于其功能的強(qiáng)大性與掌握編程各種函數(shù)和語(yǔ)句的重要性。通過(guò)了解所要編程運(yùn)行的對(duì)象的原理，我學(xué)會(huì)了熟練運(yùn)用其基本功能。我覺(jué)得近一周的課程設(shè)計(jì)對(duì)我來(lái)說(shuō)，它不僅僅是讓我們把所學(xué)的理論知識(shí)與實(shí)踐相結(jié)合起來(lái)，提高自己的實(shí)際動(dòng)手能力和獨(dú)立思考的能力，更重要的是同學(xué)間的團(tuán)結(jié)，課程設(shè)計(jì)反映的是一個(gè)從理論到實(shí)際應(yīng)用的過(guò)程，但是更遠(yuǎn)一點(diǎn)可以聯(lián)系到畢業(yè)以后從學(xué)校到踏上社會(huì)的一個(gè)過(guò)程。同組同學(xué)的方案和建議使得我們的設(shè)計(jì)得以完善，在仿真中遇到的問(wèn)題也都得到了解決。在做本次課程設(shè)計(jì)的過(guò)程中，我感觸最深的當(dāng)屬查閱大量的設(shè)計(jì)資料了。為了讓自己的設(shè)計(jì)更加完善，查閱這方面的設(shè)計(jì)資料是十分必要的，同時(shí)也是必不可少的。另外在設(shè)計(jì)過(guò)程中也用到了許多高數(shù)等其他基礎(chǔ)課的知識(shí)，對(duì)以前的內(nèi)容有了更加綜合的掌握。由于首次進(jìn)行應(yīng)用數(shù)字信號(hào)處理的設(shè)計(jì)，很多知識(shí)還未做到靈活運(yùn)用的程度，如果以后有機(jī)會(huì)，我會(huì)繼續(xù)鍛煉自己的能力。參考文獻(xiàn)1］ 張德豐《數(shù)字信號(hào)處理與應(yīng)用》清華大學(xué)2］ 高西全,丁玉美.數(shù)字信號(hào)處理(第三版)［M］.西安:西安電子科技大學(xué)出版社,2008.3］ 朱義勝，董輝等譯?信號(hào)處理濾波器設(shè)計(jì)［M］.北京：電子工業(yè)出版社,2004.4］ 薛年喜.MATLAB在數(shù)字信號(hào)處理中的應(yīng)用［M］.北京:清華大學(xué)出版社,2003.5］ 肖有平，胡霞.高階橢圓濾波器的設(shè)計(jì)與仿真［J］.電子測(cè)量技術(shù),2007,(3).6］ 王靖,李永全.橢圓濾波器Matlb設(shè)計(jì)與實(shí)現(xiàn)［J］.現(xiàn)代電子技術(shù),2007,(6).目錄1課題介紹………………12設(shè)計(jì)原理……