MATLAB代码 解线性方程组的迭代法.pdfVIP

解线性方程组的迭代法 1.rs里查森迭代法求线性方程组Ax b 的解 function [x,n] rs(A,b,x0,eps,M) if(nargin 3) eps 1.0e-6; %eps表示迭代精度 M 10000; %M表示迭代步数的限制值 elseif(nargin 4) M 10000; end I eye(size(A)); n 0; x x0; tol 1; %迭代过程 while (toleps) x (I-A)*x0+b; n n + 1; %n为最终求出解时的迭代步数 tol norm(x-x0); x0 x; if(n M) disp(Warning:迭代次数太多，可能不收敛！); return; end end 2.crs里查森参数迭代法求线性方程组Ax b 的解 function [x,n] crs(A,b