跳至主要內容

微分方程


微分方程

微分方程的概念

微分方程的通解

注意

对于求解一些微分方程通解的时候,通常会有特殊点(奇解),如果能让常数取到某一值而让方程包含这一特殊值,就让在确定常数范围时加上这一值,如果不能就直接不用管。

一阶微分方程的求解

变量可分离型

可化为变量可分离型

  1. (1) 形如 dydx=f(ax+by+c)\displaystyle \frac{\mathrm{d} y}{\mathrm{d} x}=f(a x+b y+c) 的方程, 其中常数 a,ba, b 全都不为零. 其解法为令 u=ax+by+cu=a x+b y+c, 则 dudx=a+bdydx\displaystyle \frac{\mathrm{d} u}{\mathrm{d} x}=a+b \frac{\mathrm{d} y}{\mathrm{d} x}, 代人原方程得 dudx=a+bf(u)\displaystyle \frac{\mathrm{d} u}{\mathrm{d} x}=a+b f(u).
  2. (2) 齐次型微分方程.
    如, 对于微分方程
    将其变形为

(x+ycosyx)dxxcosyxdy=0, \displaystyle \left(x+y \cos \frac{y}{x}\right) \mathrm{d} x-x \cos \frac{y}{x} \mathrm{d} y=0,

dydx=x+ycosyxxcosyx=1+yxcosyxcosyx=φ(yx). \displaystyle \frac{\mathrm{d} y}{\mathrm{d} x}=\frac{x+y \cos \frac{y}{x}}{x \cos \frac{y}{x}}=\frac{1+\frac{y}{x} \cos \frac{y}{x}}{\cos \frac{y}{x}}=\varphi\left(\frac{y}{x}\right) .

形如 dydx=φ(yx)\displaystyle \frac{\mathrm{d} y}{\mathrm{d} x}=\varphi\left(\frac{y}{x}\right) 的方程叫作齐次型微分方程. 其解法为令 u=yxu=\frac{y}{x}, 则

y=uxdydx=u+xdudx, \displaystyle y=u x \Rightarrow \frac{\mathrm{d} y}{\mathrm{d} x}=u+x \frac{\mathrm{d} u}{\mathrm{d} x},

于是原方程变为 xdudx+u=φ(u)\displaystyle x \frac{\mathrm{d} u}{\mathrm{d} x}+u=\varphi(u), 即 duφ(u)u=dxx\displaystyle \frac{\mathrm{d} u}{\varphi(u)-u}=\frac{\mathrm{d} x}{x}.

高阶线性微分方程的求解

nn 阶常系数齐次线性微分方程的解

方程 y(n)+p1y(n1)++pn1y+pny=0y^{(n)}+p_{1} y^{(n-1)}+\cdots+p_{n-1} y{\prime}+p_{n} y=0 称为 nn 阶常系数齐次线性微分方程, 其中 p1,p2,,pnp_{1}, p_{2}, \cdots, p_{n} 为常 数, 其对应的特征方程为 rn+p1rn1++pn1r+pn=0r^{n}+p_{1} r^{n-1}+\cdots+p_{n-1} r+p_{n}=0, 求出其特征根. 有如下情况需要大家牢记 (其中 大写的英文字母均为任意常数).

  1. 特征根为单实根 rr 时, 微分方程通解中对应一项 CerxC \mathrm{e}^{r x};
  2. 特征根为 kk 重实根 rr 时, 微分方程通解中对应 kk(C1+C2x++Ckxk1)erx\left(C_{1}+C_{2} x+\cdots+C_{k} x^{k-1}\right) \mathrm{e}^{r x};
  3. 特征根为单复根 α±βi(β>0)\alpha \pm \beta i(\beta>0) 时, 微分方程通解中对应两项 ear(C1cosβx+C2sinβx)\mathrm{e}^{\mathrm{ar}}\left(C_{1} \cos \beta x+C_{2} \sin \beta x\right);
  4. 特征根为 kk 重复根 α±βi(β>0)\alpha \pm \beta i(\beta>0) 时, 微分方程通解中对应 2k2 k

