annotate general/numerical/matrix/eigsel.m @ 61:eff6bddf82e3
tip
Finally implemented perceptual brightness thing.
| author |
samer |
| date |
Sun, 11 Oct 2015 10:20:42 +0100 |
| parents |
db7f4afd27c5 |
| children |
|
| rev |
line source |
|
samer@4
|
1 function [r1,r2]=eigs(T,I)
|
|
samer@4
|
2 % eigs - Return selected eigs relative to order by magnitude
|
|
samer@4
|
3 %
|
|
samer@4
|
4 % eigs :: [[N,N]], [[M]->[N]] ~'M indices between 1 and N' -> [[M]].
|
|
samer@4
|
5 % eigs :: [[N,N]], [[M]->[N]] ~'M indices between 1 and N' -> [[N,M]], [[M]].
|
|
samer@4
|
6
|
|
samer@4
|
7 if nargout==1
|
|
samer@4
|
8 L0=eig(T);
|
|
samer@4
|
9 [dummy,ord]=sort(-abs(L0));
|
|
samer@4
|
10 r1=L0(ord(I));
|
|
samer@4
|
11 else
|
|
samer@4
|
12 [V0,D0]=eig(T); L0=diag(D0);
|
|
samer@4
|
13 [dummy,ord]=sort(-abs(L0));
|
|
samer@4
|
14 r1=V0(:,ord(I));
|
|
samer@4
|
15 r2=L0(ord(I));
|
|
samer@4
|
16 end
|
