Dav*_*veP 6 python matlab numpy scipy
在我的python代码中,我想解决多项式特征值问题:
A0 + lambda*A1 + lambda^2*A2 + lambda^3*A3 + .... = 0
An密集矩阵在哪里,并且lambda是常数.在matlab中,可以使用polyeig函数解决这个问题.似乎scipy中没有相同的功能.到目前为止,我能想到的唯一方法就是形成相应的伴随矩阵.这会产生一个等效的线性特征值问题,可以给现有的scipy求解器,但是它要大得多,我相信它可能会非常恶劣.
任何人都可以建议一个现有的,开源的或免费提供的库来解决这个问题吗？我很满意fortran库,可以通过f2py或C/C++库链接到cython链接.
编辑:对于任何有兴趣在python中解决非线性特征值问题的人来说,我自己编写的用来解决这个问题的代码可以在这里找到.请注意,我处理非线性特征值问题的更一般情况(在某种意义上它与λ有非线性依赖关系).要了解该方法,请阅读代码注释中提到的文章.