- 1、本文档共25页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
5.1 引言 凡是迭代法都有一个收敛问题,有时某种方法对一类方程组迭代收敛,而对另一类方程组进行迭代时就会发散。一个收敛的迭代法不仅程序设计简单,适于自动计算,而且较直接法更少的计算量就可获得满意的解。因此,迭代法亦是求解线性方程组,尤其是求解具有大型稀疏矩阵的线性方程组的重要方法之一。 第五章 解线性方程组的迭代法 5.2 迭代法的基本思想 迭代法的基本思想是将线性方程组转化为便于迭代的等价方程组,对任选一组初始值 ,按某种计算规则,不断地 对所得到的值进行修正,最终获得满足精度要求的方程组的近似解。 设 非奇异, ,则线性方程组 有惟一解 ,经过变换构造出一个等价同解方程组 将上式改写成迭代式 选定初始向量 ,反复不断地使用迭代式逐步逼近方程组的精确解,直到满足精度要求为止。这种方法称为迭代法 如果 存在极限 则称迭代法是收敛的,否则就是发散的。 收敛时,在迭代公式 中当 时, , 则 , 故 是方程组 的解。 对于给定的方程组可以构造各种迭代公式,是否全部收敛呢? 例5.1 用迭代法求解线性方程组 解 构造方程组的等价方程组 据此建立迭代公式 取 计算得 迭代解离精确解 越来越远,迭代不收敛。 雅可比(Jacobi)迭代法 1)雅可比迭代法算法构造 将n元线性方程组 写成 写成 若 ,分离出变量 若 ,分离出变量 有 据此建立迭代公式 上式称为解方程组的Jacobi迭代公式。 例5.2 用雅可比迭代法求解方程组 解:从方程组的三个方程中分离出 和 建立迭代公式 取初始向量 进行迭代, 可以逐步得出一个近似解的序列: (k=1, 2, …) 直到求得的近似解能达到预先要求的精度, 则迭代过程终止,以最后得到的近似解作为线 性方程组的解。 当迭代到第10次有 计算结果表明,此迭代过程收敛于方程组的精 确解x*= (3, 2, 1)T。 2) 雅可比迭代法的矩阵表示 设方程组 的系数矩阵A非奇异,且主对 角元素 ,则可将A分裂成 记作 A = L + D + U 则 等价于 即 因为 ,则 这样便得到一个迭代公式 令 则有 (k = 0,1,2…) 称为雅可比迭代公式, B称为雅可比迭代矩阵 其中 在例4.2中,由迭代公式写出雅可比迭代矩阵为 雅可比迭代矩阵表示法,主要是用来讨论其收敛性,实际计算中,要用雅可比迭代法公式的分量形式。即 (k=0,1,2,…) 雅可比迭代法的算法实现 高斯-塞德尔(Gauss-Seidel)迭代法 1 )高斯-塞德尔迭代法的基本思想 在Jacobi迭代法中,每次迭代只用到前一次的迭代值,若每次迭代充分利用当前最新的迭代值,即在求 时用新分量 代替旧分量 , 就是高斯-赛德尔迭代法。其迭代法格式为: (i=1,2,…,n k=0,1,2,…) 例5.3 用Gauss?Seidel 迭代格式解方程组 精确要求为ε=0.005 解 Gauss?Seidel 迭代格式为 取初始迭代向量 ,迭代结果为: x* ≈ 2 ) Gauss—Seidel 迭代法的矩阵表示 将A分裂成A =L+D+U,则 等价于 ( L+D+U )x = b 于是,则高斯—塞德尔迭代过程 因为 ,所以 则高斯-塞德尔迭代形式为: 故 令 3) 高斯—塞德尔迭代算法实现 高斯-塞德尔迭代算法的计算步骤与流程图与雅可比迭代法大致相同,只是一旦求出变元 的某个新值
您可能关注的文档
最近下载
- 2024山西省文化旅游投资控股集团有限公司所属企业社会招聘190人笔试模拟试题及答案解析.docx
- DB36T 941-2017补充耕地质量评定技术规范.pdf
- 3500词(带音标)——完整打印版.pdf
- 轮机管理(助理管轮)考试题复习题库(460题).docx
- 体检表---三甲医院标准体检表模板.docx VIP
- 森兰HOPE510系列高性能矢量控制变频器用户手册V1.0.pdf
- 零星及紧急采购管理办法.docx VIP
- 新人教版小学数学一年级下册第八单元 《总复习》全单元备课.docx
- 河南省南阳市新野县七年级(下)期末数学试卷.pdf VIP
- 2020年新版三体系内审检查表(全条款).doc
文档评论(0)