26、MATLAB在微分方程与边界值问题求解中的应用

MATLAB在微分方程与边界值问题求解中的应用

1. MATLAB求解微分方程概述

在MATLAB中，我们开发了一系列用于求解微分方程和微分方程组的函数，这些函数对MATLAB自带的功能起到了补充作用。同时，还简要介绍了MATLAB的Simulink工具箱。对于拥有Simulink的用户，在MATLAB命令窗口输入 >> shower ，可以找到一个优秀的Simulink模型示例，这是一个使用模糊逻辑的淋浴温度控制系统的模拟，该模型有四个掩码块，例如热水阀块本身包含七个低级块。

2. 微分方程求解问题集

2.1 放射性物质衰变问题

放射性物质的衰变率与剩余量成正比，其微分方程为 $\frac{dy}{dt} = -ky$ ，当 $t = t_0$ 时，$y = y_0$ 。给定 $y_0 = 50$ ，$k = 0.05$ ，在 $t = 0$ 到 $10$ 的范围内求解该方程：
- 使用 feuler 函数，步长 $h$ 分别取 $1$、$0.1$、$0.01$ 。
- 使用 eulertp 函数，步长 $h$ 分别取 $1$、$0.1$ 。
- 使用 rkgen 函数（经典方法），步长 $h = 1$ 。
并将结果与精确解 $y = 50\exp(-0.05t)$ 进行比较。

2.2 一阶微分方程求解

求解 $y’ = 2xy$ ，初始条件为 $x_0 = 0$ 时，$y_0 = 2$ ，在 $x = 0$ 到 $2$ 的范围内。使用