4. 梯形公式与Simpson公式的误差分析
在数值积分中，复合梯形公式与复合Simpson公式是最常用的两种方法。它们都通过把积分区间分成若干小区间、在每个小区间上用低阶多项式近似被积函数来实现积分。但两者的误差来源、光滑性要求和收敛速度差别很大。理解这些差别，是选择算法和评估结果可信度的基础。
一、截断误差表达式的推导思路
先看单个小区间上的局部误差。设被积函数 f 在区间 [a,b] 上足够光滑，记 h=b−a。梯形公式用连接两端点的直线代替曲线，其局部截断误差可写成
E_T = −(h³/12)·f″(ξ), ξ∈(a,b)推导的关键是把 f 在区间中点展开成泰勒级数，再对线性插值多项式积分，二次项及更高次项无法被直线精确积分，剩余部分即给出上述误差。误差与 h 的三次方成正比，且符号由二阶导数决定：f 为凸函数时梯形公式偏大，凹函数时偏小。
Simpson公式在每个小区间上改用经过左端点、中点、右端点的二次插值多项式。其局部误差为
E_S = −(h⁵/2880)·f⁽⁴⁾(ξ), ξ∈(a,b)注意这里 h 是子区间长度，Simpson公式需要三个节点，因此实际使用时通常把 [a,b] 分成偶数个小区间，每两个小区间组成一个“双区间”套用一次。误差中出现四阶导数，说明它对函数光滑性的要求更高：只有 f 具有连续四阶导数时，这个表达式才成立。
二、光滑性假设的影响
上述误差公式都建立在“f 在区间上具有足够阶连续导数”的前提上。这个假设不是形式上的装饰，它直接决定误差估计是否可信。
- 若 f 只有连续二阶导数，梯形公式的误差表达式仍然成立，但 Simpson 公式的四阶导数可能不存在，此时 E_S 的经典表达式失效，实际收敛可能明显慢于理论值。
- 若 f 在区间内存在尖点、跳跃或导数剧烈变化（例如 |x|、分段函数、含激波或边界层的函数），高阶导数会变得很大甚至无界，误差公式中的 f⁽⁴⁾(ξ) 无法给出有意义的界。
- 此时 Simpson 公式的“高精度”优势可能消失，甚至因为过冲而比梯形公式更差；梯形公式虽然阶低，但对不光滑函数的鲁棒性通常更好。
三、代数精度与收敛阶的比较
代数精度指公式能精确积分的最高多项式次数。梯形公式对一次多项式精确，代数精度为 1；Simpson公式对三次多项式精确，代数精度为 3。这解释了为什么 Simpson 公式的局部误差从 h³ 提升到 h⁵：它“免费”多精确匹配了两阶多项式信息。
把局部误差累加到整个区间。设区间长度为 L，分成 n 个长度为 h=L/n 的小区间，复合梯形公式的总误差为
|E_T| ≤ (L·h²/12)·max|f″| → O(h²)复合 Simpson 公式（n 为偶数，每两个小区间一组）的总误差为
|E_S| ≤ (L·h⁴/180)·max|f⁽⁴⁾| → O(h⁴)因此梯形公式是二阶收敛，Simpson 公式是四阶收敛。网格加密一倍时，梯形误差约降为原来的 1/4，Simpson 误差约降为 1/16。这正是 Simpson 公式在光滑问题上更受欢迎的原因。
四、如何用误差公式做实际判断
误差公式的价值不只是给出一个界，更在于指导我们如何验证计算结果。一个实用的做法是“网格加密对比法”：用 n 和 2n 两套网格分别计算积分值 In 与 I2n，观察差值 I2n−In 随 h 缩小的衰减速度。
- 若差值大致按 1/4 衰减，说明处于梯形公式的 O(h²) 收敛区，误差公式中的二阶导数假设基本成立。
- 若差值大致按 1/16 衰减，说明处于 Simpson 公式的 O(h⁴) 收敛区，四阶导数假设可信。
- 若衰减速度明显慢于理论阶，或出现不规则的跳动，往往意味着函数在该区间不光滑，或网格尚未进入渐近收敛区，此时应缩小积分区间、拆分奇异点，或改用自适应积分。
还可以利用误差公式做“事后修正”。对 Simpson 公式，若已知 In 与 I2n，可用 Richardson 外推构造更高阶的近似：I≈(16I2n−In)/15。这一修正之所以有效，正是因为误差按 h⁴ 的规律变化；若函数不光滑，外推反而可能放大噪声。
总结来看，梯形公式与 Simpson 公式的差别不只是“精度高低”，而是误差结构对光滑性的依赖程度不同。梯形公式以较低的光滑性要求换取稳健的二阶收敛；Simpson 公式以更高的光滑性假设换取四阶收敛和更高的代数精度。实际选择时，应先判断被积函数的光滑性，再用网格加密检验实际收敛阶是否与理论一致，最后才决定是否采用外推或自适应策略。