解决Numpy中的广义特征值问题

解决Numpy中的广义特征值问题

对于对称密集矩阵，您可以

scipy.linalg.eigh()

用来解决此广义特征值问题：

from scipy.linalg import eigheigvals, eigvecs = eigh(A, B, eigvals_only=False)

您会发现这

eigvecs

很复杂

ndarray

，所以也许您必须使用

eigvecs.real

…

在同一个模块中

eigvalsh()

，对于您的情况，该模块可能会执行得更快，但不会返回特征向量。