Python SciPy linalg.fractional_matrix_power用法及代碼示例

本文簡要介紹 python 語言中 scipy.linalg.fractional_matrix_power 的用法。

用法:

scipy.linalg.fractional_matrix_power(A, t)#

計算矩陣的分數冪。

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

參數 ：：

A： (N, N) 數組

要評估其分數功效的矩陣。

t： 浮點數

分數功率。

返回 ：：

X： (N, N) 數組

矩陣的分數冪。

參考：

[1]

Nicholas J. Higham 和 Lijing lin (2011) “A Schur-Pade 矩陣的分數冪算法”。 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.]])