SOR迭代法求解线性方程组(共2页).doc

《SOR迭代法求解线性方程组(共2页).doc》由会员分享，可在线阅读，更多相关《SOR迭代法求解线性方程组(共2页).doc（2页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

精选优质文档-倾情为你奉上实验三：用SOR迭代法求解线性方程组 取初始点，松弛因子，精度要求。1，建立SOR.m函数文件，此函数文件可调用，程序源码如下：function x,n=SOR(A,b,x0,w,eps,M)if nargin=4 eps= 1.0e-6;%精度要求 M = 200;elseif nargin4 error; returnelseif nargin =5 M = 200;end if(w=2) error; return;endD=diag(diag(A); %求A的对角矩阵L=-tril(A,-1); %求A的下三角阵U=-triu(A,1); %求A的上三角阵B=inv(D-L*w)*(1-w)*D+w*U);f=w*inv(D-L*w)*b;x=B*x0+f;n=1; %迭代次数while norm(x-x0)=eps x0=x; x =B*x0+f; n=n+1; if(n=M) disp(Warning: 迭代次数太多，可能不收敛！); return; endend 2，输入矩阵。并根据要求调用函数，运行结果如下图所示：即经过7次迭代算出结果，且求得：专心-专注-专业