jacobi2.m

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

%通过Jacobi迭代求解线性方程组(直接求解x的各个分量)
function x=Jacobi2(A,b,x0,delta,max1)

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

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