分帧加窗处理增强.doc

《分帧加窗处理增强.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、分帧加窗处理增强hc=wavread('D:matlab安装程序work1.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('純语音波形');gridminor;subplot(212)plot(x);axis([04*10^4-11])title('混噪后的语音波形');gridminor;%------进行分帧，加窗处理framelength=256;%窗长framemove=80;%帧移y

2、_frame0=enframe(hc,framelength,framemove);y_frame=enframe(x,framelength,framemove);figure(2)subplot(211)plot(y_frame0);title('对纯语音分帧看结果');gridminor;subplot(212)plot(y_frame);axis([0500-11])title('对混噪后的语音分帧看结果');gridminor;%求帧数，size返回矩阵的行数和列数nf=size(y_frame,1);w=hamming(framelength)';y_window=[];

3、%忽略第一个fork=1:nfy_row=y_frame(k,:).*w;%加汉明窗y_window=[y_window;y_row];endfigure(3)subplot(211)plot(w);title('窗函数');gridminor;subplot(212)plot(y_window);axis([030001])title('对每帧加窗看结果');gridminor;%----谱减系数设定及噪声谱估计a=2;b=1;%2<=a<=4,1<=b<=7是系数的可调范围a=2,b=1时是基本谱减法n_noise=10;%取前十帧语音作为背景噪声y_noise=mean(y_w

4、indow(1:10,:));noise_amplitude=abs(fft(y_noise))';%噪声的谱估计figure(4)plot(noise_amplitude)title('前十帧噪声估计的幅度谱');gridminor;%-----基本谱减法谱减过程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('窗的幅度');gridminor;subplot(212)plot

5、(y_angle);title('窗的相位');gridminor;谱减过程y_enhanced_amplitude=[];fork=1:nfy_temp=y_amplitude(:,k).^a-b*noise_amplitude.^a;%谱相减y_temp0=y_temp.^(1/a);y_enhanced_amplitude=[y_enhanced_amplitudey_temp];%得到处理后的幅度矩阵endfigure(6)subplot(211)plot(y_temp0);%??????????title('普通谱减法处理后语音信号幅度的平方');gridminor;sub

6、plot(212)plot(y_enhanced_amplitude)title('普通谱减法处理后语音信号幅度');gridminor;y_enhanced_fft=y_enhanced_amplitude.*exp(i*y_angle);%恢复原来的相位y_enhanced=real(ifft(y_enhanced_fft)');%fft反变换求得增强后的分帧语音，取的是实部output=[];fork=1:nf-1output=[outputy_enhanced(k,1:80)];%从增强后的分帧语音中抽取出最终的语音序列endoutput=[outputy_enhanced(

7、k,:)];figure(7)subplot(211)plot(hc);title('純语音波形');gridminor;subplot(212)plot(output);title('处理后的波形');%-----做出处理前后图形，给出直观对比gridminor;sound(hc,fs)sound(x,fs)sound(output,fs)