数值分析第七章数值微分与数值积分讲述.pptVIP

上机作业 第260至261页 数值实验1题, 2题任选一题 ? 当区间[a, b] 2n等分时, 则有 在定理中以 h/2 代替 h 得 上式乘以4减去 T(h) 再除以3, 记之为 T1(h), 得 T1 (h) 逼近 I 的速度是 O ( h4 )阶, 效果比 T (h)好, 它不是别的, 就是Simpson序列. ? 类似地 上式乘以16减去 T1(h) 再除以15, 记之为 T2(h), 得 T2 (h) 逼近 I 的速度是 O ( h6 )阶, 效果比 T1 (h)好, 它不是别的, 就是Cotes公式序列. ? 对Cotes公式序列进行同样处理得到Romberg公式序列. Richardson外推加速方法 也称为Romberg求积算法 ? 收敛性说明: 如果 f (x) 充分光滑, 那么梯形公式序列, Simpson公式序列, Cotes公式序列, Romberg公式序列均收敛到所求的积分值. ? 对于 f (x)不充分光滑的函数也可用Romberg算法计算, 只是收敛慢一些. 也可以直接使用复化Simpson公式计算. 例 用Bomberg算法计算积分 解 在[0, 1]上仅是一次连续可微 用Romberg算法计算结果见下表 2 4 1 8 16 32 0.5 0.426777 0.407018 0.401812 0.400463 0.4001