计算方法(方程组的迭代法)3教程教案.ppt

北京科技大学应用学院数力系 卫鸿儒 Weihr168@;课程性质和计划（续）;二、用迭代法解线性方程组; 若在求解过程中 xk?x*(k??)，由 xk+1=?(xk)产生的迭代 xk向x*的逼近 ，在数次迭代求解之后，由于机器跳动产生的xk值误差或是有效数字产生的舍入误差，都会在第k+1次迭代计算中自动弥补过来或逐步纠正过来。因此，在迭代求解过程中产生的各种误差是可以忽略的，即迭代求解无累积误差，实际上，xk只是解的一个近似，机器的舍入误差并不改变它的此性质。;（二）、Jacobi 迭代;假设 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 ;;（三）、Seidel迭代法 为了加快收敛速度，同时为了节省计算机的内存，我们作如下的改进：每算出一个分量的近似值，立即用到下一个分量的计算中去，即用迭代格式： 这样所得的迭代法就称为Gauss-Seidel迭代法，也称为“异步迭代法”，简称为GS迭代法．利用Ax=b 及A=L+D+U,其中D为对角矩阵,L,U分别为严格下，上三角矩阵．则有，GS迭代法的矩阵形式为: 　 　　　　　　 　;Seidel迭代法的具体形式;（四）、收敛性及误差估计;;;;;;;;;;;;;;;;;;;;;;;;;;;;;（五）、例题及求解 例：用迭代法解方程组AX=b,其中 [已知该方程的解为 ] 　 解：本题分别用简单迭代法（Jacobi迭代法）和GS迭代法来解　 　(1)简单迭代法 　　　　　　　　　　　　　　　　;;表１;;;???2;;;;;;（六）、相关程序设计 原始数据（A,b)可用一个二维数组存储，也可将A用一个二维数组，b用一个一维数组分别存储，存储　　需要一个一维数组。程序中应方便地对迭代方法和终止条件的选择以及对初始向量和?值的设置。在迭代过程中，为反映迭代情况，可设置一些中间数据的输出，如迭带次数，迭代向量，迭代残向量等。当然不需要每迭代一次都作输出，这可作为收敛情况或不收敛情况的分析。作为不收敛的判定，可设置一个大的整数，当迭代次数超过该数时作为不收敛处理。 GS 迭代法的计算公式为： ;;开始;请给出用C语言或其他语言求解下面方程组的程序及结果：;（七）、方法优缺点讨论 　　由以上例题的求解过程可明显看出GS迭代法的收敛速度比简单迭代法快，但对于任意给定的一个方程组分别用简单迭代法和GS迭代法求解时，两种迭代法可能都收敛，也可能都不收敛。也有可能是GS迭代法收敛而J迭代法不收敛。但亦有相反情况，即简单迭代法收敛而GS迭代法不收敛。而且交换方程组中的方程和未知数的次序都会影响GS迭代法的计算结果，但这种交换对简单迭代法是没有影响的。　 ;（八）、SOR法介绍;;;;（九）、迭代法的特点