matlab滤波器代码要点.pdfVIP

下载本文档

1
0
约5.66千字
约 8页
2021-10-14 发布于上海
举报
版权申诉

matlab滤波器代码要点.pdf

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]) p195-14 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,-100,0]) xlabel( 频率 (kHZ));ylabel( 幅度(db));grid on P195-16 wp=2*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 on P195-15 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 切比雪夫低通滤波器 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