是否有能够解决以下特征值的 fortran 求解器？

Question

A 是一个 N 乘 N 矩阵。 I是(N-2)乘(N-2)的单位矩阵。 B 是另一个 N × N 矩阵，定义为

    B=[I 0 0;
       0 0 0;
       0 0 0]

。 x 是一个包含 N 个元素的数组。如何解决以下形式的特征值 A x=cB x，其中c是特征值， 通过使用特征值求解器？

Answer 1

您可以查看Lapack 库，它为广义矩阵的特征值问题提供了解决方案。根据您的数据类型和矩阵类型，您将需要使用不同的子例程。 在这方面看看here。另请查看 here 以了解用作矩阵类型的命名法。

最后，前一段时间我写了this module 来给出一些 Lapack 功能的实现示例，包括特征值问题。您可以在那里找到单值广义矩阵 (sgeev)。