微分方程(1-3)
第9章 微分方程与差分方程
第1节 微分方程的基本概念
我们已经知道,利用函数关系可以对客观事物的规律性进行研究. 而在许多几何,物理,经济和其他领域所提供的实际问题,即使经过分析、处理和适当的简化后,我们也只是能列出含有未知函数及其导数的关系式. 这种含有未知函数的导数的关系式就是所谓的微分方程. 求出微分方程中的未知函数的过程就叫解微分方程.
本章主要介绍微分方程的一些基本概念和几种常用的微分方程的解法.
实际问题中的数据大多数是按等时间间隔周期统计的. 因此,有关变量的取值是离散变化的,处理他们之间的关系和变化规律就是本章最后的内容——差分方程.
含有未知函数的导数或微分的方程称为微分方程. 微分方程中出现的未知函数的导数的最高阶数称为微分方程的阶.
现实世界中的许多实际问题,例如,物体的冷却,人口的增长,琴弦的振动,电磁波的传播等,都可以归结为微分方程问题. 这时微分方程也称为所研究问题的数学模型.
例9.1 质量为m 的物体只受重力作用由静止开始自由垂直降落. 根据牛顿第二定律:物体所受的力F 等于物体的质量m 与物体运动的加速度的乘积,即F 的铅垂线为x 轴,其正向向下. 下落的起点为原点. 记开始下落的时间t 距离x 与时间t 的函数关系x =
=ma . 取物体降落=0,则物体下落的
x (t ) 满足
d 2x
=g , (9.1) 2dt
其中g 为重力加速度常数. 这就是一个2阶微分方程。
例9.2 产品的月产量为x 时的边际成本
1
(9.2) c '(x ) =x +8,2
就是一个1阶微分方程.
在微分方程中,若未知函数是一元函数就称为常微分方程;若未知函数是多元函数,就称为偏微分方程. 本章只讨论常微分方程。
n 阶微分方程的一般形式是
(9.3) F (x , y , y ', y '', , y (n ) ) =0,
其中x 为自变量,y =
y (x ) 是未知函数,上式(9.3)中,y (n ) 必须出现,而其余变量(包
括低阶导数)可以不出现.
如果能从式(9.3)中解出最高阶导数得到微分方程的如下形式
y (n ) =f (x , y , y ', y '', , y (n -1) ) (9.4)
以后我们只讨论姓如式(9.4)的微分方程,并假设式(9.4)右端的函数连续.
特别地,式(9.4)中的
f
在所讨论的范围内
f
如果能写成如下形式
y (n ) +a 1(x ) y (n -1) + +a n -1(x ) y '+a n (x ) y =g (x ) (9.5)
则称式(9.5)为n 阶线性微分方程. 其中a 1(x ), , a n (x ) 和g (x ) 均为自变量x 的已知函数. 把不能表示成形如式(9.5)的微分方程称为非线性微分方程.
例9.3 试指出下列方程是什么方程,并指出微分方程的阶数. (1)
dy dy
=x 3+y (2)sin x +(cosx ) y +tan x =0 dx dx
2
3
d y d 3y dy ⎛dy ⎫
-3=5y (3)x (4)+ln x +xy =x ⎪23dx dx dx ⎝dx ⎭
解 方程(1)是一阶线性微分方程. 因为程. 因为两边除以sin x 就可看出.
dy
和y 都是一次. 方程(2)也是一阶线性微分方dx
3
⎛dy ⎫
方程(3)是2阶非线性微分方程,因为其中含有 ⎪
⎝dx ⎭
.
d 3y dy
方程(4)是3阶线性微分方程. 因为, , y 都是一次式. 3
dx dx
如果一个函数代入微分方程能使方程式为恒等式,则称这个函数为该微分方程的解. 例如,(a)x =
121
gt ,(b)x =gt 2+c 1t +c 2都是例9.1中的微分方程9.1的解,22
其中c 1, c 2为任意常数.
通常,称不含任意常数的解为微分方程的特解. 而含有相互独立的任意常数,且任意常
数的个数与微分方程的阶数相等的解为微分方程的通解(一般解).
这里所说的相互独立的任意常数,是指它们取不同的值时就得到不同的解. 从而不能通过合并而使得通解中的任意常数的个数减少.
上面的解中,(a)和(c)分别是方程(9.1)和(9.2)的特解,(b)和(d)分别是方程(9.1)和(9.2)的通解.
在实际问题中通常都要求寻找满足某些附加条件的解. 此时,这类附加条件就可以用来确定通解中的任意常数. 这类附加条件称为初始条件,也称为定解条件.
一般地,一阶微分方程y '=
f (x , y ) 的初始条件为
y x =x =y 0 (9.6)
其中x 0, y 0都是已知常数.
二阶微分方程y ''=
f (x , y , y ') 的初始条件为
y x =x =y 0, y 'x =x =y 0' (9.7)
带有初始条件的微分方程称为微分方程的初值问题. 微分方程的解的图形是一条曲线,称为微分方程的积分曲线. 例9.4 验证函数y =(x
3
+c )cos x (c 为任意常数)是方程 dy
+y tan x -3x 2cos x =0 dx
的通解,并求出满足初始条件y
x =0
=0的特解.
解 要验证一个函数是否是微分方程的通解,只要将函数代入方程,验证是否恒等,再看函数式中所含的独立的任意常数的个数是否与方程的阶数相同.
对y =(x
3
+c )cos x ,求一阶导数
dy
=3x 2cos x -(x 3+c )sin x dx
把
y 和
dy
代入方程左端,得 dx
dy
+y tan x -3x 2cos x =3x 2cos x -(x 3+c )sin x +(x 3+c )cos x tan x -3x 2cos x =0 dx
因为方程两边恒等,且
y
中含有一个任意常数,方程又是一阶的,故
. y =(x 3+c ) c o s 是题设方程的通解x 把初始条件
y x =0=0代入通解y =(x 3+c )cos x 中,得c =0. 从而所求特解为
y =x 3cos x .
习题9-1
1、 指出下列微分方程的阶数
(1)xy ''-2yy '+2x =0 (2)(y '') (3)(x
22
-y 3sin x +x 5=0
+3) dx +(4x +5y 2) dy =0
2、指出下列各题中的函数是否为所给微分方程的解. (1)xy '=2y , y =5x (2)y ''-
2
22y
'y +2=0, y =c 1x +c 2x 2 x x
λ1x
(3)y ''-(λ1+λ2) y '+λ1λ2y =0, y =c 1e 3、验证
+c 2e λ2x
y =cx +
x =0
12
(c 为任意常数)是方程x (y ') -yy '+1=0的通解,并求满足初c
始条件y =2的特解.
4、设曲线在点(x , y ) 处的切线的斜率等于该点横坐标的平方,试建立曲线所满足的微分方程,并求出通解.
习题9-1答案
1、(1)2阶 (2)2阶 (3)1阶 2、(1)是 (2)是 (3)是 3、特解为y
=
1
x +2 2
dy 143
=x ,通解为y =x +c 4、微分方程为
4dx
第2节 一阶微分方程
微分方程没有统一的解法,必须根据微分方程的不同类型,研究相应的解法. 本节我们将介绍可分离变量的微分方程以及一些可以化为这类方程的微分方程,如齐次方程等.
一、可分离变量的微分方程. 在一阶微分方程
dy
=F (x , y ) 中,如果右端函数能分解成F (x , y ) =f (x ) g (y ) ,
dx
x 与y 分离,x 的一个函数f (x ) 与y 的一个函数g (y ) 相乘的形式,即
dy
=f (x ) g (y ) (9.8) dx
其中
f (x ) ,g (y ) 都是连续函数. 根据这种方程的特点,我们可以通过积分的方法来求解.
设g (y ) ≠0. 用g (y ) 除方程(9.8)的两端,用dx 乘以方程的两端,使得未知函数y 的某已知函数及其微分与自变量x 的某已知函数及其微分置于等号的两边(又一次分离了x 与
y )得
1
dy =f (x ) dx g (y )
再对上述等式两边积分,即得
1
⎰g (y ) dy =⎰f (x ) dx (9.9)
积分出来以后就表明y 是x 的一个(隐)函数(关系),就是方程(9.8)的解. 如果g (y 0) =0,则易验证y
=y 0也是方程(9.8)的解.
上述求解可分离变量的微分方程的方法,称为分离变量法. 例9.5 求微分方程
xydx +dy =x 2dy +xdx
的通解.
解 先合并dx , dy 的各项得
x (y -1) dx =(x 2-1) dy
设y -1≠0, x
2
-1≠0, 分离变量得
dy x
=2dx y -1x -1
dy x
两端积分
⎰y -1=⎰x 2-1dx
得 ln |
11
y -1|=ln |x 2-1|+ln |c 1|
22
2
于是 (y -1)
=±c 1(x 2-1)
2
记c =±c 1,则得到题设方程的通解为(y -1)
=c (x 2-1)
dy
例9.6 求微分方程=e x y 的通解.
dx
解 分离变量后两边积分
dy x
=e ⎰y ⎰dx
得 ln |从而 y
y |=e x +ln |c 1|
=±c 1e
e x
记c =±c 1,则得到题设方程的通解为 y
=ce e
x
例9.7 一曲线通过点(3,2),它在两坐标轴间的任一切线线段均被切点所平分,求曲线的方程.
解 设曲线的方程为y =
y (x ) . 曲线上任一点(x , y ) 的切线方程为
Y -y
=y '
X -x
由假设,切点(x , y ) 的切线位于两坐标轴间的线段的两个端点分别是X
=0时,Y =2y 和
Y =0时,X =2x . 将这两个端点代入切线方程都得到曲线所满足的微分方程
dy y
=-
x dx
分离变量后积分,得到通解为xy
=c
将初始条件y |x =3=2代入通解得c =6. 从而所求的曲线方程为xy 二、齐次方程 如果一阶微分方程
=6.
dy =f (x , y )
dx
中的函数
f (x , y ) 可以写成
y ⎛y ⎫
的函数,即f (x , y ) =ϕ ⎪,于是 x ⎝x ⎭
dy ⎛y ⎫
=ϕ ⎪ (9.10) dx ⎝x ⎭
这称为齐次方程.
齐次方程可以通过引进新的未知函数的方法化成为可分离变量的微分方程. 令u
=
y dy du ,u 是x 的一个新的未知函数. 则 y =ux , =x +u , x dx dx
du
原齐次方程变成 x +u =ϕ(u )
dx
du dx
分离变量后积分得 =⎰=ln |x |+c ⎰x ϕ(u ) -u
记Φ(u ) 为
1
的一个原函数,则得通解为 Φ(u ) =ln |x |+c
ϕ(u ) -u
再以
y ⎛y ⎫
代替u , 就得所给齐次方程的通解 Φ ⎪=ln |x |+c
x x ⎝⎭
2
例9.8 求微分方程 (xy -x 解 原方程变形为
) dx -(y 2-xy ) dy =0 的通解.
y -1
dy xy -x =2=2dx y -xy ⎛y ⎫y
⎪-⎝x ⎭x
2
就是一个齐次方程
令u =
y dy du ,则 y =ux , =x +u x dx dx
du u -1
+u =2
dx u -u
u 1
du =dx 2
1-u x
代入齐次方程得 x
分离变量,u ≠0, x ≠0时,得
两边积分
u 1du =⎰1-u 2⎰x dx
1
ln |1-u 2|=ln |x |+ln |c 1| 2
得 -
y 1y
以代替u 就得到原方程的通解 -ln |1-|=ln |x |+ln |c 1| x 2x
1记c =±2
c 1
y c
得 1-=2
x x
2
从而 x -xy =c .
注. 本题也可以直接分离变量法求解.
x (x -y ) dx =y (y -x ) dy
y -x ≠0时, ydy =-xdx 12121
积分得 y =-x +c
222
即y
2
+x 2=c 为原方程的通解.
2
这样本题得到两个通解形式x 定要包含所有解!
三、一阶线性微分方程 方程
-xy =c 和y 2+x 2=c . 说明微分方程的通解并不一
dy
+p (x ) y =Q (x ) (9.11) dx
叫做一阶线性微分方程,它对于未知函数y 及其导数y '都是一次的. 如果Q (x ) ≡0,则方
程(9.11)称为齐次的,否则就称为非齐次的.
对于齐次一阶线性微分方程
dy
+p (x ) y =0 (9.12) dx
-p (x ) dx
(9.13) =Ce ⎰
通过分离变量积分,可得它的通解
y
而对于非齐次一阶线性微分方程(9.11),我们可以利用它相应的齐次一阶线性微分方程(9.12)的通解(9.13),并使用所谓常数变易法来求非齐次方程(9.11)的通解,这种方法是把齐次方程(9.12)的通解(9.13)中的任意常数C 变易换成x 的未知函数u (x ) ,即作变换
-p (x ) dx
(9.14) y =ue ⎰
假设(9.14)是非齐次方程(9.11)的解,代入(9.11)中进而求出u (x ) ,再代入(9.14)就得到非齐次方程(9.11)的解. 为此,将(9.14)对x 求导,注意u 是x 的函数,得
-⎰p (x ) dx dy du -⎰p (x ) dx
(9.15) =e -up (x ) e
dx dx
将(9.15)和(9.14)代入(9.11),得
du -⎰p (x ) dx
e =Q (x ) dx
分离变量后积分得
u
=⎰Q (x ) e ⎰
-p (x ) dx
p (x ) dx
dx +C (9.16)
p (x ) dx ⎰Q (x ) e dx (9.17)
将(9.16)代入(9.14)就得到(9.11)的通解
y =Ce ⎰+e ⎰
-p (x ) dx
⎰
易见,一阶非齐次线性方程的通解(9.17)是对应的一阶齐次线性方程的通解(9.13)与其本身的一个特解((9.17)中取C 方程也成立.
例9.9 求方程y '+
=0的解) 之和. 此后还可看到,这个结论对高阶非齐次线性
1cos x
y =的通解. x x
1cos x
解 题设方程是一阶非齐次线性方程,这时p (x ) =, Q (x ) =.
x x
于是,按公式(9.17),所求通解为
11⎰x dx cos x ⎰x dx
y =Ce +e ⎰x e dx
-ln x -ln x cos x ln x
=Ce +e ⎰e dx
x
C 1
=+⎰cos xdx x x C 1=+sin x x x
-
⎰x dx
1
-
例9.10 求方程
dy
+3y =8的通解. dx
解 这是一个非齐次线性一阶方程. 下面不利用公式(9.17),而采用常数变易法来求解. 先求解相应的齐次方程的通解. 由
dy
+3y =0 dx
-3x
分离变量后积分得相应齐次方程的通解 y =c 1e , 其中c 1为任意常数.
-3x
利用常数变易法,将c 1变易为u (x ) ,即设原非齐次方程的通解为 y =ue 求导得
dy du -3x
=e -3ue -3x dx dx
du -3x
代入原非齐次方程得 e =8
dx
分离变量后积分得 u (x ) =
83x 3x
8e dx =e +C ⎰3
8-3x
从而得到原非齐次方程的通解为 y =+Ce
3
习题9-2
1、求下列微分方程的通解 (1)x (y
2
-1) dx +y (x 2-1) dy =0
dy
=3x +y (2)dx
2、求下列微分方程的通解
(1
)xy '-
(2)y (x 2y -=0 -xy +y 2) dx +x (x 2+xy +y 2) dy =0
3、求下列微分方程的通解
(1)y '+y =e -x
y =sin x (2)xy '+
4、求下列微分方程的初值问题:
(1)cos ydx +(1+e
(2)(x +1) y '--x )sin ydy =0, y |x =0=π4 y =(x +1) 2e x , y |x =0=1
5、已知某产品生产的总成本C 由可变成本与固定成本两部分组成. 可变成本y 是产量x 的函数,且y 关于x 的变化率等于
成本函数c =c (x ) .
习题9-2答案
1、(1)(x 22xy ,当x =10时,y =1;固定成本为22x +y 100. 求总-1)(y 2-1) =C ; (2)3x +3-y =C
=Cx 2; (2)xy =Ce
-x ⎛y ⎫-arctan ⎪⎝x ⎭2、(1
)y + 3、(1)y =(x +C ) e
4、(1
)(1+e x ; (2)y 1=(C -cos x ) x )sec y = (2)y =(x +1) e x
99-1) 5
、C (x ) =100+2第3节 可降阶的二阶微分方程
本节讨论三种特殊形式的二阶微分方程的求解.
一、y ''=f (x ) 型
这种简形的方程,其解法就是多次积分.
在y ''=f (x ) 两端积分,得 y '=⎰f (x ) dx +C 1
=⎰[⎰f (x ) dx +C 1]dx +C 2=⎰⎰f (x ) dxdx +C 1x +C 2
(n ) 再次积分,得 y 注:对于n 阶微分方程y
意常数的通解.
例9.11 求方程y ''=e
解 连续积分两次,得 2x =f (x ) ,显然也可以连续积分n 次,就得到含有n 个任+sin x 的通解.
1y '=e 2x -cos x +C 12 1y =e 2x +sin x +C 1x +C 24
这就是所求通解.
二、y ''=f (x , y ') 型
这种类型的特征是不显含y ,求解方法是:
令y '=
p (x ) ,则y ''=p '(x ) ,则原二阶方程化成了一阶方程 p '=f (x , p )
dy 利用上一节的方法求出它的通解p =ϕ(x , C 1) ,再根据=y '=p =ϕ(x , C 1) 也是一dx
阶方程.
直接积分得y =⎰ϕ(x , C 1) dx +C 2,就是原二阶微分方程的通解.
注:由于一阶微分方程p '=f (x , p ) ,我们并不都会求解. 因此本类型y ''=f (x , y ') 方程的求解还不能说都可求出.
例9.12 求方程y ''=1y '+xe x 的通解. x
解 令p =y ',原方程化成
p '-1p =xe x x
的一阶线性微分方程. 从而 p =c 1e ⎰
=c 1x +1-(-) dx x +e ⎰1-(-) dx x ⎰(-) dx x ⎰x xe e dx 1x ⎰e x dx
=c 1x +xe x
x 'p =y =c x +xe 1即
因此,原方程的通解为
y =⎰(c 1x +xe x ) dx +c 2
12=c 1x +(x -1) e x +c 2 2
三、y ''=f (y , y ') 型
y 的函数,令p =y ',从这种类型的特征是不明显地含x . 这时我们把x 看成自变量
而p 也是y 的函数. 再利用复合函数的求导法则,把对x 的导数y ''化为对y 的导数,即 y ''=dp dp dy dp =⋅=⋅p dx dy dx dy
于是,y ''=f (y , y ') 就变成了 p dp =f (y , p ) dy
这样就得到一个关于y , p 的一阶微分方程.
设y '=p =ϕ(y , c 1) 是它的通解,那么分离变量再积分就得到原方程的通解为 dy ⎰ϕ(y , c 1) =x +c 2.
注. 一阶微分方程
求出解来. p dp =ϕ(y , c 1) 不一定会求解,因此本类型y ''=f (y , y ') 也不一定能dy
例9.13 求方程y ''=yy '的通解.
解 令
这时y ''=p =y ',将x 看作是y 的函数. dp dp dy dp =⋅=p dx dy dx dy 代入原方程就得到一个一阶方程 dp p =y dy
12分离变量再积分得p =y +c 1 2
再解一阶微分方程y '=
分离变量再积分得
p =12y +c 1 2
⎛⎫dy x +c 2=⎰2⎛⎫c 1+y 1+2就是原方程的通解.
习题9-3
1、 求下列方程的通解
(1)y ''=x -cos x
(2)y ''=x +y '
(3)y ''=(y +1) y '
2、求下列微分方程初始问题的特解.
(1)y ''=e
(2)y ''=3x , y |x =0=0, y '|x =0=0 1y ', y |x =1=0, y '|x =1=2 x
2(3)yy ''-(y ') -y '=0, y |x =0=2, y '|x =0=1
习题9-3答案
1、(1)y =1
6x 3+cos x +c 1x +c 2
(2)y =c x
1e -xe x +c 2
(3
)x +c =22、(1)y =1
9e 3x -1
3x -1
9
(2)y =x 2-1
(3)y =e x +1