【精选】Romberg算法.pptVIP

Romberg 算法 梯形法递推公式 梯形法递推公式 梯形法递推公式 举例 龙贝格 (Romberg) 加速 Richardson 外推加速法 梯形法的加速 Romberg 算法 举例 举例 作业 * * 4.1 Romberg 算法 数值分析及计算软件 Chap 4 数值积分与数值微分 太 大: 利用复化梯形公式、复化Simpson公式、复化Cotes公式等计算定积分时，如何选取步长 h ？ 计算精度难以保证 太 小: 计算量增大 解决办法：变步长算法 通常采取将区间不断对分的方法，即取 n = 2k ，反复采用复化求积公式，直到所得到的计算结果满足指定的精度为止。 步长折半：[xi , xi+1/2] ， [xi +1/2 , xi+1] 将 [a, b] 分成 n 等分 [xi , xi+1] ， n=1 n=2 n=4 记 解： 例：用梯形法的递推公式计算定积分 , 要求计算精度满足 Cha4.2.m 程序演示 Cha41.m: 复化N-C公式 Cha42.m: 梯形法的递推 梯形法递推公式算法简单，编程方便 解决方法：龙贝格 (Romberg) 加速 但收敛速度较 慢 复化积分公式的渐近状态 梯形公式： 辛普森公式： Cotes公式： 思想：利用余项公式与积分的定义 定义 如果一种复化求积公式 满