微分方程在MATLAB中的实现.pdfVIP

MATLAB MATLAB 微分方程在 MMAATTLLAABB中的实现 作者：吴建宏 2011.09.20 2011.09.20 时间.0099..2200 ********************************************* 作者简介： 吴建宏，男，毕业于哈尔滨工业大学（威海），现就 读于同济大学，攻读汽车电子方向，有两年的 MATLAB 实 践经验，个人喜欢建模，编程和电控方面的。真诚愿意和 各位志同道合的朋友一起探讨交流，一起搭建广阔的知识 平台。 ********************************************* PARTONE PARTONE PPAARRTTOONNEE 微分方程简介 一、微分方程基本概念 微分方程：未知函数以及它们的某些阶的导数连同自变量都由一已知方程联系在 一起的方程称为微分方程。其中，未知函数的最高阶数称为微分方程的阶。 偏微分方程：如果未知函数是多元函数，那么这个微分方程就是偏微分方程。 常微分方程：如果未知函数是一元函数，那么这个微分方程就是常微分方程。 二、微分方程的解法 高等数学里讲过，微分方程的解法有分离变量法，常数变易法，特征值法……而 在解决工程实际问题时候，利用 MATLAB 求解微分方程的方法主要有符号解法和数值 解法。其中，符号解法主要可以求解可以用特征值法求解的常系数微分方程和少数特 殊的微分方程，有很大的局限性，而且运行时间较长，优点就是可以得到精确的数学 表达式。大部分实际的微分方程不得不通过数值解法来求解，优点就是运行时间快， 可以解决复杂的线性或者非线性微分方程。缺点就是解是近似值。 PARTTWO dsolve PARTTWO dsolve PPAARRTTTTWWOO 微分方程的符号解法（ddssoollvvee） dsolve dsolve 微分方程的符号解法主要是函数 ddssoollvvee，这个函数用起来很简单，最重要的是要 dsolve dsolve 知道神马时候能够用 ddssoollvvee，神马时候不能用，当运行出错的时候，该怎么处理。接 下来进入正文。 一、语法 r= dsolve(eq1,eq2,..., cond1,cond2,..., v) r= dsolve(eq1,eq2,..., cond1,cond2,..., v) rr==ddssoollvvee((eeqq11,,eeqq22,,,, ccoonndd11,,ccoonndd22,,,,vv)) r= dsolve(eq1,eq2,...,cond1,cond2,...,v) r= dsolve(eq1,eq2,...,cond1,cond2,...,v) rr==ddssoollvvee((eeqq11,,eeqq22,,,,ccoonndd11,,ccoonndd22,,,,vv)) 二、语法详解 1. eq1,eq2……用来代表看上去能够用高等数学介绍过的手段求解出来的微分方 程，如果实在不太熟悉或者忘了，没关系，后面会提供出现错误时候的处理措施；v 代表自变量，默认值为时间 t。边界条件/初始值用 cond1,cond2,...来表示。 • •• D Dy D2y D 2. 在 eq1,eq2……中，用 DD来代表变量的微分，例如 用 表示，而 用 来