waveFile='\dataset\childSong4public\MIR-QBSH-corpus\waveFile\year2007/person00001/00011.wav';
[y, fs, nbits]=wavReadInt(waveFile);

vsdOpt=vsdOptSet;
%load bestInputIndex.mat
[inData, outData, inputName, annotation]=wave2feature(waveFile, vsdOpt);
%inData=inData(bestInputIndex, :);
%inputName=inputName(bestInputIndex);
if vsdOpt.useInputNormalize==1
	error('This is not implemneted yet!\n');
end

[feaDim, frameNum]=size(inData);

load gmmData.mat
load transProb.mat
transLogProb=log(transProb);
%transLogProb=zeros(size(transProb));
gmmcPrm=gmmData(index);
vsdResult=waveDecode(waveFile, vsdOpt, gmmcPrm, transProb, 1, inputName);

% ====== HMM recog. rate
computedClassIndex = vsdResult';
rr=sum(outData==computedClassIndex)/length(outData);
fprintf('HMM RR=%.2f%%\n', rr*100);

% ====== GMMC recog. rate
computedClassIndex = gmmcEval(inData, gmmcPrm);
rr=sum(outData==computedClassIndex)/length(outData);
fprintf('GMMC RR=%.2f%%\n', rr*100);