基本迭代方法课件.ppt

《基本迭代方法课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、5.1基本迭代方法5.1.2Jacobi迭代法和Gauss-Seidel迭代法5.1.1迭代公式的构造第五章线性方程组的迭代解法教学目的1.掌握Jacobi迭代法，G-S迭代法解大型线性方程组的方法及其收敛性的判别方法；2.掌握SOR迭代法及收敛的必要条件(0<ω<2)；3.了解三种迭代法之间的改进关系从而掌握该思想方法；4.理解迭代法基本定理。教学重点及难点重点是三种迭代法及收敛性的判别方法；难点是迭代法基本定理及三种迭代法收敛定理的证明。第5章线性方程组的迭代解法首先看一个形成大型方程组的例子。考虑下面的Poisson方程的离散逼近，其边界条件为：取进行网格剖分,用二

2、阶导数,按逐行自左至右和自下而上的自然次序离散华可得下列线性方程组其中　　　　　　　　　　　是　　　　　　　　　的近似值。这是一种特殊形状的稀疏矩阵。随着　　和　　的减少，所得到的方程组的阶数将增大。对于大型线形代数方程组，常用迭代解法。它是从某些初始向量出发，用设计好的步骤逐次算出近似解向量　　，从而得到向量序列。一般　　　的计算公式是称之为多步迭代法．若　　　只与　　有关，且　　是线性的，即其中，称为单步线性迭代法，称为迭代距阵。若和都与k无关，即称为单步定常线性迭代法。本章主要讨论具有这种形式的各种迭代方法。5.1基本迭代方法5.1.1迭代公式的构造设，，A非奇异，

3、满足方程组Ax=b。（5.1.1）如果能找到距阵，向量，使可逆，而且方程组x=Bx+f(5.1.2)的唯一解就是方程组（5.1.1）的解，则可从（5.1.2）式构造一个定常的线性迭代公式（5.1.3）给定初始向量,由(5.1.3)可以产生序列,若它有极限,显然就是(5.1.1)和(5.1.2)的解。定义5.1若对任意初始向量，迭代公式（5.1.3）产生的序列都有则称迭代法（5.1.3）是收敛的。从（5.1.1）出发，可以由不同的途径得到各种不同的等价方程组（5.1.2），从而得到不同的迭代法（5.1.3）。例如，设A可以分解为，其中M非奇异，则由（5.1.1）可得令就可以

4、得到（5.1.2）的形式。不同的分解方式，可的不同的B和f,下面给出对应不同分解方式的常用迭代计算公式。5.1.2Jacobi迭代法和Gauss-Seidel迭代法1.Jacobi迭代法记，可以把A分解为（5.1.4）其中现设D非奇异，即。方程组（5.1.1）等价于用J法计算向量序列，要用两组单元存放向量和。迭代法可以写成分量形式（5.1.8）由此构造迭代公式：（5.1.5）其中迭代距阵和向量为（5.1.6）（5.1.7）称（5.1.5）为解（5.1.1）的Jacobi迭代法，简称J法。2.Gauss-Seidel迭代法在J法中，计算时，分量已经算出，所以可考虑对J法进行

5、修改。在每个分量计算出来之后，下一个分量的计算就利用最新的计算结果。这样，在整个迭代过程中只要使用一组单元存放迭代向量，其分量形式的计算结果为（5.1.9）这就是Gauss-Seidel迭代法，简称GS法将（5.1.9）写成距阵形式经整理有（5.1.10）其中迭代距阵和向量为(5.1.11)(5.1.12)Jacobi迭代法和Gauss-Seidel迭代法的分量形式供计算编程用，它们的距阵形式供研究迭代序列是否收敛等理论分析用。例５.１　用Ｊ法和ＧＳ法分别求解方程组其准确解为　　　　　。解用J法计算，按（5.1.8）有用GS法计算，按（5.1.9）有取，J法迭代4次的计算

6、结果是GS法迭代4次的计算结果是从计算结果看，本例用GS法显然比用J法收敛快。