[工学]第三章-迭代法.ppt

* § 3.3 误差估计 * * * 注1：由 * 注2：由 * * * * §4　最速下降法与共轭梯度法 * 1.确定搜索方向 2.求极小值 § 4.1 最速下降法 * 1.确定搜索方向 § 4.2 共轭梯度法 2.求极小值 * Ex. P87 3,5(1), 6, 8(1), 10, 12. * Jacobi : 各个分量的计算顺序无关, 可用于并行计算 G-S: 减少了存储量, 要求计算顺序, SOR: 是G-S的推广, 或是G-S的加速 迭代法: 求解大型稀疏线性方程组的常用方法, 保持矩阵的稀疏性 * 分裂A 换一个角度 * * * * 例：用Jacobi和Gauss－Seidel迭代求解方程组 　　取 * 解：Jacobi迭代的计算公式为 * G－S迭代的计算公式为 * x=Gauss_Seidel_iter(A,b,20,1e-5,x0) k x1 x2 x3..... 1 1 -1 1 2 5 -5 1 3 13 -13 1 4 29 -29 1 5 61 -61 1 6 125 -125 1 7 253 -253 1 8 509 -509 9 1021 -1021 1 10 2045 -2045 1 11 4093 -409 1 12 8189 -81 1 ……….. 20 2097149 -2097149 1 Gauss-seidel method did not converge 问题：G－S迭代法， Jacobi迭代法的收敛需要 条件？由什么来决定收敛性呢？ * §3　迭代法的收敛条件 § 3.1 矩阵的谱半径 Re Im ? ? ? ? ? ? ? ? ? (A) * x 具体证明见P73 定理3.5 * § 3.2　迭代法的收敛条件 * * 充分条件 * * 例：设有迭代格式： 试证此迭代格式对任意初始向量都收敛. * 解：求B的谱半径 * 严格对角占优 非对角占优 主对角线元素的绝 对值大于同行其它 元素绝对值之和 * 充分条件 * 充分条件 * 充分条件 看系数矩阵 * 充分必要条件 看迭代矩阵 * * 例： 方程组Ax=b的系数矩阵为 严格对角占优阵. * 例：方程组Ax=b的系数矩阵为 但A为对称正定矩阵， 非严格对角占优， * 解： * * 例：设有方程组 试写出收敛的迭代式并说明理由. * 解：原方程组等价于 * (1)Jacobi 迭代公式： * (2)G-S 迭代公式： * * 第三章 解线性方程组的迭代法 §1 迭代法概述 求解大型稀疏线性方程组的常用方法,保持矩阵的稀疏性 基本思想： * § 1.1 向量序列矩阵序列的收敛性: * * * §1.2 迭代法的一般形式 * * §2　几种基本的迭代法 * §2.1　雅可比（Jacobi）迭代法 * * * 与系数矩阵的关系 * * 例1 用Jacobi迭代法求解线性方程组 解：Jacobi迭代计算公式为 * * x=Jacobi_iter(A,b,20,1e-5,x0) k x1 x2 x3.... 1.0000 7.2000 8.3000 8.4000 2.0000 9.7100 10.7000 11.5000 3.0000 10.5700 11.5710 12.4820 4.0000 10.8535 11.8534 12.8282 5.0000 10.9510 11.9510 12.9414 6.0000 10.9834 11.9834 12.9804 7.0000