《SOR迭代法求解线性方程组.docx》由会员分享，可在线阅读，更多相关《SOR迭代法求解线性方程组.docx（2页珍藏版）》请在三一办公上搜索。
1、SOR迭代法求解线性方程组实验三：用SOR迭代法求解线性方程组 0.76-0.01-0.14-0.16x10.68-0.010.88-0.030.06x21.18-0.14-0.031.01-0.12x=0.12 3-0.160.06-0.120.72x0.744 取初始点x(0)=(0,0,0,0)T，松弛因子w=1.05，精度要求e=10-6。 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;
2、 return elseif nargin =5 M = 200; end if(w=2) error; return; end D=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; end end 2，输入矩阵。并根据要求调用函数，运行结果如下图所示： 即经过7次迭代算出结果，且求得： 1.27151.2844 x=0.48581.2843
