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]])