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

下載本文檔

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

文檔簡介

1、數(shù)字信號處理： 已知通帶截止頻率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 on xlabel( 頻率(

2、 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，不是 L B1,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(abs(Hk);grid on xlabel( 頻率( kHz);

3、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(Hk);grid on xlabel( 頻率( kHz);ylabel(幅度( dB) axis(0,12,-70,5) 295

4、 二 4 WPH2*4oo0 八 WSH2*20oo0HPH0.5HSH45 八 LZWCllbusorcKWpwsdprsx pAHbuffelxzwc)八 Lhk-wkllfreqz(B- fkuwk、pi*40 八 P_of(fk20*oglo(abs(hk)八 axis(p3oxoooD x_abe_ (-酋4M(kHzmy_abe_ (-B岡(d by)込 rid on WPH2*325K5500八 WSH2*225K5500HPH 1HSH40 八 N 一 wcHe 三 pod(wpws IP Is)八 pAHe三p(Ndpswc= Lhk-wkllfreqz(B- fkuwk、p

5、i*40 八 P_of(fk20*oglo(abs(hk)八 axis(p3oxoooD x_abe_ (-酋4M(kHzmy_abe_ (-B岡(d by)込 rid on P195-15 頻率(kHZ) wp=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); 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

6、 切比雪夫低通濾波器 頻率(kHZ) 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);grid on xlabel(頻率(kHZ);ylabel( 幅度(db); axis(0,12,-70,5) 雙音頻檢測 audiofile=test.wav =

7、nlaudipfs-bifslrwawead(audiofi-e= bbHchebylpp103)八 ouflaudioHfi_(D(bbjnlaudio)八 sound(ouflaudipfskfs)八 Wavwme(ouflaudipf3bifs=esflour= xkl Hffwnlaudioull 2)八 xk2H3(ouflaudip512)八 subp-0f(21- fkuwk、pi*40 八 P_of(fk20*oglo(abs(hk)八 axis(p3oxooo2 x_abe三 s (kHZ)y_abe三B岡(dBy) 一B啓s民w-Bs吏寶ffi話藩器 一 B- WPH2*p

8、i*3000八 WSH2*pi*l 2000HPHP1Q5SH60 八 三 wp lHchebod(wpws IP bss_)八 LbkQilllchebwml mwpk運(yùn))八 subp-0f(223八 fkHpk200552=200pwkH2*pi4;?T hkHfeqs(blQ)wkx p-of(fk/_k00020*02 0(abs(hk)=gl1d on x_abe三 s (khzmy_abe三B岡(dby) axis(p270ulD 5ffikhz) WPH2*pi*3000八 WSH2*pi*l 2000HPHP1Q5SH60 八 2 一 Wp2llcheb2od(wpwsdpbs

9、s-)八 Lb262llcheby2(n2Lpwp2)八 subp-0f(223八 fkHpk200552=200pwkH2*pi4;?T h kufreqs(b2 62 一 wk)八 p-of(fk/_k00020*02 0(abs(hk)=gl1d on x_abe三 s (khzmy_abe三B岡(dby) axis(p270ulD 裁青一 WPH2*pi*3000八 WSH2*pi*l 2000HPHP1Q5SH60 八 三 wp lHchebod(wpws IP bss_)八 al blllchebwmmwpk-T-OT- subp-0f(223八 fkHpk200552=200pw

10、kH2*pi4;?T hkHfeqs(blQ)wkx p-of(fk/_k00020*02 0(abs(hk)=gl1d on x_abe三 s (khzmy_abe三B岡(dby) axis(p270ulD sf(khz) WPH2*pi*3000八 WSH2*pi*l 2000HPHP1Q5SH60 八 rn一 wpoHe三pod(wpWSJPbss_)八 rbbHe三P(PP60WPP-S) subp-0f(223八 fkHpk200552=200pwkH2*pi4;?T hkHfreqs(bawk)- p-of(fk/_k00020*02 0(abs(hk)=gl1d on x_abe三 s (khzmy_abe三B岡(dby) axis(p270ulD 菌(52) WPH2*P 二200p9000rwsH2*p 二400P70001HPA6SH20 八 Lnpweplle三pod(wpwsdpbss-)八 bs

溫馨提示

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

評論

0/150

提交評論