comparison examples/SMALL_test_coherence2.m @ 183:0d7a81655ef2 danieleb

removed cumulative coherence calculation

+%
+%
+%
+
 clc, clear, close all
 
 %% Parameteres
 nTrials   = 10;										%number of trials of the experiment
 

 	end
 end
 
 %% Evaluate coherence and snr of representation for the various methods
 sr = zeros(size(SMALL.DL));				%signal to noise ratio
-mu = zeros(nTrials,nInitDicts,nCorrLevels,nDecorrAlgs,blockSize);	%cumulative coherence
+mu1 = zeros(nTrials,nInitDicts,nCorrLevels,nDecorrAlgs,blockSize);	%cumulative coherence
+mu2 = 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

 				sr(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs) = ...
 					snr(SMALL.Problem.b,SMALL.DL(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs).D*tempSolver.solution);
 				%calculate mu
 				dic(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs) = ...
 					dictionary(SMALL.DL(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs).D);
-				mu(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs,:) = ...
+				mu1(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs,:) = ...
 					dic(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs).cumcoherence;
+				mu2(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs,:) = ...
+					dic(iTrial,iInitDicts,iCorrLevels,iDecorrAlgs).cumcoherence(2);
 			end
 		end
 	end
 end
 
 %% Plot results
 minMu = sqrt((dictSize-blockSize)/(blockSize*(dictSize-1)));	%lowe bound on coherence
 initDictsNames = {'Data','Gabor'};
-dicDecorrNames = {'IPR','INK-SVD','IP'};
+dicDecorrNames = {'IPR','INK-SVD'};
 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
+	for iDecorrAlgs=1:nDecorrAlgs-1
 		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

 	set(gca,'YLim',[0 25],'XLim',[0 1.4]);
 	legend([{'\mu_{max}'},dicDecorrNames,{'\mu_{min}'}]);
 	xlabel('\mu');
 	ylabel('SNR (dB)');
 end
+
+%% 
+mu2 = squeeze(mean(mu2,1));
+mu = squeeze(mean(mu,1));