二阶常微分方程的解法及其应用本科毕业.doc

《二阶常微分方程的解法及其应用本科毕业.doc》由会员分享，可在线阅读，更多相关《二阶常微分方程的解法及其应用本科毕业.doc（19页珍藏版）》请在三一文库上搜索。

1、 本科毕业论文 二阶常微分方程的解法及其应用 毕业论文（设计）原创性声明毕业论文（设计）原创性声明 本人所呈交的本人所呈交的毕业论毕业论文（文（设计设计）是我在）是我在导师导师的指的指导导下下进进行的研究行的研究 工作及取得的研究成果。据我所知，除文中已工作及取得的研究成果。据我所知，除文中已经经注明引用的内容外，本注明引用的内容外，本 论论文（文（设计设计）不包含其他个人已）不包含其他个人已经发经发表或撰写表或撰写过过的研究成果。的研究成果。对对本本论论文文 （ （设计设计）的研究做出重要）的研究做出重要贡贡献的个人和集体，均已在文中作了明确献的个人和集体，均已在文中作了明确说说明明 并表示

2、并表示谢谢意。意。 作者作者签签名：名： 日期：日期： 毕业论文（设计）授权使用说明毕业论文（设计）授权使用说明 本本论论文（文（设计设计）作者完全了解）作者完全了解*学院有关保留、使用学院有关保留、使用毕业论毕业论文（文（设计设计） ） 的的规规定，学校有定，学校有权权保留保留论论文（文（设计设计）并向相关部）并向相关部门门送交送交论论文（文（设计设计）的）的电电 子版和子版和纸质纸质版。有版。有权权将将论论文（文（设计设计）用于非）用于非赢赢利目的的少量复制并允利目的的少量复制并允许许 论论文（文（设计设计） ）进进入学校入学校图书馆图书馆被被查阅查阅。学校可以公布。学校可以公布论论文（文

3、（设计设计）的全部）的全部 或部分内容。保密的或部分内容。保密的论论文（文（设计设计）在解密后适用本）在解密后适用本规规定。定。 作者作者签签名：名： 指指导导教教师签师签名：名： 日期：日期： 日期：日期： 注 意 事 项 1.设计（论文）的内容包括： 1）封面（按教务处制定的标准封面格式制作） 2）原创性声明 3）中文摘要（300 字左右） 、关键词 4）外文摘要、关键词 5）目次页（附件不统一编入） 6）论文主体部分：引言（或绪论） 、正文、结论 7）参考文献 8）致谢 9）附录（对论文支持必要时） 2.论文字数要求：理工类设计（论文）正文字数不少于 1 万字（不包括图纸、程序清单等）

4、， 文科类论文正文字数不少于 1.2 万字。 3.附件包括：任务书、开题报告、外文译文、译文原文（复印件） 。 4.文字、图表要求： 1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准 请他人代写 2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符 合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒 手画 3）毕业论文须用 A4 单面打印，论文 50 页以上的双面打印 4）图表应绘制于无格子的页面上 5）软件工程类课题应有程序清单，并提供电子文档 5.装订顺序 1）设计（论文） 2）附件：按照任务书、开题报告、外文译

5、文、译文原文（复印件）次序装订 3）其它 目 录 1 引言5 2 二阶常系数常微分方程的几种解法.5 2.1 特征方程法5 2.1.1 特征根是两个实根的情形6 2.1.2 特征根有重根的情形.6 2.2 常数变易法.8 2.3 拉普拉斯变换法9 3 常微分方程的简单应用10 3.1 特征方程法11 3.2 常数变易法.13 3.3 拉普拉斯变换法14 4 总结及意义15 参考文献.16 二阶常微分方程的解法及其应用 摘要摘要：本文主要介绍了二阶常系数微分方程的三种解法：特征方程法、常数 变异法和拉普拉斯变换法，并着重讨论了特征方程根为实根、复根及重根的情形。 针对这三种解法的特点，分别将其应

