数值分析 matlab方法 线性方程组的迭代解法.pptVIP

if n~=m error(The rows and columns of matrix A must be equal!); return; end %当方程组与右端项的维数不匹配时，停止计算，并输出出错信息. if m~=nb error(The columns of A must be equal the length of b!); return; end if nargin4 it_max=100;end if nargin3 ep=1e-5;end k=0;x=zeros(n,1);y=zeros(n,1);index=1 线性方程组迭代解法 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== it_max index=0;return; end z=z/A(i,i);x(i)=z; end 线性方程组迭代解法 if norm(y-x,inf)ep break; end k=k+1; end 调用函数Gauss_Seidel.m 解例1. 线性方程组迭代解法 x = 0.9998 k = 4 index = 1 迭代成功，收敛 0.9998