基于matlab的去除工频干扰巴特沃斯滤波器

程序及波形：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);xlabel('时间/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('信号滤波后频域图');