各种滤波器程序.doc

下载文档

150
0
约3.9千字
约 8页
2017-12-12 发布于河南
举报
版权申诉
保障服务

各种滤波器程序.doc

1、本文档共8页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

各种滤波器程序

巴特沃斯低通滤波器程序代码： wp=2*pi*5000;ws=2*pi*12000;Rp=2;As=30;%设置滤波器参数 [N,wc]=buttord(wp,ws,Rp,As,s);%计算滤波器阶数N和3dB截止频率wc [B,A]=butter(N,wc,s);%计算滤波器系统函数分子和分母多项式系数 k=0:511;fk=0:14000/512:14000;wk=2*pi*fk; Hk=freqs(B,A,wk);%计算滤波器的频率响应 plot(fk/1000,20*log10(abs(Hk))); grid on;%以dB表示频率响应的幅度 title(巴特沃斯低通滤波器); xlabel(频率(kHz)); ylabel(幅度(dB)); axis([0,14,-40,5]): 运行结果为： N=5，wc=3.7792e+004，B=7.7094e+022 A=[1 1.223e+005 7.4785e+009 2.8263e+014 6.6014e+018 7.7094e+022] 巴特沃斯高通滤波器代码 %调用函数buttord和butter直接设计巴特沃斯高通滤波器 wp=2*pi*4000;ws=2*pi*1000;Rp=0.1;As=30;%设置高通滤波器参数 [N,wc]=buttord(wp,ws,Rp,As,s);%计算滤波器阶数N和3dB截止频率wc [BH,AH]=butter(N,wc,high,s);%计算高通滤波器系统函数分子分母多项式系数 k=0:511;fk=0:14000/512:14000;wk=2*pi*fk; Hk=freqs(BH,AH,wk); plot(fk/1000,20*log10(abs(Hk)));grid on; title(巴特沃斯高通滤波器); xlabel(频率(kHz));ylabel(幅度(dB)); axis([0,14,-40,5]); 运行结果为： N=4，BH=[1 0 0 0 0] AH=[1 38930 7.5778e+008 8.6405e+012 4.9261e+016 ] 巴特沃斯带通滤波器程序代码 wp=2*pi*[4000,7000];ws=2*pi*[2000,9000];Rp=1;As=20;%设置带通滤波器%参数 [N,wc]=buttord(wp,ws,Rp,As,s);%计算带通滤波器阶数N和3dB截止频率%wc [BB,AB]=butter(N,wc,s);%计算带通滤波器系统函数分子和分母多项式系%数向量BB和AB k=0:511;fk=0:14000/512:14000;wk=2*pi*fk; Hk=freqs(BB,AB,wk); plot(fk/1000,20*log10(abs(Hk)));grid on title(巴特沃斯带通滤波器); xlabel(频率(kHz));ylabel(幅度(dB)); axis([0,14,-40,5]) 运行结果： N=5;BB=[0 0 0 0 0 6.9703e+021 0 0 0 0 0] AB=[1 75625 8.3866e+009 4.0121e+014 2.2667e+019 7.0915e+023 2.5056e+028 4.9024e+032 1.1328e+037 1.1291e+041 1.6504e+045] 巴特沃斯带阻滤波器 wp=2*pi*[2000,9000];ws=2*pi*[4000,7000];Rp=1;As=20;%设置带阻滤波器%参数 [N,wc]=buttord(wp,ws,Rp,As,s);%计算带阻滤波器阶数N和3dB截止频率%wc [BSb,ASb]=butter(N,wc,stop,s);%计算带阻滤波器系统函数分子和分母多%项式系数 k=0:511;fk=0:14000/512:14000;wk=2*pi*fk; Hk=freqs(BSb,ASb,wk); plot(fk/1000,20*log10(abs(Hk)));grid on title(巴特沃斯带阻滤波器); xlabel(频率(kHz));ylabel(幅度(dB)); axis([0,14,-40,5]); 运行结果为： N=5, BSb=[1 0 5.527e+009 0 1.2219e+019 0 1.3507e+028 0 7.4652e+036 0 1.6504e+045] ASb=[1 96579 1.0191e+010 5.6622e+014 3.0252e+019 1.0394e+024 3.3441e+028 6.9186e+032 1.3764e+037 1.442e+041 1.