MATLAB数学实验报告2.docx

《MATLAB数学实验报告2.docx》由会员分享，可在线阅读，更多相关《MATLAB数学实验报告2.docx（25页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、MATLAB数学实验报告2MATLAB数学实验报告姓名：李帆班级：机械硕21学号：2120204008第一次数学实验报告线性规划问题一，实验问题1，某饲养场饲养动物出售，设每头动物天天至少需要700g蛋白质，30g矿物质，100mg维生素。现有五种饲料可供选择，各种饲料的每千克营养成分含量和单价如下表。是确定既能知足动物生长的营养需要，游客是费用最省的选用饲料方案。2，某工厂生产甲、乙、丙三种产品，单位产品所需工时分别为2、3、1个；单位产品所需原料分别为3、1、5公斤；单位产品利润分别为2、3、5元。工厂天天可利用的工时为12个，可供给的原料为15公斤。为使总利润为最大，试确定日生产计划和最

2、大利润。二，问题分析1，1该题属于采用线性规划的方式求出最优解的数学问题。该题有下面特点，1.目的函数有线性，是求目的函数的最小值;2.约束条件为线性方程组；3.未知变量都有非负限制。1，2求解该类问题的方法有图解法，理论解法和软件解法。图解法常用于解变量较少的线性规划问题。理论解法要构建完好的理论体系。目前用于解线性规划的理论解法有：单纯形法，椭球算法等。在此，我们采用单纯形法的MATLAB软件解法来求解该问题。1，3此题中，要求既要知足动物生长的营养需要，又要使费用最省，则使每种饲料的选用量为变量，以总费用的最小值为所求量，同时每种饲料的使用量要符合营养成分的要求。1，4在此，首先确定建立

3、线性规划模型。设饲料i选用量为xi公斤，i=1,2,3,4,5.则有模型：Minz=0.2x1+0.7x2+0.4x3+0.3x4+0.8x5s.t.3x1+2x2+6x4+18x5=700;x1+0.5x2+0.2x3+2x4+0.5x5=300.5x1+x2+0.2x3+2x4+0.8x5=100Xj=0,j=1,2,3,4,5解之得：x1=x2=x3=0X4=39.74359X5=25.14603Zmin=32.435902,1该问题与第一题分析步骤类似，故只在此写出其线性规划模型Z=2x+3y+5z2x+3y+z第一题：c=0.2,0.7,0.4,0.3,0.8A=3,2,1,6,18

4、;1,0.5,0.2,2,0.5;0.5,1,0.2,2,0.8;1,0,0,0,0;0,1,0,0,0;0,0,1,0,0;0,0,0,1,0;0,0,0,0,1b=700,30,100,0,0,0,0,0x,fval=linprog(c,-A,-b)c=0.20000.70000.40000.30000.8000A=3.00002.00001.00006.000018.00001.00000.50000.20002.00000.50000.50001.00000.20002.00000.80001.0000000001.0000000001.0000000001.0000000001.00

5、00b=7003010000000Optimizationterminated.x=0.0000-0.00000.000039.743625.6410fval=32.4359第二题c=-2-3-5A=231;315b=12;15lb=000x,Z,exitflag,output=linprog(c,A,b,lb,)将上述程序输入matlab。就可得出结果。结果如下：x=0.00003.21432.3571Z=-21.4286四，结果分析和结论1.由上述实验结果可知，当x1=x2=x3=0，x4=39.74539公斤，x3=25.64103公斤时，费用获得最小值32.43590元。2.由上述结果

6、能够知道当甲、乙、丙分别生产0,3.2143,2.3571公斤时，总利润最大，可到达21.4286四，总结和体会通过该次试验，我们把握了采用线性规划方式解决求最优解的数学问题的方法，同时，对于建立数学模型有了一个初步的认识。对于以后的数学学习起到了非常大的帮助。并且，我们把握了使用MATLAB软件解决问题的一些方法，学到了关于其操作的一些软件语言，开拓了解决数学相关问题的一种新方法。这对于我们来讲，是一次极大的收获。当然，由于我们本身知识的有限，对于解决一些较为复杂的问题还有一些难度，同时，对于其中的一些符号语言以及一些表示方法还不够熟悉，需要我们以后进行进一步的学习。第二次数学实验报告曲线拟

