常见的常微分方程的一般解法

常见的常微分方程的一般解法总结了常见常微分方程的通解。

如无意外，本文将不包括解的推导过程。

常微分方程，我们一般可以将其归纳为如下n类：

1.可分离变量的微分方程（一阶）

2.一阶齐次（非齐次）线性微分方程（一阶），包含伯努

利

3.二阶常系数微分方程（二阶）

4.高阶常系数微分方程（n阶），包含欧拉

1.可分离变量的微分方程（一阶）

这类微分方程可以变形成如下形式：

f ( x ) d x =

g ( y ) d y f(x)dx=g(y)dy f(x)dx=g(y)dy

函数可以通过同时整合两边来解决。难点主要在于不定积分，不定积分是最简单的微分方程。

p.s. 某些方程看似不可分离变量，但是经过换元之后，其实还是可分离变量的，不要被这种方程迷惑。

2.一阶齐次（非齐次）线性微分方程（一阶）

形如

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

的方程叫做一阶线性微分方程，若 Q ( x ) Q(x) Q(x)为0，则方程齐次，否则称为非齐次。

解法：直接套公式：

y ( x ) = e − ∫ P ( x ) d x ( ∫ e ∫ P ( x ) d x Q ( x ) d x + C ) y(x)=e^{-

\int{P(x)}dx}(\int{e^{\int{P(x)dx}}Q(x)}dx+C)

y(x)=e−∫P(x)dx(∫e∫P(x)dxQ(x)dx+C)

多套几遍熟练就好。

伯努利方程

形如

d y d x + P ( x ) y = Q ( x ) y n , n ∈R , n ≠ 1

\frac{dy}{dx}+P(x)y=Q(x)y^{n},n\in\mathbb{R},n\ne1

dxdy+P(x)y=Q(x)yn,n∈R,n=1

的方程称为伯努利方程，这种方程可以通过以下步骤化为一阶线性微分方程：

