数学建模上机实验报告4

  • 格式:doc
  • 大小:170.00 KB
  • 文档页数:7

西华大学数学建模上机实验报告课程名称: 数学建模 年级:2011级 上机实验成绩: 指导教师:蒲俊老师姓名:李国强上机实验名称:数据拟合与线性规划 学号:362011*********上机实验日期:2013.03.22上机实验编号: 组号:上机实验时间:一、实验目的直观了解插值与拟合的基本内容 掌握用数学软件求解插值、拟合问题 了解线性规划的基本内容掌握用数学软件包求解线性规划问题 二、实验内容1.已知飞机下轮廓线上数据如下 (1)作数据的点图形。

(2)确定X 和 Y 之间的近似关系。

bn an t +=23.用MA TLAB 或 Lingo 求解线性规划问题6543218121110913min x x x x x x z +++++=⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=≥≤++≤++=+=+=+6,,2,1,09003.12.15.08001.14.0500600400x ..654321635241 i x x x x x x x x x x x x t s i4.用MA TLAB 或 Lingo 求解线性规划问题X0 3 5 7 9 11 12 13 14 15 Y1.21.72.02.12.01.81.21.01.6m a x 6543216.064.072.032.028.04.0x x x x x x z +++++=85003.003.003.001.001.001.0..654321≤+++++x x x x x x t s70005.002.041≤+x x 10005.002.052≤+x x 90008.003.063≤+x x6,2,10 =≥j x j5.用MA TLAB 或 Lingo 求解线性规划问题321436min x x x z ++= 120..321=++x x x t s 301≥x5002≤≤x203≥x三、使用环境MATLAB 7.0.1四、核心代码及调试过程1. (1).>> x=[0 3 5 7 9 11 12 13 14 15]; >> Y=[0 1.2 1.72.0 2.1 2.0 1.8 1.2 1.0 1.6]; >> plot(x,y,'.')Undefined function or variable 'y'.Did you mean:>> plot(x,Y,'.')(2)、>> a=polyfit(x,Y,2) a =-0.0249 0.4416 0.0683>> a=polyfit(x,Y,3) a =0.0012 -0.0517 0.5939 -0.0541>> a=polyfit(x,Y,4) a =0.0004 -0.0123 0.0769 0.2146 0.0300 2. >> t=[0 20 40 60 80 100 120 140 160 183.5];>> n=[0 1153 2054 2800 3466 4068 4621 5135 5619 6152]; >> plot(t,n,'.') >> a=polyfit(n,t,2) a =0.0000 0.0144 0.0263 3. 用MATLAB 或 Lingo 求解线性规划问题6543218121110913min x x x x x x z +++++=⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=≥≤++≤++=+=+=+6,,2,1,09003.12.15.08001.14.0500600400x ..654321635241 i x x x x x x x x x x x x t s i4.用MA TLAB 或 Lingo 求解线性规划问题m a x 6543216.064.072.032.028.04.0x x x x x x z +++++=85003.003.003.001.001.001.0..654321≤+++++x x x x x x t s70005.002.041≤+x x 10005.002.052≤+x x 90008.003.063≤+x x6,2,10 =≥j x j5.用MA TLAB 或 Lingo 求解线性规划问题321436min x x x z ++= 120..321=++x x x t s 301≥x5002≤≤x203≥x三、使用环境Windows7MATLAB 7.0.1四、核心代码及调试过程1)……………………………………………………………………………………………………… clear clf clcx0=[0 3 5 7 9 11 12 13 14 15 ];y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 ]; x=0:0.1:15;y=interp1(x0,y0,x);y1=interp1(x0,y0,x,'spline'); y2=interp1(x0,y0,x,'nearest'); y3=interp1(x0,y0,x,'cubic'); y4=interp1(x0,y0,x,'lagrange'); subplot(2,3,1)plot(x0,y0,'*',x,y,'r') title('linear') subplot(2,3,2)plot(x0,y0,'*',x,y1,'g') title('spline') subplot(2,3,3)plot(x0,y0,'*',x,y2,'m') title('nearest') subplot(2,3,4)plot(x0,y0,'*',x,y3,'c') title('cubic') subplot(2,3,5)plot(x0,y0,'*',x,y4,'y') title('lagrange')2)…………………………………………………………………………………………………..clearclfclct=[0 20 40 60 80 100 120 140 160 183.5];y=t';n=[0000 1153 2045 2800 3466 4068 4621 5135 5619 6152];x=[(n.^2)' n' zeros(length(t),1)];plot(n,t,'*k')b=regress(y,x);format short gp=b't1=polyval(p,n)hold onplot(n,t1,'r')legend('t(原始数据)','t1(拟合后)',4)p =2.5048e-06 0.014405 0t1 =Columns 1 through 80 19.939 39.934 59.972 80.019 100.05 120.05 140.02Columns 9 through 10160.03 183.423)…………………………………………………………………………………………….clearclcC=[13 9 10 11 12 8];A=[0.4 1.1 1 0 0 0;0 0 0 0.5 1.2 1.3];b=[800;900];Aeq=[1 0 0 1 0 0;0 1 0 0 1 0;0 0 1 0 0 1];beq=[400;600;500];VLB=[0;0;0;0;0;0];VUB=[];[X,fval]=linprog(C,A,b,Aeq,beq,VLB,VUB)Optimization terminated.X =0.0000600.00000.0000400.00000.0000500.0000fval =1.3800e+004 4)………………………………………………………………………………………………. clearclcC=[-0.4 -0.28 -0.32 -0.72 -0.64 -0.6];A=[0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.05 0;0 0 0.03 0 0 0.08];b=[850;700;100;900];Aeq=[];beq=[];VLB=[0;0;0;0;0;0];VUB=[];[X,fval]=linprog(C,A,b,Aeq,beq,VLB,VUB)Optimization terminated.X =1.0e+004 *3.50000.50003.00000.00000.00000.0000fval =-2.5000e+004 5)…………………………………………………………………………………………clearclcC=[6 3 4];A=[0 1 0];b=[50];Aeq=[1 1 1];beq=[120];VLB=[30;0;20];VUB=[];[X,fval]=linprog(C,A,b,Aeq,beq,VLB,VUB)Optimization terminated.X =30.000050.000040.0000fval =490.0000五、总结基本上能解决拟合问题,遇到问题能够及时解决并掌握,对于线性规划问题也能灵活运用。