学径XUEJING · 个人学习平台通知设置

课程库 Calculus 2 S5 · 一阶常微分方程 › 第 3

一阶线性方程 —— 积分因子的魔法

普通课45 分钟

本课导读

可分离变量法专治"乘法脸",但很多方程偏偏长着"加法脸":y+y=exy' + y = e^xxy+y=3x2xy' + y = 3x^2……拆不成乘积,分离不了。本课的主角积分因子是一阶线性方程的万能钥匙:往方程两边乘一个精心挑选的函数,左边就会奇迹般地折叠成一个"完美导数",然后积分收工。课末还有一场认亲大会:第 1 课练习 2 里那个来历不明的解 y=x2+C/xy = x^2 + C/x,今天查明身世。预计 45 分钟。

说明

前置工具箱:乘积求导法则 (uv)=uv+uv(uv)' = u'v + uv'(本课的灵魂,倒着用);epdxe^{\int p\,dx} 型的积分:dxx=lnx+C\displaystyle\int \frac{dx}{x} = \ln|x| + Celnx=xe^{\ln x} = x;分部积分偶尔客串。

🔧 30 秒回顾 · 本课要用的旧手艺(熟练可跳过)

① 乘积法则 (uv)=uv+uv(uv)' = u'v + uv':本课的灵魂,而且是倒着用。微例:(e2xy)=e2xy+2e2xy(e^{2x}y)' = e^{2x}y' + 2e^{2x}y——例题 1 的"折叠"就是把这一行从右往左读。⚠️ 翻车点:漏项,尤其是 μy\mu'y 那一项。

eeln\ln 互相抵消elnx=xe^{\ln x} = x,推广 eklnx=xke^{k\ln x} = x^k。微例:μ=e2lnx=x2\mu = e^{2\ln x} = x^2——练习 3 原式。⚠️ 翻车点:把 e2lnxe^{2\ln x} 误算成 2x2x

③ 凑微分:被积式里"外面挂着内层的导数(差个常数倍)"就凑。微例:xex2dx=12ex2d(x2)=12ex2+C\displaystyle\int xe^{x^2}dx = \frac{1}{2}\int e^{x^2}\,d(x^2) = \frac{1}{2}e^{x^2} + C——练习 2 原式。⚠️ 翻车点:漏系数 12\frac{1}{2}

一阶线性方程长什么样

一阶线性方程的标准形式:

y+p(x)y=q(x)y' + p(x)\,y = q(x)

"线性"指未知数 yyyy' 都只以一次幂、不相乘地出现(没有 y2y^2yyyy'siny\sin y 这些花样)。p(x)p(x)q(x)q(x) 可以随便复杂——它们不含 yy,只是已知的系数和右端项。

