用迭代法解线性方程组.PDF

用迭代法解线性方程组 Jacobi迭代和Gauss-Seidel迭代 由于收敛速度较慢，已经越来越不 适应当前信息时代人们对计算速度 和精度的要求，所以在实际应用中 使用得并不多。但是，由于它们体 现了迭代法的最基本的思想，是学 习其它迭代法的基础。 引言 直接法 是通过有限步运算后得到线性方程组的解, 解线性方程组还有另一种解法，称为迭代法,它的基本思 想是将线性方程组Ax=b 化为 x=Bx+f (k) 再由此构造向量序列{x }: x(k+1)=Bx (k)+f (k) *, * 若{x } 收敛至某个向量x 则可得向量x 就是所求方 程组Ax=b 的准确解. 线性方程组的迭代法主要有Jocobi迭代法、Seidel迭 代法和超松弛(Sor)迭代法. 迭代法的特点 * 若在求解过程中x x (k) ，由x =(x ) 产生的迭 k k+1 k * 代x 向x 的逼近 ，在数次迭代求解之后，由于机器跳动产 k 生的x 值误差或是有效数字产生的舍入误差，都会在第k+1 k 次迭代计算中自动弥补过来或逐步纠正过来。因此，在迭 代求解过程中产生的各种误差是可以忽略的，即迭代求解 无累积误差，实际上，x 只是解的一个近似，机器的舍入 k 误差并不改变它的此性质。 如果矩阵A=(a ) 满足 ij n |a |  |a | i=1,2,……n, ii ij j=1,j i 则称方阵A是严格(行)对角占优的. U a11 a12 a13 … a1n a21 a22 a23 … a2n A= … … … … … =L+D+U an1 an3 an4 … ann D -4 2 1 L 例矩阵A= 1 -9 7 2 -6 10 Jacobi 迭代  一：设有方程组 a x +a x +····+a x =b 11 1 12 2 1n n 1 a x +a x +····+a x =b 21 1 22 2 2n n 2 . . . . . . . . . . . . .