一阶线性微分方程

                     

贡献者: addis; 零穹; int256

  • 本文需要更多参考文献。
  • 需要补充例子
预备知识 常微分方程

  1具有以下形式的微分方程叫做一阶线性微分方程(first order linear differential equation)

\begin{equation} \frac{\mathrm{d}{y}}{\mathrm{d}{x}} + p(x)y = f(x)~, \end{equation}
其中函数 $p(x)$ 和 $f(x)$ 是已知的,函数 $y(x)$ 是未知的,称为方程的解。这三个函数的自变量都是实数,函数值可以是实数或者复数。若 $f(x) \equiv 0$,则称方程为齐次的(homogeneous),否则就是非齐次的(inhomogeneous)。所以 $f(x)$ 叫做方程的非齐次项(inhomogeneous term)

1. 线性

   若干函数 $y_1(x), y_2(x), \dots$ 的线性组合2定义为他们分别乘以一个常数再相加

\begin{equation} c_1 y_1(x) + c_2 y_2(x) + \dots~ \end{equation}
由于求导操作是线性(式 6 )的,容易证明若 $y_i(x)$($i = 1, 2,\dots$)都是齐次方程的解,那么它们的线性组合也必定是齐次方程的解。微分方程的这种性质就叫做线性(linear)

2. 一阶常系数齐次线性微分方程

   当式 1 中的 $p(x) \equiv b$ 是常数时,我们说式 1 的方程是常系数的。先来看齐次的情况,方程为

\begin{equation} \frac{\mathrm{d}{y}}{\mathrm{d}{x}} + by = 0~. \end{equation}
也就是说我们想要求一个函数 $y(x)$,使得它的导数和 $y(x)$ 自己成正比。一个最显然的解就是 $y(x) \equiv 0$,但它太简单了没什么意义,我们说这个解是平凡的(trivial)。对于非零解,一个直接的猜测就是自然指数函数 $y(x) = \mathrm{e} ^{ax}$ 其中 $a$ 为常数。把它代入式 3
\begin{equation} ay(x) + by(x) = 0~, \end{equation}
所以消去 $y(x)$ 得 $b = -a$。再考虑到方程的线性,所以我们把这个解乘以任意常数 $C$ 同样是方程的解,所以式 3 的通解为
\begin{equation} y(x) = C \mathrm{e} ^{-bx}~. \end{equation}

例 1 稀释酒精

   一个容器中装有质量为 $M$(千克)的酒精,从某时刻开始以速率 $I$(千克每秒)加水,并以同样的速率从底部排水,保持容器中液体的总质量不变。假设容器中有某种快速搅拌装置,使两种液体时刻保持充分混合,求容器中酒精浓度 $x$(酒精质量比总质量)随时间的变化关系。

   解:令容器中酒精的质量为 $m_1$,水的质量为 $M - m_1$。流出容器的液体速率可以分为两部分:流出酒精的速率为 $xI$,流出水的速率为 $I - xI$。所以容器中酒精质量的变化率为

