北科大matlab数学实验作业第四次.doc.docVIP

《数学实验》报告 实验名称?? 常微分方程求解、曲线拟合及概率统计初步知识 学 院 材料科学与工程 专业班级?????????材料科学与工程 姓??? 名???????? ??? 学??? 号??????? ?? ?? 2011年10月 【实验目的】 掌握常微分方程求解和曲线拟合的方法，能够应用Matlab解决概率中的一些基本运算。 【实验任务】 第130页第8、10题，第168页第24、27题。 【实验程序】 第130页第8题: clf; x=[0.10,0.30,0.40,0.55,0.70,0.80,0.95]; y=[15,18,19,21,22.6,23.8,26]; p1=polyfit(x,y,1), p3=polyfit(x,y,3); p5=polyfit(x,y,5); poly2str(p1,x) poly2str(p3,x) poly2str(p5,x) x1=0.1:0.01:0.95; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); plot(x,y,rp,x1,y1,--,x1,y3,k-.,x1,y5,k); legend(拟合点,一次拟合,三次拟合,五次拟合) 第130页第10题: clf; x=[10,15,20,25,30]; y=[25.2,29.8,31.2,31.7,29.4]; p2=polyfit(x,y,2); p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); poly2str(p2,x) poly2str(p3,x) poly2str(p5,x) poly2str(p7,x) x1=10:1:30; y2=polyval(p2,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,rp,x1,y2,--,x1,y3,k-.,x1,y5,m--,x1,y7,k); legend(拟合点,二次拟合,三次拟合,五次拟合,七次拟合) function y=concentration(x) y= -1.0578e-005* x^5 + 0 x^4 - 0.018322 *x^3 + 0.067889 *x^2+ 3.0013 *x; end concentration(18) ans = 30.8259 concentration(26) ans = 31.6449 第168页第24题: y=dsolve(Dy=x*sin(x)/cos(y),x) 第168页第27题: 先将高阶微分方程化为一阶微分方程。令， 即原微分方程化为： 写成矩阵形式为： 放入函数exf.m中，命令如下： [t,x]=ode23(exf,[t0,tf],x0t) 其中，t0=0,tf=,x0t=。 function xdot=exf(t,x) u=1-2.*t; xdot=[0 1;1 -t]*x+[0 1]*u; end clf; t0=0;tf=pi;x0t=[0;0]; [t,x]=ode23(exf,[t0,tf],x0t) y=x(:,1), Dy=x(:,2); clf, plot(t,y,m-,t,Dy,b*) legend(y,Dy ) 【实验结果】 第130页第8题: 第130页第10题: 第168页第24题: 第168页第27题: 【实验总结】 这节课我们学习了用MATLAP作常微分方程的求解、多项式拟合等内容。MATLAP强大的运算功能大大方便了我们对数据的统计，通过形成各种图形使使用者得到直观的答案。 2