微分方程

微分方程为什么长成$\mathrm{d}y=f\mathrm{d}t$这副模样？本文从「求函数而非求数」的本质讲起，理清这套记号的莱布尼茨来源与 ODE→PDE→SDE 三百年的演化动机，最后用一张表收拢三类方程的谱系。

求的不是数，是函数

普通代数方程求的是数：$x^2=2$，解出来是一个数值。微分方程求的是函数——未知数是 $y(t)$，方程里出现它的导数 $\frac{\mathrm{d}y}{\mathrm{d}t}$、$\frac{{\mathrm{d}}^{2}y}{\mathrm{d}{t}^2}$。这一步抽象升级决定了整套理论的形态：解空间从实数轴变成了函数空间。

为什么自然规律偏要写成这种形式？因为物理、生物、经济里容易观测的往往是变化速率，而不是函数本身。传感器测的是速度和加速度，实验里读出的是反应速率，市场上看到的是价格的波动——直接拿到 $y(t)$ 表达式的场景反而罕见。

物体下落是标准例子，一个二阶常微分方程（阶数=方程中出现的最高阶导数的阶数）。牛顿第二定律不直接给位置 $x(t)$，给的是

$m\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{t}^2}=F$

已知力——也就是变化率层面的信息——反求位置函数 $x(t)$，这就是微分方程的任务。

核心思想

很多自然规律不告诉你「现在的值是多少」，只告诉你变化率是多少。已知导数约束、还原函数本身，是微分方程存在的全部理由。

$\mathrm{d}y=f\mathrm{d}t$ 这种写法从哪来

教科书的原始写法是

$\frac{\mathrm{d}y}{\mathrm{d}t}=f(y,t)$

莱布尼茨记号把导数看成两个无穷小量的比值：$\mathrm{d}y$ 是函数增量，$\mathrm{d}t$ 是时间的无穷小增量。既然是比值，交叉相乘就得到

$\mathrm{d}y=f(y,t)\mathrm{d}t$

这是无穷小视角：在极短时间 $\mathrm{d}t$ 之内，函数 $y$ 的改变量 $\mathrm{d}y$ 近似等于 $f(y,t)$ 乘以时间步长。两种写法在 ODE 里等价，但微分形式的适用范围更广——分水岭出现在引入噪声之后：$\mathrm{d}{W}_t$ 是随机无穷小，量级 $\sqrt{\mathrm{d}t}$，不是普通实数，古典微积分在此失效，详见随机微分方程与伊藤积分与伊藤公式。

实用直觉

$\mathrm{d}(\text{状态})=(\text{确定推动力量})\cdot \mathrm{d}t+(\text{随机扰动力量})\cdot \text{噪声增量}$ 读任何一个微分方程，先把它拆成这两股力量，再看各项分别是谁。

三百年演化：从决定论到噪声

ODE 与经典力学同生。17 世纪牛顿研究天体运动，牛顿第二定律加万有引力，行星轨道天然写成二阶常微分方程；牛顿本人用几何和级数方式求解，还顺手发明了早期微分工具。莱布尼茨建立微分记号 $\mathrm{d}$，才有了 $\mathrm{d}y/\mathrm{d}t$ 这套书写形式。到 18 世纪，欧拉大力发展 ODE——解法、存在唯一性、数值方法（欧拉法），大量力学、流体、单摆、振动问题被系统地转化成 ODE。实用动机很直接：我们知道「力决定加速度」，加速度是位置的二阶导数，却不知道位置的直接表达式，于是问题变成已知导数约束，还原函数。这一时期的 ODE 全部是确定性的——给定初值，未来完全确定，即拉普拉斯妖的思想。系统性的求解理论与数值方法见常微分方程。

18 世纪后期到 19 世纪，热传导、波动、电磁问题把空间放进了方程：状态不只随时间变，还随空间位置变化，未知函数成了 $u(x,y,z,t)$，方程里出现对空间的偏导数 $\partial u/\partial x$——这就是 PDE。傅里叶对热传导方程的研究是里程碑。

20 世纪轮到噪声登场。1828 年布朗观察到花粉的无规则运动；1905 年爱因斯坦用统计物理解释它并给出扩散方程；1923 年维纳给出布朗运动的严格数学构造（维纳过程 $W_t$）——但此时还没有 SDE。缺的最后一块拼图由伊藤清在 1940 年代补上：伊藤积分。布朗运动轨道处处不可微，$\frac{\mathrm{d}{W}_t}{\mathrm{d}t}$ 根本不存在，方程只能写成微分形式 $\mathrm{d}{X}_t=f\mathrm{d}t+g\mathrm{d}{W}_t$——这就是 SDE 标准形式的由来。

