土木1103水成霖41112061第6次作业.docVIP

《数学实验》报告 实验名称常微分方程求解、曲线拟合及概率统计初步知识 学 院 土木与环境工程学院 专业班级?????土木1103 姓??? 名????水成霖 ??? ? 学??? 号??? ???? ?? 2013年 6月 【实验目的】 掌握常微分方程求解和曲线拟合的方法，能够应用Matlab解决概率中的一些基本运算． 【实验任务】 p168 24 27 P190 15 16 【实验程序】 P168,24 y=dsolve(Dy=x*sin(x)/cos(y),x) 27 令，x1=y,x2=y′=x1′,y″=x2′ x2′+t*x2-x1=1-2*t x1′=x2即x1′=0*x1+x2+0 x2′=-t*x2+x1+1-2*t 即x2′=x1-t*x2+(1-2*t) 写成矩阵形式 X′=[0 1 ]x+[0]u 1 -t 1 U=[1-t^2/(2*pi)] function xdot=exf(t,x) u=1-2*t; xdot=[0,1;1,-1]*x+[0 1]*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,bo) legend(y,y一阶导数) P190第15题 load a.txt a a = 6.8000 29.6000 33.6000 35.7000 36.9000 45.2000 54.8000 65.8000 43.4000 53.8000 63.7000 69.9000 70.7000 79.5000 97.9000 139.4000 157.0000 hist(a,10) P190第16题 a a = 85 86 78 90 96 82 80 74 mean(a) ans = 83.8750 std(a) ans = 6.9783 median(a) ans = 83.5000 range(a) ans = 22 var(a) ans = 48.6964 【实验结果】 y = asin(C3 + sin(x) - x*cos(x)) P190第16题 ans = 83.8750 ans = 6.9783 ans = 83.5000 ans = 22 ans = 48.6964 五、【实验总结】 对matlab中的有关常微分的解法，尤其是数值法解常微分的方法要重点掌握，对二次微分方程化为一次方程的过程还不是很熟悉，也不习惯运用矩阵的方式解方程，在以后的学习中要加以熟练。对数据的处理掌握较好，尤其是求特征数据的命令。只是有时候因为路径或者其他问题，载入不了桌面上的数据，以后应该细心应对。 6