GaussSeidel迭代法的Matlab程序.doc

《GaussSeidel迭代法的Matlab程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、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;   %迭代过程   whilek<=N       x(1)=(b(1)-A(1,2:n)*x0(2:n)')/A(1,1);       fori=2:n         

2、  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       ifmax(abs(x-x0))<=tol           fid=fopen('G_S_iter_result.txt','wt');           fprintf(fid,'********用Gauss_Seidel迭代求解线性方程组的输出结果********');           fprintf(fid,'迭代次数:%d次',k);           f

3、printf(fid,'x的值');           fprintf(fid,'%12.8f',x);           break;       end       k=k+1;       x0=x;   end   ifk==N+1     fid=fopen('G_S_iter_result.txt','wt');     fprintf(fid,'********用Gauss_Seidel迭代求解线性方程组的输出结果********');     fprintf(fid,'迭代次数:%d次',k);     fp

4、rintf(fid,'超过最大迭代次数，求解失败！');     fclose(fid);   end