y − n d y d x + P ( x ) y 1 − n = Q ( x ) y^{-

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

+P(x)y1−n=Q(x)

1 1 − n ⋅ d y 1 − n d x + P ( x ) y 1 − n = Q ( x ) \frac{1}{1-n}·\frac{dy^{1-n}}{dx}+P(x)y^{1-n}=Q(x)

1−n1⋅dxdy1−n+P(x)y1−n=Q(x)

令 y 1 − n = u y^{1-n}=u y1−n=u，方程两边同时乘以 1 − n 1-n 1−n，得到

d u d x + ( 1 − n ) P ( x ) u = ( 1 − n ) Q ( x )

\frac{du}{dx}+(1-n)P(x)u=(1-n)Q(x) dxdu

+(1−n)P(x)u=(1−n)Q(x)

即 d u d x + P ′ ( x ) u = Q ′ ( x )

\frac{du}{dx}+P'(x)u=Q'(x) dxdu+P′(x)u=Q′(x)

这是一个可以公式化的一阶线性微分方程。

3.二阶常系数微分方程（二阶）

形如

y ′ ′ + p y ′ + q y = f ( x ) y''+py'+qy=f(x) y′′+py′+qy=f(x)

的方程称为二阶常系数微分方程，若 f ( x ) ≡ 0

f(x)\equiv0 f(x)≡0，则方程称为齐次的，反之称为非齐次的。以下默认方程是非齐次的。

求解此类方程分两步：

1.求出齐次通解

2.求出非齐次特解

方程的解是齐次通解和非齐次特解之和。

齐次通解的求法

首先假设 f ( x ) ≡ 0 f(x)\equiv0 f(x)≡0.

用特征方程法，写出对应的特征方程并且求解：

r 2 + p r + q = 0 r^{2}+pr+q=0 r2+pr+q=0

解的情况分为以下三种：

情况一：方程有两个不同的实数解

假设两个实数解分别是 r 1 、 r 2 r_{1}、r_{2} r1、r2，

此时方程的通解是

Y ( x ) = C 1 e r 1 x + C 2 e r 2 x

Y(x)=C_{1}e^{r_{1}x}+C_{2}e^{r_{2}x} Y(x)=C1er1x+C2er2 x.

情况二：方程有一个二重解

假设该解等于 r r r，

此时方程的通解是

Y ( x ) = ( C 1 + C 2 x ) e r x

Y(x)=(C_{1}+C_{2}x)e^{rx} Y(x)=(C1+C2x)erx.

情况三：方程有一对共轭复解

假设这对解是α ± i β \alpha\pm i\beta α±iβ

此时方程的通解是

Y ( x ) = e α x ( C 1 c o s ( β x ) + C 2 s i n ( β x ) ) Y(x)=e^{\alpha x}(C_{1}cos(\beta

x)+C_{2}sin(\beta x)) Y(x)=eαx(C1cos(βx)+C2sin(βx))非齐次特解的求法

对于 f ( x ) f(x) f(x)和特征根的情况，对特解的情况做如下归纳：

1. f ( x ) = P m ( x ) f(x)=P_{m}(x) f(x)=Pm(x)，其中P m ( x ) P_{m}(x) Pm(x)表示 x x x的最高次数为m的多

项式。

若0不是方程特征解

则方程有特解 y ∗ = Q m ( x ) y^{*}=Q_{m}(x) y∗=Qm(x)若0是方程的单特征解

则方程有特解 y ∗ = x Q m ( x ) y^{*}=xQ_{m}(x)

y∗=xQm(x)

若0是方程的二重特征解

则方程有特解 y ∗ = x 2 Q m ( x ) y^{*}=x^{2}Q_{m}(x)

y∗=x2Qm(x).

其中，Q m ( x ) = b 0 + b 1 x + … … + b m x m

Q_{m}(x)=b_{0}+b_{1}x+……+b_{m}x^{m} Qm(x)=b0+b1

x+……+bmxm，

b i ( i = 0 , 1 , … … m ) b_{i}(i = 0,1,……m) bi

(i=0,1,……m)是需要带回原方程来确定的系数。

2. f ( x ) = e α x P m ( x ) f(x)=e^{\alpha x}P_{m}(x) f(x)=eαx Pm(x).

若α \alpha α不是方程特征解

则方程有特解 y ∗ = e α x Q m ( x ) y^{*}=e^{\alpha x}Q_{m}(x) y∗=eαx Qm(x)

若α \alpha α是方程的单特征解

则方程有特解 y ∗ = x e α x Q m ( x )

y^{*}=xe^{\alpha x}Q_{m}(x) y∗=xeαx Qm(x)

若α \alpha α是方程的二重特征解

则方程有特解 y ∗ = x 2 e α x Q m ( x )

y^{*}=x^{2}e^{\alpha x}Q_{m}(x) y∗=x2eαx Qm(x).

3. f ( x ) = e α x ( a 1 c o s ( β x ) + a 2 s i n ( β x ) ) f(x)=e^{\alpha x}(a_{1}cos(\beta

x)+a_{2}sin(\beta x)) f(x)=eαx(a1cos(βx)+a2sin(βx))若α ± i β \alpha\pm i\beta α±iβ不是特征解

则方程有特解 y ∗= e α x ( A 1 c o s ( β x ) + A 2 s i n ( β x ) ) y^{*}=e^{\alpha x}(A_{1}cos(\beta x)+A_{2}sin(\beta x)) y∗=eαx(A1cos(βx)+A2sin(βx))

若α ± i β \alpha\pm i\beta α±iβ是特征解

则方程有特解 y ∗= x e α x ( A 1 c o s ( β x ) + A 2 s i n ( β x ) ) y^{*}=xe^{\alpha x}(A_{1}cos(\beta x)+A_{2}sin(\beta x)) y∗=xeαx(A1cos(βx)+A2sin(βx))

其中 A 1 A_{1} A1、 A 2 A_{2} A2是需要带回原方程来确定的系数。

4.高阶常系数微分方程（n阶），包含欧拉

形如 y ( n ) + p 1 y ( n − 1 ) + . . . + p n − 1

y ′ + p n y = f ( x ) y^{(n)}+p_{1}y^{(n-

1)}+...+p_{n-1}y'+p_{n}y=f(x) y(n)+p1y(n−1)+...+pn−1 y′+pny=f(x)

的方程叫做高阶常系数微分方程，若 f ( x ) ≡ 0

f(x)\equiv0 f(x)≡0，则方程是齐次的，否则是非齐次的。下面默认方程是非齐次的。

求解此类方程分两步：

1.求出齐次通解

2.求出非齐次特解

方程的解是齐次通解和非齐次特解之和。

齐次通解的求法

使用特征方程法。首先假设 f ( x ) ≡ 0 f(x)\equiv0 f(x)≡0.

列出特征方程并求解：

r n + p 1 r n − 1 + . . . + p n − 1 r + p n = 0

r^{n}+p_{1}r^{n-1}+...+p_{n-1}r+p_{n}=0 rn+p1

rn−1+...+pn−1r+pn=0

得到的解无非是以下四种情形：

1.1重实数解

2.k重实数解

3.一对1重共轭复数解

4.一对k重共轭复数解

每出现一种情形的解，通解 y ( x ) y(x) y(x)中就多出一项，各项之和就构成了整个通解。

根据下面的规则可以写出不同情况下的通解

1重实数解

设λ \lambda λ是特征方程的一个1重实数解。

则应该为它在通解中增加这样一项：

C e λ x Ce^{\lambda x} Ceλx

k重实数解

设λ \lambda λ是特征方程的一个k重实数解。

则应该为它在通解中增加这样一项（如果拆开括号就是k

项）：

e λ x ( C 1 + C 2 x + . . . + C k x k − 1 )

e^{\lambda x}(C_{1}+C_{2}x+...+C_{k}x^{k-1}) eλx(C1

+C2x+...+Ckxk−1)

一对1重共轭复数解

设α ± i β \alpha\pm i\beta α±iβ是特征方程的一

对1重共轭复数解。

则应该为它在通解中增加这样一项（如果拆开括号就是2

项）：

e α x ( C 1 c o s ( β x ) + C 2 s i n ( β x ) )

e^{\alpha x}(C_{1}cos(\beta x)+C_{2}sin(\beta x))

eαx(C1cos(βx)+C2sin(βx))

一对k重共轭复数解

设α ± i β \alpha\pm i\beta α±iβ是特征方程的一

对k重共轭复数解。

则应该为它在通解中增加这样一项（如果拆开括号就是2k 项）：

e α x [ ( C 1 + C 2 x + . . . + C k x k − 1 ) c o s ( β x ) + e^{\alpha x}[(C_{1}+C_{2}x+...+C_{k}x^{k-1})cos(\beta x)+ eαx[(C1+C2x+...+Ckxk−1)cos(βx)+

( D 1 + D 2 x + . . . + D k x k − 1 ) s i n ( β x ) ] (D_{1}+D_{2}x+...+D_{k}x^{k-1})sin(\beta x)] (D1+D2

x+...+Dkxk−1)sin(βx)]

当你写出每个特征方程的解的相应项时，把它们加起来，你

就得到齐次方程的通解。

非齐次特解的求法

1. f ( x ) = P m ( x ) f(x)=P_{m}(x) f(x)=Pm(x)，其中

P m ( x ) P_{m}(x) Pm(x)表示 x x x的最高次数为m的多

项式。

若0不是方程特征解

则方程有特解 y ∗ = Q m ( x ) y^{*}=Q_{m}(x) y∗=Qm(x)若0是方程的k重特征解

则方程有特解 y ∗ = x k Q m ( x ) y^{*}=x^{k}Q_{m}(x)

y∗=xkQm(x)

其中，Q m ( x ) = b 0 + b 1 x + … … + b m x m

Q_{m}(x)=b_{0}+b_{1}x+……+b_{m}x^{m} Qm(x)=b0+b1

x+……+bmxm，

b i ( i = 0 , 1 , … … m ) b_{i}(i = 0,1,……m) bi

(i=0,1,……m)是需要带回原方程来确定的系数。

2. f ( x ) = e α x P m ( x ) f(x)=e^{\alpha x}P_{m}(x) f(x)=eαx Pm(x).

若α \alpha α不是方程特征解

则方程有特解 y ∗ = e α x Q m ( x ) y^{*}=e^{\alpha x}Q_{m}(x) y∗=eαx Qm(x)

若α \alpha α是方程的k重特征解

则方程有特解 y ∗ = x k e α x Q m ( x )

y^{*}=x^{k}e^{\alpha x}Q_{m}(x) y∗=xkeαx Qm(x).

3. f ( x ) = e α x ( a 1 c o s ( β x ) + a 2 s i n ( β x ) ) f(x)=e^{\alpha x}(a_{1}cos(\beta

x)+a_{2}sin(\beta x)) f(x)=eαx(a1cos(βx)+a2sin(βx))若α ± i β \alpha\pm i\beta α±iβ不是特征解

则方程有特解 y ∗= e α x ( A 1 c o s ( β x ) + A 2 s i n ( β x ) ) y^{*}=e^{\alpha x}(A_{1}cos(\beta x)+A_{2}sin(\beta x)) y∗=eαx(A1cos(βx)+A2sin(βx))

若α ± i β \alpha\pm i\beta α±iβ是k重特征解

则方程有特解 y ∗= x k e α x ( A 1 c o s ( β x ) + A 2 s i n ( β x ) ) y^{*}=x^{k}e^{\alpha

x}(A_{1}cos(\beta x)+A_{2}sin(\beta x)) y∗=xkeαx(A1 cos(βx)+A2sin(βx))

其中 A 1 A_{1} A1、 A 2 A_{2} A2是需要带回原方程来确定的系数。

欧拉方程

形如

x n y ( n ) + p 1 x n − 1 y ( n − 1 ) + . . . + p n − 1 x y ′ + p n y = f ( x ) x^{n}y^{(n)}+p_{1}x^{n-1}y^{(n-1)}+...+p_{n-1}xy'+p_{n}y=f(x) xny(n)+p1

xn−1y(n−1)+...+pn−1xy′+pny=f(x)

我们的方程叫欧拉方程，可以转化为常系数微分方程。

变形过程如下：

做变换： x = e t x=e^{t} x=et, 或 t = l n ( x ) t=ln(x) t=ln(x)

于是

d y d x = d t d x d y d x = 1 x d y d t

\frac{dy}{dx}=\frac{dt}{dx}\frac{dy}{dx}=\frac{1}{x}\f rac{dy}{dt} dxdy=dxdtdxdy=x1dtdy

d 2 y d x 2 = 1 x 2 ( d 2 y d t 2 − d y d t )

\frac{d^{2}y}{dx^{2}}=\frac{1}{x^{2}}(\frac{d^2y}{dt^2 }-\frac{dy}{dt}) dx2d2y=x21(dt2d2y−dtdy)

d 3 y d x 3 = 1 x 3 ( d 3 y d t 3 − 3 d 2 y d t 2 + 2 d y d t )

\frac{d^{3}y}{dx^{3}}=\frac{1}{x^{3}}(\frac{d^3y}{dt^3 }-3\frac{d^2y}{dt^2}+2\frac{dy}{dt}) dx3d3y=x31(dt3d3y −3dt2d2y+2dtdy)

……

记 D k y = d k y d t k D^{k}y=\frac{d^{k}y}{dt^{k}} Dky=dtkdky

则x y ′ = D y xy'=Dy xy′=Dy

x 2 y ′ ′ = D ( D − 1 ) y x^{2}y''=D(D-1)y

x2y′′=D(D−1)y

x 3 y ( 3 ) = D ( D − 1 ) ( D − 2 ) y

x^{3}y^{(3)}=D(D-1)(D-2)y x3y(3)=D(D−1)(D−2)y

…

一般的， x k y ( k ) = D ( D − 1 ) . . . ( D − k + 1 ) y x^{k}y^{(k)}=D(D-1)...(D-k+1)y

xky(k)=D(D−1)...(D−k+1)y.

这样，欧拉方程就是关于 y y y和 t t t的常系数微分方程。

解出 y y y后，将 t = l n ( x ) t=ln(x) t=ln(x)带入即得 y ( x ) y(x) y(x).

如有错误欢迎指出。