常微分方程95ppt

1、第五节 高阶常系数线性微分方程一、特征方程与特征根二、二阶常系数齐次线性微分方程三、二阶常系数非齐线性微分方程四、 n 阶常系数齐线性微分方程一一. 特征方程与特征根特征方程与特征根( (1 1) ) )()(1) 1(1)(均为常数knnnnpxfypypypy称为n阶常系数非齐线性微分方程。式成为则若 ) 1 ( , 0)( xf方程（2）称为n阶常系数齐线性微分方程。xrer函数为常数时因为,和它的导数只差常数因子,xrey 代入(2)得( r 为待定常数 ),所以令(2)的解为 )( 201) 1(1)( ypypypynnnn称(3)为微分方程(2)的特征方程特征方程, 其根称为特征

2、根特征根.)( 30111 prprprnnnn0)(111xrnnnne prprpr)( 201) 1(1)(ypypypynnnn形如形如) 1 ( 0 yqypy )(常数。常数。实实为为的方程，称为二阶常系数线性齐微分方程，的方程，称为二阶常系数线性齐微分方程， qp、其中其中 得得的解，则代入方程后，的解，则代入方程后，假设方程有形如假设方程有形如xey 02，xxxeqepe即即 02。qp二、二阶常系数线性微分方程二、二阶常系数线性微分方程二阶常系数线性齐微分方程二阶常系数线性齐微分方程) 1 ( 0 yqypy的特征方程为的特征方程为 02。qp ) 121，则，则实根实根特

