實(shí)驗(yàn)四IIR數(shù)字濾波器的設(shè)計(jì)實(shí)驗(yàn)報(bào)告.doc
《實(shí)驗(yàn)四IIR數(shù)字濾波器的設(shè)計(jì)實(shí)驗(yàn)報(bào)告.doc》由會(huì)員分享,可在線閱讀,更多相關(guān)《實(shí)驗(yàn)四IIR數(shù)字濾波器的設(shè)計(jì)實(shí)驗(yàn)報(bào)告.doc(16頁(yè)珍藏版)》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
數(shù) 字 信 號(hào) 處 理 實(shí) 驗(yàn) 報(bào) 告 數(shù) 字 信 號(hào) 處 理 實(shí) 驗(yàn) 報(bào) 告 實(shí)驗(yàn)四 IIR數(shù)字濾波器的設(shè)計(jì) 學(xué)生姓名 張 志 翔 班級(jí) 電子信息工程1203班 學(xué)號(hào) 12401720522 指導(dǎo)教師 2015.4.29 實(shí)驗(yàn)四 IIR數(shù)字濾波器的設(shè)計(jì) 一、實(shí)驗(yàn)?zāi)康模? 1. 掌握雙線性變換法及脈沖響應(yīng)不變法設(shè)計(jì)IIR數(shù)字濾波器的具體設(shè)計(jì)方法及其原理,熟悉用雙線性變換法及脈沖響應(yīng)不變法設(shè)計(jì)低通、高通和帶通IIR數(shù)字濾波器的MATLAB編程。 2. 觀察雙線性變換及脈沖響應(yīng)不變法設(shè)計(jì)的濾波器的頻域特性,了解雙線性變換法及脈沖響應(yīng)不變法的特點(diǎn)。 3. 熟悉Butterworth濾波器、切比雪夫?yàn)V波器和橢圓濾波器的頻率特性。 二、實(shí)驗(yàn)原理: 1. 脈沖響應(yīng)不變法 用數(shù)字濾波器的單位脈沖響應(yīng)序列 模仿模擬濾波器的沖激響應(yīng) ,讓 正好等于 的采樣值,即 ,其中 為采樣間隔,如果以 及 分別表示 的拉式變換及 的Z變換,則 2.雙線性變換法 S平面與z平面之間滿足以下映射關(guān)系: s平面的虛軸單值地映射于z平面的單位圓上,s平面的左半平面完全映射到z平面的單位圓內(nèi)。 雙線性變換不存在混疊問(wèn)題。 雙線性變換是一種非線性變換 ,這種非線性引起的幅頻特性畸變可通過(guò)預(yù)畸而得到校正。 三、實(shí)驗(yàn)內(nèi)容及步驟: 實(shí)驗(yàn)中有關(guān)變量的定義: fc 通帶邊界頻率; fr阻帶邊界頻率;δ 通帶波動(dòng);At 最小阻帶衰減; fs采樣頻率; T采樣周期 (1) =0.3KHz, δ=0.8Db, =0.2KHz, At =20Db,T=1ms; 設(shè)計(jì)一個(gè)切比雪夫高通濾波器,觀察其通帶損耗和阻帶衰減是否滿足要求。 MATLAB源程序: wp=2*1000*tan(2*pi*300/(2*1000)); ws=2*1000*tan(2*pi*200/(2*1000)); [N,wn]=cheb1ord(wp,ws,0.8,20,s); %給定通帶(wp)和阻帶(ws)邊界角頻率,通帶波動(dòng)波動(dòng)0.8,阻帶最小衰減20dB,求出最低階數(shù)和通帶濾波器的通帶邊界頻率Wn [B,A]=cheby1(N,0.5,wn,high,s);%給定通帶(wp)和阻帶(ws)邊界角頻率,通帶波動(dòng) [num,den]=bilinear(B,A,1000); [h,w]=freqz(num,den); f=w/(2*pi)*1000; plot(f,20*log10(abs(h))); axis([0,500,-80,10]); grid;xlabel(頻率);ylabel(幅度/dB) 程序結(jié)果 num = 0.0304 -0.1218 0.1827 -0.1218 0.0304 den = 1 1.3834 1.4721 0.8012 0.2286 系統(tǒng)函數(shù): 幅頻響應(yīng)圖: 分析:由圖可知,切比雪夫?yàn)V波器幅頻響應(yīng)是通帶波紋,阻帶單調(diào)衰減的。δ=0.8,fr=0.2kHz,At=30Db,滿足設(shè)計(jì)要求 (2)fc=0.2kHz, δ=1dB,fr=0.3kHz,At=25dB,T=1ms;分別用脈沖響應(yīng)不變法及雙線性變換法設(shè)計(jì)一Butterworth數(shù)字低通濾波器,觀察所設(shè)計(jì)數(shù)字濾波器的幅頻特性曲線,記錄帶寬和衰減量,檢查是否滿足要求。比較這兩種方法的優(yōu)缺點(diǎn)。 MATLAB源程序: T = 0.001;fs = 1000;fc = 200;fr = 300; wp1 = 2*pi*fc;wr1 = 2*pi*fr; [N1,wn1] = buttord(wp1,wr1,1,25,s) [B1,A1] = butter(N1,wn1,s); [num1,den1] = impinvar(B1,A1,fs);%脈沖響應(yīng)不變法 [h1,w] = freqz(num1,den1); wp2 = 2*fs*tan(2*pi*fc/(2*fs)) wr2 = 2*fs*tan(2*pi*fr/(2*fs)) [N2,wn2] = buttord(wp2,wr2,1,25,s) [B2,A2] = butter(N2,wn2,s); [num2,den2] = bilinear(B2,A2,fs);%雙線性變換法 [h2,w] = freqz(num2,den2); f = w/(2*pi)*fs; plot(f,20*log10(abs(h1)),-.,f,20*log10(abs(h2)),-); axis([0,500,-100,10]);grid;xlabel(頻率/Hz );ylabel(幅度/dB) title(巴特沃思數(shù)字低通濾波器); legend(脈沖相應(yīng)不變法,雙線性變換法,1); 結(jié)果分析: 脈沖響應(yīng)不變法的低通濾波器系統(tǒng)函數(shù): num1 -2.3647 0.0002 0.0153 0.0995 0.1444 0.0611 0.0075 0.0002 3.6569 0 den1 1 -1.9199 2.5324 -2.2053 1.3868 -0.6309 0.2045 -0.0450 0.0060 -0.0004 雙線性變換法設(shè)計(jì)的低通濾波器系統(tǒng)函數(shù): num2 0.0179 0.1072 0.2681 0.3575 0.2681 0.1072 0.0179 den2 1 -0.6019 0.9130 -0.2989 0.1501 -0.0208 0.0025 分析:脈沖響應(yīng)不變法的頻率變化是線性的,數(shù)字濾波器頻譜響應(yīng)出現(xiàn)了混疊,影響了過(guò)渡帶的衰減特性,并且無(wú)傳輸零點(diǎn);雙線性變化法的頻率響應(yīng)是非線性的,因而消除了頻譜混疊,在f=500Hz出有一個(gè)傳輸零點(diǎn)。 脈沖響應(yīng)不變法的一個(gè)重要特點(diǎn)是頻率坐標(biāo)的變換是線性的,ω=ΩΤ,ω與Ω是線性關(guān)系:在某些場(chǎng)合,要求數(shù)字濾波器在時(shí)域上能模仿模擬濾波器的功能時(shí),如要實(shí)現(xiàn)時(shí)域沖激響應(yīng)的模仿,一般使用脈沖響應(yīng)不變法。 脈沖響應(yīng)不變法的最大缺點(diǎn):有頻譜周期延拓效應(yīng),因此只能用于帶限的頻響特性,如衰減特性很好的低通或帶通,而高頻衰減越大,頻響的混淆效應(yīng)越小,至于高通和帶阻濾波器,由于它們?cè)诟哳l部分不衰減,因此將完全混淆在低頻響應(yīng)中,此時(shí)可增加一保護(hù)濾波器,濾掉高于 的頻帶,再用脈沖響應(yīng)不變法轉(zhuǎn)換為數(shù)字濾波器,這會(huì)增加設(shè)計(jì)的復(fù)雜性和濾波器階數(shù),只有在一定要滿足頻率線性關(guān)系或保持網(wǎng)絡(luò)瞬態(tài)響應(yīng)時(shí)才采用。 雙線性變換法的主要優(yōu)點(diǎn)是S平面與Z平面一一單值對(duì)應(yīng),s平面的虛軸(整個(gè)jΩ)對(duì)應(yīng)于Z平面單位圓的一周,S平面的Ω=0處對(duì)應(yīng)于Z平面的ω=0處, Ω= ∞處對(duì)應(yīng)于Z平面的ω= π處,即數(shù)字濾波器的頻率響應(yīng)終止于折疊頻率處,所以雙線性變換不存在混迭效應(yīng)。 雙線性變換缺點(diǎn): Ω與ω成非線性關(guān)系,導(dǎo)致: a. 數(shù)字濾波器的幅頻響應(yīng)相對(duì)于模擬濾波器的幅頻響應(yīng)有畸變,(使數(shù)字濾波器與模擬濾波器在響應(yīng)與頻率的對(duì)應(yīng)關(guān)系上發(fā)生畸變)。 b. 線性相位模擬濾波器經(jīng)雙線性變換后,得到的數(shù)字濾波器為非線性相位。 c.要求模擬濾波器的幅頻響應(yīng)必須是分段恒定的,故雙線性變換只能用于設(shè)計(jì)低通、高通、帶通、帶阻等選頻濾波器。 (3)利用雙線性變換法分別設(shè)計(jì)滿足下列指標(biāo)的Butterworth型、Chebyshev型和橢圓型數(shù)字低通濾波器,并作圖驗(yàn)證設(shè)計(jì)結(jié)果:fc=1.2kHz ,δ≤0.5dB ,fr=2kHz , At≥40dB, fs=8kHz,比較這種濾波器的階數(shù)。 MATLAB源程序: clear all; wc=2*pi*1200;wr=2*pi*2000;rp=0.5;rs=40;fs=8000; w1=2*fs*tan(wc/(2*fs)); w2=2*fs*tan(wr/(2*fs)); [Nb,wn]=buttord(w1,w2,rp,rs,s) %巴特沃思 [B,A]=butter(Nb,wn,s); [num1,den1]=bilinear(B,A,fs); [h1,w]=freqz(num1,den1); [Nc,wn]=cheb1ord(w1,w2,rp,rs,s) %切比雪夫 [B,A]=cheby1(Nc,rp,wn,s); [num2,den2]=bilinear(B,A,fs); [h2,w]=freqz(num2,den2); [Ne,wn]=ellipord(w1,w2,rp,rs,s) %橢圓型 [B,A]=ellip(Ne,rp,rs,wn,low,s); [num3,den3]=bilinear(B,A,fs); [h3,w]=freqz(num3,den3); f=w/(2*pi)*fs; plot(f,20*log10(abs(h1)),-,f,20*log10(abs(h2)),--,f,20*log10(abs(h3)),:); axis([0,4000,-100,10]);grid; xlabel(Frequency in Hz); ylabel(Gain in dB); title(三種數(shù)字低通濾波器); legend(巴特沃思數(shù)字低通濾波器,切比雪夫數(shù)字低通濾波器,橢圓數(shù)字低通濾波器,3); 巴特沃思數(shù)字低通濾波器的系統(tǒng)函數(shù)系數(shù): num1= 0.0032 0.0129 0.0302 0.0453 0.0453 0.0302 0.0129 0.0032 0.0003 den1= -2.7996 4.4581 -4.5412 3.2404 -1.6330 0.5780 -0.1370 0.0197 -0.0013 切比雪夫數(shù)字低通濾波器的系統(tǒng)函數(shù)系數(shù): num2= 0.0026 0.0132 0.0264 0.0264 0.0132 0.0026 den2= 1 -2.9775 4.2932 -3.5124 1.6145 -0.3334 橢圓數(shù)字低通濾波器的系統(tǒng)函數(shù)系數(shù): num3= 0.03887 0.0363 0.0665 0.0363 0.0389 den3= 1 -2.1444 2.3658 -1.3250 0.3332 程序結(jié)果圖: 分析:設(shè)計(jì)結(jié)果表明,巴特沃思數(shù)字低通濾波器、切比雪夫數(shù)字低通濾波器、橢圓數(shù)字低通濾波器的階數(shù)分別是9、5、4階??梢?jiàn),對(duì)于給定的階數(shù),橢圓數(shù)字低通濾波器的階數(shù)最少(換言之,對(duì)于給定的階數(shù),過(guò)渡帶最窄),就這一點(diǎn)來(lái)說(shuō),他是最優(yōu)濾波器。由圖表明,巴特沃思數(shù)字低通濾波器過(guò)渡帶最寬,幅頻響應(yīng)單調(diào)下降;橢圓數(shù)字低通濾波器過(guò)渡帶最窄,并具有等波紋的通帶和阻帶響應(yīng);切比雪夫數(shù)字低通濾波器的過(guò)渡帶介于兩者之間。 (4)分別用脈沖響應(yīng)不變法及雙線性變換法設(shè)計(jì)一Butterworth型數(shù)字帶通濾波器,已知 ,其等效的模擬濾波器指標(biāo) δ<3dB,2kHz- 1.請(qǐng)仔細(xì)閱讀文檔,確保文檔完整性,對(duì)于不預(yù)覽、不比對(duì)內(nèi)容而直接下載帶來(lái)的問(wèn)題本站不予受理。
- 2.下載的文檔,不會(huì)出現(xiàn)我們的網(wǎng)址水印。
- 3、該文檔所得收入(下載+內(nèi)容+預(yù)覽)歸上傳者、原創(chuàng)作者;如果您是本文檔原作者,請(qǐng)點(diǎn)此認(rèn)領(lǐng)!既往收益都?xì)w您。
下載文檔到電腦,查找使用更方便
9.9 積分
下載 |
- 配套講稿:
如PPT文件的首頁(yè)顯示word圖標(biāo),表示該P(yáng)PT已包含配套word講稿。雙擊word圖標(biāo)可打開(kāi)word文檔。
- 特殊限制:
部分文檔作品中含有的國(guó)旗、國(guó)徽等圖片,僅作為作品整體效果示例展示,禁止商用。設(shè)計(jì)者僅對(duì)作品中獨(dú)創(chuàng)性部分享有著作權(quán)。
- 關(guān) 鍵 詞:
- 實(shí)驗(yàn) IIR 數(shù)字濾波器 設(shè)計(jì) 報(bào)告
鏈接地址:http://www.hcyjhs8.com/p-6535475.html