numpy求解方程组

线性代数中比较常见的问题之一是求解矩阵向量方程。 这是一个例子。 我们寻找解决方程的向量x

A x = b

当

1、用 np.linalg.solve(A,b) 直接求解 

我们首先构建A和b的数组。

A = np.array([[2,1,-2],[3,0,1],[1,1,-1]])
b = np.transpose(np.array([[-3,5,-2]])

为了解决这个系统 

x = np.linalg.solve(A,b)

形如A x = b，求解x的题目我们经常会遇到，因此Numpy专门提供函数用来解决这种问题。

2、通过求逆矩阵求解

A_inv=np.linalg.inv(A)
x=np.dot(A_inv,b)

这种方法能得到正确结果，但是要多求一步逆矩阵。