eigs(A,k,'sm')
返回k个最小幅度的特征值。但是,它会发出警告,并非所有特征值都收敛,也不能给出k
最小的特征值。
另一方面,使用eig(full(A))
,自己进行排序并提取最小k
个特征值确实给了我相同的答案。
通过互联网搜索很长一段时间,我听到了很多关于eigs
的负面评论。但是,因为我的矩阵A
稀疏而且很大。 eigs
似乎是速度最好的方法,有没有什么方法可以使eigs(A,k,'sm')
正确并输出我想要的正确特征值。或者是否有另一种更可靠,更快速的数值方案来寻找稀疏矩阵的特征值?