[工学]数值计算方法第二章讲义ppt.ppt

注1：使用牛顿迭代法存在从一个根跳到另一个根的情况。 注2：如果f(x)=0没有实根，则牛顿迭代序列不收敛。 全局收敛性 定理2.2.3 保证了根的存在性 函数单调，根惟一 函数图形的凸向不变 几何解释 例2.3.2 解 练习： 用Newton法解方程 在区间[0,1]上的根，要求精度为0.00005。 解 由定理2.2.3， Newton法确定的迭代序列 收敛到区间内的唯一根。取 在Matlab执行 fname=inline(‘x*exp(x)-1’); dfname=inline(‘(x+1)*exp(x)’); [x_star,k]=Newton1(fname,dfname,0.5) 迭代4次后得计算结果 k 0 1 2 3 4 xk 0.5 0 0 0 0说明：在Newton法中初始值的选取较困难： 要保证(1)成立，由一阶Taylor展开式，可以得出 定理 2.7 补充 例2.3.3 解 根据全局收敛性定理，对任何满足 上述Newton迭代产生的迭代序列 均收敛于 练习 解 在Matlab窗口执行 fname=inline(‘x*x-115’);dfname=inline(‘2*x’); [x_star,k]=Newton1(fname,dfname,10) 迭代4次后得到计算结果 k 0 1 2 3 4 xk 10 1010 1010简化Newton法 此格式称为简化Newton迭代。 迭代函数为 迭代格式 几何意义 此格式称为推广的简化Newton迭代。 迭代格式为 进一步简化 由简单迭代的局部收敛条件 收敛性与收敛速度 得 推广的简化Newton法收敛时，因为 所以，推广的简化Newton法只有线性收敛速度。 Newton下山法 在Newton迭代法中，若函数较复杂，初值的选取较困难时，为防止迭代发散，可改用如下的迭代式，以扩大初值的选取范围： 成立， 上述迭代方法称为Newton下山法。 终止准则： Newton下山法的计算步骤 要控制迭代次数？ 例2.13 解 弦割法 迭代格式 由上式定义的迭代算法称为弦割法，也称为割线法。 为什么称为弦割法？是从它的几何意义而言。 几何意义 切线 割线 还是以直线代替 曲线的近似方法啊! 收敛性与收敛速度 例2.3.6 解 两步方法 弦割法的计算步骤 例2.3.7 解 取x0=1,x1=2,代入公式，计算结果如下表所示。 k xk f(xk) 0 1 -1 1 2 5 2 1.166666667 -03 1.253112023 -04 1.337206444 0.053880579 5 1.323850096 -0.0036981168 6 1.324707936 -4.273521*10E-5 7 1.324717965 3.79*10E-8 上式称为单点弦割法。相应地，此前的迭代格式也称为双点弦割法。 单点弦割法的迭代函数： 单点弦割法的几何意义： 单点弦割法 由简单迭代的一般收敛性理论可以推出：单点弦割法是局部收敛的且一般只有线性收敛速度。 与Newton迭代法类似，可在一定条件下得到单点弦割法的全局收敛性。 比 双点 弦割法 收敛慢 单点弦割法的全局收敛性 定理2.3.2 例2.3.8 解 从而，有如下误差估计： 定理 在推论的条件下，有误差估计式 例2.2.2 解 也可化为等价方程 . 但此时定理条件不成立，迭代序列不能保证收敛。 练习： 解 注4. 由于定理中条件(1)一般难于验证，而且在大区间上，这些条件也不一定都成立。所以实际使用迭代法总是在根的邻近进行。 表明收敛性与初值的选择有关！ 定理2.2.2 上面定理所涉及的收敛性，即在根邻域的收敛性称为局部收敛性，具有局部收敛性质的迭代法通常对初值的要求很高，使用起来不太方便。因而，人们通常希望迭代算法对相对大的范围的初始点具有收敛性，这种收敛性称为全局收敛性。 构造迭代法，人们自然希望迭代算法不仅收敛，而且收敛快，