matlab雅克比迭代法

以下是使用Matlab实现雅可比迭代法求解线性方程组的示例代码：

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

3x + y - z = 10
x - 5y + z = 5
2x + y + 4z = 18

我们可以将其转化为矩阵形式：

[3 1 -1; 1 -5 1; 2 1 4] * [x; y; z] = [10; 5; 18]

然后，我们可以使用Matlab编写以下代码来实现雅可比迭代法：

% 定义系数矩阵A和常数向量b
A = [3 1 -1; 1 -5 1; 2 1 4];
b = [10; 5; 18];

% 定义初始解向量x0和迭代误差限eps
x0 = [0; 0; 0];
eps = 1e-6;

% 定义迭代次数上限max_iter和当前迭代次数iter
max_iter = 1000;
iter = 0;

% 开始迭代
x = x0;
while iter < max_iter
    % 计算下一次迭代的解向量x_new
    x_new = zeros(size(x));
    for i = 1:size(A, 1)
        x_new(i) = (b(i) - A(i, 1:i-1)*x(1:i-1) - A(i, i+1:end)*x(i+1:end)) / A(i, i);
    end
    
    % 判断是否满足迭代误差限
    if norm(x_new - x) < eps
        break;
    end
    
    % 更新解向量x和迭代次数iter
    x = x_new;
    iter = iter + 1;
end

% 输出结果
disp(['Solution: x = ', num2str(x(1)), ', y = ', num2str(x(2)), ', z = ', num2str(x(3))]);

运行以上代码，即可得到线性方程组的解。