matlab端点检测程序(完全注释版)

《matlab端点检测程序(完全注释版)》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、端点检测程序第一部分：常数设置%常数设置FrameLen=240;%指定巾贞长Framelnc=80;%指定巾贞移，每一巾贞中未重叠的部分ampl=10;%初始短时能量高门限amp2=2;%初始短时能量低门限zcrl=10;%初始过零率高门限zcr2=5;%初始过零率低门限maxsilence=8;%8^10ms=80ms%语音段中允许的最大静音长度，如果语音段中的静咅帧数未超过此值，则认为语音还没结束；如果超过了该值，则对语咅段长度count进行判断，若countminlen,则认

2、为语音段结束；minlen=15;%15*10ms=150ms%语音段的最短长度，若语音段长度小于此值，则认为其为一段噪音status=0;%初始状态为静音状态count=0;%初始语音段长度为0silence=0;%初始静音段长度为0第二部分：过零计算tmpl=enframe(x(I:end-1),FrameLen,Framelnc);%分巾贞处理，tmpl和tmp2力分巾贞后形成的二维数组tmp2=enframe(x(2:end)，FrameLen,Framelnc);%语音信号是一种典型的非平稳信号，但是语音信号具有短时平稳性，因此在处理中要对采样的语音

3、信号进行分割成一帧一帧的短吋语音序列。分帧处理有利于对语音信号进行准确的分析并且能够提高识别率这时再分别求出每帧的短时能量和短时过零率signs=(tmpl.*tmp2)<0;%当tmpl.*tmp<0的时候，说明tmpl〉0，tmp2<0或tmpl<0，tmp2〉0;即信号过零点。diffs=(tmpl-tmp2)>0.02;%当信号过零点，而tmpl与tmp2距离又过近(<0.02)的话就认为是噪音，舍去.zcr=sum(signs.*diffs，2);%定义语音信号Xn(m)的过零率为Zn，1yv-i则Zn=—ZsgnfXn(m)]-sgn[-1)111

4、m=0第三部分：计算短时能量amp=sum(abs(enframe(filter([1-0.9375]，1，x)，FrameLen，Framelnc))，2)；%过零率有两类重要的应用：用于粗略地描述信号的频谱特性，就是用将为若干个通道，对各通道进行短时平均过零率和短时能量的计算，即可粗略地估计频谱特性；用于判别清音和浊音、有话和无话。计算过零率容易受低频干扰，特别是50Hz交流干扰的影响。解决这个问题的办法，一个是做高通滤波器或带通滤波，减小随机噪声的影响(filter函数滤波器)%调整能量门限ampl=min(ampl,max(amp)/4);amp2=m

5、in(amp2,max(amp)/8);AM%若第N帧音信号Xn(m)短时能量为En，则⑽/H=0第四部分:算法说明forn=l:length(zcr)goto=0;switchstatus%status是判断语音在哪个状态的函数，当status=O时，为静音状态status=l时，为可能开始状态；当status=2时，为语音状态，当status=3时，为结束状态；case{0,1}%0=静音，1=可能开始ifamp(n)>ampl%确信进入语音段xl(end+l)=max(n-count-l，l);status=2;%进入语音段silence=0;count

6、=count+1;elseifamp(n)>amp2

7、...%可能处于语音段%从静音段开始，当amp(n)高于短时能量低门限时，可能进入语音段，如果信号长度大于最短语音长度的话，即为语音段，否则为静音段。在语音段的时候，当amp(n)高于短时能量低门限时，将保持在语音段，当amp(n)低于短时能量低门限时，进入静音段。zcr(n)>zcr2%从静音段开始，当zcr(n)高于过零率低门限时，可能进入语音段，如果信号长度大于最短语音长度的话，即为语音段，否则为静音段。在语音段的时候，当ar(n)高于过零率低门限时，将保持在语音段，当amp(n)低于过零率低门限时，

8、进入静音段。count=count+1;else%静音状态status=0;%进入静音段count=0;else%语音结束status=3;endendcase3,status=O;x2(end+1)=xl(end)+count-silence/2-1;endend%count=count-silence/2;%v_count(i)=v_count(i)+v_count(i-1);%v_silence(i)=v_count(i)+v_silence(i);iflength(x2)

9、ot(311)(将多个图画到一个平面的