可分离方程
\(dy/dx=g(x)h(y)\)
MIT 18.01 · 第 16 课
一个冷却中的物体、一棵倒下的树、一盏旋转的探照灯——变化率往往依赖于当前状态。微分方程描述这种关系,而分离变量法 (separation of variables) 是求解最简单一类微分方程的入门工具。
第 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 的函数",就可用分离变量法。
除以 \(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}).\]
这是放射性衰变、连续复利、牛顿冷却定律的统一数学模型。
物体温度 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\) 求出"煮熟"所需时间。
\[\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.\]
\[\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) 可能漏掉常数解