阅读量:8
SciPy库中提供了很多求解线性方程组的方法,其中最常用的是使用scipy.linalg.solve
函数。下面是一个简单的例子来展示如何使用SciPy求解线性方程组。
假设我们有一个线性方程组:
3x + 2y - z = 1 2x - 2y + 4z = -2 -x + 0.5y - z = 0
我们可以将其表示成矩阵形式 Ax = b,其中
A = [[3, 2, -1], [2, -2, 4], [-1, 0.5, -1]] b = [1, -2, 0]
接下来,我们可以使用scipy.linalg.solve
函数来求解该线性方程组:
import numpy as np from scipy.linalg import solve A = np.array([[3, 2, -1], [2, -2, 4], [-1, 0.5, -1]]) b = np.array([1, -2, 0]) x = solve(A, b) print(x)
运行以上代码,我们可以得到方程组的解为 x = [1, -2, -2]
。