MIT 18.01 · 第 16 课

微分方程:分离变量法

一个冷却中的物体、一棵倒下的树、一盏旋转的探照灯——变化率往往依赖于当前状态。微分方程描述这种关系,而分离变量法 (separation of variables) 是求解最简单一类微分方程的入门工具。

授课:Prof. David Jerison 日期:2006-10-13 难度:中等

这一课在课程中的位置

第 15 课引入了反导数 \(\int f(x)\,dx\)。这一课把反导数用于解未知函数是 x 的函数的方程——微分方程。分离变量法把形如 \(dy/dx=g(x)h(y)\) 的方程拆成两边各自积分,是微积分"互逆性"最直接的体现。


什么是微分方程?

微分方程(differential equation) 是包含未知函数及其导数的方程。我们关心的是:找到满足该方程的函数 \(y(x)\)。

最简单的类型是一阶方程:

\[\frac{dy}{dx}=f(x,y).\]

可分离变量方程

若能写成

分离变量形式

\[\frac{dy}{dx}=g(x)\,h(y),\]

即右边能拆成"只含 x 的函数"乘"只含 y 的函数",就可用分离变量法。

求解步骤

  1. 把方程写成 \(\dfrac{dy}{dx}=g(x)h(y)\);
  2. 分离:把含 y 的部分移到左边(连同 dy),含 x 的部分移到右边(连同 dx):\(\dfrac{dy}{h(y)}=g(x)\,dx\)(注意 \(h(y)\neq 0\));
  3. 两边积分:\(\displaystyle\int \frac{dy}{h(y)} = \int g(x)\,dx\);
  4. 解出 y(如果可能),并加上常数 C;
  5. 用初始条件(如 \(y(0)=1\))确定常数 C。

除以 \(h(y)\) 的前提:分离时两边除以 \(h(y)\),要检查 \(h(y)=0\) 的常数解是否被遗漏。例如 \(dy/dx=y\) 除了 \(y=Ce^x\) 之外,\(y=0\) 也是一个解。


例题:指数增长/衰减

\[\frac{dy}{dx}=ky\quad (k\text{ 为常数}).\]

分离:\(\dfrac{dy}{y}=k\,dx\)。两边积分:

\[\ln|y|=kx+C_1\;\Longrightarrow\; y=Ce^{kx}\quad(C=\pm e^{C_1}).\]

这是放射性衰变、连续复利、牛顿冷却定律的统一数学模型。

例题:牛顿冷却定律(复习题 11)

物体温度 T 的变化率与环境温差成正比:

\[\frac{dT}{dt}=-k(T-T_e),\quad T_e=100^\circ\text{(沸水)}.\]

设 \(u=T-T_e\),则 \(du/dt=-ku\),解得 \(u=Ce^{-kt}\),即 \(T(t)=100+Ce^{-kt}\)。

由室温 \(T(0)=20\) 得 \(C=-80\)。再由 \(T(5)=30\) 定出 k,最后解 \(T(t)=60\) 求出"煮熟"所需时间。


例题:\(dy/dx = xy + x\)(复习题 10)

\[\frac{dy}{dx}=x(y+1).\]

分离:\(\dfrac{dy}{y+1}=x\,dx\)。积分:

\[\ln|y+1|=\frac{x^2}{2}+C_1\;\Longrightarrow\; y+1=Ce^{x^2/2}.\]

由 \(y(0)=1\) 得 \(1+1=C\Rightarrow C=2\),故

\[y(x)=2e^{x^2/2}-1.\]

例题:\(dx/dt = 2t\sqrt{1-x^2}\)(复习题 5b)

\[\frac{dx}{\sqrt{1-x^2}}=2t\,dt.\]

积分:\(\arcsin x = t^2 + C\)。由 \(x(0)=1\) 得 \(\arcsin 1=\pi/2=C\),故

\[x(t)=\sin\!\left(t^2+\frac{\pi}{2}\right)=\cos(t^2).\]


常见误区与边界情况


知识检查

问题 1:下列哪个方程可以用分离变量法?

问题 2:分离变量法的第一步是什么?

问题 3:方程 \(\dfrac{dy}{dx}=ky\) 的通解是?

问题 4:分离变量时两边除以 \(h(y)\),需要注意什么?


本课要点

可分离方程

\(dy/dx=g(x)h(y)\)

分离变量

含 y 的与 dy 同侧,含 x 的与 dx 同侧

两边积分

各自求反导数,留常数 C

初始条件

确定常数 C,得到特解

检查特解

除以 h(y) 可能漏掉常数解


参考资料