数值计算课件数值分析第4章 一元线性方程的解法教案.pptVIP

图 4.9 与例2比较,同一例用一般迭代法要迭代十次才能得到满足精度ε=10-3的结果,而这里仅迭代三次便可达到ε=10-5的高精度结果。这种加速过程取得的效果极为显著。 为了避免计算导数a≈g′(x),下面介绍埃特金(Aitken)迭代方法。它也是一种加速迭代法。 (4―31) 将式(4―27)与式(4―31)联立消去a得到 可解出 (4―32) (4―33) 这样得到埃特金迭代公式 (4―33) 例6用埃特金迭代法求 ?x3-x-1=0 ?在(1,1.5)内的根。 解 前面已经提到,迭代公式 ?x k+1=x3k-1, k=0,1,2,… ?是发散的。 现用埃特金算法来求根,其迭代公式为 仍取x0=1.5,计算结果见表4―7。 表 4―7 图 4.4 迭代过程(4―8)就是在x轴取初始近似值x0,过x0作y轴的平行线交曲线y=g(x)于p0,p0的横坐标为x0,纵坐标为g(x0)(g(x0)=x1),也即 p0(x0,x1) 再在x轴上取x1作为新的近似值,过x1作y轴的平行线交曲线y=g(x)于p1,p1的横坐标为x1,纵坐标为 g(x1)(g(x1)=x2),也即 p1(x1,x2) 而这相当于过p0引平行于x轴的直线交y=x于 ? Q1(x1,x2) 再过Q1引平行于y轴的直线交曲线y=g(x)于 ? p1(x1,x2) ?仿此可得到点列 ?p0(x0,x1),p1(x1,x2),p2(x2,x3),… ?若 则迭代法收敛,见图4.4(a);否则迭代法发散,见图4.4(b)。 必须说明两点: ①要验证g(x)是否满足李氏条件一般比较困难,若g(x)可微,可用充分条件 来代替。这里q＜1是非常重要的条件,否则不能保证迭代收敛。 ②对于收敛的迭代过程,误差估计式(4―11)说明迭代值的偏差｜xk-xk-1｜相当小,就能保证迭代误差｜x-xk｜足够小。因此在具体计算时常常用条件 ｜xk-x k-1｜＜ε (4―15) 来控制迭代过程结束。 迭代法的突出优点是算法的逻辑结构简单,且在计算时,中间结果若有扰动,仍不会影响计算结果。其计算步骤为: (1)确定方程f(x)=0的等价形式x=g(x),为确保迭代过程的收敛,要求g(x)满足李普希茨条件(或｜g′(x)｜≤q＜1); (2)选取初始值x0,按公式 ? x k+1=g(xk), k=0,1,2,… ?进行迭代; (3)若｜x k+1-xk｜＜ε,则停止计算,x≈x k+1。 例2 求方程 x=e-x 在x=0.5附近的一个根。按五位小数计算,计算结果 的精度要求为ε=10-3。 解 过x=0.5以步长h=0.1计算 f(x)=x-e-x 由于 f(0.5)＜0, f(0.6)＞0 故所求的根在区间(0.5,0.6)内,且在x=0.5附近 图 4.5 表 4―3 因此用迭代公式 ? 由表可见 为方程 最后,我们给出一个说明,在将方程(4―1)化为等价形式(4―7)时,g(x)的形式是多种多样的,选取不当,迭代公式(4―8)就不会收敛。最一般的形式可以写成 ?