Simpson-Rule-Summary---辛普森法则分析和总结.docxVIP

考虑积分 I  [ a,b] ? ? b a f (x)dx ，如果在区间[a，b]内取等间隔的N 份，间隔长度为h，简述矩 形（、梯形、Simpson 法则）计算积分的i）理论、误差精度分析，和算法计算流程。解：对于缓变函数我们可以用各个区间中点上函数值作来近似该区间的平均值 其中 x i?1/ 2  ? 1 (x 2 i?1  x ) 。 i  f ? f (x ) i i?1/ 2 矩形法则：f(x)在区间[a,b]上的积分用矩形求积定义如下 第 i 个区间对积分的贡献为：  I [a,b] ? h?N i?1  f i?1/ 2 I x[ x i?1 ? ? ,x ] i x i xi?1 x f (x)dx ? hf  i?1/ 2 如果围绕该区间中点 x i?1/ 2 的邻域内对函数f(x)作泰勒级数展开， 有 1 1 1 f (x) ? f ? f ? (x ? x ) ? f ? (x ? x )2 ? f (3) (x ? x )3 ?L i?1/2 1! i?1/2 i?1/2 2! i?1/2 i?1/2 3! i?1/2 i?1/2 其 中 f ? i?1/ 2 ， f ? i?1/ 2 和 f (3) i?1/ 2 分别表示了 f(x)在 x ? x i ?1/ 2 处的一阶，二阶和三阶导数。相 应地，积分在子区间内的值可以表示为 ?ix f (x)dx ? f ? i x dx ? 1 f ? x (x ? x )dx xi?1 x i?1/ 2  x?i?i?1 x ? i ? 1! i?1/ 2  x?ii?1 x ? i i?1/ 2 ? 1 f ? x (x ? x )2 dx i2! i?1/ 2 i  x?i?1 x ? i?1/ 2 1? f (3) 1 x (x ? x )3 dx i3! i?1/ 2 i i?1/ 2 xi?1 x ?L 其中第一项是矩形积分的近似值，第二项则由于其中的积分等于零而消除。从而，矩形法则 在宽度为h 的单个子区间内的最高阶误差由第三项给出 i?I ? ?x i f (x)dx ? hf x[ x 3i?1 3 , x ] i  xi?1 x ? 1 f ? i?1/ 2 ?ix (x ? x ? i )2 dx ? h f ? 2! i?1/ 2  xi?1 x i?1/ 2 24 i?1/ 2 在整个[a,b]区间上的总误差则通过将所有N 个子区间的贡献相加得到 b b ? a 1 (b ? a)3 ?I [a,b] ? ? f (x)dx ? I a  [a,b] ? h2 f ?(?) ? f ?(?) 24 24 N 2 其中我们利用了 Nh=(b-a)，并且取 f ?(? ) 为 f(x)在[a,b]上的二阶导数的均值。梯形法则： 如果用一个梯形来近似代替每个子区间的面积，梯形的四个顶点分别位于 (x i?1 ,0) ， (x i?1 , f i?1 ) ， (x , f i i ) 与(x ,0) 。梯形的面积为 i I x[ x i?1  , x ] i ? h ( f 2  i?1 ? f ) （2-8） i 对于整个[a,b]区间，积分值由所有窄带的和给出 I [a,b] ? 1 h?N 2 i?1 ( f i ?1 ? f ) ? h( 1 f i 2 0 ? f ? f 1 2 ?L ? f ? 1 N ?1 2 f ) （2-9） N 与其他值相比， f 与 f 的贡献只具有一半的权重，这是由于它们位于区间的两个端点。根 0 N 据 Euler-Maclaurin 积分求和公式，我们有 1? b f (x)dx ? h{ f ? f ? f ?L ? f ? 1 f } 1 a 2 0 1 2 N ?1 2 N （ ） B2? 2 h2 ( f ? ? f ?) ?L ? B k B 2 h2k ( f (2 k ?1) ? f (2 k ?1) ) ?L 2-10 2! N 0 (2k )! N 0 系数 B 为伯努利数，其满足以下下规律： B (x) ? 1 ， B (x) ? nB (x) 且当 n ? 1 时， 2k 0 n n?1 ? 1 B (x)dx ? 1 。 有 B ? ? ，而其他所有奇数阶项为零。前几个偶数阶项的值如下 10 n 1 2 1 1 1 1 B ? 1 B ? B ? ? B ? 0 更一般的，它们可以通过 t et ?1 2 ? ?? B k k ?0 6 4 30 t k 来生成 k ! 6 42 梯形法则的