1. 待解方程式
已知ai,bi,ci,求解x,y,z。
a1*x + b1*y+c1*z = 1
a2*x + b2*y+c2*z = 2
a3*x + b3*y+c3*z = 3
2. 需要的模块与方法
numpy
np.linalg.solve(m, n)
3. 输入格式
其中 np.linalg.solve(m, n)
的m,n输入格式为:
m = np.array([[a1, b1, c1], [a2, b2. c2], [a3, b3, c3]])
n = np.array([1, 2, 3]) #可替换为式子右边的常数
4. 结果
solution = np.linalg.solve(m, n) #solution format: np.array([x, y, z])
完整伪代码如下:
import numpy as np
"""
a1*x + b1*y+c1*z = 1
a2*x + b2*y+c2*z = 2
a3*x + b3*y+c3*z = 3
"""
m = np.array([[a1, b1, c1], [a2, b2. c2], [a3, b3, c3]])
n = np.array([0, 0, 0]) #可替换为式子右边的常数
solution = np.linalg.solve(m, n) #solution format: np.array([x, y, z])
注意: np.linalg.solve
前提是m矩阵可逆。