数值分析--第7章德非线性方程求根.ppt

课件 课件 取x0=2, 对上式4种迭代法, 计算三步所得结果入下表. 2 1.75 1.732143 1.732051 ┆ 2 1.75 1.73475 1.732361 ┆ 2 1.5 2 1.5 ┆ 2 3 9 87 ┆ x0 x1 x2 x3 ┆ 0 1 2 3 ┆ 迭代法(4) 迭代法(3) 迭代法(2) 迭代法(1) xk k 课件 注意 ，从计算结果看到迭代法(1)及(2)均不收敛，且它们均不满足定理3中的局部收敛条件，迭代法(3)和(4)均满足局部收敛条件，且迭代法(4)比(3)收敛快，因在迭代法(4)中??(x*)=0. 为了衡量迭代法(2.2)收敛速度的快慢可给出以下定义. 定义2 设迭代过程xk+1=?(xk)收敛于方程x=?(x)的根x*，如果迭代误差ek=xk-x*当k→∞时成立下列渐近关系式 则称该迭代法是p阶收敛的. 特别地，p=1时称线性收敛，p1时称超线性收敛，p=2时称平方收敛. 课件 定理4 对于迭代过程xk+1=?(xk)，如果?(p)(x)在所求根x*的邻近连续，并且 则该迭代过程在x*的邻近是p阶收敛的. 证明 由于??(x*)=0，根据定理3立即可以断定迭代过程xk+1=?(xk)具有局部收敛性. 再将?(xk)在根x*处做泰勒展开, 利用条件(2.4), 则有 注意到?(xk)=xk+1，?(x*)= x*，由上式得 课件 因此对迭代误差，令k→∞时有 这表明迭代过程xk+1=?(xk)确实为p阶收敛. 证毕. 上述定理告诉我们，迭代过程的收敛速度依赖于迭代函数?(x)的选取. 如果x∈[a, b]但??(x)≠0时，则该迭代过程只可能是线性收敛. 对例4的讨论见书p272. 课件 的三阶方法. 假设 x0 充分靠近 x*, 求 证明 首先由泰勒展式可得 例子 证明迭代公式 xk+1=xk(xk2+3a)/(3xk2+a)是求 而1/4a≠0，故此迭代公式是三阶方法. 课件 7.3 迭代收敛的加速方法 7.3.1 埃特金加速收敛方法 对于收敛的迭代过程，只要迭代足够多次，就可以使结果达到任意的精度，但是有时迭代过程收敛较慢，从而使计算量变得很大，因此迭代过程的加速是个重要的课题. 　 设x0是根x*的某个近似值, 用迭代公式校正一次得 　　　　　 x1=?(x0) 而由微分中值定理，有 课件 　　假设??(x)改变不大, 近似地取某个近似值L, 则有 由于 x2-x*≈L(x1-x*). 　 若将校正值x1=?(x0)再校正一次，又得 　　　　　 x2=?(x1) 将它与(3.1)式联立，消去未知的L，有 由此推知 课件 在计算了x1及x2之后，可用上式右端作为x*的新近似,记作?x1，一般情形是由xk计算xk+1, xk+2，记 它表明序列{?xk}的收敛速度比{xk}的收敛速度快. (3.1)式称为埃特金(Aitken) △2加速方法. 可以证明 课件 也称为埃特金 ( Aitken ) 外推法. 可以证明: 为线性收敛,则埃特金法为平方收敛； 这个加速迭代法也可写成下面格式 若 为 p ( p 1)阶收敛， 导数连续，则埃特金法为 2p–1 阶收敛. 的 p 阶 若 课件 例题 求方程 x = e –x 在 x=0.5 附近的根. 解 取 x0=0.5, 迭代格式 x25=x26=0.5671433 若对此格式用埃特金法, 则 得 课件 仍取 x0=0.5 , 得 由此可见, 埃特金法加速收敛效果是相当显著的. 课件 7.3.2 斯蒂芬森(Steffensen)迭代法 埃特金方法不管原序列{xk}是怎样产生的，对{xk}进行加速计算，得到序列{?xk}. 如果把埃特金加速技巧与不定点迭代结合，则可得到如下的迭代法： 称为斯蒂芬森(Steffensen)迭代法. 它可以这样理解，我们要求x=?(x)的根x*，令误差ε(x)=?(x)-x，有等式ε(x*)=?(x*)-x*=0，已知x*的近似值xk及yk，其误差分别为 课件 把误差ε(x)“外推到零”，即过(xk,ε(xk))及(yk,ε(yk))两点做线性插值函数，它与x轴交点就是(3.3)中的xk+1，即方程 的解 课件 实际上(3.3)是将不定点迭代法(2.