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

第二章 非线性方程求解 第二章 非线性方程求解目录 第二章 非线性方程求解概述 求根的隔离区间的两种方法 §1 对分法 对分法的误差估计 f=x^3+10*x-20 f = x^3+10*x-20 double(solve(f)) ans = 1.5946 -0.7973 + 3.4506i -0.7973 - 3.4506i 对分法的优缺点 §2 简单迭代法 2.1 迭代法的基本思想 迭代法举例续 迭代法的几何含义 迭代法的几何含义（续） 压缩映象原理的证明 2.2 迭代法的收敛条件（三大定理） 压缩映象原理的证明 压缩映象原理的证明（续1） 压缩映象原理的证明 两个重要误差公式说明 迭代法的收敛条件之二 证明 推论1应用举例 应用举例 迭代法的收敛条件之三 2.3 Steffensen方程—简单迭代法的加速 {xn}的r 阶收敛定理 1、 Aitken加速法 Aitken加速法（续） §3 Newton法与弦截法 Newton法（续） Newton法的几何意义 Newton法举例 Newton法收敛定理 Newton法的几何意义及其优劣 3.2计算重根的牛顿迭代法 计算重根的牛顿迭代法(续1） 计算重根的牛顿迭代法(续2） 3.2 弦截法 弦截法迭代公式的几何解释 弦截法的几点说明 弦截法的几点说明(续） 则 由式（2-8）可知 故 因此可用 下式估计m 例8 用牛顿迭代法求方程 在0.95附近之根。 解 取x0=0.95，用牛顿迭代法， 按式（2-7）求得的xk见表3-3，由表中数据可见xk 收敛很慢。由 可知，所求根为m=2重根， 改用式(2-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法优点：收敛快（平方阶），固定格式； 修 正：以差商代替导数（微商） 设x, y为[a, b]上的任意两点，由微分中值定理，在 x, y之间至少存在一点?，使得: 于是： 即? (x)满足定理2.6的条件（2），故结论成立。 证毕！ 采用的三种迭代格式， 在隔根区间（1,1.2）内有： 用推论1判别简单迭代法的收敛性比定理2.6方便 如对例题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阶收敛的。 一阶收敛即为线性收敛，收敛速度较慢，下面想法加速： 证毕！ 若序