数值分析教学课件3.3.pptVIP

第三章 数值积分与数值微分 3.3 外推原理与Romberg求积法 3.3.2 Romberg 求积法 3.3.1 外推原理 3.3 外推原理与Romberg求积法 学习目标： 理解外推原理，会运用Romberg求积法。 　　在科学与工程计算中，很多算法与步长h有关，特别是数值积分、数值微分和微分方程数值解的问题。对于这些算法，我们可以通过外推技巧提高计算精度。先看一个计算π的近似值的例子，由函数sinx的Taylor展开式有 若记 则有 3.3.1 外推原理 3.3 外推原理与Romberg求积法 由此构造新的表达式： 可见，计算π的近似值的算法F(h)的截断误差是 ，而算法 的截断误差是 。外推一次，精度提高了。这就是外推法的基本 思想。 若重复以上过程，不断外推，即不断折半步长h，得到计算π的 算法序列 。随着k的增加，算法的截断误差越来越高，计算精 度越来越好。 可将外推思想推广到一般情况。设F(h) 是计算F(0)的一种近似算式， 带截断误差的表示式为 其中， 与p无关。 如果我们用h和h/q（q1）两种步长分别计算F(