matlab梯形法数值积分,数值积分（梯形法）

&nbsp

&nbsp

&nbsp

&nbsp

&nbsp

&nbsp

&nbsp

预备知识　Matlab 的函数，定积分

本词条介绍一种简单的梯形算法计算数值积分．如图 1 若要用梯形法计算定积分 $\int_a^b f(x) \,\mathrm{d}{x} $，则可将区间 $[a, b]$ 划分为 $N$ 个长度为 $h = (b-a)/N$ 的等长的小区间，区间端点从 $a$ 到 $b$ 分别为 $x_1 = a, x_2 = a + h, \dots, x_{N+1} = b$．

图 1：梯形法数值积分

接下来将每个区间的被曲线围出的面积用梯形来计算，第 $i$ 个梯形面积为(两底和乘以高除以二)

$[f(x_{i+1}) + f(x_i)]h/2$．定积分约等于所有梯形面积

\begin{equation}

\int_a^b f(x) \,\mathrm{d}{x} \approx \sum_i^N \frac h2 [f(x_{i+1}) + f(x_i)]

= h \left[\frac12 f(x_1) + \sum_{i = 2}^N f(x_i) + \frac12 f(x_{N+1}) \right]

\end{equation}

显然&#x