Chapter4_1_对分法和一般迭代法预案.ppt

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 -5 0 有根区间 中点 f 2 14 0 - 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.7734 0.7695 － 8 0.7695 － 0.7714 － 9 0.7714 － 0.7724 － 10 0.7724 － 0.7729 ＋ 取x10 0.7729，误差为| x* -x10| 1/211 。 Remark:二分法不能用来求重根 f x 0 x g x 等价变换 f x 的根 g x 的不动点 §4.2 单个方程的迭代法 f x 0化为等价方程x g x 的方式是不唯一的,有的收敛,有的发散 For example：2x3-x-1 0 xk+1 g xk 3 1 如果将原方程化为等价方程 由此可见，这种迭代格式是发散的 取初值 2 如果将原方程化为等价方程 仍取初值 依此类推,得 x3 0.9940 x4 0.9990 x5 0.9998 x6 1.0000 x7 1.0000 已经收敛,故原方程的解为 x 1.0000 同样的方程 ? 不同的迭代格式 有不同的结果 什么形式的迭代法能够收敛呢? 收敛性分析 定义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］上存在惟一不动