计算方法(方程组的迭代法)3.pptVIP

计算方法 北京科技大学应用学院数力系 卫鸿儒 Weihr168@ 课程性质和计划（续） 二、用迭代法解线性方程组 Jacobi迭代和Seidel迭代由于收敛速度较慢，已经越来越不适应当前信息时代人们对计算速度和精度的要求，所以在实际应用中使用的并不多。但是，他们体现了迭代法的最基本的思想，是学习其它迭代法的基础。 （一）引言 直接法是通过有限步运算后得到线性方程组的解，解线性方程组还有另一种解法，称为迭代法，它的基本思想是将线性方程组 Ax=b 化为 x=Bx+f 再由此构造向量序列{x (k)}: x(k+1)=Bx (k)+f 若{x (k)}收敛至某个向量x *，则可得向量x *就是所求方程组 AX=b 的准确解。 线性方程组的迭代法主要有Jocobi迭代法、Seidel迭代法和超松弛(Sor)迭代法。 （二）、Jacobi 迭代 一： 设有方程组 a11x1+a12x2+····+a1nxn=b1 a21x1+a22x2+····+a2nxn=b2 . . . . . . . . . . . . an1x1+an2x2+····+annxn=bn 用矩阵表示：Ax =b （A 为系数矩阵，非奇异；b为右端列向量，x为解向量） 把方程组写成容易迭代的形式： 假设 aii?0 令 cij = -aij /aii (i?j) gi= bi /aij , i=1,2,3?,n 则 x1(k+1)= c12x2(k)+c13x3(k)+???? +c1nxn(k)+g1 x2(k+1)=c21x1(k) +c23x3(k)+???? +c2nxn(k)+g2 。。。。。。。。。。。。。。。。。。。。。。。。。。。。 xn(k+1)=cn1x1(k) +cn2x2(k)+???? +cn(n-1)xn-1(k) + gn Jacobi迭代格式 若令 0 c12 c13 … c1n x1 c21 0 c23 … c2n x2 BJ= … … … … … x= .. cn1 cn3 cn4 … 0 xn a11 g1 a22 g= g2 易看出:BJ =D-1(D-A)=I-D-1A D= .. .. ann gn （四）、收敛性及误差估计 （九）、迭代法的特点 （1）方法简单，每次迭代都是简单的重复运算，易于编制程序；与求解线性方程的精确法相比，简单迭代法对于字长位数较少的计算机更为适用，它可以用增加迭代次数来弥补字长位数少的不足。 （2）初值可以任取，因而中间结果偶然错误不影响最后结果的获得。 （3）缺点：用计算机计算时，迭代速度较慢。 （4）就其收敛性而言，某些用Seidel迭代法不能收敛。而无法得出结果的线性代数方程组，用Jacoai迭代法却能进行收敛计算，反之已然。 （五）、例题及求解 例：用迭代法解方程组AX=b,其中 [已知该方程的解为 ] 　 解：本题分别用简单迭代法（Jacobi迭代法）和GS迭代法来解　 　(1)简单迭代法 　　　　　　　　　　　　　　　　 表１ 表2 （六）、相关程序设计 原始数据（A,b)可用一个二维数组存储，也可将A用一个二维数组，b用一个一维数组分别存储，存储　　需要一个一维数组。程序中应方便地对迭代方法和终止条件的选择以及对初始向量和?值的设置。在迭代过程中，为反映迭代情况，可设置一些中间数据的输出，如迭带次数，迭代向量，迭代残向量等。当然不需要每迭代一次都作输出，这可作为收敛情况或不收敛情况的分析。作为不收敛的判定，可设置一个大的整数，当迭代次数超过该数时作为不收敛