matlab数学实验结课作业,北科大matlab数学实验_第六次作业.doc

北科大matlab数学实验_第六次作业.doc

《数学实验》报告实验名称常微分方程学院计算机与通信工程学院专业班级姓名学号2015年6月2一、【实验目的】1、学习常微分方程符号解法和数值解法；2、学习频数直方图的作法。二、【实验任务】1.P168T24求解微分方程。yxcosin 2.P168T27用数值方法求解下列微分方程，用不同颜色和线型将y和y’画在同一个图形窗口里：tt21 初始时间：；终止时间：；初始条件：，。0tft.0|ty2.0| ty3.P190T15描绘以下数组的频数直方图：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.P168T24y=dsolve( Dy=x*sin(x)/cos(y) , x )2.P168T27先将高阶微分方程转化为一阶微分方程，令，即原微分方程化为： ,2121xyxytt 12写成矩阵形式为：)21(0t-1)(0t-10 212ttxxx3放入函数exf.m中，命令如下：utxdotu10,21x[t,x]=ode23( exf ,[t0,tf],x0t)(1)将导数表达式的右端写成exf.m函数文件：function[xdot]=exf(t,x)u=1-2*t;xdot=[01;1-t]*x+[01] *u;end(2)主程序如下：clft0=0,tf=pi,x0t=[0.1,0.2];[t,x]=ode23( exf ,[t0,tf],x0t);y=x(:,1);y1=x(:,2);plot(t,y, r- ,t,y1, b-- );legend( y , y的一阶导数 );3.P190T15A=[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]hist(A,12)四、【实验结果】1.P168T24y=asin(C3+sin(x)-x*cos(x))42.P168T273.P190T155五、【实验总结】1、学会了常微分方程符号解法和数值解法；2、学会了频数直方图的作法。