计算方法7.4-7.6剖解.ppt

* * 双点弦截迭代法较单点弦截迭代法收敛的速度快，但一般情况下比牛顿迭代法慢，是超线性收敛的，收敛阶数是 双点弦截法是牛顿迭代法的变形，实际上可以看作是用差商代替函数f(x)在xk处的导数，其最大的优点就是避免了求导数。都是线性化方程，但两者有本质区别。Newton切线法在计算xk＋1时只用到前一步的xk及f(xk)，但要计算f ′(xk)，而双点弦截法在计算xk＋1时要用前面两步的结果xk，xk－1， f(xk)， f(xk－1)，而不须计算导数。这种方法必须有两个启动值x0， x1 。 0 0.5 1 0.6 2 0.56754 0.03246 3 0.56715 –0.00036 4 0.56714 –0.00001 * 例12 用双点弦截法求方程 在 附近的根（ ） 解：取 由迭代公式求得下表 故 ，满足精度要求。 * 令 * * 本节作业 P168 （9） 用迭代法可逐步精确方程 根的近似值，但必须要找到 的等价方程 ,如果 选得不合适,不仅影响收敛速度,而且有可能造成迭代格式发散。能否找到一种迭代方法,既结构简单,收敛速度快,又不存在发散的问题。这就是本节要介绍的牛顿迭代法。 7.4.1 牛顿迭代法的基本思想 牛顿迭代法一种重要和常用的迭代法, 它的基本思想是借助泰勒展开，将非线性函数f(x)逐步线性化, 从而将非线性方程 f(x)=0 近似地转化为线性方程求解 7.4 牛顿迭代法 * * 对于方程 ,设其近似根为 , 函数f(x)可在 附近作泰勒展开 忽略高次项,用其线性部分作为函数f(x)的近似， 设 的根 ,则有 ,即 * * x1是比x0更接近于根x*的近似值。 同样将f(x)=0在x1附近作泰勒展开，取近似得到 依此类推，得到牛顿迭代公式 因此得到方程f(x)=0的近似根数列{xi} 。 * * 注：当f′(xk )=0时，由于此时在xk处切线是水平的，因此无法计算出xk+1 Newton迭代的等价方程为： 所以 若f(x)在a处为单根，则 所以，迭代格式在根a附近收敛 * * 7.4.2 牛顿迭代法的几何解释 方程f(x)=0的根x*是曲线y=f(x)与x轴交点的横坐标,设xk是根x*的某个近似值,过曲线y=f(x)的横坐标为xk的点Pk=(xk ，f (xk))引切线交x轴于xk+1 , 并将其作为x*的近似值 重复上述过程,每次通过求解切线方程来求解方程f(x)=0的近似根,所以也称为切线法。 * * * * 例7. 用牛顿迭代法求方程的根: 解: 由牛顿迭代法 x0 =0.5 x1 =0.3333333333 x2 =0.3472222222 x3 =0.3472963532 x4 =0.3472963553 迭代四次 精度达10-8 设函数 ，且满足 定理2（牛顿法收敛的充分条件） 1) f (a) f (b) 0； 2) 在整个[a, b]上 f ? 连续不变号且 f ? (x) ? 0； 3) 选取 x0 ? [a, b] 使得 f (x0) f ? (x0) 0； 则牛顿迭代法产生的序列{ xk } 收敛到f (x) 在 [a, b] 的唯一根。 有根 只有单根，根唯一 产生的序列单调有界，保证收敛。 * * 7.4.3 牛顿迭代法的收敛性 y x 0 B=x0 f′′(x)0 a y x 0 B f′′(x)0 a=x0 y x 0 B=x0 f′′(x)0 a y x 0 B f′′(x)0 a =x0 * * 例8 用迭代法求 在隔根区间[1.4,1.5] 内的根，要求准确到小数点后第4位。 （1）牛顿迭代公式为 （2） 当 时有， 因 ,故取 ,牛顿迭代法收敛。 * * 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? * * x0 ? 不满足迭代条件时，可能导致迭代值远离根的情况而找不到根或死循环的情况 * * 牛顿迭代法 x0 = 2 x1 = -3.54 x2 = 13.95 x3 = -279.34 x4 = 122017 如 x0 =1 x1 = -0.5708 x2 = 0.1169 x3 = -0.0011 x4 = 7.9631e-010 x5 = 0 收敛 发散 7.4.4 牛顿迭代法求开方（自己看