[ 常微分方程 ] 03 微分方程数值解

伊利亚家的猫

已于 2024-11-06 16:01:45 修改

阅读量1.4k

点赞数 15

分类专栏：数学文章标签：数学建模

于 2024-10-05 22:55:50 首次发布

本文链接：https://blog.csdn.net/2201_75660982/article/details/142719638

版权

今天学习常微分方程的数值解，学习视频指路：微分方程数值解，感谢本视频让我突击学会了数值解（教材真的看不懂www）。本文作为该视频的学习笔记，加入了个人的一些理解，同时参考了王高雄老师的《常微分方程》第四版，如有理解不当地方还请各位指正。

文章目录

一、欧拉方法（Euler法）
- 01 Euler法基本格式
- 02 差分形式拓展
- 03 Euler法的改进
- - （1）左矩形公式
  - （2）右矩形公式
  - （3）梯形公式
- 04 实际运行效果
- 05 欧拉法的评价
二、龙格—库塔方法（Runge—Kutta)
- 01 龙格—库塔方法基本介绍
- 02 相容性条件
- 03 龙格—库塔方法的评价

一、欧拉方法（Euler法）

01 Euler法基本格式

对于微分方程的初值问题 $\begin{cases}\frac{dy}{dx}=f(x,y)\\ y(x_0)=y_0\end{cases}$ 数值方法本质上就是希望能够找出一系列点去逼近 $y = f (x)$ 这条曲线。

微分方程的近似
假设这个微分方程的解曲线在 $[0, T]$ 上有定义，我们将区间划分为 $n$ 个区间，即 $n + 1$ 个节点，并令步长 $h=\frac{T}{n}$ ,记节点为 $t_k$ ,满足 $0<t_0<t_1<...<t_n=T$
为了方便推导符号的统一性（其实用 $y$ 也是可以的），我们记初值问题为 $\begin{cases}u(t_{0})=u_{0}\\u'(t_{0})=f(t_0,u_{0})\end{cases}$ 在下一个节点处进行泰勒展开有 $u(t_1)=u(t_0+h)=u(t_0)+u'(t_0)h+O(h^2)$ 忽略二阶小量，令一个新的符号表示近似值 $u_1=u(t_0)+u'(t_0)h$ 由于在递推过程中，精确解本身就是难以得到的，所有的逼近过程都是用近似值 $u_1=u_0+f(t_0,u_0)h\approx u(t_1)$ 以此类推，得到递推公式
$u_{i+1}=u_i+f(t_{i},u_{i})h\quad (i=0,1,...,n-1)\tag{1-1}$