3、征方程有两个不同的特征方程有两个不同的xxeyey2 1 21 ，是方程是方程 (1) 的两个线性无关的解，故方程的两个线性无关的解，故方程 (1) 的通解为的通解为 2 1 212211。xxeCeCyCyCy二阶常系数线性齐微分方程二阶常系数线性齐微分方程) 1 ( 0 yqypy的特征方程为的特征方程为 02。qp )221，则，则实重根实重根特征方程有特征方程有 ) 1 ( 1 1的一个解。是方程此时，xey 042， qp由求根公式由求根公式 22422, 1 ，pqpp021 p由刘维尔公式求另一个解：由刘维尔公式求另一个解：xeexeeeyxpxxxpxdd)()2(2d21 1

4、 1 1 021p d11 。xxexxe于是，当特征方程有重实根时，方程于是，当特征方程有重实根时，方程 ( 1 ) 的通解为的通解为 )(21211 1 1 。xCCeexCeCyxxx二阶常系数线性齐微分方程二阶常系数线性齐微分方程) 1 ( 0 yqypy的特征方程为的特征方程为 02。qp3) 特征方程有一对共轭复根：特征方程有一对共轭复根： i i21，则，则， )i (2)i (12 1 xxxxeeyeey，是方程是方程 ( 1 ) 的两个线性无关的解，其通解为的两个线性无关的解，其通解为 )i(2)i(12211。xxeCeCyCyCy 1i 。式中的虚数单位利用欧拉公式去掉

5、表达欧拉公式：欧拉公式： sinicosi。e )sini(cosi)i(1，xxeeeeyxxxx )sini(cosi)i(2。xxeeeeyxxxx由线性方程解的性质：由线性方程解的性质： cos)(21211，xeyyyx sin)(i21212xeyyyx均为方程均为方程 ( 1 ) 的解，且它们是线性无关的：的解，且它们是线性无关的：故当特征方程有一对共轭复根故当特征方程有一对共轭复根 i i2 1，时，原方程的通解可表示为时，原方程的通解可表示为 )sincos(21。xCxCeyx ) )Im( , )Re( (1 1 。二阶常系数线性齐微分方程二阶常系数线性齐微分方程 0 y

6、qypy特征方程特征方程 02。qp特特 征征 根根通通 解解 形形 式式)( 2 1 不同的实根xxeCeCy2 1 21)( 2 1 实重根)(211 xCCeyx)( i2, 1 共轭复根)sincos(21xCxCeyx解解 032 的通解。的通解。求方程求方程 yyy 032 2，特征方程特征方程 3 1 2 1 ，特征根 321。所求通解为xxeCeCy例解解 052 的通解。的通解。求方程求方程 yyy 052 2，特征方程特征方程 i21 i21 2 1 ，特征根 )2sin2cos( 21。所求通解为所求通解为xCxCeyx例解解 0 d d2 dd 22满足初始条件的解：满

7、足初始条件的解：求方程求方程ststs 012 2，特征方程特征方程 1 2 1 ，特征根 ) ( 21。所求通解为tCCest 2 d d 4 0 0 。，tttss 2 4 2 d d 4 210 0 ，得得，由初始条件由初始条件CCtsstt故所求特解为故所求特解为 ) 24(。test例解解 的弹簧从静止状态的弹簧从静止状态用手将悬挂着的质量为用手将悬挂着的质量为 m此时弹簧仅受到弹性恢复力此时弹簧仅受到弹性恢复力 f 的作用。求反映此弹的作用。求反映此弹 O 0时，时，的位移为的位移为当点当点xx 突然放手，突然放手，开始拉长，开始拉长，簧运动的规律（设其弹性系数为簧运动的规律（设其

8、弹性系数为 k ）。）。O例解解 的弹簧从静止状态的弹簧从静止状态用手将悬挂着的质量为用手将悬挂着的质量为 m此时弹簧仅受到弹性恢复力此时弹簧仅受到弹性恢复力 f 的作用。求反映此弹的作用。求反映此弹 O 0时，时，的位移为的位移为当点当点xx 突然放手，突然放手，开始拉长，开始拉长，簧运动的规律（设其弹性系数为簧运动的规律（设其弹性系数为 k ）。）。O0 xx取取 x 轴如如图所示。轴如如图所示。由力学的虎克定理，有由力学的虎克定理，有 。xkf（ 恢复力与运动方向相反恢复力与运动方向相反 ）由牛顿第二定律，得由牛顿第二定律，得 dd22。xktxm例 )0( 0 dd222。，axatx

9、它能正确描述我它能正确描述我们的问题吗？们的问题吗？ 0 ，则有初始条件：，则有初始条件：t记拉长后，突然放手的时刻为记拉长后，突然放手的时刻为 00 ，初始位移初始位移xxt 0 dd 0 。初始速度初始速度ttx 2，则有，则有移项，并记移项，并记mka 我们要找的规律是下列初值问题的解：我们要找的规律是下列初值问题的解： 0 dd222，xatx 00 ，xxt。 0 dd 0 ttx 0 22，特征方程特征方程a i 2, 1 ，特征根a sin cos 21。所求通解为taCtaCx 0100 ；，得，得由由xCxxt 0 0) cos sin( dd 220 210 。，得，得由由

10、CaCtaaCtaaCtxtt从而，所求运动规律为从而，所求运动规律为 ) ( cos0。，mkataxx) , 0 (。a形如形如)2( )( xfyqypy )(常数。常数。实实为为的方程，称为二阶常系数线性非齐微分方程，的方程，称为二阶常系数线性非齐微分方程， qp、其中其中它对应的齐方程为它对应的齐方程为) 1 ( 0 。 yqypy我们只讨论函数我们只讨论函数 f ( x ) 的几种简单情形下，的几种简单情形下，(2) 的特解。的特解。三、二阶常系数非齐线性微分方程三、二阶常系数非齐线性微分方程0 yqypy二阶常系数齐线性方程二阶常系数齐线性方程)(xfyqypy 二阶常系数非齐线

11、性方程二阶常系数非齐线性方程特征方程特征方程02qp特征根特征根 , 21 2211yCyCy通解通解 * y特解特解 * yyy通解通解 二阶常系数线性微分方程解的关系二阶常系数线性微分方程解的关系)2( )( xfyqypy ) 1 ( 0 。 yqypy )()( . 1的情形的情形xPexfnx )( 1110。其中其中nnnnnaxaxaxaxP方程方程 (2) 对应的齐方程对应的齐方程 (1) 的特征方程及特征根为的特征方程及特征根为 0 2；特征方程特征方程qp 2 1 。，特征根单根单根二重根二重根一对共轭复根一对共轭复根 你认为方程应该你认为方程应该有什么样子的特解？有什么样

12、子的特解？假设方程假设方程)2( )(xPeyqypynx 有下列形式的特解：有下列形式的特解： )(，xueyx则则 ，ueueyxx 22，ueueueyxxx 代入方程代入方程 (2) ，得，得 )()()2(2，xPeuqpupuenxx 即即 )3( )()()2(2。xPuqpupun 方程方程 (3) 的系数与方程的系数与方程 (2) 的特征根有关。的特征根有关。qpqpp)2( )(xPeyqypynx )3( )()()2(2。xPuqpupun ) 1 (不是特征根，则不是特征根，则若若 02，qp由方程由方程 (3) 及多项式求导的特点可知，应有及多项式求导的特点可知，应

13、有 )()(1110，nnnnnbxbxbxbxQxu )2( )()( 的特征根时，的特征根时，不是方程不是方程中的中的故当故当xPexfnx方程方程 (2) 有下列形式的特解有下列形式的特解: )(*。xQeynx )(xueyx )2(是单特征根，则是单特征根，则若若 02，qp由多项式求导的特点可知，应有由多项式求导的特点可知，应有 )()()(1110，nnnnnbxbxbxbxxQxxu )2( )()( 的单特征根时，的单特征根时，是方程是方程中的中的故当故当xPexfnx方程方程 (2) 有下列形式的特解有下列形式的特解: )(*。xQexynx )3( 02 2 为为。此时，

14、方程。此时，方程，即，即而而pp )()2(。xPupun )2( )(xPeyqypynx )3( )()()2(2。xPuqpupun )(xueyx )3(是二重特征根，则是二重特征根，则若若 02，qp由多项式求导的特点可知，应有由多项式求导的特点可知，应有 )()()(111022，nnnnnbxbxbxbxxQxxu )2( )()( 的二重特征根时，的二重特征根时，是方程是方程中的中的故当故当xPexfnx方程方程 (2) 有下列形式的特解有下列形式的特解: )(*2。xQexynx )3( 0 2 2 为为。此时，方程。此时，方程，即，即且且pp )(。xPun )2( )(x

15、Peyqypynx )3( )()()2(2。xPuqpupun )(xueyx当二阶常系数线性非齐方程当二阶常系数线性非齐方程)2( )( xfyqypy )()( 时，时，的右端为的右端为xPexfnx它有下列形式的特解：它有下列形式的特解： )(*，xQexynxk其中：其中： 0 ；不是特征根时，取不是特征根时，取当当k 1 ；是单特征根时，取是单特征根时，取当当k 2 。是二重特征根时，取是二重特征根时，取当当k :。可以为复数可以为复数注意注意解解 2。的通解的通解求方程求方程xxyy ) )()( 2 0 )(2xPexfnxxxfnx。，对应的齐方程的特征方程为对应的齐方程的特

16、征方程为 012，特征根为特征根为 i2, 1 。对应的齐方程的通解为对应的齐方程的通解为 sincos21。xCxCy例 0 ，原方程有特解，原方程有特解不是特征根，故取不是特征根，故取由于由于k *2120，bxbxby将它代入原方程，得将它代入原方程，得 2221200，xxbxbxbb xxyy 2比较两边同类项的系数，得比较两边同类项的系数，得 10，b 11，b 0220，bb 10，b 11 ，b 2 2，b故原方程有一特解为故原方程有一特解为 2*2。xxy综上所述，原方程的通解为综上所述，原方程的通解为 2sincos*221。xxxCxCyyy解解 32 。的通解的通解求方

17、程求方程xeyyy ) )()( 0 1 )(xPexfnexfnxx。，对应的齐方程的特征方程为对应的齐方程的特征方程为 0322，特征根为特征根为 1 32 1 。，对应的齐方程的通解为对应的齐方程的通解为 231。xxeCeCy 1 ，原方程有特解，原方程有特解是单特征根，故取是单特征根，故取由于由于k *0，bexyx将它代入原方程，得将它代入原方程，得 3)1 (2)2(000，xxeexbxbxb例上式即上式即 140， b 410，b故原方程有一特解为故原方程有一特解为 41*。xexy综上所述，原方程的通解为综上所述，原方程的通解为 41*231。xxxexeCeCyyy将它代

18、入原方程，得将它代入原方程，得 3)1 (2)2(000，xxeexbxbxb解解 1332 。的通解的通解求方程求方程 xeyyyx 1332 xeyyyx 32xeyyy 1332 xyyy 41*1xexy31*2xy对应的齐方程的通解为对应的齐方程的通解为 231。xxeCeCy综上所述，原方程的通解为综上所述，原方程的通解为 3141*231。xexeCeCyyyxxx例请看上例请看上例 . )( ),( , sin)( cos)()( . 2次多项式次和分别是是常数，这里nlxQxPxxQxxPexfnlnlxsin)(cos)(*xxBxxAexymmxk特解：其中k; , 0不

19、是特征根i. , 1是特征根i.,max )()(nlmmxBxAmm次待定多项式，是，y+4y=cos2xii2, 2, 0特征方程 r2+4=0, 得r1,2=2i)2sin2cos(* xBxAxy故设41 , 0 BA解得xxy2sin41* 例解解 xxQxxPexfnlxsin)(cos)()(y+y=exsinxii1 , 1 , 1特征方程 r2+1=0, 得r1,2=i)sincos(e* xBxAyx故设51 ,52 BA解得)cos52sin51(e* xxyx故例解解xyyyxsine2652 求解特征方程 r25r60得 r1=2, r2=3xxCCy3221ee (

20、2) 再求y*：xyyy2e265 .sin65xyyy .sine2)( 2xxfx有(1)(2)(1) 先求y：例解解 求方程(1)的 y1*: 设y1*=Axe2x代入方程(1)得 A2xxy21e2*xyyy2e265 (1)r1=2, r2=3r25r60, 求方程(2)的y2*:xBxAysincos* 2设代入方程(2)得101 BA)sin(cos101* 2xxy* 21yyy)sin(cos10122xxxex)sin(cos101e2ee*23221xxxCCyyyxxx.sin65xyyy (2)r25r60, r1=2, r2=3解解 cos 的一个特解。的一个特解。

21、求方程求方程xyy 01 2，特征方程特征方程 i 2, 1 ，特征根 1 0 i ，且有，且有，故取，故取是特征根，是特征根，由于由于kn ，)sincos(*xbxaxy代入上述方程，得代入上述方程，得原方程有一特解为原方程有一特解为 。xxsin21*y例解解 sincos 的一个特解。的一个特解。求方程求方程xxxyy )cossin(41sin21*221xxxxxxyyy cos41sin432。xxxx思思 考考形如形如) 1 ( 01)1(1)(ypypypynnnn )(常数。常数。实实为为的方程，称为的方程，称为 n 阶常系数线性齐微分方程，阶常系数线性齐微分方程， , 1

22、npp 其中其中四、四、n 阶常系数齐线性微分方程阶常系数齐线性微分方程n 阶常系数线性齐微分方程的特征方程为阶常系数线性齐微分方程的特征方程为 单实根单实根xCe 1 项项 实重根k)( 121kkxxCxCCek项项 一对共轭复根一对共轭复根)sincos( 221xCxCex项项 011 1 nnnnpppi 2, 1 重复根重复根一对共轭一对共轭 ki 2, 1 2 项项k cos)(121xxCxCCekkx sin)(121xxDxDDkk特特 征征 根根通通 解解 中中 的的 对对 应应 项项解解 0dd3dd3dd 2233的通解。求方程yxyxyxy 0133 23，特征方程

23、特征方程 1 321，特征根特征根 ) ( 2321。所求通解为所求通解为xCxCCeyx例解解在研究弹性地基梁时，遇到一个微分方程在研究弹性地基梁时，遇到一个微分方程 )0( 0dd444。，x试求此方程的通解。试求此方程的通解。 0 44，特征方程特征方程 i)1 (2 i)1 (2 43 2, 1 ，特征根， 所求通解为所求通解为 ) 2sin2cos(212xCxCeyx ) 2sin2cos(432。xCxCex 2)(2222244例第六节 欧拉方程形如形如 )(1)1(11)(xfypyxpyxpyxnnnnnn的方程，称为的方程，称为 n 阶欧拉方程，其中阶欧拉方程，其中 )

