第八章：线性方程组迭代法.PPT

* 《数值分析》第八讲 理学院 张立杰 方程组求解的 迭代方法 §8.1 迭代法 第八章：线性方程组迭代法 记为 1、 一般迭代过程 第八章：线性方程组迭代法 则 构造迭代格式为 例 第八章：线性方程组迭代法 迭代格式 或 第八章：线性方程组迭代法 取 得 第八章：线性方程组迭代法 例2 迭代格式一 取 860 645 6 320 285 5 20 15 2 50 45 3 140 105 4 5 5 1 0 x2 0 x1 0 k 计算程序 x1(1)=0; x2(1)=0; for k=1:6 x1(k+1)=2*x2(k)+5; x2(k+1)=3*x1(k)+5; end 精确解 第八章：线性方程组迭代法 迭代格式二 取 -3.9815 -2.9861 6 -3.9583 -2.963 5 -3.3333 -2.5 2 -3.75 -2.7778 3 -3.8889 -2.9167 4 -2.5 -1.6667 1 0 x2 0 x1 0 k 计算程序 x1(1)=0; x2(1)=0; for k=1:6 x1(k+1)=(x2(k)-5)/3; x2(k+1)=(x1(k)-5)/2; end 第八章：线性方程组迭代法 已知 2、Jacobi迭代过程 令 其中 第八章：线性方程组迭代法 由 得 则有迭代格式为 即 第八章：线性方程组迭代法 由 得 迭代格式为 如 迭代格式 第八章：线性方程组迭代法 如 迭代格式 Gauss-Seidel迭代格式 第八章：线性方程组迭代法 x3 x2 x1 x3 x2 x1 K 1 1 1 1.0001 1.0002 0.9959 1.0041 1.2273 0 2 2 2 2.0001 1.9997 1.9968 2.0289 2.0909 0 3 3 3 2.9998 2.9998 3.0098 2.9773 2.5 0 Gauss-Seidel 1.0031 2.0026 2.9990 7 0.9998 2.0006 3.0002 8 1.0038 2.0000 2.9938 6 0.9205 1.9478 3.0241 4 1.0010 1.9840 3.0003 5 0.9716 2.0455 3.1364 3 1 2.3636 2.8750 2 3 3 2.5 1 0 0 0 0 Jacobi 第八章：线性方程组迭代法 x1(1)=0;x2(1)=0;x3(1)=0; for k=1:8 x1(k+1)=(20+3*x2(k)-2*x3(k))/8; x2(k+1)=(33-4*x1(k+1)+x3(k))/11; x3(k+1)=(36-6*x1(k+1)-3*x2(k+1))/12; end x1(1)=0;x2(1)=0;x3(1)=0; for k=1:8 x1(k+1)=(20+3*x2(k)-2*x3(k))/8; x2(k+1)=(33-4*x1(k)+x3(k))/11; x3(k+1)=(36-6*x1(k)-3*x2(k))/12; end Jacobi Gauss-Seidel 第八章：线性方程组迭代法 §8.2 迭代法的收敛性 如 定义1 收敛于 成立，则称 如果 及矩阵 设矩阵序列 记为 第八章：线性方程组迭代法 迭代法的基本定理 如 860 645 6 320 285 5 20 15 2 50 45 3 140 105 4 5 5 1 0 x2 0 x1 0 k 由 得 设有方程组 对于 及 迭代格式 收敛 其中， 第八章：线性方程组迭代法 又如 得 -3.9815 -2.9861 6 -3.9583 -2.963 5 -3.3333 -2.5 2 -3.75 -2.7778 3 -3.8889 -2.9167 4 -2.5 -1.6667 1 0 x2 0 x1 0 k 由 第八章：线性方程组迭代法 其中 定理：如果方程组系数矩阵严格对角占优或不可约弱对角占优 ，则对任意初值，Jacobi迭代法，Gauss-Seidel迭代法均收敛。 1、严格对角占优 3、不可约 2、弱对角占优 如果存在置换矩阵 使得 称矩阵A可约，否则称矩阵A不可约。 第八章：线性方程组迭代法 如 -3.9815 -2.9861 6 -3.9583 -2.963 5 -3.3333 -2.5 2 -3.75 -2.7778 3 -3.8889 -2.9167 4 -2.5 -1.6667 1 0 x2 0 x1 0 k 860 645 6 320 285 5 20 15 2 50 45 3 140 105 4 5 5 1 0 x2 0 x1 0 k 第八章：线性方程组