清华第五版数值分析第6章ppt课件.pptVIP

本题考查学生综合知识解决问题的能力。最早研制成功的火车、飞机分别是在第一、二次工业革命中。火车动力来源于蒸汽机，飞机动力来源于内燃机 第五章 线性代数方程组的数值解法 迭代法 求解 一. 迭代法的基本思想 B：迭代矩阵 向量和矩阵收敛性定义 定义：向量收敛 范数为任一种向量范数。 定义：矩阵收敛 范数为矩阵的任一种算子范数。 定理： 定理：设 则下面三命题等价。 （1） （2） （3）至少存在一种从属的矩阵范数 有 五 迭代法的收敛性 定理：设有线性方程组 及迭代格式 ， 若有 的某种算子范数 ，则 （1）迭代法收敛，即对任意初值 ，有 且 （2） （3） （4） 迭代法研究的主要问题 1）迭代格式的构造； 2）迭代的收敛性分析； 3）收敛速度分析； 4）复杂性分析。（计算工作量） 二 . 雅可比（Jacobi）迭代法 建立与式(1)相等价的形式： 例如：考虑方程组 据此建立迭代公式： 一般地，设方程组 AX=b , 通过分离变量的过程建立Jacobi迭代公式，即 由此我们可以得到 Jacobi 迭代公式： 或 将系数矩阵分裂为： 其中 ? 雅可比迭代法的矩阵表示 BJ fJ Jacobi 迭代阵 Jacobi迭代矩阵 例1：设方程组为 解：Jacobi迭代格式为 试写出其Jacobi迭代格式。 Jacobi迭代算法描述 1 输入 2 if , then 2.1 for 2.1.1 s=0, 2.1.2 for 2.1.3 2.1.4 if then 2.2 k=k+1 2.3 if then 2.3.1 2.3.2 goto 2 else 输出 结束。 else 2.4 输出 迭代次数太大。 3 结束 注意到利用Jacobi迭代公式计算 时，已经计算好了 的值，而Jacobi迭代公式并不利用这些最新的近似值计算， 仍用 这启发我们可以对其加以改进,即在每个分量的计算中尽量利用最新的迭代值，得到 上式称为 Gauss-Seidel 迭代法. 三. 高斯-塞德尔（Gauss-Seidel）迭代法 … … … … 即： BG-S Gauss-Seidel 迭代阵 写成矩阵形式： fG-S 例1的Jacobi 迭代公式为 例1的Gauss-Seidel 迭代公式为 G-S迭代算法描叙 1 输入 2 if , then 2.1 for 2.1.1 s=0， 2.1.2 for 2.1.3 2.1.4 if then 2.2 k=k+1 2.3 if ,输出结果，结束。 else 2.4 输出迭代次数太大。 3 结束 Remark: Gauss-Seidel迭代法的计算过程比Jacobi迭代法更简单。计算过程中只需用一个一维数组存放迭代向量。 Gauss-Seidel迭代不一定比Jacobi迭代收敛快。 思想 类似于G-S迭代法的改进方法，利用第k次 迭代值和第k +1次的G-S迭代值作加权平均 G-S迭代法的计算公式 ： 作加权平均 四、逐次超松弛迭代法（SOR方法） ?超松弛(SOR)迭代法的分量形式： 利用高斯-塞德尔迭代法求得: 例1的SOR迭代公式为 例1