comparison examples/SMALL_test_coherence2.m @ 174:dc2f0fa21310 danieleb

multiple trials with error bars

 	end
 end
 
 %% Evaluate coherence and snr of representation for the various methods
 sr = zeros(size(SMALL.DL));				%signal to noise ratio
-mu = zeros(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs,blockSize);	%cumulative coherence
+mu = zeros(nTrials,nInitDicts,nCorrLevels,nDecorrAlgs,blockSize);	%cumulative coherence
 dic(size(SMALL.DL)) = dictionary;		%initialise dictionary objects
 for iTrial=1:nTrials
 	for iInitDicts=1:nInitDicts
 		for iCorrLevels=1:nCorrLevels
 			for iDecorrAlgs=1:nDecorrAlgs

 end
 
 %% Plot results
 minMu = sqrt((dictSize-blockSize)/(blockSize*(dictSize-1)));	%lowe bound on coherence
 initDictsNames = {'Data','Gabor'};
-dicDecorrNames = {'Iter. Proj. + Rotation','INK-SVD','Iter. Proj.'};
-lineStyles     = {'ks-','kd-','ko-'};
+dicDecorrNames = {'IPR','INK-SVD','IP'};
+lineStyles     = {'k.-','r*-','b+-'};
 for iInitDict=1:nInitDicts
 	figure, hold on, grid on
 	title([initDictsNames{iInitDict} ' Initialisation']);
 	plot([1 1],[0 25],'k-');
 	for iDecorrAlgs=1:nDecorrAlgs
-		plot(mu(1,iInitDict,:,iDecorrAlgs,1),sr(1,iInitDict,:,iDecorrAlgs),...
-			lineStyles{iDecorrAlgs});
+		coherenceLevels = squeeze(mean(mu(:,iInitDict,:,iDecorrAlgs,1),1));
+		meanSNRs		= squeeze(mean(sr(:,iInitDict,:,iDecorrAlgs),1));
+		stdSNRs			= squeeze(std(sr(:,iInitDict,:,iDecorrAlgs),0,1));
+		errorbar(coherenceLevels,meanSNRs,stdSNRs,lineStyles{iDecorrAlgs});
 	end
 	plot([minMu minMu],[0 25],'k--')
 	
 	set(gca,'YLim',[0 25],'XLim',[0 1.4]);
 	legend([{'\mu_{max}'},dicDecorrNames,{'\mu_{min}'}]);