7、合一，实验问题1，下表中，X是华氏温度，Y是一分钟内一只蟋蟀的鸣叫次数，试用多项式模型拟合这些数据，画出拟合曲线，分析你的拟合模型能否很好2、(1)在下列数据中，W表示一条鱼的重量，l表示它的长度，使用最小二乘准则拟合模型W=kl3长度l英寸14.512.517.2514.512.62517.7514.12512.625重量w盎司2717412617492316(2)*在下列数据中，g表示一条鱼的身围，使用最小二乘准则拟合模型W=klg2长度l英寸14.512.517.2514.512.62517.7514.12512.625身围g英寸9.758.37511.09.758.512.59.08.

8、5重量w盎司2717412617492316(3)*两个模型哪个拟合数据较好？为何？二，问题分析1,1该题属于曲线拟合问题。最佳曲线拟合问题等价于确定一条平面曲线使它和实验数据点最接近。并不要求曲线严格通过每个已知数据点，但要求曲线在各数据点处的取值与已知观测值的总体误差最小。1,2首先设定某一类型的函数y=f(x),，然后确定函数中的参数，使得在各点处的偏差ri=f(xi)-yi的平方和最小。1,3此题中采用最小二乘法的基本思想。运用MATLAB中命令polyfit编写拟合程序。其中x,y为已知的数据点横坐标向量和总坐标向量，m为要拟合的多项式次数，结果p为m次多项式系数向量，且从高次到低次

9、存放在向量p中。2.1该题属于曲线拟合问题。最佳曲线拟合问题等价于确定一条平面曲线使它和实验数据点最接近。并不要求曲线严格通过每个已知数据点，但要求曲线在各数据点处的取值与已知观测值的总体误差最小。2.2首先设定某一类型的函数y=f(x),，然后确定函数中的参数，使得在各点处的偏差ri=f(xi)-yi的平方和最小。三，程序设计流程图第一题x=4649515254565758596061626364666768717271;y=40505563727077709093968899110113120217137132137;plot(x,y,k.,markersize,20);axis(4080

10、35145);p=polyfit(x,y,8);t=40:1:80;s=polyval(p,t);holdonplot(t,s,k-,linewidth,2)grid第二题1l=14.5;12.5;17.25;14.5;12.625;17.75;14.125;12.625;w=27;17;41;26;17;49;23;16;c=0;d1=0;fori=1:8c=c+l(i)4;d1=d1+l(i)*w(i);endA=cD=d1ab=inv(A)*Dfori=1:8x(i)=ab(1)*l(i)3;endplot(l,w,r*-,l,x,b.-);2l=14.512.517.2514.512.

11、62517.7514.12512.625g=9.758.37511.09.758.512.59.08.5w=2717412617492316plot3(l,g,w,k.,markersize,25)axis(10207.5121555)m=l.*(g.2)k=inv(m*(m.)*(m)*(w.)x=11:0.1:20y=8:0.1:13X,Y=meshgrid(x,y)Z=k*X.*Y.2surf(X,Y,Z)shadingflat四，结果分析和结论运行结果如下列图所示：第一题：第二题1l=14.500012.500017.250014.500012.625017.750014.125012

12、.6250g=9.75008.375011.00009.75008.500012.50009.00008.5000w=2717412617492316m=1.0e+003*1.37840.87682.08731.37840.91222.77341.14410.9122k=0.01873经过对于实验题1和实验题2的数据拟合曲线的比拟，我以为第二个数学模型拟合数据较好。由于鱼的身体重不仅遭到其长度的影响，也与他的身围有一定的关系。通过综合数据的拟合，能够愈加准确的判定鱼的体重的影响关系。五，总结和体会通过这次试验，我们对于数据拟合模型有了一定的认识，同时对于其基本原理和使用方法有了一定的把握。我们发现，该数学模型，在我们的日常生活中，具有极大的使用范围，能够为我们带来相当大的便利。