scipy.linalg.

fractional_matrix_power#

scipy.linalg.fractional_matrix_power(A, t)[source]#

计算矩阵的分数幂。

按照 [1] 第 (6) 章中的讨论进行。

参数 :

A(N, N) array_like: 用于计算分数幂的矩阵。
tfloat: 分数幂。

返回 :

X(N, N) array_like: 矩阵的分数幂。

引用

[1]

Nicholas J. Higham 和 Lijing lin（2011）“矩阵分数幂的舒尔-帕德算法”。SIAM 矩阵分析与应用杂志，第 32 卷 (3)。第 1056-1078 页。ISSN 0895-4798

示例

>>> import numpy as np
>>> from scipy.linalg import fractional_matrix_power
>>> a = np.array([[1.0, 3.0], [1.0, 4.0]])
>>> b = fractional_matrix_power(a, 0.5)
>>> b
array([[ 0.75592895,  1.13389342],
       [ 0.37796447,  1.88982237]])
>>> np.dot(b, b)      # Verify square root
array([[ 1.,  3.],
       [ 1.,  4.]])