wave2mfccTest

 測試不同的程式所產生的 MFCC 的差異

0001 % 測試不同的程式所產生的 MFCC 的差異
0002 
0003 waveFile='主人下馬客在船.wav';
0004 waveFile='ws2_16khz.wav';
0005 [y, fs, nbits]=wavread(waveFile);
0006 %sound(y, fs);
0007 %y=y*2^nbits/2;
0008 feature=wave2mfccMex(y, fs);
0009 mesh(feature);
0010 axis tight; box on;
0011 
0012 dos(['hcopy\HCopy -C hcopy\htk.cf ', waveFile, ' file.fea']);
0013 feature2=readFeaFile('file.fea');
0014 
0015 fprintf('Difference between Vincent & hcopy = %g\n', max(max(abs(feature-feature2))));