用MATLAB求解常微分方程实验报告.docxVIP

实验五 利用Matlab求解常微分方程(组)的实验报告 学院：数计学院 班级：1003班 姓名：黄晓丹 学号：1051020144 一.实验目的: 熟悉Matlab软件中关于求解常微分方程的各种命令. 掌握利用Matlab软件进行常微分方程的求解。 二.相关知识 在MATLAB中，由函数dsolve()解决常微分方程（组）的求解问题，其具体格式如下： X=dsolve(‘eqn1’,’eqn2’,…) 函数dsolve用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解. 用字符串方程表示，自变量缺省值为t。导数用D表示，2阶导数用D2表示，以此类推。S返回解析解。在方程组情形，s为一个符号结构。 [tout,yout]=ode45(‘yprime’,[t0,tf],y0) 采用变步长四阶Runge-Kutta法和五阶Runge-Kutta-Felhberg法求数值解，yprime是用以表示f(t,y)的M文件名，t0表示自变量的初始值，tf表示自变量的终值，y0表示初始向量值。输出向量tout表示节点(t0,t1, …,tn)T,输出矩阵yout表示数值解，每一列对应y的一个分量。若无输出参数，则自动作出图形. 三.实验内容: 例1? 求下列微分方程的解析解 （1） （2） （3） 方程（1）求解的MATLAB代