原创力文档- 1、本文档共22页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
*第四章解线性方程组的迭代法§1迭代法的构造设A?Rn?n，用迭代法求解线性方程组 Ax=b 首先是构造迭代序列，通常可找到等价的方程组 x=Bx+f 从而可构造迭代法 x(k+1)=Bx(k)+fk=0,1,?B?Rn?n称为迭代矩阵，当给出初始近似x(0)?Rn可逐次求出x(k+1)只依赖于前一步x(k)的值，故称为单步定常迭代。若{x(k)}收敛于x*，则x*即为解。定义2若迭代公式产生的序列{x(k)}满足?x(0)?Rn则称迭代法是收敛的。为构造迭代法，可将A分裂为 A=M–N 其中M非奇异，且容易求逆，则可得 x=M-1Nx+M-1b 若令 B=M-1N=I–M-1A f=M-1b 便可得到x=Bx+f的形式，从而可构造迭代法。通常为使M容易求逆，可取M为对角阵、三角阵或三对角阵等不同形式，可构造不同的迭代法。常将方程A=(aij)?Rn?n分裂为 A=D–L–U 其中D=diag(a11,a22,…,ann)为A的对角部分，若aii?0(i=1,…,n)，取M=D,N=L+U则得x(k+1)=BJx(k)+f k=0,1,… BJ=D-1(L+U)=I-D-1A f=D-1b称为Jacobi迭代法。若取M=D–L，则得Gauss-Seidel迭代法: x(k+1)=BGx(k)+f 其中:BG=(D–L)-1U=I–(D–L)-1A,f=(D-L)-1b按Jacobi迭代法与G-S迭代法计算时可改写为分量形式，分别为:(i=1,…,n)k=0,1,…(i=1,…,n)k=0,1,…两者之差别仅在于已得分量用新值代替旧值。§3Jacobi法和Gauss-Seidel法BJacobi迭代法Gauss-Seidel迭代法…………写成矩阵形式：B*
文档评论（0）