Jacobi迭代算法

%雅克比迭代法：建立函数文件Jacobi.m
function Jacobi（AbXnmw_error）
%用雅克比迭代法求解方程组Ax=b
%输入：A为方程组的系数矩阵，b为方程组右端的列向量，X为迭代初值构成的列向量，nm为最大迭代次数，w_error为误差精度
%输出：x为求得的方程组的解构成的列向量，n为迭代次数
n=1;
m=length（A）;
D=diag（diag（A））;        %令A=D-L-U计算矩阵D
L=tril（-A）+D;           %令A=D-L-U计算矩阵L
U=triu（-A）+D;          %令A=D-L-U计算矩阵U
M=inv（D）*（L+U）;       %计算迭代矩阵
g=inv（D）*b;            %计算迭代格式中的常数项
%