数学建模数值分析模型.pptVIP

6.3 迭代法 黑 龙 江 科 技 学 院 数 学 建 模 理学院 将方程 ( ) 0 f x ? 等价变形为 ( ) x x ? ? ，若要求满足 * ( ) 0 f x ? 的根 * x ，等价的是求 * x 使得 * * ( ) x x ? ? ，称 ( ) 0 f x ? 与 ( ) x x ? ? 同解；反之亦然。这时的 * x 称为是函数 ( ) x ? 的 一个不动点 。求方程 ( ) 0 f x ? 的根等价于求 ( ) x ? 的不动点。 不动点迭代关系式 ( 也称简单迭代法 ) 为 1 ( ), 0,1,2,... k k x x k ? ? ? ? (6.3.1) 其中函数 ( ) x ? 称为迭代函数 . 如果对任意 1 ( ), 0,1,2,... k k x x k ? ? ? ? 由式（ 6.3.1 ）产生的序列 { } n x 有极限 * lim k k x x ?? ? 则称 不动点迭代法（ 6.3.1 ）收敛 . 黑 龙 江 科 技 学 院 数 学 建 模 1 、 简单迭代法 x y o ) ( x y ? ? x y ? 是否对于任意的等价形式 该迭代法 都是收敛的？什么情况下收敛？ ) ( x x ? ? 黑 龙 江 科 技 学 院 数 学 建 模 理学院 x y y = x x y y = x x y y = x x y y = x x* x* x* x* x 0 p 0 x 1 p 1 ? x 0 p 0 x 1 p 1 ? x 0 p 0 x 1 p 1 ? x 0 p 0 x 1 p 1 ? ) ( x y ? ? ) ( x y ? ? ) ( x y ? ? ) ( x y ? ? 黑 龙 江 科 技 学 院 数 学 建 模 理学院 定理 6.3.1 ( 不动点存在性定理 ) 设 ( ) [ , ] x C a b ? ? 满足以下两个条件 : （ 1 ）对任意 [ , ] x a b ? ，有 ( ) [ , ] x a b ? ? （ 2 ）存在正常数 1 L 使对任意 , [ , ] x y a b ? 都有 ( ) ( ) x y y x ? ? ? ? ? 则 ( ) x ? 在 [ , ] a b 上存在惟一的不动点 * x 定理 6.3.2 （不动点迭代法的全局收敛性定理） 设 ( ) [ , ] x C a b ? ? 满足定理 6.3.1 中的两个条件 , 则对任意 0 [ , ] x a b ? 得到的迭代序列 由 (6.3.1) 式 { } n x 收敛到 ( ) x ? 的不动点，并有 误差估计 式 * 1 1 k k k L x x x x L ? ? ? ? ? * 1 0 1 k k L x x x x L ? ? ? ? 和 黑 龙 江 科 技 学 院 数 学 建 模 理学院 定理 6.3.3 （不动点迭代法的局部收敛性定理） 设 * x 为 ( ) x ? 的不动点 , ( ) x ? 在 * x 的某个邻域连续，且 ( ) 1 x ? ，则迭代法 (6.3.1) 局部收敛 . 黑 龙 江 科 技 学 院 数 学 建 模 例 6.3.2 求方程 0 1 3 3 ? ? ? x x 要求结果精确到 10-5 。 在 [0, 0.5] 内的根， 解 将方程变形 ) ( ) 1 ( 3 1 3 x x x ? ? ? ? , 因为 [0,0.5] x ? 0 ) ( 2 ? ? x x ? 在 [0, 0.5] 内为增函数，所以 1 25 . 0 5 . 0 ) ( max 2 ? ? ? x L ? 满足收敛条件，取 x 0 = 0.25 ，用迭代公式 ) 1 ( 3 1 3 1 ? ? ? n n x x 计算步骤如下 x 1 = ? (0.25) = 0.3385416 x 2 = ? ( x 1 ) = 0.3462668 x 3 = ? ( x 2 ) =0.3471725 x 4 = ? ( x 3 ) =0.3472814 x 5 = ? ( x 4 ) =0.3472945 x 6 = ? ( x 5 ) =0.3472961 x 7 = ? ( x 6 ) =0.3472963 取方程 0 1 3 3 ? ? ? x x 的近似根为 x * = 0.347296 。 黑 龙 江 科 技 学 院 数 学 建 模 2 、牛顿迭代法 基本思想 ：牛顿法是将非线性方程线性化的一种近似 方法。它是将 在初始点 附近展开成泰勒级数 ，取其线性部分 ( 即前两项 ) ，作为非线性方程 的近似方程，则 有： 。 ) ( x f 0 x ? ? ? ? ? ? ? ? ? ? ) ( ! 2 ) ( )