Matlab根据特征值排序特征向量

介绍

因为matlab中像[X,D]=eig(X)这样的函数计算特征值和特征向量之后是分别返回的，尽管可以使用sort()函数对diag(D)进行排序，但是特征向量并不能随着特征值的排序而排序，为了能够排序特征值的同时，使与特征值相对应的特征向量也排序，我们可以利用sortrows(X,col,direction)按某列排序，同时整行随着排序的特性来实现以上功能。思路仅供参考。

代码

function [sorted_X,sorted_D]=sorted(X,direction)
    [X1,D]=eig(X); %对矩阵求特征值和特征向量
    col_vector=diag(D); %将特征值矩阵转换为只有一列的向量
    X1=[X1',col_vector] %将特征值整列追加到特征向量最后
    [row,col]=size(X1); 
    sorted_result=sortrows(X1,col,direction) %按照最后一列排序
    sorted_X=sorted_result(:,1:end-1)'; %返回排序后特征向量
    sorted_D=sorted_result(:,end);  %返回排序后的特征值
end

测试

测试代码如下

X=[1 -2 0;-2 5 0;0 0 2];
[X1,D]=eig(X)
[X,D]=sorted(X,'decend')

运行结果如下：
下面两张图是用eig()计算的结果


下图是将特征值追加到特征向量末尾的结果

下图是根据最后一列（也就是特征值）降序排列的结果

最后分别返回排序后的特征向量和特征值

上图中，X矩阵的第1列是最大特征值对应特征向量。即，5.8284对应的特征向量是(-0.3827,0.9239,0)