查找特征值的代码

时间:2014-03-14 10:50:07

标签: c matrix eigenvector eigenvalue

您好我一直在尝试编码以查找n * n矩阵的特征值。但是我无法想到它应该是什么算法。

步骤1:找到det(A-(lamda)* I)= 0

一般矩阵的算法应该是什么,用于查找lamda?

我已经编写了查找矩阵行列式的代码,可以在我们的算法中使用。

请帮助。非常感谢。

1 个答案:

答案 0 :(得分:0)

[假设你的矩阵是hermitian(简单地说是对称),所以特征向量是实数]

在计算中,你不能使用行列式来求解eignenvectors和特征值。它在数字上太慢而且不稳定。

您所做的是应用转换(户主缩减)以将矩阵缩减为三对角形式。

之后,您可以应用所谓的 QL 算法。

作为起点,请查看数字食谱(www.nr.com)中的tred2tqli。这些是我刚刚描述的算法。

请注意,这些例程还可以恢复候选特征向量