山东建筑大学语音信号处理实验四实验报告.doc

山东建筑大学实验报告学院:信息与电气工程学院 班级: 电信111 姓名: 学号: 课程: 语音信号处理 实验日期: 2014 年 4 月 23 日 成绩: 实验四 减谱法语音增强技术研究一、实验目的本实验要求掌握减谱法语音增强的原理，会利用已学的相关语音特征，构建语音特征矢量，然后自己设计减谱法语音增强程序（也可参考相关文献），能显示干净语音和加噪语音信号及处理后的结果语音信号波形，分析实验结果，写出报告。二、实验原理谱相减方法是基于人的感觉特性，即语音信号的短时幅度比短时相位更容易对人的听觉系统产生影响，从而对语音短时幅度谱进行估计，适用于受加性噪声污染的语音。由于语音是短时平稳的，所以在短时谱幅度估计中认为它是平稳随机信号，假设、和分别代表语音、噪声和带噪语音，、和分别表示其短时谱。假设噪声是与语音不相关的加性噪声。于是得到信号的加性模型： (4-1)对功率谱有 (4-4)原始语音的估值为 (4-5)只要在频域用(4-5)式得到纯净语音的谱估计，就可以根据(4-6)式得到增强后的语音。 (4-6)根据前面分析，我们可以给出谱相减算法的整个算法流程，如图4-1所示：增强后的语音带噪语音啊FFTIFFT相位信息减去估计噪声图4-1 谱减法的算法流程三、实验程序及结果分析1、不同强度的噪声与纯净信号叠加后的带噪信号及它们的频谱比较%在噪声环境下语音信号的增强%语音信号为读入的声音文件%噪声为正态随机噪声clear;input=wavread(H:soundb1.wav);count=length(input);noise1=0.1*randn(1,count);signal=input;for i=1:countvoice1(i)=signal(i)+noise1(i);endnoise2=0.01*randn(1,count);for i=1:countvoice2(i)=signal(i)+noise2(i);endnoise3=randn(1,count);signal=input;for i=1:countvoice3(i)=signal(i)+noise3(i);endn=1:count;figure %对比纯净语音信号，噪音信号和带噪语音信号subplot(3,1,1);plot(n,signal);title(纯净信号)subplot(3,1,2);plot(n,noise1);title(噪音信号)subplot(3,1,3);plot(n,voice1);title(带噪信号)figure %对比纯净语音信号频谱，噪音信号和带噪信号频谱Fss=fft(signal);subplot(3,1,1);plot(n,abs(Fss);title(纯净信号频谱)Fss1=fft(noise1);subplot(3,1,2);plot(n,abs(Fss1);title(噪音信号频谱)Fv1=fft(voice1);subplot(3,1,3)plot(n,abs(Fv1);title(带噪信号的频谱)figure %对比纯净语音信号，噪音信号和带噪语音信号subplot(3,1,1);plot(n,signal);title(纯净信号)subplot(3,1,2);plot(n,noise2);title(噪音信号)subplot(3,1,3);plot(n,voice2);title(带噪信号)figure %对比纯净语音信号频谱，噪音信号和带噪信号频谱Fss=fft(signal);subplot(3,1,1);plot(n,abs(Fss);title(纯净信号频谱)Fss2=fft(noise2);subplot(3,1,2);plot(n,abs(Fss2);title(噪音信号频谱)Fv2=fft(voice2);subplot(3,1,3)plot(n,abs(Fv2);title(带噪信号的频谱)figure %对比纯净语音信号，噪音信号和带噪语音信号subplot(3,1,1);plot(n,signal);title(纯净信号)subplot(3,1,2);plot(n,noise3);title(噪音信号)subplot(3,1,3);plot(n,voice3);title(带噪信号)figure %对比纯净语音信号频谱，噪音信号和带噪信号频谱Fss=fft(signal);subplot(3,1,1);plot(n,abs(Fss);title(纯净信号频谱)Fss3=fft(noise3);subplot(3,1,2);plot(n,abs(Fss3);title(噪音信号频谱)Fv3=fft(voice3);subplot(3,1,3)plot(n,abs(Fv3);title(带噪信号的频谱)实验结果：（1）noise1=0.1*randn(1,count);（2）noise2=0.01*randn(1,count);（3）noise3=randn(1,count);分析以上三种情况可知，不同强度的噪声叠加到纯净信号上时，对纯净信号的干扰是不同的。噪声加强，会使原始信号被噪声淹没；噪声减弱，对原始信号的影响较小。2、 用减谱法给语音信号降噪 %在噪声环境下语音信号的增强%语音信号为读入的声音文件%噪声为正态随机噪声clear;input=wavread(H:soundb1.wav);count=length(input);noise=0.1*randn(1,count);signal=input;for i=1:countvoice(i)=signal(i)+noise(i);endFv=fft(voice);anglev=angle(Fv);Fn=fft(noise);power1=(abs(Fv).2;power2=(abs(Fn).2;power3=power1-power2;power4=sqrt(power3);Fs=power4.*exp(j*anglev);sound=ifft(Fs);n=1:count;%纯净语音信号频谱Fss=fft(signal);figure %对比纯净语音信号和输出信号subplot(3,1,1);plot(n,signal);title(纯净信号)subplot(3,1,2);plot(n,noise);title(噪声信号)subplot(3,1,3)plot(n,sound);title(输出信号)figure %对比纯净语音信号频谱和输出语音信号频谱subplot(3,1,1)plot(n,abs(Fss);title(纯净信号频谱)subplot(3,1,2)plot(n,abs(Fn);title(噪声信号频谱)subplot(3,1,3)plot(n,abs(Fs);title(输出信号频谱)max_v=max(voice); %对带噪信号抽样值点进行归一化处理re_voice=voice/max_v;%对输出信号抽样点值进行归一化处理max_s=max(sound);re_sound=sound/max_s;%读出带噪语音信号，存为1001.wavwavwrite(re_voice,5500,16,1001);%读出处理后语音信号，存为1002.wavwavwrite(re_sound,5500,16,1002)实验结果：噪声为0.1*randn(1,coun)纯净信号为输入信号分析上面的实验结果可知，用减谱法対带噪信号进行降噪处理后，信号质量明显提升，噪声的干扰减小，从而得到了较纯净的语音信号。3、用改进的减谱法给语音信号降噪当噪声为0.1*randn(1,coun) , 纯净信号为输入信号时，程序如下：%在噪声环境下语音信号的增强%语音信号为读入的声音文件%噪声为正态随机噪声clear;input=wavread(H:soundb1.wav);count=length(input);noise=0.1*randn(1,count);signal=input;for i=1:countvoice(i)=signal(i)+noise(i);endFv=fft(voice);anglev=angle(Fv);Fn=fft(noise);power1=(abs(Fv).2;power2=(abs(Fn).2;a=3;b=0.01;for i=1:count if(power1(i)=a*power2(i) power3(i)=power1(i)-a*power2(i); else power3(i)=b*power2(i); endendpower4=sqrt(power3);Fs=power4.*exp(j*anglev);sound=ifft(Fs);n=1:count;%纯净语音信号频谱Fss=fft(signal);figure %对比纯净语音信号和输出信号subplot(2,1,1)plot(n,signal);title(纯净信号)subplot(2,1,2)plot(n,sound);title(输出信号)figure %对比纯净语音信号频谱和输出语音信号频谱subplot(2,1,1)plot(n,abs(Fss);title(纯净信号频谱)subplot(2,1,2)plot(n,abs(Fs);title(输出信号频谱)实验结果：（1）a=3，b=0.01时：（2） a=9,b=0.01时：（3） a=0.5,b=0.01时：（4） a=3,b=0.001时：（5）a=3,b=0.05(6) a=3,b=0.1时：分析以上取不同的a、b值得到的实验结果可知，用改进的减谱法対带噪信号进行语音增强处理后，语音信号质量进一步提高，信号更加地纯净。参数a、b的选择对结果影响很大，因为本实验所用的是辅音，所以选择a=3,b=0.01左右较合适。a和b太大或太小，都不能取得较好的降噪和抑制纯音噪声的效果。4、 实验总结通过本次实验，我掌握了减谱法语音增强的原理，并能够自己设计减谱法语音增强程序。减谱法能够有效地降低包含在语音中