Chapter4对分法和般迭代法.ppt

Proof: (1)由g′(x*)=0?必存在?0，当x0∈ ［x* - ? ，x*+ ?］ ? U(x)时，由迭代格式(3)产生的序列{xk}收敛于x*,并有xk∈ ［x* - ? ，x*+ ?］ (2)由泰勒公式有xk+1=g(xk)=g(x* ）＋g′(x*)(xk- x*)+…+g (p-1) (x*)(xk-x*) p-1/(p-1)! + g (p)(x*+ ?(xk-x*))(xk-x*) p /p! ，0?1. 利用g在x*的各阶导数条件及g(x*)=x*,上式可改写成 (11) (3)由于g在x*处p阶连续可微且g(p)(x*)≠0，知必存在x*的某邻域Ｕ(x*)，当x∈U(x*)时，有g (p) (x)≠0. 由于x*+ ?(xk-x*) ∈ ［x* - ? ，x*+ ?］ ?U(x*)，故 g (p)(x*+ ?(xk-x*)) ≠0，k=0，1,2,…. 可见，当初值x0≠x*时，由(11)式可推出xk≠x* 于是由(11)式有 上式令k→∞取极限. 即{xk}有p阶收敛速度. * 方程是在科学研究中不可缺少的工具.方程求解是科学计算中一个重要的研究对象. 几百年前就已经找到了代数方程中二次至五次方程的求解公式.但是,对于更高次数的代数方程目前仍无有效的精确解法.对于无规律的非代数方程的求解也无精确解法.因此,研究非线性方程的数值解法成为必然. 第四章 非线性方程的数值解法  简介（Introduction） 我们知道在实际应用中有许多非线性方程的例子，例如 （1）在光的衍射理论（the theory of diffraction of light)中,我们需要求x-tanx=0的根 （2）在行星轨道（ planetary orbits）的计算中,对任意的a和b，我们需要求x-asinx=b的根 （3） 在数学中，需要求n次多项式xn + a1 xn-1+...+an-1 x + an ＝0的根 求f(x)=0的根 §4.1 对分区间法 (Bisection Method ) 原理：若 f(x) ?C[a, b]，且 f (a) · f (b) 0，则f(x) 在 (a, b) 上必有一根。 a b x1 x2 a1 b2 x* b1 a2 停机条件（termination condition ）： 或 误差 分析： 第1步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： 例1 用二分法求 在(1,2)内的根，要求绝对误差不超过 解： f(1)=-50 有根区间 中点 f(2)=140 -(1,2)+ f(1.25)0 (1.25,1.5) f(1.375)0 (1.25,1.375) f(1.313)0 (1.313,1.375) f(1.344)0 (1.344,1.375) f(1.360)0 (1.360,1.375) f(1.368)0 (1.360,1.368) f(1.5)0 (1,1.5) 12 例2，求方程f(x)= x 3 –e-x =0的一个实根。 因为 f(0)0,f(1)0。 故f(x)在(0,1)内有根 用二分法解之，(a,b）=（0,1)’计算结果如表： k a bk xk f(xk)符号 0 0 1 0.5000 － 1 0.5000 － 0.7500 － 2 0.7500 － 0.8750 ＋ 3 － 0.8750 0.8125 ＋ 4 － 0.8125 0.7812 ＋ 5 － 0.7812 0.7656 － 6 0.7656 － 0.7734 ＋ 7 － 0.7