SDE 写不成导数形式

$\frac{\mathrm{d}{X}_t}{\mathrm{d}t}=\dots$ 这种写法对 SDE 不成立——布朗运动没有导数，只能保留 $\mathrm{d}{W}_t$ 的微分形式。这不是记号偏好，是数学上的硬约束。

SDE 的实用动机同样朴素：真实系统永远有扰动——分子碰撞、热噪声、市场冲击。ODE 只能描述理想的无噪声系统，而观测数据永远带涨落，于是需要 SDE 框架。整条时间线可以压缩成一行：ODE（1680s，力学，确定）→ PDE（1800s，场，时空）→ SDE（1940s，带随机噪声）。

三类方程对比

三类方程按未知函数的类型切分：依赖单个自变量的 ODE、依赖多个自变量的 PDE、取值为随机过程的 SDE。细节展开在各自笔记，一张表收拢全貌：

方程	未知函数	标准形式	一句话本质	深入
常微分方程 ODE	单变量函数 $y(t)$	$\frac{\mathrm{d}y}{\mathrm{d}t}=f(t,y)$	确定性轨迹：给定初值，未来唯一锁定	常微分方程
偏微分方程 PDE	多变量函数 $u(x,t)$	$u_t=D{u}_{xx}$	场的时空演化：一点的变化耦合整个空间	—
随机微分方程 SDE	随机过程 $X_t$	$\mathrm{d}{X}_t=\mu \mathrm{d}t+\sigma \mathrm{d}{W}_t$	轨道随机而密度确定：噪声驱动的动力学	随机微分方程

三者不是三门孤立的课：ODE 是 SDE 扩散项 $\sigma =0$ 的退化特例；SDE 轨道集合的概率密度又服从 Fokker‑Planck 方程——一个描述随机轨道、一个描述密度演化的 PDE，同一系统的两副面孔（见刘维尔方程 Liouville Equation的确定性对照）；PDE 做空间离散化后，每个格点各自变成一个 ODE。掌握这几条换算关系，三类方程就串成一张网。

PDE 三大经典方程

PDE 尚无独立笔记，经典款留在这里：

方程	形式	描述对象
波动方程	$u_{tt}=c^2{u}_{xx}$	波、弦振动
热方程	$u_t=D{u}_{xx}$	扩散、热传导
拉普拉斯方程	$u_{xx}+u_{yy}=0$	稳态静电场

PDE 同样分线性与非线性——纳维‑斯托克斯方程就是著名的非线性 PDE。

旁支：DDE、差分方程、SPDE

延迟微分方程 DDE：$\frac{\mathrm{d}y}{\mathrm{d}t}=f(y(t),y(t-\tau ))$，变化率依赖过去的状态，带时间延迟，用于传染病模型和控制系统。

差分方程不是微分方程——时间离散，$y_{n+1}=F(y_n)$，出现在仿真与迭代中；微分方程离散化之后得到的就是差分格式（如欧拉法）。

随机偏微分方程 SPDE：PDE 加上空间噪声，即场加随机扰动，用于湍流、随机场。

整张地图收成一棵树：

微分方程
├─ 确定性
│  ├─ ODE 常微分方程：单自变量（通常 t）
│  │   ├─ 线性 ODE（常系数 / 变系数，齐次 / 非齐次）
│  │   └─ 非线性 ODE
│  └─ PDE 偏微分方程：多个自变量（空间 + 时间）
│      ├─ 线性 PDE
│      └─ 非线性 PDE
└─ 随机
   ├─ SDE 随机微分方程：时间连续，状态随机（伊藤 / 斯特拉托诺维奇两种积分版本）
   └─ SPDE 随机偏微分方程：时空都带随机噪声

收尾：几条贯穿全文的线索

微分方程的本质是已知变化率，求函数——现实世界容易观测变化速率、难以直接得到函数表达式，这套工具因此诞生。$\mathrm{d}y=f\mathrm{d}t$ 来自莱布尼茨的无穷小记号；SDE 因为布朗运动轨道不可微而写不成导数形式，只能保留 $\mathrm{d}{W}_t$。历史顺序也映射着世界观的演进：ODE 出自牛顿力学的确定世界观，PDE 处理空间分布，SDE 面对无处不在的随机扰动。