数值分析上机实验——高斯—赛德尔迭代求解线性方程

 

 高斯赛德尔

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

结果

 

雅克比

求个关注