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

2744557306

MATLAB代码展示了如何通过数值积分和符号积分计算函数的定积分，并对不同步长的数值积分结果进行比较。以下是逐步解析：
# |, T* W0 c# E; h: `! K0 _
% |3 D  s, i9 u### 1. 创建自变量和函数
```matlab
x = [0:0.01:3*pi/2, 3*pi/2];  % 包含 3π/2 这一点
y = cos(15*x);
( L- X; x8 P5 S+ `- {7 v* F( fplot(x, y)                     % 绘制 y = cos(15*x) 的图形
/ K% w, e% P/ P0 D7 H  E```
! p+ E2 V0 N( y- `x` 创建了一个从 0 到 \( \frac{3\pi}{2} \) 的向量，并确保 \( \frac{3\pi}{2} \) 这个点被包含在内。
- `y` 计算了在 `x` 每个点上的 `cos(15*x)` 值。
5 j% ~! p: a6 D  r6 h: t: G8 P+ |- 使用 `plot(x, y)` 绘制了函数 `y = cos(15*x)` 的图形。
( m0 n- y7 M0 b: ~* |
### 2. 计算理论值
```matlab
syms x, A = int(cos(15*x), 0, 3*pi/2);   % 使用符号积分求取定积分
```
- h4 _& W/ _# E  J: \& h! r- `syms x` 定义符号变量 `x`。
7 q  f: l; d6 z$ n, L+ J- `int(cos(15*x), 0, 3*pi/2)` 计算 `cos(15*x)` 从 0 到 \( \frac{3\pi}{2} \) 的定积分 `A`，这个结果为理论值。

### 3. 定义步长并进行数值积分
, s' `- @, S! g3 w( B, F```matlab
. R0 \% X) [2 W# xh0 = [0.1, 0.01, 0.001, 0.0001, 0.00001, 0.000001];
: S8 ^+ A% W: e( q. I1 vv = [];
. T) v2 B6 u; S8 C4 D7 Nfor h = h0,
    x = [0:h:3*pi/2, 3*pi/2];
" Y& E0 {- g1 D" D    y = cos(15*x);
    I = trapz(x, y);            % 使用梯形法进行数值积分
0 n7 _1 k3 k, v9 a    v = [v; h, I, 1/15 - I];    % 将步长、数值积分结果和错误存储到 v 中
, p+ U  I7 v- t! E( ?) D  b- Fend
# H8 \* ^7 m0 Z5 M% G5 q) R```
  ?3 C. }$ v% P6 R- `h0` 是一个数组，包含了不同的步长值，用于进行数值积分。
( i2 z# t( p$ `, f% \# G$ g- `v = []` 初始化一个空数组 `v`，用于存储每种步长的结果。
6 i) z0 h& {+ f- 对于每个步长 `h`，代码：
; ]5 t' W( m: f4 ~1 y0 I7 ~% H  - 创建新的 `x` 向量，从 0 到 \( \frac{3\pi}{2} \)，步长为 `h`，并确保包含 \( \frac{3\pi}{2} \)。
. [3 l* I/ b8 H8 i  b$ c  - 计算 `y = cos(15*x)`。
  - 使用 `trapz(x, y)` 进行数值积分，得到的结果存储在 `I` 中。
$ ]. T3 o0 F! A3 }  - 将当前步长 `h`、计算得到的数值积分 `I` 和与理论值 \( \frac{1}{15} \) 的误差 \( 1/15 - I \) 存储到 `v` 中。

9 `2 k3 E$ d8 N; ]$ m### 总结
这段代码通过图形显示、符号积分和数值积分来综合比较结果，展示了函数 `y = cos(15*x)` 的行为。通过不同的步长进行数值积分，并对结果与理论值进行比较，可以分析数值积分的精度随步长变化的情况。这是一种在数值分析中检验数值方法有效性的常见手段。

0 t2 R6 E8 V# n8 m- y+ S, g' A1 i6 t
9 o* a3 N* l4 j6 C" o$ @# \
- i+ j  F. x4 X! z! e