一阶线性方程:把左边凑成一个导数

从乘积求导出发理解积分因子,并用一个初值问题检查计算中的每一步。

为什么要乘一个函数

考虑区间 II 上的一阶线性方程

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

其中 p,qp,q 连续。我们希望乘上一个非零函数 μ\mu,使左边变成乘积的导数。比较

(μy)′=μy′+μ′y(\mu y)'=\mu y'+\mu' y

可知,只要 μ′=pμ\mu'=p\mu,目标就能实现。积分因子因此来自乘积求导法则,而不是需要单独背下来的技巧。

用定积分固定常数

设初值为 y(x0)=y0y(x_0)=y_0,记

P(x)=∫x0xp(s) ds,μ(x)=eP(x).P(x)=\int_{x_0}^{x}p(s)\,\mathrm{d}s, \qquad \mu(x)=e^{P(x)}.

这样 μ(x0)=1\mu(x_0)=1,积分常数的位置便很清楚。对 (μy)′=μq(\mu y)'=\mu q 从 x0x_0 积分到 xx,得到

y(x)=e−P(x)[y0+∫x0xeP(t)q(t) dt].y(x)=e^{-P(x)} \left[y_0+\int_{x_0}^{x}e^{P(t)}q(t)\,\mathrm{d}t\right].

公式中,积分变量 tt 与外部变量 xx 各有职责;把它们写清楚,能减少换元和代入时的混淆。

一个可以直接检验的例子

求解 y′+2y=e−xy'+2y=e^{-x},且 y(0)=0y(0)=0。这里 P(x)=2xP(x)=2x,所以

y(x)=e−2x∫0xet dt=e−x−e−2x.y(x)=e^{-2x}\int_0^x e^t\,\mathrm{d}t =e^{-x}-e^{-2x}.

求导后,y′=−e−x+2e−2xy'=-e^{-x}+2e^{-2x},于是 y′+2y=e−xy'+2y=e^{-x};同时 y(0)=0y(0)=0。把答案代回原方程与初值条件,比只检查最后一行的形式更可靠。

别忽略方程所在的区间

若原方程写成 a(x)y′+b(x)y=c(x)a(x)y'+b(x)y=c(x),需要先在 a(x)≠0a(x)\ne0 的区间内除以 a(x)a(x)。遇到零点时,不能直接沿用除法后的方程跨过去。积分因子本身始终为正,因此在系数连续的区间内不会引入新的除零问题。先确定区间,再做计算,是这类题目很有用的起点。

这是一篇用于展示网站功能的示例,可替换为自己的学习笔记或真实项目记录。