testthree.m

提供了集中语音去噪的matlab代码

% if nargin<3 
% nbits=8; 
% end 
% if nargin<2 
% fs=8000; 
% end 
% if nargin<1 
% [filename,pathname]=uigetfile('speech_clean.wav'); 
% [y,fs,nbits]=wavread([pathname,filename]); 
% end 
[filename,pathname]=uigetfile('*.wav','请选择语音文件：');
[y,fs,nbits]=wavread([pathname filename]);
%------进行分帧，加窗处理 
framelength=256; 
step=80; 
y_frame=enframe(y,framelength,step); 
% [m,n]=size(y_frame); % 求帧数 
nf=347;
w=hamming(framelength)'; 
y_window=[]; 
for k=1:nf 
y_row=y_frame(k,:).*w; % 加汉明窗 
y_window=[y_window;y_row]; 
end 

%----谱减系数设定 及噪声谱估计 
a=2; 
b=4; % 2<=a<=4,1<=b<=7 是系数的可调范围 a=2,b=1时是基本谱减发 

n_noise=10; % 取前十帧语音作为背景噪声 
y_noise=mean(y_window(1:10,:)); 
noise_amplitude=abs(fft(y_noise))'; % 噪声的谱估计 

%-----基本谱减法谱减过程 
y_fft=fft(y_window'); 
y_amplitude=abs(y_fft); % ff后y的幅度 
y_angle=angle(y_fft); % ff后y的相位 

y_enhanced_amplitude=[]; 
for k=1:nf 
y_temp=y_amplitude(:,k).^a-b*noise_amplitude.^a; % 谱相减 
y_temp=y_temp.^(1/a); 
y_enhanced_amplitude=[y_enhanced_amplitude y_temp];%得到处理后的幅度矩阵 
end 
y_enhanced_fft=y_enhanced_amplitude.*exp(i*y_angle); % 恢复原来的相位 
y_enhanced=real(ifft(y_enhanced_fft)'); % fft反变换求得增强后的分帧语音，取的是实部 
output=[]; 
for k=1:nf-1 
output=[output y_enhanced(k,1:80)]; % 从增强后的分帧语音中抽取出最终的语音序列 
end 
output=[output y_enhanced(k,:)]; 


%-----做出处理前后图形，给出直观对比 
subplot(2,1,1) 
plot(y); 
title('处理前的波形'); 


subplot(2,1,2) 
plot(output); 
title('处理后的波形');