python 求解多元函数方程编码

目录

方法一

方法二

---

方法一

矩阵的方法

此方法只适用于多云系数分离的方程

 一般地，我们设解线性方程组形如 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)