辛普森(Simpson)公式是牛顿-科特斯公式当n=2时的情形,也称为三点公式。利用区间
二等分的三个点来进行积分插值。其科特斯系数分别为1/6,4/6,1/6。
应用
公式内容
设
拟柱体的高(两底面α,β间的距离)为H,如果用平行于底面的平面γ去截该图形,所得到的
截面面积是平面γ与平面α之间距离h的不超过3次的函数,那么该拟柱体的体积V为
式中,和是两底面的面积,是
中截面的面积(即平面γ与平面α之间距离h=H/2时得到的截面的面积)。
事实上,不光是拟柱体,其他符合条件(所有顶点都在两个平行平面上、用平行于底面的平面去截该图形时所得到的截面面积是该平面与一底之间距离的不超过3次的函数)的
立体图形也可以利用该公式
求体积。
计算实例
解:此题中S_1 = S_0 = S_2 = S,H = h,所以V = H (S_1 + 4S_0 + S_2) /6 = h (S + 4S + S) /6 = S h。
解:此题中S_1 = S,S_0 = S /4,S_2 = 0,H = h,所以V = H (S_1 + 4S_0 + S_2) /6 = h (S + 4S /4 + 0) /6 = S h /3。
例3:计算半径为r的球的体积。
解:此题中S_1 = S_2 = 0,S_0 = πr^2,H = 2r,所以V = H (S_1 + 4S_0 + S_2) /6 = 2r (0 + 4πr^2 + 0) /6 = 4πr^3 /3。
公式证明
只需要证明根据公式算出来的体积和用积分算出来的体积相等即可。
设
截面面积是截面高h的不超过3次的函数:f(h)= ah^3 + bh^2 + ch + d。
那么,
利用积分计算体积,可以得到(积分限为0~h):
V = ∫ f(x) dx
= ah^4 /4 + bh^3 /3 + ch^2 /2 +dh;
利用公式计算体积,可以得到:
V = H (S_1 + 4S_0 + S_2) /6
= h ( f(0) + 4f(h/2) + f(h) ) /6
= h [ d + 4 (ah^3 /8 + bh^2 /4 + ch /2 + d) + (ah^3 + bh^2 + ch + d) ]/6
= ah^4 /4 + bh^3 /3 + ch^2 /2 +dh。
因此两式相等,公式得证。
Remark:当函数f(h)次数高于或等于4次时,公式一般不成立。这只需验证f(h)=h^4时公式不成立即可。
C++代码实现