语音信号在不同信噪比条件下的端点检测matlab仿真+仿真录像

function [voiceseg,vsl,SF,NF]=vad_ezrm(x,wlen,inc,NIS) x=x(:); % 把x转成列数组 maxsilence = 15; % 初始化 minlen = 5; status = 0; count = 0; silence = 0; %计算短时能量 y=enframe(x,wlen,inc)'; % 分帧 fn=size(y,2); % 帧数 amp=sum(y.^2); % 求取短时平均能量 zcr=zc2(y,fn); % 计算短时平均过零 ampm = multimidfilter(amp,5); % 中值滤波平滑处理 zcrm = multimidfilter(zcr,5); ampth=mean(ampm(1:NIS)); % 计算初始无话段区间能量和过零率的平均值 zcrth=mean(zcrm(1:NIS)); amp2=1.2*ampth; amp1=1.5*ampth; % 设置能量和过零率的阈值 zcr2=0.8*zcrth; %开始端点检测 xn=1; for n=1:fn switch status case {0,1} % 0 = 静音, 1 = 可能开始 if amp(n) > amp1 % 确信进入语音段 x1(xn) = max(n-count(xn)-1,1); status = 2; silence(xn) = 0; count(xn) = count(xn) + 1; elseif amp(n) > amp2 | ... % 可能处于语音段 zcr(n) < zcr2 status = 1; count(xn) = count(xn) + 1; else % 静音状态 status = 0; count(xn) = 0; x1(xn)=0; x2(xn)=0; end case 2, % 2 = 语音段 if amp(n) > amp2 | ... % 保持在语音段 zcr(n) < zcr2 count(xn) = count(xn) + 1; else % 语音将结束 silence(xn) = silence(xn)+1; if silence(xn) < maxsilence % 静音还不够长，尚未结束 count(xn) = count(xn) + 1; elseif count(xn) < minlen % 语音长度太短，认为是噪声 status = 0; silence(xn) = 0; count(xn) = 0; else % 语音结束 status = 3; x2(xn)=x1(xn)+count(xn); end end case 3, % 语音结束，为下一个语音准备 status = 0; xn=xn+1; count(xn) = 0; silence(xn)=0; x1(xn)=0; x2(xn)=0; end % fprintf('%4d \n',status); end el=length(x1); if x1(el)==0, el=el-1; end % 获得x1的实际长度 if x2(el)==0 % 如果x2最后一个值为0，对它设置为fn fprintf('Error: Not find endding point!\n'); x2(el)=fn; end SF=zeros(1,fn); % 按x1和x2，对SF和NF赋值 NF=ones(1,fn); for i=1 : el SF(x1(i):x2(i))=1; NF(x1(i):x2(i))=0; end speechIndex=find(SF==1); % 计算voiceseg voiceseg=findSegment(speechIndex); vsl=length(voiceseg);