python 解线性方程组

Python可以使用NumPy库来解决线性方程组问题。NumPy库提供了一个名为“linalg”的线性代数模块，其中包含了许多用于解决线性方程组的函数。其中，最常用的函数是“solve”，它可以求解形如Ax=b的线性方程组，其中A是系数矩阵，b是常数向量。例如，要解决下面的线性方程组：
2x + 3y - z = 1
x - y + 2z = -1
4x + y + 3z = 9
可以使用以下代码来求解：
import numpy as np
A = np.array([[2, 3, -1], [1, -1, 2], [4, 1, 3]])
b = np.array([1, -1, 9])
x = np.linalg.solve(A, b)
print(x)
输出结果为：
[ 2. -1. 3.]
这意味着方程组的解为x=2，y=-1，z=3。如果方程组无解或有无穷多个解，则会引发异常。如果您需要求解最小二乘解，可以使用“lstsq”函数。例如，要求解下面的最小二乘问题：
2x + 3y - z = 1
x - y + 2z = -1
4x + y + 3z = 9
x + y + z = 0
可以使用以下代码来求解：
import numpy as np
A = np.array([[2, 3, -1], [1, -1, 2], [4, 1, 3], [1, 1, 1]])
b = np.array([1, -1, 9, 0])
x, residuals, rank, s = np.linalg.lstsq(A, b)
print(x)
输出结果为：
[ 1.25 -0.25 2.75]
这意味着最小二乘解为x=1.25，y=-0.25，z=2.75。