1、用matlap编写求解线性方程组的Jacobi迭代法。

好的，我可以为您提供一份 MATLAB 代码来实现 Jacobi 迭代法求解线性方程组。以下是代码：

function [x, iter] = Jacobi(A, b, x0, tol, max_iter)
% Jacobi 迭代法求解线性方程组 Ax=b
% 输入：系数矩阵 A，右端向量 b，初值 x0，容许误差 tol，最大迭代次数 max_iter
% 输出：解向量 x，迭代次数 iter

n = length(b); % 系数矩阵 A 的维数

D = diag(diag(A)); % 提取系数矩阵 A 的对角线元素
L = tril(A, -1); % 提取系数矩阵 A 的下三角元素
U = triu(A, 1); % 提取系数矩阵 A 的上三角元素
M = D; % 构造迭代矩阵 M
N = L + U; % 构造迭代矩阵 N
invM = inv(M); % 计算矩阵 M 的逆矩阵

x = x0; % 初值
iter = 0; % 迭代次数

while (iter < max_iter)
    x = invM * (b - N * x); % Jacobi 迭代公式
    if (norm(A * x - b) < tol) % 判断是否满足容许误差 tol
        return;
    end
    iter = iter + 1; % 迭代次数加 1
end

error('Jacobi 迭代法未收敛'); % 达到最大迭代次数仍未收敛，输出错误信息

您可以按照以下步骤使用 Jacobi 迭代法求解线性方程组：
1. 将上面的代码复制到 MATLAB 编辑器中，并保存为 Jacobi.m 文件。
2. 在 MATLAB 命令窗口中输入系数矩阵 A、右端向量 b、初值 x0、容许误差 tol 和最大迭代次数 max_iter，例如：

A = [4 -1 0; -1 4 -1; 0 -1 4];
b = [15; 10; 10];
x0 = [0; 0; 0];
tol = 1e-6;
max_iter = 1000;

3. 调用 Jacobi 函数求解线性方程组，例如：

[x, iter] = Jacobi(A, b, x0, tol, max_iter);

4. 得到解向量 x 和迭代次数 iter。

希望这份代码能够帮助您解决问题！