[x1, fs,]=audioread('quzao.wav');
n=length(x1);
y_p=fft(x1,n);
f=fs*(0:n/2-1)/n;
figure(1);
subplot(2,1,1);
plot(x1);xlabel('t');ylabel('y');title(' 时域波形 ');
subplot(2,1,2);
plot(f,abs(y_p(1:n/2)));
title('原始语音信号采样后频谱图');
xlabel('频率Hz');
L=length(x1);
noise=0.1randn(L,1);
y_z=x1+noise;
n=length(x1);
y_zp=fft(y_z,n);
f=fs(0:n/2-1)/n;
figure(2)
subplot(2,1,1);
plot(y_z);
title('加噪语音信号时域波形');
xlabel('时间轴')
ylabel('幅值 A')
subplot(2,1,2);
plot(f,abs(y_zp(1:n/2)));
title('加噪语音信号频谱图');
xlabel('频率Hz');
ylabel('频率幅值');
fp=1000;
fr=1200;
ap=1;
as=40;
[n,fn]=cheb1ord(fp/(fs/2),fr/(fs/2),ap,as,'z');
[b,a]=cheby1(n,ap,fn);
freqz(b,a);
[h,w]=freqz(b,a);
subplot(221);plot(wfs/(2pi),abs(h));title('cheby1低通滤波');
signal=filter(b,a,y_z);
X=fft(signal,n);
f=fs*(0:n/2-1)/n;
figure(3);
subplot(2,1,1);
plot(f,abs(y_zp(1:n/2)));
title('滤波前信号的频谱');
subplot(2,1,2);
plot(X);
title('滤波后信号频谱');
figure(4);
subplot(2,1,1);
plot(y_z);
title('滤波前信号的波形')