scipy.linalg.logm

scipy.linalg.logm(A, disp=True)[源代码]

计算矩阵对数。

矩阵对数是exPM:exPM(Logm)的逆 (A ))== A

参数
A(n,N)类数组

要计算其对数的矩阵

disp布尔值,可选

如果结果中的误差估计较大,则打印警告,而不是返回估计误差。(默认值:TRUE)

退货
logm(N,N)ndarray

的矩阵对数 A

errest浮动

(如果disp==false)

1-估计误差的范数, ||err||_ 1/ ||A||_ 1

参考文献

1

Awad H.Al-Mohy和Nicholas J.Higham(2012)“改进了矩阵对数的逆缩放和平方算法。暹罗科学计算学报,34(4)。C152-C169。ISSN 1095-7197

2

尼古拉斯·J·希格汉姆(2008年)“矩阵函数:理论与计算”ISBN978-0-898716-46-7

3

尼古拉斯·J·希格姆和林丽静(2011),“求矩阵分数次幂的Schur-Pade算法”。“暹罗矩阵分析与应用学报”,32(3)。1056-1078页。ISSN 0895-4798

示例

>>> from scipy.linalg import logm, expm
>>> a = np.array([[1.0, 3.0], [1.0, 4.0]])
>>> b = logm(a)
>>> b
array([[-1.02571087,  2.05142174],
       [ 0.68380725,  1.02571087]])
>>> expm(b)         # Verify expm(logm(a)) returns a
array([[ 1.,  3.],
       [ 1.,  4.]])