北科大MATLAB第六次作业.doc

《数学实验》报告

实验名称??MATLAB绘图

学院东凌经济管理学院

专业班级?????工管111

姓??名?????李伊童

学??号???

??

2013年5月

【实验目的】

运用MATLAVB求解常微分方程

【实验任务】

第七章、24、27；第八章、15、16

【实验程序】

24、y=dsolve(Dy=x*sin(x)/cos(y),x)

27、编写函数文件

functionxdot=exf(t,x)

u=1-2*t;

xdot=[01;1-1]*x+[01]*u;

编写另一个主文件

t0=0;

tf=pi;

x0t=[0.1;0.2];

[t,x]=ode23(exf,[t0,tf],x0t);

y1=x(:,1);

y2=x(:,2);

plot(t,y1,r-,t,y2,b-)

legend(y,y一阶导数)

15、A=[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];

hist(A,10)

16、A=[85,86,78,90,96,82,80,74];

M=[mean(A),median(A),range(A),var(A),std(A)]

M=

83.875083.500022.000048.69646.9783

【实验结果】

24、y=

-asin(-sin(x)+x*cos(x)-C1)

27、令：x1=y,x2=y’=x1’,y’’=x2’

X2’+tx2-x1=1-2t

X1’=x2

X2’=-tx2+x1+1-2t

即：x1’=0x1+x2+0

X2’=x1-tx2+(1-2t)

U=1-2t

15、

16、

【实验总结】

MATLAB中由于解常微分方程比较中综合，所以要充分理解各种命令，仔细分析方程，再写程序。