线性代数方程组的迭代解法详解演示文稿_第1页
线性代数方程组的迭代解法详解演示文稿_第2页
线性代数方程组的迭代解法详解演示文稿_第3页
线性代数方程组的迭代解法详解演示文稿_第4页
线性代数方程组的迭代解法详解演示文稿_第5页
已阅读5页,还剩15页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

线性代数方程组的迭代解法详解演示文稿目前一页\总数二十页\编于八点(优选)线性代数方程组的迭代解法目前二页\总数二十页\编于八点如果原方程组可化为其中相应的迭代格式上述方法称为Jacobi迭代法,简称J法或简单迭代法分量形式:目前三页\总数二十页\编于八点二、Gauss-Seidel迭代法G-S迭代法是J迭代法的一种改进在J迭代公式中,计算时,利用已经算出来的新的值,从而得到G-S迭代法。G-S迭代法的分量形式:目前四页\总数二十页\编于八点例1:利用Jacobi和Gauss-Seidel迭代法求解方程组解:Jacobi迭代格式目前五页\总数二十页\编于八点G-S迭代格式计算结果取初值Jacobi迭代法要求精度迭代次数

0.0019(1.00025071.00006941.0002507)0.000110(0.99995411.00012530.9999541)0.0000114(0.99999811.00000200.9999981)方程组的近似解目前六页\总数二十页\编于八点G-S迭代法的迭代矩阵:计算结果Gauss-Seidel迭代法要求精度迭代次数

0.0015(0.99979160.99984791.0000664)0.00017(0.99999290.99999491.0000022)0.000018(1.00000131.00000090.9999996)方程组的近似解取初值由迭代公式迭代矩阵目前七页\总数二十页\编于八点三、Jacobi和Gauss-Seidel迭代法的收敛性Jacobi迭代法收敛的充要条件是Gauss-Seidel迭代法收敛的充要条件是推论1:Jacobi迭代法收敛的充分条件是Gauss-Seidel迭代法收敛的充分条件是如例1:利用J和G-S迭代法求解方程组目前八页\总数二十页\编于八点Jacobi迭代矩阵系数矩阵目前九页\总数二十页\编于八点Gauss-Seidel迭代矩阵目前十页\总数二十页\编于八点设满足称为严格对角占优矩阵如果且至少有一个严格不等式成立,则称为弱对角占优矩阵。设,如果能找到排列阵,使得其中与均为方阵,称为可约的否则称为不可约的目前十一页\总数二十页\编于八点例如:矩阵是可约的若系数矩阵是可约的,则可通过行与列重排化为(*)式,从而可以将方程组简化为低阶方程组。目前十二页\总数二十页\编于八点(补充:可约矩阵的等价定义)是可约矩阵,当且仅当存在一个下标的非空子集,使得例如:矩阵矩阵不可约目前十三页\总数二十页\编于八点如果严格对角占优,则,且非奇异。如果不可约且弱对角占优,则,且非奇异。自己看证明:首先证明设由条件:是弱对角占优,交换的第k、n行与k、n列,则矩阵变为与不可约矛盾!目前十四页\总数二十页\编于八点其次证明是非奇异的设则存在非零向量满足定义下标的集合且令对某个j显然J非空,否则目前十五页\总数二十页\编于八点对,有由此可知,当时,但对于都有所以否则与弱对角占优矛盾!与不可约矛盾目前十六页\总数二十页\编于八点如果为严格对角占优或为不可约且弱对角占优矩阵,则求解方程组的J法和G-S法均收敛。证明:仅给出不可约且弱对角占优矩阵G-S法的证明只要证明,其中设有一个特征值,满足,且有是不可约且弱对角占优矩阵,由定理6.8:目前十七页\总数二十页\编于八点因此注意到和的零元素和非零元素的位置完全一样,故是不可约也是弱对角占优矩阵矛盾!如果为严格对角占优矩阵,易证其中为J法的迭代矩阵目前十八页\总数二十页\编于八点如果是对称矩阵,且有正的对角元,则求解方程组的J法收敛的充要条件是矩阵和均为正定的,其中证明:记其中迭代矩阵矩阵和相似,故有相同的特征值;且、、对称目前十九页\总数二十页\编于八点必要性设J法收敛,则记

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论