matlab求矩阵前k小的特征值和对应的特征向量

最新推荐文章于 2024-10-09 10:12:01 发布

原创

最新推荐文章于 2024-10-09 10:12:01 发布 · 1w 阅读

13 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #function #测试

在MATLAB中，使用eig函数得到的特征值是随机排列的，若要获取绝对值前k小的特征值及其对应的特征向量，标准的eigs函数是用于求解最大特征值。可以参考论坛中zqx917提出的方法来实现这一目标。

matlab中eig函数求特征值和特征向量，但特征值是随机排列的。要想取得前k个特征值（绝对值前k小的特征值）和对应的特征向量，还得另想办法。

matlab自带的eigs是用来求绝对值前k大的特征值。

实现了http://www.ilovematlab.cn/thread-50763-1-1.html 的zqx917提到的方法

function [firstkVector,firstKlemda]=firstKEigenVector(A,k)
%输入一个方阵
%将特征值和特征向量重排序，按照特征值的绝对值升序来排
%取前k个特征值和对应的特征向量

%测试用
%A=magic(5);

matrixSize=size(A);
[v,d] = eig(A) 
v = [v ; diag(d)' ; abs(diag(d))'  ]' ;%将特征值的信息加入到矩阵便于按特征值排序
v = sortrows(v,matrixSize+2) ;%升序排列
firstKlemda = v(1:k,matrixSize(1)+1) ;
firstkVector = v(1:k,1:matrixSize)';

以运行幻方矩阵为例