切比雪夫低通濾波器設(shè)計(jì)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上課程設(shè)計(jì)任務(wù)書20102011學(xué)年第一學(xué)期專業(yè)： 通信工程 學(xué)號(hào)： 姓名： 課程設(shè)計(jì)名稱： 數(shù)字信號(hào)處理課程設(shè)計(jì) 設(shè)計(jì)題目： 切比雪夫型低通濾波器的軟件實(shí)現(xiàn) 完成期限：自 年 月 日至 年 月 日共 周一設(shè)計(jì)目的1鞏固所學(xué)的理論知識(shí)。2提高綜合運(yùn)用所學(xué)理論知識(shí)獨(dú)立分析和解決問(wèn)題的能力。3更好地將理論與實(shí)踐相結(jié)合。4掌握信號(hào)分析與處理的基本方法與實(shí)現(xiàn)。5熟練使用MATLAB語(yǔ)言進(jìn)行編程實(shí)現(xiàn)。二設(shè)計(jì)內(nèi)容用MATLAB編程設(shè)計(jì)切比雪夫I型低通濾波器，各參數(shù)要求如下：fp=100 Hz，p=2 dB，fs=120 Hz，s=60 dB。

2、三、設(shè)計(jì)要求給出所設(shè)計(jì)濾波器的幅度及幅度衰減特性并分析是否滿足設(shè)計(jì)需要。四設(shè)計(jì)條件計(jì)算機(jī)、MATLAB語(yǔ)言環(huán)境五、參考資料1數(shù)字信號(hào)處理（第三版），丁玉美，高西全.西安電子科技大學(xué)出版社，2000.2MATLAB及在電子信息課程中的應(yīng)用，陳懷堔，吳大正，高西全.電子工業(yè)出版社，2006.3MATLAB 7.0從入門到精通，求是科技.人民郵電出版社，2006.4數(shù)字信號(hào)處理（第三版）學(xué)習(xí)指導(dǎo)，高西全，丁玉美.西安科技大學(xué)出版社，2001.指導(dǎo)教師（簽字）： 教研室主任（簽字）： 批準(zhǔn)日期： 年 月 日摘要隨著信息和數(shù)字時(shí)代的到來(lái)，數(shù)字信號(hào)處理已成為當(dāng)今一門極其重要的學(xué)科和技術(shù)領(lǐng)域。在現(xiàn)代通信系統(tǒng)

3、中，由于信號(hào)中經(jīng)常混有各種復(fù)雜成分，因此很多信號(hào)的處理都是基于濾波器而進(jìn)行的。所以，數(shù)字濾波器在數(shù)字信號(hào)處理中起著舉足輕重的作用。而數(shù)字濾波器的設(shè)計(jì)都要以模擬濾波器為基礎(chǔ)的，這是因?yàn)槟M濾波器的理論和設(shè)計(jì)方方法都已發(fā)展的相當(dāng)成熟， 且有典型的模擬濾波器供我們選擇。，如巴特沃思濾波器、切比雪夫?yàn)V波器等。本次課程設(shè)計(jì)將運(yùn)用MATLAB設(shè)計(jì)一個(gè)基于切比雪夫低通濾波器，并出所設(shè)計(jì)濾波器的幅度及幅度衰減特性。關(guān)鍵詞：模擬低通 濾波 切比雪夫目錄專心-專注-專業(yè)1課題描述數(shù)字濾波器是數(shù)字信號(hào)處理的重要工具之一，它通過(guò)數(shù)值運(yùn)算處理改變輸入信號(hào)所含頻率成分的相對(duì)比例或者濾出某些頻率成分的數(shù)字器件或程序，而數(shù)字