24、, 2 , 1 ( 为常数。为常数。nipi tex 令令关于变量关于变量 t 的常系数线性微分方程的常系数线性微分方程 。引入算子记号：引入算子记号： ddkkktyyDkkktD dd ) , 2 , 1 (。k ，则有，则有令令tex dd1dddddd，tyxxttyxyy dddd1dd22222， tytyxxyy dd2dd3dd1dd2233333， tytytyxxyyyDyxyDDyx) 1(2 yDDDyx)2)(1(3 由数学归纳法可以证明：由数学归纳法可以证明： ) 1()2)(1()(。ynDDDDyxnn 例解解 34 223的通解。的通解。求方程求方程xyxyx

25、yx 这是三阶欧拉方程，这是三阶欧拉方程， ，原方程化为，原方程化为令令tex 34) 1()2)(1( 2，teDyyDDyDDD作代数运算后，得作代数运算后，得 332 223，teDyyDyD即即 (1) 3dd3dd2dd 22233，tetytyty这是一个三阶常系数线性非齐微分方程，且这是一个三阶常系数线性非齐微分方程，且 032 23，特征方程特征方程 , 3 1 0 3 2 1 ，特征根方程方程 (1) 对应的齐方程的通解为对应的齐方程的通解为 33 21。tteCeCCy 2 0 2 3)( 2特征根，故特征根，故不是不是，且，且，由于由于netftteby 20*为方程为方

26、程 (1) 特解形式，代入方程特解形式，代入方程 (1) 中，得中，得 ，tteebbb220003)688(从而从而 21* 21 20。，teyb故原欧拉方程的通解为故原欧拉方程的通解为ttteeCeCCyyy 2 33 2121* 21123321。xxCxCC 0 k故取.22的通解求方程xxyxyy 将方程化为xyyxyx22 (欧拉方程) ,ddtD 记则方程化为,tex 令teyDDD2)1) 1(即teyDD2) 12(2特征根:, 121 rr设特解:,2 tetAy 代入 解得 A = 1,ttetetCCy221)(xxxxCC221ln)ln(所求通解为 例解解满足设函数)(xyy 1,ln5d)(321 xxttytyyxx,01xy且.