§5.3 牛顿迭代法.ppt

5.3 牛顿迭代法 * * * * * * * §5.3 牛顿迭代法 §5.3 牛顿迭代法 过曲线 作为初始近似值， 上的点 作切线 设方程 (图形如下)。首先在根 的附近取一点 ? 切线方程 以上迭代方法－牛顿法（切线法） 若 切线与 轴的交点为 ? 作为根 的第二次近似值。 上的点 作切线，若 以切线与 轴交点 以 作为根 的第一次近似值，然后又过曲线 ? 仿此不断作下去，得到一般迭代公式： ? §5.3 牛顿迭代法 §5.3 牛顿迭代法 由于牛顿法仍是迭代法，其迭代函数为 牛顿法收敛性分析 为常数, 则牛顿法收敛。 若 在根 的某个邻域内不为零，且 存在， 问题：如何保证牛顿迭代法收敛？ 由上节定理可知，若 使在 上， 上保持严格单调和凹向不变。 都不变号，即 在 用逐步搜索法，找到方程的有根区间 (尽量小） 有以下四种情况： §5.3 牛顿迭代法 可以看出，用牛顿法求得的序列 均是单调地趋于 故牛顿法是收敛的。 均可作为初始值。 凡满足关系式 ，图(2)，(3)取 例如图(1)，(4)取 结 论 §5.3 牛顿迭代法 定理 设函数 在 上存在二阶导数，且满足下列条件： (1) (2) 在 上不为零 (3) 在 上不变号 (4) 对 ，有 则牛顿法迭代序列 收敛于方程 在 内的 唯一根 §5.3 牛顿迭代法 1．选初值 （为保证迭代收敛，可用定理或上式选定初值） 4．以 替代 回到2继续，直到 （或 ）为止。 牛顿法编程步骤 否则到4。 3．若 （或 ），迭代终止， 即为根， 计算 。 2．迭代：按公式 ，计算 。 §5.3 牛顿迭代法 例11：用牛顿迭代法编程求方程 在 附近的一个实根。 解： real*8 x0, x1, f, f1 open(1, file=nt.dat) read (*,*) x0 N=1 10 f=x0**3-2.*x0**2+4.*x0+1 f1=3.*x0**2-4.*x0+4. x1=x0-f/f1 write(*,*) N, x0, x1 write(1,*) N, x0, x1 If(abs(x1-x0). gt. 1e-6) then x0=x1 N=N+1 goto 10 end if end 输出结果： 1 1.000000000000000 -3.333333333333333E-001 2 -3.333333333333333E-001 -2.287581699346405E-001 3 -2.287581699346405E-001 -2.225152446820112E-001 4 -2.225152446820112E-001 -2.224945143483070E-001 5 -2.224945143483070E-001 -2.224945141207864E-001 输入：1.0 （经验证满足初值条件） (取初值 也可以，满足条件) §5.3 牛顿迭代法 §5.3 牛顿迭代法 例12：用牛顿迭代法求方程 在区间 内根的近似值。 解：设 又 为初值。 牛顿法比一般的函数迭代法(8次）收敛速度要快。 §5.3 牛顿迭代法 例13：计算 在 附近的实根。 解： ∴ 可取 为初值。 算得： 精确到小数点后四位。 §5.3 牛顿迭代法 弦截法 弦截法(割线法)是它的一种改进，在牛顿迭代公式中， 代入牛顿迭代公式可得 牛顿法计算时要用到函数的导数，很多情况下难以使用。 说明 迭代格式中没有用到函数的导数，计算方便，但收敛速度较牛顿法要慢，开始时要用到两个不同的根的近似值作初值。 特点 用差商代替导数，即 §5.3 牛顿迭代法 例18：在相距100m的两个塔（高度相等的点）上悬挂 一根电缆(如下图所示)，允许电缆中间下垂10m。 试确定悬链线方程 中的参数 解：由于曲线最低点和最高点相差10m，有 ，先构造函数 要确定参数 弦截法计算结果如下表： 上表中 ，所求悬链线方程为 根的两个初始值为 。 由图可知，根的区间在[120，150]范围内，故选取以