您好,欢迎访问三七文档
当前位置:首页 > 电子/通信 > 综合/其它 > SOR迭代法求解线性方程组
实验三:用SOR迭代法求解线性方程组74.012.018.168.072.012.006.016.012.001.103.014.006.003.088.001.016.014.001.076.04321xxxx取初始点Tx)0,0,0,0()0(,松弛因子05.1,精度要求610。1,建立SOR.m函数文件,此函数文件可调用,程序源码如下:function[x,n]=SOR(A,b,x0,w,eps,M)ifnargin==4eps=1.0e-6;%精度要求M=200;elseifnargin4error;returnelseifnargin==5M=200;endif(w=0||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;%迭代次数whilenorm(x-x0)=epsx0=x;x=B*x0+f;n=n+1;if(n=M)disp('Warning:迭代次数太多,可能不收敛!');return;endend2,输入矩阵。并根据要求调用函数,运行结果如下图所示:即经过7次迭代算出结果,且求得:1.27151.28440.48581.2843x
本文标题:SOR迭代法求解线性方程组
链接地址:https://www.777doc.com/doc-4721641 .html