建模、初值问题与解的几何图像

很多微分方程题目看起来从公式开始，其实真正的起点通常是一句现实语言：水温怎样下降，人口怎样增长，药物浓度怎样变化，电路电流怎样响应。把这句话写成方程以后，我们还没有得到一个完整预测。方程只告诉我们“在每个状态下怎样变”，初始条件告诉我们“现在在哪里”，图像则帮助我们检查这种变化是否说得通。

本章先不急着求解析公式。我们的目标是看懂一阶初值问题的结构，并学会用方向场判断解曲线的大致行为。后面几章会训练具体求解方法；这里先建立一个更稳的眼光：即使没有公式，也能读出增长、衰减、平衡和长期趋势。

---

一条方程还不是一个预测

一阶常微分方程常写成

$$\frac{dy}{dt}=f(t,y)$$

这里 $t$ 是自变量，通常可以理解为时间；$y(t)$ 是随时间变化的量；右边 $f(t,y)$ 给出当前时刻和当前状态下的变化率。

这个式子本身描述的是一整族可能的运动。比如温度模型

$$\frac{dT}{dt}=-0.08(T-22)$$

说的是温度 $T$ 和室温 $22$ 的差越大，降温或升温越快。可是咖啡一开始是 $90^\circ\mathrm{C}$，还是刚从冰箱拿出来只有 $5^\circ\mathrm{C}$，会得到完全不同的温度曲线。同一个变化规律，需要配上一个起点，才会确定一条具体曲线。

在实际建模中，我们常把一个初值问题写成

$$\frac{dy}{dt}=f(t,y), \qquad y(t_0)=y_0$$

它读作：函数 $y(t)$ 的变化率由 $f(t,y)$ 决定，并且在 $t=t_0$ 时，状态等于 $y_0$。

---

变量、单位与变化率假设

建模最容易出错的地方不是微积分计算，而是变量还没说清楚就开始写方程。一个简洁的建模过程至少要回答四个问题。

以牛顿冷却定律为例，如果环境温度是 $T_s$，物体温度是 $T(t)$，常见模型是

$$\frac{dT}{dt}=-k(T-T_s), \qquad k>0$$

负号不是装饰。若 $T>T_s$，右边为负，温度下降；若 $T<T_s$，右边为正，温度上升；若 $T=T_s$，右边为零，温度保持不变。一个符号就包含了图像趋势。

---

初值问题：把变化规律钉在一个起点上

考虑方程

$$y'=x-y$$

它没有告诉我们 $y(0)$ 是多少。因此它对应许多条解曲线。若补上 $y(0)=1$，就得到一个初值问题：

$$y'=x-y, \qquad y(0)=1$$

这句话的几何意思是：在整张平面上，每一点的斜率都按 $x-y$ 计算；我们要找的那条曲线必须经过点 $(0,1)$，并且在经过每个点时都服从当地斜率。

初值不一定总在 $t=0$ 给出。比如 $y(2)=5$ 也完全可以，它表示我们在 $t=2$ 这个时刻知道状态。数学上，初值点是 $(t_0,y_0)$；现实中，它就是开始预测时手里已有的一次测量。

初值问题的关键词是“局部”。一阶方程先告诉你当前点处的斜率，然后你沿着斜率前进一点，再用新位置重新计算斜率。这个想法后来会发展成 Euler 方法和 Runge-Kutta 方法；在本章，它先以方向场的形式出现。

---

解曲线与积分曲线

如果函数 $y=\phi(t)$ 满足微分方程

$$\phi'(t)=f(t,\phi(t))$$

那么它的图像称为这条微分方程的一条解曲线，也常叫积分曲线。名字里的“积分”提醒我们：从变化率恢复状态，本质上是在累积变化。

对初值问题来说，解曲线还要经过初值点：

$$\phi(t_0)=y_0$$

解曲线不是随手画出来的光滑曲线。它在每个点的切线斜率都必须等于微分方程右边给出的数值。比如方程 $y'=x-y$ 中，点 $(2,1)$ 处的斜率是 $1$；若一条候选曲线经过 $(2,1)$，它在那里就必须向右上方倾斜，斜率为 。

一个有用的自检是把“曲线、导数、方程”放在一起看。给定候选函数 $\phi(t)$，先求 $\phi'(t)$，再把 $\phi(t)$ 代入方程右边。如果两边一致，并且初值也满足，它才是这个初值问题的解。

---

方向场：先看斜率，再看公式

方向场也叫斜率场。对一阶方程

$$y'=f(x,y)$$

我们在平面上的许多点 $(x,y)$ 计算 $f(x,y)$，然后画一小段斜率为 $f(x,y)$ 的短线。短线不是解本身，而是给解曲线准备的“局部路标”。

看方向场时，可以按三个层次读图。

下面的交互可以直接移动初值点，观察同一张方向场中不同积分曲线如何改变。

---

从方向场读长期趋势

长期趋势通常从四类现象开始看：增长、衰减、平衡、靠近或远离某个平衡。对自治方程

$$y'=g(y)$$

右边不显含时间，斜率只由当前状态 $y$ 决定。若存在常数 $c$ 使 $g(c)=0$，那么

