第二章2解线性方程组的迭代法讲解.ppt

例： 试分别用雅可比、高斯-塞德尔、超松弛法迭代（取 ）解线性方程组。 解： 取 雅可比迭代公式为 迭代24次后得近似解 高斯-塞德尔迭代公式为 迭代9次后得近似解 超松弛迭代公式为 迭代8次后得近似解 对于迭代法，我们指出两点 （1）从理论上讲，迭代法可以得到任意精度的近似解。但 由于受机器字长的限制，不可能达到任意精度。故用 误差估计式来控制迭代时， 要选取恰当，否则可能 出现死循环。 （2）当给定的方程组不满足迭代的收敛条件时，可适当调 整方程组中方程的次序或作一定的线性组合，就可能 得到满足迭代收敛条件的同解方程组。 例如： 线性方程组 可以验证，无论用雅可比迭代或高斯-塞德尔迭代 均不能满足迭代收敛条件。 但交换其方程的次序为 则该方程组的系数矩阵为严格对角占优。因此， 无论用雅可比迭代或高斯-塞德尔迭代均满足迭代收 敛条件。 例如： 线性方程组 （1） （2） （3） （2）+（1）×2 （2） （3）+（1） （3） 从而与原方程同解的线性方程组 该方程组的系数矩阵为严格对角占优。因此，无论用雅可比迭代或高斯-塞德尔迭代均满足迭代收敛条件。 不尚贤， 使民不争。 不贵难得之货，使民不为盗。不见可欲，使民心不乱。 是以圣人之治，虚其心，实其腹，弱其志，强其骨；常使民无知、无欲，使夫智者不敢为也。 为无为，则无不治。 Thank you！ School of mathematics and statistics 无论是插值公式与求积公式的建立，还是常微分方程差分格式的构造，其基本思想都是将其转化为代数问题来处理，最后归结为解线性方程组。工程技术的科学计算中，线性方程组也会经常遇到。因此，线性方程组的解法在数值分析中占有极其重要的地位。 线性方程组的解法大致分为直接法和迭代法两大类。将介绍迭代法，这类算法一个突出优点就是算法简单，因而编制程序比较容易。但是迭代法也有缺点,它要求方程组的系数矩阵具有某种特殊性质，以保证迭代过程的收敛性。发散的迭代过程是没有实用价值的。 2.4线性方程组的迭代法 第2章 线性代数方程组 在用直接法解线性方程组时要对系数矩阵不断变换，如果方程组的阶数很 2.4 解线性方程组的迭代法 高，则运算量将会很大，并且大量占用计算机资源。因此对线性方程组 要求找寻更经济、适用的数值解法。 对大型的线性代数方程组，常用迭代解法。它是从某些初始向量出发，用 设计好的步骤逐次算出近似解的向量 ，从而得到向量序列 . 一般 的计算公式是 称之为多步迭代法。 若 只与 有关，且 是线性的，即 其中 ，称为单步线性迭代法， 称为迭代矩阵。若 和 都与k无关，即 称为单步定常线性迭代法。下面主要讨论具有这种形式的各种迭代法。 第2章 线性代数方程组 2.4 解线性方程组的迭代法 迭代公式的建立 1、雅可比迭代公式 解线性方程组迭代法的基本思想是将联立方程组的求解归结为重复计算一组彼此独立的线性表达式，这就使问题得到了简化。 求解一般方程 的迭代法的基本想法是： 将 改写为等价的 ，从而得到迭代格式： ，当 时 或很小时，就得到原方程的解。 类似地，可以将这种思想用到方程组的求解上，从而得 到线性方程组的迭代解法。 例： 精确解为 解： 把 ， ， ，从三个方程中分离出来得： 记 得迭代格式： 取初值 ， ， ，…… 从迭代结果看，随着迭代次数的 增加，迭代值 越来越接近精确解 则： 一般地，设 是与 等价的方程组，适 当地取 后，得到迭代格式： 如果当 时， 序列 则称该迭代法收敛，否则称迭代法发散。 由于这里只用到上一步的计算结果 所以，这类迭代又称为单步迭代法或简单迭代法。 这里B称为迭代矩阵 如果，迭代用到前面多步结果 即 则称这类迭代为多步迭代法。 与方程求根的迭代格式设计一样 将方程组AX=b等价转换为X=BX+f的形式也是不唯一的。 例如： 因此，简单迭代法的迭代格式不唯一，而由此产生的 迭代序列 也就不同，有的收敛，有的可能不收敛。对 于收敛的序列，当 k 充分