eax[(C1+C2x++Ckxk1)cosβx+(D1+D2x++Dkxk1)sinβx]. \mathrm{e}^{\mathrm{ax}}\left[\left(C_{1}+C_{2} x+\cdots+C_{k} x^{k-1}\right) \cos \beta x+\left(D_{1}+D_{2} x+\cdots+D_{k} x^{k-1}\right) \sin \beta x\right] .

一、一阶线性微分方程

方程

dydx+P(x)y=Q(x) \frac{dy}{dx}+P(x)y=Q(x)

叫做一阶线性微分方程,如果Q(x)0Q(x)\equiv0,则方程称为齐次,反之则为非齐次。

1.齐次线性方程

dydx+P(x)y=0 \displaystyle\frac{dy}{dx}+P(x)y=0

分离变量

dyy=P(x)dx \displaystyle\frac{dy}{y}=-P(x)dx

两端积分,整理可得

y=CeP(x)dx \displaystyle y=\displaystyle Ce^{-\displaystyle\int \displaystyle P(x)dx}

2.非齐次线性方程

dydx+P(x)y=Q(x) \frac{dy}{dx}+P(x)y=Q(x)

由常数易变法,可得通解

y=eP(x)dx(Q(x)eP(x)dxdx+C) y=e^{-{\displaystyle\int} P(x)dx}(\int Q(x)e^{\displaystyle\int P(x)dx}dx+C)

推导过程

在方程y+p(x)y=q(x)y'+p(x)y=q(x) 两边同时乘以 ep(x)dx\mathrm{e}^{\displaystyle\int p(x) \mathrm{d} x}, 得

ep(x)dxy+ep(x)dxp(x)y=ep(x)dxq(x), \mathrm{e}^{\int p(x) \mathrm{d} x} \cdot y{\prime}+\mathrm{e}^{\int p(x) \mathrm{d} x} p(x) \cdot y=\mathrm{e}^{\int p(x) \mathrm{d} x} \cdot q(x),

[ep(x)dxy]=ep(x)dxq(x) {\left[\mathrm{e}^{\int p(x) \mathrm{d} x} \cdot y\right]{\prime}=\mathrm{e}^{\int p(x) \mathrm{d} x} \cdot q(x)}

ep(x)dxy=ep(x)dxq(x)dx+C \mathrm{e}^{\int p(x) \mathrm{d} x} \cdot y=\int \mathrm{e}^{\int p(x) \mathrm{d} x} \cdot q(x) \mathrm{d} x+C

y=ep(x)dx[ep(x)dxq(x)dx+C] y=\mathrm{e}^{-\int p(x) \mathrm{d} x}\left[\int \mathrm{e}^{\int p(x) \mathrm{d} x} \cdot q(x) \mathrm{d} x+C\right]

二、伯努利方程(仅数一)

方程

dydx+P(x)y=Q(x)yn(n0,1) \frac{dy}{dx}+P(x)y=Q(x)y^n\quad(n\neq0,1)

叫做伯努利方程.

方程两边同除yny^n,再令z=y1nz=y^{1-n},可化为线性方程:

dzdx+(1n)P(x)z=(1n)Q(x) \frac{dz}{dx}+(1-n)P(x)z=(1-n)Q(x)

三、可降阶的高阶微分方程

1.y(n)=f(x)y^{(n)}=f(x)

连续积分即可.

