1、求矩阵的逆
代码
#coding:utf8
import numpy as np
A=np.mat("2 3 4; 4 2 6;10 -4 18")
print "A\n",A
inverse=np.linalg.inv(A) #求矩阵的逆
print "inverse of A\n",inverse
print "Check\n",A*inverse #检验相乘是否为单位矩阵
print "Error\n",A*inverse -np.eye(3)
运行结果
A
[[ 2 3 4]
[ 4 2 6]
[10 -4 18]]
inverse of A
[[-1. 1.16666667 -0.16666667]
[ 0.2 0.06666667 -0.06666667]
[ 0.6 -0.63333333 0.13333333]]
Check
[[ 1.00000000e+00 0.00000000e+00 0.00000000e+00]
[ 4.44089