高阶微分方程的解法及应用
哈尔滨学院本科毕业论文(设计)
题目:院(系)专 业 年 级 姓 名 指导教师
高阶微分方程的解法及应用
理学院 数学与应用数学 2009级 刘晓辉 学 号 09031212 徐亚兰
职 称 副教授
2013年6月1日
哈尔滨学院本科毕业论文(设计)
目 录
摘 要 ............................................................................................................................................. 1 Abstract ......................................................................................................................................... 2 前 言 ............................................................................................................................................. 3 第一章 高阶微分方程的理论与结构 ........................................................................................... 4 第二章 高阶常系数线性微分方程 ............................................................................................. 6 2.1 高阶常系数线性齐次微分方程 ........................................................................................ 6 2.1.1 特征根是单根的情况 ................................................................................................. 6 2.1.2 特征根是重根的情况 ................................................................................................. 7 2.2 高阶常系数线性非齐次方程 ............................................................................................ 8 2.2.1 常数变易法 ................................................................................................................. 8 2.2.2 比较系数法 ............................................................................................................... 10 2.2.3 拉普拉斯变换法 ....................................................................................................... 11 2.3 Euler方程 ........................................................................................................................ 13 第三章 可降阶的高阶微分方程的解法 ................................................................................... 15
dny
3.1 形如nf(x)的高阶方程 .......................................................................................... 15
dx
3.2 形如F(x,y(k),y(k1),,y(n))0的高阶方程 ................................................................. 16 3.3 形如F(y,y,,y(n))0的高阶方程 ............................................................................. 17 3.4 恰当导数方程 .................................................................................................................. 19 第四章 高阶微分方程的应用 ................................................................................................... 21 参考文献 ....................................................................................................................................... 25 致 谢 ........................................................................................................................................... 26
摘 要
本文首先介绍了高阶微分方程的一些理论与结构。进而介绍了高阶齐次线性微分方程的求解方法和高阶非齐次线性微分方程的求解方法,在求解齐次线 性微分方程里主要采用了特征根法;在求解非齐次线性微分方 程里主要采用了比较系数法、拉普拉斯变换法和常
dny
数变易法。其次又介绍了几类可降阶的微分方程的解法,主要有形如nf(x),
dx
F(x,y(k),y(k1),,y(n))0,F(y,y,,y(n))0,恰当导数方程和Euler方程的降阶方法,并且研究了几类较为复杂的高阶微分方程的降阶问题。最后通过一些在现实生活中例子对这些方法的具体应用做了介绍。
关键词:高阶常微分方程;常数变易法;特征根法;降阶法
Abstract
This paper introduces some of the theories and higher order differential structure. Then introduce higher-order homogeneous linear differential equation methods and high-order non-homogeneous linear differential equation method for solving homogeneous linear differential equation where the main use of the eigenvalue method; in solving inhomogeneous linear differential equations in mainly uses the comparison coefficient method, Laplace transform method and the constant variation. And secondly describes several types of differential equations can be reduced for the solution, the main tangible eg, appropriate derivative equations and Euler equations reduction method, and studied several types of more complex higher order differential equations reduction problem. Finally some real life examples of specific applications of these methods have been described.
Key words: Higher Order Ordinary Differential Equations; constant variation; eigenvalue method; reduction method
前 言
常微分方程作为数学系重要专业的一门基础课程,对学习好其他的科目起到了至关重要的作用。它的形成与发展是和力学、天文学、物理学,以及其他科学技术的发展密切相关的。数学的其他分支的新发展,如复变函数、拓扑学等,都对常微分方程的发展产生了深刻的影响,当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具。而高阶微分方程是常微分方程中的一个重要的组成部分,在现实的生活中也有着广泛的应用,比如工程问题。常系数线性微分方程的解法,高阶微分方程的降阶问题又是高阶微分方程的重中之重。
常微分方程是在生产实践和科学技术中产生的。目前,常微分方程在很多学科领域内有着重要的应用,自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等。这些问题都可以化为求常微分方程的解,或者化为研究解的性质的问题。人们对于二阶以及简单的高阶微分方程求解的方法有了很多理论成果,而高阶常微分方程并没有固定的解法,例如,高阶常系数线性齐次微分方程,我们可以运用特征根的方法进行求解,高阶常系数线性非齐次微分方程,我们可以运用常数变易法,比较系数法,拉普拉斯变换法进行求解。而对于可以降阶的高阶微分方程,我们通常采用降阶法,也就是通过一定的变换把高阶微分方程求解的问题转化成低阶微分方程的求解问题。
dny
本篇论文我总结了形如nf(x),F(x,y(k),y(k1),,y(n))0,F(y,y,,y(n))0,
dx
恰当导数方程和Euler方程的降阶方法,并且研究了几类较为复杂的高阶微分方程的降阶问题,进而介绍此类问题在科学技术中的应用。
第一章 高阶微分方程的理论与结构
定义1(方程的阶) 在一个常微分方程里,未知函数的最高阶导数的阶数叫做方程的阶。
n阶隐式方程的一般形式为
F(x,y,y',,y(n))0
n阶显式方程的一般形式为
y(n)f(x,y,y',,y(n1))
定义2(解) 设函数y(x)在区间I上有直到n阶的导数。如果把y(x)代入到
方程F(x,y,y',,y(n))0
得到在区间I上关于x的恒等式是
F(x,(x),'(x),,(n)(x))0
则称y(x)是方程F(x,y,y',,y(n))0在区间I上的一个解。
微分方程的解可以包括任意的常数,其中任意常数的个数可以多到和方程的阶数相等,当然也可以不包括任意常数。我们把方程
F(x,y,y',,y(n))0
的含有n个独立的任意常数C1,C2,,Cn的解
y(x,C1,C2,,Cn)
称做该方程的通解。如果方程的解y(x)不包含任意常数,则把它叫做特解。
方程
y(n)an1(x)y(n1)a1(x)y'a0(x)yf(x) (1-1)
称做n阶线性微分方程,它关于未知函数y以及各阶导数y',y'',,y(n)都是线性的。在这 里,我们通常假设ai(x)(i0,1,,n1)和f(x)是区间(a,b)上的连续函数。如果
ai(x)(i0,1,,n1)都是常数,则把方程(1-1)叫做n阶常系数线性方程。如果方程的右端项f(x)0,即
y(n)an1(x)y(n1)a1(x)y'a0(x)y0
则称方程(1-1)是齐次的,否则为非齐次的。所以对于方程(1-1)的齐次方程是
y(n)an1(x)y(n1)a1(x)y'a0(x)y0 (1-2)
定理1(叠加原理) 设y1(x)和y2(x)是齐次方程(1-2)的解,则对于任意常数c1和
c2,c1y1(x)c2y2(x)也是方程(1-2)的解。
定理2 设yi(x)(i1,2)是方程
y(n)an1(x)y(n1)a1(x)y'a0(x)yfi(x)
的解,则y1(x)y2(x)也是方程
y(n)an1(x)y(n1)a1(x)y'a0(x)yf1(x)f2(x)
的解。
定理3 设y1(x),y2(x),,yn(x)是齐次方程(1-2)的n个线性无关的特解,则
yciyi(x)
i1n
是方程(1-2)的通解,其中c1,c2,,cn是任意常数。
定理4 设y*(x)是非齐次线性方程(1-1)的任意一个确定的解,Y(x)是(1-1)对应的齐次线性方程(1-2)的通解。则
y(x)Y(x)y*(x)
是(1-1)的通解。
第二章 高阶常系数线性微分方程
2.1 高阶常系数线性齐次微分方程
对于n阶常系数线性齐次方程
ynan1yn1a1y'a0y0 (2-1)
其中y是关于x的未知函数,系数an1,,a1,a0是实常数。如果yex是方程的根,把他代入到方程中,得
(nan1n1a1a0)ex0
因为ex0,因此
nan1n1a1a00 (2-2) 反之,如果满足等式(2-2),则yex是方程(2-1)的解。式子(2-2)是关于的n次代数方程,则把他叫做微分方程(2-1)的特征方程,它的根就称做特征根。下面根据特征根的不同情形分别进行讨论方程解的情况。
2.1.1 特征根是单根的情况
定义 我们把P()na1n1an1an0称为方程
y(n)a1y(n1)an1y'any0
的特征方程,它的根叫做特征根。在这里把叫做待定系数。
定理 如果特征方程P()na1n1an1an0有n个互异的根
1,2,n,则
y1e1x,y2e2x,,ynenx
是方程
y(n)a1y(n1)an1y'any0
的一个基本解组。
特征方程P()na1n1an1an0可能有复根,由于他的系数是实的,他的复根一定是共轭成对的出现。即此时在相异特征根1,2,,n中有复数。例如
k1aib(a,b为实数),则k1aib也是P()na1n1an1an0的根。这
两个特征根所对应的解是实变量复值函数
yke(aib)xeaxcosbxieaxsinbxyk1e
(aib)x
ecosbxiesinbx
axax
d4x
例1 求方程4x0的通解。
dt
解 特征方程410的根是11,21,3i,4i,其中有两个实根和两个复根,但他们都是单根,所以所求方程的通解是
xc1etc2etc3costc4sint
在这里c1,c2,c3,c4是任意的常数。
2.1.2 特征根是重根的情况
定理 假设方程y(n)a1y(n1)an1y'any0有互异的特征根1,2,p,他们的重数分别是m1,m2,,mp,mi1,并且m1m2mpn,则与他们相对应的
y(n)a1y(n1)an1y'any0的特解是
e1x,xe1x,,xm11e1xe2x,xe2x,,xm21e2x
,
xxm1xep,xep,,xpep
并且该特解构成y(n)a1y(n1)an1y'any0在区间(,)上的基本解组。
例2 解初值问题
y(4)y0
''''''
y(0)y(0)y(0)0,y(0)1
解 特征方程是41(21)(21)0,特征根是1,i.所以方程的通解是
yc1exc2exc3cosxc4sinx.
又因为
y'c1exc2exc3sinxc4cosx,
y''c1exc2exc3cosxc4sinx, y'''c1exc2exc3sinxc4cosx,
根据初始条件,得
c1c2c30ccc0124
ccc0231c1c2c41
再解方程组,得
c1
111,c2,c30,c4 4421x1
(eex)sinx 42
于是初值问题的解是
y
2.2 高阶常系数线性非齐次方程
对于n阶常系数线性非齐次方程
ynan1yn1a1y'a0yf(x) (2-3)
他的通解等于齐次方程的通解再加上加其对应的非齐次方程的一个特解。在上一节中我们知道了怎样求解齐次方程的通解,下面我们主要来研究求解非齐次方程的特解的方法。
2.2.1 常数变易法
常数变易法 实际上是一种变量变换的方法,在这里我们 简单的介绍一下 在n阶方程中的应用。
可以设方程(2-3)的特 解是:
x(t)c1(t)x1(t)c2(t)x2(t)cn(t)xn(t) (2-4) ~
其中ci,i1,2,,n是待定的 常函数。并且把它代入到方程(2-3)中,再附加上n-1个条件,就可以得到 方程组
''
x1c1'(t)x2c2(t)xncn(t)0''''''
xc(t)xc(t)xcn(t)01122n
(2-5)
x(n2)c'(t)x(n2)c'(t)x(n2)c'(t)0
122nn1
(n1)'(n1)'(n1)'x1c1(t)x2c2(t)xncn(t)f(t)
(t),c2(t),,cn(t)的表 达式,把它们分别进行 积分进解方程 组(2-5)就会得到关于c1
x(t)。而得 到ci,i1,2,,n,再把它 们代入到(2-4)式中,继而求得方程(2-3)的一个特解~
由于这种方法 对于自由项f(t)的形式没有任何的限制,因此使用的范围会比较广,但是求解的工作量相对来说会大一些。
例3 求解 方程
x''x
1
cost
的通解,已知 它所对应的齐次线性微 分方程的基本解组是cost,sint。
解 运用常数变易法,设
xc1(t)costc2(t)sint
并且把它代入到方程里,就可以得到关于c'1(t)和c'2(t)的两个方程
c'1(t)costc'2(t)sint0
和
'
sintc1'(t)costc2(t)
1
cost
解得
c1'(t)
据此得到
sint
,c'2(t)1 cost
c1(t)lncost1,c2(t)t2
所以原方程的通解是
x1cost2sintcostlncosttsint
其中1,2是 任意的常数。
2.2.2 比较系数法
对于常系数非线性方程(2-3),我们通常采用的方法是比较系数法,它是把所要求解的微分方程的问题转化成代数问题,在自由项是
f(t)pm(t)et或f(t)pn(t)costps(t)sintet
(其中pm(t),pn(t),ps(t)分别是m次,n次,s次的多项式。,,都是实常数)时,就可以确定
x的形式,即分别令 特解~
~xtkQm(t)et
(Qm(t)是一个待定的m次的多项式,k是方程x(n)p1x(n1)p2x(n2)pnxf(t)的特征方程有根时的次数)或者
(1)(2)~xtkQm(t)costQm(t)sintet
(1)(2)
(t),Qm(t)是两个待定的m次多项式,k是方程含有根t的次(其中mmaxn,s.Qm
数)然后把它代入到方程(2-3)中,再进行比较等式的左右两边同次幂的系数来确定待定系数多项式。再根据线性微分方程解的结构便可以求解出方程的通解。
例4 求方程x'''3x''xet(t5)的通解。 解 特征方程
33231(1)30
有三重根1,2,31,所对应的齐次方程的通解是
x(c1c2tc3t2)et
并且方程有~xt3(ABt)et的特解,将它代入到方程中得
(6A24Bt)etet(t5)
再比较两边的系数求得
51
A,B
624
进而
13~xt(t20)et 24
所以所求的方程的通解是
x(c1c2tc3t2)et
13
t(t20)et 24
其中c1,c2,c3是任意的常数。
2.2.3 拉普拉斯变换法
根据积分
F(s)estf(t)dt
所定义的确定在复平面(Res)上的复变数s的函数F(s),叫做函数f(t)的拉普拉斯变换,其中f(t)在t0上有定义,并且满足不等式
f(t)Met
在这里M,是某两个正常数。我们把f(t)称为原函数,而把F(s)称为象函数。
设所给定的微分方程
dnxdn1x
a1n1anxf(t) (2-6) ndtdt
和初始条件
'(n1)
x(0)x0,x'(0)x0,,x(n1)(0)x0
其中a1,a2,,an是常数,而f(t)是连续的并且满足原函数的条件。
可以证明,假如x(t)是方程(2-1)的任意解,则x(t)以及它的各阶的导数
x(k)(t)(k1,2,,n)都是原函数。记
F(s)f(t)estf(t)dt
0
X(s)x(t)estx(t)dt
那么,根据原函数的微分性质就有
x'(t)sX(s)x0
'(n1) x(n)(t)snX(s)sn1x0sn2x0x0
于是,再对方程(2-6)的两边进行拉普拉斯变换,并且运用线性性质就可以得到
'(n2)(n1)
snX(s)sn1x0sn2x0sx0x0'(n2)a1sn1X(s)sn2x0sn3x0x0
an1sX(s)x0anX(s)F(s)
即
(sna1sn1an1san)X(s)F(s)(sn1a1sn2an1)x0
'(n1)(sn2a1sn3an2)x0x0
或者
A(s)X(s)F(s)B(s)
其中A(s),B(s)和F(s)都是已知的多项式,由此得到
X(s)
F(s)B(s)
A(s)
这就是方程(2-6)所满足所给初始条件的解x(t)的象函数。而x(t)可以直接查表或者根据反变换公式计算求解出来。
例5 求方程x'''3x''3x'x1的满足初始条件x(0)x'(0)x''(0)0的解。 解 对方程的左右两边进行拉普拉斯变换得到
(s33s23s1)X(s)
1 s
由此得到
X(s)
1
s(s1)3
再把上面式子的右面分解成为部分分式
11111
323
ss1(s1)s(s1)(s1)
对上面式子的右端的各项分别求出或者查表得出他们的原函数,则他们的和就是X(s)的原函数
11
x(t)1ettett2et1(t22t2)et
22
这就是所要求的解。
2.3 Euler方程
定义:形如
anxny(n)an1xn1y(n1)a1xy'a0yf(x)
的方程叫做Euler方程,其中an,an1,,a1,a0是实常数,并且an0。它的特点就是包含y的k阶导数项的系数是akxk。当x0时,各阶导数项的系数是0,所以我们令x0。在现实里,我们仅需要考虑x0的这种情况,因为在x0的时候,在上述的方程里做自变量变换tx,则方程就化成
n1
dnyydyn1dantatata0yf(t) n11nn1
dtdtdtn
求出他的解,再用x替换t就可以得出方程关于x0的解。 再做自变量变换
lnx
则xe,
dydyd1dy
dxdtdxxd
d2y1d2ydy
2(2) 2
ddxxdd3y1d3yd2ydy
(32) 3332
ddxxdd
一般的,假如
dky1dkydk1ydy
(bb),其中bk1,,b1是常数,则 k11kkkk1
ddxxdddk1y1
dxk1xk1
dk1ydkyd2ydy
(bk)(bkb)kbk1121k1 k2
dddd
k
dkydydky
,,k的常系数的线性组合,进而把Euler方程 所以,对于每一个正整数k,x是k
ddxd
anxny(n)an1xn1y(n1)a1xy'a0yf(x)
化成了常系数线性方程。
例5 求解方程
x3y'''5x2y''7xy'8y0,x0
解 设tlnx,方程化做
d3yd2ydy
248y0 32
dtdtdt
他的特征方程是
32248(2)(24)0
他的特征值是-2和2i,方程的通解是 1
yc1e2tc2cos2tc3sin2tc12c2cos(2lnx)c3sin(2lnx)
x
第三章 可降阶的高阶微分方程的解法
本部分我将介绍4类比较常见的高阶微分方程的解法,在这些解法里有一个比较类似的思路,就是把这些的高阶微分方程通过某些变换降成比较低阶的微分方程再进行求解。所以,我们把这种方法称做“降阶法”。
3.1 形如dny
dx
nf(x)的高阶方程
方程
dny
dx
nf(x) 这种类型的方程比较简单,通常令y(n1)z,则zf(x).
积分得zf(x)dxC(n1)1,也就得到yf(x)dxC(n2)1.同理可以令yz,得到y(n2)f(x)dxC1dxC2 (f(x)dx)dxC1xC2
.
如此继续下去,再通过n次积分就可以求出(3-1)的通解是
ydxdxf(x)dx
C1(n1)!xn1C2
(n2)!
xn2Cn1xCn. 例1 求解微分方程ymsinxcosx的通解。 解 对原方程的左右两边依次进行积分,得
ycosxsinxC1,
ysinxcosxC1xC2.
再次进行积分,求解出原方程的通解是
ycosxsinx
C12
x2
C2xC3. 例2 求方程y'''e2xcosx的通解
3-1)(
解 y''
12x1
esinxC1,y'e2xcosxC1xC2 24
所以所求原方程的通解是
11
ye2xsinxC1x2C2xC3
82
3.2 形如F(x,y(k),y(k1),,y(n))0的高阶方程
方程
(k)k(1)n()
F(x,y (3-2a) ,y,,y) 0
这种方程的特点极其容易看出,方程不显含y或者y,y,y,,y(k1).在这时我们只要把
y(k)z,代入到上述的方程中,原方程就可以化作
F(x,z,z,,z(nk))0. (3-2b)
如果方程(3-2b)可以求解出通解
zz(x,C1,,Cnk).
则再对方程
y(k)z(x,C1,,Cnk)
积分k次,便可以求出y了。在这里需要注意的是每积分一次,就需要增加一个独立的任意常数.
例1 求解方程
(1x2)y2xy.
解 设yp(x),则yp,代入到上面的方程中得
(1x2)p2xp.
再积分得
lnpln(1x2)lnC1,
即
pC1(1x2).
积分四次就可以求解出原方程的通解
yC2x3C3xC4.
例2 求微分方程xy''y'0的通解 解 设y'p则y''
dp
。所以原方程就可以写成 dx
x
dpdpdxp0, dxpx
左右两边进行积分得
lnplnxlnC1
所以p
C1C
,也就是y'1,两边再次进行积分得出 xx
yC1lnxC2
3.3 形如F(y,y,,y(n))0的高阶方程
方程
F(y,y,,y(n))0 (3-3)
这类方程也有一定的特点,就是不显含自变量x,这时,我们总可以利用代换yp,使这类方程降低一阶。以二阶方程
F(y,y,y)0
为例,设yp,于是便有
y
dpdpdydp
p. dxdydxdy
代入到原方程中,就有
F(y,p,p
dp
)0. dy
这是一个关于未知函数p的一个一阶方程。如果用它可求出pp(y,C),就有
ypp(y,C).
这是一个关于x,y的变量可分离方程,进而可以求解出通积分。
1y2
例1 求解方程y.
2y
解 根据上面的分析,我们可以令yp,则yp
dp
,代到原方程里得 dy
dp1p2p. dy2y
即
2pdpdy
.
1p2y
左右两边进行积分得
ln(1p2)lnylnC1.
求解出p得
p
积分后得
xC2. 于是便有
y
C12C3
xxC4. 48
例2 求解方程a2y2(1y2)3.
解 首先令yp,则yp,于是原方程就化成了
3dp122
(1p). dxa
再令ptant,则
dp
dx
px
sec2tdt1
(1tant2t)a
32
acostdt.
即
xasintC1.
进而
p即
p进而可以积分求解出通积分
(yC2)2(xC1)2a2.
3.4 恰当导数方程
假如方程
F(x,y,y,,y(n))0 (3-4)
的左面正好是某一个函数(x,y,y,,y(n1))对x的导数,则(3-4)就可化为
d
(x,y,y,,y(n1))0. dx
于是我们就把(3-4)称作恰当导数方程。
其实这类方程的解法和全微分方程的解法很类似,可以降低一阶,化成
(x,y,y,,y(n1))C
之后,再想办法求解这个方程。
例1 求解方程 yyy20. 解 我们可以把方程写成
d
(yy)0,所以就有yyC1.即 dx
ydyC1dx.
积分后就可以得出通积分
y2C1xC2.
这样的问题虽然简单,但是需要具有很强的观察能力和比较牢固的基础才可以观察出来。下面有一个关于这方面的例子,解法技巧很高明,关键还是配导数的方法。
例2 求解方程yyy20.
解 经过观察我们可以先把等式的两边同时乘以一个不是0的因子
1
,便有 所以
从而通解是
y
2yyy2y2
ddx(y
y
)0. yC1y.
yC2eC1x.
第四章 高阶微分方程的应用
要利用微分方程解决实际问题,首先必须要根据物理和几何关系规律来建立微分方程,然后再对进一步的问题进行分析与微分方程的建设,并且考虑初始条件,边界条件,收敛条件来确定定解的条件,这是数学建模过程。
模型建立好了就有了微分方程,我们就可以根据前面的内容来解除方程,因为解决的是实际问题,我们还要用解出来的结果来分析问题。这部分内容因为实际应用相对比较强,所以我用三个简单的例子来简单的介绍一下。
例1 设质量是m的物体自由悬挂在一个一端固定的弹簧上,当重力跟弹簧力相互抵消的时候,物体就会处在一个平衡的状态,若用手向下拉物体使它离开平衡的位置以后放开,物体在弹性力与阻力的作用下将做往复运动,阻力的大小与运动的速度成正比,方向相反。
(1) 建立位移所满足的微分方程
解 设时刻t物体的位移是x(t)。
1. 自由振动的情况,物体所受到的力有弹力恢复力fcx, 阻力R根据牛顿第二定律得到
d2xdx
m2cx. dtdt
dx
. dt
令2n
m
,k2
c
,就可以得到阻尼自由振动方程 m
d2xdx2
2nkx0. 2dtdt
2. 强迫振动情况,如果物体在运动的过程里还受到铅直外力FHsinpt的作用,设
h
得到强迫振动的方程
H
, m
d2xdx2
2nkxhsinpt. 2dtdt
(2)在没有外力的作用下做自由运动,设t0时物体的位置是xx0,初始的速度为v0,
求物体的运动规律xx(t).
解 由(1)知,位移满足的定解问题是
d2xdx2
2nkx0,dt2dt
x|x,dx|v.t00t00dt
1. 无阻尼自由振动的情况(n0).
d2x
方程 2k2x0.
dt
解得方程的通解是
xC1cosktC2sinkt.
再利用初始条件可以得到
C1x0,C2
所以所求的特解是
v0
. k
xx0coskt
v0
sinktk
kxtan0).v0
Asin(kt)(Ad2xdx
1有阻尼自由振动情况,方程22nk2x0,特征方程是 r22nrk2
0,特征值
dtdtr1,2n
在这个时候我们需要三种情况来进行讨论 小阻尼nk,则
xent(C1costC2sint).(
大阻尼nk,则
xC1er1tC2er2t.
临界阻尼nk,则
x(C1C2t)ent.
例2 人类将要向宇宙发射一颗人造地球卫星,为了让她摆脱地球引力,初始速度应该不少于第二宇宙速度,试求该速度。
解 在物理问题中,关键是要通过建立模型,把物理问题转化成数学问题,在这个题目里设人造地球卫星的质量是m,地球的质量是M,卫星的质心到地心的距离是h,根据牛顿第二定律得
d2hGMm
m22. (G是引力系数) dth
再设卫星的初速度是v0,已知地球的半径R63105,于是就有初值问题
d2hGM
,dt22
h.
hR,dhv.
0t0
dtto
于是以上的物理问题就转化成了求二阶常微分方程的特解的问题,设方程组的第二个式子就可以得到
dh
v(h),代入到上述dt
v
从而就有
dvGM2, dhh
GM
dh. 2h
vdv
两边进行积分得到
12GMvC. 2h
再利用初始条件得
12GMCv0.
2R
所以
12121
vv0GM(R). 22h
注意到
12121
vv0GM. h22Rlim
为了让v0,v0应该满足
v0
(4-1)
因为在h=R(在地面上)时,引力跟重力是相等的,即
GMm2
1s/ )mg. (g9.8m2
R
所以
GMR2g.
代入到方程(4-1)中得
v011.210(m/s).
3
这就说明第二宇宙速度是11.2km/s.
例3 在船上向海里沉放某一种探测器,按照探测的要求,需要确定仪器的下沉深度y和下沉的速度v之间的函数关系。假设仪器在重力的作用下在海平面由静止开始往下沉,在下沉的过程中还受到了阻力和浮力的作用,我们设仪器的质量是m,体积是B,海水的比重是,仪器所受到的阻力跟下沉的速度成正比,比例系数是k(k0),试着建立y与v所满足的微分方程,并且求出函数关系式yy(v).
解 同样也是首先把实际的问题转化成常微分方程,根据题目中的条件和牛顿第二定律可以将问题转化成求解初始条件是v|y00的二阶微分方程
d2y
m2mgBkv. dt
的特解的问题。我们有
d2ydvdvdydv
v. 2dtdtdydtdy
得
mv
dv
mgBkv. dy
初始条件是
v|y00,
再利用分离变量法解上述初值问题得
y
mm(mgB)mgBkv
vln. 2kkmgB
参考文献
[1]金福临,李训经.常微分方程.上海科技出版社, 1997.(书籍)
[2]同济大学应用数学系.高等数学(下册, 第五版) .高等教育出版社,2002.(书籍) [3]华北师范大学常微分方程教研室.常微分方程(第二版).高等教育出版社,2005.(书籍) [4]史金麟,张剑峰.常微分方程分类原理.科学出版社,2003.(书籍) [5]张伟年,杜正东,徐冰.常微分方程.高等教育出版社,2006.(书籍)
[6]时宝,张德存,盖明久.常微分方程理论及应用.国防工业出版社,2010.(书籍) [7]叶彦谦.常微分方程讲义(第二版).人民出版社,1982.(书籍)
[8]化存才,赵奎奇等.常微分方程解法与建模应用选讲.科学出版社,2010.(书籍) [9]丁同仁,李承治.常微分方程.高等教育出版社,1985.(书籍)
[10]斯米尔诺夫著,孙念增译.高等数学教程.人民教育出版社,1955.(书籍) [11]周义仓等.常微分方程及其应用.科学出版社,2010.(书籍) [12]武卓群,李勇.常微分方程.高等教育出版社, 2004.(书籍)
[13]王高雄,周之铭,朱思铭,王寿松.常微分方程.高等教育出版社,1983.(书籍)
[14] Hughes,A. Testing for Language Teachers. Cambridge:Cambridge University Press,1989(期刊) [15] Nunan,D. The Learner-Centred Curriculum. Cambridge:Cambridge University Press,1988(期刊)
致 谢
作为即将从哈尔滨学院理学院毕业的我,在四年的大学生活中,认真学习了各科专业知识,积极参加了社会实践活动。在大四的师范实习期间,我的教学方面有了显著的提高,特别是在教态、教学方法、教学过程中与学生的沟通技能方面有了明显的改进。回首大学四年时光,匆匆而过,我要诚挚的感谢教育和培养我的老师们,感谢徐亚兰老师对我完成论文的选题、撰写方面给予的指导和帮助.在撰写论文期间,徐老师对我的谆谆教导、和她的严谨的治学态度使我终身受益,对我未来参加工作必将产生深远的影响。我真诚感谢哈尔滨学院理学院的各位领导和老师等给我长期的传道授业解惑,对本论文在撰写过程中给我的知识指导、帮助和启迪。借此机会向在大学四年里给予我的帮助和关怀的同学们表示感谢。论文中还有诸多的问题和不足之处,敬请大家给予批评指证。