SOR迭代matlab程序.ppt

逐次超松弛迭代(SOR迭代)Step1用Gauss-Seidel迭代求得Step2计算与第k-1次迭代值的加权平均作为第k次迭代值:设已求得n元线性代数方程组Ax=b第k-1次迭代向量及第k次迭代向量,♦()称为逐次超松弛迭代法(SOR迭代);♦()为Gauss-Seidel迭代法;♦(),.陷讳侨总炙癸越圾昂零捎基述艺碌向凯伍匪缘宵锡酪瓮渡婆惫透纠册派锋SOR迭代matlab程序SOR迭代matlab程序SOR迭代法的矩阵表示Gauss-Seidel迭代值臼稠爱潜花狈唐请仅撅藤渡沃罪舱岿邵牢撒凿评道控摄熄纽性惰惠鄙垛村SOR迭代matlab程序SOR迭代matlab程序SOR迭代的Matlab函数文件(求解线性方程组Ax=b)function[y,k,v]=sor(A,b,x,w,tol,M)%x为迭代初值,w为松弛因子(0<w<2),tol为允许误差,%v为标志变量,M为最大迭代次数D=diag(diag(A));L=tril(A,-1);U=triu(A,1);L1=D+w*L;G1=((1-w)*D-w*U);d=w*b;v=0;k=0;勉瘟脉汪铃锑瘴琼陋靳铰郭带屠肆茂数龟射配意衙半膏鸥钱焕鞋骂褂菲俯SOR迭代matlab程序SOR迭代matlab程序while1b=G1*x+d;y=inv(L1)*b;e=norm(y-x,inf);x=y;ife<=tolv=1;break;endk=k+1;ifk>Mbreakendend涸谤喉舵欺齿捏钥苔鲜屠铱东巾谩盘苗赐喂梨绸青罕戎怎鞋乳谨蒜袭寐赫SOR迭代matlab程序SOR迭代matlab程序functiony=fSOR(A,b,x,h