第七章数值微分与数值积分初次修改稿1.pptVIP

上机作业 第260至261页 数值实验1题，2题任选一题 2 4 1 8 0.9207355 0.9397933 0.9445135 0.9456909 0.9461459 0.9560869 0.9460833 0.9460830 0.9460831 0.9460831 ? 利用只有两三位有效数字的T1, … ,T8 经过三次外推得到7位有效数字. 可见加速的效果十分显著. 用Romberg算法计算如下 ? 理论依据: 复化梯形公式的余项展开. 记 ?定理 设 则 其中系数 ?k ( k=0, 1, …)是与 h 无关的常数. ? T (h) 逼近 I 的速度是 O ( h2 )阶. ? 当区间[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外推加速方法 也称为Romb