基于Matlab编写的语音端点检测

基于Matlab编写的语音端点检测专业：班级：姓名：指导教师：2011 年6月18 日一、实验目的1．学会MATLAB的使用，掌握MATLAB的程序设计方法；3．掌握语音处理的基本概念、基本理论和基本方法；4．掌握基于MATLAB编程实现带噪语音信号端点检测；5．学会用MATLAB对信号进行分析和处理。

二、实验内容简介：可利用时域分析（短时能量、短时过零率、短时自相关）方法的某一个特征或某几个特征的结合，判定某一语音信号的端点，尤其在有噪声干扰时，如何准确检测语音信号的端点，这在语音处理中是富有挑战性的一个课题。

要求：（1）录制语音，读入文件，绘制波形（2）分帧，绘制能量曲线和短时过零率曲线（3）根据上述端点检测原理，实现端点检测（4）界面三．课程设计原理端点检测是语音信号处理过程中非常重要的一步，它的准确性直接影响到语音信号处理的速度和结果，因此端点检测方法的研究一直是语音信号处理中的热点。

本设计使用传统的短时能量和过零率相结合的语音端点检测算法利用短时过零率来检测清音．用短时能量来检测浊音，两者相配合便实现了信号信噪比较大情况下的端点检测。

算法对于输入信号的检测过程可分为短时能量检测和短时过零率检测两个部分。

算法以短时能量检测为主，短时过零率检测为辅。

根据语音的统计特性，可以把语音段分为清音、浊音以及静音（包括背景噪声）三种。

在本算法中，短时能量检测可以较好地区分出浊音和静音。

对于清音，由于其能量较小，在短时能量检测中会因为低于能量门限而被误判为静音；短时过零率则可以从语音中区分出静音和清音。

