scipy.sparse.linalg.

expm#

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

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

参数:
A(M,M) array_like 或稀疏矩阵

要进行指数运算的二维数组或矩阵(稀疏或稠密)

返回:
expA(M,M) ndarray

矩阵 A 的指数

注释

这是算法 (6.1),它是算法 (5.1) 的简化版本。

Added in version 0.12.0.

参考文献

[1]

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

示例

>>> from scipy.sparse import csc_matrix
>>> from scipy.sparse.linalg import expm
>>> A = csc_matrix([[1, 0, 0], [0, 2, 0], [0, 0, 3]])
>>> A.toarray()
array([[1, 0, 0],
       [0, 2, 0],
       [0, 0, 3]], dtype=int64)
>>> Aexp = expm(A)
>>> Aexp
<Compressed Sparse Column sparse matrix of dtype 'float64'
    with 3 stored elements and shape (3, 3)>
>>> Aexp.toarray()
array([[  2.71828183,   0.        ,   0.        ],
       [  0.        ,   7.3890561 ,   0.        ],
       [  0.        ,   0.        ,  20.08553692]])