matlab雅克比迭代法算例

雅克比迭代法算例子函数：m文件

%雅克比Jacobi.m

function x=Jacobi(A,b,x0,eps,N)%用雅克比迭代法解方程

n=length(b);x=ones(n,1);k=0;

while k<=N

for i=1:n

x(i)=(b(i)-A(i,[1:i-1,i+1:n])*x0([1:i-1,i+1:n]))/A(i,i);

end

k=k+1;

if norm(x-x0,inf)

x0=x;

end

if k>N

warning('算法超出迭代的次数');

else

disp(['迭代次数= ',num2str(k)])

x

end

主程序：（调用以上程序）。M文件

clc

format long

A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];

b=[6;25;-11;15];

x0=[0;0;0;0];eps=1e-6;N=30;

x=Jacobi(A,b,x0,eps,N);