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]
。