用MATLAB解常微分方程.docxVIP

用MATLAB解常微分方程 实验四求微分方程的解 一、问题背景与实验目的 用通数学建模所而获得的方程，大多半都是微分方程，真实能获得代数方程的时机很少．另一方面，能求解的微分方程也是十分有 限的，特是高方程和偏微分方程（）．就要求我必研究微分方程（）的解法，既要研究微分方程（）的解析解法（精准解），更要研究微分方程（）的数解法（近似解）． 微分方程（）的解析解法(精准解)，Matlab有的函数能够用，本将作一定的介． 本将主要研究微分方程( )的数解法（近似解），重点介Euler折法． 二、有关函数（命令）及简介 1．dsolve(equ1,equ2,?)：Matlab求微分方程的解析解．equ1、equ2、? 方程（或条件）．写方程（或条件）用Dy表示y对于自量的一数，用用D2y表示y对于自量的二数，依此推． 2．simplify(s)：表达式s使用maple的化行化． 比如： symsx simplify(sin(x)^2+cos(x)^2) ans=1 3．[r,how]=simple(s)：由于Matlab提供了多种化，simple命令就是表达式s用各样行化，然后用r返回最形式，how返回形成种形式所用的． 比如： symsx [r,how]=simple(cos(x)^2-sin(x)^2) r=cos(2*x) how=combine 4．[T,Y]=solver(o