matlab濾波器代碼_第1頁(yè)
matlab濾波器代碼_第2頁(yè)
matlab濾波器代碼_第3頁(yè)
matlab濾波器代碼_第4頁(yè)
matlab濾波器代碼_第5頁(yè)
已閱讀5頁(yè)，還剩3頁(yè)未讀， 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說(shuō)明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、數(shù)字信號(hào)處理：已知通帶截止頻率fp=5kHz，通帶最大衰減ap=2dB，阻帶截止頻率fs=2kHz，阻帶最小衰減as=30dB，按照以上技術(shù)指標(biāo)設(shè)計(jì)巴特沃斯低通濾波器：wp=2*pi*5000;ws=2*pi*12000;Rp=2;As=30;N,wc=buttord(wp,ws,Rp,As,'s');B,A=butter(N,wc,'s');k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;Hk=freqs(B,A,wk);subplot(2,2,1);plot(fk/1000,20*log10(abs(Hk);grid onxl

2、abel('頻率（kHz）');ylabel('幅度（dB）')axis(0,14,-40,5)切比雪夫1型低通濾波器：wp=2*pi*3000;ws=2*pi*12000;Rp=0.1;As=60;N1,wpl=cheb1ord(wp,ws,Rp,As,'s'); %cheb1ord,里面的是1，不是LB1,A1=cheby1(N1,Rp,wpl,'s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;Hk=freqs(B1,A1,wk);plot(fk/1000,20*log10

3、(abs(Hk);grid onxlabel('頻率（kHz）');ylabel('幅度（dB）')axis(0,12,-70,5)橢圓模擬低通濾波器：wp=2*pi*3000;ws=2*pi*12000;Rp=0.1;As=60;N,wpo=ellipord(wp,ws,Rp,As,'s');B,A=ellip(N,Rp,As,wpo,'s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;Hk=freqs(B1,A1,wk);plot(fk/1000,20*log10(abs(H

4、k);grid onxlabel('頻率（kHz）');ylabel('幅度（dB）')axis(0,12,-70,5)p195-14wp=2*4/80;ws=2*20/80;rp=0.5;rs=45;N,wc=buttord(wp,ws,rp,rs);B,A=butter(N,wc);hk,wk=freqz(B,A);fk=wk/pi*40;plot(fk,20*log10(abs(hk);axis(0,30,-100,0)xlabel('頻率(kHZ)');ylabel('幅度(db)');grid onP195-16wp=2

5、*325/2500;ws=2*225/2500;rp=1;rs=40;N,wc=ellipord(wp,ws,rp,rs);B,A=ellip(N,rp,rs,wc);hk,wk=freqz(B,A);fk=wk/pi*40;plot(fk,20*log10(abs(hk);axis(0,30,-100,0)xlabel('頻率(kHZ)');ylabel('幅度(db)');grid onP195-15wp=2*4/80;ws=2*20/80;rp=0.5;rs=45;N,wc=cheb1ord(wp,ws,rp,rs);B,A=cheby1(N,rp,wc)

6、;hk,wk=freqz(B,A);fk=wk/pi*40;plot(fk,20*log10(abs(hk);axis(0,30,-100,0)xlabel('頻率(kHZ)');ylabel('幅度(db)');grid on切比雪夫低通濾波器wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;N1,wp1=cheb1ord(wp,ws,rp,as,'s');B1,A1=cheby1(N1,rp,wp1,'s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*

7、fk;hk=freqs(B1,A1,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel('頻率(kHZ)');ylabel('幅度(db)');axis(0,12,-70,5)雙音頻檢測(cè)audiofile='test.wav'in_audio,fs,bits=wavread(audiofile);b,a=cheby1(5,0.1,0.3);out_audio=filter(b,a,in_audio);sound(out_audio,fs,bits);wavwrite(out_audio,fs,bits,

8、'test_out');xk1=fft(in_audio,512);xk2=fft(out_audio,512);subplot(2,1,1);stem(abs(xk1);subplot(2,1,2);stem(abs(xk2);巴特沃斯模擬高通濾波器。wp=2*pi*4000;ws=2*pi*1000;rp=0.1;as=40;n,wc=buttord(wp,ws,rp,as,'s');bh,ah=butter(n,wc,'high','s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*p

9、i*fk;hk=freqs(bh,ah,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel('頻率(khz)');ylabel('幅度(db)')axis(0,12,-70,5)巴特沃斯模擬低通濾波器。wp=2*pi*5000;ws=2*pi*12000;rp=2;as=30;n,wc=buttord(wp,ws,rp,as,'s');b,a=butter(n,wc,'s');k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;hk=freqs(b,a,wk)

10、;subplot(2,2,1);plot(fk/1000,20*log10(abs(hk);grid onxlabel('頻率(khz)');ylabel('幅度(db)')axis(0,14,-40,5)巴特沃斯模擬帶通濾波器wp=2*pi*4000,7000;ws=2*pi*2000,9000;n,wc=buttord(wp,ws,rp,as,'s');bb,ab=butter(n,wc,'s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(bb,ab,wk)

11、;plot(fk/1000,20*log10(abs(hk);grid onxlabel('頻率(khz)');ylabel('幅度(db)')axis(0,12,-70,5)巴特沃斯帶阻濾波器wp=2*pi*2000,9000;ws=2*pi*4000,7000;rp=1;as=20;nb,wc=buttord(wp,ws,rp,as,'s');bsb,asb=butter(nb,wc,'stop','s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=fr

12、eqs(bsb,asb,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel('頻率(khz)');ylabel('幅度(db)')axis(0,12,-70,5)第14題。巴特沃斯低通濾波器（采樣頻率80，同代便捷頻率4，最大衰減0.5）wp=2*4/80;ws=2*20/80;rp=0.5;rs=45; N,wc=buttord(wp,ws,rp,rs);B,A=butter(N,wc);hk,wk=freqz(B,A);fk=wk/pi*40;plot(fk,20*log10(abs(hk);axis(0,30,-

13、100,0);xlabel('頻率(kHz)');ylabel('幅度(dB)')設(shè)計(jì)切比雪夫I型和切比雪夫II型模擬低通濾波器。 I型：wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;n1,wp1=cheb1ord(wp,ws,rp,as,'s');b1,a1=cheby1(n1,rp,wp1,'s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(b1,a1,wk);plot(fk/1000,20*log10(abs(hk);gr

14、id onxlabel('頻率(khz)');ylabel('幅度(db)')axis(0,12,-70,5)II型：wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;n2,wp2=cheb2ord(wp,ws,rp,as,'s');b2,a2=cheby2(n2,rp,wp2,'s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(b2,a2,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabe

15、l('頻率(khz)');ylabel('幅度(db)')axis(0,12,-70,5)切比雪夫I型高通濾波器wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;n1,wp1=cheb1ord(wp,ws,rp,as,'s');b1,a1=cheby1(n1,rp,wp1,'HIGH','s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(b1,a1,wk);plot(fk/1000,20*log10(abs(hk);

16、grid onxlabel('頻率(khz)');ylabel('幅度(db)')axis(0,12,-70,5)設(shè)計(jì)橢圓模擬低通濾波器。wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;n,wpo=ellipord(wp,ws,rp,as,'s');b,a=ellip(n,rp,as,wpo,'s');subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(b,a,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel('頻率(khz)');ylabel('幅度(db)')axis(0,12,-70,5)橢圓帶阻濾波器wp=2*pi*2000,9000;ws=2*pi*4000,7000;rp=1;as=20;ne,wep=e

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間，僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請(qǐng)與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論