\begin{equation} m'_1 = -xI~. \end{equation}
酒精浓度定义为 $x = m_1/M$,$M$ 是常数,所以
\begin{equation} x' = \left(\frac{m_1}{M} \right) ' = \frac{m'_1}{M} = -\frac{I}{M}x~. \end{equation}
该微分方程符合式 3 的结构,可以直接写出方程的解为
\begin{equation} x(t) = C \mathrm{e} ^{-It/M}~. \end{equation}
但这里的 $C$ 到底是什么呢?如果我们令开始的时间为 $t = 0$,那么根据题目条件必须满足 $x(0) = 1$,代入得 $C = 1$。

3. 一阶常系数非齐次线性微分方程

   现在,来看 $p(x) \equiv b$ 是常数时非齐次的情况,方程为

\begin{equation} \frac{\mathrm{d}{y}}{\mathrm{d}{x}} + by = f(x)~, \end{equation}
通解为
\begin{equation} y(x)=C \mathrm{e} ^{-bx}+ \mathrm{e} ^{-bx}\int f(x) \mathrm{e} ^{bx} \,\mathrm{d}{x} ~. \end{equation}

推导

   我们用所谓的常数变易法来求常系数非齐次线性方程式 9 的通解。这方法是把式 5 的通解中的 $C$ 换为 $x$ 的未知函数 $u(x)$,即作变换

\begin{equation} y(x)=u \mathrm{e} ^{-bx}~, \end{equation}
于是
\begin{equation} \frac{\mathrm{d}{y}}{\mathrm{d}{x}} =u' \mathrm{e} ^{-bx}-bu \mathrm{e} ^{-bx}~. \end{equation}
式 11 式 12 代入式 9
\begin{equation} u' \mathrm{e} ^{-bx}=f(x)~, \end{equation}
\begin{equation} \,\mathrm{d}{u} =f(x) \mathrm{e} ^{bx} \,\mathrm{d}{x} ~, \end{equation}
两边积分,得
\begin{equation} u(x)=\int f(x) \mathrm{e} ^{bx} \,\mathrm{d}{x} +C~. \end{equation}
式 15 代入式 11 ,便得常系数非齐次线性方程式 9 的通解
\begin{equation} y(x)= \mathrm{e} ^{-bx} \left[\int f(x) \mathrm{e} ^{bx} \,\mathrm{d}{x} +C \right] ~. \end{equation}

   将上式写成两项之和,就得到了式 10

\begin{equation} y(x)=C \mathrm{e} ^{-bx}+ \mathrm{e} ^{-bx}\int f(x) \mathrm{e} ^{bx} \,\mathrm{d}{x} ~. \end{equation}
右端第一项对应常系数齐次线性方程式 3 的通解,第二项则是常系数非齐次线性方程式 9 的一个特解(在通解式 16 中取 $C=0$ 便得)。由此可知,一阶常系数非齐次线性方程的通解等于对应的齐次方程的通解与非齐次方程的一个特解之和。

4. 齐次的常见情况

   我们先来看式 1 对应的齐次方程

\begin{equation} \frac{\mathrm{d}{y}}{\mathrm{d}{x}} + p(x)y = 0~. \end{equation}
这是一个可分离变量的方程(分离变量是指把所有含 $x$ 的项放到一边、含 $y$ 的项放到等式另一边,变为 $P(y) \,\mathrm{d}{y} = Q(x) \,\mathrm{d}{x} $ 的形式),分离变量得
\begin{equation} \frac{ \,\mathrm{d}{y} }{y} = -p(x) \,\mathrm{d}{x} ~, \end{equation}
分离变量使得两边可以分别对各自的变量积分。两边分别积分得
\begin{equation} \ln \left\lvert y \right\rvert = -\int p(x) dx + C~. \end{equation}
两边取自然指数得
\begin{equation} y = \pm \mathrm{e} ^C \mathrm{e} ^{-\int p(x) \,\mathrm{d}{x} }~. \end{equation}
把 $\pm \mathrm{e} ^C $ 整体看做一个任意常数 $C$,上式变为。
\begin{equation} y = C \mathrm{e} ^{-\int p(x) \,\mathrm{d}{x} }~. \end{equation}
这就是一阶线性齐次微分方程式 18 的通解,也叫式 1 齐次解

5. 非齐次的一般形式

   重新考虑回式 1 的情况:

\begin{equation} y'+ P(x) y = Q(x) ~. \end{equation}
我们已经知道其中 $y' + P(x) y = 0$ 的通解是 $C \exp\left(- \int P(x) \,\mathrm{d}{x} \right) $(这里 $ \exp\left(f\right) = e^f$)。而对于现在的这个方程我们只需要找到一个 “特解” 即可。也就是俗称的 “解集=通解+特解”。

   做一个常见假设:考虑 $y$ 还可以表示为 “分离变量” 的形式,即 $y=(C+\xi) \exp\left(- \int P(x) \,\mathrm{d}{x} \right) $。其中 $\xi=\xi(x)$ 贡献特解项。那么代回式 23 可以得到:

\begin{equation} \xi = \int \left[Q(x) \exp\left(\int P(x) \,\mathrm{d}{x} \right) \right] \,\mathrm{d}{x} ~. \end{equation}

   所以有,对于这一般形式的微分方程,有通解:

\begin{equation} y = \exp\left(- \int P(x) \,\mathrm{d}{x} \right) \left[\int\left(Q(x) \exp\left[ \int P(x) \,\mathrm{d}{x} \right] \,\mathrm{d}{x} \right) + C\right]~. \end{equation}

6. 常数变易法

   现在我们用常数变易法来解非齐次方程式 1 。为书写方便,式 22 中令 $y_0(x) = \exp\left(-\int p(x) \,\mathrm{d}{x} \right) $。假设上式中的 $C$ 是一个函数 $C(x)$ 而不是常数,代入式 1

\begin{equation} C'y_0 + C[y_0' + p(x)y_0] = f(x)~. \end{equation}
由于 $y_0$ 是齐次解,上式方括号中求和为 0,分离变量得
\begin{equation} \,\mathrm{d}{C} = \frac{f(x)}{y_0} \,\mathrm{d}{x} ~, \end{equation}
两边积分得
\begin{equation} C(x) = \int \frac{f(x)}{y_0} \,\mathrm{d}{x} ~, \end{equation}
所以一阶线性非齐次微分方程的通解为
\begin{equation} y = y_0 \int \frac{f(x)}{y_0} \,\mathrm{d}{x} ~, \end{equation}
其中
\begin{equation} y_0(x) = \mathrm{e} ^{-\int p(x) \,\mathrm{d}{x} }~. \end{equation}
注意待定常数包含在式 29 的不定积分中,式 30 中的不定积分产生的待定常数在代入式 29 后可消去。


1. ^ 本文参考 Wikipedia 相关页面[1]
2. ^ 线性组合是线性代数中的概念,线性代数中,线性空间都可以做线性组合,例如几何矢量(式 14 )。特定函数的集合可以看作线性空间(例 2 )。


[1] ^ 同济大学数学系. 高等数学 (上下册) 高等教育出版社 (2014) 第七版

                     

© 小时科技 保留一切权利