Python中两个平方矩阵的本征分解

Question

在matlab中，无论是否存在对称或非对称乘积，我们都可以选择找到两个矩阵的本征分解

A = [1 3; 4 9];
B = [4 7; 9 16];
[Vec,Val]=eig(A,B)

向量是

`[-1,-1;0.54,0.85]`  

和值是

[-3.79,0;0,0.79]

我已经签入python numpy.linalg，但是没有这样的选项。所有eig变体仅接受一个参数。有没有办法在python中处理这个问题

Answer 1

您可以使用scipy.linalg.eig：

from scipy import linalg
linalg.eig(A, B)

其中A = [[1,3],[4,9]]和B = [[4,7], [9,16]]是您的两个矩阵。