scipy.linalg.

fractional_matrix_power#

scipy.linalg.fractional_matrix_power(A, t)[原始碼]#

計算矩陣的分數次方。

根據 [1] 的第 (6) 節中的討論進行。

參數:

A(N, N) array_like: 要評估分數次方的矩陣。
tfloat: 分數次方。

回傳:

X(N, N) array_like: 矩陣的分數次方。

參考文獻

[1]

Nicholas J. Higham and Lijing lin (2011) “A Schur-Pade Algorithm for Fractional Powers of a Matrix.” SIAM Journal on Matrix Analysis and Applications, 32 (3). pp. 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.]])