CODEWITHSUNDEEP
matlabmatrixlinear-algebraeigenvectoreigenvalue
kamaci
kamaci

Reputation: 75227

Does Matlab eig always returns sorted values?

I use a function at Matlab:

[V,D] = eig(C);

I see that V and D are always sorted ascending order. Does it always like that or should I sort them after I get V and D values?

Upvotes: 16

Views: 36485

Answers (2)

Rody Oldenhuis
Rody Oldenhuis

Reputation: 38042

If you want to guarantee sorted-ascending values, just do an extra 

if ~issorted(diag(D))
    [V,D] = eig(A);
    [D,I] = sort(diag(D));
    V = V(:, I);
end

to sort them the way you want.

Alternatively, use eigs: 

[V,D] = eigs(A,size(A,1)-1)

Upvotes: 33

user85109
user85109

Reputation:

V is NOT sorted in any order, except to correspond to the order of the associated eigenvalues. But perhaps you did not mean that.

The eigenvalues TEND to be in descending order, but this is not assured at all. They tend to be in order because the largest tend to trickle out of the algorithm on top. Eig has no sort at the end to ensure that fact.

I might point out the eigenshuffle tool, designed to take a sequence of eigenproblems, then resorting the eigenvalues (and the corresponding eigenvectors) so they are consistent along the sequence.

If you really need them certainly in decreasing order, then do a sort to ensure that fact. Make sure you also sort the vectors in the same order.

Upvotes: 25

Related Questions