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

2744557306

MATLAB代码展示了如何通过数值积分和符号积分计算函数的定积分，并对不同步长的数值积分结果进行比较。以下是逐步解析：
. n9 w0 Z% i, W* B. o2 \
2 ?. T9 P4 x9 ^; H7 l3 X### 1. 创建自变量和函数
2 D6 p/ ~2 ~1 i, P```matlab
x = [0:0.01:3*pi/2, 3*pi/2];  % 包含 3π/2 这一点
' M! K: Z  L0 p$ Cy = cos(15*x);
9 ]+ c0 U- Z; {plot(x, y)                     % 绘制 y = cos(15*x) 的图形
```
- `x` 创建了一个从 0 到 \( \frac{3\pi}{2} \) 的向量，并确保 \( \frac{3\pi}{2} \) 这个点被包含在内。
5 L. r" A  ]5 j- `y` 计算了在 `x` 每个点上的 `cos(15*x)` 值。
. t1 x; D' D' \: Y- 使用 `plot(x, y)` 绘制了函数 `y = cos(15*x)` 的图形。

  h* m) M: A' e0 [, c### 2. 计算理论值
- _3 q: }) R( M2 C( q( N2 _```matlab
syms x, A = int(cos(15*x), 0, 3*pi/2);   % 使用符号积分求取定积分
$ j* O; s8 h: X2 t4 T" y```
) d$ }1 B8 n4 R, |) p! g% G! I- `syms x` 定义符号变量 `x`。
- `int(cos(15*x), 0, 3*pi/2)` 计算 `cos(15*x)` 从 0 到 \( \frac{3\pi}{2} \) 的定积分 `A`，这个结果为理论值。

  i. N2 L, Y9 R3 ]; K, k& T' v### 3. 定义步长并进行数值积分
```matlab
# D/ m' }! K/ J8 Q) `) K. t1 {9 Vh0 = [0.1, 0.01, 0.001, 0.0001, 0.00001, 0.000001];
! L9 v( r) v2 E- c5 wv = [];
% f( G( e( O9 h+ [2 n" }3 \for h = h0,
    x = [0:h:3*pi/2, 3*pi/2];
% }$ s0 z, ~" b* Q( N% j/ {7 u" E# s    y = cos(15*x);
    I = trapz(x, y);            % 使用梯形法进行数值积分
    v = [v; h, I, 1/15 - I];    % 将步长、数值积分结果和错误存储到 v 中
end
```
- `h0` 是一个数组，包含了不同的步长值，用于进行数值积分。
2 Y3 c' i; b7 x1 m6 k* R' m- `v = []` 初始化一个空数组 `v`，用于存储每种步长的结果。
- 对于每个步长 `h`，代码：
7 B9 a. @. k& ?# i  N: ^  - 创建新的 `x` 向量，从 0 到 \( \frac{3\pi}{2} \)，步长为 `h`，并确保包含 \( \frac{3\pi}{2} \)。
  - 计算 `y = cos(15*x)`。
* A8 ^. h$ |$ z" `7 Z  - 使用 `trapz(x, y)` 进行数值积分，得到的结果存储在 `I` 中。
  - 将当前步长 `h`、计算得到的数值积分 `I` 和与理论值 \( \frac{1}{15} \) 的误差 \( 1/15 - I \) 存储到 `v` 中。

3 N5 X: C4 D4 l3 l8 o$ O; s- K### 总结
% F: }8 T9 ^$ _4 t! t6 }这段代码通过图形显示、符号积分和数值积分来综合比较结果，展示了函数 `y = cos(15*x)` 的行为。通过不同的步长进行数值积分，并对结果与理论值进行比较，可以分析数值积分的精度随步长变化的情况。这是一种在数值分析中检验数值方法有效性的常见手段。
& v; @4 b2 M1 r% m8 ]1 a

5 ?1 O+ G$ ]+ q7 F
1 L1 ]8 |" K/ Z  Z8 e) S( x; ~3 N* g