$$y(t)=c$$

就是一条平衡解。图像上它是一条水平线；方向场上沿着这条线的短线都是水平的。

以逻辑斯蒂型方程为例：

$$y'=y\left(1-\frac{y}{4}\right)$$

它有两个平衡解：$y=0$ 和 $y=4$。当 $0<y<4$ 时，$y'>0$，解向上增长；当 时，，解向下回落。于是从正初值出发的解会趋向 ，这条水平线可以解释为长期容量或稳定水平。

这里有一个细节：如果 $y$ 表示人口，负初值在现实中没有意义；即使方程在数学上可以计算，我们也不应把负人口的曲线当成现实预测。几何图像能提醒我们模型行为，但不能替代对变量意义的检查。

---

例题：冷却模型的图像解释

一杯热咖啡放在 $22^\circ\mathrm{C}$ 的房间中。设 $T(t)$ 表示 $t$ 分钟后的温度，假设冷却速度与温差成正比，并取 $k=0.08$。若初始温度为 $90^\circ\mathrm{C}$，模型为

$$\frac{dT}{dt}=-0.08(T-22), \qquad T(0)=90$$

我们先不求公式，只用方程读图像。

这个例子展示了本章的核心思路：不求解析解，也能判断曲线的方向、凹平趋势和长期水平。求解公式当然有用，但它应该和图像解释互相校验。

---

解的合理性检查

解出公式以后，很多错误并不会自动暴露。一个表达式可能看起来很整齐，却在单位、极端情况或图像趋势上出问题。下面三种检查特别适合一阶 ODE 建模。

单位检查

方程两边单位必须一致。若 $T$ 用摄氏度，$t$ 用分钟，则 $dT/dt$ 的单位是摄氏度每分钟。在冷却模型

$$\frac{dT}{dt}=-k(T-T_s)$$

中，$T-T_s$ 的单位是摄氏度，所以 $k$ 的单位必须是每分钟。若把 $k$ 当成一个没有单位的普通数字，模型在换时间单位时就会混乱。

极端情况

把参数取到很简单的值，观察模型是否符合常识。冷却模型中若 $k=0$，方程变成

$$\frac{dT}{dt}=0$$

这表示物体温度不变。它对应“没有热交换”的极端情况，是合理的。若 $k$ 很大，温度应当很快靠近环境温度，也符合直觉。

图像趋势

看导数符号和斜率大小是否匹配现实。热咖啡高于室温时应下降，低于室温时应升温，等于室温时应保持。若一个模型让热咖啡在普通房间里越放越热，除非题目中有加热源，否则符号就错了。

---

常见误区

把通解当成初值问题的答案

方程 $y'=y$ 的通解是

$$y=Ce^t$$

但初值问题 $y'=y,\ y(0)=3$ 的答案必须把 $C$ 定下来。代入初值可得 $C=3$，所以具体解是 $$。如果答案里还留着任意常数，就没有完成初值问题。

把方向场里的短线当成速度箭头

在 $x$-$y$ 平面中，方向场短线的作用是表示解曲线的切线斜率。它不是粒子在平面中以箭头方向移动的速度向量。对一阶方程 $y'=f(x,y)$，我们通常把 $x$ 当作自变量向前推进，短线告诉我们 $y$ 随 的升降趋势。

忽略变量的现实范围

如果 $y$ 表示人口、浓度、质量，它通常不能为负。数学方程可能在负数区域仍有方向场，但那部分图像未必有现实意义。建模时要把变量允许范围和方程一起记录下来。

只相信公式，不看图像

公式可能因为积分常数、符号、代数化简出错而偏离模型。方向场和简单趋势检查可以快速发现异常：该增长的地方是否增长？该趋近的平衡是否被穿过？初值点是否真的在曲线上？

---

小结

本章建立了三件事之间的联系。微分方程给出局部变化率，初始条件选定起点，解曲线是在每一点都服从斜率规则的图像。方向场把这套规则可视化，让我们在还没有解析公式时也能判断解的大致形状。

接下来学习可分离变量方程时，你会开始真正求出一些显式解。请保留本章的检查习惯：先看变量和单位，再看导数符号，最后把公式放回图像中验证。

---

练习

练习一：设一阶方程为

$$y'=2-y$$

不求通解，判断 $y=2$ 是否为平衡解；若 $y(0)=5$，解曲线一开始上升还是下降？长期可能靠近哪里？

练习二：某药物在血液中的含量 $A(t)$ 以与当前含量成正比的速度减少。若时间单位是小时，比例常数为 $0.3$，初始含量是 $40$ 毫克，写出初值问题，并说明 $0.3$ 的单位。

练习三：考虑冷却模型

$$\frac{dT}{dt}=-k(T-20), \qquad k>0$$

若初值为 $T(0)=10$，图像一开始应上升还是下降？它能否表示“从冰箱拿出的饮料在室温下回温”？

练习四：方程

$$y'=x-y$$

在点 $(1,3)$ 处的方向场短线斜率是多少？若一条解曲线经过这个点，它在那里应向上还是向下？