tInit=clock;
fprintf('Loading output/mat/gmmparam.mat\n');
load output/mat/gmmParam.mat
fprintf('Loading output/mat/speakerSet.mat\n');
load output/mat/speakerSet.mat
if maxSpeakerNum<length(speakerSet1)
	speakerSet1(maxSpeakerNum+1:end)=[];
	speakerSet2(maxSpeakerNum+1:end)=[];
end

speakerSet=speakerSet2;		% Use the test data
speakerNum=length(speakerSet);

fprintf('Compute the log prob. of each sentence to every speakers...\n');
for i=1:speakerNum
	fprintf('%d/%d\n', i, speakerNum);
	sentenceNum=length(speakerSet(i).sentence);
	for j=1:sentenceNum
		for k=1:speakerNum
			logProb=gmmEvalMex(speakerSet(i).sentence(j).mfcc, gmm(k).mean, gmm(k).covariance, gmm(k).weight);
			speakerSet(i).sentence(j).logProb(k)=mean(logProb);
		end
	end
end

for i=1:speakerNum
	sentenceNum=length(speakerSet(i).sentence);
	for j=1:sentenceNum
		if length(speakerSet(i).sentence(j).logProb)~=speakerNum
			fprintf('i=%d, j=%d\n', i, j);
			error('Wrong');
		end
	end
end

for i=1:speakerNum
%	fprintf('%d/%d\n', i, speakerNum);
	sentenceNum=length(speakerSet(i).sentence);
	for j=1:sentenceNum
%		fprintf('\t%d/%d\n', j, sentenceNum);
		% ====== Positive cases
		logProb=speakerSet(i).sentence(j).logProb;
		speakerSet(i).sentence(j).positiveLogProb1=logProb(i);		% logProb for 本身
		logProb(i)=[];
		speakerSet(i).sentence(j).positiveLogProb2=max(logProb);	% 除掉本身之外的最大機率值
		% ====== Negative cases (i: index of the right person, which should be deleted since this is a negative case, index1: index of the imitated person)
		logProb=speakerSet(i).sentence(j).logProb;
		while 1,
			index1=ceil(rand*speakerNum);
			if index1~=i, break; end
		end
		speakerSet(i).sentence(j).negativeLogProb1=logProb(index1);
		logProb([i, index1])=[];
		speakerSet(i).sentence(j).negativeLogProb2=max(logProb);
	end
end

positiveData=[];
negativeData=[];
for i=1:speakerNum
	positiveData=[positiveData, [speakerSet(i).sentence.positiveLogProb1]-[speakerSet(i).sentence.positiveLogProb2]];
	negativeData=[negativeData, [speakerSet(i).sentence.negativeLogProb1]-[speakerSet(i).sentence.negativeLogProb2]];
end
figure; roc(negativeData, positiveData, 1);
eval(['print -dpng output\png\', mfilename, '01']);
figure; rocPlot(negativeData, positiveData);
eval(['print -dpng output\png\', mfilename, '02']);
figure; rocPlot(negativeData, positiveData, 'log');
eval(['print -dpng output\png\', mfilename, '03']);
fprintf('%s ===> %g seconds\n', mfilename, etime(clock, tInit));