一阶线性微分方程的初等积分法例题

变量分离方程与变量变换

可分离变量的微分方程


求解方程
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' lnsinu=lnx+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=xy2xy+5(1.3.1)


u = x − y u=x-y u=xy
则方程(1.3.1)可化为
1 − d u d x = u + 5 u − 2 1-\frac{du}{dx}=\frac{u+5}{u-2} 1dxdu=u2u+5
整理后得到
− ( u − 2 ) d u = 7 d x -(u-2)du=7dx (u2)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 (xy)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+y3xy+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. {xy+1=0x+y3=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+YXY(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=12uu21+udu
两边求积分,整理,得:
X 2 ( 1 − 2 u − u 2 ) = c ′ X^2(1-2u-u^2)=c' X2(12uu2)=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 x2y22xy+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)dxdyny=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=eP(x)dx(Q(x)eP(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=6xyxy2(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=6xy1x(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=y2x2+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=2xw1(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=ex2(cex2dx)
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+cex2dxex2

恰当方程与积分因子

恰当方程

不定积分法


求解方程
( 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+(x2siny)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)=x2siny
因为
∂ M ( x , y ) ∂ y = 1 = ∂ N ( x , y ) ∂ x \frac{\partial M(x,y)}{\partial y} = 1 =\frac{\partial N(x,y)}{\partial x} yM(x,y)=1=xN(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} xu=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} yu=x2siny(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+(x21)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) NyMxN=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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值