我知道你可以像这样逆矩阵:
。I, A**-1在python中,我想知道是否有什么不同。我问的原因是,我正在考虑使用mpmath高精度,它只支持**-1实现。
谢谢编辑:
添加一些说明:
我想知道.I
是**-1
的简写,还是调用np.linalg.inv
或其他函数。
我打赌你可以在numpy矩阵代码中找到答案:
https://github.com/numpy/numpy/blob/master/numpy/matrixlib/defmatrix.pydef matrix_power....
'''
If ``n < 0``, the inverse is computed and then raised to the ``abs(n)``
...
matrix Provides an equivalent function as the exponentiation operator (``**``, not ``^``).
'''
from numpy.linalg import inv
...
elif n<0:
M = inv(M)
n *= -1
为matrix
类
def __pow__(self, other):
return matrix_power(self, other)
M.I
的查看M.getI
(IPython
是一种方便的查看代码的方式,而不需要拉出源文件,例如M.getI??
)