numpy - 矩阵基本操作

# -*- coding: utf-8 -*-

import numpy as np

# 二维随机数组
array = np.random.rand(4, 4)
print array

# 二维数组转成矩阵
randmat = np.mat(array)
print randmat

# 矩阵求逆
invrandmat = randmat.I
print invrandmat

# 创建4*4的单位矩阵
eye4 = np.eye(4)
print eye4

# 矩阵乘法
myeye = randmat * invrandmat
# 由于误差原因， 得到的可能并不是4*4的单位矩阵
print myeye

# 矩阵减法，误差大小
print myeye - eye4

print eye4 + eye4