常微分方程（Ordinary Differential Equations, ODE）是描述未知函数及其导数之间关系的方程，其中未知函数仅依赖于一个自变量。ODE 是数学建模的核心工具，广泛应用于物理、工程、生物、经济等各个领域。本文系统梳理 ODE 的基础理论、求解方法、定性分析及数值计算，从入门到深入，涵盖核心概念与实用技巧。
**常微分方程（ODE）**是包含一个自变量 x、未知函数 y(x) 及其导数 y′,y′′,…,y(n) 的方程：
F(x,y,y′,y′′,…,y(n))=0
方程中出现的最高阶导数的阶数称为方程的阶。
- 一阶 ODE：y′=f(x,y)
- 二阶 ODE：y′′+p(x)y′+q(x)y=f(x)
- n 阶 ODE：y(n)=f(x,y,y′,…,y(n−1))
线性 ODE：可写成以下形式的 ODE
an(x)y(n)+an−1(x)y(n−1)+⋯+a1(x)y′+a0(x)y=g(x)
特征：未知函数及其导数均为一次，无乘积项 y⋅y′ 或 (y′)2 等。
非线性 ODE：不满足线性条件的方程。常见的非线性类型包括：
- 伯努利方程：y′+p(x)y=q(x)yn
- 黎卡提方程：y′=p(x)y2+q(x)y+r(x)
- 一阶隐式方程：F(x,y,y′)=0
非线性 ODE 通常没有通用解析解，数值方法和定性分析成为主要手段。
- 齐次：g(x)=0
- 非齐次：g(x)=0
- 初值问题（IVP，Initial Value Problem）：给定 x=x0 处 y 及其各阶导数的值
- 边值问题（BVP，Boundary Value Problem）：给定区间端点处 y 或其导数的条件
形式：dxdy=g(x)h(y)
解法：将变量分离后两边积分
h(y)dy=g(x)dx⇒∫h(y)dy=∫g(x)dx+C
典型例子：指数增长模型
dtdP=rP⇒P(t)=P0ert
形式：dxdy=f(xy)
解法：令 u=y/x，则 y=ux，y′=u+xu′，代入化简为可分离变量方程
形式：y′+p(x)y=q(x)
解法（积分因子法）：
- 计算积分因子 μ(x)=e∫p(x)dx
- 两边乘以 μ(x)：μ(x)y′+μ(x)p(x)y=μ(x)q(x)
- 左边是 (μy)′：(μy)′=μq
- 积分得：y=μ(x)1[∫μ(x)q(x)dx+C]
公式总结：
y=e−∫p(x)dx[∫e∫p(x)dxq(x)dx+C]
典型应用：RC 电路中的电流
RdtdI+C1I=dtdV
形式：y′+p(x)y=q(x)yn，n=0,1
解法：令 u=y1−n，化为 u′+(1−n)p(x)u=(1−n)q(x)，即一阶线性方程
形式：M(x,y)dx+N(x,y)dy=0
恰当条件：∂y∂M=∂x∂N
解法：存在 ψ(x,y) 使得 dψ=Mdx+Ndy，方程的解为 ψ(x,y)=C
非恰当方程的积分因子：若恰当条件不满足，可寻找 μ(x,y) 使 μMdx+μNdy=0 成为恰当方程。
形式：y=xy′+ψ(y′)
解法：令 p=y′，两边对 x 求导得 p=p+xp′+ψ′(p)p′，即 [x+ψ′(p)]p′=0
- p′=0 得通解 y=Cx+ψ(C)
- x+ψ′(p)=0 得奇解（包络线）
n 阶线性 ODE 的一般形式：
y(n)+p1(x)y(n−1)+⋯+pn−1(x)y′+pn(x)y=f(x)
解的结构：
- 齐次解（补函数）yh：对应齐次方程的通解，由 n 个线性无关的解构成
- 特解yp：非齐次方程的一个特解
- 通解：y=yh+yp
解的存在唯一性定理：若 pi(x) 和 f(x) 在区间 I 上连续，则对任意 x0∈I 和给定的初值，IVP 存在唯一解。
形式：y(n)+a1y(n−1)+⋯+an−1y′+any=0
解法：设 y=erx，代入得特征方程 rn+a1rn−1+⋯+an−1r+an=0
根据特征根的情况：
- 单实根 r：对应 erx
- k 重实根 r：对应 erx,xerx,…,xk−1erx
- 单复根 α±iβ：对应 eαxcosβx，eαxsinβx
- k 重复根：各乘以 1,x,…,xk−1
二阶 ODE 的完整公式：
特征方程 r2+ar+b=0：
- Δ>0：y=C1er1x+C2er2x
- Δ=0：y=(C1+C2x)erx
- Δ<0：y=eαx(C1cosβx+C2sinβx)，其中 r=α±iβ
适用于 f(x) 为多项式、指数函数、正余弦函数或其组合。
| f(x) 的形式 |
特解形式 yp |
| Pm(x) |
xkQm(x) |
| Pm(x)eαx |
xkQm(x)eαx |
| Pm(x)eαxcosβx 或 Pm(x)eαxsinβx |
xkeαx[Qn(x)cosβx+Rn(x)sinβx] |
其中 k 为 α（或 α+iβ）作为特征根的重数，Qm 为 m 次待定多项式，n=max(m,l)。
更通用的方法，不需要猜测特解形式。对于二阶方程 y′′+p(x)y′+q(x)y=f(x)：
设 y1,y2 为齐次方程的两个线性无关解，特解为：
yp=−y1∫Wy2fdx+y2∫Wy1fdx
其中 W=y1y2′−y1′y2 为朗斯基行列式。
形式：anxny(n)+an−1xn−1y(n−1)+⋯+a1xy′+a0y=f(x)
解法：令 x=et（x>0）化为常系数方程，或设 y=xr 代入齐次方程得指标方程。
二阶欧拉方程：ax2y′′+bxy′+cy=0，设 y=xr：
ar(r−1)+br+c=0
- 两实根 r1=r2：y=C1xr1+C2xr2
- 重根 r：y=(C1+C2lnx)xr
- 复根 r=α±iβ：y=xα[C1cos(βlnx)+C2sin(βlnx)]
- 勒让德方程：(1−x2)y′′−2xy′+n(n+1)y=0，解为勒让德多项式
- 贝塞尔方程：x2y′′+xy′+(x2−ν2)y=0，解为贝塞尔函数
这些是数学物理方程中常见的特殊函数方程，在偏微分方程的分离变量法中频繁出现。
形式：y′=A(x)y+f(x)
其中 y=(y1,…,yn)T，A(x) 为 n×n 矩阵。
y′=Ay，A 为常数矩阵。
解法：设 y=veλt，代入得特征值问题 Av=λv
- A 有 n 个线性无关特征向量时，通解为 y=∑civieλit
- 若特征向量不足，使用广义特征向量
矩阵指数法：
y(t)=eAty(0)
其中 eAt=∑k=0∞k!Aktk，可通过对角化、若尔当标准型或拉普拉斯变换计算。
对于二阶自治系统 x′=f(x)，平衡点 x0 满足 f(x0)=0。
在平衡点附近的线性化：设 J 为雅可比矩阵 ∂x∂f(x0)
根据 J 的特征值类型，平衡点分为：
| 特征值情况 |
平衡点类型 |
稳定性质 |
| 两负实根 |
稳定结点 |
渐近稳定 |
| 两正实根 |
不稳定结点 |
不稳定 |
| 一正一负实根 |
鞍点 |
不稳定 |
| 实部为负的共轭复根 |
稳定焦点 |
渐近稳定 |
| 实部为正的共轭复根 |
不稳定焦点 |
不稳定 |
| 纯虚根 |
中心点 |
稳定（Lyapunov稳定） |
对于一阶 IVP：y′=f(x,y)，y(x0)=y0
若 f 在矩形域 R=[x0−a,x0+a]×[y0−b,y0+b] 上连续且关于 y 满足利普希茨条件：
∣f(x,y1)−f(x,y2)∣≤L∣y1−y2∣
则存在 δ>0，使得 IVP 在区间 [x0−δ,x0+δ] 上有唯一解。
皮卡迭代（证明构造）：
y0(x)=y0
yn+1(x)=y0+∫x0xf(t,yn(t))dt
迭代序列在足够小的区间上一致收敛于方程的解。
若 f 在 R 上连续（但不一定满足利普希茨条件），则解存在（但不一定唯一）。
典型非唯一例子：
y′=3y2,y(0)=0
解有 y(x)≡0 和 y(x)=(x/3)3 等无穷多个解。
若 f 在整个 x 轴上连续且有界（或满足线性增长条件 ∣f(x,y)∣≤a(x)∣y∣+b(x)），则解可延拓至整个 x 轴。
L{f(t)}=F(s)=∫0∞e−stf(t)dt
拉普拉斯变换将微分方程转化为代数方程，特别适合处理具有间断或脉冲输入的 IVP。
| 原函数 f(t) |
变换 F(s) |
| 1 |
s1 |
| tn |
sn+1n! |
| eat |
s−a1 |
| sinat |
s2+a2a |
| cosat |
s2+a2s |
| eatsinbt |
(s−a)2+b2b |
| eatcosbt |
(s−a)2+b2s−a |
| δ(t−a) |
e−as |
| u(t−a) |
se−as |
线性性：L{af+bg}=aF(s)+bG(s)
微分性质：L{f′(t)}=sF(s)−f(0)
卷积性质：L{f∗g}=F(s)G(s)，其中 (f∗g)(t)=∫0tf(τ)g(t−τ)dτ
平移性质：
- 在 s 域：L{eatf(t)}=F(s−a)
- 在 t 域：L{f(t−a)u(t−a)}=e−asF(s)
- 对 ODE 两边取拉普拉斯变换，利用初值条件
- 解代数方程得到 Y(s)
- 对 Y(s) 进行拉普拉斯逆变换得到 y(t)（通常使用部分分式展开和查表）
例子：求解 y′′+3y′+2y=e−t，y(0)=1，y′(0)=0
L{y′′+3y′+2y}=L{e−t}
[s2Y−sy(0)−y′(0)]+3[sY−y(0)]+2Y=s+11
(s2+3s+2)Y−s−3=s+11
Y(s)=(s+1)(s+2)s+3+(s+1)2(s+2)1
通过部分分式展开和逆变换得到 y(t)。
若 p(x) 和 q(x) 在 x0 处解析，则 y′′+p(x)y′+q(x)y=0 的解可表示为：
y(x)=n=0∑∞an(x−x0)n
代入方程后比较系数可逐项确定 an。
若 x0 是正则奇点（(x−x0)p(x) 和 (x−x0)2q(x) 在 x0 处解析），设：
y(x)=(x−x0)rn=0∑∞an(x−x0)n
代入得指标方程确定 r，再接代求解系数。
贝塞尔函数的级数解：对 ν 阶贝塞尔方程，指标方程为 r2−ν2=0，解得 r=±ν。r=ν 对应的解为第一类贝塞尔函数 Jν(x)：
Jν(x)=k=0∑∞k!Γ(k+ν+1)(−1)k(2x)2k+ν
自治 ODE 系统 x′=f(x) 的右端不显含时间。其解 x(t) 在相空间中的轨迹称为轨线。
- 平衡点：f(x0)=0
- Lyapunov 稳定：对任意 ϵ>0，存在 δ>0，使得若初始偏差 <δ，则对所有 t≥0，偏差 <ϵ
- 渐近稳定：Lyapunov 稳定且 limt→∞∥x(t)−x0∥=0
- Lyapunov 第二方法：构造正定函数 V(x)（势函数），若 V′(x)=∇V⋅f(x)≤0，则平衡点稳定；若严格 <0，则渐近稳定
- 极限环：非线性系统中孤立的闭轨线，周围轨线螺旋趋近或远离
- 庞加莱-本迪克森定理：在二维相平面中，若有界区域内轨线始终不离开且无平衡点，则必存在极限环
- 范德波尔振子：x′′+μ(x2−1)x′+x=0，μ>0 时存在稳定极限环
- 鞍结分岔：平衡点产生与湮灭 x˙=r+x2
- 跨临界分岔：平衡点交换稳定性 x˙=rx−x2
- 叉形分岔：对称破缺 x˙=rx−x3
- 霍普夫分岔：平衡点失稳产生极限环
- 洛伦兹系统：x˙=σ(y−x),y˙=x(ρ−z)−y,z˙=xy−βz
- 对初值敏感依赖（蝴蝶效应）
- 奇异吸引子：相空间中分形维数的吸引集
yn+1=yn+hf(xn,yn)
局部截断误差 O(h2)，全局误差 O(h)。简单但精度低。
y~n+1=yn+hf(xn,yn)
yn+1=yn+2h[f(xn,yn)+f(xn+1,y~n+1)]
k1=f(xn,yn)
k2=f(xn+2h,yn+2hk1)
k3=f(xn+2h,yn+2hk2)
k4=f(xn+h,yn+hk3)
yn+1=yn+6h(k1+2k2+2k3+k4)
局部截断误差 O(h5)，工程中最常用的通用 ODE 求解器。
- 亚当斯-巴什福思法：显式，利用前几步的值外推
- 亚当斯-莫尔顿法：隐式，精度更高但需迭代求解
- 预估-校正法：显式预测 + 隐式校正的组合
某些 ODE 包含快变和慢变分量，快变分量迅速衰减但主导了时间步长的约束。
特点：
- 特征值实部为负且差距悬殊
- 显式方法需要极小的步长保持稳定
- 隐式方法（如隐式欧拉、后向微分公式 BDF）可允许大步长
吉尔（Gear）方法：面向刚性方程的 BDF 多步法族，1-6 阶可用。
- 打靶法：将 BVP 转化为含未知参数的 IVP，通过迭代调整初值满足边界条件
- 差分法：将微分算子差分化，转化为线性/非线性代数方程组
- 伽辽金法：用基函数逼近解，转化为弱形式的积分方程
- 牛顿第二定律：mdt2d2x=F(t,x,v) — 二阶 ODE 建模质点运动
- 谐振子：mx′′+kx=0 — 简谐运动
- 阻尼振子：mx′′+cx′+kx=0
- 受迫振动：mx′′+cx′+kx=F0cosωt — 共振现象
- RC 电路：RdtdQ+CQ=V(t)
- RLC 电路：Ldt2d2I+RdtdI+CI=dtdV
- 一阶 ODE 描述 RC/RL 瞬态响应，二阶 ODE 描述 RLC 振荡
- Logistic 增长：dtdP=rP(1−KP)
- 洛特卡-沃尔泰拉（捕食-猎物）模型：
- dtdx=ax−bxy
- dtdy=−cy+dxy
- 传染病模型（SIR 模型）：
- dtdS=−βSI
- dtdI=βSI−γI
- dtdR=γI
- 化学反应动力学：dtd[A]=−k[A]n（一级、二级反应速率方程）
- 传质扩散：菲克第二定律在一维稳态下的 ODE 形式
- 连续搅拌反应器（CSTR）：物料和能量平衡方程组
- 索洛增长模型：k˙=sf(k)−δk
- 新古典增长模型：人均资本积累的 ODE
- 存货动态模型：需求-供给调整过程的微分方程
dxd[p(x)dxdy]+q(x)y+λw(x)y=0
配上适当的边界条件构成 Sturm-Liouville 问题（S-L 问题）。
- 特征值均为实数
- 特征函数正交：∫abym(x)yn(x)w(x)dx=0，m=n
- 特征值可数且无聚点：λ1<λ2<⋯→∞
- 特征函数完备性：任意分段连续函数可按特征函数展开
| 方程 |
p(x) |
q(x) |
w(x) |
区间 |
特征函数 |
| y′′+λy=0 |
1 |
0 |
1 |
[0,L] |
sin(λnx) |
| 勒让德 |
1−x2 |
0 |
1 |
[−1,1] |
Pn(x) |
| 贝塞尔 |
x |
−ν2/x |
x |
[0,a] |
Jν(αnx) |
| 埃尔米特 |
e−x2 |
0 |
e−x2 |
(−∞,∞) |
Hn(x) |
- 偏微分方程的分离变量法（热传导、波动、拉普拉斯方程）
- 量子力学中的薛定谔方程（通解以 S-L 特征函数展开）
- 信号处理中的傅里叶级数（S-L 问题 y′′+λy=0 的特例）
- 正交多项式理论
本手册从基本概念出发，系统梳理了常微分方程的核心内容：
| 主题 |
核心要点 |
难度 |
| 一阶方程 |
分离变量、线性方程、恰当方程 |
⭐ |
| 高阶线性方程 |
常系数齐次/非齐次、待定系数法 |
⭐⭐ |
| 线性方程组 |
特征值方法、矩阵指数 |
⭐⭐⭐ |
| 存在唯一性 |
Picard-Lindelöf 定理 |
⭐⭐⭐ |
| 拉普拉斯变换 |
解 IVP、卷积、传递函数 |
⭐⭐ |
| 幂级数解法 |
常点、奇点、Frobenius 方法 |
⭐⭐⭐ |
| 定性理论 |
相平面、稳定性、分岔、混沌 |
⭐⭐⭐⭐ |
| 数值方法 |
RK4、多步法、刚性方程 |
⭐⭐⭐ |
| Sturm-Liouville |
特征函数展开、正交多项式 |
⭐⭐⭐⭐ |
进阶方向：
- 动力系统：从 ODE 定性理论延伸至全局分岔、正规形理论
- 微分代数方程（DAE）：混合微分和代数约束的系统
- 延迟微分方程（DDE）：导数依赖于过去时刻状态
- 随机微分方程（SDE）：含随机噪声项的微分方程
- 遍历理论：动力系统的长期统计性质
本文为数学知识库的一部分，更多内容请参阅：数学知识库索引
相关页面：偏微分方程 | 数值分析 | 实分析