专题一(二阶常微分方程解法)

专题一(二阶常微分方程解法)
专题一(二阶常微分方程解法)

二阶微分方程:

时为非齐次

时为齐次,0)(0)()()()(22≠≡=++x f x f x f y x Q dx dy x P dx y d 二阶常系数齐次线性微分方程及其解法:

2

122,)(2,,(*)0)(1,0(*)r r y y y r r q pr r q p qy y p y 式的两个根、求出的系数;

式中的系数及常数项恰好是,,其中、写出特征方程:求解步骤:

为常数;

,其中?'''=++?=+'+''式的通解:出的不同情况,按下表写、根据(*),321r r

为常数;

型,为常数

,]sin )(cos )([)()()(,)(x x P x x P e x f x P e x f q p x f qy y p y n l x m x ωωλλλ+===+'+''

二阶常系数非齐次线性微分方程的一般形式是 ''+'+=y py qy f x () (1)

其中p q ,是常数。

方程(1)的通解为对应的齐次方程

0=+'+''qy y p y (2)

[

的通解Y 和方程(1)的一个特解*y 之和。即 *y Y y +=.我们已解决了求二阶常系数齐

次线性方程通解的问题,所以,我们只需讨论求二阶常系数非齐次线性微分方程的特解*

y 的方法。

下面我们只介绍当方程(1)中的)(x f 为如下两种常见形式时求其特解*y 的方法。 一、

f x e P x x m ()()=?λ型

由于方程(1)右端函数f x ()是指数函数e x λ?与m 次多项式P

x m ()的乘积,而指数函数与多项式的乘积的导数仍是这类函数,因此,我们推测:

方程(1)的特解应为

y e Q x x *?=λ()( Q x ()是某个次数待定的多项式 ) y

e Q x e Q x x x *??'=+'λλλ()()

y e Q x Q x Q x x *?"=?+'+''λλλ[()()()]22

代入方程(1),得

e Q x p Q x p q Q x e P x x x m λλλλλ???''++'+++≡?[()()()()()]()22 消去e

x λ?,得

【 ''++'+++≡Q x p Q x p q Q x P x m ()()()()()()22λλλ (3)

讨论 01、如果λ不是特征方程

r pr q 20++=的根。 即

02≠++q p λλ 由于P x m

()是一个m 次的多项式,欲使(3)的两端恒等,那未Q x ()必为一个m 次多项式,设为

Q x b x b x b x b m m m m m ()=++++--0111

将之代入(3),比较恒等式两端x 的同次幂的系数,就得到以b b b b

m m 01

1,,,, -为未知数的m +1个线性方程的联立方程组,解此方程组可得到这m +1个待定的系数,并得到特解 y e Q x x m *?=λ()

02、如果λ是特征方程

r pr q 20++=的单根。 即

λλ20++=p q ,但 20λ+≠p ^

欲使(3)式的两端恒等,那么'Q x ()必是一个m 次多项式。

因此,可令 Q x x Q x m ()()=?

并且用同样的方法来确定)(x Q 的系数b b b b m m 0

11,,,, -。 03、如果λ是特征方程r pr q 20++=的二重根。

λλ20++=p q ,且 20λ+=p 。 欲使(3)式的两端恒等,那么''Q x ()必是一个m 次多项式

因此, 可令

Q x x Q x m ()()=?2

并且用同样的方法来确定)(x Q 的系数b b b b m m 0

11,,,, -。

综上所述,我们有结论 如果

f x e P x x m ()()=?λ,则方程(1)的特解形式为 》

y x Q x e k m x *?=()λ

其中Q x m ()是与P x m ()同次的多项式,k 的取值应满足条件 k =?????012λλλ不是特征方程的根是特征方程的单根

是特征方程的二重根

例1求 ''-'+=y y y xe x

562 的通解。

解 特征方程为 0652=+-r r

特征根为 3,221==r r

齐次方程的通解为 x x e C e C Y 3221+=

因为2=λ是特征单根,所以,设非齐次方程的特解为

y x b x b e x *=+()012

~

*'y =+++[()]222020112b x b b x b e x

*''y =++++[()]484240201012b x b b x b b e x

将上述三式代入原方程,得

()-+-≡2200122b x b b e xe x x ,

比较恒等式两端的系数,得

-=-=???2120001b b b

解得

21

0-

=b , 11-=b 因此 x

e x x y 2)121(*+-= 所以方程的通解为

y c e

c e x x e x x x =+-+12232121() 】

二、f x e P x x P x x x l n ()[()cos ()sin ]=+λωω型

由于方程(1)右端函数为

[]x x p x x p e n l x ωωλsin )(cos )(+,这种形式得到非齐次方程的特解*y 的过程稍微复杂些,所以我们这里就只给出结论

y x e R x x R x x k x m m *=+λωω[()cos ()sin ]()()12

其中,R x m ()()1、R x m ()()2是两个m 次多项式,m l n =max{,},

???++=是特征方程的根若不是特征方程的根若ωλωλi i k 10 例2求方程

''+=y y x x cos2 的通解。 解 特征方程 r 210+=

特征根

i r ±=2,1 齐次方程的通解为 x C x C Y sin cos 21+=

这里1,2,0===m ωλ,由于i i 2=+ωλ不是特征方程的根,所以设方程的特解为 &

y ax b x cx d x

*=+++()cos ()sin 22

代入原方程,得 x x x a d Cx x C b ax 2cos 2sin )433(2cos )433(=++-+--

比较两端同类项的系数,得

???????=--=-=+-=-0430

304313a d C C b a

解得

94,0,0,31===-=d C b a 于是

y x x x *=-+132492cos sin

所以非齐次方程的通解为 y c x c x x x x =+-+12132492cos sin cos sin

相关主题
相关文档
最新文档