DTW语音识别算法（部分）(2)

Status=1;

HoldTime=HoldTime+1; else

Status=0; HoldTime=0; end case 2,

if amp(n)>AmpLow | zcr(n)>ZcrLow HoldTime=HoldTime+1; else

SilenceTime=SilenceTime+1; if SilenceTime

elseif (HoldTime-SilenceTime)

Status=3; end end case 3, break; end

if Status==3 break; end end

HoldTime =HoldTime-SilenceTime; EndPoint=StartPoint+HoldTime; disp('显示端点??');

figure,subplot(3,1,1); plot(k);

axis([1,length(k),min(k),max(k)]); title('(I) 语音信号');

xlabel('Sample'); ylabel('Speech');

line([StartPoint*FrameInc,StartPoint*FrameInc],[min(k),max(k)],'color','red'); line([EndPoint*FrameInc,EndPoint*FrameInc],[min(k),max(k)],'color','red'); subplot(3,1,2);

plot(zcr);

axis([1,length(zcr),0,max(zcr)]); title('(II) 短时过零率');

xlabel('Frame'); ylabel('ZCR');

line([StartPoint,StartPoint],[0,max(zcr)],'Color','red'); line([EndPoint,EndPoint],[0,max(zcr)],'Color','red'); subplot(3,1,3); plot(amp);

axis([1,length(amp),0,max(amp)]); title('(III) 短时能量');

xlabel('Frame'); ylabel('Energy');

line([StartPoint,StartPoint],[0,max(amp)],'Color','red'); line([EndPoint,EndPoint],[0,max(amp)],'Color','red');

附录5:DTW算法 function dist = dtw(test, ref) global x y_min y_max global t r global D d global m n

t = test; r = ref;

n = size(t,1); m = size(r,1);

d = zeros(m,1);

D = ones(m,1) * realmax; D(1) = 0;

% 如果两个模板长度相差过多，匹配失败 if (2*m-n<3) | (2*n-m<2) dist = realmax; return end

% 计算匹配区域 xa = round((2*m-n)/3);

xb = round((2*n-m)*2/3);

if xb>xa

%xb>xa, 按下面三个区域匹配 % 1 :xa % xa+1:xb % xb+1:N for x = 1:xa y_max = 2*x; y_min = round(0.5*x); warp end

for x = (xa+1):xb y_max = round(0.5*(x-n)+m); y_min = round(0.5*x); warp end

for x = (xb+1):n y_max = round(0.5*(x-n)+m); y_min = round(2*(x-n)+m); warp end elseif xa>xb

%xa>xb, 按下面三个区域匹配 % 0 :xb % xb+1:xa % xa+1:N for x = 1:xb y_max = 2*x; y_min = round(0.5*x); warp end

for x = (xb+1):xa y_max = 2*x; y_min = round(2*(x-n)+m); warp end

for x = (xa+1):n y_max = round(0.5*(x-n)+m); y_min = round(2*(x-n)+m);

warp end

elseif xa==xb

%xa=xb, 按下面两个区域匹配 % 0 :xa % xa+1:N for x = 1:xa y_max = 2*x; y_min = round(0.5*x); warp end

for x = (xa+1):n y_max = round(0.5*(x-n)+m); y_min = round(2*(x-n)+m); warp end end

%返回匹配分数

dist = D(m); function warp

global x y_min y_max global t r global D d global m n d = D;

for y = y_min:y_max D1 = D(y); if y>1 D2 = D(y-1); else

D2 = realmax; end if y>2 D3 = D(y-2); else

D3 = realmax; end

d(y) = sum((t(x,:)-r(y,:)).^2) + min([D1,D2,D3]); end D = d;

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库DTW语音识别算法（部分）(2)在线全文阅读。