view DL/Majorization Minimization DL/ExactDicoRecovery/dict_update_MAP_cn.m @ 234:c96880c0c47c

renamed file.
author luisf <luis.figueira@eecs.qmul.ac.uk>
date Thu, 19 Apr 2012 17:21:05 +0100
parents b14209313ba4
children
line wrap: on
line source
function [Phiout,unhatnz] = dict_update_MAP_cn(Phi,x,unhat,mu,maxIT,eps,cvset)
%% Maximum A Posteriori Dictionary Update with the constraint on the column norms %%%%%



%%

K = Phi;
B = zeros(size(Phi,1),size(Phi,2));
i = 1;
%%
while (sum(sum((B-K).^2))>eps)&&(i<=maxIT)
    B = K;
    for j = 1:size(K,2),
        K(:,j) = K(:,j) - mu*(eye(size(K,1))-K(:,j)*K(:,j)')*(K*unhat-x)*unhat(j,:)';           
    end
    i = i+1;
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