数学建模社区-数学中国

标题: 通过泰勒级数展开绘制了正弦函数的逼近 [打印本页]

作者: 2744557306 时间: 2024-8-27 11:05
标题: 通过泰勒级数展开绘制了正弦函数的逼近

x0=-2*pi:0.01:2*pi; y0=sin(x0); syms x; y=sin(x);& X8 D9 C. n, y5 h: F9 s
plot(x0,y0), axis([-2*pi,2*pi,-1.5,1.5]); hold on
' K" |% s1 r2 A- O7 v8 W
for n=[8:2:16]
- F+ y1 T( i l# T$ e
p=taylor(y,x,n), y1=subs(p,x,x0); line(x0,y1)
$ G, C. K7 [$ U2 E) B5 j
end% D" q9 `0 u; ?% r

复制代码

这段 MATLAB 代码通过泰勒级数展开绘制了正弦函数的逼近。下面是代码的逐步解释：

### 代码解释

1. **定义 x 范围**：
```matlab
x0 = -2*pi:0.01:2*pi;
```
- 这行代码创建了一个从 $-2\pi$ 到 $2\pi$ 的向量 `x0`，步长为 0.01。这个向量将用于计算和绘图。

2. **计算 sin(x0)**：
```matlab
y0 = sin(x0);
```
- 计算 `x0` 中每个值的正弦，并将结果存储在 `y0` 中。`y0` 将是 `sin` 函数的实际值，用于绘图。

3. **定义符号变量和函数**：
```matlab
syms x;
y = sin(x);
```
- 使用 `syms` 创建符号变量 `x`，然后定义符号函数 $ y = \sin(x) $。这个函数用于后续的泰勒级数展开。

4. **绘制 sin(x) 图形**：
```matlab
plot(x0, y0), axis([-2*pi, 2*pi, -1.5, 1.5]); hold on
```
- 使用 `plot` 函数绘制 `y0` 关于 `x0` 的图形，即实际的正弦波。
- `axis` 函数设置坐标轴的范围为 $[-2\pi, 2\pi]$ 和 $[-1.5, 1.5]$。
- `hold on` 使得后续绘图不会覆盖当前的图形。

5. **进行泰勒级数展开和绘图**：
```matlab
for n = [8:2:16]
   p = taylor(y, x, n);
   y1 = subs(p, x, x0);
   line(x0, y1)
end
```
- 使用 `for` 循环遍历 `n` 的值，从 8 到 16，步长为 2（即分别为 8、10、12、14 和 16）。
- 在循环内部：
   - `p = taylor(y, x, n)` 计算在点 0 附近的 $ n $ 次泰勒级数展开，得到多项式 $ p $。
   - `y1 = subs(p, x, x0)` 将泰勒展开多项式 $ p $ 替换中 `x` 的值为 `x0`，以计算对应的 `y1`（即泰勒多项式的值）。
   - `line(x0, y1)` 在当前图中绘制泰勒级数的结果。

### 效果

- 代码运行后，会得到一幅包含原始正弦函数图像和不同阶次的泰勒多项式的图形。每个泰勒多项式的图形与正弦函数重合得越近，表示这一级数的逼近效果越好。

### 知识点总结

1. **泰勒级数**：
- 泰勒级数是表示函数的一种多项式近似，适用于在某一点附近的函数描述。

2. **符号计算**：
- 使用 MATLAB 的符号工具箱，能够对符号函数进行解析计算并获取多项式形式。

3. **绘图与数据可视化**：
- `plot` 和 `line` 函数用于展示函数图像，`hold on` 功能允许在同一图中叠加多个图形。

4. **遍历与替换**：
- 通过循环和 `subs` 函数，可以对多次定义和计算的函数值进行有效处理。

### 结论

这段代码展示了如何利用 MATLAB 对正弦函数进行泰勒级数展开，并通过可视化的方式展示其近似效果。可以通过这个示例了解泰勒级数的适用性和效率，同时为函数逼近与数值计算提供了直观的理解。

examp3_18.m

211 Bytes, 下载次数: 0, 下载积分: 体力 -2 点

售价: 2 点体力 [记录] [购买]

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)