pujianfa.m

基于谱减法的语音增强算法

clear all;  
close all; 
if nargin<3 
    nbits=8; 
end 
if nargin<2 
    fs=8000; 
end 
if nargin<1 
    [y,fs,nbits]=wavread('Pool.wav'); 
  end 
 
%------进行分帧，加窗处理 
framelength=256; 
step=80; 
y_frame=enframe(y,framelength,step); 
nf=size(y_frame,1);              % 求帧数 
w=hamming(framelength); 
y_window=[]; l.
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('处理后的波形');