在我的 python 代码中，我想解决多项式特征值问题:
A0 + lambda*A1 + lambda^2*A2 + lambda^3*A3 + .... = 0
其中 An 是稠密矩阵，lambda 是常量。在 matlab 中，可以使用 polyeig function 解决这个问题。 .似乎 scipy 中没有等效的功能。到目前为止，我能想到的唯一方法就是形成相应的伴随矩阵。这产生了一个等效的线性特征值问题，可以将其提供给现有的 scipy 求解器，但它要大得多，我相信它可能是病态的。
谁能推荐一个现有的、开源的或免费提供的库来解决这个问题？我对可以通过 f2py 或 C/C++ 库链接到 cython 链接的 fortran 库感到非常满意。
编辑:对于任何有兴趣在 python 中解决非线性特征值问题的人，可以找到我自己编写的解决此问题的代码 here .请注意，我处理的是非线性特征值问题的更一般情况(从某种意义上说，它对 lambda 具有非线性依赖性)。要了解该方法，请阅读代码注释中提到的论文。
最佳答案
This discussion指出将多项式特征值问题转化为广义特征值问题的一般方法，稍后可以使用 scipy's linear algebra 求解职能。希望这对您有所帮助!
关于python - 如何在python中求解多项式特征值？，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/8252428/