非方程的数值解法要素.ppt

由高等数学知识知, 设f (x)为区间[a,b]上的单值连续, 如果f (a)·f (b)0 , 则[a,b]中至少有一个实根。如果f (x)在[a,b]上还是单调地递增或递减，则仅有一个实根。 记y1=x , y2=φ(x) , 它们交点的横坐标α即为方程的根 Steffensen迭代格式 已知 ，则 ， 改成 Steffensen迭代格式 也可以改写成 其中迭代函数 Steffensen迭代法收敛的充要条件 定理7.4.1 Steffensen迭代法收敛的充要条件 证明：必要性 Steffensen迭代法收敛的充要条件 充分性 Steffensen算法的收敛速度 定理7.4.2 在定理7.4.2假设下，若 产生的序列 至少平方收敛到 。 Steffensen算法的收敛速度 Steffensen算法的收敛速度 Steffensen算法的收敛速度 Steffensen算法的收敛速度 由定理7.4.2知 至少以平方速度收敛到 。 也就是说：简单迭代法是线性收敛；Steffensen迭代至少平方以上收敛（加速收敛）。 例题 例7.9试用Steffensen算法求解方程 解法一、取 ，由 例题 取初值 ，计算结果如下： 例题 解法二、取 ，由 对于该迭代函数在一般迭代法中是发散的，而Steffensen格式却是收敛的。 例题 取初值 ，计算结果如下： Steffensen迭代格式几何解释 Steffensen迭代算法 Steffensen迭代算法 练习 练习1 用Newton法求 的近似解。 解：由零点定理。 练习 练习 练习2 用Newton法计算 。 解： 1.理解收敛性、收敛阶的概念及二分法思想方法。 2.会求用二分法解非线性方程时的执行次数k 。 为给定的误差界. 3.理解简单迭代法的思想方法，几何意义，压缩不动点定理。 4. 掌握简单迭代法的收敛(局部)定理（定理证明，会判断简单迭代法是否收敛）。 7.4.2 割线法与抛物线法 7.4.1 Newton迭代法 7.4 一元方程的常用迭代法 用迭代法可逐步精确方程 根的近似值，但必须要找到 的等价方程 ,如果 选得不合适,不仅影响收敛速度,而且有可能造成迭代格式发散。能否找到一种迭代方法,既结构简单,收敛速度快,又不存在发散的问题。这就是本节要介绍的牛顿迭代法. 1 牛顿迭代法的基本思想 牛顿迭代法一种重要和常用的迭代法, 它的基本思想是将非线性函数f(x)逐步线性化, 从而将非线性方程f(x)=0近似地转化为线性方程求解。 7.4.1 Newton迭代法 对于方程 ,设其近似根为 , 函数f (x)可在 附近作泰勒展开 忽略高次项,用其线性部分作为函数 f (x)的近似， 设 的根 ,则有 ,即 将右端取为 ,即 是比 更接近于 的近似值 这就是著名的牛顿迭代公式 （7.4.2） （7.4.1） 3 牛顿迭代法的几何解释 任取初始值 ， 上过点 的切线方程为： 与 轴交于点 过点 的切线方程为 与 轴交于点 如此下去得牛顿迭代公式: 用切线代替曲线，用 线性函数的零点作为 f(x)的零点的近似值。 牛顿迭代法也称切线法 因此牛顿法产生的序列{xk}如下图所示。 x0 x2 x1 过P0的切线 过P1的切线 将(7.4.2)写成一般的不动点迭代(7.3.3)的形式,有 所以有 ， Newton迭代法是超线性收敛的。更准确地，从(7.4.1)和(7.4.2)可得下面的定理. （收敛的充分条件）设 f ?C2[a, b]，若 (1) f (a) f (b) 0；(2) 在整个[a, b]上 f 〞不变号且 f′(x) ? 0； (3) 选取 x0 ? [a, b] 使得 f (x0) f 〞(x0) 0； 则Newt