寻找对应于最小特征值的特征向量

Question

我想找到矩阵U最小幅度特征值的相应特征向量。最简单的方法是什么？

目前我正在使用算法

[evecs, D] = eigs(U);
evals = diag(D);
smallesteig = inf;
for k = 1:length(evals)
    if (evals(k) < smallesteig)
        smallesteig = evals(k);
        vec = evecs(:, k);
    end
end

有更有效的方法吗？

Answer 1

有一个非常简单的简写：[V,D] = eigs(U,1,'SM')。

如果查看eigs文档，请说明：

EIGS(A,K,SIGMA) and EIGS(A,B,K,SIGMA) return K eigenvalues. If SIGMA is:
   'LM' or 'SM' - Largest or Smallest Magnitude
For real symmetric problems, SIGMA may also be:
   'LA' or 'SA' - Largest or Smallest Algebraic
   'BE' - Both Ends, one more from high end if K is odd
For nonsymmetric or complex problems, SIGMA may also be:
   'LR' or 'SR' - Largest or Smallest Real part
   'LI' or 'SI' - Largest or Smallest Imaginary part

因此，当{strong> S mallest M agnitude排序时，[V,D] = eigs(U,1,'SM')会返回U的 1 st特征值的特征向量和值。