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:matrixSi