分帧加窗处理增强

分帧加窗处理增强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(純语音波形); grid minor;subplot(212)plot(x);axis(0 4*104 -1 1)title(混噪后的语音波形); grid minor;%-进行分帧，加窗处理framelength=256; %窗长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;% 求帧数，size返回矩阵的行数和列数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(21