线性代数中比较常见的问题之一是求解矩阵向量方程。 这是一个例子。 我们寻找解决方程的向量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)
这种方法能得到正确结果,但是要多求一步逆矩阵。