计算方法 第2章 一元线性方程的解法概要.pptVIP

上述迭代法需要计算导数a≈g′(x)， 下面介绍埃特金(Aitken)迭代方法。它也是一种加速迭代法。 埃特金迭代法 将上两式联立消去a得到 可解出 这样得到埃特金迭代公式 例6 用埃特金迭代法求 ?x3-x-1=0 ?在(1,1.5)内的根。 解 前面已经提到,迭代公式 ?xk+1=x3k-1, k=0,1,2,… ?是发散的。 现用埃特金算法来求根,其迭代公式为 仍取x0=1.5,计算结果见下表。 定义2.2.1 设序列 收敛到 ，若有实数 和非零常数C，使得 其中 ， ，则称该序列是p 阶收敛的，C 称为渐进误差常数。 迭代法收敛的阶 当p=1时，称为线性收敛； 当p1时，称为超线性收敛； 当p=2时，称为平方收敛或二次收敛。 一般迭代法是线性收敛 线性收敛到 。 迭代法的收敛速度 } { ) ( 0 ) ( 0 ) ( ... ) ( ) ( ) ( , ) ( * 0 1 0 * ) ( * ) 1 ( * * * x p x x x x x x x x x x x x n k k p p j j j j j j j = D ? 1 = = = ￠ ￠ = ￠ = ￥ + - 阶收敛速度收敛到 ，以 列 产生的序 ，由 ，则 而 邻域内连续，且有 的不动点 是 设 定理 ) ( p 在x*的 证： ! ) ( ) ( lim * ) ( * * 1 p x x x x x p p n n n j = - - + ￥ ? 将 (x ) 在x*附近做泰勒展式： j 因为有： ，由上式有： 定理：由埃特金迭代公式 产生的序列 至少平方收敛到 。 n=0，1，2，… n n+1 n+1 n n+1 n n n+1 n+1 n n+1 x x x x x x x x x x x + - - - = = = + 2 ) ( ) ( ) ( 2 1 j j 还有： 由埃特金迭代公式至少是2阶（平方）收敛的。 由线性收敛知 当n充分大时有 即 加速迭代法的另种推导 展开有： 已知 ，则 ， 改成 n=0，1，2，… n n+1 n+2 n n+1 n n n+1 n+2 n n+1 x x x x x x x x x x x + - - - = = = + 2 ) ( ) ( ) ( 2 1 j j 也可以改写成 其中迭代函数 习题 习题二:2,3,4,5 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * 若从任何可取的初值出发都能保证收敛，则称它为大范围收敛。如若为了保证收敛性必须选取初值充分接近于所要求的根，则称它为局部收敛。 通常局部收敛方法比大范围收敛方法收敛得快。因此，一个合理的算法是先用一种大范围收敛方法求得接近于根的近似值（如对分法），再以其作为新的初值使用局部收敛法（如迭代法）。 这里讨论迭代法的收敛性时，均指的是局部收敛性。 则方程在(a,b)内有唯一的根； 定理/收敛定理/压缩映像原理 设方程x=g(x)在(a,b)内有根x* g(x)满足李普希茨(Lipschitz)条件:即对(a,b)内任意的x1 和x2都有： q为某个确定的正数，q＜1 条件 结果 3. 还有误差估计式 2. 且迭代公式 xk+1=g(xk) 对任意初始近似值x0均收敛于方程的根x； 由已知条件知，x*为方程x=g(x)的根，即x*=g(x*) 设 也是方程的根，即 于是，由李普希茨条件得 因为q＜1，所以上式矛盾，故必有 证 1. 有唯一根 考虑迭代公式 x k+1=g(xk) , k=0,1,2,… 由李普希茨条件 证 2. 迭代公式收敛 因为q＜1，当k→∞