第4章解非线性方程-2概念.ppt

取初始值(1,1,1), 计算如下 N x y z 0 1.0000000 1.0000000 12.1893260 1.5984751 1.3939006 1.8505896 1.4442514 1.2782240 1.7801611 1.4244359 1.2392924 1.7776747 1.4239609 1.2374738 1.7776719 1.4239605 1.2374711 1.7776719 1.4239605 1.2374711 练习： 3. Newton 迭代法是如何推出的? 它若在单根附近收敛,是几阶收敛?在重根附近是几阶收敛？求方程重根时,能达到2阶收敛的改进 Newton 迭代公式是什么 用牛顿法求方程 在区间 [1，2] 内 的一个实根，要求 2. 导出求立方根 的迭代公式，并讨论其收敛性。 首先导出求根方程 ，再对 使用牛顿法 得迭代公式 ,用全局收敛性定理或局部收 敛性定理讨论其收敛性。 §4.4 非线性方程的牛顿法 (Newton Method of Nonlinear Equations ) 内容提纲（Outline) 牛顿法及其几何意义 收敛性及其收敛速度 计算实例及其程序演示 取x0作为初始近似值，将 f (x)在x0做Taylor展开: 重复上述过程 ? 作为第一次近似值 一、牛顿法及其几何意义 Newton 迭代公式 基本思路：将非线性方程 f (x)=0 线性化 牛顿法的几何意义 x y x* x0 x 1 x 2 牛顿法也称为切线法, 迭代函数 （局部收敛性定理） 设 f (x)?C2[a, b]，若 x* 为 f (x) 在[a, b]上的根,且 f ?(x*) ? 0，则存在 x* 的邻域 使得任取初始值 ，Newton 法产生的序列{ xk } 收敛到 x*，且满足 至少平方收敛 二、牛顿法的收敛性与收敛速度 定理4.4.1 在x*的附近收敛 由Taylor 展开： 令k?? ，由 f ?(x*) ? 0，即可得结论。 证明：Newton法实际上是一种特殊的迭代法 迭代函数为： 思考题1 若　　　 ，Newton法是否仍收敛？ 设 x* 是 f 的 m 重根，则令： 且 Answer1: 有局部收敛性 Answer2: 线性收敛 思考题2 当x* 是 f (x)=0的m重根, 是否平方收敛？ 结论：Newton法的收敛性依赖于x0 的选取。 局部收敛定理对初始值 x0 要求较高。 x* x0 ? x0 ? x0 ? 有根 根唯一 (全局收敛性定理): 设 f (x)?C2[a, b], 若 f (a) f (b) 0； 在整个[a, b]上 f ?(x) ? 0； f ??(x)在 [a, b]上不变号 选取初始值x0 ? [a, b] 使得 f ??(x0) f (x0) 0； 则由Newton法产生的序列{ xk } 单调地收敛到 f (x)=0 在 [a, b] 的唯一根x*,且收敛速度至少是二阶的 保证产生的序列{xk}单调有界 保证Newton迭代函数将[a , b]映射于自身 定理4.4.2 将f (x* )在 xk 处作Taylor展开 对迭代公式两边取极限，得 证明：以 为例证明 说明数列{ xk }有下界 故{xk}单调递减, 从而{xk}收敛.令 ？ 定理4.4.3 (全局收敛性定理): 设 f (x)?C2[a, b], 若 (1) f (a) f (b) 0； (2) 在整个[a, b]上 f ?(x) ? 0, f ??(x) ? 0 ； (3) 则对任何 , Newton 迭代格式产生的序列 都收敛于f (x)=0 的根 x* . 注:定理的条件(3) 保证了从 x*两侧任取 x0 , 所得到的数列{ xk }均在[a , b]内. 三、计算实例及其程序演示 辅助工具: VC程序设计语言 Matlab数学软件 (1) 选定初值x0 ,计算f (x0) , f ?(x0) 计算步骤 (2) 按公式