第二章 非线性方程求解.ppt

第二章 非线性方程求解 第二章 非线性方程求解目录 第二章 非线性方程求解概述 求根的隔离区间的两种方法 §1 对分法 对分法的误差估计 对分法的优缺点 §2 简单迭代法 2.1 迭代法的基本思想 迭代法举例续 迭代法的几何含义 迭代法的几何含义（续） 压缩映象原理的证明 2.2 迭代法的收敛条件（三大定理） 压缩映象原理的证明 压缩映象原理的证明（续1） 压缩映象原理的证明 两个重要误差公式说明 迭代法的收敛条件之二 证明 推论1应用举例 应用举例 迭代法的收敛条件之三 2.3 Steffensen方程—简单迭代法的加速 {xn}的r 阶收敛定理 定理2.4（续） 1、 Aitken加速法 Aitken加速法（续） §3 Newton法与弦截法 Newton法（续） Newton法的几何意义 Newton法举例 Newton法（续） Newton法收敛定理 Newton法的几何意义及其优劣 3.2计算重根的牛顿迭代法 计算重根的牛顿迭代法(续1） 计算重根的牛顿迭代法(续2） 3.2 弦截法 弦截法迭代公式的几何解释 弦截法的几点说明 弦截法的几点说明(续） 则 由式（8-8）可知 故 因此可用 下式估计m 例8 用牛顿迭代法求方程 在0.95附近之根。 解 取x0=0.95，用牛顿迭代法， 按式（3-7）求得的xk见表3-3，由表中数据可见xk 收敛很慢。由 可知，所求根为m=2重根， 改用式(8-9)迭代格式，得： 收敛速度大大快于直接用牛顿迭代公式（8-6）. 表3-3 k xk ?k 0 1 2 3 4 5 6 0.95 0.9744279 0.9870583 0.9934878 0.9967328 0.9983576 0.9991901 ? ? ? 0.5090 0.5047 0.5007 0.5125 ? ? ? 2.0369 2.0190 2.0028 2.0511 m 不 足 之 处：需要计算导数值，较难； ——这就是弦截法迭代公式 Newton法优点：收敛快（平方阶），固定格式； 修 正：以差商代替导数（微商） 采用的三种迭代格式， 在隔根区间（1,1.2）内有： 用推论1判别简单迭代法的收敛性比定理8.1方便 如对例题5： 第一种迭代格式发散，第二、三种迭代格式收敛且第三种迭代格式比第二种迭代格式中的L要小，因而收敛要快得多，这与实际迭代结果完全吻合。 故可取n = 7，只需迭代7次就可达到所要求的精度。 根据推论1可知， 对第三种迭代格式，为使与方程近似根的误差不超过 10-6，可估计迭代次数： Leonardo于1225年研究了方程 曾经轰动一时，因为没有人知道他用的是什么方法。 我们现在可用迭代法求解： 还可用Newton法，弦截法求解 定理2.7 定理2.3强调迭代初值x0应取在根x*的邻域中。如果对任意给定的x0，迭代格式均收敛，则称此格式具有全局收敛性，但这样的格式是极其稀少的。如果对根x*的某邻域内的任一点x0，迭代格式均收敛，则此格式具有局部收敛性。 即可保证对其中任取的一点x0迭代收敛。事实上，在用迭代法求解方程（2-1）时，常常先用对分区间求得较好的初值，然后再进行迭代。 本定理给出的就是局部收敛性条件。具体解题时，虽然无法判别隔根区间是否为以x*为中心的邻域，但只要它足够小，且在邻域中满足： 收敛速度（收敛速度的阶）： 成立，则称{xn}是r 阶收敛的，或称{xn}的收敛阶为r，收敛阶r 的大小刻划了序列{xn}的收敛速度： r 越大，收敛越快: r =1 线性收敛 r 1 超线性收敛 r =2 平方收敛 设序列{xn}收敛于x*，若存在正数r和a使得： 定理2.8 设迭代函数? (x)在x*邻近有r阶连续导数，且 x* = ?(x*)，并且有 证明： 1） ?(x)满足收敛定理的条件?{xn}?x*； 紧接下屏 2）利用Taylor公式将?(x)在x*附近展开： 这表明：{xn}是r阶收敛的。 一阶收敛即为线性收敛，收敛速度较慢，下面想法加速： 证毕！ 若序列{xn}线性收敛于x*，可按式： 当n充分大时，有： 紧接下屏 由此式可推导出： 由此可得比值： 3.1 Newton法 将非线性方程线性化，以线性方程的解逐步逼近非线性方程的解，这就是Newton法的基本思想。 设已知方程f (x) = 0的近似根x0，f (x)在其零点x*邻近一阶连续可微， 且f ?(x) ? 0，当x