4、濾波器處理精度高、體積小、穩(wěn)定、重量輕、靈活、不存在阻抗匹配問(wèn)題，可以實(shí)現(xiàn)模擬濾波器無(wú)法實(shí)現(xiàn)的特殊功能。故本課題使用MATLAB信號(hào)處理箱和運(yùn)用切比雪夫法設(shè)計(jì)數(shù)字低通濾波器。 2設(shè)計(jì)原理2.1切比雪夫?yàn)V波器介紹在巴特沃茲濾波器中，幅度響應(yīng)在通帶和阻帶內(nèi)都是單調(diào)的。因此，若濾波器的技術(shù)要求是用最大通帶和阻帶的逼近誤差來(lái)給出的話，那么，在靠近通帶低頻端和阻帶截止頻率以上的部分都會(huì)超出技術(shù)指標(biāo)。一種比較有效的途徑是使逼近誤差均勻地分布于通帶或阻帶內(nèi)，或同時(shí)在通帶和阻帶內(nèi)都均勻分布，這樣往往可以降低所要求的濾波器階次。通過(guò)選擇一種具有等波紋特性而不是單調(diào)特性的逼近方法可以實(shí)現(xiàn)這一點(diǎn)。切比雪夫型濾波器就

5、具有這種性質(zhì)：其頻率響應(yīng)的幅度既可以在通帶中是等波紋的，而在阻帶中是單調(diào)的（稱為I型切比雪夫?yàn)V波器），也可以在通帶中是單調(diào)的，而在阻帶中是等波紋的（稱為II型切比雪夫?yàn)V波器）。I型切比雪夫?yàn)V波器的幅度平方函數(shù)是 = (2.1)式中為N階切比雪夫多項(xiàng)式，定義為 (2.2)從定義切比雪夫多項(xiàng)式可以直接得出由和求的遞推公式。將三角恒等式代入 (2.2)式，得 =2x (2.3)從 (2.2)式我們注意到，當(dāng)0<x<1時(shí)，在0和1之間變化；當(dāng)x>1時(shí)，是虛數(shù)，所以像雙曲余弦一樣單調(diào)地增加。參考(2.1)，對(duì)于01呈現(xiàn)出在1和1/（）之間的波動(dòng)；而對(duì)于1單調(diào)地減小。需要用三個(gè)參量來(lái)確定

6、該濾波器：，和N。在典型的設(shè)計(jì)中，用容許的通帶波紋來(lái)確定，而用希望的通帶截止頻率來(lái)確定。然后選擇合適的階次N，以便阻帶的技術(shù)要求得到滿足。定義允許的通帶最大衰減用下式表示： 為了求切比雪夫?yàn)V波器在橢圓上極點(diǎn)的位置，我們首先要這樣確定，在大圓和小圓上以等角度 等間隔排列的那些點(diǎn)：這些點(diǎn)對(duì)于虛軸呈對(duì)稱分布，并且沒(méi)有一個(gè)點(diǎn)落在虛軸上；但當(dāng)N為奇數(shù)時(shí)要有一個(gè)點(diǎn)落在實(shí)軸上，而當(dāng)N為偶數(shù)時(shí)，就都不會(huì)落在實(shí)軸上。切比雪夫?yàn)V波器的極點(diǎn)落在橢圓上，起縱坐標(biāo)由相應(yīng)的大圓上點(diǎn)的縱坐標(biāo)來(lái)表示，起橫坐標(biāo)由相應(yīng)的小圓上點(diǎn)的橫坐標(biāo)來(lái)表示。1.1 濾波器的分類（1）從功能上分；低、帶、高、帶阻。 （2）從實(shí)現(xiàn)方法上分:FIR

