上海交大数值分析课件数值分析.pptVIP

一、Jacobi迭代法 第六章 解线性方程组的迭代法 §2 基本迭代法 二、Gauss-Seidel迭代法 三、超松弛迭代法 一、Jacobi迭代法 1. Jacobi迭代法举例 例：求解方程组 其中 精确解是x*=(3, 2, 1)T 解：将原方程组改写为 则迭代公式为： 若选 x(0) =(0, 0, 0)T, 则迭代10次有 x(10) = (3.000032, 1.999838, 0.9998813)T 这就是Jacobi迭代法！ 2. Jacobi迭代法一般形式 由方程组 的系数矩阵A非奇异，不妨设 aii≠0，方程组 变形为 对应上述的方程组，可得迭代公式为 其中 x(k)为第 k 次迭代向量. Jacobi迭代法的一般公式 3. Jacobi迭代法的矩阵形式 将方程组记为 Ax = b 其中A非奇异且aii ≠0 (I=1, 2, …, n). 将A分裂为 A = D－Ｌ－Ｕ 其中 由此可将变形过程用矩阵表示为 Dx = (L + U)x +b 即 x = D-1(L + U)x + D-1b 简记为 x = B0x + f 故Jacobi迭代公式的矩阵形式为 二、Gauss-Seidel迭代法 1. Gauss-Seidel迭代法举例 例：求解方程组 精确解是x*=