高斯-塞德尔迭代法ppt课件.ppt

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

1、第6章解线性方程组的迭代法16.1引言考虑线性方程组（1.1）其中为非奇异矩阵，当为低阶稠密矩阵时，第5章所讨论的选主元消去法是有效方法.但对于的阶数很大，零元素较多的大型稀疏矩阵方程组，例如求某些偏微分方程数值解所产生的线性方程组来说，利用迭代法求解则更为合适.迭代法通常都可利用中有大量零元素的特点.2例1（1.2）记为,方程组的精确解是.求解方程组其中现将(1.2)改写为3（1.3）或写为,其中4将这些值代入(1.3)式右边(若(1.3)式为等式即求得方程组的解，但一般不满足).任取初始值，例如取.再将分量代入(1.3)式右边得到，反复利用这个计

2、算程序，得到一向量序列和一般的计算公式(迭代公式)得到新的值5（1.4）简写为其中表示迭代次数迭代到第10次有6从此例看出，由迭代法产生的向量序列逐步逼近方程组的精确解.对于任何由变形得到的等价方程组，迭代法产生的向量序列不一定都能逐步逼近方程组的解.如对方程组7构造迭代法则对任何的初始向量，得到的序列都不收敛.对于给定方程组，设有唯一解，（1.5）又设为任取的初始向量，（1.6）其中表迭代次数.则按下述公式构造向量序列8定义1(1)对于给定的方程组，逐步代入求近似解的方法称为迭代法(或称为一阶定常迭代法，这里与无关).(2)如果存在(记为)，显然就

3、是方程组的解，否则称此迭代法发散.用公式(1.6)称此迭代法收敛，研究的收敛性.引进误差向量由(1.6)减去(1.5)式，得，9要考察的收敛性,就要研究在什么条件下有亦即要研究满足什么条件时有递推得106.2基本迭代法设有（2.1）其中，为非奇异矩阵.将分裂为（2.2）其中，为可选择的非奇异矩阵，且使容易求解，一般选择为的某种近似，称为分裂矩阵.11于是，求解转化为求解，即求解可构造一阶定常迭代法（2.3）其中称为迭代法的迭代矩阵.12选取阵，就得到解的各种迭代法.设,并将写为三部分（2.4）136.2.1雅可比迭代法由，选取为的对角元素部分，解的雅

4、可比(Jacobi)迭代法即选取(对角阵)，，（2.5）其中称为解的雅可比迭代法的迭代阵.由(2.3)式得到14研究雅可比迭代法(2.5)的分量计算公式.记由雅可比迭代公式(2.5),有或于是，解的雅可比迭代法的分量计算公式为15（2.6）由(2.6)式可知，雅可比迭代法计算公式简单，每迭代一次只需计算一次矩阵和向量的乘法且计算过程中原始矩阵始终不变.16(下三角阵)，6.2.2高斯-塞德尔迭代法选取分裂矩阵为的下三角部分，即选取于是由(2.3)式得到解（2.7）其中称为解的高斯-塞德尔迭代法的迭代阵.的高斯-塞德尔(Gauss-Seidel)迭

5、代法17研究高斯-塞德尔迭代法的分量计算公式.由(2.7)式有或即记18于是解的高斯-塞德尔迭代法计算公式为（2.8）或（2.9）19而由高斯-塞德尔迭代公式可知，雅可比迭代法不使用变量的最新信息计算，计算的第个分量时，利用了已经计算出的最新分量.由(2.8)可知，高斯-塞德尔迭代法每迭代一次只需计算一次矩阵与向量的乘法.高斯-塞德尔迭代法可看作雅可比迭代法的一种改进.算法1(高斯-塞德尔迭代法)设，其中为非奇异矩阵且本算法用高斯-塞德尔迭代法解，20迭代一次，这个算法需要的运算次数至多与矩阵的非零元素的个数一样多.数组开始存放,后存放为最

6、大迭代次数.21例2按高斯-塞德尔迭代公式迭代7次，得，（1.2）用高斯-塞德尔迭代法解线性方程组(1.2).取，22由此例可知，用高斯-塞德尔迭代法，雅可比迭代法解线性方程组(1.2)(且取)均收敛.而高斯-塞德尔迭代法比雅可比迭代法收敛较快(即取相同，达到同样精度所需迭代次数较少).但这结论只当满足一定条件时才是对的.且236.2.3解大型稀疏线性方程组的逐次超松弛迭代法选取分裂矩阵为带参数的下三角阵其中为可选择的松弛因子.于是，由(2.3)可构造一个迭代法，其迭代矩阵为从而得到解的逐次超松弛迭代法(SuccessiveOverRelaxati

7、onMethod,简称SOR方法).24解的SOR方法为（2.10）其中研究解的SOR迭代法的分量计算公式.记25或由(2.10)式可得由此，得到解的SOR方法的计算公式（2.11）26或（2.12）关于SOR迭代法，有(1)显然，当时，SOR方法即为高斯-塞德尔迭代法.27(2)SOR方法每迭代一次主要运算量是计算一次矩阵与向量的乘法.(3)当时，称为超松弛法；当时，称为低松弛法.(4)在计算机实现时可用控制迭代终止，或用控制迭代终止.SOR迭代法是高斯-塞德尔迭代法的一种修正.28设已知及已计算的分量(1)首先用高斯-塞德尔迭代法定义辅助量（

8、2.13）(2)再由与加权平均定义，将(2.13)代入(2.14)得到解的SOR迭代(2.11)式.即（2.