巴特沃兹带阻滤波器.doc

下载文档 降价啦

7
0
约2.1千字
约 5页
2020-07-02 发布于湖北
举报
版权申诉
保障服务

巴特沃兹带阻滤波器.doc

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

clf; ft=2;fp1=0.35;fp2=0.55;fs1=0.4;fs2=0.5; wp1=2*pi*fp1/ft;wp2=2*pi*fp2/ft;ws1=2*pi*fs1/ft;ws2=2*pi*fs2/ft; omegp1=tan(wp1/2);omegp2=tan(wp2/2);omegs1=tan(ws1/2);omegs2=tan(ws2/2); BW=omegs2-omegs1;omeg0=omegs1*omegs2;omeg00=sqrt(omeg0); omegs=1;omegp=omegp1*BW/(omeg0-omegp1*omegp1); [n,wn]=buttord(omegs, omegp, 1, 40,'s'); [b,a] = butter(n, wn, 's'); [num1,den1]=lp2bs(b,a,omeg00,BW); [num,den]=bilinear(num1,den1,0.5) [z,p,k]=tf2zp(num,den); disp('分子系数是');disp(num); disp('分母系数是');disp(den); [h,w]=freqz(num,den,512); g = 20*log10(abs(h)); figure(1) plot(w/pi,g);grid axis([0 1 -200 1]); xlabel('频率，Hz'); ylabel('增益， dB'); title('巴特沃兹带阻滤波器'); figure(2) zplane(z,p); xlabel('z'); ylabel('p'); title('传输函数零极点'); k=0:1:1001; f1=2*pi*450/2000; f2=2*pi*800/2000; x1=sin(f1*k);%正弦波信号 figure(3); subplot(2,2,1) plot(x1);grid on;%500Hz正弦波波形 axis([0,100*pi,-5,5]); xlabel('t(s)'); ylabel('x1'); title('x1正弦波信号'); x2=sin(f2*k);%正弦波信号 subplot(2,2,2) plot(x2);grid on;%500Hz正弦波波形 axis([0,100*pi,-5,5]); xlabel('t(s)'); ylabel('x2'); title('x2正弦波信号'); x=sin(f1*k)+sin(f2*k);%正弦波信号 subplot(2,2,3) plot(x);grid on;%500Hz正弦波波形 axis([0,100*pi,-5,5]); xlabel('t(s)'); ylabel('x'); title('x正弦波信号'); y=filter(num,den,x) ; subplot(2,2,4);%数字滤波函数输出 plot(y) ;grid on;%数字滤波器输出波形 axis([0,100*pi,-5,5]); xlabel('f(hz)'); ylabel('y'); title('数字滤波器输出波形 ');