Maltab 数学机实验报告 2.docVIP

Maltab 数学上机实验报告 Part 1 实验问题 1蟋蟀鸣叫次数与华氏温度之间的关系 2（1）鱼的体重与身长的三次方之间的拟合关系（2）鱼的体重和身长与胸围二次方之积间的拟合关系 Part 2 问题的分析（涉及的理论知识、数学建模与求解的方法等） Part3 程序设计（粘贴过来） 1蟋蟀鸣叫次数与华氏温度之间的关系 x=[46 49 51 52 54 56 57 58 59 60 61 62 63 64 66 67 68 71 72 71]; y=[40 50 55 63 72 70 77 73 90 93 96 88 99 110 113 120 127 137 132 137]; h1=subplot(1,2,1); plot(x,y,k.,markersize,25); axis([45 73 40 140]); p1=polyfit(x,y,1); p2=polyfit(x,y,2); p3=polyfit(x,y,3); t=45:1:73; s1=polyval(p1,t); s2=polyval(p2,t); s3=polyval(p3,t); hold on plot(t,s1,r-,linewidth,2) plot(t,s2,m-,linewidth,2) plot(t,s3,k-,linewidth,2) title(华氏温度与蟋蟀鸣叫次数拟合模型) xlabel(华氏温度X) ylabel(蟋蟀鸣叫次数Y) grid h2=subplot(1,2,2); plot(x,y,k.,markersize,25); axis([45 73 40 140]); p4=polyfit(x,y,4); p5=polyfit(x,y,5); p6=polyfit(x,y,6); t=45:1:73; s1=polyval(p4,t); s2=polyval(p5,t); s3=polyval(p6,t); hold on plot(t,s4,r-,linewidth,2) plot(t,s5,m-,linewidth,2) plot(t,s6,k-,linewidth,2) title(华氏温度与蟋蟀鸣叫次数拟合模型) xlabel(华氏温度X) ylabel(蟋蟀鸣叫次数Y) grid 2（1）鱼的体重与身长的三次方之间的拟合关系 a=[14.5 12.5 17.25 14.5 12.625 17.75 14.125 12.625]; b=[27 17 41 26 17 49 23 16]; plot(a,b,k.,markersize,25); axis([10 20 15 50]); n=8; a1=a.^6; b1=b.*a.^3; c1=b.^2; a2=sum(a1(1:n)); b2=sum(b1(1:n)); c3=sum(c1(1:n)); k=b2/a2; q=a2*k^2-2*b2*k+c3; x=10:20; y=k*(x.^3); hold on plot(x,y,k-,linewidth,2) xlabel(鱼的长度l); ylabel(鱼的重量w); title(鱼的重量与长度w=kl^3拟合模型) Grid Q1 Q1= 12.1683 2（2）鱼的体重和身长与胸围二次方之积间的拟合关系 l=[14.5 12.5 17.25 14.5 12.625 17.75 14.125 12.625]; g=[9.75 8.375 11.0 9.75 8.5 12.5 9.0 8.5]; a=l.*g.^2; b=[27 17 41 26 17 49 23 16]; n=8; a1=a.^2; b1=b.*a; c1=b.^2; a2=sum(a1(1:n)); b2=sum(b1(1:n)); c3=sum(c1(1:n)); k=b2/a2; q=a2*k^2-2*b2*k+c3; plot3(l,g,b,r*); hold on x=10:20; y=8.5:12.5; [X,Y]=meshgrid(x,y); Z=k*X.*Y.^2; surf(X,Y,Z); shading flat; titl