matlab程序使用动态时间规整语音识别.zip

function [x1,x2,Voiceframe]=frameDetect(Y,FS) %归一化幅度至【－1，1】 Y=double(Y); Y=Y-mean(Y); Y=Y/max(abs(Y)); %帧长10ms 帧内交织5ms FrameLength=round(FS*20/1000); % FrameLength=240; FrameInc=round(FS*10/1000); % FrameInc=80; % Amp1=50; % Amp2=10; Zcr1=10; Zcr2=5; Delta =0.01;%过零率监测域值 MaxSilence=3; % 30ms MinLength=15; % 150ms Status=0; % 0静音,1开始,2语音段,3结束 Count=0; Silence=0; %计算过零率 frametemp=enframe(Y,FrameLength,FrameInc); tmp1=enframe(Y(1:length(Y)-1),FrameLength,FrameInc); tmp2=enframe(Y(2:length(Y)),FrameLength,FrameInc); signs=(tmp1.*tmp2)<0; diffs=abs((tmp1-tmp2))>Delta; zcr=sum(signs.*diffs,2); %计算短时能量 %滤波 filter_Y=filter([1 -0.9375],1,Y); %分帧 %短时能量 frame=enframe(filter_Y,FrameLength,FrameInc); %三种能量 amp1=sum(abs(frame),2); amp2=sum(frame.*frame,2); amp3=sum(log(frame.*frame+eps),2); amp=amp1; %采用第1种幅度检测 %amp=amp-min(amp); %调整幅度门限 Max=max(amp); Min=min(amp); B=(Max-Min)/Max;%比例因子调整域值 factor=1.0; %Amp1=Min+(Max-Min)/3; %Amp2=Min+(Max-Min)/60; Amp1=Min+B.^factor*(Max-Min)/3; Amp2=Min+B.^factor*(Max-Min)/60; %开始端点检测 x1=0; x2=0; Count=0; Silence=0; for n = 1:length(zcr) switch Status case {0,1} if amp(n) > Amp1 x1=max(n-Count-1,1); Status =2; Silence=0; Count=Count+1; elseif amp(n)> Amp2 | zcr(n) > Zcr2 Status =1; Count=Count+1; else Status =0; Count =0; end case 2, if amp(n)> Amp2 | zcr(n) > Zcr2 Count=Count+1; else Silence=Silence+1; if Silence< MaxSilence Count=Count+1; elseif Count< MinLength Status=0; Silence=0; Count=0; else Status=3; end end case 3 break; end end Count=Count-Silence; x2=x1+Count-1; Voiceframe=frametemp(x1:x2,:); % x2=x2-mod(x2,1); figure(1) subplot(3,1,1); plot(Y); axis([1 length(Y) -1 1]); ylabel('Speach'); line([x1*FrameInc x1*FrameInc],[-1 1],'Color','red'); line([x2*FrameInc x2*FrameInc],[-1 1],'Color','red'); subplot(3,1,2); plot(amp); axis([1 length(amp) 0 max(amp)]); ylabel('Energy'); line([x1 x1],[0 max(amp)],'Color','red'); line([x2 x2],[0 max(amp)],'Color','red'); subplot(3,1,3); plot(zcr); axis([1 length(zcr) 0 max(zcr)]); ylabel('ZCR'); line([x1 x1],[0 max(zcr)],'Color','red'); line([x2 x2],[0 max(zcr)],'Color','red'); % disp('[x1,x2,Amp1,Amp2,mean(amp(1:x1)),mean(amp)]'); %[x1,x2,Amp1,Amp2,mean(amp(1:x1)),mean(amp)]