《用SOR方法求解线性方程组.docx》由会员分享，可在线阅读，更多相关《用SOR方法求解线性方程组.docx（2页珍藏版）》请在三一办公上搜索。
用SOR方法求解线性方程组第七次作业：编写用SOR方法求解线性方程组Ax=B的标准程序，并求下列方程组的解，并比较松弛因子取1.0、1.25、1.5时所需迭代的次数。 可取初始向量 代终止条件x(0)=(1,1,1)(k)T ，迭-8|x(k+1)-x|10430x12434-1x=3020-14x-243function x,k =sor(A,b,w) % A 为方程组的系数矩阵 % b 为方程组的右端项 % x 为方程组的解 % N为迭代次数 n,m=size(A);n=length(b);e=10-8;x1=ones(1,n);x=zeros(1,n);M=100;s=zeros(1,n) y=zeros(1,n); r=max(abs(b); k=0; while re for i=1:n sum=0; for j=1:n if ji sum=sum+A(i,j)*x1(j); end end s(i)=w*(b(i)-sum)/A(i,i); y(i)=(1-w)*x1(i)+s(i); end r=max(abs(s-w*x1); x1=y;k=k+1; if kM warning(不收敛); return; end end x=y; SOR方法迭代30次 SOR方法 迭代17次 w=1w=1.25结时结果如果如：图：时图 SOR方法w=1.5时结果如图迭代31次 ：
