微积分-常微分方程解题方法
北京理工大学
微积分-常微分方程解法
常微分方程各种解题方法
程功 2011/2/16
1.几个基本定义
(1)微分方程:凡含有未知函数的导数或微分的方程叫微分方程.
实质: 联系自变量, 未知函数以及未知函数的某些导数(或微分) 之间的关系式.
分类1: 常微分方程: 未知函数为一元函数 偏微分方程: 未知函数为多元函数
分类2:
微分方程的阶: 微分方程中出现的未知函数的最高阶导数的阶数称之. 一阶微分方程F (x , y , y ') =0, y '=f (x , y );
高阶(n )微分方程F (x , y , y ', , y (n ) ) =0, y (n ) =f (x , y , y ', , y (n -1) ).
分类3: 线性与非线性微分方程. y '+P (x ) y =Q (x ), x (y ') 2-2yy '+x =0;
⎧dy
=3y -2z , ⎪⎪dx
分类4: 单个微分方程与微分方程组. ⎨
⎪dz =2y -z , ⎪⎩dx
(2)微分方程的解:代入微分方程能使方程成为恒等式的函数称之.
微分方程的解的分类:
① 通解: 微分方程的解中含有任意常数, 且任意常数的个数与微分方程的阶数相同.
例y '=y , 通解y =Ce x ;
y ''+y =0, 通解y =C 1sin x +C 2cos x ;
② 特解: 确定了通解中任意常数以后的解. (3)解的图象: 微分方程的积分曲线. 通解的图象: 积分曲线族.
(4)初始条件: 用来确定任意常数的条件. 初值问题: 求微分方程满足初始条件的解的问题.
⎧⎪y '=f (x , y ) 一阶:⎨过定点的积分曲线;
y =y 0⎪⎩x =x 0
⎧⎪y ''=f (x , y , y ')
二阶:⎨过定点且在定点的切线的斜率为定值的积分曲线.
'y =y , y =y 00x =x 0⎪⎩x =x 0
2. 可分离变量的微分方程
-dy 25
可分离变量微分方程的形式g (y ) dy =f (x ) dx 例如=2x y ⇒y 5dy =2x 2dx ,
dx
4
4
解法:设函数g (y ) 和f (x ) 是连续的, ⎰g (y ) dy =⎰f (x ) dx 设函数G (y ) 和F (x ) 是依次为g (y ) 和f (x ) 的原函数, G (y ) =F (x ) +C 为微分方程的解.
3. 齐次方程
dy y
=f () 的微分方程称为齐次方程. 形如dx x
y dy du
=u +x , 解法:作变量代换u =, 即y =xu , ∴
x dx dx
du du f (u ) -u =f (u ), 即=. (可分离变量的方程) 代入原式u +x dx dx x
(1)当f (u ) -u ≠0时, 得⎰
du
=ln C 1x ,
f (u ) -u
y
即x =Ce , (ϕ(u ) =⎰
ϕu )
ϕ() du y
将u =代入, 得通解x =Ce x , )
x f (u ) -u
(2)当∃u 0, 使f (u 0) -u 0=0, 则u =u 0是新方程的解, 代回原方程,
得齐次方程的解y =u 0x . 4.可化为齐次的方程 定义形如
dy ax +by +c
=f () 的微分方程 dx a 1x +b 1y +c 1
当c =c 1=0时, 为齐次方程. 否则为非齐次方程. 解法:令x =X +h ,
y =Y +k ,(其中h 和k 是待定的常数)dx =dX , dy =dY
dY aX +bY +ah +bk +c ⎧ah +bk +c =0,
=f () ⎨
dX a 1X +b 1Y +a 1h +b 1k +c 1⎩a 1h +b 1k +c 1=0,
(1)
a 1b 1⎧X =x -h , dY aX +bY
≠有唯一一组解. =f () 得通解代回⎨a 2b 2dX a 1X +b 1Y ⎩Y =y -k ,a 1b 1dy ax +by +c ==λ, 方程可化为=f (), 令z =ax +by , a b dx λ(ax +by ) +c 1
(2)
则
dz dy 1dz z +c =a +b (-a ) =f (). 可分离变量. dx dx b dx λz +c 1
5.其它类型:通过变量代换化为可分离变量方程
(1)f (x ±y )(dx ±dy ) =g (x ) dx 令u =x ±y , du =dx ±dy , 方程化为f (u ) du =g (x ) dx (2)f (xy )(xdy +ydx ) =g (x ) dx 令u =xy , du =xdy +ydx , 代入方程得f (u ) du =g (x ) dx
y y xdy -ydx g (x ) (3)f ()(xdy -ydx ) =g (x ) dx 令u =, 则du =, f (u ) du =dx 代入方程得22
x x x x
(4)f (x 2+y 2)(xdx +ydy ) =g (x ) dx 令u =x 2+y 2, 则du =2xdx +2ydy , 代入方程得
f (u ) du =2g (x ) dx
6.线性方程
一阶线性微分方程的标准形式:
dy
+P (x ) y =Q (x ) dx
当Q (x ) ≡0, 上方程称为齐次的. 当Q (x ) 0, 上方程称为非齐次的. 例如
dy dx =y +x 2, =x sin t +t 2, 线性的; dx dt
yy '-2xy =3, y '-cos y =1, 非线性的。
一阶线性微分方程的解法
dy
(1)线性齐次方程+P (x ) y =0. (使用分离变量法)
dx
dy dy =-P (x ) dx , ⎰=-⎰P (x ) dx , ln y =-⎰P (x ) dx +ln C , y y
-⎰P (x ) dx
y =Ce . 齐次方程的通解为
(2)线性非齐次方程讨论
dy
+P (x ) y =Q (x ). dx
⎤dy ⎡Q (x )
=⎢-P (x ) ⎥dx , y ⎣y ⎦
Q (x ) Q (x )
-⎰P (x ) dx , 设⎰dx 为v (x ), ∴ln y =v (x ) -⎰P (x ) dx , y y
两边积分ln y =⎰
即y =e v (x ) e -⎰P (x ) dx . 非齐次方程通解形式与齐次方程通解相比:C ⇒u (x ) ⎰作变换y =u (x ) e
-P (x ) dx
-⎰P (x ) dx -⎰P (x ) dx
''y =u (x ) e +u (x )[-P (x )]e ,
-P (x ) dx
=Q (x ), 将y 和y '代入原方程得u '(x ) e ⎰
积分得u (x ) =Q (x ) e ⎰
⎰
P (x ) dx
dx +C ,
一阶线性非齐次微分方程的通解为:
-⎰P (x ) dx -⎰P (x ) dx P (x ) dx P (x ) dx -⎰P (x ) dx ⎰⎰=Ce +e ⋅⎰Q (x ) e dx y =[⎰Q (x ) e dx +C ]e
对应齐次方程通解非齐次方程特解
7. 伯努利方程
伯努利(Bernoulli ) 方程的标准形式
dy
+P (x ) y =Q (x ) y n (n ≠0,1) dx
方程为线性微分方程. 当n ≠0,1时,方程为非线性微分方程. 当n =0,1时,
解法:需经过变量代换化为线性微分方程.
dy dz dy +P (x ) y 1-n =Q (x ), 令z =y 1-n , 则=(1-n ) y -n , 两端除以y n ,得y -n
dx dx dx
dz
代入上式 +(1-n ) P (x ) z =(1-n ) Q (x ),
dx 求出通解后,将z =y
1-n
代入即得y
1-n
=z =e ⎰
-(1-n ) P (x ) dx
(1-n ) P (x ) dx ⎰(⎰Q (x )(1-n ) e dx +C )
8. 可降阶高阶微分方程 (1)y ''=f (x ) 型的微分方程 方程特点:方程右端仅含有自变量x
只需积分两次就可以得到方程的通解.
y '=⎰f (x ) d x +C 1y =⎰[⎰f (x ) d x +C 1]dx +C 2=⎰[⎰f (x ) d x ]d x +C 1x +C 2
这种逐次积分的方法可以用于解更高阶的方程y (n ) =f (x ) 依次通过n 次积分, 可得含n 个任意常数的通解 . (2)y ''=f (x , y ') 型的微分方程 方程特点:方程右端不显含未知函数y
设y '=p (x ) , 则y ''=p ', 原方程化为p '=f (x , p ) (一阶微分方程) 设其通解为p =ϕ(x , C 1) 则得y '=ϕ(x , C 1)
再一次积分, 得原方程的通解y =⎰ϕ(x , C 1) d x +C 2 (3)y ''=f (y , y ') 型的微分方程 方程特点:右端不显含自变量x
则y ''=令y '=p (y ), 即把y 作为自变量,P 为新的未知函数,
d p d p d y d p
=⋅=p
d x d y d x d y
故方程化为p
d p
=f (y , p ) (一阶微分方程)设其通解为p =ϕ(y , C 1), 即得y '=ϕ(y , C 1) d y
d y
=x +C 2
φ(y , C 1)
分离变量后积分, 得原方程的通解⎰9. 线性方程解的结构
(n -1)
'n (x ) y =f (x ). (1)n 阶线性非齐次微分方程y (n ) +P (x ) y + +P 1n -1(x ) y +P
(n -1)
'n (x ) y =0 (2)相应的齐次微分方程为y (n ) +P + +P 1(x ) y n -1(x ) y +P
10. 二阶齐次方程解的结构:
y ''+P (x ) y '+Q (x ) y =0
(1)
定理1. 如果函数y 1(x ) 与y 2(x ) 是方程(1)的两个解, 那末y =C 1y 1+C 2y 2
此性质称为解的线性性质 也是(1)的解.(并不是通解() C 1, C 2是常数)
为解决通解的判别问题。下面引入函数的线性相关与线性无关概念.
定义:设y 1, y 2, , y n 为定义在区间I 内的n 个函数.如果存在n 个不全为零的常数,使得当x 在该区间内有恒等式成立k 1y 1+k 2y 2+ +k n y n =0, 那么称这n 个函数在区间I 内线性相关.否则称线性无关。 特别的:若在I 上有
y 1(x )
≠常数,则函数y 1(x ) 与y 2(x ) 在I 上线性无关 y 2(x )
定理2:如果y 1(x ) 与y 2(x ) 是方程(1)的两个线性无关的特解, 那么y =C 1y 1+C 2y 2就是方程(1)的通解.
例如y ''+y =0, y 1=cos x , y 2=sin x , 且
y 2
=tan x ≠常数, y =C 1cos x +C 2sin x . y 1
则解的实部(1)的解,
定理3. 若y 1(x ) +iy 2(x ) 是齐次方程y ''+P (x ) y '+Q (x ) y =0
y 1(x ) 和虚部y 2(x ) 也是方程的解.
11. 二阶非齐次线性方程的解的结构:
定理4. 设y *是二阶非齐次线性方程y ''+P (x ) y '+Q (x ) y =f (x ) (2)的一个特解, Y 是与(2)对应的齐次方程(1)的通解, 那么y =Y +y *是二阶非齐次线性微分方程(2)的通解.
定理5. 设非齐次方程(2)的右端f (x ) 是几个函数之和,
**
如y ''+P (x ) y '+Q (x ) y =f 1(x ) +f 2(x ) 而y 1与y 2分别是方程,
y ''+P (x ) y '+Q (x ) y =f 1(x ) y ''+P (x ) y '+Q (x ) y =f 2(x )
**的特解, 那么y 1就是原方程的特解. +y 2
定理6. 若y 1(x ) +iy 2(x ) 是线性非齐次方程y ''+p (x ) y '+q (x ) y =f 1(x ) +if 2(x ) 的解,则解的实部y 1和虚部y 2分别是方程y ''+p (x ) y '+q (x ) y =f 1(x ) 和y ''+p (x ) y '+q (x ) y =f 2(x ) 的
解,其中p , q , f 1, f 2, y 1, y 2都是实数。
12. 常数变易法
设对应齐次方程通解为y =C 1y 1+C 2y 2 (3)
'(x ) y 2+c 1(x ) y 1'+c 2(x ) y 2' 设非齐次方程通解为y =c 1(x ) y 1+c 2(x ) y 2y '=c 1'(x ) y 1+c 2''(x ), 设为使y ''中不含C 1''(x ), C 2''(x ), 为使y ''中不含C 1''(x ), C 2
'(x ) y 2'+c 1(x ) y 1''+c 2(x ) y 2'' y ''=c 1'(x ) y 1'+c 2
将y , y ', y ''代入方程(2),得
'(x ) y 2'+c 1(x )(y 1''+P (x ) y 1'+Q (x ) y 1) +c 2(x )(y 2''+P (x ) y 2'+Q (x ) y 2) =f (x ) c 1'(x ) y 1'+c 2
'(x ) y 2'=f (x ) (5) 即c 1'(x ) y 1'+c 2
'(x ) y 2=0y ⎧c 1'(x ) y 1+c 2
(4),(5)联立方程组⎨ 系数行列式w (x ) =1
'(x ) y 2'=f (x ) y 1'⎩c 1'(x ) y 1'+c 2
∴c 1'(x ) =-
y 2f (x ) y f (x ) '(x ) =1, c 2,
w (x ) w (x )
y 2f (x ) y f (x )
dx , c 2(x ) =C 2+⎰1dx , w (x ) w (x )
y 2f (x ) y f (x )
dx +y 2⎰1dx . w (x ) w (x )
y 2
≠0 'y 2
积分可得c 1(x ) =C 1+⎰-
非齐次方程通解为y =C 1y 1+C 2y 2-y 1⎰13. 定义
(n -1)
'n y =f (x ) n 阶常系数线性微分方程的标准形式y (n ) +P + +P 1y n -1y +P
二阶常系数齐次线性方程的标准形式y ''+py '+qy =0 二阶常系数非齐次线性方程的标准形式y ''+py '+qy =f (x )
14. 二阶常系数齐次线性方程解法y ''+py '+qy =0(特征方程法)
设y =e rx , 将其代入上方程, 得(r 2+pr +q ) e rx =0 e rx ≠0, 故有r 2+pr +q =0
特征根r 1,2=
① 有两个不相等的实根(∆>
0) 特征根为r 1=r 2=
r x r x
两个线性无关的特解y 1=e 1, y 2=e 2, 得齐次方程的通解为y =C 1e 1+C 2e 2;
r x r x
② 有两个相等的实根(∆=0) 特征根为r 1=r 2=-
p , 2
',y 2''代入原方程并化简,一特解为y 1=e r 1x , 设另一特解为y 2=u (x ) e r 1x , 将y 2,y 2
u ''+(2r 1+p ) u '+(r 12+pr 1+q ) u =0, 知u ''=0, 取u (x ) =x , 则y 2=xe r 1x ,
得齐次方程的通解为y =(C 1+C 2x ) e r 1x ;
③ 有一对共轭复根(∆
y 1=e (α+j β) x , y 2=e (α-j β) x ,
下面讨论实数形式的解:(欧拉(Euler ) 公式e ix =cos x +i sin x ) 重新组合1=
11
(y 1+y 2) =e αx cos βx , 2=(y 1-y 2) =e αx sin βx , 22j
得齐次方程的通解为y =e αx (C 1cos βx +C 2sin βx ).
15.n 阶常系数齐次线性方程解法
(n -1) n -1
'n y =0特征方程为r n +Pr y (n ) +P + +P + +P 1y n -1y +P 1n -1r +P n =0
(1)首先求相应齐次方程的特征根,确定Q (x ) 的形式。 (2)设非齐方程特解为=Q (x ) e λx 代入原方程
Q ''(x ) +(2λ+p ) Q '(x ) +(λ2+p λ+q ) Q (x ) =P m (x ) 根据多项式恒等,确定待定系数。
特别地y ''+py '+qy =Ae λx
A ⎧λx
e ⎪λ2+p λ+q , λ不是特征方程的根
⎪
A ⎪
=⎨xe λx λ是特征方程的单根
⎪2λ+p ⎪A 2λx
x e λ是特征方程的重根⎪
2⎩
,
17. 线性常系数非齐次方程 (1)f (x ) =e λx P m (x ) 型
y ''+py '+qy =f (x ) 二阶常系数非齐次线性方程
对应齐次方程y ''+py '+qy =0, 通解结构y =Y +,
λx λx λx
常见类型P m (x ), P cos βx , P sin βx , m (x ) e m (x ) e , P m (x ) e
解法:设非齐方程特解为=Q (x ) e λx 代入原方程
Q ''(x ) +(2λ+p ) Q '(x ) +(λ2+p λ+q ) Q (x ) =P m (x )
①若λ不是特征方程的根,λ2+p λ+q ≠0, 可设Q (x ) =Q m (x ), =Q m (x ) e λx ; ②若λ是特征方程的单根,λ2+p λ+q =0, 2λ+p ≠0,
可设Q (x ) =xQ m (x ), =xQ m (x ) e λx ;
③若λ是特征方程的重根,λ2+p λ+q =0, 2λ+p =0,
可设Q (x ) =x 2Q m (x ), =x 2Q m (x ) e λx .
⎧0λ不是根
⎪
综上讨论设=x k e λx Q m (x ) k =⎨1λ是单根,
⎪2λ是重根⎩
上述结论可推广到n 阶常系数非齐次线性微分方程(k 是重根次数). (2)f (x ) =e λx [P l (x )cos ωx +P n (x )sin ωx ]型
e i ωx +e -i ωx e i ωx -e -i ωx
+P n ] 方法一:f (x ) =e [P l cos ωx +P n sin ωx ]利用欧拉公式e [P l
22i
λx
λx
P l P n (λ+i ω) x P l P n (λ-i ω) x
=(+) e +(-) e =P (x ) e (λ+i ω) x +(x ) e (λ-i ω) x ,
22i 22i 设y ''+py '+qy =P (x ) e (λ+i ω) x , 1=x k Q m e (λ+i ω) x ,
k (λ-i ω) x
, y ''+py '+qy =(x ) e (λ-i ω) x , 2=x m e
(1)(2)
∴=x k e λx [Q m e i ωx +m e -i ωx ]=x k e λx [R m (x )cos ωx +R m (x )sin ωx ],
⎧0λ±i ω不是根其中R (x ), R (x ) 是m 次多项式,m =max {l , n }k =⎨,
1λ±i ω是单根⎩
(1)
m
(2)m
注意:上述结论可推广到n 阶常系数非齐次线性微分方程.
方法二:f (x ) =e [P l (x )cos ωx +P n (x )sin ωx ]=e P l (x )cos ωx +e P n (x )sin ωx
λx
λx
λx
=Re[P l (x ) e λ+i ωx ]+Im[P n (x ) e λ+i ωx ] f 1(x ) f 2(x )
**
y ''+py '+qy =f 1(x ), y ''+py '+qy =f 2(x ) 分别求下面两个方程的特解y 1, y 2:**
则y *=Re(y 1) +Im(y 2) 为原方程的特解。
(3)n 阶线性常系数非齐次方程
二阶方程的待定系数法可以推广用来求解n 阶方程,只需根据λ与特征根的重数关系把特解形式作相应的修改。 (4)欧拉方程
形如x n y (n ) +p 1x n -1y (n -1) + +p n -1xy '+p n y =f (x ) 的方程(其中p 1, p 2 p n 为常数) 叫欧拉方程
特点:各项未知函数导数的阶数与乘积因子自变量的方次数相同. 解法:欧拉方程是特殊的变系数方程,通过变量代换可化为常系数微分方程 作变量变换x =e t 或t =ln x ,
1⎛d 2y dy ⎫dy dy dt 1dy d 2y
=2 2-⎪, ==, 将自变量换为t
dt ⎭x ⎝dt dx dt dx x dt dx 2
d 3y 1⎛d 3y d 2y dy ⎫
=-3+2 ⎪,
dt ⎭dx 3x 3⎝dt 3dt 2用D 表示对自变量t 求导的运算
d
, dt
d 2y dy =(D 2-D ) y =D (D -1) y 上述结果可以写为xy '=Dy , x y ''=2-dt dt 2
d 3y d 2y dy x y '''=3-32+2=(D 3-3D 2+2D ) y =D (D -1)(D -2) y dt dt dt 3
一般地,x k y (k ) =D (D -1) (D -k +1) y . 将上式代入欧拉方程,则化为以t 为自变量的常系数线性微分方程. 求出这个方程的解后,把t 换为ln x ,即得到原方程的解.
18. 常系数线性微分方程组
(1)微分方程组:由几个微分方程联立而成的方程组称为微分方程组.
注意:这几个微分方程联立起来共同确定了几个具有同一自变量的函数.
(2)常系数线性微分方程组 微分方程组中的每一个微分方程都是常系数线性微分方程叫做常系数线性微分方程组.
(3)常系数线性微分方程组的解法
步骤:
①从方程组中消去一些未知函数及其各阶导数,得到只含有一个未知函数的高阶常系数线性微分方程.
②解此高阶微分方程, 求出满足该方程的未知函数.
③把已求得的函数带入原方程组, 一般说来,不必经过积分就可求出其余的未知函数.
10