matlab-denoise-of-voice.rar_denoise_signal denoise_speech denois

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频谱图')