Mercurial > hg > smallbox
diff DL/Majorization Minimization DL/ExactDicoRecovery/mapfn.m @ 159:23763c5fbda5 danieleb
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/mapfn.m Wed Aug 31 10:43:32 2011 +0100 @@ -0,0 +1,16 @@ +%% Maximum A Posteriori Dictionary Learning with the constraint on the column norms %%%%% +function [Phiout,unhatnz] = mapfn(Phi,x,unhat,mu,maxIT,eps,phim,res) +K = Phi; +B = zeros(size(Phi,1),size(Phi,2)); +i = 1; +while (sum(sum((B-K).^2))>eps)&&(i<=maxIT) + B = K; + E = x-K*unhat; + K = K+mu*(phim*E*unhat'-trace(unhat*E'*K)*K); + i = i+1; +end +%%% depleted atoms cancellation %%% +[Y,I] = sort(sum(K.^2),'descend'); +RR = phim; +Phiout = K(:,I(1:RR)); +unhatnz = unhat(I(1:RR),:); \ No newline at end of file