fp=1500;fc=1700;As=100;Ap=1; %(以上为低通滤波器的性能指标)
wc=2*pi*fc/fs;wp=2*pi*fp/fs;
wdel=wc-wp;
beta=0.112*(As-8.7);
N=ceil((As-8)/2.285/wdel);
wn=kaiser(N+1,beta);
ws=(wp+wc)/2/pi;
b=fir1(N,ws,wn);figure(3);
freqz(b,1);
%(此前为低通滤波器设计阶段)——接下来为去除噪声信号的程序 x=fftfilt(b,y,z);
X=fft(x,n);
figure(4);
subplot(2,2,1);
plot(f,abs(y_zp(1:n/2)));
axis([0 20000 0 2000]);
title('滤波前信号的频谱');
subplot(2,2,2);
plot(f,abs(X(1:n/2)));
axis([0 20000 0 2000]);
title('滤波后信号频谱');
subplot(2,2,3);
plot(y_z);
title('滤波前信号的波形')
axis([0 200000 -15 15]);
subplot(2,2,4);
plot(x);
axis([0 200000 -1 1]);
title('滤波后信号的波形')
Magnitude (dB)
评论2