Mercurial > hg > smallbox

diff DL/Majorization Minimization DL/ExactDicoRecovery/dict_update_MOD_cn.m @ 159:23763c5fbda5 danieleb

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge
author Daniele Barchiesi <daniele.barchiesi@eecs.qmul.ac.uk>
date Wed, 31 Aug 2011 10:43:32 +0100
parents b14209313ba4
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/DL/Majorization Minimization DL/ExactDicoRecovery/dict_update_MOD_cn.m	Wed Aug 31 10:43:32 2011 +0100
@@ -0,0 +1,20 @@
+function [Phiout,unhatnz] = dict_update_MOD_cn(x,unhat,cvset)
+%% Method of Optimized Direction (MOD) Dictionary Update followed by the
+%% constraint on the column norms %%
+
+%%
+eps = .01;
+K = x*unhat'*inv(unhat*unhat'+eps*eye(size(unhat,1)));   
+if cvset == 1,
+    K = K*diag(min(sum(K.^2).^(-.5),1)); % Projecting on the convex constraint set
+else
+    K = K*diag(sum(K.^2).^(-.5)); % Projecting on the Unit-Norm constraint set
+end
+
+%% depleted atoms cancellation %%%
+[Y,I] = sort(sum(K.^2),'descend');
+RR = sum(Y>=.01);
+Phiout = K(:,I(1:RR));
+unhatnz = unhat(I(1:RR),:);
+
+end
\ No newline at end of file