数值分析李庆扬第七章.pptVIP

( p ) 定理 4: 定理：对于迭代过程 x k ? 1 ? ? ( x k ), 如果 ? ( x ) 在所 求根 x 的邻近连续，并且 ? ( x ) ? ? ( x ) ? L ? ? * * * ( p ? 1) * ( x ) ? 0; ? * ( p ) ( x ) ? 0 * 则该迭代过程在点 x 邻近是 p 阶收敛的。 证：由于 ? ( x ) ? 0 ? 1, 故 x k ? 1 ? ? ( x k ) 具有局部收敛性。 将 ? ( x k ) 在根 x 处展开，由条件有 ? ( x k ) ? ? ( x ) ? ? x k ? 1 ? x ? * * * * ? ( p ) ( ? ) * p ( x k ? x ) p ! ( p ) * ? ( p ) e k ? 1 ( ? ) ? ( x ) * p ( x k ? x ) ? ? p p ! e k p ! § 3. 迭代收敛的加速方法 埃特金 (Aitken) 方法 设 x 0 是根 x 的某个预测值，通过两次迭代校正 有 x 1 ? ? ( x 0 ) x 2 ? ? ( x 1 ) 由微分中值定理，有 x 1 ? x ? ? ( x 0 ) ? ? ( x ) ? ? ( ? 1 )( x 0 ? x ) x 2 ? x ? ? ( x 1 ) ? ? ( x ) ? ? ( ? 2 )( x 1 ? x ) 假定 ? ( x ) 改变不大，近似地取某个近似值 L , 则由 x 1 ? x ? L ( x 0 ? x ) x 0 ? x x 1 ? x ? ? * * x 2 ? x x 1 ? x * * * * * * * * * * * x 2 ? x ? L ( x 1 ? x ) ( x 2 ? x 1 ) x ? x 2 ? x 0 ? 2 x 1 ? x 2 * 2 * * ? 此种加速需用两次迭代值进行加工。 (P273) 如果将一次改进值作为一步， 则计算公式如下： % 校正 x k ? 1 ? ? ( x k ) % 再校正 x k ? 1 ? ? ( x k ? 1 ) 改进 x k ? 1 ? x k ? 1 即为斯蒂芬森 (Steffensen) 迭代法 2 % ( x k ? 1 ? x k ? 1 ) ? % x k ? 1 ? 2 x k ? 1 ? x k (P273) 可以证明 表明 x k ? 1 ? x lim ? 0 * k ? ? x ? x k * x k 的收敛速度比 x k 的收敛速度快 , 且为 2 阶收敛的． § 4. 牛顿 (Newton) 法 非线性问题的最简单解法是线性近似 . 将非线性方程线性化，以线性方程的 解逐步逼近非线性方程的解，这就是 Newton 法的基本思想 f ( x ) ? 0 在真根附近 x 0 点展开成 Taylor 级数： f(x) ? f( x 0 ) ? (x - x 0 ) f ( x 0 ) ? ? x ? x 0 ? 2 f ( x 2 ! ) 0 ? L 取线性部分近似代替 f(x) ? 0 : f(x ) ? f (x ) ? x ? x ? ? 0 0 0 0 解出 x 作为近似根 x ： 1 x ? x ? f( x ) f ( x ) ( f ( x ) ? 0 ) 1 0 0 0 0 依次产生迭代格式，称 Newton 法： x n ? 1 ? x n ? f ( x n ) f ( x ), n ? 0 , 1 , 2 , ? n Newton 法的几何解释 当 x 0 在取后（在真根 ? 附近），过 ? x 0 , f ( x 0 ) ? 作 f ( x ) 的切线，则切线方程： y ? f ( x 0 ) ? f ( x )( x ? x ) 0 0 y ? f ( x ) x x 0 * x 1 x 2 1 与 x 轴交点 ( y ? 0 ) 为第二个近似根 x ： x ? 1 x ? f ( x ) f ( x ) 0 1 0 故也称切线法 迭代过程的收敛速度依赖于迭代函数 ? ( x ) 的选 取。如果当 x ? [ a , b ] 时 ?