【高等数学】第七章 微分方程——第四节 一阶线性微分方程
上一节:【高等数学】第七章 微分方程——第三节 齐次方程
总目录:【高等数学】 目录
文章目录
- 1. 线性方程
- 2. 伯努利方程
1. 线性方程
- 方程
dydx+P(x)y=Q(x)\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y = Q(x)dxdy+P(x)y=Q(x)
叫做一阶线性微分方程,因为它对于未知函数yyy及其导数是一次方程. 如果Q(x)≡0Q(x) \equiv 0Q(x)≡0,那么方程是齐次的
如果Q(x)≢0Q(x) \not\equiv 0Q(x)≡0,那么方程是非齐次的. - 解法
令Q(x)=0Q(x)=0Q(x)=0,写出对应的齐次线性方程
dydx+P(x)y=0\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y = 0dxdy+P(x)y=0
分离变量后,两边积分得
ln∣y∣=−∫P(x)dx+C1\ln \vert y \vert = -\int P(x) \mathrm{d}x + C_1ln∣y∣=−∫P(x)dx+C1
即y=Ce−∫P(x)dxy = C\mathrm{e}^{-\int P(x) \mathrm{d}x}y=Ce−∫P(x)dx
使用常数变易法,求非齐次线性方程的通解
令u=Cu=Cu=C,y=ue−∫P(x)dxy = u\mathrm{e}^{-\int P(x) \mathrm{d}x}y=ue−∫P(x)dx
dydx=u′e−∫P(x)dx−uP(x)e−∫P(x)dx\displaystyle\dfrac{\mathrm{d}y}{\mathrm{d}x} = u'\mathrm{e}^{-\int P(x) \mathrm{d}x} - uP(x)\mathrm{e}^{-\int P(x) \mathrm{d}x}dxdy=u′e−∫P(x)dx−uP(x)e−∫P(x)dx
代入dydx+P(x)y=Q(x)\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y = Q(x)dxdy+P(x)y=Q(x)得
u′e−∫P(x)dx−uP(x)e−∫P(x)dx+P(x)ue−∫P(x)dx=Q(x)u'\mathrm{e}^{-\int P(x) \mathrm{d}x} - uP(x)\mathrm{e}^{-\int P(x) \mathrm{d}x} + P(x)u\mathrm{e}^{-\int P(x) \mathrm{d}x} = Q(x)u′e−∫P(x)dx−uP(x)e−∫P(x)dx+P(x)ue−∫P(x)dx=Q(x)
即u′=Q(x)e−∫P(x)dxu'=Q(x)\mathrm{e}^{-\int P(x) \mathrm{d}x}u′=Q(x)e−∫P(x)dx
u=∫Q(x)e−∫P(x)dxdx+C\displaystyle u=\int Q(x)\mathrm{e}^{-\int P(x) \mathrm{d}x} \mathrm{d}x+ Cu=∫Q(x)e−∫P(x)dxdx+C
因此非齐次线性方程的通解
y=e−∫P(x)dx(∫Q(x)e∫P(x)dxdx+C)=Ce−∫P(x)dx+e−∫P(x)dx∫Q(x)e∫P(x)dxdx\begin{aligned}y &= \mathrm{e}^{-\int P(x) \mathrm{d}x} \left( \int Q(x) \mathrm{e}^{\int P(x) \mathrm{d}x} \mathrm{d}x + C \right)\\ &= C\mathrm{e}^{-\int P(x) \mathrm{d}x} + \mathrm{e}^{-\int P(x) \mathrm{d}x} \int Q(x) \mathrm{e}^{\int P(x) \mathrm{d}x} \mathrm{d}x\end{aligned}y=e−∫P(x)dx(∫Q(x)e∫P(x)dxdx+C)=Ce−∫P(x)dx+e−∫P(x)dx∫Q(x)e∫P(x)dxdx
右边是齐次线性方程的通解,左边是非齐次线性方程的一个特解
2. 伯努利方程
- 方程
dydx+P(x)y=Q(x)yn(n≠0,1)\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y = Q(x)y^n \quad (n \neq 0,1)dxdy+P(x)y=Q(x)yn(n=0,1)
叫做伯努利(Bernoulli)方程.
当n=0n = 0n=0或n=1n = 1n=1时,这是线性微分方程.
当n≠0n \neq 0n=0,n≠1n \neq 1n=1时,这方程不是线性的,但是通过变量的代换,便可把它化为线性的. - 解法
两边同除以yny^nyn得
y−ndydx+P(x)y1−n=Q(x).y^{-n}\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y^{1 - n} = Q(x).y−ndxdy+P(x)y1−n=Q(x).
令z=y1−nz=y^{1-n}z=y1−n,dzdx=(1−n)y−ndydx\dfrac{\mathrm{d}z}{\mathrm{d}x}=(1-n)y^{-n}\dfrac{\mathrm{d}y}{\mathrm{d}x}dxdz=(1−n)y−ndxdy,则可得线性方程
dzdx+(1−n)P(x)z=(1−n)Q(x)\dfrac{\mathrm{d}z}{\mathrm{d}x} + (1-n)P(x)z = (1-n)Q(x)dxdz+(1−n)P(x)z=(1−n)Q(x)
求得通解后,再将z=y1−nz=y^{1-n}z=y1−n回代即可得伯努利方程通解
下一节:
总目录:【高等数学】 目录