下载地址: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 =