数学实验 常微分方程.docVIP

实验六 常微分方程的Matlab解法 一、实验目的 了解常微分方程的解析解。 了解常微分方程的数值解。 学习掌握MATLAB软件有关的命令。 二、实验内容 一根长的无弹性细线，一段固定，另一端悬挂一个质量为的小球，在重力的作用下小球处于垂直的平衡位置。若使小球偏离平衡位置一个角度，让它自由，它就会沿圆弧摆动。在不考虑空气阻力的情况下，小球会做一定周期的简谐运动。利用牛顿第二定律得到如 下的微分方程 问该微分方程是线性的还是非线性的？是否存在解析解？如果不存在解析解，能否求出其近似解？ 三、实验准备 MATLAB中主要用dsolve求符号解析解，ode45,ode23,ode15s求数值解。 s=dsolve(‘方程1’, ‘方程2’,…,’初始条件1’，’初始条件2’ …,’自变量’) 用字符串方程表示，自变量缺省值为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表示数