Romberg算法幻灯片.pptVIP

第三节 Romberg算法 * ? 2009, Henan Polytechnic University * §3 Romberg算法 第四章 数值积分与数值微分 * 4.3.1 步长的选取: 通过上述截断误差可以看出用复化的求积公式提高了计算精度，但是对给定的误差限，应将[a,b]分成多少份？由截断误差可以大概算出步长的取值情况，但是比较麻烦，接下来我们考虑步长的选取。 x 0 0.125 0.25 0.375 0.5 0.625 0.75 0.875 1 f(x) 1 0.9973978 0.9896158 0.9767267 0.958851 0.9361556 0.9088516 0.8771925 0.8414709 例：用变步长的梯形法计算 解： 注：若精度达不到，则应继续二分下去。 类似于梯形公式，也可得到变步长的Simpson公式及Cotes公式。 4.3.2 Romberg 求积公式 根据复化梯形法的误差公式，我们推得 希望得到的结果更精确。 类似的，可得 用(1) (2) (3)来计算积分,得到的结果应更精确。 精度如此高的积分公式其实就是复化的Simpson公式。 上述公式称为Romberg公式。 综上所述，用对分区间的方法，先算出T2n，再由(1)(2)(3)很容易算出S2n,C2n及R2n，通过这种方法，能将粗糙的积分值逐步加工成精度较高的积分值 流程图如下：