jacobi1.m

Gauss-seidel迭代与两种Jacobi迭代的matlab实现

%通过Jacobi迭代求解线性方程组(转化为迭代矩阵x=Bx+g求解)
function x=Jacobi1(A,b,x0,delta,max1)

%Input  - A要求解的线性方程组的系数矩阵
%       - b是线性方程组标准形式Ax=b中的b列向量
%       - x0是初始迭代值
%       - delta可以是接受的误差
%       - max1为最大迭代次数
%Output - x为所求的解

n=length(b);
for i=1:n
    b(i)=b(i)/A(i,i);
    A(i,1:n)=-A(i,1:n)/A(i,i);
    A(i,i)=0;
end
for i=1:max1
    x=A*x0+b;
    err=norm(x-x0);
    relerr=err/norm(x);
    x0=x;
    if (err<delta)|(relerr<delta)
        break;
    end
end