9_YQJ_Ch4非线性方程求根Q.pptVIP

迭代过程的加速 对于收敛的迭代过程，从理论上说，只要迭代足够多次，总可得到满意的精度，但有时迭代过程过于缓慢，计算量极大，实际上得不到满意的结果。因此，迭代过程的加速便成了一个重要的研究课题。 ? 加速方案一 ? 加速方案二 ? Steffensen 加速：把 xk+1, xk+2 看作中间值, 分别记为 yk , zk 定理4.5 设 s为g x 的不动点, g x 在含 s的某个开区间内有连续的二阶导数, 且 g s ≠1, 则存在 ? 0, 对于 x0?B? [s- ?, s+ ?] 但 x0 ≠ s时，由 Steffensen 迭代 4.6 产生的序列 xk 至少以二阶速度收敛于 s. 证 构造函数 ψ x , x1 2.座机电话号码845905 x2 2.0座机电话号码69250 x3 1.0座机电话号码73631 x4 0.座机电话号码673049 x5 0.1座机电话号码04091 x6 0.1座机电话号码75348 x7 0.1座机电话号码79268 x8 0.1座机电话号码48468 x9 0.1座机电话号码67211 x10 0.1座机电话号码72488 * 当迭代函数 g x 在不动点 s 处导数不为零, 迭代 xk+1 g xk 仅是线性收敛的. 现在介绍 加速法. 设由迭代公式计算出 迭代值 ，若 x 在 s 附近变化, g x 连续, 则 g x 变化不大, 近似地记为常数L，则用微分中值定理，可以得到 4.1.4 Steffensen（斯蒂芬森）迭代法 把它作为 的一种补偿 解出s： 即记： 那么x1 应比 近似得更好。 具体算法如下： 例 用简单迭代法求方程 x e-x 在x 0.5附近的一个根, 要求 解：方程 x e-x ,迭代函数为 g x e-x , 在区间 [0.5, 0.6] 内有根, 根据定理4.2知, 在 [0.5, 0.6] 内 g x 有不动点 s; 且迭代公式 xk+1 e-xk产生的序列 xk 局部收敛于 s. 计算结果: f 0.5 0.5-e-0.5 0.5-0.6065 0 f 0.6 0.6-e-0.6 0.6-0.5488 0 function y JDDDL31 x0,epsilon ek 1;xk x0; %n 1; for n 1:16 yk g xk ; xk yk; ek abs yk-xk /abs xk ; fprintf \n y%d %15.14f \n,n,yk end fprintf \n ek |yk-xk|/|xk| %f \n,ek y yk; %%%%%%%%%%%%%%% function Y g x Y exp -x ; x0 0.5 x0 0.5000 epsilon 0epsilon 1.0000e-008 y JDDDL31 x0,epsilon y1 0.座机电话号码971263 y2 0.座机电话号码189261 y3 0.座机电话号码487807 y4 0.座机电话号码793890 y5 0.座机电话号码897722 y6 0.座机电话号码698032 y7 0.座机电话号码757007 y8 0.座机电话号码274692 y9 0.座机电话号码426224 y10 0.座机电话号码293547 y11 0.座机电话号码597078 y12 0.座机电话号码185373 y13 0.座机电话号码008764 y14 0.座机电话号码425699 y15 0.座机电话号码370764 y16 0.座机电话号码365927 ek |yk-xk|/|xk| 0.000000 y 0.5671 例 用加速方案一求x e-x 在x 0.5附近的一个根, 解: 取L -0.6，则加速公式为： 0.56715 0.56714 3 0.56746 0.56712 2 0.60653 0.56658 1 0.50000 0 xk k 计算结果: 当迭代函数 g x 在不动点 s 处导数不为零, 迭代 xk+1 g xk 4.3 仅是线性收敛的. 现在介绍 Steffensen 加速法. 设由 4.3 计算出 3 个迭代值 xk , xk+1, xk+2 , 则用微分中值定理, 可有 若 x 在 s 附近变化, g