《计算方法》第二章.pptVIP

- n xn 0 0 1 0.5 2 03 04 05 0提纲 §2.1 对半分法 §2.2 一般迭代法 §2.3 加速迭代收敛的埃特金(Aitken)方法 §2.4 牛顿法 §2.5 弦截法 §2.5 弦截法 用牛顿法解方程f(x)=0，虽然在单根附近具有较快的收敛速度，但它有个明显的缺点，就是需要计算导数f’(x)，当f(x)比较复杂时，计算f’(x)可能有困难。 此式为解方程f(x)=0的弦截迭代格式。 弦截法与牛顿迭代法的比较： 与牛顿法一样，当在根的领域内有直至二阶的连续导数且 时，弦截法具有局部收敛性，弦截法的收敛速度比牛顿法稍慢，但也是超线性收敛，其收敛阶为1.618。 n xn 0 0 1 1 2 0.612699836 3 0.572181411 4 0.567703214 5 0.567205552 n xn 0 0 1 1 2 0.612699836 3 0.563838389 4 5 第二章 完 定理3 对于迭代过程xk+1=φ(xk)，如果迭代函数φ(x)在所求根x*的邻近有连续二阶导数，且|φ’(x*)| 1，则有 (1)当φ’(x*)≠0时，迭代过程为线性收敛； (2)当φ’(x*)=0时，而φ’’(x*)≠0时迭代过程为平方收敛。 证明： (1)对于在根x*邻近收敛的迭代公式xk+1= φ(xk)， 由于xk+1–x*=φ(xk)-φ(x*)=φ’(ξ) (xk –x*) ξ为x*与xk之间的某一点，因而当xk在根x*的附近时，将有xk+1–x* ≈φ’(x*) (xk–x*) 即 这样，若φ’(x*)≠0，则迭代过程仅为线性收敛。 (2)若φ’(x*)=0，将φ(xk)在根x*处作泰勒展开有（取n=1） 而 ，由上式得 因此，对迭代误差，当k→∞时，有 这表明迭代过程为平方收敛。 收敛阶是衡量算法的重要标准。对于一个迭代过程，只要迭代次数足够多，从理论上讲就能得到满足任意精度的结果。 但这里有一个重要问题——收敛速度问题。若迭代过程的收敛速度太慢，则计算工作量就很大，因而有必要研究加速收敛的方法。 提纲 §2.1 对半分法 §2.2 一般迭代法 §2.3 加速收敛的埃特金(Aitken)方法 §2.4 牛顿法 §2.5 弦截法 从定理3，当迭代函数φ(x)在根x*处有φ’(x*)≠0 时，迭代过程仅为线性收敛，如何能使收敛加速？ §2.3 加速收敛的埃特金(Aitken)法 埃特金算法的几何解释： n 0 1 2 3 4 5 6 7 8 x(1) 1.5 1.35721 1.33086 1.32588 1.32494 1.32476 1.32473 1.32472 1.32472 x(2) 1.5 2.375 12.39 埃x(1) 1.5 1.325 1.32472 埃x(2) 1.5 1.41629 1.35565 1.32895 1.32480 1.32472 提纲 §2.1 对半分法 §2.2 一般迭代法 §2.3 加速迭代收敛的埃特金(Aitken)方法 §2.4 牛顿法 §2.5 弦截法 求f(x)=0的解。 设xk是第k次近似解，将f(x)在xk处作泰勒展开 取前两项来近似代替f(x)，则得近似方程 设f’(xk)≠0 ，令其解为xk+1得牛顿迭代格式 §2.4 牛顿法（切线法） 牛顿迭代法的几何意义 条件1保证了根的存在; 条件2表明了函数的单调，根唯一; 条件3表示f(x)的图形在[a,b]上凹向不变; 条件4保证了x?[a,b]时，?(x)?[a,b]。 如(a)、(d)，若f’’(x)0，f(x)0，即(a)时，x0取b开始，(d)时，x0取a 开始，则f(x)f’’(x)0，，这时求出x1的就不一定在[a,b]之中，反之情况亦然。 下面四种情况都满足定理条件: 备注： 一阶导数f’(x)表示曲线的斜率，若f’(x) ≥0，则函数单调递增，反之递减。 二阶导数f’’(x)表示曲线的凹凸，当f’’(x) 0，曲线上凹； f’’(x)0，曲线上凸。 第二章 一元非线性代数方程的求解 求方程 f(x)=0的实根的几何意义是寻找曲线 y= f(x)与x轴的交点。但一般来讲，不可能找到一个数学式子来表