Gauss-Seidel迭代法的Matlab程序-(共2页).doc

Gauss_Seidel迭代法的Matlab程序function x=Gauss_Seidel_iterative(A,b)% 用Gauss_Seidel迭代求解线性方程组,矩阵A是方阵 x0=zeros(1,length(b); % 赋初值 tol=10(-2); % 给定误差界 N=1000; % 给定最大迭代次数 n,n=size(A); % 确定矩阵A的阶 k=1; % 迭代过程 while k=N x(1)=(b(1)-A(1,2:n)*x0(2:n)/A(1,1); for i=2:n x(i)=(b(i)-A(i,1:i-1)*x(1:i-1)-A(i,i+1:n)*x0(i+1:n)./A(i,i); end if max(abs(x-x0)=tol fid = fopen(G_S_iter_result.txt, wt); fp