用数值积分的两种方法来计算函数的曲线下的面积

2744557306

这段MATLAB代码展示了如何使用数值积分的两种方法来计算函数的曲线下的面积，并确保它们得到一致的结果。以下是代码的逐步解析：

5 J+ t: a7 O; A& Q( h" r( W/ ^### 1. 定义区间和函数值
& W+ K/ h' h$ }# J% r' D& r```matlab
x1 = [0:pi/30:pi]';
y = [sin(x1) cos(x1) sin(x1/2)];
```
- `x1` 创建了一个列向量，从 0 到 π，以 π/30 的步长划分。
- `y` 是由三个函数（`sin(x1)`、`cos(x1)` 和 `sin(x1/2)`）组成的矩阵，每一列对应一个函数在 `x1` 上的值。

### 2. 准备进行数值积分的变量
```matlab
# U0 U+ p+ J! G4 F; ?x = [x1 x1 x1];
3 i3 W" L2 L2 ^" ~2 e- A# @```
- `x` 将列向量 `x1` 复制三次，形成一个三列的矩阵。每列表示不同的 `x` 值，目的可能是为了与 `y` 矩阵相对应。
( o: K* [9 m6 p# ?
4 |% t; a& f, p- m### 3. 使用梯形法计算面积
```matlab
2 e8 ?; ^9 P2 r, RS = sum((2*y(1:end-1, + diff(y)) .* diff(x)) / 2;
  k( t0 F: W7 }- i7 Z+ H```
- `diff(y)` 计算 `y` 的相邻差分，得到每个函数在相邻采样点间的变化。
' Q3 P+ V# A6 o% f$ O( N. z- `y(1:end-1,` 取 `y` 的前 n-1 行，作为计算梯形法时的基准。
: m( h, u+ O' R; U- `2*y(1:end-1, + diff(y)` 计算了每个小梯形的高度，`diff(x)` 计算了相邻 `x` 的差分，表示每个小段的宽度。
( v# u. W# J- `3 J- 最后，通过求和并除以 2 来得出总面积 \( S \)。

### 4. 使用 `trapz` 函数进行数值积分
$ X( [2 V& z  `% S0 R# r5 d" N```matlab
1 l% G, Y+ r  B/ m" yS1 = trapz(x1, y);
; Q; ]% n. s# |  W0 ^```
- `trapz(x1, y)` 使用 MATLAB 的内置函数 `trapz` 来计算由 `x1` 和 `y` 定义的曲线下的面积。这也采用了梯形法进行积分。
! {, z0 n! h0 L  |
### 5. 一致性检查
- 变量 `S` 和 `S1` 都存储了通过不同积分方法得到的面积值。由于使用了相同的函数和区间，理论上这两个值应该是相等的。
! G6 ?3 v9 k7 A$ G. i2 O# W+ k
6 q3 _1 o/ m; f( _# h### 总结
' V4 }* a1 |- J- w' c' C+ F这段代码有效地演示了如何使用手动实施的梯形法和 MATLAB 内置的 `trapz` 函数来计算两个或多个函数下的区域面积。通过比较两个结果，用户可以验证数值积分的正确性和精度。在数值分析中，这种方法用于确保计算的可靠性非常重要。
$ R. K# m/ E" q5 B: s