将两种检测结合起来，就可以检测出语音段（清音和浊音）及静音段1、短时能量计算定义n 时刻某语音信号的短时平均能量n E 为：∑∑--=+∞-∞=-=-=n N n m m n m n w m x m n w m x E )1(22)]()([)]()([式中N 为窗长，可见短时平均能量为为一帧样点值的甲醛平方和。

特殊地，当窗函数为矩形窗时，有∑--==n N n m n m x E )1(2)(2、短时过零率过零就是指信号通过零值。

过零率就是每秒内信号值通过零值的次数。

对于离散时间序列，过零则是指序列取样值改变符号，过零率则是每个样本的改变符号的次数。

对于语音信号，则是指在一帧语音中语音信号波形穿过横轴（零电平）的次数。

可以用相邻两个取样改变符号的次数来计算。

如果窗的起点是n=0，短时过零率Z 为波形穿过横轴（零电平）的次数短时过零可以看作信号频率的简单度量浊音的短时平均幅度最大，无声的短时平均幅度最小，清音的短时过零率最大，无声居中，浊音的短时过零率最小。

3、短时自相关函数1、是偶函数；2、s(n)是周期的，那么R(k)也是周期的；3、可用于基音周期估计和线性预测分析4、判断语音信号的起点和终点利用短时平均幅度和短时过零率可以判断语音信号的起点和终点。

语音端点检测方法可采用测试信号的短时能量或短时对数能量、联合过零率等特征参数，并采用双门限判定法来检测语音端点，即利用过零率检测清音，用短时能量检测浊音，两者 配合。

首先为短时能量和过零率分别确定两个门∑-=--=100))1(())((21N n w w n S Sgn n S Sgn Z ∑--=+=10)()()(k N n w w w k n s n s k R ⎩⎨⎧<-≥=0,10,1)sgn(x x x限，一个是较低的门限数值较小，对信号的变化比较敏感，很容易超过;另一个是比较高的门限，数值较大。

低门限被超过未必是语音的开始，有可能是很短的噪声引起的，高门限被超过并且接下来的自定义时间段内的语音超首先，可根据浊语音情况下短时平均幅度M的概率密度函数P(M/V)确定一个阈值参数MH, MH的值定得比较高。

当一帧输入信号的M值超过MH 时，就可以十分肯定该帧语音信号不是无声，而有相当大的可能性是浊音由N1’向前和N2’向后继续用短时过零率Z进行搜索。

根据无声情况下短时过零率Z的均值，设置一个参数Z0，如果由N1’向前搜索时Z始终大于Z0的3倍，则认为这些信号仍属于语音段，直至Z突然下降到低于3Z0值时，这时可以确定语音的精确起点（为了保证可靠，由N1’向前搜索时间不超过25ms）。

对终点做同样的处理。

采取这一算法的原因在于，N1’以前可能是一段清辅音段（如f,s），它的能量相当弱，依靠能量不可能把他们与无声段区别开，而他们的过零率明显高于无声，因而能用这个参数来精确的判断二者的分割点，也就是语音真正的起点。

四、实验步骤及分步结果1.界面设计2.程序设计（1）.①打开语音文件function openbutton_Callback(hObject, eventdata, handles) clc;axes(handles.wavaxes);cla reset;boxon;set(gca,'XTickLabel',[],'YTickLabel',[]);axes(handles.Energeaxes);cla reset;boxon;set(gca,'XTickLabel',[],'YTickLabel',[]);axes(handles.Zerorateaxes);cla reset;boxon;set(gca,'XTickLabel',[],'YTickLabel',[]);set(handles.T1edit,'string',0);set(handles.T2edit,'string',0);set(handles.T3edit,'string',0);%´ò¿ª´ý´¦ÀíµÄÓïÒôÎļþ[filename,pathname]=uigetfile({'*.wav','All WavFiles'},'Ñ¡ÔñÓïÒôÎļþ');if filename == 0return;%Èç¹ûûÓÐÑ¡ÔñеÄÎļþ£¬Ôò·µ»Øendfile=fullfile(pathname,filename);%ÎļþÃû[signal,fs,bit]=wavread(file);%¶ÁÈ¡Ñ¡ÔñµÄÓïÒôÊý¾Ýaxes(handles.wavaxes)plot(signal)%ÏÔʾ²¨ÐÎhandles.wavsignal=signal;%update handles structureguidata(hObject,handles);②预处理预处理程序function prebutton_Callback(hObject, eventdata, handles)% hObject handle to prebutton (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)signal=handles.wavsignal;signal=filter([1,-0.9375],1,signal);%Ô¤¼ÓÖØÂ˲¨axes(handles.wavaxes)plot(signal)%ÏÔʾԤ¼ÓÖغó²¨ÐÎtitle('Ô¤¼ÓÖØÂ˲¨ºóÓïÒô²¨ÐÎ');handles.wavsignal=signal;framelength=256;%´°³¤£»framenumber=fix(length(signal)/framelength);%×ÜÖ¡Êýfor i=1:framenumber;framesignal(i,1:framelength)=signal((i-1)*framelength+1:i*framelength);%·ÖÖ¡´¦Àíendhandles.framesignal=framesignal;%Update handles structureguidata(hObject,handles);③短时能量短时能量function Energebutton_Callback(hObject, eventdata, handles)% hObject handle to Energebutton (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) framesignal=handles.framesignal;framenumber=size(framesignal,1);%»ñµÃÊý¾ÝèåÊýfor i=1:framenumber;E(i)=0;E(i)=sum(framesignal(i,:).^2);%¼ÆËã¶ÌʱÄÜÁ¿endaxes(handles.Energeaxes)%c=[1 2 3 4 5 6 4 3 3 3 3];plot(E);%½«ÄÜÁ¿Ö»ÏÔʾÔÚÄÜÁ¿´°¿Úhandles.E=E;%Update handles structureguidata(hObject,handles);% --- Executes on button press in Zeroratebutton.function Zeroratebutton_Callback(hObject, eventdata, handles)% hObject handle to Zeroratebutton (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA) framesignal=handles.framesignal;framenumber=size(framesignal,1);%ÓïÒôÊý¾ÝÖ¡Êýframelength=size(framesignal,2);%»ñÈ¡ÓïÒôÖ¡³¤for i=1:framenumber;Z(i)=0;for j=2:framelength-1;Z(i)=Z(i)+abs(sign(framesignal(i,j))-sign(framesignal(i,j-1)));%¼ÆËã¶Ìʱ¹ýÁãÂÊendendaxes(handles.Zerorateaxes)%ÏÔʾ֡¹ýÁãÂÊÇúÏßplot(Z)handles.Z=Z;handles.framelength=framelength;%Update handles structureguidata(hObject,handles);④短时过零率短时过零率function Zeroratebutton_Callback(hObject, eventdata, handles)% hObject handle to Zeroratebutton (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA) framesignal=handles.framesignal;framenumber=size(framesignal,1);%ÓïÒôÊý¾ÝÖ¡Êýframelength=size(framesignal,2);%»ñÈ¡ÓïÒôÖ¡³¤for i=1:framenumber;Z(i)=0;for j=2:framelength-1;Z(i)=Z(i)+abs(sign(framesignal(i,j))-sign(framesignal(i,j-1)));%¼ÆËã¶Ìʱ¹ýÁãÂÊendendaxes(handles.Zerorateaxes)%ÏÔʾ֡¹ýÁãÂÊÇúÏßplot(Z)handles.Z=Z;handles.framelength=framelength;%Update handles structureguidata(hObject,handles);⑤设定门限1.2阈值设置门限1.2（编辑框2）function T2edit_Callback(hObject, eventdata, handles)% hObject handle to T2edit (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)clcaxes(handles.Energeaxes);cla reset;box on;set(gca,'XTicklabel',[],'YTicklabel',[]);%ÔÚÊäÈëеÄa,bֵʱ£¬ÒªÇå¿ÕEn ergeaxesÀïµÄͼÏñTL=get(handles.T1edit,'string');%»ñµÃTLÖµTH=get(handles.T2edit,'string');%»ñµÃTHÖµTL=str2num(TL);%½«×Ö·ûÀàÐÍת»¯³ÉÊýÖµÀàÐÍTH=str2num(TH);E=handles.E;MeanE=mean(E(1:10));%ÇóÇ°10Ö¡µÄÄÜÁ¿¾ùÖµaxes(handles.Energeaxes)plot(E)hold online([1 length(E)],[TL TL],'Color','r');%ÓúìÏß»-³öTLÖµhold online([1 length(E)],[TH,TH],'Color','b');%ÓÃÀ¶Ïß»-³öTHÖµhandles.TL=TL;handles.TH=TH;%Update handles structureguidata(hObject,handles);⑥编辑框3门限3function T3edit_Callback(hObject, eventdata, handles)% hObject handle to T3edit (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)clcaxes(handles.Zerorateaxes);cla reset;box on;set(gca,'XTickLabel',[],'YtickLabel',[]);v=get(handles.T3edit,'string');v=1*str2num(v);Z=handles.Z;MeanZ=mean(Z(1:10));VarZ=std(Z(1:10));ZT=min(v,MeanZ+VarZ);axes(handles.Zerorateaxes)plot(Z)hold online([1,length(Z)],[ZT,ZT],'Color','b'); handles.ZT=ZT;guidata(hObject,handles);⑦端点检测端点检测function Detectpointbutton_Callback(hObject, eventdata, handles) TL=handles.TL;%¼ÓÔØÊý¾ÝTH=handles.TH;ZT=handles.ZT;Z=handles.Z;E=handles.E;startflag=1;startframe=11;mins=min(handles.wavsignal);maxs=max(handles.wavsignal);axes(handles.wavaxes);plot(handles.wavsignal);accustartflag=1;accustartframe=startframe-15;i=accustartframe;k=0;i=11;while(startflag)%ÓöÌʱÄÜÁ¿³õ²½¼ÆËãÆðʼµãif((E(i)>TL)&(E(i+1)>TH))startframe=i;startflag=0;%±ê־λÇåÁãÍ˳öÑ-»·while(accustartflag)%ÓöÌʱ¹ýÁãÂʽøÒ»²½¼ÆËãÆðʼµãif(Z(i)>ZT);k=k+1;if(k==3)accustartflag=0;accustartframe=i-3startframe=accustartframe;endelsek=0;endi=i+1;endelsei=i+1;endendendflag=1;endframe=length(E);i=endframe;accuendflag=1;accuendframe=endframe+15;i=accuendframe;k=0;j=size(E);for i=11:j-1;%ÓɶÌʱÄÜÁ¿´ÖÅÐÖÕÖ¹µãif(i==1)break;endif((E(i)>TL)&(E(i-1)>TH))endframe=i;while(accuendflag)%ÓɶÌʱ¹ýÁãÂʽøÒ»²½¼ÆËãÖÕÖ¹µãif(i>length(Z))break;endif(Z(i)>ZT)k=k+1;if(k==3)accuendflag=0;accuendframe=i;endframe=accuendframe;endflag=0;endelsek=0;endi=i+1;endelse i=i-1;endendhold online([startframe*256 startframe*256],[-0.50.5],'Color','b');%»-³öÆðʼµãhold online([endframe*256 endframe*256],[-0.5 0.5],'Color','r');%»-³öÖÕÖ¹µã五、实验结论分析：（1）从图中明显可以看出，浊音信号的具有明显的周期性，其自相关函数和平均幅度差函数也表现出周期性。