testdtw.m

patten regnization source从1-14章能运行

disp('1窜拌:' )
disp('澜己牢侥饶焊客 矫氰侩 澜己俊 措茄 傈贸府')
for i=1:5
	fname = sprintf('%dr.wav',i);
	x = wavread(fname);
	x = filter([1 -0.9375], 1, x);
	m = melcepst(x,16000,'M',12,24,256,80);
	ref(i).mfcc = m;
end

disp('矫氰侩 澜己俊 措茄 漂隆眠免...')
for i=1:5
	fname = sprintf('%dt.wav',i);
	x = wavread(fname);
	x = filter([1 -0.9375], 1, x);
	m = melcepst(x,16000,'M',12,24,256,80);
	test(i).mfcc = m;
end

disp('2窜拌: DTW 概莫...')
dist = zeros(5,5);
for i=1:5
for j=1:5
	dist(i,j) = dtw(test(i).mfcc, ref(j).mfcc);
end
end


disp('3窜拌: 搬苞...')

dist 
for i=1:5
	[d,j] = min(dist(i,:));
	fprintf('矫氰 澜己 %d ===> %d 肺 牢侥\n', i, j);
end