通过数值积分和符号积分计算函数的定积分

2744557306

MATLAB代码展示了如何通过数值积分和符号积分计算函数的定积分，并对不同步长的数值积分结果进行比较。以下是逐步解析：
4 z; }1 H4 r: L. Z  e! }
### 1. 创建自变量和函数
) R' e2 j% k1 p- v: F3 {7 a' [```matlab
& L$ |" d& j5 r4 j7 t9 w2 ^* `x = [0:0.01:3*pi/2, 3*pi/2];  % 包含 3π/2 这一点
y = cos(15*x);
' q. s; e9 E  Q& K2 @plot(x, y)                     % 绘制 y = cos(15*x) 的图形
$ B3 Q- r) p' j8 j* T5 v```
3 V" r$ A2 _( W3 f2 \# `- `x` 创建了一个从 0 到 \( \frac{3\pi}{2} \) 的向量，并确保 \( \frac{3\pi}{2} \) 这个点被包含在内。
- `y` 计算了在 `x` 每个点上的 `cos(15*x)` 值。
- 使用 `plot(x, y)` 绘制了函数 `y = cos(15*x)` 的图形。
0 i) m7 I/ k- ?' g6 d
### 2. 计算理论值
```matlab
3 b+ m4 }8 a% ysyms x, A = int(cos(15*x), 0, 3*pi/2);   % 使用符号积分求取定积分
7 P7 D4 W6 y- K0 f( f( w$ P1 g& L7 }$ o```
- `syms x` 定义符号变量 `x`。
- `int(cos(15*x), 0, 3*pi/2)` 计算 `cos(15*x)` 从 0 到 \( \frac{3\pi}{2} \) 的定积分 `A`，这个结果为理论值。
  P- g; \( J5 Y- Z! i- S* J$ e$ y5 L
### 3. 定义步长并进行数值积分
```matlab
2 V/ i, \2 k! e1 `5 V( O5 Vh0 = [0.1, 0.01, 0.001, 0.0001, 0.00001, 0.000001];
5 D, s* k% i* R) Fv = [];
for h = h0,
$ y; l" ]4 b- S    x = [0:h:3*pi/2, 3*pi/2];
    y = cos(15*x);
    I = trapz(x, y);            % 使用梯形法进行数值积分
    v = [v; h, I, 1/15 - I];    % 将步长、数值积分结果和错误存储到 v 中
end
```
( a8 @, t. L" j; d3 J+ L- `h0` 是一个数组，包含了不同的步长值，用于进行数值积分。
- `v = []` 初始化一个空数组 `v`，用于存储每种步长的结果。
- 对于每个步长 `h`，代码：
& x1 z9 @. _; S7 p7 A# `/ x  - 创建新的 `x` 向量，从 0 到 \( \frac{3\pi}{2} \)，步长为 `h`，并确保包含 \( \frac{3\pi}{2} \)。
" ?0 H' K1 b. ?8 _+ Y  x  - 计算 `y = cos(15*x)`。
  - 使用 `trapz(x, y)` 进行数值积分，得到的结果存储在 `I` 中。
  - 将当前步长 `h`、计算得到的数值积分 `I` 和与理论值 \( \frac{1}{15} \) 的误差 \( 1/15 - I \) 存储到 `v` 中。
( M# r+ S8 C, M9 o9 l, p! W+ x
! w+ L; ^0 h% n5 W' Q### 总结
; D% t  A7 t9 E, ]( I# u6 C这段代码通过图形显示、符号积分和数值积分来综合比较结果，展示了函数 `y = cos(15*x)` 的行为。通过不同的步长进行数值积分，并对结果与理论值进行比较，可以分析数值积分的精度随步长变化的情况。这是一种在数值分析中检验数值方法有效性的常见手段。
: D0 t# A6 m- W  T4 l. b1 K% y$ v
! c7 m2 Z  Y* C: c. x# h3 D: b+ }7 {
+ ~: b, ^1 l& E, [9 }3 R5 ~