高阶微分方程的解法及应用毕业论文

哈尔滨学院本科毕业论文（设计）本科毕业论文(设计)题目：高阶微分方程的解法及应用哈尔滨学院本科毕业论文（设计）毕业论文（设计）原创性声明本人所呈交的毕业论文（设计）是我在导师的指导下进行的研究工作及取得的研究成果。据我所知，除文中已经注明引用的内容外，本论文（设计）不包含其他个人已经发表或撰写过的研究成果。对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中作了明确说明并表示谢意。作者签名：日期：毕业论文（设计）授权使用说明本论文（设计）作者完全了解**学院有关保留、使用毕业论文（设计）的规定，学校有权保留论文（设计）并向相关部门送交论文（设计）的电子版和纸质版。有权将论文（设计）用于非赢利目的的少量复制并允许论文（设计）进入学校图书馆被查阅。学校可以公布论文（设计）的全部或部分内容。保密的论文（设计）在解密后适用本规定。

作者签名：指导教师签名：日期：日期：

注意事项1.设计（论文）的内容包括：1）封面（按教务处制定的标准封面格式制作）2）原创性声明3）中文摘要（300字左右）、关键词4）外文摘要、关键词5）目次页（附件不统一编入）6）论文主体部分：引言（或绪论）、正文、结论7）参考文献8）致谢9）附录（对论文支持必要时）2.论文字数要求：理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。4.文字、图表要求：1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印4）图表应绘制于无格子的页面上5）软件工程类课题应有程序清单，并提供电子文档5.装订顺序1）设计（论文）2）附件：按照任务书、开题报告、外文译文、译文原文（复印件）次序装订3）其它

目录摘要 摘要本文首先介绍了高阶微分方程的一些理论与结构。进而介绍了高阶齐次线性微分方程的求解方法和高阶非齐次线性微分方程的求解方法，在求解齐次线性微分方程里主要采用了特征根法；在求解非齐次线性微分方程里主要采用了比较系数法、拉普拉斯变换法和常数变易法。其次又介绍了几类可降阶的微分方程的解法，主要有形如，，，恰当导数方程和Euler方程的降阶方法，并且研究了几类较为复杂的高阶微分方程的降阶问题。最后通过一些在现实生活中例子对这些方法的具体应用做了介绍。关键词：高阶常微分方程；常数变易法；特征根法；降阶法AbstractThispaperintroducessomeofthetheoriesandhigherorderdifferentialstructure.Thenintroducehigher-orderhomogeneouslineardifferentialequationmethodsandhigh-ordernon-homogeneouslineardifferentialequationmethodforsolvinghomogeneouslineardifferentialequationwherethemainuseoftheeigenvaluemethod;insolvinginhomogeneouslineardifferentialequationsinmainlyusesthecomparisoncoefficientmethod,Laplacetransformmethodandtheconstantvariation.Andsecondlydescribesseveraltypesofdifferentialequationscanbereducedforthesolution,themaintangibleeg,appropriatederivativeequationsandEulerequationsreductionmethod,andstudiedseveraltypesofmorecomplexhigherorderdifferentialequationsreductionproblem.Finallysomereallifeexamplesofspecificapplicationsofthesemethodshavebeendescribed.Keywords:HigherOrderOrdinaryDifferentialEquations;constantvariation;eigenvaluemethod;reductionmethod

