第7章 非线性方程与方程组的数值解法.ppt

数值分析 李小林 重庆师范大学数学学院 若令 则得到斯蒂芬森(Steffensen)迭代法的另一种形式： Steffensen迭代法的优点：可以改进收敛速度，有时也能把不收敛的迭代法改进为收敛的二阶方法. 埃特金（Aitken）加速方法： 例5： 已知方程 在 上有一个根（正根） 下面选取3种迭代格式： 1、 即 2、 即 3、 即 取 计算结果如下： 法2 原迭代次数29 法3 原来不收敛 法1 原来不收敛 设不动点迭代 的迭代函数 在其不动点 的某邻域内具有二阶连续导数， 则斯蒂芬森(Steffensen)的迭代技术是二阶收敛的，而 且其极限仍为 。 证明： 设斯蒂芬森(Steffensen)的迭代为 其中 首先证明有 相同的不动点 设 则 ? 反之，设 型 洛比塔法则 ? 其次证明Steffensen迭代是二阶收敛的 由Taylor公式： 在 处展开 则 上式中用 替换 即证 代入上述结果： ? 注意：对本来就是p(1)阶收敛的方法， 改用Stefensen迭代方法优点不多。 §4 牛顿法 /* Newton - Raphson Method */ 一、牛顿迭代公式的推导 1、待定参数法 不动点迭代的关键是构造满足收敛条件的迭代函数 一种自然的选择是令 为了加速不动点迭代的收敛过程，应尽可能使迭代函数 在 处有更多阶导数等于零。 令 现设 取 满足 ? 因此，选取迭代函数 Newton – Raphson迭代格式 称之为牛顿—拉夫森方法，简称牛顿法 原理：将非线性方程线性化 取 x0 ? x*，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间 2、Taylor展开法/* Taylor’s expansion Method */ 将 (x* ? x0)2 看成高阶小量，则有： x y x* x0 只要 f ?C1，每一步迭代都有 而且 ，则 x*就是 f 的根。 与x轴交点的横坐标 无开方运算，又无除法运算。 例1：?写出求 的Newton迭代格式； ?写出求 的Newton迭代格式,要求公式中既 解： ?等价于求方程 的正根 ?解法一： 等价于求方程 的正根 ? 解法二： 等价于求方程 的正根 ? （局部收敛性） 设 x* 为方程 f (x) =0的根，在包含x*的某个开区间内 连续， 且 ，则存在 x* 的邻域 ，使得任取初值 ，由Newton’s Method产生的序列 以不低于二阶的收敛速度收敛于x*，且 Th 证明：Newton’s Method 事实上是一种特殊的不动点迭代 其中 ，则 至少平方收敛 由 Taylor 展开： 在单根 /*simple root */ 附近收敛快 ? 只要 ，则令 可得结论。 有根 根唯一 产生的序列单调有界保证收敛 证明： 因为f ?C2[a, b]，由（1）和（2）知f (x) 在[a, b]内有唯一根 下面由条件（1）、（2）分4种情况讨论： ? ? ? ? 仅证明第一种情况，其它情况类似讨论 （收敛的充分条件）设f (x) =0 且f ?C2[a, b]，若 (1) f (a) f (b) 0；(2) 在整个[a, b]上 不变号且 ； (3) 选取 x0 ? [a, b] 使得 ； 则Newton’s Method产生的序列{ xk } 收敛于方程的根 , Th 且 由中值定理，