7、、IIR （3）從設(shè)計(jì)方法上來(lái)分：Chebyshev(切比雪夫）,Butterworth（巴特沃斯） （4）從處理信號(hào)分：經(jīng)典濾波器、現(xiàn)代濾波器 1.2 模擬濾波器的設(shè)計(jì)指標(biāo)設(shè)是一個(gè)模擬濾波器的頻率響應(yīng)，則基于平方幅度響應(yīng)的低通濾波器技術(shù)指標(biāo)為： （2-2-1） （2-2-2）其中為通帶波動(dòng)系數(shù)，和是通帶和阻帶邊緣頻率。A為阻帶衰減系數(shù)。這些指標(biāo)如圖所示。從圖知必須滿足 （2-2-3）其中參數(shù)和A是數(shù)字濾波器指標(biāo)。1.3 切比雪夫I型濾波器1.3.1 切比雪夫低通濾波器的設(shè)計(jì)原理切比雪夫?yàn)V波器的幅頻特性具有等波紋特性。它有兩種形式：振幅特性在通帶內(nèi)是等波紋的，在阻帶內(nèi)是單調(diào)遞減的切比雪夫I型濾

8、波器，振幅特性在阻帶內(nèi)是等波紋的，在通帶內(nèi)是單調(diào)遞減的切比雪夫II型濾波器，如圖所示分別畫出了濾波器的幅頻特性和衰減函數(shù)。以切比雪夫I型為例介紹其設(shè)計(jì)原理幅度平方函數(shù)用表示 （2-3-1）式中，為小于1的正數(shù)，表示通帶內(nèi)幅度波動(dòng)的程度，越大，波動(dòng)幅度也越大。稱為通帶截止頻率。令 ，稱為對(duì)的歸一化頻率。定義允許的通帶內(nèi)最大衰減用下式表示 （2-3-2） 式中 （2-3-3） （2-3-4） 因此 （2-3-5） （2-3-6） 這樣，可以根據(jù)通帶內(nèi)最大衰減，可求出參數(shù)。 階數(shù)N影響過(guò)渡帶的寬度，同時(shí)也影響通帶內(nèi)波動(dòng)的疏密，因?yàn)镹等于通帶內(nèi)的最大值和最小值的總個(gè)數(shù)。設(shè)阻帶的起點(diǎn)頻率為，則有 （2-

9、3-7） 令，由>1，有 （2-3-8） 可以解出 （2-3-9） （2-3-10）3dB截止頻率用表示， （2-3-11） 按照（2-3-1）式，有 （2-3-12） 經(jīng)過(guò)一系列推論得歸一化系統(tǒng)函數(shù)為（2-3-13）去歸一化的系統(tǒng)函數(shù)為 （2-3-14） 1.3.2 切比雪夫低通濾波器的設(shè)計(jì)步驟（1）確定低通濾波器的技術(shù)指標(biāo)：邊帶頻率，通帶最大衰減、阻帶最大衰減、阻帶截至頻率，它們滿足 （2-3-15）（2-3-16）（2）求濾波器階數(shù)N和參數(shù) （2-3-17） （2-3-18） （2-3-19）這樣，先由（2-18）式求出，代入（2-3-19），求出階數(shù)N，最后取大于或等于N的最小整

10、數(shù)。（3）求歸一化系統(tǒng)函數(shù) （2-3-20）（4）將去歸一化，得到實(shí)際的 （2-3-21）2.3.3用MATLAB設(shè)計(jì)切比雪夫低通濾波器 （1）該格式用于計(jì)算N階切比雪夫I型歸一化模擬低通濾波器系統(tǒng)的零極點(diǎn)和增益因子。返回長(zhǎng)度為N的列向量Z和P，分別給出N個(gè)零點(diǎn)和極點(diǎn)的位置。Rs是阻帶最小衰減。（2）該格式用于計(jì)算切比雪夫I型數(shù)字濾波器的階數(shù)N和阻帶截止頻率wso。調(diào)用參數(shù)分別為數(shù)字濾波器的通帶頻率和阻帶邊界頻率的歸一化值。（3） 該格式用于計(jì)算切比雪夫I型模擬濾波器的階數(shù)N和阻帶截止頻率wso。wp和ws是實(shí)際模擬角頻率。（4） 該格式用于計(jì)算N階切比雪夫I型數(shù)字濾波器系統(tǒng)函數(shù)的分子和分母多

