clear
clc
[x1,fs,bits]=wavread('xinhao.wav') ; %x为数据点，fs为采样频率，bits为位数
sound(x1,fs,bits) %播放原始语音
pause(3);%暂停三秒，等待语音播放完毕
f=fs*linspace(0,1,size(x1,1));
figure(1)
y1=fft(x1,size(x1,1));%y1为x的频谱
subplot(2,1,1)
plot(x1);
title('原始语音信号')
subplot(2,1,2)
plot(f,abs(y1));
title('语音fft频谱图 ')
%加入噪声
t=0:1/fs:(size(x1,1)-1)/fs; %将所加噪声信号的点数调整到与原始信号相同
Au=0.5;
d=[Au*cos(2*pi*10000*t)+Au*cos(3*pi*10000*t)]'; %噪声为10khz与15khz的余弦
x2(:,1)=x1(:,1)+d;
x2(:,2)=x1(:,2)+d;
sound(x2,fs,bits) ; %播放加噪声后的语音
pause(3);%暂停三秒，等待语音播放完毕
figure(2)
y2=fft(x2,size(x1,1)); %加噪声后的fft频谱图
subplot(2,1,1)
plot(t,x2);
title('加噪声语音时域图');
subplot(2,1,2)
plot(f,abs(y2));
title('加噪语音fft频域图');
%设计切比雪夫滤波器
[N,Wn]=cheb1ord(9000*2/fs,9800*2/fs,1,40,'s');
[B,A]=cheby1(N,1,Wn);
figure(3)
freqz(B,A)
title('滤波器频谱图')
x3=filter(B,A,x2); %x3为滤波完的时域数据
sound(x3,fs,bits) ; %播放滤波后的语音
pause(3);%暂停三秒，等待语音播放完毕
y3=fft(x3,size(x1,1));
figure(4)
subplot(2,1,1)
plot(t,x3)
title('滤波后的语音时域图')
subplot(2,1,2)
plot(f,abs(y3))
title('滤波后的语音fft频谱图')
