级牛顿法.pptVIP

* §5 牛顿迭代法 广泛应用的 高效计算方法 基本思想 非线性方程局部线性化 (化繁为简) 特点 单根附近具有较高的收敛速度 5.1 Newton迭代公式 ( Newton法或切线法或牛顿切线法 ) 1 公式 设非线性方程 ,其精确解或真解为 二次连续可导, xk是(4.1)的k次近似解. 介于x与xk之间 用泰勒展开式的线性部分近似 f(x),即有 并设f (x) 则f(x)在点xk的泰勒展开式为 Newton迭代公式 把 作为第k+1次近似解，即得 说明 Newton法是一个逐次线性化方法. 2 几何意义 设曲线y =f(x)的图象如图示, 的k 次近似, 曲线 y =f(x) 上对应的点为 过该点作曲线 y=f (x)的切线Lk 交x 轴于点xk+1 近似曲线y =f(x), 则Lk与x轴的交点 作为f(x)=0的第k+1次近似解. 即由(4.2)式令y =0得 该式就是Newton法的迭代公式, 因此Newton法也称为切线法. 例1 用牛顿切线法求方程x = e-x在0.5附近的根 . 解 由 x=e-x 得 xex=1， 取 f (x)=xex-1， 取x0=0.5计算结果如表5-1所示 . 0.5 0 0.567155569 0.567143291 0 0 1 2 3 4 表5-1 说明 牛顿切线法比一般迭代法收敛要快得多 . 牛顿公式 3 牛顿切线法在数学上的应用实例 简单迭代法的例4 迭代18次 解方程、求算术平 方根、造倒数表等 3. 牛顿切线法在数学上的应用实例 解方程、求算术平 方根、造倒数表等 例2 不用除法运算，构造求 解 的值的迭代公式. 迭代格式 令 取 取 例3 导出计算 解 的牛顿公式，并计算 计算结果如表5-2所示 . 10 10.750000 10.723837 10.723805 10.723805 0 1 2 3 4 表5-2 当k=3时，得8位有效数字10.723805. 计算 是求x2-c = 0的正根， 解出x即得到 取 f (x)=x2-c ， 的值在10和11之间， 取初值x0=10，则由牛顿迭代公式 注 牛顿法也存在收敛问题，并不 是对所有函数牛顿法都收敛.如图示: 因此讨论收敛 定理. 牛顿法迭代会越跑越远. 5.2 Newton法收敛定理 定理9 （Newton法局部超线性收敛性） 如果f(x)在解x*邻近连续 可导,且 则存在 只要 收敛于x*即 Newton序列超线性 分析 Newton迭代公式为 若证其收敛性，可 看作简单迭代，即 则迭代函数为 由f(x)在解x*邻近连续可导， 则g(x)在不动点x*处可微，若 则利用推论2知迭代序列{xk}收敛于x*，而且 成立. 证明 其中 介于 与 之间. 即 若令 则 再由 f (x)的连续性知 则 即 所以由推论2得 ＃ 中值定理 牛顿法二阶收敛 分析 由定理9知Newton法超线性收敛，因此该定理只证(5.6) 式,因(5.6)式中有f(x)的二阶导数及 因此试用含二阶导数 项的泰勒展开式. 定理10 （Newton法局部二阶收敛性） 如果f(x)在解 邻近二次 则存在 只要 Newton序列二阶 即 连续可导,且 收敛于 证明 首先 将 在点 进行泰勒展开: 代入上式,得 其中 介于 与 之间. 说明 当 时,有 Newton法超二阶收敛. 定理9和定理10都要求满足 对初值要求高,初值 选取较困难. 在实际应用中可以放宽对初值的要求,即有以下的非 局部收敛定理. # 当 时, 并由 的连续性,则 其中 介于 与 之间. 注 （1）掌握Newton迭代公式并会用Newton迭代法编程求方程的近似解. （2）掌握Newton法局部二阶收敛定理的证明. 作业 2. 证明计算 的Newton程序为 (并用此迭代 编程计算 精确到10-6 ) . 3.导出计算 的牛顿迭代程序，要求公式中既无开方， 又无除法运算. 上机作业 本节例题1、2及作业1、2 1. 用Newton法迭代两次求方程x3-3x-1=0在x0=2附近的根，取4 位有效数字（准确解x*=1）. *