import numpy as np
 
a = np.array([[1, 2], [3, 4]]) # 初始化一个非奇异矩阵(数组)
print(np.linalg.inv(a)) # 对应于MATLAB中 inv() 函数
 
# 矩阵对象可以通过 .I 更方便的求逆
A = np.matrix(a)
print(A.I)

import numpy as np
 
# 定义一个奇异阵 A
A = np.zeros((4, 4))
A[0, -1] = 1
A[-1, 0] = -1
A = np.matrix(A)
print(A)
# print(A.I) 将报错，矩阵 A 为奇异矩阵，不可逆
print(np.linalg.pinv(a))  # 求矩阵 A 的伪逆（广义逆矩阵），对应于MATLAB中 pinv() 函数