4.1 解非线性方程及迭代解法.pptVIP

4.1.5　Newton迭代法 第四章 非线性方程(组)的迭代解法 §4.1非线性方程的迭代解法 §4.1.1 对分法 f (x)在区间[a,b]上只有一根， 且 设 则可用对分法求方程 的近 似根。 算法4.1(对分法) 输入 及精度水平 和 对于k =0,1,…执行 (1) 计算 (2) 若 或 则停止计算， 取 否则执行(3)； (3) 若 则令 若 则令 解 例1 §4.1.2 简单迭代法及其收敛性 简单迭代法的基本步骤： 将方程 化为等价方程 选定方程的解 s 的一个初始近似值 按递推式 产生序列 其极限必为方程(4.2)的解s， (4.2) (4.3) 收敛， 若 则当k 足够大时可用 作为根的近似值。 称 为迭代函数， 为 的不动点； 简单迭代法(4.3)又称为不动点迭代法。 迭代过程的几何表示 例2 构造求解方程 的迭代公式。 解： 由此可构造出如下多种简单迭代法: 为求同一个根，它们所产生的迭代序列{xk},有的收敛 , 有的发散，有的收敛的快，有的收敛的慢。 可将方程等价地改写成以下几种形式 (4.4) (4.5) (4.6) (4.7) 1.324717957 0.000038542 -8.494e+231 14 1.324717957 673222706.25 -∞ 15 … … … -2.040e+77 … 1.324717957 1.324701749 3.959001389 - 5.88e+25 6 1.324717957 1.324632625 0.644502197 -730.0000000 5 1.324718174 1.324268745 2.243252653 -9.000000000 4 1.325200399 1.322353819 0.926780890 -2.000000000 3 1.347826087 1.312293837 1.710000000 -1.000000000 2 1.500000000 1.259921050 1.111111111 0.000000000 1 1.000000000 1.000000000 1.500000000 1.000000000 0 k 计算结果 收敛定理 定理4.1 设函数 在(a,b)内可导，满足 (1) 当 时， (2) 当 时， 其中L为常数。 则有如下结论 产生的序列 收敛于 s ； ③ 成立误差估计式 ② 则对任取的 由迭代公式 ① 方程 在区间[a,b]上有唯一的根 s ； (4.8) (4.9) 定理4.2 设 在包含s的某个开区间内连 续。如果 则存在 当 时，由简单迭代法 产生的序 列 列 收敛。 例3 用简单迭代法求方程 在区间 内的根，要求 解 记 由于 在区间 内严格单调递增，因而原方程 在区间 内有且仅有一根 s,且 将方程化为以下同解方程 则 由定理4.1，对任意的 迭代公式 产生的序列 必收敛于s。 计算结果 其中： 所以，取s≈3.146193215作为方程的近似根。 另一根的计算 其中： 所以，取s≈0为方程的近似根。 §4.1.3 简单迭代法的收敛速度 定义4.1 设序列 收敛于s，并且 若存在常数 和常数c0,使 成立，或者存在正整数K,当kK时，有 成立，则称序列 收敛于 s 具有 r 阶收敛速度， 简称序列 是 r 阶收敛的。 线性收敛： 超线性收敛： 平方收敛： 定理4.3 设函数 且满足 (1) 当 时， (2) 当 时， 其中L为常数。 则对任取的 由迭代公式 产生的序 列 收敛于方程 在[a,b]内的 唯一根 s，且当 时， 是线性收敛的。 定理4.4 设 在包含s的某个开区间内连 续(m1)。如果 则存在 当 代法 产生的序列 的速度收敛于 s 。 以 m 阶 时，由简单迭 例4 对于方程 和 都收敛于方程的唯一解s，且 s 0。 对于(4.6) ， 的导数 在 s 附近连续，且 由例2，我们已经知道 迭代公式 所以，当x0充分 接近于s 时，迭代公式(4.6)产生的序列{xk}线性收敛于s。 (4.6) (4.7) 迭代函数 对于(4.7) ， 它的导数 在s附近连续，且 迭代公式(4.7)产生的序列{xk}至少以二阶收敛速度收敛于s。 迭代函数为 所以，只要x0充分接近于s ， §4.1.4 Steffensen加速收敛方法 设迭代公式 产生的序列{xk}线性收 敛于s。则由微分中值定理得 其中 介于 xk+1 和 s 之间， 介于 xk+2 和 s 之间，假定当 在s附近变化时， 变化不大， 由此解出 则有 上式右端的值可能比xk+1、xk+2更接近于