6、用到求解弹簧振子系统的振子的运动方程。 关键词：关键词：二阶常微分方程;特征根法；常数变异法；拉普拉斯变换 METHODS FOR TWO ORDER ORDINARY DIFFERENTIAL EQUATION AND ITS APPLICATION Abstract:This paper mainly introduces three kinds of solution for two order differential equation with constant coefficients: the characteristic equation method, the method

7、of variation of constant and Laplasse transform method, and discusses the characteristics of Fang Chenggen is the real root, complex roots and root. According to the characteristics of the three solution, were applied to the equations of motion of vibrator for spring oscillator system. Keywords:seco

8、nd order ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform 1 引言引言 数学发展的历史告诉我们，300年来数学分析是数学的首要分支，而微分方 程又是数学分析的心脏，它还是数学分析里大部分思想和理论的根源。人所共知， 常微分方程从它产生的那天起，就是研究自然界变化规律、研究人类社会结构、 生态结构和工程技术问题的强有力工具。常微分方程已有悠久的历史，而且继续 保持着进一步发展的活力，主要原因是它的根源深扎在各种实际问题之中。常微 分

9、方程在很多学科领域内有着重要的应用，自动控制、各种电子学装置的设计、 弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等。 二阶常系数常微分方程在常微分方程理论中占有重要地位，在工程技术及力学和 物理学中都有十分广泛的应用。关于它的解结构己有十分完美的结论，但其求解 方法却各有不同，因此.二阶常系数线性微分方程的求解方法成为常微分方程研 究的热点问题之一。而本文正是在这一背景下对于二阶常系数常微分方程的解法 和应用做出研究。 2 二阶常系数常微分方程的几种解法二阶常系数常微分方程的几种解法 通常来说，纵观二阶常系数常微分方程的解法来看，其中比较有代表性的是 特征方程法、常数变

10、易法、拉普拉斯变换法这三种解法，因为篇幅和个人能力有 限，本文则选取这三种具备代表性的解法进行分析。 2.1 特征方程法 所谓特征方程，实际上就是为研究相应的数学对象而引入的一些等式，它因 研究对象的不同而不同，包括数列特征方程，矩阵特征方程，微分方程特征方程， 积分方程特征方程等等。 求微分方程 2 2 0 d xdx pqx dtdt 的通解. 解 特征方程 0 2 qp 的根 21, , (1)若这是两个不等实根,则该方程有两个实值解 12 , tt ee ,故通解为 12 12 tt xc ec e （ 21,c c 为任意常数）. (2)若这两个根相等,则该方程有二重根,因此方程的通

11、解具有形状 11 12 tt xc ec te （ 21,c c 为任意常数）. (3)若这两个根为共轭复根z abi ,则该方程的通解具有形状 12 (sincos) at xecbtcbt （ 21,c c 为任意常数）. 数学的许多公式与定理都需要证明,下面本文给出上面前两个解答的理论依据. 2.1.1 特征根是两个实根的情形 设 12 , 是上面特征方程的两个不相等的实根,从而相应的方程有如下两个解 12 , tt ee , 我们指出这两个解在a tb 上线性无关,从而它们能够组成方程的基本解组.事 实上,这时 12 12 12 () 12 12 11 ( ) tt t tt ee w

12、 te ee , 而最后一个行列式是著名的范德蒙德（Vandermonde）行列式,它等于 21 () .由 于假设 21 ,故此行列式不等于零,从而 ( )0w t ,于是 12 , tt ee 线性无关,这就是 所要证明的.而此方程的通解可表示为 12 12 tt xc ec e （其中 12 ,c c 为任意数）. 如果特征方程有复根,则因方程的系数是实常数,复根将成对共轭出现.设 1 i 是一特征根,则 2 i 也是特征根,因而与这对共轭复根对应的,方 程有两个复值解 () (cossin) itt eetit , () (cossin) itt eetit . 根据定理可知,复值解的

13、实部和虚部也是方程的解.这样一来,对应于特征方 程的一对共轭复根 i ,我们可求的方程 2 2 0 d xdx pqx dtdt 的两个实值解 cos,sin tt et et . 2.1.2 特征根有重根的情形 设特征方程有k重根 1, 则众所周知 (1) 111 ()()()0, k FFF ( ) 1 ()0 k F , 先设 1 0 ,即特征方程有因子 k ,于是 11 0 nnn k aaa , 也就是特征方程的形状为 1 1 0 nnk n k aa , 而对应的方程 1 11 1 0 nn nn nn d xdx L xaaa x dtdt 变为 1 1 1 0 nnk n k

