微分方程
约 2302 字大约 8 分钟
2025-11-21
微分方程(Differential Equation) 是含有未知函数及其导数的方程。微分方程是数学建模的重要工具,广泛应用于物理、工程、生物、经济等领域。
按阶数分类:
- 一阶微分方程:F(x,y,y′)=0
- 二阶微分方程:F(x,y,y′,y′′)=0
- 高阶微分方程:阶数 ≥3
按线性性分类:
- 线性微分方程:an(x)y(n)+⋯+a1(x)y′+a0(x)y=f(x)
- 非线性微分方程:不能表示为线性形式
一阶微分方程
1. 变量可分离方程
基本形式:dxdy=f(x)g(y)
解法:分离变量 g(y)dy=f(x)dx,然后两边积分
2. 一阶线性微分方程
标准形式:dxdy+P(x)y=Q(x)
解法:使用积分因子 u(x)=e∫P(x)dx
通解公式:
y=e−∫P(x)dx(∫Q(x)e∫P(x)dxdx+C)
3. 齐次微分方程
基本形式:dxdy=f(xy)
解法:令 u=xy,则 y=ux,dxdy=u+xdxdu
特殊类型的一阶微分方程
1. 变量可分离方程
基本形式:M(x)dx+N(y)dy=0
例题1:求微分方程 y′tanx=ylny 的通解。
解:
将方程改写为变量可分离形式:
dxdytanx=ylny
分离变量:
ylnydy=tanxdx=cotxdx
两边积分:
∫ylnydy=∫cotxdx
左边:令 u=lny,则 du=ydy:
∫udu=ln∣u∣+C1=ln∣lny∣+C1
右边:
∫cotxdx=ln∣sinx∣+C2
因此:
ln∣lny∣=ln∣sinx∣+C
化简得:
∣lny∣=C∣sinx∣⇒lny=Csinx
注意特殊情况:当 ylny=0 时,需要单独考虑:
- 若 y=1,则 lny=0,代入原方程验证可知 y≡1 也是解
2. 齐次微分方程
基本形式:dxdy=f(xy)
换元方法:令 u=xy,则 y=ux,y′=u+xu′
例题1
求解 (y+x2+y2)dx=xdy
解:
将方程化为标准形式:
dxdy=xy+x2+y2=xy+1+(xy)2
令 u=xy,则 y=ux,dxdy=u+xdxdu:
u+xdxdu=u+1+u2
简化得:
xdxdu=1+u2
分离变量:
1+u2du=xdx
两边积分:
∫1+u2du=∫xdx
左边积分:ln∣u+1+u2∣+C1 右边积分:ln∣x∣+C2
因此:
ln∣u+1+u2∣=ln∣x∣+C
化简:
u+1+u2=Cx
将 u=xy 代回:
xy+1+(xy)2=Cx
例题2
求解 (1−e−yx)ydx+(y−x)dy=0
解:
令 u=yx,则 x=uy,dx=udy+ydu
代入原方程:
(1−e−u)y(udy+ydu)+(y−uy)dy=0
展开并整理:
(1−e−u)y2du+[(1−e−u)uy+(1−u)y]dy=0
化简:
y2(1−e−u)du+y[1−ue−u]dy=0
分离变量:
yy2[1−e−u]du+[1−ue−u]dy=0
即:
y(1−e−u)du+(1−ue−u)dy=0
重新排列:
(1−e−u)du=yue−u−1dy
积分后得到:
u−e−u=yC
将 u=yx 代回:
yx−e−yx=yC
二阶微分方程
1. 常系数线性微分方程
齐次方程:y′′+ay′+by=0
特征方程:r2+ar+b=0
解的情况:
- 两个不同实根 r1,r2:y=C1er1x+C2er2x
- 重根 r:y=(C1+C2x)erx
- 复根 α±βi:y=eαx(C1cosβx+C2sinβx)
非齐次方程:y′′+ay′+by=f(x)
通解:y=yh+yp,其中 yh 是齐次解,yp 是特解
2. 可降阶的二阶微分方程
类型1:y′′=f(x)
- 解法:直接积分两次
类型2:y′′=f(x,y′)(不显含 y)
- 解法:令 p=y′,则 y′′=p′,化为一阶方程
类型3:y′′=f(y,y′)(不显含 x)
- 解法:令 p=y′,则 y′′=pdydp
可降阶的高阶微分方程
例题1
求微分方程 yy′′−y′2=0 的通解及满足 y(0)=1,y′(0)=21 的特解。
解:
令 p=y′,则 y′′=dxdp=dydpdxdy=pdydp
代入原方程:
y⋅pdydp−p2=0
整理得:
ypdydp=p2
当 p=0 时:
p1dp=y1dy
两边积分:
ln∣p∣=ln∣y∣+C1
因此:
p=C2y
代回 p=y′:
dxdy=C2y
分离变量:
ydy=C2dx
积分得:
ln∣y∣=C2x+C3
因此通解为:
y=CeC2x
利用初始条件 y(0)=1,y′(0)=21:
C=1,C2=21
所以特解为:
y=e2x
全微分方程
基本理论:
全微分方程:P(x,y)dx+Q(x,y)dy=0
充要条件:∂y∂P=∂x∂Q
解法:若存在函数 u(x,y) 使得 du=Pdx+Qdy,则 u(x,y)=C
例题1
求微分方程 (4−x+y)dx+(x+y−2)dy=0 的通解
解:
设 P(x,y)=4−x+y,Q(x,y)=x+y−2
验证是否为全微分方程:
∂y∂P=1,∂x∂Q=1
因为 ∂y∂P=∂x∂Q,所以是全微分方程。
求 u(x,y):
u(x,y)=∫P(x,y)dx+ϕ(y)=∫(4−x+y)dx+ϕ(y)=4x−2x2+xy+ϕ(y)
对 y 求导:
∂y∂u=x+ϕ′(y)
与 Q(x,y) 比较:
x+ϕ′(y)=x+y−2⟹ϕ′(y)=y−2
积分得:
ϕ(y)=2y2−2y+C
因此:
u(x,y)=4x−2x2+xy+2y2−2y+C
通解为:
4x−2x2+xy+2y2−2y=C
例题2
求微分方程 (y2+1)dx=y(y−2x)dy 的通解
解:
将方程重写为标准形式:
(y2+1)dx+(2x−y)ydy=0
设 P(x,y)=y2+1,Q(x,y)=(2x−y)y=2xy−y2
验证:
∂y∂P=2y,∂x∂Q=2y
是全微分方程。
求 u(x,y):
u(x,y)=∫(y2+1)dx+ϕ(y)=xy2+x+ϕ(y)
对 y 求导:
∂y∂u=2xy+ϕ′(y)
与 Q(x,y) 比较:
2xy+ϕ′(y)=2xy−y2⟹ϕ′(y)=−y2
积分得:
ϕ(y)=−3y3+C
因此通解为:
xy2+x−3y3=C
可化为齐次方程的微分方程
形式:(a1x+b1y+c1)dx+(a2x+b2y+c2)dy=0
解法:通过坐标变换化为齐次方程
例题1
求微分方程 (x−y+3)y′=(x+y−1) 的通解
解:
将方程化为标准形式:
dxdy=x−y+3x+y−1
求交点:
{x+y−1=0x−y+3=0⟹{x=−1y=2
作坐标变换:
{u=x+1v=y−2⟹{x=u−1y=v+2
则 dxdy=dudv,方程变为:
dudv=(u−1)−(v+2)+3(u−1)+(v+2)−1=u−vu+v
这是齐次方程,令 w=uv,则 v=uw,dudv=w+ududw:
w+ududw=u−uwu+uw=1−w1+w
因此:
ududw=1−w1+w−w=1−w1+w−w+w2=1−w1+w2
分离变量:
1+w21−wdw=udu
积分:
∫1+w21−wdw=∫udu
左边:
∫1+w21dw−∫1+w2wdw=arctanw−21ln(1+w2)
右边:ln∣u∣+C
因此:
arctanw−21ln(1+w2)=ln∣u∣+C
代回原变量:
arctanx+1y−2−21ln[1+(x+1y−2)2]=ln∣x+1∣+C
伯努利方程与特殊技巧
1. 反向代理技巧
思想:将 y 视为自变量,x 视为因变量
例题1
求微分方程 dxdy=x4+y22x3y 满足 y(1)=1 的特解
解:
将方程反向:
dydx=2x3yx4+y2=2yx+2x3y
整理:
dydx−2y1x=2yx−3
这是关于 x 的伯努利方程。令 z=x4,则 dydz=4x3dydx:
4x31dydz−2y1x=2yx−3
两边乘以 4x3:
dydz−2yx4=2y
即:
dydz−2yz=2y
这是一阶线性微分方程,积分因子为 e−y2:
dyd(ze−y2)=2ye−y2
积分得:
ze−y2=−e−y2+C
因此:
z=−1+Cey2
代回 z=x4:
x4=−1+Cey2
利用初始条件 y(1)=1:
1=−1+Ce1⟹C=e2
特解为:
x4=−1+e2ey2=−1+2ey2−1
2. 参数变换技巧
例题1
求微分方程 y′′+(x+siny)y′3=0 的通解
解:
利用 x′=y′1(当 y′=0 时):
y′′=dxd(x′1)=−(x′)3x′′
代入原方程:
−(x′)3x′′+(x+siny)(x′)31=0
乘以 (x′)3:
−x′′+x+siny=0⟹x′′−x=siny
这是关于 x 的二阶线性非齐次方程,但要注意 siny 中的 y 是 x 的函数。
设通解为 x=xˉ+x∗,其中 xˉ 是齐次解,x∗ 是特解。
齐次方程:x′′−x=0,通解为 xˉ=C1ey+C2e−y
实际上,这个问题需要更深入的技巧,通常在这种情况下,需要寻找特殊的方法或数值解。
