南开大学计算物理讲义4.ppt

* 2. The Gauss-Seidel Iteration In the Jacobi Iteration * In the Gauss-Seidel Iteration * 高斯――赛得尔（Gauss-Seidel）迭代法 （i = 1, 2,…, n） * Jacobi 迭代法 B = I – D-1A F = D-1 b Gauss-Seidel迭代法 * 用高斯――赛得尔（Gauss-Seidel）迭代法求解方程组 前三次叠代的结果，要求取 * * * Gauss-Seidel Iteration Jacobi Iteration * 从此例看出,高斯—塞德尔迭代法比雅可比迭代法收敛快(达到同样的精度所需迭代次数少),但这个结论,在一定条件下才是对的,甚至有这样的方程组,雅可比方法收敛，而高斯—塞德尔迭代法却是发散的. * 定义1：如果矩阵的每一行中，不在主对角线上的所有元素绝对值之和小于主对角线上元素的绝对值，即 则称矩阵A按行严格对角占优，类似地，也有按列严格对角占优。 迭代收敛的充分条件 * 定理：若线性方程组AX = b的系数矩阵A按行严格对角占优，则雅克比迭代法和高斯―赛得尔迭代法对任意给定初值均收敛。 如果矩阵A严格对角占优，那么高斯――赛得尔迭代法的收敛速度快于雅克比迭代法的收敛速度。 * 用 Gauss-