jocabi_1.m

来自「本代码是本人在做作业时所编」· M 代码 · 共 33 行

function x=jacobi_1(a,b)
a=[10 1 2 3 4;1 9 -1 2 -3;2 -1 7 3 -5;3 2 3 12 1;4 -3 -5 -1 15];
b=[12 -27 14 -17 12];
tol=1e-6;
x0=[0;0;0;0;0];
max=100;
[n,m]=size(a);
xold=x0;
B=-a;
for i=1:n
B(i,i)=0;
B(i,:)=B(i,:)/a(i,i);
g(i,1)=b(i)/a(i,i);
end
i=1;
while i<max;
xnew=B*xold+g;
if norm(xnew-xold)<=tol
   x=xnew;
   disp('迭代收敛');
   return;
else
   xold=xnew;
  end
  disp([xnew']);
  i=i+1;
end
disp('不收敛');
disp('jiguo')
x=xnew;