scipy.linalg.expm

scipy.linalg.expm(A)[源代码]

使用Pade近似计算矩阵指数。

参数
A(n,N)类阵列或稀疏矩阵

要求幂的矩阵。

退货
expm(N,N)ndarray

的矩阵指数 A

参考文献

1

Awad H.Al-Mohy和Nicholas J.Higham(2009)“矩阵指数的一种新的缩放和平方算法”。暹罗矩阵分析与应用杂志。31(3)。970-989页。ISSN 1095-7162

示例

>>> from scipy.linalg import expm, sinm, cosm

公式exp(0)=1的矩阵版本:

>>> expm(np.zeros((2,2)))
array([[ 1.,  0.],
       [ 0.,  1.]])

欧拉的身份(EXP(I theta) = cos(theta) + i sin(Theta))应用于矩阵的:

>>> a = np.array([[1.0, 2.0], [-1.0, 3.0]])
>>> expm(1j*a)
array([[ 0.42645930+1.89217551j, -2.13721484-0.97811252j],
       [ 1.06860742+0.48905626j, -1.71075555+0.91406299j]])
>>> cosm(a) + 1j*sinm(a)
array([[ 0.42645930+1.89217551j, -2.13721484-0.97811252j],
       [ 1.06860742+0.48905626j, -1.71075555+0.91406299j]])