数值分析08非线性方程求解.ppt

阜师院数科院 第八章 非线性方程求解 第八章 非线性方程求解目录 第八章 非线性方程求解概述 非线性方程求解概述(续） 求方程根近似解的几个问题： 求根的隔离区间的两种方法 例2（续） §1 对分法 对分法（续） 对分法的误差估计 对分法举例 对分法的优缺点 §2 迭代法 2.1 迭代法的基本思想 迭代法举例 例4（续） 迭代法举例续 迭代法的几何含义 迭代法的几何含义（续） 2.2 迭代法的收敛条件（三大定理） 压缩映象原理的证明 压缩映象原理的证明（续1） 压缩映象原理的证明（续2） 两个重要误差公式说明 两个重要误差公式说明（续） 迭代法的收敛条件之二 定理8.2证明 定理2.2应用举例 定理2.2应用举例（续） 应用举例 迭代法的收敛条件之三 2.3 Steffensen方程——简单迭代法的加速 {xn}的r 阶收敛定理 定理2.4（续） 1、 Aitken加速法 Aitken加速法（续） 2、 Steffensen加速收敛式 Steffensen加速收敛举例 §3 Newton法与弦截法 Newton法（续） Newton法的几何意义 Newton法举例 Newton法收敛定理 定理8.5（续） Newton法的几何意义及其优劣 3.2计算重根的牛顿迭代法 计算重根的牛顿迭代法(续1） 计算重根的牛顿迭代法(续2） 计算重根的牛顿迭代法(续3） 例8（续） 3.2 弦截法 弦截法迭代公式的几何解释 弦截法的几点说明 弦截法的几点说明(续） 弦截法举例 表8-4 §4 抛物线法 抛物线法（续） 抛物线法的几点说明 抛物线法举例 第八章 结 束 定理8.3 定理8.3强调迭代初值x0应取在根x*的邻域中。如果对任意给定的x0，迭代格式均收敛，则称此格式具有全局收敛性，但这样的格式是极其稀少的。如果对根x*的某邻域内的任一点x0，迭代格式均收敛，则此格式具有局部收敛性。 即可保证对其中任取的一点x0迭代收敛。事实上，在用迭代法求解方程（8-1）时，常常先用对分区间求得较好的初值，然后再进行迭代。 本定理给出的就是局部收敛性条件。具体解题时，虽然无法判别隔根区间是否为以x*为中心的邻域，但只要它足够小，且在邻域中满足： 定理8.3 （续） 收敛速度（收敛速度的阶）： 成立，则称{xn}是r 阶收敛的，或称{xn}的收敛阶为r，收敛阶r 的大小刻划了序列{xn}的收敛速度： r 越大，收敛越快: r =1 线性收敛 r 1 超线性收敛 r =2 平方收敛 设序列{xn}收敛于x*，若存在正数r和a使得： 定理8.4 设迭代函数? (x)在x*邻近有r阶连续导数，且 x* = ?(x*)，并且有 证明： 1） ?(x)满足收敛定理的条件?{xn}?x*； 紧接下屏 2）利用Taylor公式将?(x)在x*附近展开： 这表明：{xn}是r阶收敛的。 一阶收敛即为线性收敛，收敛速度较慢，下面想法加速： 证毕！ 若序列{xn}线性收敛于x*，可按式： 当n充分大时，有： 紧接下屏 由此式可推导出： 由此可得比值： 将Aitken加速法与简单迭代格式xn+1 = ? (xn)相结合就得到 Steffensen加速收敛式 ： 当n ? ?时，比值中分子趋于0的速度比分母趋于0的速度快，亦即分子是比分母高阶的无穷小，这表明{xn}比{xn}更快地收敛于x*。 例6 的根，取x0 =1.5，误差精度? = 10?6。 于是由x0 = 1.5， 可计算： 继续下去 ，在此可 求x2, x3,… 由此例题可见：Steffensen方法收敛很快 ，达到了加快收敛的目的。 用Steffensen方法求方程： 3.1 Newton法 将非线性方程线性化，以线性方程的解逐步逼近非线性方程的解，这就是Newton法的基本思想。 设已知方程f (x) = 0的近似根x0，f (x)在其零点x*邻近一阶连续可微， 且f ?(x) ? 0，当x0充分接近x*时，f (x)可用Taylor公式近似表示为 ： 则方程f (x) = 0可用线性方程近似代替，即： 解此线性方程得： 取此x作为原方程的新近似值x1，重复以上步骤， 于是得迭代公式： 按式（8-7）求方程f (x) = 0近似解称为Newton法。 如此继续下去，xn+1为曲线上点(xn, f (xn))处的切线与x轴的交点。因此Newton法是用曲线的切线与x轴的交点作为曲线与x轴交点的近似，故