有没有可以解决以下特征值的Fortran解算器?
问题描述:
A是N乘N矩阵。 I是(N-2)乘以(N-2)的单位矩阵。 B是另一个N乘N矩阵,其被定义为有没有可以解决以下特征值的Fortran解算器?
B=[I 0 0;
0 0 0;
0 0 0]
。 x是一个有N个元素的数组。如何解决以下形式的本征值 甲X = C 乙X,其中c是本征值, 通过使用特征值求解器?
答
您可以查看Lapack库,它提供了广义矩阵的特征值问题的解决方案。根据您的数据类型和矩阵类型,您将需要使用不同的子程序。 在这方面看看here。也可以看看here作为用作矩阵类型的术语。
最后,我之前写了this module来给出一些Lapack功能的实现示例,其中包括特征值问题。在那里你可以找到一个单值的广义矩阵(sgeev
)。
另外谷歌为[广义特征值问题](https://en.wikipedia.org/wiki/Eigendecomposition_of_a_matrix#Generalized_eigenvalue_problem) –