分帧加窗处理增强.docVIP

分帧加窗处理增强 hc wavread D:\matlab安装程序\work\1.wav ; fs 22050; t 0:1/22050: length hc -1 /22050; h 0.1*randn size t ; s hc; x s+h; figure 1 subplot 211 plot hc ; title 純语音波形 ; grid minor; subplot 212 plot x ; axis [0 4*10^4 -1 1] title 混噪后的语音波形 ; grid minor; % framemove 80; %帧移 y_frame0 enframe hc,framelength,framemove ; y_frame enframe x,framelength,framemove ; figure 2 subplot 211 plot y_frame0 ; title 对纯语音分帧看结果 ; grid minor; subplot 212 plot y_frame ; axis [0 500 -1 1] title 对混噪后的语音分帧看结果 ; grid minor; % 求帧数nf size y_frame,1 ; w hamming framelength ; y_window []; %忽略第一个 for k 1:nf y_row y_frame k,: .*w ; % 加汉明窗 y_window [y_window;y_row] ; end figure 3 subplot 211 plot w ; title 窗函数 ; grid minor; subplot 212 plot y_window ; axis [0 300 0 1] title 对每帧加窗看结果 ; grid minor; %----谱减系数设定 及噪声谱估计a 2; b 1; % 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 ; % 噪声的谱估计figure 4 plot noise_amplitude title 前十帧噪声估计的幅度谱 ; grid minor; %-----基本谱减法谱减过程y_fft fft y_window ; y_amplitude abs y_fft ; %fft后y的幅度 y_angle angle y_fft ; % fft后y的相位 figure 5 subplot 211 plot y_amplitude ; title 窗的幅度 ; grid minor; subplot 212 plot y_angle ; title 窗的相位 ; grid minor; 谱减过程 y_enhanced_amplitude []; for k 1:nf y_temp y_amplitude :,k .^a-b*noise_amplitude.^a; % 谱相减y_temp0 y_temp.^ 1/a ; y_enhanced_amplitude [y_enhanced_amplitude y_temp];%得到处理后的幅度矩阵end figure 6 subplot 211 plot y_temp0 ; % title 普通谱减法处理后语音信号幅度的平方 ; grid minor; subplot 212 plot y_enhanced_amplitude title 普通谱减法处理后语音信号幅度 ; grid minor; 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,: ]; figure 7 subplot 211 plot hc ; title 純语音波形 ; grid minor; subplot 212 plot output ; title 处理后的波形 ;%-----做出处理前后图形，给出直观对比 grid minor; sound hc,fs s