14、nnk d ydyd y aa dxdxdx . 易见它有k个解1, 21 , k t tt ,而且它们是线性无关的.这样一来,特征方程的 k重零根就对应方程的k个线性无关的解1, 21 , k t tt .如果这个k重根 1 0 ,我 们作变量变换 1t xye ,注意到 11 ()()()(1)2(2) 111 (1) () 2! ttmmmmmm m m xyeeymyyy ， 可得 111 1 11 1 () nn ttt n nn d ydy L yebb y eLy e dtdt , 于是对应方程化为 1 11 1 0 nn n nn d ydy Lybb y dtdt , 其中

15、123 , n b b bb 仍为常数,而相应的特征方程为 1 11 ( )0 nn nn Gbbb , 直接计算易得 1111 ()()() 11 ()( ) ttttt FeL eLeeGe , 因此 1 ()( )FG , 从而 1 ()( ) jj FG ， 1,2,jk , 这样,问题就化为前面讨论过的情形了. 2.2 常数变易法 常数变易法是求解微分方程的一种很重要的方法,常应用于一阶线性微分方 程的求解。在常数变易法中，通过将常数 C 放入当中就可以得到非齐次线 XU 性方程的通解。它是拉格朗日十一年的研究成果，我们所用仅是他的结论，并无 过程。它是连接非齐次线性微分方程与相应的

16、齐次线性微分方程的桥梁。 对于二阶常系数非线性常微分方程的解法,只要先求出其一个特解,再运用特 征方程法求得方程的通解. 求常微分方程 2 2 ( ) d xdx pqxf t dtdt 的通解. 解 方程 2 2 ( ) d xdx pqxf t dtdt 对应齐次方程为 2 2 0 d xdx pqx dtdt , 其特征方程为 0 2 qp . 由于方程 2 2 ( ) d xdx pqxf t dtdt 的通解等于其对应的齐次线性微分方程的通解 与其自身的一个特解之和，而二阶常系数齐次线性微分方程的通解我们已经研究 过了,所以此处只需求出其一个特解. 若为上面方程的实根,则 t xe

17、是方程 2 2 0 d xdx pqx dtdt 的解.由常数变易 法设 2 2 ( ) d xdx pqxf t dtdt 的一个解为 * ( ) t xc t e ,代入原方程并化简得 “ ( )(2) ( )( ) t c tp c tef t , 这是关于 ( ) c t 的一阶线性微分方程,其一个特解为 (2)() ( )( ) p tp t c teef tdt dt , 从而得上面方程的一个特解为 *(2)() ( ) tp tp t xeeef t dt dt . 若为上面方程的复根,我们可以设 , ,abi a bR 且 0b ，则 * sin at xebt 是方程 2 2

