看了一些大牛的博客,终于对自适应Simpson积分有了初步的了解。
Q:自适应Simpson积分是用来做什么的?
A:用来求积分啊!
Q:它和一般的微积分有什么区别?
A:不能求出精确解,但是可以近似地求出一些微积分不能够求解的图像的面积。
首先是Simpson积分公式:
(感谢维基百科)
由于直接套这个公式很容易得出差得离谱的解,所以要稍加改进,使用自适应Simpson积分:
1.取三个点a, (a+b)/2, b
2.利用Simpson积分公式分别计算原图像在[a, b], [a, (a+b)/2], [(a+b)/2. b]的面积(分别记为S0, S1, S1),若S0与S1+S2的值相差无几,则可以认为S0为原图像在[a, b]内的面积。
另外需要提及的是,自适应Simpson积分不仅可以求“正常”函数的积分,还可以用来求不规则图像的面积。这个时候,我们就不能狭隘地认为f[a]表示原函数在自变量为a的时候的值,而应该视为直线x=a被图像覆盖的长度(这也就是微积分的本质)。最后求出来的面积是整个图像的总面积,