Mercurial > hg > smallbox

comparison util/classes/dictionaryMatrices/iterativeprojections.m @ 184:8fc38e8df8c6 danieleb

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Minor edits
author Daniele Barchiesi <daniele.barchiesi@eecs.qmul.ac.uk>
date Fri, 27 Jan 2012 13:18:50 +0000
parents 68fb71aa5339
children
comparison
equal deleted inserted replaced
183:0d7a81655ef2 184:8fc38e8df8c6
30 srs(iIter) = SNR(dic); 30 srs(iIter) = SNR(dic);
31 31
32 % calculate gram matrix 32 % calculate gram matrix
33 G = dic'*dic; 33 G = dic'*dic;
34 34
35 % project into the structural constraint set 35 % project onto the structural constraint set
36 H = zeros(size(G)); %initialise matrix 36 H = zeros(size(G)); %initialise matrix
37 ind1 = find(abs(G(:))<=mu); %find elements smaller than mu 37 ind1 = find(abs(G(:))<=mu); %find elements smaller than mu
38 ind2 = find(abs(G(:))>mu); %find elements bigger than mu 38 ind2 = find(abs(G(:))>mu); %find elements bigger than mu
39 H(ind1) = G(ind1); %copy elements belonging to ind1 39 H(ind1) = G(ind1); %copy elements belonging to ind1
40 H(ind2) = mu*sign(G(ind2)); %threshold elements belonging to ind2 40 H(ind2) = mu*sign(G(ind2)); %threshold elements belonging to ind2