18、 ( ) d xdx pqxf t dtdt 的解，根据常数变易法可设其一个特解 为 * ( )sin at xc t ebt ，与情形 1 的解法类似得方程 2 2 ( ) d xdx pqxf t dtdt 的一个特解 为 (2 )(2 ) * 2 ( )sin sin. sin papa t at ef t ebtdt xebtdt bt 由于 * x 是特解,则积分常量可以都取零. 2.3 拉普拉斯变换法 拉普拉斯变换法是工程数学中常用的一种积分变换法，又名拉氏转换法。拉 氏变换法是一个线性变换法，可将一个有因数实数的函数转换为 )0(tt 一个因数为复数 s 的函数。有些情形下一个实

19、变量函数在实数域中进行一些运算 并不容易，但若将实变量函数作拉普拉斯变换，并在复数域中作各种运算，再将 运算结果作拉普拉斯反变换来求得实数域中的相应结果，往往在计算上容易得多。 拉普拉斯变换的这种运算步骤对于求解线性微分方程尤为有效，它可把微分方程 化为容易求解的代数方程来处理，从而使计算简化。在经典控制理论中，对控制 系统的分析和综合，都是建立在拉普拉斯变换的基础上的。引入拉普拉斯变换的 一个主要优点，是可采用传递函数代替常系数微分方程来描述系统的特性。这就 为采用直观和简便的图解方法来确定控制系统的整个特性、分析控制系统的运动 过程，以及提供控制系统调整的可能性。常系数线性微分方程可以应用

20、拉普拉斯 变换法进行求解，这往往比较简单。 由积分 ( )( ) 0 st F sef t dt . 所定义的确定于复平面（Re ）上的复变数s的函数 ( )F s ,称为函数 ( )f t 的拉普 拉斯变换,我们称 ( )f t 为原函数,而 ( )F s 称为像函数. 拉普拉斯变换法主要是借助于拉普拉斯变换把常系数线性微分方程转换成复 平面s的代数方程.通过一些代数运算,一般地再利用拉普拉斯变换表,即可求出微 分方程的解.方法十分简单方便,为工程技术工作者所普遍采用.当然,方法本身有 一定的局限性,它要求所考察的微分方程的右端函数必须是原函数。 求解方程 2 2 2, (1)(1)0 t

21、d xdx xexx dtdt . 解 先使 1t ，将问题化为 2 (1) 2 2, (0)(0)0 t d xdx xexx dtdt , 再对新方程两边作拉普拉斯变换,得到 2 11 ( )2( )( ) 1 s X ssX sX s se , 因此 3 11 ( ) (1) X s se , 查拉普拉斯变换表可得 21 1 ( ) 2 xe , 从而 2 1 ( )(1) 2 t x tte , 这就是所要求的解. 当然,求解二阶或者更高阶的常微分方程的方法还有很多,这里我们不能一一 列出.然而我们利用上面的一些结论就可以解决下面的几个物理问题了。 3 常微分方程的简单应用常微分方程的

22、简单应用 为直观的了解常微分方程的简单应用，本文特选取动力学方程当中简单应用 常微分方程。通常来说，对于物理问题进行求解主要应该分为以下三个步骤内容： 第一步是对问题进行分析从而做到对方程的建立并且对定解条件进行明确；第二 步是对解的性质进行讨论或者求出方程以便满足初始条件的特解；第三步是定性 分析对解，对原来问题反着进行解释，其中最为关键的因素就是要将方程列出， 而列出方程的方法主要有：微元分析法和瞬时变化法。而在对阻尼振动进行研究 的过程当中，对运动方程所进行的求解这一问题显得比较复杂，以下就分别使用 特征值法、常数变异法以及拉普拉斯变换法来求动力学方程。 3.1 特征方程法 例如在弹簧振

23、子系统当中，测试出物体的阻尼系数 1 10.0s ，物体质量 1.0mkg ,该弹簧所具备的劲度系数 1 75kN m ,在此背景下，假设整个质点 从静止状态开始逐步运动，求解弹簧振子的位移方程。 解：按照牛顿的第二运动定律的结果可以得到 kxcvma , （1） 或 2 2 0 d xdx mckx dtdt , （2） 相对来说振动系统这是之前给定的，其中的常量为 , ,m k c，如果可以确定 2 0, 2k mc m ,那么以上的方程式可以转变为： 2 2 0 2 20 dd dtdt , （3） 那么把所得到的数据代入公式（3）就可以得到 2 2 20750 d xdx x dtdt

24、 . （4） 通过对以上公式的细致观察和研究则可以得到对其进行求解能够使用特征值 法，那么在这里的特征方程可以表述为： 2 20750 ,并且在这一特征方程 当中包含有两个分别根 12 15,5 ,这样相对应的则(4)的两个根分别为 515 12 , tt ee （5） 那么按照公式（5）进行计算可以得到固有角频率数值为 0 5 2k m ， 在这时候阻尼系数值为 10 ,也就是说 22 0 ,则方程(5)的解可以表述为 515tt AeBe （初始条件觉得 ,A B数值）. （6） 在公式（6）当中，所保持的属于一个非振动状态，在如此背景之下，所存 在的质点也只是在原先的不平衡位置逐步恢复到

25、平衡状态当中，质点并不具备周 期振动的特征。而我们的关注点是在基于 0 此种情况下,质点呈现出逐渐衰减 的振动。可是正是由于受到阻尼作用的影响，不能够长久的维持这种自由振动系 统的振动，通常都会经历着从振动的逐渐衰减延续至振动停止，为了保持震荡持 续不停的状态，就必须不断的从外界当中获得必要的能量，学术界将这种因为受 到外部持续作用而产生的振动归纳成为强迫振动。 又例如案例：假如在以上的振动系统当中受到某个外力 100cos(30 )Ft N 的 作用，在公式当中 100 A F 表示为驱动力所具备的幅度值， 30 则表示为驱动 力所拥有的圆频率， f 也就是驱动力所保持的频率。 解：在质点振

26、动系统当中受到驱动力的作用，那么就可以得到关于系统振动 的方程为： 2 2 d xdx mckxF dtdt , (7) 或者还可以将上述公式改成 2 2 0 2 2cos(30 ) d xdx xHt dtdt . (8) 在以上的公式当中 A F H m 表示为在单位质量上面所受到的外力幅值。 （7）与 （8）这两个方程式都属于质点强迫振动方程。从本质上来看，这种强迫振动方 程属于二阶的非齐次常微分方程，这个方程所得到的一般解也就是这个方程所得 到的某一个特解和相对应的齐次方程一般解两者之和。由于在之前的篇幅当中已 经得到相对应的自由振动方程的一般解，这就导致其在的关键问题就是对于 （8）

27、当中的一个特解进行寻找，把所得到的数据代入到（8）当中就可以得到： 2 2 2075100cos(30 ) d xdx xt dtdt , （9） 在这里可以通过假设（9）有着 1 sin30cos30xAtBt 这样的特解，将这个特别往 （9）当中进行替代并且将其进行简化之后得到 (3324 )sin30(2433 )cos304cos30ABtABtt , 按照比较同类项系数可以得到 3244 , 555555 AB ,这样就可以进一步得到 1 3244 sin30cos30 555555 xtt ,根据以上所得到的结果没那么原方程所存的通解就可 以表述为 515 3244 ( )sin3

28、0cos30 555555 tt x tAeBett . 在以上的公式当中，初始条件决定 ,A B的数值，而其中的瞬态解是之前的两 项，瞬态项能够对于整个系统的自由衰减振动进行有效描述，而所能够起作用的 只是在震动的开始阶段，而当经历比较长的时间之后，瞬态解所起到的影响则会 逐渐的减弱并且在最后阶段消失。稳态解则是之后的两项，稳态解则是对于系统 受到驱动力的作用之下进行强制振动的状态进行描述，这主要是由于立足于恒定 的幅值条件下，从而将这种状态称之为稳定振动。从以上的公式可以得到，如果 质点振动系统受到外力作用之后，整个系统有着比较复杂的振动状态，这属于稳 态振动和自由衰减振动两者的有机合成体

29、，在这样的振动状态之下对于强迫振动 当中逐步建立稳态振动的过程进行有效描述。如果经历一定时间之后，就会消失 瞬态振动，使得整个系统保持着稳态振动的状态。 3.2 常数变易法 从之前的分析当中可以了解到 5t xe 这属于特征方程 2 20750 的实根， 那么就可以得到 5t xe 这个属于方程（9）当中的一个根，然后通过常数变异法 设置 *5 ( ) t xc t e ，那么在这一过程当中也可以得到方程的一个解为 * x ，把数值代 入到（9）当中并且进行简化之后可以得到 “5 ( ) 10 ( )100cos30 t c tc tet . 以上属于 ( ) c t 的一阶线性微分方程,并且

30、在方程当中一个特解为 55 1 84 ( )sin30cos30 33 tt c tetetc , 从而得出（9）的一个特解为（取 12 0cc ） *555 12 84 ( )( (sin30cos30 ) 33 ttt x teetet dtc 3244 sin30cos30 555555 tt , 从而可得（9）的通解 515 3244 ( )sin30cos30 555555 tt x tAeBett . 由之前可知 2 2 d xdx mckxF dtdt . （10） 将数据代入公式中可以得到 2 2 20400cos(2 ) d xdx xt dtdt . （11） 按照自己所做

31、的观察可以发现，在进行求解的过程当中使用常数变异法，首要就 是必须得出公式（11） ，而在之前的研究当中可以得到公式（11）齐次线性微分 方程的特征方程为 2 204000 。这样就可以进一步的假设特征方程的根为 10 10 3i ，那么 10 ( )sin(10 3 ) t x tet 这就是公式（11）的一个解。由常数 变易法可设为 *10 ( )( )sin(10 3 ) t x tc t et . 与情形 1 中的解法类似，将 *( ) x t 代入（12）并化简得 * 1099 ( )sin(2 )cos(2 ) 3960439604 x ttt . 由于 * x 是特解,则积分常量

32、可以都取零。 3.3 拉普拉斯变换法 依然使用之前的例子，由牛顿第二运动定律可以得到以下的公式 2 2 d xdx mckxF dtdt , 将这一公式代入数据之后可以得到 2 2 20400cos(2 ) d xdx xt dtdt , (12) 由于质点通过开设的静止状态逐步运动，那么就可以得到以下的公式 0 0,0 t dx x dt , 对方程(12)进行拉普拉斯变换,得到 2 2 ( )20( )400( ) 4 s s X ssX sX s s , 即 22 1 ( ) 420400 s X s sss , 把上式右端分解为部分分式 22 10299 ( ) 39604439604

33、4 s X s ss 2222 101 310 39910 11881239604(10)(10 3)(10)(10 3) s ss , 由拉普拉斯变换表可得 1099 ( )sin(2 )cos(2 ) 3960439604 x ttt 1010 101 399 sin(10 3 )cos(10 3 ) 11881239604 tt etet 。 4 总结及意义总结及意义 总而言之，现在常微分方程在很多学科领域内有着重要的应用，自动控制、 各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反 应过程稳定性的研究等。现今对于二阶常微分方程解法的研究已经取得了不少成 就，尤其在

34、二阶常系数线性微分方程的求解问题方面卓有成效。而幂级数解法作 为求解二阶变系数齐次线性微分方程的一种方法，其过程还是比较繁琐的，计算 量偏大，且需要考虑函数是否解析，幂级数在某个区间是否收敛等。另外，对于 二阶常系数非齐次线性微分方程，目前还尚有通用的求解方法，只有一些特殊类 型是可以求解的。应该说，应用常微分方程理论已经取得了很大的成就，但是， 它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理 论更加完善。 参考文献参考文献 1 (瑞典)L.戈丁(Lars.Garding)著,胡作玄译.数学概观M.科学出版社, 2001：121-147 2 赵慈庚,朱鼎勋主编.大学数学

35、自学指南M.中国青年出版社,1984:74- 91 3 王高雄等编.常微分方程M.高等教育出版社,1978:39-53 4 李瑞遐,何志庆编著.微分方程数值方法M.华东理工大学出版社, 2005：41-58 5 余德浩,汤华中编著.微分方程数值解法M.科学出版社,2003：14-21 6 胡燧林.一阶方程初值问题解的存在与唯一性定理的几点注记J.韶关 学院学报.1988(02)：38-47 7 弭鲁芳,纪在秀.论一类常微分方程解的最大存在区间J.聊城大学学报 (自然科学版).2006(04):24-26+28 8 赵慧娟,陈伟丽,赵晨霞,袁书娟.关于常微分方程初值问题数值解法的分 析J.中国科教创新导刊.2011(08):83 9 李孝诚,刘兆丽.常微分方程解题模式的构建J.高等数学研究.2009(04)： 96-99 10 韦程东,高扬,陈志强.在常微分方程教学中融入数学建模思想的探索与 实践J.数学的实践与认识 2008(20)：228-233 11 舒小保.一类二阶常微分方程边值问题的无穷多个解J.系统科学与数 学.2008(01)：91-98