作者热门文章
- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
在我的 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/
我是一名优秀的程序员,十分优秀!