Matlab数学实验报告05.docVIP

《数学实验》报告 实验名称?? Matlab高等数学相关运算 ? 学 院 机械工程学院 专业班级??????工业101 ???? 姓??? 名?????于洋??? 学??? 号??? ?? 2011年10月26日 【实验目的】 掌握Matlab高等数学相关的运算； 初步运用概率统计相关的功能。 【实验任务】 1、在钢线碳含量对于电阻效应的研究中，得到以下数据。分别用一次，三次，五次多项式拟合曲线来拟合这组数据并画出图形。 碳含量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 2、已知在某实验中测得某质点的位移s和速度v随时间t变化如下： t 0 0.5 1.0 1.5 2.0 2.5 3.0 v 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411 s 1 1.5 2 2.5 3 3.5 4 3、由方程所确定的隐函数的导数。 4、求下列定积分： I II III 5、用三种方法求下列积分的数值。 I II III 6、绘制下列数组的频数直方图： 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 7、若样本为85,86,78,90,96,82,80,74，求样本均值、标准差、中位数、极差和方差。 【实验程序】 1、 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); disp(1),f1=poly2str(p1,x)； x1=[0.10,0.30,0.40,0.55,0.70,0.80,0.95]; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); plot(x1,y1,--,x1,y3,k-,x1,y5); legend(一次拟合,三次拟合,五次拟合) 2、 t=[0,0.5,1.0,1.5,2.0,2.5,3.0]; v=[0,0.4794,0.8415,0.9975,0.9093,0.5985,0.1411]; s=[1,1.5,2,2.5,3,3.5,4]; pvt=polyfit(t,v,1); pst=polyfit(t,s,1); psv=polyfit(v,s,1); t1=[0,0.5,1.0,1.5,2.0,2.5,3.0]; v1=[0,0.4794,0.8415,0.9975,0.9093,0.5985,0.1411]; y1=polyval(pvt,t1); y2=polyval(pst,t1); y3=polyval(psv,v1); polt(t1,y1,--,t1,y2,k-,v1,y3)； plot(t1,y1,--,t1,y2,k-,v1,y3); legend(v-t,s-t,s=v) 3、 syms x y F=exp(y)+x*y-exp(1); dF_dx=diff(F,x); dF_dy=diff(F,y); pretty(-dF_dx/dF_dy) 4、 I syms x int(x/(1+sqrt(1+x)),1/2,2) II syms x int((sin(x))^7,0,pi/2) III syms x int((1-x^2)^6,1,0) 5、 I x=0:0.01:3; y=exp(-0.5*x)*sin(x+pi/6)； y=exp(-0.5*x).*sin(x+pi/6); t=length(x); s1=sum(y(1:(t-1)))*0.01; s2=sum(y(2:t))*0.01; s3=trapz(x,y); s4=quad(jifen,0,3) II clear all x=0:0.01:2.5; y=exp(-x); t=length(x); s1=sum(y(1:(t-1)))*0.01; s2=sum(y(2:t))*0.01; s3=trapz(x,y); s4=quad(jifen,1,2.5)； III x=0:0.01:2; y=exp(3*x).*sin(2*x); t=length(x); s1