数值分析解线性方程组的迭代法.ppt

§ 1 引言 第 6 章 解线性代数方程组的迭代法 考虑线性方程组 也就是 AX=b. (1.1) ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? n n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a ? ? ? ? ? ? 2 2 1 1 2 2 2 22 1 21 1 1 2 12 1 11 低阶稠密的线性方程组用直接法 ( 如高斯消去法和三 角分解法 ) 。 大型稀疏非带状的线性方程组 (n 很大 , 且零元素很多 . 如偏微方程数值解产生的线性方程组,n≥10 4 ) 的求解 问题？ 零元素多，适合用迭代法。 我们将介绍迭代法的一般理论及雅可比迭代法、高 斯 — 塞德尔迭代法、超松弛迭代法，研究它们的收 敛性。 例 1 求解线性方程组 (1.2) . 36 12 3 6 , 33 11 4 , 20 2 3 8 3 2 1 3 2 1 3 2 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? x x x x x x x x x 记为 Ax=b, 即 . 36 33 20 12 3 6 1 11 4 2 3 8 3 2 1 x x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 精确解 x * =(3,2,1) T . 改写 (1.2) 为 (1.3) ). 36 3 6 ( ), 33 4 ( ), 20 2 3 ( 2 1 12 1 3 3 1 11 1 2 3 2 8 1 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? x x x x x x x x x 或写为 x=B 0 x+f, 即 . 0 0 0 12 36 11 33 8 20 3 2 1 12 3 12 6 11 1 11 4 8 2 8 3 3 2 1 x x x x x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 任取初值 , 如 x (0) =(0,0,0) T , 代入 (1.3) 得到 x (1) = (2.5,3,3) T . 反复迭代 (1.4) . 12 / ) 36 3 6 ( , 11 / ) 33 4 ( , 8 / ) 20 2 3 ( ) ( 2 ) ( 1 ) 1 ( 3 ) ( 3 ) ( 1 ) 1 ( 2 ) ( 3 ) ( 2 ) 1 ( 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? k k k k k k k k k x x x x x x x x x 即 x (k+1) =B 0 x (k) +f ， (k=0,1,2,…) 解。 逐步逼近此方程的精确 序列 代法产生的向量 从此例可以看出，由迭 其中 ) ( ) 10 ( ) 10 ( ) 10 ( ) 10 ( x . * , 000187 . 0 , ) 9998813 . 0 , 999838 . 1 , 000032 . 3 ( k T x x ε ε x ? ? ? ? ? . f Bx x b Ax ? ? ? 变形得到等价的 一般地，由 (1.5) * * *, f Bx x x ? ? 则 设有唯一解 (1.6) , ) ( ) 1 ( ) 0 ( f Bx x x ? ? ? k k 则可构造迭代序列 又设任取初值 . ) 6 . 1 ( (1) 1 ) 迭代法( 定义 一阶定常迭代法 近似解的方法称为 逐步代入求 ，用公式 对于方程组 f Bx x ? ? . , * , *) ( lim ) 2 ( ) ( 否则迭代法发散 是解 显然 ，则称此 记为 存在 若 x x x 迭代法收敛 k k ? ? .