93~数值分析(阅读).pptVIP

§3.8 椭圆型差分方程的迭代法 一、迭代法的基本理论 1、迭代公式 线性方程组 （3.45） 其中A是N阶非奇异矩阵，x，b均为N维列向量。 将方程组改写为 的迭代公式 （3.47） 2、迭代的收敛性判断 定理3.4 解方程组（3.45）的迭代格式（3.47） 对任意右端 c 及任意初始向量 X（0） 收敛的充分必 要条件为 推论1 迭代格式（3.47）收敛的充分条件为 3、收敛速度 迭代公式 准确解满足 记解的误差向量 误差向量满足 设迭代矩阵 G 有 N 个线性无关的特征向量 vi，分别 对应于特征值λi，且 由 得 对于相当大的 n ，有 则 要求从第 n 步迭代到第（n+p）步的误差减少为 e（n） 的10–q ，即要求 两边去对数得 或 越小， 越大，称 为迭代法的渐近收敛速度。 二、Jacobi 迭代和 Gauss-Seidel 迭代 方程组 A = D – L – R 方程组写成 若 ，Jacobi 迭代格式为 或 即 Gauss-Seidel 迭代法 或 即 迭代矩阵 三、椭圆型差分格式的 Jacobi 迭代和Guass-Seidel 迭代 一般二阶线性椭圆型方程的五点差分格式 1、 Jacobi 迭代格式（3.63） 2、 Guass-Seidel 迭代格式(3.64) 3、收敛速度 考虑Laplace方程的Dirichlet条件问题 Jacobi 迭代矩阵的谱半径为 收敛速度为 由于 所以 Guass-Seidel 迭代矩阵的谱半径为 收敛速度为 例1 利用五点差分格式近似Laplcce方程的Dirichlet问题 取步长h=1，试用Jacobi 迭代法求解，迭代误差限ε=0.1。 解 Jacobi 迭代法格式 dirichjac.m lapjacobi 四、超松弛迭代法 1、逐次超松弛迭代法 方程组 当 时，逐次超松弛迭代格式为 （1）计算 (3.75.1) （2）加权平均 (3.75.2) 或 即 矩阵形式 记逐次超松弛迭代矩阵为 则SOR迭代格式为 2、收敛性的判断 定理3.5 SOR方法收敛的充分必要条件为 方程 的根 η按模小于1。 定理3.6 对所有ω成立 证 即 故 由定理3.5 SOR收敛的必要条件为 ， 若A对称 正定 ，则SOR收敛的充分条件为 。 为低超松弛因子， 为超松弛因子， ，SOR是 G-S 迭代。 3、最佳松弛因子的确定 定义 使