获取对称矩阵的特征值和特征向量的简单 Lanczos 算法代码答案

【问题标题】：Simple Lanczos algorithm code to obtain eigenvalues and eigenvectors of a symmetric matrix获取对称矩阵的特征值和特征向量的简单 Lanczos 算法代码
【发布时间】：2013-01-11 09:36:00
【问题描述】：

我想使用 Lanczos 算法编写一个简单的程序（用 C 语言）。我遇到了一个 Matlab 示例，它帮助我进一步理解了算法，但是从这段代码中我找不到获取特征值和特征向量的方法。我可以遵循算法，但我认为我一定遗漏了一些东西。有人可以指导从这个示例中获取特征值，以便我可以理解该方法，然后用 C 对其进行编码吗？

% Create a random symmetric matrix 
D=6
for i=1:D,
    for j=1:i,
        A(i,j)=rand; 
        A(j,i)=A(i,j);
    end 
end

% Iteration with j=0 
r0 = rand(D,1); 
b0 = sqrt(r0'*r0); 
q1 = r0/b0; 
a1 = q1'*A*q1

%Iteration with j=1
r1 = A*q1 - a1*q1
b1 = sqrt(r1'*r1)
q2 = r1/b1;
a2 = q2'*A*q2

%Iteration with j=2
r2 = A*q2 - a2*q2 - b1*q1;
b2 = sqrt(r2'*r2)
q3 = r2/b2
a3 = q3'*A*q3

% Create Matrix Q
Q = [q1 q2 q3];

%Check orthogonality
EYE = Q'*Q
T = Q'*A*Q

【问题讨论】：

标签： matlab linear-algebra eigenvector eigenvalue lanczos

【解决方案1】：

在初始Lanczos方法中，首先要计算矩阵A的最大特征值。然后计算与该特征值对应的特征向量。计算完这两个对象后，您可以减少您在一个对象上使用的矩阵的维数，然后找到新矩阵的最大特征值。并且您将迭代这 m 次，其中 m 是初始矩阵 A 的维数。

但是如果你想同时计算所有的特征值，你要使用Paige iterative procedure (see in middle) 首先你建立三对角矩阵。然后你可以使用一种众所周知的快速算法来计算它的特征值，因为这样的矩阵非常稀疏，通过上面文章中指定的公式，你可以很容易地计算初始矩阵的特征值及其对应的特征向量。

【讨论】：