雅可比（Jacobi）迭代法求解线性方程组matlab程序设计

目录

算法步骤

 例题

 程序代码

运行结果

---

算法步骤

---

 例题

---

 程序代码

function x=Jacobidiedai(A,b,x0,eps,N)
n=length(b);
if nargin<5
    N=500;
end
if nargin<4
    eps=1e-6;
end
if nargin<3
    x0=zeros(n,1);
end
k=0;
x=zeros(n,1);
while k<N
    for i=1:n
        x(i)=(b(i)-A(i,1:i-1)*x0(1:i-1)-A(i,i+1:n)*x0(i+1:n))/A(i,i);
    end
    if norm(x-x0,inf)<eps
        break;
    end
    x0=x;k=k+1;
end
if k==N
    warning('已达到迭代次数上限');
end
fprintf('k= %d',k)
end

---