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

2744557306

MATLAB代码展示了如何通过数值积分和符号积分计算函数的定积分，并对不同步长的数值积分结果进行比较。以下是逐步解析：
, ]; V& t' G- j8 |
' r2 ~5 @- b5 r7 ~### 1. 创建自变量和函数
% g; z+ s8 i8 X/ A4 B5 ?7 ^, i2 H```matlab
2 Y% @0 B/ H  E4 }& U, Sx = [0:0.01:3*pi/2, 3*pi/2];  % 包含 3π/2 这一点
8 g& B8 Z1 z! N0 Q* ~+ wy = cos(15*x);
- ~9 y3 T+ M; l: V' U, Gplot(x, y)                     % 绘制 y = cos(15*x) 的图形
9 ~1 n' l* y/ {% m& a& P```
- `x` 创建了一个从 0 到 \( \frac{3\pi}{2} \) 的向量，并确保 \( \frac{3\pi}{2} \) 这个点被包含在内。
+ Y, c5 K% f# {0 ]- `y` 计算了在 `x` 每个点上的 `cos(15*x)` 值。
- 使用 `plot(x, y)` 绘制了函数 `y = cos(15*x)` 的图形。
. J' p) N8 U/ f9 r8 d
### 2. 计算理论值
```matlab
syms x, A = int(cos(15*x), 0, 3*pi/2);   % 使用符号积分求取定积分
```
; D( N. B: ^( e; b# r& I0 K- X- `syms x` 定义符号变量 `x`。
& K# o* ]6 I0 o9 h0 C- `int(cos(15*x), 0, 3*pi/2)` 计算 `cos(15*x)` 从 0 到 \( \frac{3\pi}{2} \) 的定积分 `A`，这个结果为理论值。
2 `; Q, I3 }; I  x0 z& A
) U4 `/ |5 S' A, o% B### 3. 定义步长并进行数值积分
```matlab
h0 = [0.1, 0.01, 0.001, 0.0001, 0.00001, 0.000001];
9 P# F: F4 b2 e. P) U, D1 _. ]v = [];
. a2 S4 G8 M% c# x; Pfor h = h0,
  [6 N, t5 a8 _( N( T1 D# d    x = [0:h:3*pi/2, 3*pi/2];
    y = cos(15*x);
    I = trapz(x, y);            % 使用梯形法进行数值积分
4 m% y- s8 R: D$ y    v = [v; h, I, 1/15 - I];    % 将步长、数值积分结果和错误存储到 v 中
5 Q! I- }1 D! j9 k4 J# C6 ~) h3 o6 }end
```
- `h0` 是一个数组，包含了不同的步长值，用于进行数值积分。
% p% @* }; t  Z) H0 `% w" I2 Z9 ?- `v = []` 初始化一个空数组 `v`，用于存储每种步长的结果。
- 对于每个步长 `h`，代码：
  - 创建新的 `x` 向量，从 0 到 \( \frac{3\pi}{2} \)，步长为 `h`，并确保包含 \( \frac{3\pi}{2} \)。
  - 计算 `y = cos(15*x)`。
  - 使用 `trapz(x, y)` 进行数值积分，得到的结果存储在 `I` 中。
3 c, `& ^0 Y4 L; |5 D6 [/ G2 ^- v  - 将当前步长 `h`、计算得到的数值积分 `I` 和与理论值 \( \frac{1}{15} \) 的误差 \( 1/15 - I \) 存储到 `v` 中。
2 m9 E2 x4 l3 B5 p: E# {& H& Q
8 S9 `8 }; z( `! I/ A### 总结
这段代码通过图形显示、符号积分和数值积分来综合比较结果，展示了函数 `y = cos(15*x)` 的行为。通过不同的步长进行数值积分，并对结果与理论值进行比较，可以分析数值积分的精度随步长变化的情况。这是一种在数值分析中检验数值方法有效性的常见手段。
/ f2 E" [: u  y. x% ~4 p

# Z- ~  \( t1 O# |2 g