数值分析第7章4 6节.pptVIP

若取 f ( x ) ? ( x ) ? x ? m , f ? ( x ) 则 ? ? ( x *) . ? 0 用迭代法 x k ? 1 f ( x k ) ? x k ? m f ? ( x k ) ( k ? 0 , 1 , ? ) （ 4.13 ） 求 m 重根，则具有 2 阶收敛，但要知道 x 的重数 . m * 构造求重根的迭代法，还可令 若 x 是 * 重根，则 f ( x 的 ) ? 0 m ( x ? x *) g ( x ) ? ( x ) ? , mg ( x ) ? ( x ? x *) g ? ( x ) 21 ， ? ( x ) ? f ( x ) / f ? ( x ) 故 x 是 * ? ( x 的单根 ) ? 0 . 对它用牛顿法，其迭代函数为 ? ( x ) f ( x ) f ? ( x ) ? ( x ) ? x ? ? x ? . 2 ? ? ( x ) [ f ? ( x )] ? f ( x ) f ? ? ( x ) 从而可构造迭代法 x k ? 1 f ( x k ) f ? ( x k ) ? x k ? 2 ? [ f ( x k )] ? f ( x k ) f ? ? ( x k ) ( k ? 0 , 1 , ? ), （ 4.14 ） 它是 2 阶收敛的 . 22 例 9 方程 x 4 ? 4 x 2 ? 的根 4 ? 0 是二重根， x * ? 2 用上述三种方法求根 . 解 f ( x ) ? x ? 4 x ? 4 , 先求出三种方法的迭代公式： 4 2 （ 1 ） 牛顿法 x ? 4 x ? 4 x ? 2 x k ? 1 ? x k ? ? x k ? . 3 4 x k ? 8 x k 4 x k x ? 2 m ? 2 , x k ? 1 ? x k ? . 2 x k 2 x k ( x k ? 2 ) x k ? 1 ? x k ? . 2 ( k ? 0 , 1 , ? ) x k ? （ 2 4.13 ） 2 k 4 k 2 k 2 k （ 2 ） 用（ 4.13 ）式 （ 3 ） 用（ 4.14 ）式 f ( x k ) x k ? 1 ? x k ? m f ? ( x k ) 取初值 x 0 ? ，计算结果如表 7-7. 1 . 5 f ( x k ) f ? ( x k ) x k ? 1 ? x k ? [ f ? ( x k )] ? f ( x k ) f ? ? ( x k ) 2 ( k ? 0 , 1 , ? ), 23 表7 ? 7 三种方法数值结果 k 1 2 3 x k x 1 x 2 x 3 方法（1） 1 . 458333333 1 . 436607143 1 . 425497619 方法（2） 1 . 416666667 1 . 414215686 1 . 414213562 方法（3） 1 . 411764706 1 . 414211438 1 . 414213562 从结果看出，经过三步计算，方法（ 2 ）及（ 3 ）均达到 10 位有效数字，而由于牛顿法只有线性收敛，所以要达到同 样精度需迭代 30 次 . 24 7.5 弦截法与抛物线法 用牛顿法求方程 还要算 f ? ( x . k ) 的根，每步除计算 f ( x ) ? 0 外， f ( x k ) 当函数 f 比较复杂时，计算 ( x ) 往往较困难， f ? ( x ) 为此可以利用已求函数值 f ( x k ), f ( x 来回避导数 k ? 1 ), ? 值 f ? ( x 的计算 . k ) 25 7.5.1 弦截法 设 x k 是 , x k ? 1 的近似根， f ( x ) ? 0 利用 f ( x k ), f ( x k ? 1 ) p 1 ( 的根作为新的 x ) ? 0 构造一次插值多项式 p 1 ( ，并用 x ) 近似根 x k ? . 1 由于 f ( x k ) ? f ? x k ? 1 ? p 1 ( x ) ? f ( x k ) ? ( x ? x k ) ， x k ? x k ? 1 （ 5.1 ） 因此有 x k ? 1 f ( x k ) ? x k ? ( x k ? x k ? 1 ) ， f ( x k ) ? f ( x k ? 1 ) （ 5.2 ） 26 （ 5.2 ）可以看作将牛顿公式 x k ? 1 f ( x k ) ? x k ? f ? ( x k ) f ( x k ) 取代的结果 ? f ? x k ? 1 ? . x k ? x k ? 1 中的导数 f ? ( x 用差商 k ) 接着讨论几何意义 . 曲线 y ? 上横坐标为 f ( x )