matlab实验七__常微分方程.docVIP

实验七 常微分方程 【实验目的】 1.了解常微分方程的基本概念。 2.了解常微分方程的解析解。 3.了解常微分方程的数值解。 4. 学习、掌握MATLAB软件有关的命令。 【实验内容】 如下图所示，一根长l的无弹性细线，一端固定，另一端悬挂一个质量为m的小球。在重力的作用下小球处于平衡位置。若使小球偏离平衡位置一定角度θ，放开它，它就会沿圆弧摆动。在不考虑空气阻力的情况下，小球会做一定周期的简谐运动。利用牛顿第二定理得到如下的微分方程 问该微分方程是线性还是非线性的？是否存在解析解？如果不存在解析解，能否求出其近似解？ 【实验准备】 1.微分方程的概念 2.常微分方程的解析解 3.微分方程的数值解法 4.解微分方程的MATLAB命令 MATLAB中主要用dsolve求符号解析解，ode45,ode23,ode15s求数值解。 matlab求解微分方程命令dsolve，调用格式为： dsolve(‘微分方程’)给出微分方程的解析解，表示为t的函数. dsolve(‘微分方程’，‘初始条件’)给出微分方程初值问题的解，表示为t的函数. dsolve(‘微分方程’，‘变量x’)给出微分方程的解析解，表示为x的函数. dsolve(‘微分方程’，‘初始条件’，‘变量x’)给出微分方程初值问题的解，表示为x的函数. 【实验重点】 1.常微分方程的解析解与数值解的计算 【实验难点】 1.常微