前言常微分方程作为数学系重要专业的一门基础课程，对学习好其他的科目起到了至关重要的作用。它的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的。数学的其他分支的新发展，如复变函数、拓扑学等，都对常微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具。而高阶微分方程是常微分方程中的一个重要的组成部分，在现实的生活中也有着广泛的应用，比如工程问题。常系数线性微分方程的解法，高阶微分方程的降阶问题又是高阶微分方程的重中之重。常微分方程是在生产实践和科学技术中产生的。目前，常微分方程在很多学科领域内有着重要的应用，自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等。这些问题都可以化为求常微分方程的解，或者化为研究解的性质的问题。人们对于二阶以及简单的高阶微分方程求解的方法有了很多理论成果，而高阶常微分方程并没有固定的解法，例如，高阶常系数线性齐次微分方程，我们可以运用特征根的方法进行求解，高阶常系数线性非齐次微分方程，我们可以运用常数变易法，比较系数法，拉普拉斯变换法进行求解。而对于可以降阶的高阶微分方程，我们通常采用降阶法，也就是通过一定的变换把高阶微分方程求解的问题转化成低阶微分方程的求解问题。本篇论文我总结了形如，，，恰当导数方程和Euler方程的降阶方法，并且研究了几类较为复杂的高阶微分方程的降阶问题，进而介绍此类问题在科学技术中的应用。第一章高阶微分方程的理论与结构定义1（方程的阶）在一个常微分方程里，未知函数的最高阶导数的阶数叫做方程的阶。n阶隐式方程的一般形式为n阶显式方程的一般形式为定义2（解）设函数在区间上有直到阶的导数。如果把代入到方程得到在区间上关于的恒等式是则称是方程在区间上的一个解。微分方程的解可以包括任意的常数，其中任意常数的个数可以多到和方程的阶数相等，当然也可以不包括任意常数。我们把方程的含有个独立的任意常数的解称做该方程的通解。如果方程的解不包含任意常数，则把它叫做特解。方程（1-1）称做n阶线性微分方程，它关于未知函数以及各阶导数都是线性的。在这里，我们通常假设和是区间上的连续函数。如果都是常数，则把方程（1-1）叫做n阶常系数线性方程。如果方程的右端项，即则称方程（1-1）是齐次的，否则为非齐次的。所以对于方程（1-1）的齐次方程是（1-2）定理1（叠加原理）设和是齐次方程（1-2）的解，则对于任意常数和，也是方程（1-2）的解。定理2设是方程的解，则也是方程的解。定理3设是齐次方程（1-2）的n个线性无关的特解，则是方程（1-2）的通解，其中是任意常数。定理4设是非齐次线性方程（1-1）的任意一个确定的解，是（1-1）对应的齐次线性方程（1-2）的通解。则是（1-1）的通解。第二章高阶常系数线性微分方程2.1高阶常系数线性齐次微分方程对于n阶常系数线性齐次方程（2-1）其中是关于的未知函数，系数是实常数。如果是方程的根，把他代入到方程中，得因为，因此（2-2）反之，如果满足等式（2-2），则是方程(2-1)的解。式子（2-2）是关于的n次代数方程，则把他叫做微分方程（2-1）的特征方程，它的根就称做特征根。下面根据特征根的不同情形分别进行讨论方程解的情况。2.1.1特征根是单根的情况定义我们把称为方程的特征方程，它的根叫做特征根。在这里把叫做待定系数。定理如果特征方程有个互异的根，则是方程的一个基本解组。特征方程可能有复根，由于他的系数是实的，他的复根一定是共轭成对的出现。即此时在相异特征根中有复数。例如，则也是的根。这两个特征根所对应的解是实变量复值函数例1求方程的通解。解特征方程的根是，其中有两个实根和两个复根，但他们都是单根，所以所求方程的通解是在这里是任意的常数。2.1.2特征根是重根的情况定理假设方程有互异的特征根，他们的重数分别是，并且，则与他们相对应的的特解是，并且该特解构成在区间上的基本解组。例2解初值问题解特征方程是，特征根是所以方程的通解是又因为根据初始条件，得再解方程组，得于是初值问题的解是2.2高阶常系数线性非齐次方程对于n阶常系数线性非齐次方程（2-3）他的通解等于齐次方程的通解再加上加其对应的非齐次方程的一个特解。在上一节中我们知道了怎样求解齐次方程的通解，下面我们主要来研究求解非齐次方程的特解的方法。2.2.1常数变易法常数变易法实际上是一种变量变换的方法,在这里我们简单的介绍一下在n阶方程中的应用。可以设方程（2-3）的特解是：(2-4)其中是待定的常函数。并且把它代入到方程(2-3)中,再附加上n-1个条件,就可以得到方程组（2-5）解方程组(2-5)就会得到关于的表达式,把它们分别进行积分进而得到,再把它们代入到(2-4)式中,继而求得方程（2-3）的一个特解。由于这种方法对于自由项的形式没有任何的限制,因此使用的范围会比较广,但是求解的工作量相对来说会大一些。例3求解方程的通解,已知它所对应的齐次线性微分方程的基本解组是。解运用常数变易法,设并且把它代入到方程里,就可以得到关于和的两个方程和解得,据此得到所以原方程的通解是其中是任意的常数。2.2.2比较系数法对于常系数非线性方程（2-3）,我们通常采用的方法是比较系数法,它是把所要求解的微分方程的问题转化成代数问题,在自由项是(其中分别是次,次,次的多项式。都是实常数)时,就可以确定特解的形式,即分别令是一个待定的次的多项式,是方程的特征方程有根时的次数)或者(其中是两个待定的次多项式,是方程含有根的次数)然后把它代入到方程（2-3）中,再进行比较等式的左右两边同次幂的系数来确定待定系数多项式。再根据线性微分方程解的结构便可以求解出方程的通解。例4求方程的通解。解特征方程有三重根,所对应的齐次方程的通解是并且方程有的特解,将它代入到方程中得再比较两边的系数求得进而所以所求的方程的通解是其中是任意的常数。2.2.3拉普拉斯变换法根据积分所定义的确定在复平面（）上的复变数的函数，叫做函数的拉普拉斯变换，其中在上有定义，并且满足不等式在这里是某两个正常数。我们把称为原函数，而把称为象函数。设所给定的微分方程（2-6）和初始条件其中是常数，而是连续的并且满足原函数的条件。可以证明，假如是方程（2-1）的任意解，则以及它的各阶的导数都是原函数。记那么，根据原函数的微分性质就有于是，再对方程（2-6）的两边进行拉普拉斯变换，并且运用线性性质就可以得到即或者其中和都是已知的多项式，由此得到这就是方程（2-6）所满足所给初始条件的解的象函数。而可以直接查表或者根据反变换公式计算求解出来。例5求方程的满足初始条件的解。解对方程的左右两边进行拉普拉斯变换得到由此得到再把上面式子的右面分解成为部分分式对上面式子的右端的各项分别求出或者查表得出他们的原函数，则他们的和就是的原函数这就是所要求的解。2.3Euler方程定义：形如的方程叫做Euler方程，其中是实常数，并且。它的特点就是包含的阶导数项的系数是。当时，各阶导数项的系数是0，所以我们令。在现实里，我们仅需要考虑的这种情况，因为在的时候，在上述的方程里做自变量变换，则方程就化成求出他的解，再用替换就可以得出方程关于的解。再做自变量变换则，一般的，假如，其中是常数，则所以，对于每一个正整数是的常系数的线性组合，进而把Euler方程化成了常系数线性方程。例5求解方程解设，方程化做他的特征方程是他的特征值是-2和，方程的通解是第三章可降阶的高阶微分方程的解法本部分我将介绍4类比较常见的高阶微分方程的解法，在这些解法里有一个比较类似的思路，就是把这些的高阶微分方程通过某些变换降成比较低阶的微分方程再进行求解。所以，我们把这种方法称做“降阶法”。3.1形如的高阶方程方程（3-1）这种类型的方程比较简单，通常令，则积分得也就得到同理可以令得到如此继续下去，再通过次积分就可以求出（3-1）的通解是例1求解微分方程的通解。解对原方程的左右两边依次进行积分，得再次进行积分，求解出原方程的通解是例2求方程的通解解所以所求原方程的通解是3.2形如的高阶方程方程（3-2a）这种方程的特点极其容易看出，方程不显含或者在这时我们只要把代入到上述的方程中，原方程就可以化作（3-2b)如果方程（3-2b）可以求解出通解.则再对方程积分次，便可以求出了。在这里需要注意的是每积分一次，就需要增加一个独立的任意常数.例1求解方程解设则代入到上面的方程中得再积分得即积分四次就可以求解出原方程的通解例2求微分方程的通解解设则。所以原方程就可以写成左右两边进行积分得所以也就是，两边再次进行积分得出3.3形如的高阶方程方程（3-3）这类方程也有一定的特点，就是不显含自变量，这时，我们总可以利用代换使这类方程降低一阶。以二阶方程为例，设于是便有代入到原方程中，就有这是一个关于未知函数的一个一阶方程。如果用它可求出就有这是一个关于的变量可分离方程，进而可以求解出通积分。例1求解方程解根据上面的分析，我们可以令则代到原方程里得即左右两边进行积分得求解出得积分后得于是便有例2求解方程解首先令则于是原方程就化成了再令则即进而即进而可以积分求解出通积分3.4恰当导数方程假如方程（3-4）的左面正好是某一个函数对的导数，则（3-4）就可化为于是我们就把（3-4）称作恰当导数方程。其实这类方程的解法和全微分方程的解法很类似，可以降低一阶，化成之后，再想办法求解这个方程。例1求解方程解我们可以把方程写成所以就有即积分后就可以得出通积分这样的问题虽然简单，但是需要具有很强的观察能力和比较牢固的基础才可以观察出来。下面有一个关于这方面的例子，解法技巧很高明，关键还是配导数的方法。例2求解方程解经过观察我们可以先把等式的两边同时乘以一个不是的因子便有所以从而通解是第四章高阶微分方程的应用要利用微分方程解决实际问题，首先必须要根据物理和几何关系规律来建立微分方程，然后再对进一步的问题进行分析与微分方程的建设，并且考虑初始条件，边界条件，收敛条件来确定定解的条件，这是数学建模过程。模型建立好了就有了微分方程，我们就可以根据前面的内容来解除方程，因为解决的是实际问题，我们还要用解出来的结果来分析问题。这部分内容因为实际应用相对比较强，所以我用三个简单的例子来简单的介绍一下。例1设质量是的物体自由悬挂在一个一端固定的弹簧上，当重力跟弹簧力相互抵消的时候，物体就会处在一个平衡的状态，若用手向下拉物体使它离开平衡的位置以后放开，物体在弹性力与阻力的作用下将做往复运动，阻力的大小与运动的速度成正比，方向相反。（1)建立位移所满足的微分方程解设时刻物体的位移是。1.自由振动的情况，物体所受到的力有弹力恢复力阻力根据牛顿第二定律得到令就可以得到阻尼自由振动方程2.强迫振动情况，如果物体在运动的过程里还受到铅直外力的作用，设得到强迫振动的方程(2)在没有外力的作用下做自由运动，设时物体的位置是，初始的速度为，求物体的运动规律解由（1）知，位移满足的定解问题是1.无阻尼自由振动的情况方程解得方程的通解是再利用初始条件可以得到所以所求的特解是1有阻尼自由振动情况,方程特征方程是特征值在这个时候我们需要三种情况来进行讨论小阻尼则大阻尼则临界阻尼则例2人类将要向宇宙发射一颗人造地球卫星，为了让她摆脱地球引力，初始速度应该不少于第二宇宙速度，试求该速度。解在物理问题中，关键是要通过建立模型，把物理问题转化成数学问题，在这个题目里设人造地球卫星的质量是，地球的质量是,卫星的质心到地心的距离是，根据牛顿第二定律得(是引力系数）再设卫星的初速度是，已知地球的半径于是就有初值问题.于是以上的物理问题就转化成了求二阶常微分方程的特解的问题，设代入到上述方程组的第二个式子就可以得到从而就有两边进行积分得到再利用初始条件得所以注意到为了让应该满足(4-1)因为在h=R（在地面上）时，引力跟重力是相等的，即所以代入到方程(4-1)中得这就说明第二宇宙速度是例3在船上向海里沉放某一种探测器，按照探测的要求，需要确定仪器的下沉深度和下沉的速度之间的函数关系。假设仪器在重力的作用下在海平面由静止开始往下沉，在下沉的过程中还受到了阻力和浮力的作用，我们设仪器的质量是，体积是，海水的比重是，仪器所受到的阻力跟下沉的速度成正比，比例系数是试着建立与所满足的微分方程，并且求出函数关系式解同样也是首先把实际的问题转化成常微分方程，根据题目中的条件和牛顿第二定律可以将问题转化成求解初始条件是的二阶微分方程的特解的问题。我们有得初始条件是再利用分离变量法解上述初值问题得参考文献[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.TestingforLanguageTeachers.Cambridge:CambridgeUniversityPress,1989（期刊）[15]Nunan,D.TheLearner-CentredCurriculum.Cambridge:CambridgeUniversityPress,1988（期刊）致谢作为即将从哈尔滨学院理学院毕业的我，在四年的大学生活中，认真学习了各科专业知识，积极参加了社会实践活动。在大四的师范实习期间，我的教学方面有了显著的提高，特别是在教态、教学方法、教学过程中与学生的沟通技能方面有了明显的改进。回首大学四年时光，匆匆而过，我要诚挚的感谢教育和培养我的老师们，感谢徐亚兰老师对我完成论文的选题、撰写方面给予的指导和帮助.在撰写论文期间，徐老师对我的谆谆教导、和她的严谨的治学态度使我终身受益，对我未来参加工作必将产生深远的影响。我真诚感谢哈尔滨学院理学院的各位领导和老师等给我长期的传道授业解惑，对本论文在撰写过程中给我的知识指导、帮助和启迪。借此机会向在大学四年里给予我的帮助和关怀的同学们表示感谢。论文中还有诸多的问题和不足之处，敬请大家给予批评指证。基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究基于TCP/IP协议的单片机与Internet互联的研究与实现变频调速液压电梯单片机控制器的研究基于单片机γ-免疫计数器自动换样功能的研究与实现基于单片机的倒立摆控制系统设计与实现单片机嵌入式以太网防盗报警系统基于51单片机的嵌入式Internet系统的设计与实现单片机监测系统在挤压机上的应用MSP430单片机在智能水表系统上的研究与应用基于单片机的嵌入式系统中TCP/IP协议栈的实现与应用单片机在高楼恒压供水系统中的应用HYPERLINK