一阶线性微分方程的初等积分法
变量分离方程与变量变换
可分离变量的微分方程
例
求解方程
d
y
d
x
=
−
x
y
(1.1.1)
\frac{dy}{dx}=-\frac{x}{y}\tag{1.1.1}
dxdy=−yx(1.1.1)
解
显然
y
≠
0
y\ne 0
y=0,则方程(1.1.1)可化为
y
d
y
=
−
x
d
x
ydy=-xdx
ydy=−xdx
两边求积分,得
∫
y
d
y
=
−
∫
x
d
x
\int ydy =-\int xdx
∫ydy=−∫xdx
求得通解
x
2
+
y
2
=
c
x^2+y^2=c
x2+y2=c
齐次方程
例
求解方程
d
y
d
x
=
y
x
+
tan
y
x
(1.2.1)
\frac{dy}{dx} = \frac{y}{x}+\tan\frac{y}{x}\tag{1.2.1}
dxdy=xy+tanxy(1.2.1)
解
令
u
=
y
x
u=\frac{y}{x}
u=xy
则方程(1.2.1)可化为
x
d
u
d
x
+
u
=
u
+
tan
u
x\frac{du}{dx}+u=u+\tan u
xdxdu+u=u+tanu
整理得到
d
u
tan
u
=
d
x
x
\frac{du}{\tan u}=\frac{dx}{x}
tanudu=xdx
两边积分,得到
ln
∣
sin
u
∣
=
ln
∣
x
∣
+
c
′
\ln|\sin u|=\ln|x|+c'
ln∣sinu∣=ln∣x∣+c′
整理得到
sin
u
=
c
x
(1.2.2)
\sin u=cx\tag{1.2.2}
sinu=cx(1.2.2)
另外,特解
u
=
0
u=0
u=0也被包含在方程(1.2.2)中。
综上,方程的通解为
sin
y
x
=
c
x
\sin \frac{y}{x} =cx
sinxy=cx
可化为齐次方程的类型
例
d
y
d
x
=
x
−
y
+
5
x
−
y
−
2
(1.3.1)
\frac{dy}{dx}=\frac{x-y+5}{x-y-2}\tag{1.3.1}
dxdy=x−y−2x−y+5(1.3.1)
解
令
u
=
x
−
y
u=x-y
u=x−y
则方程(1.3.1)可化为
1
−
d
u
d
x
=
u
+
5
u
−
2
1-\frac{du}{dx}=\frac{u+5}{u-2}
1−dxdu=u−2u+5
整理后得到
−
(
u
−
2
)
d
u
=
7
d
x
-(u-2)du=7dx
−(u−2)du=7dx
两边求积分,得到
−
u
2
+
4
u
=
14
x
+
c
-u^2+4u=14x+c
−u2+4u=14x+c
得到通解
(
x
−
y
)
2
+
4
y
+
10
x
=
c
(x-y)^2+4y+10x=c
(x−y)2+4y+10x=c
例
d
y
d
x
=
x
−
y
+
1
x
+
y
−
3
(1.3.2)
\frac{dy}{dx}=\frac{x-y+1}{x+y-3}\tag{1.3.2}
dxdy=x+y−3x−y+1(1.3.2)
解
求解方程组:
{
x
−
y
+
1
=
0
x
+
y
−
3
=
0
\left\{\begin{aligned} x-y+1=0\\ x+y-3=0 \end{aligned} \right.
{x−y+1=0x+y−3=0
得
x
=
1
,
y
=
2
x=1,y=2
x=1,y=2
令
{
x
=
X
+
1
y
=
Y
+
2
\left\{\begin{aligned} x = X+1\\ y = Y+2 \end{aligned} \right.
{x=X+1y=Y+2
则方程
(
1.3.2
)
(1.3.2)
(1.3.2)可化为:
d
Y
d
X
=
X
−
Y
X
+
Y
(1.3.3)
\frac{dY}{dX}=\frac{X-Y}{X+Y}\tag{1.3.3}
dXdY=X+YX−Y(1.3.3)
再令
u
=
Y
X
u=\frac{Y}{X}
u=XY
则方程
(
1.3.3
)
(1.3.3)
(1.3.3)可化为
d
X
X
=
1
+
u
1
−
2
u
−
u
2
d
u
\frac{dX}{X}=\frac{1+u}{1-2u-u^2}du
XdX=1−2u−u21+udu
两边求积分,整理,得:
X
2
(
1
−
2
u
−
u
2
)
=
c
′
X^2(1-2u-u^2)=c'
X2(1−2u−u2)=c′
把
X
,
u
X,u
X,u重新用
x
,
y
x,y
x,y表示,得:
x
2
−
y
2
−
2
x
y
+
2
x
+
6
y
=
c
x^2-y^2-2xy+2x+6y=c
x2−y2−2xy+2x+6y=c
线性方程与常数变易法
一阶线性齐次微分方程
即 可分离变量的微分方程
一阶线性非齐次微分方程
例
求解方程
(
x
+
1
)
d
y
d
x
−
n
y
=
e
x
(
x
+
1
)
n
+
1
(2.1.1)
(x+1)\frac{dy}{dx}-ny=e^x(x+1)^{n+1}\tag{2.1.1}
(x+1)dxdy−ny=ex(x+1)n+1(2.1.1)这里
n
n
n为常数.
解
%%%%%%%%%%%%方法一%%%%%%%%%%%
整理方程得到
d
y
d
x
=
n
x
+
1
y
+
e
x
(
x
+
1
)
n
\frac{dy}{dx}=\frac{n}{x+1}y+e^x(x+1)^n
dxdy=x+1ny+ex(x+1)n
首先,求解微分方程
d
y
d
x
=
n
x
+
1
y
(2.1.2)
\frac{dy}{dx}=\frac{n}{x+1}y\tag{2.1.2}
dxdy=x+1ny(2.1.2)
不难求得方程(2.1.2)的通解为:
y
=
c
1
(
x
+
1
)
n
y =c_1 (x+1)^n
y=c1(x+1)n
利用常数变易法
设
y
=
c
(
x
)
(
x
+
1
)
n
(2.1.3)
y = c(x)(x+1)^n\tag{2.1.3}
y=c(x)(x+1)n(2.1.3)
是方程(2.1.1)的通解
将方程(2.1.3)代入方程(2.1.1)中,化简得:
c
′
(
x
)
=
e
x
c'(x)=e^x
c′(x)=ex
所以
c
(
x
)
=
e
x
+
c
c(x)=e^x+c
c(x)=ex+c
所以方程(2.1.1)的通解为
y
=
(
e
x
+
c
)
(
x
+
1
)
n
y=(e^x+c)(x+1)^n
y=(ex+c)(x+1)n
%%%%%%%%%%%%方法二%%%%%%%%%%%
一阶线性非齐次微分方程
d
y
d
x
=
P
(
x
)
y
+
Q
(
x
)
\frac{dy}{dx}=P(x)y+Q(x)
dxdy=P(x)y+Q(x)的通解为:
y
=
e
∫
P
(
x
)
d
x
(
∫
Q
(
x
)
e
−
∫
P
(
x
)
d
x
d
x
+
c
)
(2.1.4)
y = e^{\int P(x)dx}(\int Q(x)e^{-\int P(x)dx}dx+c)\tag{2.1.4}
y=e∫P(x)dx(∫Q(x)e−∫P(x)dxdx+c)(2.1.4)
在本题中,
P
(
x
)
=
n
x
+
1
P(x)=\frac{n}{x+1}
P(x)=x+1n,
Q
(
x
)
=
e
x
(
x
+
1
)
n
Q(x)=e^x(x+1)^n
Q(x)=ex(x+1)n
代入公式,即可得到通解
y
=
(
e
x
+
c
)
(
x
+
1
)
n
y=(e^x+c)(x+1)^n
y=(ex+c)(x+1)n
伯努利方程
例
求解方程
d
y
d
x
=
6
y
x
−
x
y
2
(2.2.1)
\frac{dy}{dx}=6\frac{y}{x}-xy^2\tag{2.2.1}
dxdy=6xy−xy2(2.2.1)
解
情形一:
y
=
0
y=0
y=0,此时
y
=
0
y=0
y=0为方程的一个特解
情形二:
y
≠
0
y\ne 0
y=0
等式两边除以
y
2
y^2
y2,得到:
1
y
2
d
y
d
x
=
6
1
x
y
−
x
(2.2.2)
\frac{1}{y^2}\frac{dy}{dx} = 6\frac{1}{xy}-x\tag{2.2.2}
y21dxdy=6xy1−x(2.2.2)
令
z
=
1
y
z = \frac{1}{y}
z=y1
方程(2.2.2)可化为:
d
z
d
x
=
−
6
x
z
+
x
\frac{dz}{dx}=-\frac{6}{x}z+x
dxdz=−x6z+x
这是一个一阶线性非齐次方程,代入公式(2.1.4),得到:
z
=
x
2
8
+
c
x
6
(2.2.3)
z=\frac{x^2}{8}+\frac{c}{x^6}\tag{2.2.3}
z=8x2+x6c(2.2.3)
将
z
z
z用
y
y
y替换回来,得到:
1
z
=
x
2
8
+
c
x
6
\frac{1}{z}=\frac{x^2}{8}+\frac{c}{x^6}
z1=8x2+x6c
黎卡提方程
例
求解方程
y
′
=
y
2
−
x
2
+
1
(2.3.1)
y'=y^2-x^2+1\tag{2.3.1}
y′=y2−x2+1(2.3.1)
解
观察得到,
y
~
=
x
\tilde{y}=x
y~=x是方程的一个特解
设
y
=
x
+
u
y=x+u
y=x+u是方程的通解,将
y
y
y带入(2.3.1)中,得:
d
u
d
x
=
2
x
u
+
u
2
(2.3.2)
\frac{du}{dx}=2xu+u^2\tag{2.3.2}
dxdu=2xu+u2(2.3.2)
方程(2.3.2)是一个伯努利方程
情形一:
u
=
0
u=0
u=0,得到特解:
y
=
x
y=x
y=x
清醒二:
u
≠
0
u\ne 0
u=0:
方程(2.3.2)等式两端除以
u
2
u^2
u2,得到:
1
u
2
d
u
d
x
=
2
x
u
+
1
(2.3.3)
\frac{1}{u^2}\frac{du}{dx}=2\frac{x}{u}+1\tag{2.3.3}
u21dxdu=2ux+1(2.3.3)
令
w
=
1
u
w=\frac{1}{u}
w=u1,则方程(2.3.3)可化为:
d
w
d
x
=
−
2
x
w
−
1
(2.3.4)
\frac{dw}{dx}=-2xw-1\tag{2.3.4}
dxdw=−2xw−1(2.3.4)
方程(2.3.4)是一个一阶线性非齐次微分方程,带入公式(2.1.4),求得
w
=
e
−
x
2
(
c
−
∫
e
x
2
d
x
)
w=e^{-x^2}(c-\int e^{x^2}dx)
w=e−x2(c−∫ex2dx)
将
w
w
w用
x
,
y
x,y
x,y表示,得到:
y
=
x
+
e
x
2
c
−
∫
e
x
2
d
x
y=x+\frac{e^{x^2}}{c-\int e^{x^2}dx}
y=x+c−∫ex2dxex2
恰当方程与积分因子
恰当方程
不定积分法
例
求解方程
(
e
x
+
y
)
d
x
+
(
x
−
2
sin
y
)
d
y
=
0
(3.1.1)
(e^x+y)dx+(x-2\sin y)dy =0\tag{3.1.1}
(ex+y)dx+(x−2siny)dy=0(3.1.1)
解
设:
M
(
x
,
y
)
=
e
x
+
y
M(x,y)=e^x+y
M(x,y)=ex+y,
N
(
x
,
y
)
=
x
−
2
sin
y
N(x,y)=x-2\sin y
N(x,y)=x−2siny
因为
∂
M
(
x
,
y
)
∂
y
=
1
=
∂
N
(
x
,
y
)
∂
x
\frac{\partial M(x,y)}{\partial y} = 1 =\frac{\partial N(x,y)}{\partial x}
∂y∂M(x,y)=1=∂x∂N(x,y)
所以方程(3.1.1)是恰当方程
设方程的通解为
u
(
x
,
y
)
u(x,y)
u(x,y)
则
∂
u
∂
x
=
e
x
+
y
(3.1.2)
\frac{\partial u}{\partial x}=e^x+y \tag{3.1.2}
∂x∂u=ex+y(3.1.2)
∂
u
∂
y
=
x
−
2
sin
y
(3.1.3)
\frac{\partial u}{\partial y}=x-2\sin y \tag{3.1.3}
∂y∂u=x−2siny(3.1.3)
方程(3.1.2)对
x
x
x 求积分,得到:
u
(
x
,
y
)
=
∫
(
e
x
+
y
)
d
x
+
ψ
(
y
)
=
e
x
+
x
y
+
ψ
(
y
)
(3.1.4)
u(x,y) =\int(e^x+y)dx+\psi(y)=e^x+xy+\psi(y)\tag{3.1.4}
u(x,y)=∫(ex+y)dx+ψ(y)=ex+xy+ψ(y)(3.1.4)
将(3.1.4)代入(3.1.3)中,得到:
ψ
′
(
x
)
=
−
2
sin
y
(3.1.5)
\psi'(x)=-2\sin y \tag{3.1.5}
ψ′(x)=−2siny(3.1.5)
对(3.1.5)求积分,得到:
ψ
(
x
)
=
2
cos
y
(3.1.6)
\psi(x) = 2\cos y \tag{3.1.6}
ψ(x)=2cosy(3.1.6)
将(3.1.6)代入(3.1.4)中,得到:
u
(
x
,
y
)
=
e
x
+
x
y
+
2
cos
y
u(x,y)=e^x+xy+2\cos y
u(x,y)=ex+xy+2cosy
所以方程(3.1.1)的通解为:
e
x
+
x
y
+
2
cos
y
=
c
e^x+xy+2\cos y=c
ex+xy+2cosy=c
分组凑微分法
例
求解方程
3
x
2
d
x
+
4
y
3
d
y
+
6
x
y
2
d
x
+
6
x
2
y
d
y
(3.2.1)
3x^2dx + 4y^3dy+6xy^2dx+6x^2ydy\tag{3.2.1}
3x2dx+4y3dy+6xy2dx+6x2ydy(3.2.1)
解
不难看出方程(3.2.1)可化为:
d
(
x
3
+
y
4
)
+
d
(
3
x
2
y
2
)
d(x^3+y^4)+d(3x^2y^2)
d(x3+y4)+d(3x2y2)
即:
d
(
x
3
+
y
4
+
3
x
2
y
2
)
=
0
d(x^3+y^4+3x^2y^2)=0
d(x3+y4+3x2y2)=0
所以方程的通解为:
x
3
+
y
4
+
3
x
2
y
2
=
c
x^3+y^4+3x^2y^2=c
x3+y4+3x2y2=c
线积分法
例
(
y
cos
x
+
2
x
e
y
)
d
x
+
(
sin
x
+
x
2
e
y
+
2
)
d
y
=
0
(3.3.1)
(y\cos x+2xe^y)dx+(\sin x+x^2e^y+2)dy=0\tag{3.3.1}
(ycosx+2xey)dx+(sinx+x2ey+2)dy=0(3.3.1)
解
u
(
x
,
y
)
=
∫
(
0
,
0
)
(
x
,
y
)
(
y
cos
x
+
2
x
e
y
)
d
x
+
(
sin
x
+
x
2
e
y
+
2
)
d
y
=
∫
0
x
2
x
d
x
+
∫
0
y
(
sin
x
+
x
2
e
y
+
2
)
d
y
=
x
2
+
y
sin
x
+
(
x
2
−
1
)
e
y
+
2
y
=
y
sin
x
+
x
2
e
y
+
2
y
\begin{aligned} u(x,y)&=\int_{(0,0)}^{(x,y)}(y\cos x+2xe^y)dx+(\sin x+x^2e^y+2)dy\\ &=\int_0^x2xdx+\int_0^y(\sin x+x^2e^y+2)dy\\ &=x^2+y\sin x+(x^2-1)e^y+2y\\ &=y\sin x+x^2e^y+2y \end{aligned}
u(x,y)=∫(0,0)(x,y)(ycosx+2xey)dx+(sinx+x2ey+2)dy=∫0x2xdx+∫0y(sinx+x2ey+2)dy=x2+ysinx+(x2−1)ey+2y=ysinx+x2ey+2y
所以方程的通解为
c
=
y
sin
x
+
x
2
e
y
+
2
y
c=y\sin x+x^2e^y+2y
c=ysinx+x2ey+2y
积分因子
例
求解微分方程
(
y
2
2
+
2
y
e
x
)
d
x
+
(
y
+
e
x
)
d
y
=
0
(3.4.1)
(\frac{y^2}{2}+2ye^x)dx+(y+e^x)dy=0\tag{3.4.1}
(2y2+2yex)dx+(y+ex)dy=0(3.4.1)
其中
M
(
x
,
y
)
=
y
2
2
+
2
y
e
x
M(x,y)=\frac{y^2}{2}+2ye^x
M(x,y)=2y2+2yex,
N
(
x
,
y
)
=
y
+
e
x
N(x,y)=y+e^x
N(x,y)=y+ex
由于
∂
M
∂
y
−
∂
N
∂
x
N
=
y
+
e
x
y
+
e
x
=
1
=
ψ
(
x
)
\frac{\frac{\partial M}{\partial y}-\frac{\partial N}{\partial x}}{N}=\frac{y+e^x}{y+e^x}=1=\psi(x)
N∂y∂M−∂x∂N=y+exy+ex=1=ψ(x)
仅与
y
y
y无关,所以方程有一个只与
x
x
x有关的积分因子
μ
(
x
)
=
e
∫
ψ
(
x
)
d
x
=
e
x
\mu(x)=e^{\int\psi(x)dx}=e^x
μ(x)=e∫ψ(x)dx=ex
对方程(3.4.1)两边乘以
e
x
e^x
ex得到:
(
y
2
2
e
x
+
2
y
e
2
x
)
d
x
+
(
y
e
x
+
e
2
x
)
d
y
(3.4.2)
(\frac{y^2}{2}e^x+2ye^{2x})dx+(ye^x+e^{2x})dy\tag{3.4.2}
(2y2ex+2ye2x)dx+(yex+e2x)dy(3.4.2)
求解恰当方程(3.4.2),得到通解:
y
2
2
e
x
+
y
e
2
x
=
c
\frac{y^2}{2}e^x+ye^{2x}=c
2y2ex+ye2x=c