- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
第1页,共22页,星期日,2025年,2月5日§2Jacobi和Gauss-Seidel迭代法一、Jacobi迭代法设方程组将系数矩阵分裂为:其中第2页,共22页,星期日,2025年,2月5日如果原方程组可化为其中相应的迭代格式上述方法称为Jacobi迭代法,简称J法或简单迭代法分量形式:第3页,共22页,星期日,2025年,2月5日二、Gauss-Seidel迭代法G-S迭代法是J迭代法的一种改进在J迭代公式中,计算时,利用已经算出来的新的值,从而得到G-S迭代法。?G-S迭代法的分量形式:第4页,共22页,星期日,2025年,2月5日例1:利用Jacobi和Gauss-Seidel迭代法求解方程组解:Jacobi迭代格式第5页,共22页,星期日,2025年,2月5日G-S迭代格式计算结果取初值Jacobi迭代法要求精度迭代次数0.0019(1.00025071.00006941.0002507)0.000110(0.99995411.00012530.9999541)0.0000114(0.99999811.00000200.9999981)方程组的近似解第6页,共22页,星期日,2025年,2月5日?G-S迭代法的迭代矩阵:计算结果Gauss-Seidel迭代法要求精度迭代次数0.0015(0.99979160.99984791.0000664)0.00017(0.99999290.99999491.0000022)0.000018(1.00000131.00000090.9999996)方程组的近似解取初值由迭代公式迭代矩阵第7页,共22页,星期日,2025年,2月5日三、Jacobi和Gauss-Seidel迭代法的收敛性Jacobi迭代法收敛的充要条件是Gauss-Seidel迭代法收敛的充要条件是推论1:Jacobi迭代法收敛的充分条件是Gauss-Seidel迭代法收敛的充分条件是如例1:利用J和G-S迭代法求解方程组第8页,共22页,星期日,2025年,2月5日Jacobi迭代矩阵系数矩阵第9页,共22页,星期日,2025年,2月5日Gauss-Seidel迭代矩阵第10页,共22页,星期日,2025年,2月5日设满足称为严格对角占优矩阵如果且至少有一个严格不等式成立,则称为弱对角占优矩阵。设,如果能找到排列阵,使得其中与均为方阵,称为可约的否则称为不可约的第11页,共22页,星期日,2025年,2月5日例如:矩阵是可约的若系数矩阵是可约的,则可通过行与列重排化为(*)式,从而可以将方程组简化为低阶方程组。第12页,共22页,星期日,2025年,2月5日(补充:可约矩阵的等价定义)是可约矩阵,当且仅当存在一个下标的非空子集,使得例如:矩阵矩阵不可约第13页,共22页,星期日,2025年,2月5日如果严格对角占优,则,且非奇异。如果不可约且弱对角占优,则,且非奇异。自己看证明:?首先证明设由条件:是弱对角占优,交换的第k、n行与k、n列,则矩阵变为与不可约矛盾!第14页,共22页,星期日,2025年,2月5日?其次证明是非奇异的设则存在非零向量满足定义下标的集合且令对某个j显然J非空,否则第15页,共22页,星期日,2025年,2月5日对,有由此可知,当时,但对于都有所以否则与弱对角占优矛盾!与不可约矛盾第16页,共22页,星期日,2025年,2月5日因为
文档评论(0)