4. 一阶线性微分方程

4. 一阶线性微分方程

4.1 线性方程

4.1.1 基本概念

对于方程 $\frac{dy}{dx}+P(x)y=Q(x)$（这是一个一阶线性微分方程）

（这个方程对于未知函数y及其导数是 一次方程）

① $Q(x)\equiv 0$ 时，此方程为 齐次方程

② $Q(x)$ 不恒为0时，此方程为 非齐次方程

4.1.2 如何解此类方程？

4.1.2.1 齐次的情况

此时方程为 $\frac{dy}{dx}+P(x)y=0$，显然这时候方程是一个可分离变量的方程。

我们分离变量得到 $\frac{dy}{y}=-P(x)dx$

对两端积分得到 $\ln |y|=-\int P(x)dx+C$

变成左端只剩 y 的样子就是 $y=C{e}^{-\int P(x)dx}$（这里动手试一下，是可以化成这种样子的）

此时得到的 $y=C{e}^{-\int P(x)dx}$ 就是方程 $\frac{dy}{dx}+P(x)y=0$ 的通解

4.1.2.2 非齐次的情况

此时我们要求解的式子是 $\frac{dy}{dx}+P(x)y=Q(x)$

为了简洁，我们将式子写成这个样子 $y^{{}^{\prime }}+py=q$

不难发现，左侧式子类似于乘法求导公式，只不过是少了一个因子。

不妨两端同乘一个 u(x) ，(这里同样为了简便写成 u )，即 $u{y}^{{}^{\prime }}+puy=q$

又因为 $(yu{)}^{\prime }=y^{{}^{\prime }}u+u^{{}^{\prime }}y$

此时只需要找到一个 $u^{{}^{\prime }}$ 使之等价于 $pu$，即得解。

此时解 $u^{{}^{\prime }}=pu$

即 $\frac{du}{u}=P(x)dx$

由上面讨论的齐次情况，我们可以得到 $u=e^{\int P(x)dx}$

此时 $uy=\int qu+C$

那么就得到 $y=\frac{\int qu+C}{u}$

同济版高等数学上的“常数变易法”是这样的：将齐次方程的通解 $y=C{e}^{-\int P(x)dx}$ 中的 C 换成关于 x 的未知函数 u(x) ，作变换 $y=u{e}^{-\int P(x)dx}$

4.2 伯努利方程

4.2.1 基本概念

对于方程 $\frac{dy}{dx}+P(x)y=Q(x)y^n$，$(n\ne 0,1)$

①当 $n=0$ 或 $1$ 的时候，方程为线性微分方程

②当 $n\ne 0,1$ 的时，方程显然不为线性，但是可以通过变量代换将其变换为线性的。

4.2.2 如何解?

$\frac{dy}{dx}+P(x)y=Q(x)y^n$，$(n\ne 0,1)$

方程两端同时除以 $y^n$

$\frac{dy}{dx}\cdot y^{-n}+P(x)y^{1-n}=Q(x)$

设 $z=y^{1-n}$ ，且 $\frac{dz}{dx}=(1-n)y^{-n}\cdot \frac{dy}{dx}$

那么此时原方程变为 $\frac{dz}{dx}+(1-n)P(x)\cdot z=(1-n)Q(x)$

此时解这个线性方程得到解之后，再用 $z=y^{1-n}$ 反代就得到了方程的通解。

4.3 例题

1. 求方程 $\frac{dy}{dx}-\frac{2y}{x+1}=(x+1{)}^{\frac{5}{2}}$ 的通解

运用上面得到的解法， $P(x)=-\frac{2}{x+1}$

那么 $u=e^{\int P(x)dx}=e^{-\int \frac{2}{x+1}dx}=(x+1{)}^{-2}$

则 $$\begin{gathered} y=\frac{\int qu+C}{u} \\ =\frac{\int (x+1{)}^{\frac{5}{2}}\cdot (x+1{)}^{-2}dx+C}{(x+1{)}^{-2}} \\ =(x+1{)}^2[\int (x+1{)}^{\frac{1}{2}}dx+C] \\ =(x+1{)}^2[\frac{2}{3}(x+1{)}^{\frac{3}{2}}+C] \end{gathered}$$

2. 解方程 $\frac{dy}{dx}=\frac{1}{x+y}$

此时令 $u=x+y$ ，那么

$\frac{du}{dx}=\frac{dy}{dx}+1$

即 $\frac{du}{dx}=\frac{1}{u}+1$

分离变量就得到 $\frac{u}{u+1}du=dx$

两端同时积分得到 $u-\ln |u+1|=x+C$

此时反代 $u=x+y$ 得到通解。

3. 求方程 $\frac{dy}{dx}+\frac{y}{x}=alnx\cdot y^2$ 的通解

方程两端同除 $y^2$，方程变为 $\frac{dy}{dx}\cdot \frac{1}{y^2}+\frac{1}{xy}=alnx$

此时令 $z=\frac{1}{y}$ ，则 $\frac{dz}{dx}=-\frac{1}{y^2}\cdot \frac{dy}{dx}$

代回原式就有 $\frac{dz}{dx}-\frac{z}{x}=-alnx$

此时解此一阶线性微分方程即可。

$$\begin{gathered} u=e^{\int -\frac{1}{x}dx}=\frac{1}{x} \\ z=\frac{\int alnx\cdot \frac{1}{x}dx+C}{\frac{1}{x}}=x[C-\frac{a}{2}(lnx{)}^2] \end{gathered}$$

又 $z=\frac{1}{y}$ ，再反代即得解。