第五章线性方程组的迭代法(23-24).pptVIP

* 重庆大学数理学院 数 值 分 析 第十二讲 主讲教师： 谭 宏 第五章 线性方程组的迭代法 1．教学内容： 首先通过例子介绍解线性方程组的迭代法的基本思想；然后介绍雅可比迭代公式及其程序设计；介绍高斯-塞德尔迭代公式；超松驰迭代法及其程序设计；以及迭代公式的矩阵表示。 ? 2．重点难点： 雅可比迭代法、高斯—塞德尔迭代法、超松驰迭代法 ? ? 3．教学目标： 掌握三种迭代公式，能利用这三种迭代公式进行线性方程组的迭代求解，并编制相应的应用程序。 前几章研究过的几个课题，无论是插值公式与求积公式的建立，还是常微分方程差分格式的构造，其基本思想都是将其转化为代数问题来处理，最后归结为解线性方程组。工程技术的科学计算中，线性方程组也会经常遇到。因此，线性方程组的解法在数值分析中占有极其重要的地位。 线性方程组的解法大致分为直接法和迭代法两大类。本章将先介绍迭代法，这类算法一个突出优点就是算法简单，因而编制程序比较容易。但是迭代法也有缺点,它要求方程组的系数矩阵具有某种特殊性质，以保证迭代过程的收敛性。发散的迭代过程是没有实用价值的。 §5、1 迭代公式的建立 1、雅可比迭代公式 解线性方程组迭代法的基本思想是将联立方程组的求解归结为重复计算一组彼此独立的线性表达式，这就使问题得到了简化。 求解一般方程 的迭代法的基本想法是： 将 改写为等价的 ，从而得到迭代格式： ，当 时 或很小时，就得到原方程的解。 类似地，可以将这种思想用到方程组的求解上，从而得 到线性方程组的迭代解法。 例： 精确解为 解： 把 ， ， ，从三个方程中分离出来得： 记 得迭代格式： 取初值 ， ， ，…… 从迭代结果看，随着迭代次数的 增加，迭代值 越来越接近精确解 则： 一般地，设 是与 等价的方程组，适 当地取 后，得到迭代格式： 如果当 时， 序列 则称该迭代法收敛，否则称迭代法发散。 由于这里只用到上一步的计算结果 所以，这类迭代又称为单步迭代法或简单迭代法。 这里B称为迭代矩阵 如果，迭代用到前面多步结果 即 则称这类迭代为多步迭代法。 与方程求根的迭代格式设计一样 将方程组AX=b等价转换为X=BX+f的形式也是不唯一的。 例如： 因此，简单迭代法的迭代格式不唯一，而由此产生的 迭代序列 也就不同，有的收敛，有的可能不收敛。对 于收敛的序列，当 k 充分大的时候， 就可以作为方程 组的近似解。 另外，同一简单迭代法也可能关于某个初始向量 收 敛而对于另一个初始向量可能不收敛，即简单迭代法的收敛 性有时也与初始向量的选取有关。 考察一般形式的线性方程组 设有n 阶方程组 若系数矩阵非奇异，且 (i = 1, 2,…, n)，将上述方程组 改写成 然后写成迭代格式 即得到解方程组的雅可比迭代公式: 写成矩阵形式： A = L U D B Jacobi 迭代阵 定理3： 对给定矩阵 B ，若||B||1，则矩阵I-B为非奇异。 证： 用反证法，若||I-B||为奇异阵，则存在非零向量X，使 即有： 所以： 而||B||1，故||X||||X|| 矛盾 定理4： 若迭代矩阵 B满足 则雅可比迭代格式 对于任意初值 均收敛。 证： 因为 据定理3，I-B为非奇异 故方程 有唯一解 由雅可比迭代格式 可知 所以 而 故迭代过程 收敛 2、高斯—塞德尔迭代 例： 解： 把 ， ， ，从三个方程中分离出来得： 由此可建立雅可比迭代公式 （1） （2） （3） 设已得到近似值 则由方程（1）得 对于收敛的迭代过程， 所以，把 代入（2） 则由此计算得到的计算结果 要比用 的好的多。 再把 ， 代入（3） 取初值 按上述格式进行迭代计算得 1.29997 1.19993 1.09989 5 1.30000 1.19999 1.09999 6 1.29972 1.19947 1.09913 4 1.29777 1.19572 1.09313 3 1.28205 1.16719 1.04308 2 1.16440 0.9020