python - 如何将 numpy.matrix 提升为非整数幂？

Question

numpy.matrix 的** 运算符不支持非整数幂:

>>> m
 matrix([[ 1. ,  0. ],
    [ 0.5,  0.5]])

>>> m ** 2.5
TypeError: exponent must be an integer

我想要的是

octave:14> [1 0; .5 .5] ^ 2.5
ans =

   1.00000   0.00000
   0.82322   0.17678

我可以使用 numpy 或 scipy 来实现吗？

备注:

这不是逐元素操作。正如 this post 中所述，它是一个矩阵(在线性代数中)的某个幂次。 .

Answer 1

你可以使用 scipy.linalg.fractional_matrix_power :

>>> m
matrix([[ 1. ,  0. ],
        [ 0.5,  0.5]])
>>> scipy.linalg.fractional_matrix_power(m, 2.5)
array([[ 1.       ,  0.       ],
       [ 0.8232233,  0.1767767]])