目录
方法一
矩阵的方法
此方法只适用于多云系数分离的方程
一般地,我们设解线性方程组形如 Ax=b,其中 A 是系数矩阵,b 是一维(n 维也可以,这个下面会提到),x 是未知变量. 再拿上面地最简单的二元一次方程组为例,我们用 numpy.linalg.solve 可以这样写:
import numpy as np
A = np.mat('1,2; 4,5') # 构造系数矩阵 A
b = np.mat('3,6').T # 构造转置矩阵 b (这里必须为列向量)
r = np.linalg.solve(A,b) # 调用 solve 函数求解
print (r)
方法二
SymPy 通吃一切
先设置变量、然后写方程即可
x = Symbol('x')
y = Symbol('y')
z = Symbol('z')
r = solve([
(x + y) / 2 - (y + z) / 3,
(x + y) / 2 - (x + z) / 4,
(y + z) / 3 - (x + z) / 4,
x + 2 * y + z - 12
], [x, y, z])
print(r)