北科大_matlab第六次作业_常微分方程.doc

PAGE PAGE 4 《数学实验》报告 实验名称? 常微分方程 学 院 计算机与通信工程学院 专业班级??? ?? 姓?? 名???? ? 学?? 号??? ? ?? ?? 2015年 6月 【实验目的】 1、学习常微分方程符号解法和数值解法； 2、学习频数直方图的作法。 【实验任务】 1.P168 T24 求解微分方程。 2.P168 T27 用数值方法求解下列微分方程，用不同颜色和线型将y和y’画在同一个图形窗口里： 初始时间：；终止时间：；初始条件： ，。 3.P190 T15 描绘以下数组的频数直方图： 6.8,29.6,33.6,35.7,36.9,45.2,54.8,65.8,43.4,53.8,63.7,69.9,70.7,79.5,97.9,139.4,157.0 【实验程序】 1.P168 T24 y=dsolve(Dy=x*sin(x)/cos(y),x) 2.P168 T27 先将高阶微分方程转化为一阶微分方程，令，即原微分方程化为： 写成矩阵形式为： 放入函数exf.m中，命令如下： [t,x]=ode23(exf,[t0,tf],x0t) 将导数表达式的右端写成exf.m函数文件： function [xdo