在SciPy中,可以使用 scipy.linalg.solve
函数来实现高斯消除法(也称为高斯消元法)解线性方程组。以下是一个简单的示例:
import numpy as np
from scipy.linalg import solve
# 定义系数矩阵和常数矢量
A = np.array([[2, 1, -1], [4, 1, 3], [1, -1, 2]])
b = np.array([8, 10, 3])
# 使用 solve 函数求解线性方程组
x = solve(A, b)
print("解为:", x)
在这个例子中,我们先定义了一个系数矩阵 A
和一个常数矢量 b
,然后使用 solve
函数来解线性方程组 Ax = b。最后打印出解向量 x
。
通过这种方法,可以很容易地使用SciPy中的高斯消除法来解决线性方程组。