49、数值积分方法与MATLAB实现

数值积分方法与MATLAB实现

1. 数值积分规则

1.1 Simpson 1/3 规则

Simpson 1/3 规则的公式及误差表达式如下：
[
I = \frac{h}{3}\left[f(a) + 4\sum_{i = 1, i\text{ odd}}^{n - 1}f(x_i)+2\sum_{i = 2, i\text{ even}}^{n - 2}f(x_i)+f(b)\right]+E_{net}
]
[
E_{net}=-\frac{(b - a)h^4}{180}f^{(4)}(\xi)
]
其中 (h=\frac{b - a}{n})，全局截断误差（GTE）与 (h^4) 成正比，即 (O(h^4))。

1.2 Simpson 3/8 规则

可以推导 Simpson 3/8 规则的多重应用公式，有兴趣的读者可阅读数值技术的高级文献，了解使用更多节点以提供更高精度的高阶 Newton - Cotes 积分公式。

1.3 梯形规则、Simpson 1/3 规则和 Simpson 3/8 规则的比较

为了比较这三种方法，选择 (n = 6)，则 (h=\frac{1}{6})，这对应于梯形规则应用 6 次，Simpson 1/3 规则应用 3 次，Simpson 3/8 规则应用 2 次。具体操作步骤如下：
1. 计算 (h=\frac{b - a}{n})；
2. 生成向量 (xVec=[a:h:b])；
3. 计算函数值 (fVec = fun4Int(xVec))；
4. 使用相应公