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

《北科大_matlab第六次作业_常微分方程.docx》由会员分享，可在线阅读，更多相关《北科大_matlab第六次作业_常微分方程.docx（4页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、北科大_matlab第六次作业_常微分方程当前位置：文档视界北科大_matlab第六次作业_常微分方程北科大_matlab第六次作业_常微分方程一、【实验目的】1、学习常微分方程符号解法和数值解法；2、学习频数直方图的作法。二、【实验任务】1.P168T24求解微分方程yxxycossin=。2.P168T27用数值方法求解下列微分方程，用不同颜色和线型将y和y画在同一个图形窗口里：tytyy21-=-+初始时间：00=t；终止时间：=ft；初始条件：1.0|0=ty，2.0|0=ty。3.P190T15描绘下面数组的频数直方图：6.8,29.6,33.6,35.7,36.9,45.2,54.

2、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先将高阶微分方程转化为一阶微分方程，令,2121xyxyxyx=?=，即原微分方程化为：?-+-=txtxxxx2112221写成矩阵形式为：)21(10t-110)21(10t-1102121ttxxxxx-?+?=-?+?=?=xutxdottu?+?-=-=10110,21x放入函数exf.m中，命令如下：t,x=ode23(exf,t0,tf,x0t)(1)将导数表达式的

3、右端写成exf.m函数文件：functionxdot=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.0hist(A,12)四、【实验结果】1.P168T24y=asin(C3+sin(x)-x*cos(x)当前位置：文档视界北科大_matlab第六次作业_常微分方程北科大_matlab第六次作业_常微分方程当前位置：文档视界北科大_matlab第六次作业_常微分方程北科大_matlab第六次作业_常微分方程