下載地址： http://download.csdn.net/detail/fzxy002763/4082023

算法不是原創的，由于網上流傳版本很多，不過都是獨立給出一些m文件，不完整，故筆者整理了一下，調試通過，方便大家學習下dtw算法。

其中dtw_test是測試文件，其他皆為function文件

    disp('正在計算參考模板的參數...')
for i=1:5
	fname = sprintf('%da.wav',i);
	x=fname;
    [x,fs]=wavread(x);
	[x1 x2] = vad(x);
	m = mfcc(x);
	m = m(x1-2:x2-2,:);
	ref(i).mfcc = m;
   % soundview(x);
end

disp('正在計算測試模板的參數...')
for i=1:1
    %[x,fs]=wavread('E:\\3.wav')
	fname = sprintf('%db.wav',i);
	x=fname;
    [x,fs]=wavread(x);    
	[x1 x2] = vad(x);
	m = mfcc(x);
	m = m(x1-2:x2-2,:);
	test(i).mfcc = m;
end

disp('正在進行模板匹配...')
dist = zeros(1,5);
for i=1:1
for j=1:5
	dist(i,j) = dtw(test(i).mfcc, ref(j).mfcc);
end
end

disp('正在計算匹配結果...')
for i=1:1
	[d,j] = min(dist(i,:));
	fprintf('測試模板 %d 的識別結果為：%d\n', i, j);
end

  

1.vad

2.識別

故是4a.wav匹配1b.wav，可以自行設定理解，PS： 錄音最好時間差不多，且大于1S，要不然可能會報index exceeds matrix dimensions這種錯誤。

matlab DTW語音識別算法調試