11、項(xiàng)式系數(shù)向量B和A。調(diào)用參數(shù)N和wso分別為切比雪夫I型數(shù)字濾波器的階數(shù)和阻帶截止頻率的歸一化值。（5） 該格式用于計(jì)算N階切比雪夫I型模擬濾波器系統(tǒng)函數(shù)的分子和分母多項(xiàng)式系數(shù)向量B和A。調(diào)用參數(shù)N和wso分別為切比雪夫I型模擬濾波器的階數(shù)和阻帶截止頻率的歸一化值。 2 設(shè)計(jì)內(nèi)容2.1 用MATLAB編程實(shí)現(xiàn)編程如下：wp=2*pi*100;ws=2*pi*120;Rp=2;As=60;%設(shè)置指標(biāo)參數(shù)N1,wp1=cheb1ord(wp,ws,Rp,As,'s')%計(jì)算切比雪夫I型模擬低通濾波器的階數(shù)和通帶邊%界頻率B1,A1=cheby1(N1,Rp,wp1,'s&#

12、39;); %計(jì)算切比雪夫I型模擬低通濾波器系統(tǒng)函數(shù)系數(shù)fk=0:120;wk=2*pi*fk;Hk=freqs(B1,A1,wk);%模擬低通濾波器的系統(tǒng)函數(shù)subplot(2,2,1);plot(fk/10,abs(Hk);%畫出幅頻特性xlabel('頻率(Hz)');ylabel('幅度');subplot(2,2,2);plot(fk/10,20*log10(abs(Hk);grid on;%畫出衰減函數(shù)xlabel('頻率(Hz)');ylabel('幅度(dB)');axis(0,15,-50,5);圖形如下：2.2

13、 設(shè)計(jì)結(jié)果分析（1）求得階數(shù)N1 =13，通帶邊界頻率為wp1=628.3185，表明階數(shù)越大，模擬低通濾波器越理想，且接近矩形。（2）從圖上可以看出：切比雪夫I型濾波器是在通帶頻率響應(yīng)幅度等波紋的濾波器，通帶內(nèi)具有等波紋起伏特性，阻帶內(nèi)單調(diào)下降且衰減更大；幅度在低頻處通過(guò)，隨著頻率的增加，逐漸衰減至零；綜上所看：本設(shè)計(jì)符合要求。3 總結(jié)本次課程設(shè)計(jì)，讓我對(duì)設(shè)計(jì)濾波器的原理和步驟有了更深入的理解。濾波器在數(shù)字信號(hào)處理中占據(jù)著重要的地位，因此掌握模擬低通濾波器的設(shè)計(jì)對(duì)以后的學(xué)習(xí)有著不菲的收獲。首先，設(shè)計(jì)濾波器前一定要選好設(shè)計(jì)的方法。這次課程設(shè)計(jì)采用切比雪夫I型來(lái)設(shè)計(jì)一個(gè)模擬低通濾波器。其次，確定

14、好參數(shù)。最后，用MATLAB軟件編程式設(shè)計(jì)過(guò)程大為簡(jiǎn)化。4 參考文獻(xiàn)1數(shù)字信號(hào)處理（第三版），丁玉美，高西全.西安電子科技大學(xué)出版社，2000.2MATLAB及在電子信息課程中的應(yīng)用，陳懷堔，吳大正，高西全.電子工業(yè)出版社，2006.3MATLAB 7.0從入門到精通，求是科技.人民郵電出版社，2006.4數(shù)字信號(hào)處理（第三版）學(xué)習(xí)指導(dǎo)，高西全，丁玉美.西安科技大學(xué)出版社，2001.T=0.1;Fs=1/T;fp=3;fs=4;Rp=1;As=15;wp=fp/Fs*2*pi;ws=fs/Fs*2*pi;Wp=(2/T)*tan(wp/2);Ws=(2/T)*tan(ws/2);%設(shè)計(jì)切比雪夫低通濾波器原型ep=sqrt(10(Rp/10)-1);A=10(As/20);Wc=Wp;Wr=Ws/Wp;g=sqrt(A*A-1/ep);N=cei