21、常微分方程的数值解法探究

常微分方程的数值解法探究

1. 物理学中的微分方程

微分方程在物理学中占据核心地位，许多物理方程都是微分方程，例如牛顿第二定律、薛定谔方程和麦克斯韦方程组等，其中大部分属于偏微分方程，但在某些情况下可简化为常微分方程。

1.1 常见的常微分方程

• 牛顿第二定律（一维运动） ：$F = m\frac{d^2x}{dt^2}$，这是一个二阶常微分方程。
• 小振幅单摆 ：$\ddot{\theta}=-\frac{mg}{l}\theta$，同样是二阶常微分方程。
• 一维定态薛定谔方程 ：$-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2}+V(x)\psi = E\psi$，也是二阶常微分方程。

这些方程大多是线性微分方程，但实际上自然界并非完全线性。以单摆为例，其实际方程为$\ddot{\theta}=-\frac{mg}{l}\sin\theta$，右侧包含因变量的正弦函数，是高度非线性的。只有在小振幅条件下，该方程才简化为线性方程。

1.2 线性与非线性方程

传统上，物理学家和数学家更关注线性方程，因为它们易于求解且具有一些非线性方程所没有的便利性质，如叠加原理仅适用于线性微分方程。这也导致非线性动力学和混沌领域直到20世纪后半叶才得到发展，因为大多数非线性微分方程只能用数值方法求解，需要高速计算机的支持。

1.3 非线性方程的应用

如果不局限于线性方程，许多