Как пометить комплексные собственные значения матрицы в порядке возрастания их мнимой части в Octave?

Question

Когда используется функция eig(a), Octave, похоже, не упорядочивает комплексные собственные значения в каком-либо конкретном порядке (тогда как для реальных собственных значений Octave хранит собственные значения и их собственные векторы в порядке возрастания).

Мне нужно чтобы упорядочить комплексные собственные значения (а также их соответствующие собственные векторы) матрицы в порядке возрастания их мнимой части. Есть ли встроенная процедура для этого?

Если нет, как это реализовать?

Answer 1

Функция sort имеет второй выходной аргумент, который задает индексы входных данных в отсортированном порядке. В вашем случае, если d является вектором собственных значений:

[~,I] = sort(imag(d));
d = d(I);

Вы можете использовать тот же I для сортировки соответствующих собственных векторов V в том же порядке:

V = V(:,I);