773 字
4 分钟
数值积分:当解析解「摆烂」,计算机如何算出定积分?

一张温暖的俯拍特写——旧木桌上摊开着一本数值分析教材,书页上手写着积分公式和梯形法则示意图,旁边搁着一支老式钢笔,清晨的阳光透过窗户洒在纸面上,投下柔和的阴影,背景书架被虚化成温暖的焦外光斑

宝贝们好呀~今天我们来聊一个非常实用的话题:数值积分 🧮✨

大家应该都有这种体验——翻开微积分课本,里面的定积分总是算得漂漂亮亮:01x2dx=13\int_0^1 x^2 dx = \frac{1}{3}0πsinxdx=2\int_0^\pi \sin x \, dx = 2……但一到现实世界就发现,大部分函数的原函数根本写不出来!比如 01ex2dx\int_0^1 e^{-x^2} dx,不好意思,初等函数表示不了。这时候就需要数值积分登场啦~

核心思想:用简单图形替代曲线#

定积分的几何意义是曲线下的面积。如果我们没法精确算出这个面积,能不能用一些简单的几何图形来近似它呢?

梯形法则(Trapezoidal Rule)#

最直观的想法:把区间 [a,b][a, b] 分成 nn 个小区间,每个小区间 [xi,xi+1][x_i, x_{i+1}] 上用梯形的面积代替曲线下的面积:

abf(x)dxh2[f(x0)+2f(x1)+2f(x2)++2f(xn1)+f(xn)]\int_a^b f(x) dx \approx \frac{h}{2} \left[ f(x_0) + 2f(x_1) + 2f(x_2) + \cdots + 2f(x_{n-1}) + f(x_n) \right]

其中 h=banh = \frac{b-a}{n} 是步长。梯形法则的误差是 O(h2)O(h^2)——把步长减半,误差大约变为原来的 1/41/4

辛普森法则(Simpson’s Rule)#

如果我们不满足于用直线(梯形)来近似曲线,而是用抛物线呢?这就是辛普森法则的核心思想——每三个点过一条二次曲线:

abf(x)dxh3[f(x0)+4f(x1)+2f(x2)+4f(x3)++f(xn)]\int_a^b f(x) dx \approx \frac{h}{3} \left[ f(x_0) + 4f(x_1) + 2f(x_2) + 4f(x_3) + \cdots + f(x_n) \right]

注意系数的规律是 1, 4, 2, 4, 2, …, 4, 1——非常优雅!辛普森法则的误差是 O(h4)O(h^4),远优于梯形法则。

高斯求积(Gaussian Quadrature)#

梯形和辛普森都是等距节点的。但等等——谁说节点一定要等距?如果我们允许自由选择节点位置,就可以用 nn 个节点达到 2n12n-1 次多项式的精确积分!

高斯求积的基本形式是:

11f(x)dxi=1nwif(xi)\int_{-1}^1 f(x) dx \approx \sum_{i=1}^n w_i f(x_i)

其中 xix_i 是勒让德多项式的根,wiw_i 是计算出的权重。例如两点高斯求积:

11f(x)dxf(13)+f(13)\int_{-1}^1 f(x) dx \approx f\left(-\frac{1}{\sqrt{3}}\right) + f\left(\frac{1}{\sqrt{3}}\right)

什么时候用哪个?#

方法精度适用场景
梯形法则O(h2)O(h^2)快速粗糙估计,周期性函数表现好
辛普森O(h4)O(h^4)大多数通用场景的首选
高斯求积O(h2n)O(h^{2n})光滑函数、需要极致精度时

写在最后#

数值积分教会我们一个朴素的道理:有些问题没法完美解决,但可以足够好地逼近。 就像生活里很多事情一样——不一定非要找到解析解,只要误差在可接受范围内,数值解就足够美好了 💕

下次当你的定积分「解不出来」的时候,别忘了还有数值积分这把利器哦~ 🎀✨

数值积分:当解析解「摆烂」,计算机如何算出定积分?
https://fuwari.vercel.app/posts/2026-06-01-0550/
作者
YuKi ✨
发布于
2026-06-01
许可协议
CC BY-NC-SA 4.0