Simpson(辛普森)积分公式

造bug的小歪猫

已于 2023-09-25 00:49:39 修改

阅读量1.4w

点赞数 14

分类专栏：路径规划文章标签：算法

于 2023-09-23 10:50:34 首次发布

本文链接：https://blog.csdn.net/qq_41186865/article/details/133198282

版权

路径规划专栏收录该内容

2 篇文章

订阅专栏

文章详细介绍了微积分中的牛顿-莱布尼茨公式，并重点讨论了辛普森积分公式，通过二次曲线逼近提高积分精度。文章还提到如何通过二分法和递归设计Simpson自适应算法来控制计算精度，但未深入数学误差分析。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

一、牛顿 - 莱布尼茨公式（微积分的基本公式）：
$\int_{a}^{b}f\left(x\right)dx=F(b)-F(a)$
其中，函数F(x) 是连续函数f(x)在区间[a,b]上的原函数。
二、辛普森（Simpson）积分公式
1.Simpson积分公式是将区间端点和区间中点三个点近似看成抛物线上对应的三个点，以二次曲线逼近的方式取代矩形或梯形积分公式，以求得定积分的数值近似解。
$\int_{a}^{b}f\left(x\right)dx\approx\ \ \frac{(b-a)}{6}\bullet\left[f\left(a\right)+4f\left(\frac{a+b}{2}\right)+f(b)\right]$
2．Simpson积分公式可以利用微积分基本公式推导出来，推导过程如下：
设置一个一元二次函数
$f(x)=ax^2+bx+c$
求积分为
$F\left(x\right)=\int_{0}^{x}{f(x)}=\ \frac{a}{3}x^3+\frac{b}{2}x^2+cx+d$
那么
$\begin{align} \int_{L}^{R}{f\left(x\right)\ }dx \nonumber &=\ F\left(R\right)-F\left(L\right)\\ \nonumber &=\frac{a}{3}{(R}^3-L^3)+\frac{b}{2}{(R}^2-L^2)+c(R-L)\\ \nonumber &=\frac{a}{3}{(R}^3-\ L^3)+\frac{b}{2}{(R}^2-L^2)+c(R-L)\\ \nonumber &=\ (R-L)\left[\frac{a}{3}{(R}^2+\ L^2+LR)+\frac{b}{2}(R+L)+c\right]\\ \nonumber &=\frac{\left(R-L\right)}{6}(2aR^2+2aL^2+2aRL+3bR+3bL+6c)\\ \nonumber &=\frac{\left(R-L\right)}{6}\left\{\left(aL^2+bL+c\right)+\left(aR^2+bR+c\right)+4\left[a\left(\frac{L+R} {2}\right)^2+b\left(\frac{L+R}{2}\right)+c\right]\right\}\\ \nonumber &=\frac{\left(R-L\right)}{6}[f(L)+4f(\frac{L+R}{2})+f(R)]\\ \nonumber \end{align}$
有了 Simpson积分公式，一个自然的想法是把积分区间拆成多个小区间后求和，但是分成区间的个数和长度因积分区间和精度要求甚至被积函数而异。
本文不在从数学角度进行误差分析，不是本文的重点。而从计算机运算的角度，精度的控制可以通过设定一个能够容忍的误差作为最终的返回情况，基于二分设计出递归的计算流程，即所说的Simpson自适应算法。