常见一阶微分方程@可降阶的二阶微分方程@Bernoulli方程
可用变量代换法求解的一阶微分方程
可分离变量的方程
-
d y d x = f ( x ) g ( y ) g ( y ) ≠ 0 \frac{\mathrm{d}y}{\mathrm{d}x}=f(x)g(y)\quad g(y)\neq{0} dxdy=f(x)g(y)g(y)=0可以分离变量为 d y g ( y ) = f ( x ) d x \frac{\mathrm{d}y}{g(y)}=f(x)\mathrm{d}x g(y)dy=f(x)dx
-
两边积分: ∫ d y g ( y ) = ∫ f ( x ) d x + C \int{\frac{\mathrm{d}y}{g(y)}}=\int{f(x)\mathrm{d}x}+C ∫g(y)dy=∫f(x)dx+C
齐次方程
- 形如或可化为
d
y
d
x
=
f
(
y
x
)
\frac{\mathrm{d}y}{\mathrm{d}x}=f(\frac{y}{x})
dxdy=f(xy),
(
x
≠
0
)
(x\neq{0})
(x=0)
(0)
的微分方程称为齐次方程- 令函数
u
=
y
x
u=\frac{y}{x}
u=xy
(1)
,则 y = x u y=xu y=xu,(2)
对其两侧求关于 x x x的导数: d d x y \frac{\mathrm{d}}{\mathrm{d}x}y dxdy= d d x ( x u ) \frac{\mathrm{d}}{\mathrm{d}x}(xu) dxd(xu),即 d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= u + x d u d x {u}+x{\frac{\mathrm{d}u}{\mathrm{d}x}} u+xdxdu(3)
- 其中式(1)用一个新的未知函数 u u u代替了原未知函数 y y y,
- 由于 y y y是关于 x x x得函数,所以由式(1)可知,函数 u u u也是关于 x x x的函数
- 对(1)变形得到(2)式,再对(2)式两边求导,得到(3)式
- 将(1)代入(0),得
d
y
d
x
=
f
(
u
)
\frac{\mathrm{d}y}{\mathrm{d}x}=f(u)
dxdy=f(u)
(4)
,将(4)代入(3),得 f ( u ) = u + x d u d x f(u)=u+x\frac{\mathrm{d}u}{\mathrm{d}x} f(u)=u+xdxdu(5)
,这就将原微分方程中的 y y y用 u u u代替 - 再对式(5)进行分离变量,得
f
(
u
)
−
u
d
u
=
x
d
x
\frac{f(u)-u}{\mathrm{d}u}=\frac{x}{\mathrm{d}x}
duf(u)−u=dxx
(6)
,两边取倒数,得 d u f ( u ) − u = d x x \frac{\mathrm{d}u}{f(u)-u}=\frac{\mathrm{d}x}{x} f(u)−udu=xdx,两边积分: ∫ d u f ( u ) − u \int\frac{\mathrm{d}u}{f(u)-u} ∫f(u)−udu= ∫ d x x + C \int{\frac{\mathrm{d}x}{x}}+C ∫xdx+C= ln ∣ x ∣ + C \ln|x|+C ln∣x∣+C(7)
- 设(7)左边的积分结果为 Φ ( u ) \Phi(u) Φ(u),代回式(1),得微分方程(0)得通解 Φ ( y x ) \Phi(\frac{y}{x}) Φ(xy)= ln ∣ x ∣ + C \ln|x|+C ln∣x∣+C
- 令函数
u
=
y
x
u=\frac{y}{x}
u=xy
可齐次化的方程
- 形如
d
y
d
x
\frac{\mathrm{d}y}{\mathrm{d}x}
dxdy=
f
(
a
x
+
b
y
+
c
a
1
x
+
b
1
y
+
c
1
)
f(\frac{ax+by+c}{a_1x+b_1y+c_1})
f(a1x+b1y+c1ax+by+c)
(0)
的方程可以通过变换化为齐次方程 - 分三种情形讨论
- 当 c = c 1 = 0 c=c_1=0 c=c1=0,方程(0)就是齐次方程
- 当
a
a
1
=
b
b
1
\frac{a}{a_1}=\frac{b}{b_1}
a1a=b1b时,设此比值为
λ
\lambda
λ,则方程(0)可以表示为
d
y
d
x
\frac{\mathrm{d}y}{\mathrm{d}x}
dxdy=
f
(
λ
(
a
1
x
+
b
1
y
)
+
c
a
1
x
+
b
1
y
+
c
1
)
f(\frac{\lambda(a_1x+b_1y)+c}{a_1x+b_1y+c_1})
f(a1x+b1y+c1λ(a1x+b1y)+c)
(0-1)
- 引入新变量(函数)
u
=
a
1
x
+
b
1
y
u=a_1x+b_1y
u=a1x+b1y
(1)
,对其两边关于 x x x求导: d u d x \frac{\mathrm{d}u}{\mathrm{d}x} dxdu= a 1 + b 1 d y d x a_1+b_1\frac{\mathrm{d}y}{\mathrm{d}x} a1+b1dxdy(2)
,变形可得 d y d x \frac{\mathrm{d}y}{\mathrm{d}{x}} dxdy= 1 b 1 ( d u d x − a 1 ) \frac{1}{b_1}(\frac{\mathrm{d}u}{\mathrm{d}x}-a_1) b11(dxdu−a1)(2-1)
- 将(1)代入(0-1),得
d
y
d
x
\frac{\mathrm{d}y}{\mathrm{d}x}
dxdy=
f
(
λ
u
+
c
u
+
c
1
)
f(\frac{\lambda{u}+c}{u+c_1})
f(u+c1λu+c)
(3)
;比较(2-1),(3),得 1 b 1 ( d u d x − a 1 ) \frac{1}{b_1}(\frac{\mathrm{d}u}{\mathrm{d}x}-a_1) b11(dxdu−a1)= f ( λ u + c u + c 1 ) f(\frac{\lambda{u}+c}{u+c_1}) f(u+c1λu+c)(4)
,这是可分离变量( u , x u,x u,x)的方程,问题得到转化- 将(4)分离变量,可得 b 1 f ( λ u + c u + c 1 ) + a 1 b_1f(\frac{\lambda{u}+c}{u+c_1})+a_1 b1f(u+c1λu+c)+a1= d u d x \frac{\mathrm{d}u}{\mathrm{d}x} dxdu,即 d x \mathrm{d}x dx= d u b 1 f ( λ u + c u + c 1 ) + a 1 \frac{\mathrm{d}u}{b_1f(\frac{\lambda{u}+c}{u+c_1})+a_1} b1f(u+c1λu+c)+a1du
- 引入新变量(函数)
u
=
a
1
x
+
b
1
y
u=a_1x+b_1y
u=a1x+b1y
- 当
a
a
1
≠
b
b
1
\frac{a}{a_1}\neq{\frac{b}{b_1}}
a1a=b1b,且
c
,
c
1
≠
0
c,c_1\neq{0}
c,c1=0
- 借助平面直线方程性质来求解
- 方程
a
x
+
b
y
+
c
=
0
ax+by+c=0
ax+by+c=0;
a
1
x
+
b
1
y
+
c
=
0
a_1x+b_1y+c=0
a1x+b1y+c=0表示平面上的两条直线
l
1
,
l
2
l_1,l_2
l1,l2
- 设两直线交点为 P ( α , β ) P(\alpha,\beta) P(α,β),即 α , β \alpha,\beta α,β满足 a α + b β + c = 0 a\alpha+b\beta+c=0 aα+bβ+c=0, a 1 α + b 1 β + c 1 = 0 a_1\alpha+b_1\beta+c_1=0 a1α+b1β+c1=0
- 由直线的一般式方程的性质可知, c = 0 c=0 c=0表示直线过原点(0,0)
- 显然 l 1 , l 2 l_1,l_2 l1,l2的均不过原点,其交点 P P P不可能为原点,即 α , β \alpha,\beta α,β不同时为0( α ≠ 0 \alpha\neq{0} α=0或 β ≠ 0 \beta\neq{0} β=0),否则就是情形(1)的直线方程表示
- 为了解决情形3,参考情形1,可以将坐标系( x O y xOy xOy)平移到 P P P,这就将问题转化为情形(1)
- 而平移后的坐标记为
X
O
′
Y
XO'Y
XO′Y其中
O
′
O'
O′就是
P
P
P,即新坐标为
X
P
Y
XPY
XPY
- 同一位置的点的新旧坐标转换公式
(5)
: X = x − α X=x-\alpha X=x−α; Y = y − β Y=y-\beta Y=y−β- 或表示为 x = X + α x=X+\alpha x=X+α; y = Y + β y=Y+\beta y=Y+β
- l 1 , l 2 l_1,l_2 l1,l2分别表示为 a X + b Y = 0 aX+bY=0 aX+bY=0, a 1 X + b 1 Y = 0 a_1X+b_1Y=0 a1X+b1Y=0;
- 在新坐标系上,方程(0)表示为
d
Y
d
X
\frac{\mathrm{d}Y}{\mathrm{d}X}
dXdY=
f
(
a
X
+
b
Y
a
1
X
+
b
1
Y
)
f(\frac{aX+bY}{a_1X+b_1Y})
f(a1X+b1YaX+bY)
(6)
- 这就将问题转化为情形(1)(齐次化)
- 同一位置的点的新旧坐标转换公式
- 解出(6)的通解后,将(5)代入得到方程(0)的通解
Bernoulli方程
- 方程
y
′
+
p
(
x
)
y
=
q
(
x
)
y
n
y'+p(x)y=q(x)y^{n}
y′+p(x)y=q(x)yn
(1)
, ( n ≠ 0 , n ≠ 1 ) (n\neq{0},n\neq{1}) (n=0,n=1)称为Bernoulli方程
伯努利方程求解步骤
- 原式化为
y
−
n
d
y
d
x
+
p
(
x
)
y
1
−
n
=
q
(
x
)
y^{-n}\frac{\mathrm{d}y}{\mathrm{d}x}+p(x)y^{1-n}=q(x)
y−ndxdy+p(x)y1−n=q(x)
(2)
,确定出其中的 p ( x ) , q ( x ) , n p(x),q(x),n p(x),q(x),n三要素 - 令
z
=
y
1
−
n
z=y^{1-n}
z=y1−n
(3)
,列出关于函数 z ( x ) z(x) z(x)的一阶线性微分方程(标准形式): z ′ + ( 1 − n ) p ( x ) z = ( 1 − n ) q ( x ) z'+(1-n)p(x)z=(1-n)q(x) z′+(1−n)p(x)z=(1−n)q(x)(3-0)
- 代入线性微分方程的通解公式,再代回
y
y
y,得出
y
=
y
(
x
)
y=y(x)
y=y(x)
(3-1)
(得原微分方程的通解
分析和推导
- 当 n = 0 n=0 n=0或 n = 1 n=1 n=1时,方程(1)退化为普通线性微分方程;
- 当
n
≠
0
,
n
≠
1
n\neq{0},n\neq{1}
n=0,n=1时,此方程虽然不是线性的,但可以通过变量代换,化为线性方程
- 通常以 y n y^{n} yn除方程(1)两边,得方程(2)
- 式(2)左端第一项
y
−
n
d
y
d
x
y^{-n}\frac{\mathrm{d}y}{\mathrm{d}x}
y−ndxdy
(3-2)
与 y 1 − n y^{1-n} y1−n的导数 d y 1 − n d x \frac{\mathrm{d}y^{1-n}}{\mathrm{d}x} dxdy1−n= ( 1 − n ) y − n ⋅ d y d x (1-n)y^{-n}\cdot\frac{\mathrm{d}y}{\mathrm{d}x} (1−n)y−n⋅dxdy(4)
仅差一个常数因子 ( 1 − n ) (1-n) (1−n),变形式(4),可得 y − n d y d x y^{-n}\frac{\mathrm{d}y}{\mathrm{d}x} y−ndxdy= 1 1 − n d y 1 − n d x \frac{1}{1-n}{\frac{\mathrm{d}y^{1-n}}{\mathrm{d}x}} 1−n1dxdy1−n(4-1)
- Note: y y y是 x x x函数, y 1 − n y^{1-n} y1−n对 x x x求导是复合函数求导: ( y 1 − n ) x ′ (y^{1-n})_{x}' (y1−n)x′= ( 1 − n ) y − n ⋅ y ′ (1-n)y^{-n}\cdot{y'} (1−n)y−n⋅y′
- 因此引入新的因变量
z
=
y
1
−
n
z=y^{1-n}
z=y1−n,即式(3),从而式(4-1)改写为
y
−
n
d
y
d
x
y^{-n}\frac{\mathrm{d}y}{\mathrm{d}x}
y−ndxdy=
1
1
−
n
d
z
d
x
\frac{1}{1-n}{\frac{\mathrm{d}z}{\mathrm{d}x}}
1−n1dxdz
(4-2)
- 由式(4-2),式(2)写成: 1 1 − n d z d x + p ( x ) z \frac{1}{1-n}\frac{\mathrm{d}z}{\mathrm{d}x}+p(x)z 1−n1dxdz+p(x)z= q ( x ) q(x) q(x)
- 即
d
z
d
x
{\frac{\mathrm{d}z}{\mathrm{d}x}}
dxdz+
(
1
−
n
)
p
(
x
)
z
(1-n)p(x)z
(1−n)p(x)z=
(
1
−
n
)
q
(
x
)
(1-n)q(x)
(1−n)q(x)
(5)
- 方程(5)显然是一阶线性方程,这就将Bernoulli方程转换为一阶线性方程问题
例
- 求一阶常系数微分方程
y
′
−
y
+
2
x
y
y'-y+2\frac{x}{y}
y′−y+2yx=
0
0
0
(1)
- 解
- 这是一个伯努利方程,将其化为标准形:
y
′
−
y
y'-y
y′−y=
−
2
x
⋅
y
−
1
-2x\cdot{y^{-1}}
−2x⋅y−1
(2)
- 令 p ( x ) p(x) p(x)= − 1 -1 −1; q ( x ) = − 2 x q(x)=-2x q(x)=−2x; n = − 1 n=-1 n=−1,即 1 − n = 2 1-n=2 1−n=2
- 令
z
=
y
1
−
n
z=y^{1-n}
z=y1−n=
y
2
y^2
y2
(3)
- 由Bernoulli方程的解公式得
z
z
z的一阶线性微分方程:
z
′
+
2
(
−
1
)
z
=
2
(
−
2
x
)
z'+2(-1)z=2(-2x)
z′+2(−1)z=2(−2x),即
z
′
−
2
z
=
−
4
x
z'-2z=-4x
z′−2z=−4x
(4)
- 代入一阶线性微分方程解公式
- g ( x ) g(x) g(x)= ∫ − 2 d x \int{-2}\mathrm{d}x ∫−2dx= − 2 x -2x −2x
- 从而
z
z
z=
e
2
x
(
C
+
∫
−
4
x
⋅
e
−
2
x
d
x
)
e^{2x}(C+\int{-4x}\cdot{e^{-2x}}\mathrm{d}x)
e2x(C+∫−4x⋅e−2xdx)=
C
(
e
2
x
+
2
x
+
1
)
C(e^{2x}+2x+1)
C(e2x+2x+1)
(5)
,由分部积分可以算出该结果
- 回代(3)到(5),得
y
2
y^2
y2=
C
(
e
2
x
+
2
x
+
1
)
C(e^{2x}+2x+1)
C(e2x+2x+1),
(6)
这就是原方程的隐式通解
- 这是一个伯努利方程,将其化为标准形:
y
′
−
y
y'-y
y′−y=
−
2
x
⋅
y
−
1
-2x\cdot{y^{-1}}
−2x⋅y−1
可降阶的二阶微分方程
- 二阶及以上的微分方程称为高阶微分方程
- 某些高阶微分方程可以通过代换法转化为较低阶的方程
类型0
-
y
′
′
=
f
(
x
)
y''=f(x)
y′′=f(x)
(0)
型的微分方程- 方程右端仅为 x x x的函数,作2次积分即可
- 第1次:
y
′
y'
y′=
∫
f
(
x
)
d
x
+
C
1
\int{f(x)\mathrm{d}x}+C_1
∫f(x)dx+C1
(1)
- 第2次:
y
y
y=
∫
(
∫
f
(
x
)
d
x
)
d
x
+
C
1
x
+
C
2
\int(\int{f(x)}\mathrm{d}x)\mathrm{d}x+C_1x+C_2
∫(∫f(x)dx)dx+C1x+C2
(2)
- Note:直接作2次积分:
∫
(
∫
f
(
x
)
d
x
)
d
x
+
C
\int{(\int{f(x)\mathrm{d}x})\mathrm{d}x}+C
∫(∫f(x)dx)dx+C会少一个自由变量
- 因此,要们分开作2次积分,要么直接用公式(2)
类型1
-
y
′
′
=
f
(
x
,
y
′
)
y''=f(x,y')
y′′=f(x,y′)
(1)
型的微分方程- 该方程右端不含有未知函数 y y y
- 令 y ′ = p y'=p y′=p( y y y的一阶导数代换为某个函数 p p p)
- 那么
y
′
′
=
d
p
d
x
=
p
′
y''=\frac{\mathrm{d}p}{\mathrm{d}x}=p'
y′′=dxdp=p′;而方程(1)改写成
p
′
=
f
(
x
,
p
)
p'=f(x,p)
p′=f(x,p)
(2)
,这就将二阶方程以代换的方法降阶为一阶方程 - 设(2)的通解为
p
=
ϕ
(
x
,
C
1
)
p=\phi(x,C_1)
p=ϕ(x,C1),即
d
y
d
x
\frac{\mathrm{d}y}{\mathrm{d}x}
dxdy=
ϕ
(
x
,
C
1
)
\phi(x,C_1)
ϕ(x,C1)
(3)
,这是一个关于函数 y y y的一阶微分方程,继续解这个方程 - 对(3)两边积分,得通解:
y
=
∫
ϕ
(
x
,
C
1
)
d
x
+
C
2
y=\int\phi(x,C_1)\mathrm{d}x+C_2
y=∫ϕ(x,C1)dx+C2
(4)
类型2
-
y
′
′
=
f
(
y
,
y
′
)
y''=f(y,y')
y′′=f(y,y′)
(1)
型的微分方程- 方程右端不明显得含有自变量 x x x
- 同样令
y
′
=
p
y'=p
y′=p,利用复合函数求导法,将
y
′
′
y''
y′′表示为
p
p
p关于
y
y
y的导数:
d
p
d
y
\frac{\mathrm{d}p}{\mathrm{d}y}
dydp(微商的角度描述导数)
-
y
′
′
y''
y′′=
d
p
d
x
\frac{\mathrm{d}p}{\mathrm{d}x}
dxdp=
d
p
d
x
\frac{\mathrm{d}p}{\mathrm{d}x}
dxdp
d
y
d
y
\frac{\mathrm{d}y}{\mathrm{d}y}
dydy=
d
p
d
y
\frac{\mathrm{d}p}{\mathrm{d}y}
dydp
d
y
d
x
\frac{\mathrm{d}y}{\mathrm{d}x}
dxdy=
d
p
d
y
p
\frac{\mathrm{d}p}{\mathrm{d}y}p
dydpp
(2)
,即公式 y ′ ′ y'' y′′= d p d y p \frac{\mathrm{d}p}{\mathrm{d}y}p dydpp,此公式是核心,配项手法是重点
-
y
′
′
y''
y′′=
d
p
d
x
\frac{\mathrm{d}p}{\mathrm{d}x}
dxdp=
d
p
d
x
\frac{\mathrm{d}p}{\mathrm{d}x}
dxdp
d
y
d
y
\frac{\mathrm{d}y}{\mathrm{d}y}
dydy=
d
p
d
y
\frac{\mathrm{d}p}{\mathrm{d}y}
dydp
d
y
d
x
\frac{\mathrm{d}y}{\mathrm{d}x}
dxdy=
d
p
d
y
p
\frac{\mathrm{d}p}{\mathrm{d}y}p
dydpp
- 从而方程(1)化为
p
d
p
d
y
p\frac{\mathrm{d}p}{\mathrm{d}y}
pdydp=
f
(
y
,
p
)
f(y,p)
f(y,p)
(3)
,这是一个一阶微分方程 - 设(3)的通解为
y
′
=
p
=
ϕ
(
y
,
C
1
)
y'=p=\phi(y,C_1)
y′=p=ϕ(y,C1)
(4)
,即 d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= ϕ ( y , C 1 ) \phi(y,C_1) ϕ(y,C1)分离变量得 d y ϕ ( y , C 1 ) = d x \frac{\mathrm{d}y}{\phi(y,C_1)}=\mathrm{d}x ϕ(y,C1)dy=dx(5)
, - 两边积分,得方程的 ∫ d y ϕ ( y , C 1 ) \int{\frac{\mathrm{d}y}{\phi(y,C_1)}} ∫ϕ(y,C1)dy= x + C 2 x+C_2 x+C2,可求得通解
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
2022-12-17 math_基本导数公式@积分公式@部分推导@数学工具@公式识别@计算推演工具
2022-12-17 算法导论答案网Solutions to Introduction to Algorithms Third Edition_开源免费完整