f1=filter(bz,az,z); %利用上面已做好的滤波器(已bz和az为特征的
滤波器)对x2信号进行滤波,赋值给f1 figure(10)
subplot(2,1,1) %创建两行一列绘图区间的第1绘图区间
plot(t,z) %做加噪后的信号时域图形
title('滤波前的时域波形');
subplot(2,1,2) %创建两行一列绘图区间的第2绘图区间
plot(t,f1);
title('滤波后的时域波形');
sound(f1,22050); %播放语音信号
p=length(f1);
F0=fft(f1,p);
f=0:fs/p:fs*(s-1)/p;
figure(11)
y2=fft(z,1024); %对信号y2做1024点FFT变换
subplot(2,1,1); %创建两行一列绘图区间的第1绘图区间
plot(f,abs(y4));
title('滤波前的频谱')
xlabel('Hz'); %x轴的名字是“Hz”
ylabel('fuzhi'); %y轴的名字是“fuzhi”
subplot(2,1,2)
F1=plot(f,abs(F0));
title('滤波后的频谱')
xlabel('Hz');
ylabel('fuzhi');