我正在使用迭代子空间算法(dsrrit)来获得特征向量方程的特征值
在哪里是通常的拉普拉斯算子。该算法返回一组特征值和一个子空间跨越的特征向量.
文献称
“程序不会产生一组与计算的特征值相对应的特征向量”......“如果需要明确的特征向量,它们可以通过评估特征向量来获得并申请(2)。”
但我已经在返回的子空间中绘制了向量它们确实与预期的特征向量相匹配. 是否有特定条件可以保证这一点?(我的特征向量是实数,我的算子用对称矩阵表示)
我正在使用迭代子空间算法(dsrrit)来获得特征向量方程的特征值
文献称
“程序不会产生一组与计算的特征值相对应的特征向量”......“如果需要明确的特征向量,它们可以通过评估特征向量来获得并申请(2)。”
但我已经在返回的子空间中绘制了向量它们确实与预期的特征向量相匹配. 是否有特定条件可以保证这一点?(我的特征向量是实数,我的算子用对称矩阵表示)
我没有详细阅读它,但来自您的链接文件(第 2 页):
SRRIT 试图计算一个嵌套的正交基序列. 具体来说,如果一切顺利,子程序会生成一个矩阵 Q,其正交列具有以下属性:然后跨度.
这意味着如果您的矩阵只有实特征值,则, 然后将是一个特征向量,利用矩阵是对称的并且特征向量将彼此正交的事实。
如果是两个复共轭特征值,那么对应的列不会是特征向量,因为它们是真实的。
原因是是原始矩阵的特征向量实际上是由于该矩阵的对称性。是的,由于对称性,特征值确实是实数。但也由于对称性,从返回的矩阵dsrrit
不仅是上三角矩阵,而且本身也是对称的。因此它必须是一个对角矩阵。这矩阵可以被认为是表示简化特征系统的特征向量。通常,原始矩阵的特征向量计算为. 在这种情况下,因为是对角线,它只是从矩阵。
特征值都是实数不足以使特征向量成为. 构造具有实特征值的非对称矩阵很容易。在这种情况下,将不是对角线,因此特征向量必须计算为.