高斯赛德尔
function x = Gauss(a,b,e)
a[4,-1,0,-1,0,0;-1,4,-1,0,-1,0;0,-1,4,-1,0,-1;-1,0,-1,4,-1,0;0,-1,0,-1,4,-1;0,0,-1,0,-1,4];
b =[0;5;-2;5;-2;6];
e=0.0001;
n=length(b);
N=100;
x0=zeros(n,1);
x=x0;
x0=x+2*e;
k=0;
a1=tril(a);
a2=inv(a1);
while norm(x0-x)>=e&&k<N
k=k+1;
x0=x;
x=-a2*(a-a1)*x0+a2*b;
end
if k==N
warning('已达最大迭代次数');
disp('迭代次数已到达上限!');
end
disp(['迭代次数 k=',num2str(k)])
end
结果
雅克比迭代法
function x=jacobi_fun(a,b,n,x0,tol,N)
x0=zeros(6,1);
k=0;
n=6;
a=[4,-1,0,-1,0,0;-1,4,-1,0,-1,0;0,-1,4,-1,0,-1;-1,0,-1,4,-1,0;0,-1,0,-1,4,-1;0,0,-1,0,-1,4];
b =[0;5;-2;5;-2;6];
tol=0.0001;
N=500;
while k<N
k=k+1;
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
if norm(x-x0)<=tol
break;
end
x0=x;
disp(['when k=',num2str(k)])
disp('x=');
disp(x); %输出计算的中间结果
end
if k==N
disp('迭代次数已到达上限!');
end
disp(['迭代次数 k=',num2str(k)])
end
法2
function x=Jacobi_iteration(n,a,b,tol,N,x0)
x=zeros(6,1);
k=0;
a=[4,-1,0,-1,0,0;-1,4,-1,0,-1,0;0,-1,4,-1,0,-1;-1,0,-1,4,-1,0;0,-1,0,-1,4,-1;0,0,-1,0,-1,4];
b =[0;5;-2;5;-2;6];
x0=zeros (6,1);
n=6;
tol=0.0001;
N=500;
while k<N
k=k+1;
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
if norm(x-x0)<tol
break;
end
%未达到精度,继续进行运算
x0=x;
disp(['迭代次数 k=',num2str(k)])
disp('x=');
disp(x);
end
if k==N
disp('迭代次数已经达到上限!');
end
disp(['迭代次数 k=',num2str(k)])
end
结果
![](https://img-blog.csdnimg.cn/direct/bb8510758e0844bd9465bb810795f046.png)
求个关注