第二章非线性程的数值解法2.pptVIP

§4 牛顿法 一、牛顿迭代公式的推导(Taylor展开法) 思想：非线性方程线性化（以直代曲） 取 x0 ? x*，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间 将 (x* ? x0)2 看成高阶小量，则有： 令 注：牛顿法是一种特殊的迭代法。 Newton – Raphson迭代格式 称之为牛顿—拉夫森方法，简称牛顿法. x y x* x0 与x轴交点的横坐标 无开方运算，又无除法运算。 例1：?写出求 的Newton迭代格式； ?写出求 的Newton迭代格式,要求公式中既 解： ?等价于求方程 的正根 ?解法一： 等价于求方程 的正根 ? 解法二： 等价于求方程 的正根 ? Th2.7 （局部收敛性） 设 x* 为方程 f (x) =0的根，在包含x*的某个开区间内 连续， 且 ，则存在 x* 的邻域 ，使得任取初值 ，由Newton’s Method产生的序列 以不低于二阶的收敛速度收敛于x*，且 证明：Newton’s Method 事实上是一种特殊的不动点迭代 其中 ，则 收敛 由 Taylor 展开： ? 只要 ，则令 可得结论。 Th2.5 有根 根唯一 产生的序列单调有界保证收敛 证明省略。 Th2.8 （收敛的充分条件）设f (x) =0 且f ?C2[a, b]，若 (1) f (a) f (b) 0；(2) 在整个[a, b]上 不变号且 ； (3) 选取 x0 ? [a, b] 使得 ； 则Newton’s Method产生的序列{ xk } 收敛于方程的根 . Th2.9 （收敛的另一充分条件）设 在[a, b]上连续， (1) f (a) f (b) 0；(2) 在整个[a, b]上 且 ； (3) ， 则对 ，Newton’s Method产生的序列{ xk } 收敛于方 程 在[a, b]内的唯一实根 。 且 证明省略。 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 改进与推广（补充） ? 重根 加速收敛法： Q1: 若　　　 ，Newton’s Method 是否仍收敛？ 设 x* 是 f 的 n 重根，则： 且 。 因为 Newton’s Method 事实上是一种特殊的不动点迭代,其 中 ，则 A1: 有局部收敛性，但重数 n 越高，收敛越慢。 Q2: 如何加速重根情况时的收敛速度？ A2: 将求 f 的重根转化为求另一函数的单根。 ? 令　　　　　，则 f 的重根 = ? 的单根。 ? 求复根 —— Newton 公式中的自变量可以是复数 记 z = x + i y, z0 为初值，同样有 设 代入公式，令实、虚部对应相等，可得 §5 弦割法与抛物线法 x0 x1 割线 切线斜率　?　割线斜率 需要2个初值 x0 和 x1。 Newton’s Method每一步要计算 ，为了避免计算导数值，现用 的近似，得到弦割法（割线法）。 x2 一、弦割法 Th2.10 局部收敛性 设 表示区间 ， x*为方程 f (x) =0的根， 函数f (x)在 中有足够阶连续导数， 且 满足