comparison examples/SMALL_test_coherence.m @ 167:8324c7ea6602 danieleb

Added symmetric de-correlation function, modified target de-correlation in test function.

 percActiveAtoms = 5; %percentage of active atoms
 
 % Dependent parameters
 nActiveAtoms = fix(blockSize/100*percActiveAtoms); %number of active atoms
 epsilon      = 1/sigma; %error constraint for sparse representation step (corresponds to noise applied to signals)
-minCoherence = sqrt((dictSize-blockSize)/(blockSize*(dictSize-1))); %target coherence (based on coherence lower bound)
+%minCoherence = sqrt((dictSize-blockSize)/(blockSize*(dictSize-1))); %target coherence (based on coherence lower bound)
+minCoherence = 0.4;		%target coherence
 
 % Initial dictionaries
 dctDict = dictionary('dct',blockSize,dictSize);
 dctDict = dctDict.phi;
 gaborParam = struct('N',blockSize,'redundancyFactor',2,'wd',@rectwin);

 
 %% Plot results
 minMu = sqrt((dictSize-blockSize)/(blockSize*(dictSize-1)));
 maxSNR = max(sigNoiseRatio);
 
-figure, subplot(2,2,1)
+figure, subplot(1,2,1)
 snrMat = buffer(sigNoiseRatio(1:9),3);
 bar(snrMat');
-title('SNR - KSVD Update')
+title('Signal to noise ratio')
 xlabel('Initial dictionary')
 ylabel('SNR (dB)')
-set(gca,'XTickLabel',{'data','dct','gabor'},'YLim',[0 maxSNR+1]);
+set(gca,'XTickLabel',{'data','dct','gabor'});
 legend('none','Mailhe','Tropp')
 grid on
 
-subplot(2,2,2), grid on
-snrMat = buffer(sigNoiseRatio(10:18),3);
-bar(snrMat');
-title('SNR - Mailhe Update')
-xlabel('Initial dictionary')
-ylabel('SNR (dB)')
-set(gca,'XTickLabel',{'data','dct','gabor'},'YLim',[0 maxSNR+1]);
-legend('none','Mailhe','Tropp')
-grid on
-
-subplot(2,2,3), hold on, grid on
-title('Coherence - KSVD Update')
+% subplot(2,2,2), grid on
+% snrMat = buffer(sigNoiseRatio(10:18),3);
+% bar(snrMat');
+% title('SNR - Mailhe Update')
+% xlabel('Initial dictionary')
+% ylabel('SNR (dB)')
+% set(gca,'XTickLabel',{'data','dct','gabor'},'YLim',[0 maxSNR+1]);
+% legend('none','Mailhe','Tropp')
+% grid on
+
+subplot(1,2,2), hold on, grid on
+title('Coherence')
 muMat = buffer(mu(1:9),3);
 line([0.5 3.5],[1 1],'Color','r');
 bar(muMat');
 line([0.5 3.5],[minMu minMu],'Color','k');
 set(gca,'XTick',1:3,'XTickLabel',{'data','dct','gabor'},'YLim',[0 1.05])
 legend('\mu_{max}','none','Mailhe','Tropp','\mu_{min}')
 ylabel('\mu')
 xlabel('Initial dictionary')
 
-subplot(2,2,4), hold on, grid on
-title('Coherence - Mailhe Update')
-muMat = buffer(mu(10:18),3);
-line([0.5 3.5],[1 1],'Color','r');
-bar(muMat');
-line([0.5 3.5],[minMu minMu],'Color','k');
-set(gca,'XTick',1:3,'XTickLabel',{'data','dct','gabor'},'YLim',[0 1.05])
-legend('\mu_{max}','none','Mailhe','Tropp','\mu_{min}')
-ylabel('\mu')
-xlabel('Initial dictionary')
+% subplot(2,2,4), hold on, grid on
+% title('Coherence - Mailhe Update')
+% muMat = buffer(mu(10:18),3);
+% line([0.5 3.5],[1 1],'Color','r');
+% bar(muMat');
+% line([0.5 3.5],[minMu minMu],'Color','k');
+% set(gca,'XTick',1:3,'XTickLabel',{'data','dct','gabor'},'YLim',[0 1.05])
+% legend('\mu_{max}','none','Mailhe','Tropp','\mu_{min}')
+% ylabel('\mu')
+% xlabel('Initial dictionary')