认脸练习:y+2y=4y' + 2y = 4 ✓;xy+y=3x2xy' + y = 3x^2 ✓(两边除以 xx 就站好队形:y+1xy=3xy' + \frac{1}{x}y = 3x);y=y2yy' = y^2 - y ✗(y2y^2 出局——不过下一课它有别的路)。

核心思想:让左边折叠成完美导数

观察标准形式的左边 y+p(x)yy' + p(x)y,再看乘积求导法则:

(μ(x)y)=μy+μy\big(\mu(x)\, y\big)' = \mu y' + \mu' y

两个式子长得像不像?如果给方程两边同乘一个函数 μ(x)\mu(x),左边变成 μy+μpy\mu y' + \mu p\, y——只要 μ\mu 满足 μ=μp(x)\mu' = \mu\, p(x),它就和 (μy)(\mu y)' 一模一样。而 μ=μp\mu' = \mu p 是个可分离方程(上一课的手艺):dμμ=pdx\int \frac{d\mu}{\mu} = \int p\,dx,解得

μ(x)=ep(x)dx\mu(x) = e^{\int p(x)\,dx}

这个 μ\mu 就叫积分因子(integrating factor)。于是整个解法四步:

  1. 整理成标准形式,认出 p(x)p(x)
  2. 算积分因子 μ=epdx\mu = e^{\int p\,dx}(这一步不用带 +C+C,取一个就够)
  3. 两边同乘 μ\mu,左边必定折叠为 (μy)(\mu y)'——如果折不出来,说明前面算错了,这是免费的自检点
  4. 两边积分,解出 yy(这时才出现 CC

提示

类比记忆:积分因子像拉链头——y+pyy' + py 是两排张开的链齿,看着没法收拾;乘上 μ\mu 这个拉链头一拉,"唰"地合成一条完整的 (μy)(\mu y)',接下来积分就是顺着拉链滑到底。

例题 1 · 第一次拉拉链

解方程 y+2y=4y' + 2y = 4

第 1-2 步 · 认出 p,算出 μ

已是标准形式,p(x)=2p(x) = 2q(x)=4q(x) = 4。积分因子:

μ=e2dx=e2x\mu = e^{\int 2\,dx} = e^{2x}
第 3 步 · 两边同乘,见证折叠e2xy+2e2xy=4e2xe^{2x}y' + 2e^{2x}y = 4e^{2x}

左边正是 (e2xy)(e^{2x}y)'——自检:把 (e2xy)(e^{2x}y)' 用乘积法则展开,e2xy+2e2xye^{2x}y' + 2e^{2x}y,分毫不差。于是

(e2xy)=4e2x\big(e^{2x}\,y\big)' = 4e^{2x}
第 4 步 · 积分收工

两边积分:e2xy=4e2xdx=2e2x+Ce^{2x}y = \displaystyle\int 4e^{2x}dx = 2e^{2x} + C。除以 e2xe^{2x}

y=2+Ce2xy = 2 + Ce^{-2x}

读答案:tt \to \infty 时所有解都奔向 y=2y = 2——又一个"平衡值"(y2y \equiv 2 代回方程:0+4=40 + 4 = 4 ✓)。指数项 Ce2xCe^{-2x} 是初值留下的"记忆",随时间被遗忘。

例题 2 · 认亲大会:y = x² + C/x 的身世

xy+y=3x2xy' + y = 3x^2x>0x > 0)的通解。——第 1 课练习 2 让你"验证"过它的解,现在正面推导。

第 1-2 步 · 标准化,算 μ

两边除以 xxy+1xy=3xy' + \dfrac{1}{x}y = 3x,故 p(x)=1xp(x) = \dfrac{1}{x}

μ=e1xdx=elnx=x\mu = e^{\int \frac{1}{x}dx} = e^{\ln x} = x

积分因子就是 xx 本身(🔧 工具箱 ②:elnx=xe^{\ln x} = x)!

第 3-4 步 · 折叠、积分

两边乘 xx——方程变回 xy+y=3x2xy' + y = 3x^2(原来出题人早就乘好了!),左边折叠:

(xy)=3x2xy=x3+Cy=x2+Cx(x\,y)' = 3x^2 \quad\Rightarrow\quad xy = x^3 + C \quad\Rightarrow\quad y = x^2 + \frac{C}{x}

身世大白:第 1 课那个"来历不明"的解,是积分因子两步推出来的。你现在拥有的不只是验证能力,而是生产能力

顺手再看一位老朋友:咖啡方程 T=k(T20)T' = -k(T - 20) 移项后是 T+kT=20kT' + kT = 20k——它也是线性方程!用 μ=ekt\mu = e^{kt} 走一遍流程,得到的和上一课分离变量的结果一字不差:T=20+CektT = 20 + Ce^{-kt}。一题两解不是浪费,是互相验算:两条路通到同一个终点,答案几乎不可能是错的。

✍️ 练习

练习 1 · 完整流程 + 初值

解初值问题 y+3y=6y' + 3y = 6y(0)=1y(0) = 1

✅ 参考答案

μ=e3x\mu = e^{3x}(e3xy)=6e3x(e^{3x}y)' = 6e^{3x}e3xy=2e3x+Ce^{3x}y = 2e^{3x} + C;通解 y=2+Ce3xy = 2 + Ce^{-3x}。代入 (0,1)(0,1)1=2+C1 = 2 + CC=1C = -1

y=2e3xy = 2 - e^{-3x}

验证:y=3e3xy' = 3e^{-3x}y+3y=3e3x+63e3x=6y' + 3y = 3e^{-3x} + 6 - 3e^{-3x} = 6 ✓,y(0)=21=1y(0) = 2 - 1 = 1 ✓。解从 1 出发单调爬向平衡值 2。

练习 2 · 变系数热身

y+2xy=xy' + 2xy = x 的通解。(提示:2xdx=x2\int 2x\,dx = x^2;右边积分请凑微分。)

✅ 参考答案

μ=ex2\mu = e^{x^2}(ex2y)=xex2(e^{x^2}y)' = xe^{x^2};右边凑微分:xex2dx=12ex2+C\displaystyle\int xe^{x^2}dx = \frac{1}{2}e^{x^2} + C。故 ex2y=12ex2+Ce^{x^2}y = \frac{1}{2}e^{x^2} + C

y=12+Cex2y = \frac{1}{2} + Ce^{-x^2}

验证:y=2xCex2y' = -2xCe^{-x^2}y+2xy=2xCex2+x+2xCex2=xy' + 2xy = -2xCe^{-x^2} + x + 2xCe^{-x^2} = x ✓。

练习 3 · 像例题 2 那样先站队

xy+2y=x2xy' + 2y = x^2x>0x > 0)的通解。

✅ 参考答案

标准化:y+2xy=xy' + \dfrac{2}{x}y = xp=2xp = \dfrac{2}{x}μ=e2lnx=x2\mu = e^{2\ln x} = x^2。两边乘 x2x^2

(x2y)=x3x2y=x44+Cy=x24+Cx2(x^2 y)' = x^3 \quad\Rightarrow\quad x^2 y = \frac{x^4}{4} + C \quad\Rightarrow\quad y = \frac{x^2}{4} + \frac{C}{x^2}

验证:y=x22Cx3y' = \dfrac{x}{2} - \dfrac{2C}{x^3}xy+2y=x222Cx2+x22+2Cx2=x2xy' + 2y = \dfrac{x^2}{2} - \dfrac{2C}{x^2} + \dfrac{x^2}{2} + \dfrac{2C}{x^2} = x^2 ✓(CC 项漂亮相消——好兆头)。

📝 随堂测验

随堂测验0 / 5 题正确

1. 下列哪个是一阶线性方程?

2. 积分因子 $\mu = e^{\int p\,dx}$ 的设计目标是?

3. 解 $xy' + y = 3x^2$ 时第一步应该做什么?

4. 乘完积分因子后发现左边折叠不成 $(\mu y)'$,最可能说明?

5. 咖啡方程既能用分离变量又能用积分因子解,这件事的价值是?

本课小结

  • 一阶线性标准形式 y+p(x)y=q(x)y' + p(x)y = q(x)yyyy' 一次幂不相乘;认 pp 前必须先站队(yy' 系数化为 1)
  • 积分因子 μ=epdx\mu = e^{\int p\,dx}:乘上去左边必定折叠为 (μy)(\mu y)'——折不动 = 报错信号
  • 四步流程:标准化 → 算 μ\mu(不带 CC)→ 折叠 → 积分(此时才有 CC
  • 第 1 课悬案告破:y=x2+C/xy = x^2 + C/xμ=x\mu = x 两步推出;从"会验证"升级为"会生产"
  • 既线性又可分离的方程(如咖啡方程):两种解法互为验算

延伸阅读(可选,跳过不影响主线)

  • Paul's Online Notes · Linear Equations——从头推导积分因子的来历,然后连做五道由易到难的例题,第四、五题的系数比本课更凶悍,适合练手感(英文,跳过不影响主线)

可分离管"乘法脸",线性管"标准队形"——那 y=yy2y' = y - y^2 这种谁都不服的怎么办?下一课学换元法:给陌生方程做一场"整容手术",让它变成这两张熟脸之一。顺便剧透:yy2y - y^2 正是人口增长的 Logistic 方程的雏形,第 7 课的主角。

划选正文任意文字可高亮、批注或加入复习卡

讨论

载入中…