二阶变系数线性微分方程的一些解法
第九节 二阶变系数线性微分方程
的一些解法
常系数线性齐次方程和某些特殊自由项的常系数线性非齐次方程的解法已在第七节中介绍,而对于变系数线性方程,要求其解一般是很困难的。本节介绍处理这类方程的二种方法
§9.1 降阶法
在第五节中我们利用变量替换法使方程降阶,从而求得方程的解,这种方法也可用于二阶变系数线性方程的求解。
考虑二阶线性齐次方程
dyd2y
2+p(x) +q(x)y=0 (9.1)
dxdx
设已知其一个非零特解y1,作变量替换,令 y=uy1 (9.2)
其中u=u(x)为未知函数,求导数有
dydudy1
=y1+u
dxdxdx
dudy1d2yd2ud2y1
求二阶导数有2=y12+2+u2
dxdxdxdxdx
代入(9.1)式得
dy1dudy1d2ud2y1
y12+(2+p(x)y1)+(2+p(x)
dxdxdxdxdx
+q(x)y1)u=0 (9.3)
这是一个关于u的二阶线性齐次方程,各项系数是x的已知函数,因为y1是(9.1)的解,所以其中
dy1d2y1
2+p(x) +q(x)y1≡0
dxdx
故(9.3)式化为 dy1dud2u
y12+(2+p(x)y1) =0
dxdxdx
dy
再作变量替换,令=z得
dx
dzdy1
y1+(2+p(x)y1)z=0
dxdx
21
分离变量 dz=-[+p(x)]dx
y1z
两边积分,得其通解
C2-∫p(x)dx
z=2e 其中C2为任意常数
y1
1-∫p(x)dx
积分得u=C2∫2edx+C1代回原变量得(9.1)
y1
的通解
1-∫p(x)dx
y=y1[C1+C2∫2edx]
y1
此式称为二阶线性方程的刘维尔(Liouville)公式。
综上所述,对于二阶线性齐次方程,若已知其一个非零特解,作二次变换,即作变换y=y1∫zdx可将其降为一阶线性齐次方程,从而求得通解。 对于二阶线性非齐次方程,若已知其对应的齐次方程的一个特解,用同样的变换,因为这种变换并不影响方程的右端,所以也能使非齐次方程降低一阶。
sinxd2y2dy
例1. 已知y1=是方程2++y=0的
xxdxdx
一个解,试求方程的通解
解 作变换 y=y1∫zdx
dydy1
则有 =y1z+∫zdx
dxdx
22
dzdydydy11
=y1+2z+2∫zdx 2
dxdxdxdx
代入原方程,并注意到y1是原方程的解,有
dzdy1dy1
y1+(2+)z=0
dxdxdx
dz
即 =-2ctanx·z
dx
C1
积分得 z=2
sinx
sinxC1
于是 y =y1∫zdx=[∫2dx+C2]
xsinx
sinx= (-C1ctanx+C2)
x1
= (C2sinx-C1cosx) x
这就是原方程的通解。
§9.2 常数变易法
在第三节求一阶线性非齐方程通解时,我们曾对其对应的齐次方程的通解,利用常数变易法求得非齐次方程的通解。对于二阶线性非齐次方程
dyd2y
2+p(x) +p(x)y=f(x) (9.4)
dxdx
其中p(x),q(x),f(x)在某区间上连续,如果其对应的齐次方程
dyd2y
2+p(x) +q(x)y=0
dxdx
的通解 y=C1y1+C2y2已经求得。 那么也可通过如下的常数变易法求得非齐次方程的通解。
设非齐次方程(9.4)具有形式 y=u1y1+u2y2 (9.5)
的特解,其中u1=u1(x),u2=u(x)是两个待定函数,对y求导数得
~
~
y'=u1y′1+u2y′2+y1u′1+y2u′2 由于用(9.5)代入(9.4),可确定u1,u2的一个方程,为了同时确定这两个函数,还须添加一个条件,为计算方便,我们补充一个条件:y1u′1+y2u′2=0 这样 y'=u1y′1+u2y′2 y
u′1y′1+u′2y′2=f(x)
~
~
~
y1u'1y2u'20与补充条件联列得方程组
y'1u'1y'2y'2u'2f(x)因为y1,y2线性无关,即
y2y2y1y'2y2y'1
≠常数,所以()′=≠0 2
y1y1y1
设w(x)=y1y′2-y2y′1,则有w(x)≠0所以上述方程组有唯一解。 解得
y2f(x)y2f(x)
u'1y1y'2y2y'1w(x)
u'2y1f(x)y1f(x)y1y'2y2y'1w(x)
积分并取其一个原函数得
y2f(x)
u1=-∫dx
w(x)y1f(x)
u2=∫dx
w(x)
y2f(x)
则所求特解为 y=y1∫dx+y2∫
w(x)
~
y1f(x)
dx w(x)
所求方程的通解 y=Y+y=C1y1+C2y2+y1∫
y2f(x)y1f(x)dx+y2∫dx w(x)w(x)
~
上述求特解的方法也适用于常系数非齐次方程情形。
d2y1dy例1. 求方程2-=x的通解
xdxdx
解 先求对应的齐次方程 d2y1dy 2-=0
xdxdx
2
dy1dy
的通解,由 2=
xdxdx
11dy
·d()=dx
dyxdxdx
dy
得 ln||=ln|x|+ln
dx
|C|
dy2
即 =Cx得通解y=C1x+C2
dx
2
所以对应齐次方程的两个线性无关的特解是x和
1。
为求非齐次方程的一个解y将C1,C2换成待定函数u1,u2,且u1,u2满足下列方程
~
x2u'11u'20
2xu'10u'2x
112
解上述方程得 u′1= u′2=-x
22
1x3
积分并取其一原函数得 u1=x,u2=-
26
于是原方程的一个特解为
333xxx2
y=u1·x+u2·1=-=
263
从而原方程的通解为
~
x
y=C1x+C2+
3
2
3
第十节 数学建模(二)——微分方
程在几何、物理中的应用举例
一、镭的衰变
例1. 镭、铀等放射性元素因不断地放出各种射线而逐渐减少其质量,称为放射性物的衰变。由实验得知,衰变速度与现存物质的质量成正比,求放射性元素在时刻t的质量。
解 用x表示该放射性物质在时刻t的现存物质,dx
则表示x在时刻t的衰变速度,于是“衰变速度与dt
现存质量成正比”可表示为
dx
=-kx
dt
这是一个以x为未知函数的一阶方程,它就是放射性元素衰变的数学模型。其中k>0是比例常数,称为衰变常数,因元素的不同而异。方程右端的负号表示
dx
当时间t增加时,质量x减少,即t>0时,<0。
dt
解这个方程得通解 x=Ce-kt
若已知当t=t0时,x=x0,即x|tt0=x0 代入方程可得 C=x0ekt0 得特解 x=x0ek(tt0)
它反映了某种放射性元素衰变的规律。 二、正交轨线
已知曲线族方程F(x,y,C)=0,其中包含了一个参数C,当C固定时就得到一条曲线,当C改变就得整族曲线,称为单参数曲线族。例如y=Cx2为一抛物线族。
图6-3
如果存在另一族曲线G(x,y,C)=0,其每一条曲线都与曲线族F(x,y,C)=0的每条曲线垂直相交,即不同族中的曲线在交点处的切线互相垂直。则称G(x,y,C)=0为F(x,y,C)=0的正交轨线。 将曲线族方程F(x,y,C)=0对x求导与F(x,y,C)=0联列并消去常数C,得曲线族上任一点的坐标(x,y)和曲线在该点的斜率y′所满足的微分方程 f(x,y,y′)=0
这就是曲线族F(x,y,C)=0所满足的微分方程。 因为正交轨线过点(x,y),且在该点与曲线族中过该点的曲线垂直,故正交轨线在点(x,y)处的斜率
1 k=-
y'
于是可知曲线族F(x,y,C)=0的正交轨线满足方程
1
f(x,y,-)=0
y'
这是正交轨线的数学模型,其积分曲线族(通解),就是所要求的正交轨线。
例2 求抛物线族y=Cx2的正交轨线。 解 对y=Cx关于x求导,得y′=2Cx与原方程联列
2
yCx2
消去C
y'2Cx
2y
得微分方程 y′=
x
1
将-代入y′得所求抛物线的正交轨线微分方程
y'
图6-4
12y
-=
y'x
x
即 ydy=-dx
2
x2y22
积分得 +=C
42
即抛物线族 y=Cx的正交轨线是一个椭圆族,如图6-4。 三、追迹问题
例3. 开始时,甲、乙水平距离为1单位,乙从A点沿垂直于OA的直线以等速v0向正比行走;甲从乙的左侧O点出发,始终对准乙以nv0(n>1)的速度追赶,
并问乙行多远时,被甲追到。 图6-5 解 如图6-5建立坐标系,设所求追迹曲线方程为
y=y(x)
经过时刻t,甲在追迹曲线上的点为p(x,y),乙在点B(1,v0t)。于是有
v0ty
tanθ=y′= (10.1)
1x
由题设,曲线的弧长OP为 ∫0y'2dx=nv0t
解出v0t代入(10.1)得
1x
(1-x)y′+y=∫0y'2dx
n
两边对x求导,整理得
x
2
1
(1-x)y″=y'2
n
这就是追迹问题的数学模型。
这是一个不显含y的可降阶的方程,设y′=p,y″=p′代入方程得
1
(1-x)p′=p2
n
dxdp
或 = 2
n(1x)p
1
两边积分得 ln(p+p)=-ln|1-x|
n
+ln|C1|
C12
即 p+p=
将初始条件 y′|x=0=p|x=0=0代入上式,得C1
2
=1,于是
1
y′+y'= (10.2)
x
两边同乘 y′-y'2,并化简得
2
y′-y'2=- (10.3) (10.2)与(10.3)两式相加,得
11
y′= (-x)
2x
n1
1nnn
积分,得 y=[- (1-x)+ (1-
n1n12
x)
n1n
]+C2
n
代入初始条件 y|x=0=0得C2=2,所求追迹
n1
曲线方程为
(1x)n(1x)n
y=[-]+2
n1n12n1
(n>1)
n1n
n1n
甲追到乙时,即曲线上点P的横坐标x=1,此时
n
y=2
n1
n
即乙行走至离A点2个单
n1
位距离时即被甲追到。 四、弹簧振动
下面我们讨论机械振动的简单模型——弹簧振动问题,研究
图6-6
悬挂重物的弹簧的振动,并假定弹簧的质量与重物的质量相比较可以忽略不计。
如图6-6,一弹簧上端固定,下端与一质量为m的物体连接,弹簧对物体的作用力(恢复力)与弹簧的伸长度成正比(比例常数为k);
物体在运过程中所受的
阻力与速度成正比(比例常数为λ)。此外,物体还与一个连杆连接,连杆对物体的作用力(强迫力)为F(t)。下面建立物体运动方程(数学模型)。 如图6-6,物体的平衡位置为原点,向下方向为Ox轴的正向,以x=x(t)表示物体在时刻t的位置,因为物体共受到三个力的作用。
(1)恢复力:一kx (负号表示恢复力与位移x方向相反);
dxdx(2)阻力:-λ (负号表示阻力与速度的方向
dtdt
相反);
(3)强迫力:F(t) 由牛顿第二定律 F=ma
dxd2x
得 m2=F(t)-kx-λ
dtdt
F(t)d2xdxk
或 2++x=
mmdtmdt
这就是物体运动的数学模型——振动方程。
k
为方便起见,记=2β (β>0),=ω2
mmF(t)
(ω>0),=f(t),则上述方程可写成
m
2
dxdx2
2+2β+ωx=f(t) (10.4)
dtdt
1.自由振动,当f(t)≡0时称为自由振动。
分两种情况讨论
(1)当β=0时称为无阻尼自由振动,其运动方程为
dx2
2+ωx=0
dt
其通解 x=C1cosβt+C2sinβt =Asin(ωt+φ)
C1
(其中A=CC,tanφ=)
C2
21
22
2
图6-7
这是简谐振动,如图6-7,这里振幅A及初相角φ,可由物体的初始位置和初始速度决定。
(2)当β≠0时称为有阻尼自由振动,其运动方程为
dxd2x
2+2β+ω2x=0
dtdt
其特征方程为 r2+2βr+ω2=0 下面就其根的三种情形分别讨论: (ⅰ)β>ω(大阻尼情形),其根为 r=-β±
22特征方程有两个不相等的实根,由于它们都是负数,可令r1=-η1,r2=-η2,(η1>0,η2>0)所以方程的通解为
x=C1e1t+C2e2t
图6-8 图6-9
这里的位移x不是周期函数,因而物体不作任何振动,当t→+∞时x→0,即随时间的无限增加而趋于平衡位置,如图6-8(当C1+C2>0,η1C1+η2C2<0的情形)
(ⅱ)β=ω(临界阻尼情形),特征方程有二重根,r1=r2=-β,此时通解为 x=(C1+C2t)e-βt
这是位移x也不是周期函数,物体也不作任何振动,当t→+∞时x→0,即随时间无限增加而趋于平衡位置,如图6-9(当C1>0,C2-βC1<0的情况)。 (ⅲ)0<β<ω(小阻尼情形),特征方程有一对共轭复根 r=-β±i22 此时通解为 x=Ae
这里A,φ都是任意常数,
可由振动的初始条件决定。由
-βt
cos(22t+φ)
上式看到,振幅Ae-βt随时间的增加而减少,
其减少的快慢程度由系数β=决定。当t→+∞
2m
时,振幅Ae-βt→0,于是x→0,即随时间t无限增加而趋于平衡位置。这种情形称为有阻尼的衰减振动,如图6-10所示 2.强迫振动
设外力 f(t)=asinω0t
我们只考虑无阻尼的强迫振动,其振动方程为
图6-10
dx2
2+ωx=asinω0t
dt
a
它的通解为 x=Asin(ωt+φ)+2·sin2
0
2
ω0t,当ω≠ω0时
2t
x=Asin(ωt+φ)-cosωt,当ω=ω0时。
2
由解的形式可以看出,振动由两种运动所合成,一种是自由振动也称固有振动;另一种是由外力所致的振动,称为强迫振动。前一种情况(当ω≠ω0时)强迫
a
振动的振幅为2,当ω与ω0很接近时,振幅就2
0
at
很大;后一种情况(当ω=ω0时)强迫振幅为,当t
2
at
→+∞时,振幅→+∞,这就是共振现象。因此当
2
外力a0sinωt同系统处于共振状态时,将会引起振幅无限增大的振动,这在机械和建筑中一般是必须严格避免的。
五、R、L、C电路中的电振荡
如图6-11所示的简单串联电路,在电路中,电阻为R,电感为L,电容为C及电动势E(t)。由电学知识,在电阻R上的电压降为RI
dI
在电感L上的电压降为L,在电
dt1
容C上的电压降为∫Idt,根据克
C
希霍夫第二定律,得到
图6-11
dI1
L+RI+∫Idt=E(t)
dtC
将方程的两边关于t求导数得。
dI1d2I
L2+R+
dtCdtdE(t)I=
dt
这就是R、L、C电路中的电振荡方程,它与四中所
指述的弹簧振动的运动方程(10.4)形式上完全一样,类似于弹簧振动的情况,在电路中也会发生共振,但是与机械系动不同的是,在电路系统中共振现象大有用处,例如收音机中的调谐电路。