Jacobi迭代求解MATLAB

进过18次Jacobi迭代后，其相邻迭代解间无穷范数误差小于：1.0e-8 此时Jacobi迭代解如下： x = 1.099999996412137 1.199999996412137 1.299999995744652

%*********ZIBAI*********
%Jacobi迭代求解

clear all;
format long;

%设置方程参数
A = [10-1-2;-110-2;-1-15];
b = [7.2;8.3;4.2];
x0 = [1.1;1.2;1.3];
%设置精度，最大迭代次数，
eps = 1.0e-8;
M = 200;

%Jacobi迭代求解
D = diag（diag（A））%求A的对角矩阵
L = -tril（A-1）%求A的下三角矩阵
U = -triu（A1）%求A的上三角矩阵
B = D\（L+U）;
f = D\b;
x00 = [0;0;0];
x = B*x00+f;
n = 1;%迭代次数
while norm（（x-x0） ‘inf‘）>=eps
    xk = x;
    x = B*xk+f;
    n = n+1;
    if（n>=M）
        disp（‘Warning:迭代次数太多可能不收敛!‘）
        break;
    end
end

fprintf （‘%d次迭代后，相邻无穷范数误差小于：1.0e-8\n‘ n）
fprintf （‘迭代解：‘）;x

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-10-12 20:12  Jacobi迭代\
     文件         658  2018-10-12 20:11  Jacobi迭代\Jacobi.m