常微分方程考研讲义第四章 高阶微分方程 下载本文

为了获得更明显的物理意义,令

M?H*sin?*,N?H*cos?*

即令

H*?M2?N2??2np

?2?p2H(??p)?4np22222 (4.54)

及 tg?*?这时(4.53)可以写成

??H*sin?*cospt?H*cos?*sinpt?H*sin(pt??*) ?因此,(4.52)的通解为

??Ae?ntsin(?1t??)?(4.55)

H(?2?p2)2?4n2p2sin(pt??*)

从(4.55)可以看出,摆的运动由两部分叠加而成,第一部分就有阻尼的自由振动,它是系统本身的固有振动,它随时间的增长而衰减,第二部分是由外力而引起的强迫振动项,它的振幅不随时间的增长而衰减。因此,考虑强迫振动时主要就考虑后一项

H(?2?p2)2?4n2p2sin(pt??*),它与外力的频率一样,但相位和振幅都不同了。

我们现在来研究外力的圆频率p取什么值时所引起的强迫振动项的振幅H*达到最大值。

从(4.54)看出,只需讨论当p取何值时(?2?p2)2?4n2p2达到最小值即可。为此,记?(p)?(?2?p2)2?4n2p2,将它对p求导数,并令导数等于零,得到

??(p)??4p(?2?p2)?8n2p

因此,只要2n2??2,即只要阻尼很小时,就解得

p??2?2n2 (4.56)

而当p取此值时,我们有???(p)?8p2?0,因而?(p)在

p??2?2n2 时达到最小值。

把(4.56)代入(4.54),得到相应的最大振幅值为

*Hmax?H4n?4n(??2n)4222?H2n??n22 就是说,当外力的圆频率p??2?2n2时,强迫振动项的振幅达到最大值,这时的圆频率称为共振频率,所产生的现象也叫共振现象。

§4.3 高阶方程的降阶和幂级数解法

4.3.1 可降阶的一些方程类型

n阶微分方程一般地可写为

F(t,x,x?,?,x(n))?0

下面讨论三类特殊方程的降阶问题:

1)方程不显含未知函数x,或更一般地,设方程不含x,x?,?,x(k?1),即方程呈形状:

F(t,x(k),x(k?1),?,x(n))?0 (1?k?n) (4.57)

若令x(k)?y,则方程即降为关于y的n?k阶方程

F(t,y,y?,?,y(n?k))?0 (4.58)

如果能够求得方程(4.58)的通解

y??(t,c1,c2,?,cn?k)

即 x(k)??(t,c1,c2,?,cn?k) 再经过k次积分得到

x??(t,c1,c2,?,cn)

其中c1,c2,?,cn为任意常数。可以验证,这就是方程(4.57)的通解。

特别地,若二阶方程不显含x(相当于n?2,k?1的情形)。则用变换x??y便把方程化为一阶方程。

d5x1d4x?0的解。 例1 求方程5?4dttdtd4x解 令4?y,则方程化为

dtdy1?y?0 dxtd4x这是一阶方程,积分后得y?ct,即4?ct,于是

dtx?c1t5?c2t3?c3t2?c4t?c5

其中c1,c2,?,c5为任意常数,这就是原方程的通解。 2)不显含自变量t的方程

F(x,x?,?,x(n))?0 (4.59)

我们指出,若令x??y,并以它为新未知函数,而视x为新自变量,则方程就可降低一阶。

事实上,在所作假定下,x??y,x???22dydydy?x??y,dtdxdxdy?dy?x????y???y22,…,采用数学归纳法不难证明,可用表出(k?n)。将这些

dx?dx?表达式代入(4.59)就得到

?dydn?1y?G?x,y,,?,n?1??0

dxdx??这是关于x,y的n?1阶方程,比原方程(4.59)低一阶。 例2 求解方程xx???(x?)2?0。 解 令x??y,直接计算可得x???ydydy?y2?0得到y?0或,于是原方程化为xydxdxdyccx?y?0,积分后得y?,即x??,所以x2?c1t?c2(c1?2c),这就是原dxxx方程的通解。

例3 求数学摆的运动方程

d2?g??sin? dt2l满足初始条件:t?0时,???0?0,

d??0的解。 dtd2?dpd??p,则2?p解 令,这时,方程变为 dtdtd?p积分之,得到

dpg??sin? d?l12gp?(cos??c1) 2l1?d??g或者 ??(cos??c1) (4.60) ?2?dt?l这里c1是任意常数。用初始条件代入(4.60),得到c1??cos?0。于 是(4.60)变为

2?d??2g(cos??cos?0) ???dtl??将上式开方得到

2d?2g??cos??cos?0 (4.61) dtl我们先讨论摆从最大的正偏离角???0到最大的负偏离角????0之间的第一次摆动的情况,这时

d??0,(4.61)的右端取负号,得到 dtd?2g??cos??cos?0 (4.62) dtl将方程(4.62)分离变量,然后积分,并计及初始条件即得

??令

?0td?2g2g???dt??t (4.63)

0llcos??cos?0t0?l?0d?

2g?0cos??cos?0则(4.63)可些写为

l?d?t0?t? (4.64)

2g?0cos??cos?0