线性方程组迭代求解法.docVIP

线性方程组的迭代求解法 题目 用雅可比迭代法、高斯-塞德尔迭代法、逐次超松弛迭代法求解4元线性方程组（*）的近似解。 （*） 方法 雅可比迭代法、高斯-塞德尔迭代法、逐次超松弛迭代法 程序 Jacobi.M（雅可比迭代法的程序） function[x,k,flag]=jacobi(A,b,delta,max1) if nargin4 max1=100;end if nargin3 delta=1e-4;end n=length(A);k=0; x=zeros(n,1);y=zeros(n,1);flag=ok; while 1 for i=1:n y(i)=b(i); for j=1:n if j~=i y(i)=y(i)-A(i,j)*x(j); end end if abs(A(i,i))1e-10|k==max1 flag=fail;return; end y(i)=y(i)/A(i,i); end if norm(y-x,inf)delta break; end x=y;k=k+1; end gao.M（高斯-塞德尔迭代法的程序） function[x,k,flag]=gao(A,b,delta,max1) if nargin4 max1=100;end if nargin3 delta=1e-4;end n=length(A);k=0; x=zeros(n,1);y=zeros(n,1);flag=ok; while 1 y=x; for i=1:n z=b(i); for j=1:n if j~=i z=z-A(i,j)*x(j); end end if abs(A(i,i))1e-10|k==max1 flag=fail;return; end z=z/A(i,i); x(i)=z; end if norm(y-x,inf)delta break; end k=k+1; end sor.M（逐次超松弛迭代法的程序） function[x,k,flag]=sor(A,b,delta,w,max1) if nargin5 max1=100;end if nargin4 w=1;end if nargin3 delta=1e-4;end n=length(A);k=0; x=zeros(n,1);y=zeros(n,1);flag=ok; while 1 y=x; for i=1:n z=b(i); for j=1:n if j~=i z=z-A(i,j)*x(j); end end if abs(A(i,i))1e-10|k==max1 flag=fail;return; end z=z/A(i,i); x(i)=(1-w)*x(i)+w*z; end if norm(y-x,inf)delta break; end k=k+1; end 结果 format compact A=[2 -1 0 0;-1 2 -1 0;0 -1 2 -1;0 0 -1 2] A = 2 -1 0 0 -1 2 -1 0 0 -1 2 -1 0 0 -1 2 b=[1 0 1 0] b = 1 0 1 0 [x,k,flag]=jacobi(A,b) x = 1.1999 1.3997 1.5998 0.7998 k = 41 flag = ok [x,k,flag]=gao(A,b) x = 1.1999 1.3998 1.5999 0.7999 k =