One - One Code All

Blog Content

python科学计算之线性方程组求解sympy

Python   2013-07-13 19:27:21

解方程可以使用scipy库,也可以使用sympy库。


闲话少叙,直接上代码:

'''
# 方程组
3a + b - 2c = 1.8
a - b + 0.9c = -2.2
0.8a + 1.3c = 1.5
'''

import numpy as np
from scipy.linalg import solve
a = np.array([[3, 1, -2], [1, -1, 0.9], [0.8, 0, 1.3]])
b = np.array([1.8, -2.2, 1.5])
x = solve(a, b)
print(x)

输出结果:

[0.18585526 3.32138158 1.03947368]


例子2:

'''
# 方程组
2 * x - y = 8
3 * x + y = 28
'''

from sympy import *
x = Symbol('x')
y = Symbol('y')
print(solve([2 * x - y - 8, 3 * x + y - 28],[x, y]))

输出结果:

{x: 36/5, y: 32/5}


上一篇:pandas乘法DataFrame.multiply()
下一篇:python科学计算之非线性方程组求解scipy

The minute you think of giving up, think of the reason why you held on so long.