matlab用辛普森公式求积分_数值积分常用方法

数值积分的基本思想

由积分中值定理可知，在积分区间

内存在一点

，成立

式的几何意义即为：底为

而高为

的矩形的面积恰等于所求曲边梯形的面积

。因此，要想求出

式左端积分，我们只需要知道三个值：

即可，这里

是显然的，问题在于

的具体位置一般是不清楚的，从而

未知。我们暂且将

称为区间

上的平均高度，我们的目标就是寻求一种求出平均高度

的算法，这样我们的问题也就解决了。

---

首先我们想到的是用区间两端点的“高度”

与

的算术平均值作为平均高度

的近似值，从而导出下述求积公式

式即为我们我们熟知的

梯形公式。写出MATLAB函数为

function

而若用区间中点

的“高度”

近似取代平均高度

，则可导出

中矩形公式（简称 矩形公式）

其MATLAB函数为

function

总结：更一般地，我们可以在区间

上适当选取某些节点

，然后用

的加权平均得到平均高度