Chapter3_1_对分法和一般迭代法..pptVIP

数值分析第三章 非线性方程的数值解法  数学与统计学院 简介（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的根 Remark1：求奇数个根 Find solutions to the equation Remark2:要区别根与奇异点 Remark3:二分发不能用来求重根 收敛性分析 定义2 若存在常数?(0≤ ?1),使得对一切x1,x2∈［a,b］, 成立不等式 |g(x1)-g(x2)|≤ ? |x1-x2|， (5) 则称g(x)是［a,b］上的一个压缩映射， ?称为压缩系数 考虑方程 x = g(x), g(x)?C[a, b], 若 ( I ) 当 x?[a, b] 时， g(x)?[a, b]； ( II )在［a,b］上成立不等式：|g(x1)-g(x2)|≤ ? |x1-x2| 。 则（1）g在［a,b］上存在惟一不动点x* （2）任取 x0?[a, b]，由 xk+1 = g(xk) 得到的序列 {xk}（?［a,b】） 收敛于x* 。 （3）k次迭代所得到的近似不动点xk与精确不动点x*有有误差估计式： Remark： 定理条件非必要条件，而且定理3.2.1中的压缩条件不好验证，一般来讲， 若知道迭代函数g(x)∈C1『a,b］,并且满足|g′(x)|≤ ? ≤1,对任意的x∈［a,b］, 则g（x)是［a,b］上的压缩映射 例题 已知方程2x-7-lgx＝0，求方程的含根区间，考查用迭代法解此方程的收敛性。 在这里我们考查在区间[3.5,4]的迭代法的收敛性 很容易验证：f(3.5)0,f(4)0 将方程变形成等价形式：x＝(lgx+7)/2 局部收敛性定理 定理3.2.2设x*为g的不动点，g(x)与g′(x)在包含x*的某邻域U(x*) (即开区间)内连续，且|g′(x*)|1,则存在?0，当x0∈［x* - ? ，x*+ ?］时，迭代法(3)产生的序列{xk}? ［x* - ? ，x*+ ?］且收敛于x*. 证明略（作为练习） 举例 用一般迭代法求x3-x-1=0的正实根x* 例题 用一般迭代法求方程x-lnx＝2在区间（2，?）内的根，要求|xk-xk-1|/|xk|=10-8 将方程化为等价方程：x＝2＋lnx k xi 0 3.000000000 1 3.098612289 2 3.130954362 3 3.141337866 4 3.144648781 5 3.145702209 6 3.146037143 另一种迭代格式： 程序演示 由此可见，对同一个非线性方程的迭代格式，在收敛的情形下，有的收敛快，有的收敛慢。 定理3.2.3 设x*为g的不动点，p≥2为正整数，g在x*的某邻域Ｕ(x*)内p阶连续可微，且 g′(x*)=g″(x*)=…=g (p-1)(x*)=0，而g(p)(x*)≠0， 则存在?0，当x0∈ ［x* - ? ，x*+ ?］ (x0≠x*)时，由迭代法(3)产生的序列｛xk｝以p阶收敛速度收敛于x*. Prove: (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*,上式可改写成 (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)式有 作业（homework) P159：第二题 P160：第三题