2. 从实际问题到微分方程

建立微分方程模型，本质上不是“套公式”，而是把一句关于现实世界的自然语言，翻译成一句关于“变化”的数学语言。翻译过程通常围绕三个抓手展开：变化率（某量随时间的瞬时变化速度）、微元法（在极短时间或极小空间上做局部平衡）、守恒律（总量变化 = 流入 − 流出 + 内部生成 − 内部消耗）。掌握这三者，人口增长、放射性衰变、牛顿冷却这些经典案例就不再是三个孤立公式，而是同一套思想的三种表现。

一、变化率：把“速度”写成导数

建模的第一步是明确“谁是未知函数、谁在变、对谁变”。设未知量为 y(t)，则它在 t 时刻的瞬时变化率就是 y′(t) = dy/dt。关键在于：题目里“增长快慢”“冷却快慢”“衰变快慢”这类描述，往往不是直接给出导数，而是给出导数与当前状态的依赖关系。

以人口增长为例。若假设“单位时间内新增人数与当前人口成正比”，比例系数为 r（内禀增长率），则

dP/dt = r·P,   P(0) = P₀

这里 r 的量纲是“1/时间”，P 的量纲是“人”，乘积 rP 的量纲是“人/时间”，与 dP/dt 一致——量纲一致性检查是建模后必做的第一道校验。该方程的解为 P(t) = P₀e^rt，即指数增长。但现实资源有限，于是把常数 r 改为随人口下降的 r(1 − P/K)，得到逻辑斯谛方程 dP/dt = rP(1 − P/K)，其中 K 为环境容纳量。这一步改动体现的正是“变化率依赖关系”的建模自由度。

二、微元法：在 [t, t+Δt] 上做局部平衡

微元法是守恒律落地为方程的具体操作。取一个极短时间区间 [t, t+Δt]，考察未知量的增量 Δy = y(t+Δt) − y(t)，把它拆成“流入、流出、生成、消耗”四项之和，再两边除以 Δt 并令 Δt→0，就得到微分方程。这个“先写差分、再取极限”的路径，比直接写导数更不容易漏项。

以放射性衰变为例。设 N(t) 为未衰变核数。在 Δt 内衰变的核数近似正比于当前核数与时间长度：ΔN ≈ −λNΔt，其中 λ 为衰变常数。除以 Δt 取极限得

dN/dt = −λN,   N(0) = N₀  ⟹  N(t) = N₀e^(−λt)

负号表示“减少”，这是符号约定的核心：先定义正方向（此处“核数增加为正”），再让负号自然出现，而不是事后凭感觉补。半衰期 T_1/2 = ln2/λ ≈ 0.693/λ，可作为参数换算的检验点：若 λ = 0.0231 /年，则半衰期约 30 年，与碳-14 的常用数据量级吻合，说明参数取值合理。

三、守恒律：牛顿冷却的“热流平衡”

牛顿冷却定律常被直接写成 dT/dt = −k(T − T_a)，但用微元法推一遍更能看清物理含义：物体在 Δt 内散失的热量正比于它与环境的温差和时长，即 ΔQ ≈ −h(T − T_a)Δt；又 ΔQ = mcΔT，联立得 mcΔT ≈ −h(T − T_a)Δt，取极限即得上述方程，其中 k = h/(mc)。

求解时令 u = T − T_a，方程化为 du/dt = −ku，得 T(t) = T_a + (T₀ − T_a)e^−kt。举例：一杯 90 ℃ 的水放在 20 ℃ 房间，测得 10 分钟后为 60 ℃，则 40 = 70e^−10k，解得 k ≈ 0.056 /min。据此可预测 20 分钟后温度约 20 + 70e^−1.12 ≈ 42.8 ℃。这类“用一组观测反推参数、再预测另一时刻”的流程，是参数辨识的最小闭环，也是后续数值求解章节的常见起点。

四、建模流程中的几个易错点

第一，未知量定义不清。人口模型里 P 是“总人数”还是“人口密度”？两者方程形式可能相同，但参数含义与量纲不同，后续代入数据会出错。第二，正方向与符号约定不统一。衰变中若把“减少”定义为正，则方程变为 dN/dt = λN，解为指数增长，物理意义完全相反。第三，忽略初始条件。微分方程给出的是“一族解”，只有配上 y(0) 才确定唯一解；建模时必须同时交代初始状态。第四，参数未做合理性检验。用半衰期反算 λ、用稳态值检验 K、用观测点反推 k，都是低成本的有效校验。

五、从解析解到数值求解的衔接

上述三个案例都能求出解析解，但它们的价值不止于“会解”。一方面，解析解给出了定性结论：指数增长无上界、衰变趋于零、冷却趋于环境温度，这些是判断数值结果是否合理的基准。另一方面，当模型加入变系数、非线性项或外部输入（如 r(t) 随季节变化、冷却中环境温度 T_a(t) 变化）时，解析解往往不存在，必须转向数值方法。此时，本节建立的“变化率 + 微元 + 守恒”框架依然是数值格式设计的出发点：欧拉法本质上就是把 dy/dt ≈ (y_n+1 − y_n)/Δt 代回方程，正是微元平衡的离散版本。

因此，学习本节的目标不是记住三个公式，而是形成一种可迁移的建模习惯：先问“什么在变、对谁变、变化率由什么决定”，再用微元法把守恒关系写成差分式，取极限得到方程，最后用初始条件、量纲和参数合理性做三重校验。这套习惯将直接支撑后续的数值求解与模型分析。