基于matlab的去除工频干扰巴特沃斯滤波器(共3页).doc

《基于matlab的去除工频干扰巴特沃斯滤波器(共3页).doc》由会员分享，可在线阅读，更多相关《基于matlab的去除工频干扰巴特沃斯滤波器(共3页).doc（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上程序及波形：f_N=200;f_p=30,70;f_s=49,51;R_p=3;R_s=20;Ws=f_s/(f_N/2);Wp=f_p/(f_N/2);n,Wn=buttord(Wp,Ws,R_p,R_s);b,a=butter(n,Wn,stop);f1=30;f2=50;f3=70;fs=200;figure,freqz(b,a,512,200);title(带阻滤波器)t=0:1/fs:1;s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t); figure,subplot(2, 1, 1),plot(t,s);xlab

2、el(时间/s);ylabel(幅度);title(信号滤波前时域图);Fs=fft(s,512);AFs=abs(Fs);f=(0:255)*fs/512;figure,subplot(2, 1, 2),plot(f,AFs(1:256);xlabel(频率/Hz);ylabel(幅度);title(信号滤波前频域图);sf=filter(b,a,s);figure,subplot(2, 1, 1),plot(t,sf);xlabel(时间/s);ylabel(幅度);title(信号滤波后时域图);subplot(2, 1, 2)Fsf=fft(sf,512);AFsf=abs(Fsf);f=(0:255)*fs/512;plot(f,AFsf(1:256);xlabel(频率/Hz);ylabel(幅度);title(信号滤波后频域图);专心-专注-专业