第 3 章 微分方程与动力系统：随时间演化的数学
第 3 章 · 微分方程与动力系统：随时间演化的数学
章节摘要：凡是"下一刻的状态由当前状态决定"的系统，都是微分方程的地盘。本章 3.1 讲常微分方程的建模套路与数值解法（欧拉法到龙格-库塔），3.2 用相平面与李雅普诺夫函数判断平衡点的稳定性，3.3 沿着参数变化看分岔与混沌，并顺路接上控制理论的入口。读完本章，你能对一个演化系统回答三个问题：轨迹怎么算、平衡稳不稳、参数变了会发生什么。
学习目标
阅读完本章，你应当能够：
- 把一个连续变化过程写成状态方程组并选择合适求解器；
- 在相平面上画轨线、线性化平衡点并用特征值判稳定性；
- 解释倍周期分岔通向混沌的路线，计算李雅普诺夫指数的含义；
- 说出闭环反馈改变系统行为的基本机制。
核心概念速览
微分方程的解不是一条公式，而是一族轨迹；建模者真正要回答的往往不是"轨迹长什么样"，而是"它最终去哪、去得稳不稳"。
3.1 ODE 建模与数值解
从化学反应釜降温工单出发：状态方程怎么列、欧拉法为什么精度不够、龙格-库塔与自适应步长怎么用。
3.2 相平面与稳定性判据
捕食者-猎物系统的相平面解读、平衡点线性化、李雅普诺夫直接法的直觉与一个构造实例。
3.3 分岔、混沌与控制入门
以虫口模型演示倍周期分岔，讨论混沌的初值敏感性，最后用倒立摆说明反馈如何改写稳定性。
子章节之间的逻辑关系
3.1 会算轨迹 ──► 3.2 会判稳定 ──► 3.3 会看参数与控制 │ │ └── 数值解是后两节的实验平台 ──┘
3.1 提供计算能力，3.2 提供定性判断力，3.3 把视野从"固定系统"扩展到"变参数系统"与"可干预系统"。三节合起来构成动力系统方法的完整闭环。
门诊记录：本章算例一览
| 工单 | 用到的方法 | 关键产出 |
|---|---|---|
| 反应釜降温 | 状态方程、欧拉法、RK45 | 误差阶对比与刚性判别 |
| 渔业种群管理 | 相平面、闭合轨道 | 振荡管理含义 |
| 阻尼摆稳定证明 | 李雅普诺夫能量函数 | 单调性数值验证 |
| 虫口种群预测 | 分岔图、李雅普诺夫指数 | 长期预测失效的机理 |
| 倒立摆稳定 | 状态反馈配置特征值 | 开环闭环特征值对照 |
本章方法在 7.2 节以 SIR 模型的形态正面出场；步长与稳定性的权衡则与 5.1 节的数值误差分析互为表里。
本章的实践建议分两层。基础层：把 3.1 的求解器调用练熟，能对任意一阶方程组给出数值解与误差量级判断，这是动力系统分析的入场券。进阶层：用 3.2 的相平面工具分析你所在行业的一个振荡问题（库存、舆情、心率皆可），再用 3.3 的分岔视角问一句"哪个参数是行为的开关"。分岔图建议亲手扫一遍参数轴，那张图比任何文字都更能传达"渐变引发突变"的冲击。
出科小考
自测三题：其一，RK45 跑同一方程步数异常多，你的第一个怀疑与对策是什么？其二，线性化后特征值是一对纯虚数，为什么结论要打引号？其三，为什么说混沌面前"第 50 步的值"是错误交付物？（答案分别在 3.1、3.2、3.3。）
若三题都能脱口而出，说明动力系统的"算轨迹、判稳定、看参数"三重能力已经立起来了。
前置知识与后续延伸
- 前置知识：导数的物理含义、第 2 章的特征值概念；
- 后续延伸：第 7 章的流行病模型是本章方法在生物系统的直接应用；第 5 章的数值稳定性概念与本章的步长选择互为镜像。