jacobi.asv

来自「Numerical analysis （Gauss、Bisection、jaco」· ASV 代码 · 共 24 行

function s=jacobi(a,b,x0,eps,N)
%Jacobi迭代法解线性方程组
%a为系数矩阵，b为方程组ax=b的右端项，x0为初值
if nargin==4
    eps=1.0e-6;
elseif nargin<4
%    error
    return
end
D=diag(diag(a));%求对角矩阵
D=inv(D);%求对角矩阵的逆
L=tril(a,-1);%求严格的下三角矩阵
U=triu(a,1);%求严格的上三角矩阵
B=-D*(L+U);
f=D*b;
s=B*x0+f;
i = 0;
while ((norm(s-x0)>=eps)&&(i <= N))
    x0=s;
    s=B*x0+f;
    i = i+1
end
return