工程数学ch03.pdf

工程数学 Engineering Mathematics 第三章 线性方程组的数值解法 • 设有线性方程组 Ax = b ，其中 ⎛a11 a12 L a1n ⎞ ⎡x1 ⎤ ⎡b1 ⎤ ⎜ ⎟ ⎢ ⎥ ⎢ ⎥ A ⎜⎜a21 a22 L a2n ⎟⎟为非奇异阵， x ⎢x2 ⎥ b ⎢b2 ⎥ L L O L ⎢ M ⎥ ⎢ M ⎥ ⎜⎜⎝an 1 an 2 L ann ⎠⎟⎟ ⎢⎣xn ⎦⎥ ⎢⎣bn ⎦⎥ 线性方程组的数值解法一般有两类：直接法与 迭代法 – 直接法 – 就是经过有限步算术运算，可求得方程组精确解的方 法(若计算过程中没有舍入误差)。 – 但实际计算中由于舍入误差的存在和影响，这种方法 也只能求得线性方程组的近似解。 – 这类算法中最基本的是高斯消元法及其某些变形。 – 这类方法是解低阶稠密矩阵方程组的有效方法，近十 几年来直接法在求解具有较大型稀疏矩阵方程组方面 取得了较大进展。 –迭代法 –迭代法是用某种极限过程去逐步逼近线性方程组精确 解的方法。 –迭代法具有需要计算机的存贮单元较少、程序设计简 单、原始系数矩阵在计算过程中始终不变等优点，但 存在收敛性及收敛速度问题。 –迭代法是解大型稀疏矩阵方程组(尤其是由微分方程离 散后得到的大型方程组)的重要方法。 线性方程组的迭代求解方法 3.1 几个常用的迭代方法 3.2 迭代法的收敛性与误差估计 3.3 方程组的性态与条件数 3.1 几个常用的迭代方法 一、简单迭代法－－Jacobi迭代法 二、Gauss-Seidel迭代法 三、松驰迭代法－－SOR迭代法 1. Jacobi迭代法 1. Jacobi迭代法 设n阶线性方程组 Ax=b 或 + + + a x a x L a x b ⎧ 11 1 12 2 1n n 1 ⎪ + + + a x a x L a x b ⎪ 21 1 22 2 2n n 2 ⎨ (3.1) ⎪L L ⎪ + + + a x a x L a x b ⎩ n1 1 n 2 2 nn n n 其中 A (a ) , x (x , x , L, x )T , b (b , b , L, b )T ij n×n 1 2 n 1 2 n 假设 aii ≠0 ，令 bij −aij / aii , i ≠j f b / a , i 1,2,