matlab滤波器代码要点.docxVIP

下载本文档

8
0
约4.69千字
约 8页
2021-04-14 发布于天津
举报
版权申诉

matlab滤波器代码要点.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

数字信号处理：已知通带截止频率fp=5kHz，通带最大衰减ap=2dB,阻带截止频率fs=2kHz，阻带最小衰减as=30dB,按照以上技术指标设计巴特沃斯低通滤波器： 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( 频率( 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));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 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 三 po「d(wpws IP Is)八 pAHe三p(Ndpswc= Lhk-wkllfreqz(B- fkuwk、pi*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 切比雪夫低通滤波器频率(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 =nlaudipfs-bifslrwawead(audiofi-e= 〔bbHchebylpp103)八 ouflaudioHfi_(D「(bbjnlaudio)八 sound(ouflaudipfskfs)八 Wavwme(ouflaudipf3bifs=esflour= xk