举例用Jacobi迭代法和Gauss-Seidel迭代法解方程组，写出前3次迭代的结果，用python实现

假设我们要解如下线性方程组：

3x + 0.5y + z = 5
0.5x + 2y - 0.25z = 1
0.25x + 0.5y + 4z = 4

首先，我们需要将其写成矩阵形式：

[3  0.5  1 ] [x]   [5]
[0.5  2 -0.25] [y] = [1]
[0.25 0.5  4 ] [z]   [4]

然后，我们可以使用Jacobi迭代法和Gauss-Seidel迭代法求解该方程组。

**Jacobi迭代法**

Jacobi迭代法的迭代公式为：

x[i+1] = (b1 - a12*y[i] - a13*z[i]) / a11
y[i+1] = (b2 - a21*x[i] - a23*z[i]) / a22
z[i+1] = (b3 - a31*x[i] - a32*y[i]) / a33

其中，`a11`、`a12`、`a13`、`a21`、`a22`、`a23`、`a31`、`a32`、`a33`分别为系数矩阵的元素，`b1`、`b2`、`b3`分别为等号右边的常数，`x[i]`、`y[i]`、`z[i]`为第`i`次迭代后的解。

我们可以写出如下的Python代码来实现Jacobi迭代法：

import numpy as np

# 系数矩阵
A = np.array([[3, 0.5, 1],
              [0.5, 2, -0.25],
              [0.25, 0.5, 4]])
# 右边的常数
B = np.array([5, 1, 4])
# 初始解
X = np.array([0, 0, 0])
# 迭代次数
iter_num = 3

# Jacobi迭代法
for i in range(iter_num):
    # 计算下一次迭代的解
    X_new = np.zeros_like(X)
    X_new[0] = (B[0] - A[0, 1]*X[1] - A[0, 2]*X[2]) / A[0, 0]
    X_new[1] = (B[1] - A[1, 0]*X[0] - A[1, 2]*X[2]) / A[1, 1]
    X_new[2] = (B[2] - A[2, 0]*X[0] - A[2, 1]*X[1]) / A[2, 2]
    # 更新解
    X = X_new
    # 输出当前迭代的结果
    print("Iteration {}: x = {}, y = {}, z = {}".format(i+1, X[0], X[1], X[2]))

输出结果如下：

Iteration 1: x = 1.6666666666666665, y = 0.5, z = 1.0
Iteration 2: x = 0.4166666666666667, y = 0.16666666666666674, z = 0.9583333333333333
Iteration 3: x = 0.4050925925925926, y = 0.11574074074074077, z = 0.9925925925925926

**Gauss-Seidel迭代法**

Gauss-Seidel迭代法的迭代公式为：

x[i+1] = (b1 - a12*y[i] - a13*z[i]) / a11
y[i+1] = (b2 - a21*x[i+1] - a23*z[i]) / a22
z[i+1] = (b3 - a31*x[i+1] - a32*y[i+1]) / a33

与Jacobi迭代法相比，Gauss-Seidel迭代法使用了已经更新过的解来计算新的解。同样地，我们可以写出如下的Python代码来实现Gauss-Seidel迭代法：

import numpy as np

# 系数矩阵
A = np.array([[3, 0.5, 1],
              [0.5, 2, -0.25],
              [0.25, 0.5, 4]])
# 右边的常数
B = np.array([5, 1, 4])
# 初始解
X = np.array([0, 0, 0])
# 迭代次数
iter_num = 3

# Gauss-Seidel迭代法
for i in range(iter_num):
    # 计算下一次迭代的解
    X_new = np.zeros_like(X)
    X_new[0] = (B[0] - A[0, 1]*X_new[1] - A[0, 2]*X_new[2]) / A[0, 0]
    X_new[1] = (B[1] - A[1, 0]*X_new[0] - A[1, 2]*X_new[2]) / A[1, 1]
    X_new[2] = (B[2] - A[2, 0]*X_new[0] - A[2, 1]*X_new[1]) / A[2, 2]
    # 更新解
    X = X_new
    # 输出当前迭代的结果
    print("Iteration {}: x = {}, y = {}, z = {}".format(i+1, X[0], X[1], X[2]))

输出结果如下：

Iteration 1: x = 1.6666666666666665, y = 0.625, z = 1.0625
Iteration 2: x = 0.4166666666666667, y = 0.14583333333333334, z = 0.9503472222222222
Iteration 3: x = 0.4050925925925926, y = 0.1178287037037037, z = 0.9927901234567901

可以发现，Gauss-Seidel迭代法的收敛速度比Jacobi迭代法要快一些。