2.y=f(x,y)y^{''}=f(x,y^{'})

y=py^{'}=p,则y=dpdx=py^{''}=\displaystyle\frac{dp}{dx}=p^{'}

可化为一阶微分方程.

3.y=f(y,y)y^{''}=f(y,y^{'})

y=py^{'}=p,则y=dpdx=dpdydydx=pdpdyy^{''}=\displaystyle\frac{dp}{dx}=\frac{dp}{dy}\cdot\frac{dy}{dx}=p\frac{dp}{dy}

可化为一阶微分方程.

四、高阶线性微分方程

1.线性微分方程解的结构

先讨论二阶齐次线性方程:

y+P(x)y+Q(x)y=0 y^{''}+P(x)y^{'}+Q(x)y=0

定理1

如果函数y1(x)y_1(x)y2(x)y_2(x)是上述方程的两个解,那么

y=C1y1(x)+C2y2(x) y=C_1y_1(x)+C_2y_2(x)

也是解,其中C1,C2C_1,C_2为任意常数.

定理2

如果函数y1(x)y_1(x)y2(x)y_2(x)是上述方程的两个线性无关的解,那么

y=C1y1(x)+C2y2(x) y=C_1y_1(x)+C_2y_2(x)

就是方程的通解,其中C1,C2C_1,C_2为任意常数.

推论

如果函数y1(x)y_1(x),y2(x)y_2(x),...yn(x)y_n(x)nn阶齐次线性方程的nn个线性无关的解,那么

y=C1y1(x)+C2y2(x)++Cn(x) y=C_1y_1(x)+C_2y_2(x)+\cdots+C_n(x)

就是方程的通解.

定理3

y(x)y^{*}(x)是二阶非齐次线性方程

y+P(x)y+Q(x)y=f(x) y^{''}+P(x)y^{'}+Q(x)y=f(x)

的一个特解.Y(x)Y(x)是对应的齐次方程的通解,则

y=Y(x)+y(x) y=Y(x)+y^{*}(x)

是二阶非齐次线性方程的通解.

定理4(叠加原理)

设非齐次线性方程的右端f(x)f(x)是两个函数之和,即

y+P(x)y+Q(x)y=f1(x)+f2(x) y^{''}+P(x)y^{'}+Q(x)y=f_1(x)+f_2(x)

y1(x)y_1^{*}(x)y2(x)y_2^{*}(x)分别是方程

y+P(x)y+Q(x)y=f1(x) y^{''}+P(x)y^{'}+Q(x)y=f_1(x)

y+P(x)y+Q(x)y=f2(x) y^{''}+P(x)y^{'}+Q(x)y=f_2(x)

的特解,则y1(x)+y2(x)y_1^{*}(x)+y_2^{*}(x)就是原方程的特解.

五、常数变易法

1.对于一阶线性非齐次微分方程:

dydx+P(x)y=P(x) \frac{dy}{dx}+P(x)y=P(x)

对应的齐次方程的通解为

y=CeP(x)dx y=Ce^{-{\displaystyle\int} P(x)dx}

将其中的常数CC换成xx的未知函数u(x)u(x),即

y=ueP(x)dx y=ue^{-{\displaystyle\int} P(x)dx}

可得

dydx=ueP(x)dxuP(x)eP(x)dx \frac{dy}{dx}=u^{'}e^{-{\displaystyle\int} P(x)dx}-uP(x)e^{-{\displaystyle\int} P(x)dx}

最终整理可得到

y=eP(x)dx(Q(x)eP(x)dxdx+C) y=e^{-{\displaystyle\int} P(x)dx}(\int Q(x)e^{\displaystyle\int P(x)dx}dx+C)

2.对于二阶线性非齐次微分方程

y+P(x)y+Q(x)y=f(x) y^{''}+P(x)y^{'}+Q(x)y=f(x)

其对应的齐次方程的通解:

y=C1y1(x)+C2y2(x) y=C_1y_1(x)+C_2y_2(x)

将其中的常数换为关于xx的函数v1(x),,v2(x)v_1(x),,v_2(x),即

y=v1y1(x)+v2y2(x) y=v_1y_1(x)+v_2y_2(x)

求导可得:

y=y1v1+y2v2+y1v1+y2v2 y^{'}=y_1v_1^{'}+y_2v_2^{'}+y_1^{'}v_1+y_2^{'}v_2