Simpson Rule Summary 辛普森法则.docVIP

考虑积分，如果区间[a，b]内取等间隔的N份，间隔长度为h，简述Simpson法则缓变函数区间中点上值近似平均值 其中f(x)在区间[a,b]上的积分 第i个区间对积分的贡献为： 如果围绕该区间中点的邻域内对函数f(x)作泰勒级数展开 其中，和分别表示f(x)在处的一阶，二阶和三阶导数。积分在子区间内的值可以表示为 其中第一项是近似，第二项则由于其中的积分等于零而消除。从而，矩形法则在宽度为h的单个子区间内的最高阶误差由第三项给出 在整个[a,b]区间上的总误差则通过将所有N个子区间的贡献相加得到 其中我们利用了Nh=(b-a)，并且取为f(x)在[a,b]上的二阶导数的均值。 梯形法则用一个梯形来近似代替每个子区间的面积，梯形的四个顶点分别位于，，与。梯形的面积为 （2-8） 对于整个[a,b]区间，积分值由所有窄带的和给出 （2-9） 与其他值相比，与的贡献只具有一半的，这是由于它们位于区间的两个端点。Euler-Maclaurin积分求和公式， （2-10） 系数为伯努利数，，且当时，。有 ，而其他所有奇数阶项为零。前几个偶数阶项的值如下 更一般的，它们来生成 梯形法则的优点在f(x)的值只在格点（网点）上，即处求出。，此方法中由有限步长引起的误差的预期值是矩形法则中的2倍。而另一方面，梯形法则的简洁性则有助于进行进一步的改进。此方法可以与外插法一起使用而构成一种非常高效的数值积分算法。框图2-1列出了梯形法则算法的主要步骤。该算法还包含了一种方法，即由一个小的子区间数开始，通过循环迭代逐步增加数目，而无需在格点上重复计算被积函数的值。Box 2-1 Subroutine TRAPZ(A,B,KEY,RSLT,H) Trapezoidal rule integration Arguments: A: Lower integration limit. B: Upper integration limit. KEY=0: Initiate a new integration; 0: subsequent iterations. RSLT: Value of the integral. H: Integration step size. Initialization: Define the integrand as an external function. First call (KEY=0): Calculate the integral using only two subintervals. Define the number of subintervals to be 2. Let SUM be the average of the two end points plus the midpoint. Define the step size H to be half of the distance from A to B. Change the value of KEY to 1. Return the product of H and SUM. Subsequent calls (KEY0): Calculate the contribution from points halfway between old ones: Start from the middle of the first subinterval. Add the value of the integrand at this point to SUM. Move to the next subinterval and repeat steps (i) and (ii). Double the number of subintervals and halve the step size H. Return H times SUM. Simpson法则 当积分的上下限相对于展开的中心点对称积分泰勒展开式中含有f(x)的奇数阶导数的项都将等于零。利用这一性质，我们可以在相邻的两个子区间内对面积作泰勒级数展开可得 积分的精确度达到。中心差分，可将f(x)的二阶导数近似表为 将此结果代入式，则子区间[xi-1,xi+1]内的积分近似为 对于整个[a,b]区间，结果为 假设区间的数目N为偶数奇数格点的贡献是偶数格点贡献的2倍。这种权重上的差异来自我们为修正基本方法所得的一级结果所引入的f(x)二阶导数的贡献。两端点和的权重则仅为偶数点的一半。Simpson法则的另一优点在于其自然引出了一