Mercurial > hg > smallbox
comparison util/classes/dictionaryMatrices/grassmannian.m @ 169:290cca7d3469 danieleb
Added dictionary decorrelation functions and test script for ICASSP paper.
author | Daniele Barchiesi <daniele.barchiesi@eecs.qmul.ac.uk> |
---|---|
date | Thu, 29 Sep 2011 09:46:52 +0100 |
parents | 1495bdfa13e9 |
children | 68fb71aa5339 |
comparison
equal
deleted
inserted
replaced
168:ff866a412be5 | 169:290cca7d3469 |
---|---|
46 | 46 |
47 % A = normc(A); %normalise dictionary | 47 % A = normc(A); %normalise dictionary |
48 | 48 |
49 [V_gram Sigma_gram] = svd(G); %calculate svd decomposition of gramian | 49 [V_gram Sigma_gram] = svd(G); %calculate svd decomposition of gramian |
50 Sigma_new = sqrt(Sigma_gram(1:n,:)).*sign(Sigma); %calculate singular values of dictionary | 50 Sigma_new = sqrt(Sigma_gram(1:n,:)).*sign(Sigma); %calculate singular values of dictionary |
51 A = Uinit*Sigma_new*V_gram'; %update dictionary | 51 A = Uinit*Sigma_new*V_gram'; %update dictionary |
52 | |
53 % param.step = 0.01; | |
54 % param.reg = 0.01; | |
55 % param.nIter = 20; | |
56 % A = rotatematrix(initA,A,'linesearchlie',param); | |
52 | 57 |
53 % %% Debug visualization function | 58 % %% Debug visualization function |
54 % function plotcart2d(A) | 59 % function plotcart2d(A) |
55 % compass(A(1,:),A(2,:)); | 60 % compass(A(1,:),A(2,:)); |