用matlab设计一个带阻滤波器,信号有5Hz,40Hz,50Hz,要求把40Hz的滤除,程序是这样的,滤波器和产生的信

用matlab设计一个带阻滤波器,信号有5Hz,40Hz,50Hz,要求把40Hz的滤除,程序是这样的,滤波器和产生的信w1=38/100;w2=42/100;[B,A]=butter(1,[w1,w2],'stop');[h,w]=freqz(B,A);f=w/pi*100;figure(1),plot(f,20*log10(abs(h)));axis([0,80,-30,10]);grid;xlabel('频率/Hz')ylabel('幅度/dB')title('频率响应曲线')fs=200;t=0:1/fs:1;x1=sin(2*pi*5*t)+sin(2*pi*50*t)+sin(2*pi*40*t);Y1=fft(x1,128);Y2=abs(Y1);f=(0:64)*100/64;figure(2),subplot(2,1,1),plot(f,Y2(1:65)),title('原始信号');y=filter(B,A,Y2);subplot(2,1,2),plot(y),title('滤波后信号');