数值分析第7章—非线性方程与方程组的数值解法.ppt

只是收敛速度将大大减慢。 1、当 为单根时，牛顿迭代法在根 的附近 是二阶收敛的； 2、当 为重根时，设为m重根，则 f (x)可表为 其中 此时用牛顿迭代法求 仍然收敛， 证方程 之解存在且唯一. 由于 在[a , b]上存在, f (x) 在[a , b]上连续。 作函数 由条件 连续。 所以 证 使 即 则 (1) f (a) ≤0 , f (b) ≥0 , 故存在 , 则由微分中值定理及条件值定理及条件(2)有 此式仅当 才能成立， 因此 则由微分中值定理及条件(2) 有 设方程 还有一不动点 定理 2 (1) 当x∈[a , b]时， (2) 存在正数L1，使对任意的 x∈[a , b], 对任意迭代初值 x0∈[a , b],迭代序列 ,且满足条件： 设 收敛于 。 则 且有下列误差估计式 即迭代过程收敛， 且 反复用此不等式，并注意 0 L 1 , 因此 先证迭代格式 收敛 任取 x0∈[ a, b ]，由微分中值定理，有 提示：定理的证明利用定理1以及微分中值定理。 则任取 x0∈ U , 迭代格式 均收敛于 , 定理 3 若方程 之根的某邻域 L 1，使 上面的定理在理论上十分重要, 但是条件(1)却不容易判别. 如果仅在根的邻域中考察迭代格式, 则下述定理可避免条件(1)的判别。 内 存在,且存在正常数 则迭代必发散。 反之，若在根 的邻域 U 内 称为局部收敛定理 证明 例1中采用的三种迭代格式 ,在有根区间 (1, 1.2)内有 例如 迭代收敛。 迭代收敛。 迭代发散。 例3 考察例2中四中迭代法在根附近的收敛情况，取根的近似值为 解 不收敛 不收敛 收敛 收敛 上例说明 值越小,收敛速度就越快. 四、迭代法的收敛速度 则称迭代格式 是 p 阶收敛的. p = 1时称为线性收敛, p1时称为超线性收敛. 利用微分中值定理及泰勒展式可得下面的定理3. 显然,收敛阶越大,收敛越快. p = 2 时称为二阶(平方)收敛, 特别地, 令 若 则迭代过程在 的邻近为 p 阶收敛。 (1) 若 为线性收敛； 则迭代过程在 的邻近 (2) 若 定理 4 之根,在 的邻域 U内 有连续的 p 阶导数，则 设 为 由泰勒展式可得 解 的三阶方法。假设 x0 充分靠近 , 求 证明迭代公式 xk+1=xk(xk2+3a)/(3xk2+a) ，试求 例 2 加速迭代法 松弛法迭代公式： 为松弛因子 7.3 迭代收敛的加速方法 称为斯蒂芬森 ( Steffensen ) 加速法. 则埃特金法为平方收敛； 若 为 p ( p 1)阶收敛， 导数连续, 的 p 阶 可以证明: 若 为线性收敛, 则埃特金法为 2p –1 阶收敛。 迭代格式 求方程 x = e –x 在 x = 0.5 附近的根. x25 = x26 = 0.5671433 若对此格式用斯蒂芬森法, 则 取 x0 = 0.5, 迭代格式 得 解 例3 仍取 x0 = 0.5 , 得 由此可见,斯蒂芬森法加速收敛效果是相当显著的. 例4 分别用松弛法、斯蒂芬森法求方程 在初值 附近的一个根，取迭代格式 解 用松弛法计算，取 因此松弛法的迭代公式为 列表计算如下： 1.365230013 1.365230012 1.364953916 1.5 0.887130869 0.887130869 0.890803686 3 2 1 0 n 用斯蒂芬森方法计算，迭代格式为： 列表计算如下： 1.365230583 1.367376372 1.365225534 1.348399725 1.365230013 1.365265224 1.5 2 1 0 n 7.4 牛顿法 一、牛顿法的基本思想 设已知方程 f (x) = 0 的近似根 x0，且在 x0附近 可用一阶泰勒多项式近似，表示为 方程 f (x) = 0 可用线性方程近似代替，即 解此线性方程得 取此 x 作为原方程的新近似根 x1，重复以上步骤， 得迭代公式 此式称为牛顿(Newton)迭代公式。 例 1 用牛顿法求方程 在 内一个实根，取初始近似值 解 所以迭代公式为： 列表计算如下： 111.3733333 1.5