常微分和matlab.docx

本科毕业论文（设计）模板本科中期论文(设计)论文题目：用MATLAB解常微分方程学生姓名： 李涛 学 号：0804280134专 业：数学与应用数学班 级：数学0801 指导教师： 朱少平 完成日期：2011年04月 24 日目 录内容摘要………………………………………………………………1序言……………………………………………………………………2常微分方程………………………………………………………2（一）常微分方程的基本概念………………………………………2（二）常微分方程的解 ………………………………………………3（三）常微分方程的几种典型形式…………………………………4二、MATLAB概述 ………………………………………………………5三、用MATLAB求解常微分方程………………………………………5（一）常微分方程的符号解……………………………………………5（二）常微分方程的数值解法…………………………………………7四、总结………………………………………………………………20参考文献………………………………………………………………21用MATLAB解常微分方程内容摘要在本文中，我们讨论了MATLAB在微分方程中的应用，包括常微分方程的符号求解方法和数值求解方法以及偏微分方程的有限元方法。我们介绍了MATLAB求解微分方程的有关函数命令，为了求解方程根据算法编写了实现算法的M文件，运行结果与理论估计是一致的。关键词： 常微分方程；MATLAB； 欧拉方法； 泰勒方法。序言现代的科学、技术、工程中的大量数学模型都可以用微分方程来描述，很多近代自然科学的基本方程本身就是微分方程。从微积分理论形成以来，人们就一直用微分方程来描述、解释或预见各种自然现象，不断地取得显著成效。遗憾的是，绝大多数微分方程（特别是偏微分方程）定解问题的解不能以实用的解析形式来表示，这就产生了理论与应用的矛盾：一方面，人们建立了大量实用的数学模型，列出了反映客观现象的微分方程；另一方面，人们又无法得到这些方程的准确解以定量地描述客观过程。随着电子计算机的出现和发展，解决上述矛盾的一门学科---微分方程的数值方法得到了前所未有的发展和应用。今天，需要掌握和应用微分方程数值方法的已不再限于数学专业人员和专门从事科学和工程计算的人员，大量从事力学、物理学、天文学方面的科研人员，航空、土木、船舶、机械、电机、电子等领域的工程技术人员，甚至连金融工程、风险投资等非科技领域都把这门学科作为自己领域的一种重要研究手段和方法。从20世纪80年代以来，世界各国所开发的数学类科技软件多达几十种，在我国流行的数学软件主要有四种：MATLAB、Mathematica、Maple和MathCAD。其中，MATLAB有着其它几种数学软件无法比拟的优势和适用面，近几年，MATLAB已成为科学工作者首选的数学软件。目前关于MATLAB语言和应用的书籍在国际上数以千计，但从其覆盖面和应用水平来说，往往难以达到日益增长的MATLAB使用者的要求。国内外出版的著作缺乏全面系统介绍MATLAB语言在应用数学的各个分支的应用、计算和求解的书籍。本文旨在填补就微分方程方面的空白，重点是微分方程数值解法方面，以帮助微分方程求解工作者掌握MATLAB在这方面的功能，避免重复性的劳动，尽快地“站在巨人肩上”开展工作。本文把微分方程求解和MATLAB有机的结合起来，全面介绍了微分方程的求解在MATLAB中的实现，使得让数学基础不深厚的读者同样能轻易利用MATLAB解决较高深的微分方程问题。常微分方程（一）常微分方程的基本概念定义1凡含有未知函数导数（或微分）的方程，称为微分方程，有时简称为方程，未知函数是一元函数的微分方程称做常微分方程，未知函数是多元函数的微分方程称做偏微分方程。这里仅讨论常微分方程，并简称为微分方程。例1已知一曲线通过点（1，3），且在该曲线上任意点M处的切线的斜率为，求这曲线的方程。解根据导数的几何意义，可知所求曲线应满足方程． (１)将方程(1)两边积分，得，(２)其中C是任意常数。因为曲线通过点，将条件代入（3）式，得．故得所求曲线方程为．（３）例2 一物体以初速垂直上抛，设此物体的运动只受重力的影响，试确定该物体运动的路程与时间的函数关系。解因为物体运动的加速度是路程对时间的二阶导数，且题设只受重力的影响，所以由牛顿第二定律有即 (４)其中物体的质量为，重力加速度为，且垂直向上的方向为正方向。因为物体的运动速度，所以 (４) 式可写为或两边积分，得或　　 (５)再对上式积分，得(６)其中为任意常数．如果假设物体开始上抛时的路程为，则由题意有，代入（5）、(６)式，得。于是所求的函数关系为:上述两个例子中的方程（