matlab在求解常微分方程中的应用.doc

《matlab在求解常微分方程中的应用.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、《MATLAB语言》课程论文MATLAB在求解常微分方程中的应用姓名：袁学婷学号：12010245278专业：通信工程班级：通信一班指导老师：汤全武学院：物理电气信息学院完成日期：2012年12月10日MATLAB在求解常微分方程中的应用（袁学婷120102452782010通信班）【摘要】MATLAB成为许多学科的解题工具，将MATLAB融入其它课程的学习中，可以大大提高运算效率和准确性。随着计算机的普及和国民整体素质的提高，科学计算将会更加的普及。MATLAB在矩阵及数值计算、多项式和线形代数、符号数学的基本方法等

2、方面都有较好的应用。本文概括地介绍了MATLAB在求解常微分方程中的应用。【关键词】MATLAB求解常微分方程解析解数值解一、问题的提出自20世纪80年代以来，出现了多种科学计算语言，亦称数学软件，比较流行的有MATLAB、Mathematica、Maple等。因为他们具有功能强、效率高、简单易学等特点，在在许多领域等到广泛应用。MATLAB便是一种影响大、流行广的科学计算语言。MATLAB的语法规则简单，更加贴近人的思维方式。MATLAB是英文MATrixLABoratory(矩阵实验室)的缩写。自1984年由美国M

3、athWorks公司推向市场以来，得到了广泛的应用和发展。在欧美各高等院校MATLAB已经成为线性代数、自动控制理论、数字信号处理、时间序列分析、动态系统仿真、图像处理等诸多课程的基本教学工具，成为大学生、硕士生以及博士生必须掌握的基本技能。在设计研究单位和工业部门，MATLAB已被广泛的应用于研究和解决各种具体的工程问题。近年来，MATLAB在我国也开始流行，应用MATLAB的单位和个人急剧增加。可以预见，MATLAB将在我国科学研究和工程应用中发挥越来越大的作用。现在简单的介绍一下MATLAB在解微分方程中的应用。

4、众所周知，只有对一些典型的常微分方程，才能求出它们的一般解表达式并用初始条件确定表达式中的任意常数。所以能够求解的微分方程是十分有限的，特别是高阶方程（组）．这就要求我们必须研究微分方程（组）的解法，既要研究微分方程（组）的解析解法（精确解），更要研究微分方程（组）的数值解法（近似解）．对微分方程（组）的解析解法(精确解)，Matlab有专门的函数可以用，本文将作一定的介绍．二、常微分方程的概述1、微分方程的概念未知的函数以及它的某些阶的导数连同自变量都由一已知方程联系在一起的方程称为微分方程。如果未知函数是一元函数，

5、称为常微分方程。常微分方程的一般形式为0),,",',,()(=nyyyytF...(1)2、常微分方程的解析解有些微分方程可直接通过积分求解.例如,一解常系数常微分方程可化为,两边积分可得通解为.其中为任意常数.有些常微分方程可用一些技巧,如分离变量法,积分因子法,常数变异法,降阶法等可化为可积分的方程而求得解析解.线性常微分方程的解满足叠加原理,从而他们的求解可归结为求一个特解和相应齐次微分方程的通解.一阶变系数线性微分方程总可用这一思路求得显式解。高阶线性常系数微分方程可用特征根法求得相应齐次微分方程的基本解，再

6、用常数变异法求特解。一阶常微分方程与高阶微分方程可以互化，已给一个阶方程(2)设，可将上式化为一阶方程组(3)反过来，在许多情况下，一阶微分方程组也可化为高阶方程。所以一阶微分方程组与高阶常微分方程的理论与方法在许多方面是相通的，一阶常系数线性微分方程组也可用特征根法求解。3、微分方程的数值解法除常系数线性微分方程可用特征根法求解，少数特殊方程可用初等积分法求解外，大部分微分方程无限世界，应用中主要依靠数值解法。考虑一阶常微分方程初值问题(4)其中所谓数值解法，就是寻求在一系列离散节点上的近似值称为步长，通常取为常量。

7、最简单的数值解法是Euler法。Euler法的思路极其简单：在节点出用差商近似代替导数(5)这样导出计算公式（称为Euler格式）(6)他能求解各种形式的微分方程。Euler法也称折线法。4、解微分方程的MATLAB命令MATLAB中主要用dsolve求符号解析解,0de45,0de23,ode15s求数值解。在MATLAB中，由函数dsolve()解决常微分方程（组）的求解问题，其具体格式如下：r=dsolve('eq1,eq2,...','cond1,cond2,...','v')'eq1,eq2,...'为微分方

8、程或微分方程组，'cond1,cond2,...',是初始条件或边界条件，'v'是独立变量，默认的独立变量是't'。函数dsolve用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解。基于龙格-库塔法，MATLAB提供了求常微分方程的数值解的函数，一般调用格式为：[t,y]=ode23(file