微分方程初步：从导数到动力系统
本文基于模型知识整理（生成时未联网核对），解法与稳定性建议对照 Boyce & DiPrima 复核。
一句话定义
微分方程是"含未知函数导数的方程"——描述"变化规律"而非"状态本身"：$\frac{dy}{dt}=ky$（增长率正比当前量 → 指数解 $Ce^{kt}$）、$\frac{dy}{dt}=-ay+b$（趋向平衡 → 稳态）；解微分方程=从局部变化规则重建全局轨迹，解族靠初始条件定成员（一阶方程一个常数）。
为什么重要
一切连续演化的科学定律都是微分方程：牛顿第二定律、放射性衰减、传染病 SIR、神经 ODE——"写下方程"是建模，"解方程"是预测。它也是 kp-005（e 的微分方程身份）、kp-013（原函数=最简微分方程）、kp-028（梯度流=优化动力学）三线的自然交汇。
前置知识
kp-006/007（导数与链式）、kp-013（原函数）。
核心概念
- 阶与解族：最高阶导数的阶数=方程阶；n 阶方程的通解含 n 个自由常数——初始条件（x(0)=?、x'(0)=?）从中挑出特解。
- 可分离变量：$\frac{dy}{dx}=g(x)h(y)$ → $\int\frac{dy}{h(y)}=\int g(x)dx$——两边各自积分（kp-014 的直接应用），最基础的手算法。
- 一阶线性：$\frac{dy}{dx}+p(x)y=q(x)$ → 积分因子法（乘 $e^{\int p}$ 化为全微分）。
- 两个范式解：
- 指数增长/衰减 $\frac{dy}{dt}=ky$ ⇒ $y=Ce^{kt}$（k>0 爆炸增长：人口/复利；k<0 衰减：放射性/冷却）； - 趋向平衡 $\frac{dy}{dt}=k(A-y)$ ⇒ $y=A+Ce^{-kt}$——"变化率正比于差距"（牛顿冷却、电容充电、学习率 warmup）。
- 稳定性：平衡点 y* 处扰动后的去向——$\left.\frac{\partial f}{\partial y}\right|_{y^*}<0$ 扰动衰减（稳定）、>0 发散（不稳定）；线性化判据（非线性系统的局部稳定性）。
- 数值解预告：绝大多数微分方程无解析解——欧拉法（$y_{n+1}=y_n+h\cdot f(y_n)$，kp-017 的泰勒一阶思想）与 RK4 是机器的世界（现代神经网络时代的 Neural ODE 直接把 ODE 积分器做成层）。
原理与机制
为什么 $y'=ky$ 的解必是指数：kp-005 的等价说法——"导数正比自身"的函数正是 $Ce^{kt}$（分离变量或唯一性定理皆可证）；e 的微分方程身份在此兑现：自然增长的形状就是指数，因为自然的规律写成正比式。
平衡与稳定的判据从哪来：$\frac{dy}{dt}=f(y)$ 的平衡点 y*（f(y*)=0）附近线性化 $f(y)\approx f'(y^*)(y-y^*)$——扰动 $\delta=y-y^*$ 满足 $\delta'=f'(y^*)\delta$，解为 $\delta_0e^{f'(y^*)t}$：导数为负=扰动衰减=稳定。用一阶泰勒（kp-008）把非线性问题线性化，是整个动力系统理论的起手式。
为什么"积分常数"有物理身份：一阶方程通解的 C 由初始条件确定（t=0 时代入）；丢 C 丢掉的是整个解族中除一条外的所有轨迹——kp-013 警告的概率版实操。
图示
y' = k y → y = Ce^{kt} (正比自身 → 指数)
y' = k(A−y) → y = A + Ce^{−kt} (趋向平衡 → 稳态)
分离变量: dy/h(y) = g(x)dx → 两边积分
稳定性: y* 平衡, f′(y*)<0 稳定 / >0 不稳定 (线性化判据)
数值: 欧拉法 yₙ₊₁ = yₙ + h·f(yₙ) (无解析解时的机器解)直观类比
微分方程像"交规而非照片"：它不告诉你"现在在哪"（状态），而告诉你"每个位置该往哪走、走多快"（方向场）。解=按交规从起点走出的一条完整轨迹；初始条件=起点，积分常数=不同起点的家族。
实例或案例
- 放射性测年：$\frac{dN}{dt}=-\lambda N$ → $N=N_0e^{-\lambda t}$，半衰期 $t_{1/2}=\frac{\ln2}{\lambda}$——碳十四测年的数学核。
- Newton 冷却：$\frac{dT}{dt}=k(T_{env}-T)$ → 温度指数逼近室温——法医估死亡时间。
- 人口逻辑斯蒂：$\frac{dP}{dt}=rP(1-\frac PK)$——指数增长加环境上限（K），S 型曲线；分层 kp-025 的平衡点稳定性分析素材。
常见误区
- 误区一："解微分方程=找公式"。绝大多数（尤其高维非线性）无解析解；数值积分（欧拉/RK4）与定性分析（稳定性）才是主战场——"解不出"不等于"学不了"。
- 误区二："分离变量时丢解"。$\frac{dy}{dx}=y^2$ 除以 $y^2$ 时 y=0 这个（恒定）解被除没了——除法前检查"除数可为零"的常值解。
- 误区三："平衡点都是稳定的"。$\frac{dy}{dt}=k(y-A)$ 的平衡 A 处 f'=k>0——不稳定（偏离即加速远离）；稳定性要算导数符号不是数平衡点。
与其他知识点的关系
自测题
- 解 $\frac{dy}{dt}=3y$，y(0)=5。
答：$y=5e^{3t}$（指数增长；C 由初始条件定）。
- 解牛顿冷却 $\frac{dT}{dt}=-k(T-20)$，T(0)=90，k=0.1。
答：$T=20+70e^{-0.1t}$（趋向 20 的稳态解）。
- 平衡点 $y^*=2$（$f(y)=(y-2)(y-3)$ 的两个平衡）稳定吗？
答：f'(y)=2y−5，f'(2)=−1<0 稳定；y*=3 处 f'=1>0 不稳定——两个平衡一稳一不稳。
延伸阅读
- Boyce & DiPrima《Elementary Differential Equations》（标准教材）。
- 3Blue1Brown《微分方程系列》（方向场与欧拉法动画）。
- Strogatz《Nonlinear Dynamics and Chaos》（动力系统直觉神书）。