大学数学建模实习报告作业汇总

一、线性规划问题评分：题目1： 一车队有8辆车，这8辆车存放在不同地点，队长要派其中5辆到个不同的工地去运货。各车从存放处调到装货地点所需费用列于表2.5.问应选用哪五辆车调到何处去运货，才能使各车从所在地点调到装货地点所需的费用最少（要求分别用matlab和Lingo编程求解）？ 表2.5 调车费用表-装货地点车号123456781302518322719222622931191821203019328293019192223264293019242519182152120181716141618（1）解：记表示第j号车调到装货地点i所需的费用。引进0-1变量，建立如下的0-1整数规划模型：（1）编写的matlab程序如下clc,clearc=30 25 18 32 27 19 22 26 29 31 19 18 21 20 30 19 28 29 30 19 19 22 23 26 29 30 19 24 25 19 18 21 21 20 18 17 16 14 16 18;c=c(:);a=zeros(8,40);for j=1:8 a(j,(j-1)*5+1:5*j)=1; endb=ones(8,1);d=zeros(5,40);for i=1:5 d(i,i:5:40)=1;end e=ones(5,1);x,fval=bintprog(c,a,b,d,e);x=reshape(x,5,8) fval运行结果如下： x = (1,3) 1 (2,4) 1 (3,5) 1 (5,6) 1 (4,7) 1 fval = 87所以最优方案为Lingo程序如下model:sets:col/1.5/:i;var/1.8/:j;links(col,var):c,x;endsetsdata:c=30 25 18 32 27 19 22 26 29 31 19 18 21 20 30 19 28 29 30 19 19 22 23 26 29 30 19 24 25 19 18 21 21 20 18 17 16 14 16 18;enddatamin=sum(links:c*x);for(col(i):sum(var(j):x(i,j)=1);for(var(j):sum(col(i):x(i,j)=1);for(links:bin(x);end运行结果为Objective value: 87.00000 X( 1, 1) 0.000000 30.00000 X( 1, 2) 0.000000 25.00000 X( 1, 3) 1.000000 18.00000 X( 1, 4) 0.000000 32.00000 X( 1, 5) 0.000000 27.00000 X( 1, 6) 0.000000 19.00000 X( 1, 7) 0.000000 22.00000 X( 1, 8) 0.000000 26.00000 X( 2, 1) 0.000000 29.00000 X( 2, 2) 0.000000 31.00000 X( 2, 3) 0.000000 19.00000 X( 2, 4) 1.000000 18.00000 X( 2, 5) 0.000000 21.00000 X( 2, 6) 0.000000 20.00000 X( 2, 7) 0.000000 30.00000 X( 2, 8) 0.000000 19.00000 X( 3, 1) 0.000000 28.00000 X( 3, 2) 0.000000 29.00000 X( 3, 3) 0.000000 30.00000 X( 3, 4) 0.000000 19.00000 X( 3, 5) 1.000000 19.00000 X( 3, 6) 0.000000 22.00000 X( 3, 7) 0.000000 23.00000 X( 3, 8) 0.000000 26.00000 X( 4, 1) 0.000000 29.00000 X( 4, 2) 0.000000 30.00000 X( 4, 3) 0.000000 19.00000 X( 4, 4) 0.000000 24.00000 X( 4, 5) 0.000000 25.00000 X( 4, 6) 0.000000 19.00000 X( 4, 7) 1.000000 18.00000 X( 4, 8) 0.000000 21.00000 X( 5, 1) 0.000000 21.00000 X( 5, 2) 0.000000 20.00000 X( 5, 3) 0.000000 18.00000 X( 5, 4) 0.000000 17.00000 X( 5, 5) 0.000000 16.00000 X( 5, 6) 1.000000 14.00000 X( 5, 7) 0.000000 16.00000 X( 5, 8) 0.000000 18.00000所以最优方案为第2题某单位需要加工制作100套钢架，每套用长为2.9m，2.1m和1m的圆钢各一根。已知原料长6.9m,问应如何下料，使用的原材料最省。解：列出所有可行性方案如下表方案长度一X1二X2三X3四X4五X5六X6七X72.92.11 11114 21 32 2146合计（米）剩余（米）60.96.906.80.16.3066.20.76.10.860.9第一种方案x1：切割2.9米一截,2.1米一截,1米一截，浪费0.9米；第二种方案x2:切割2.9米一截，1米4截，浪费0米；第三中方案x3：切割2.9米2截,1米1截，浪费0.1米；第四种方案x4：切割2.1米3截，浪费0.6米；第五种方案x5:切割2.1米2截，1米2截，浪费0.7米；第六种方案x6：切割2.1米1截，1米4截，浪费0.8米；第七种方案x7：切割1米6截，浪费0.9;设：使用第i方案xi次（i=,3,4,5,6,7）Min=x1+x2+x3+x4+x5+x6+x7；model:sets:row/1.3/:b;col/1.7/:c,x;links(row,col):a;endsetsdata:c=1 1 1 1 1 1 1;a=-1 -1 -2 0 0 0 0-1 0 0 -3 -2 -1 0-1 -4 -1 0 -2 -4 -6;b=-100 -100 -100;enddatamin=sum(col:c*x);for(row(i):sum(col(j):a(i,j)*x(j)=b(i);for(col(j):gin(x);end运行结果：Objective value: 91.00000 X( 1) 0.000000 1.000000 X( 2) 12.00000 1.000000 X( 3) 44.00000 1.000000 X( 4) 33.00000 1.000000 X( 5) 1.000000 1.000000 X( 6) 0.000000 1.000000 X( 7) 1.000000 1.000000则最少根数为91根，选择的方案为第二种方案12根，第三种方案44根，第四种方案33根，第五种方案1根，第七种方案1根。二、图与网络问题评分：第一题 已知有6个村庄，各村庄的小学生人数如表4.11所列，各村庄间的距离如表4.5所示。现在计划建造一所医院和一所小学，问医院应建在哪个村庄才能使最远村庄的人到医院看病所走的路最短？又问小学建在哪个村庄使得所有小学生走的总路程最短？ 表4.11 各村小学生人数村庄小学生504060207090解：（1）建立赋权完全图G=（V，E，W），其中V=（），权重邻接矩阵为： 利用Floyd算法可以求出任意两个村庄的最短距离矩阵：其中表示第i村庄到第j村庄的最短距离。A的第j列表示其他各村庄到该村的最短距离，第j列最大值最远的村到该村的最短距离，由矩阵su可知，每一列的最大值为11,9,6,7,8,11；所以医院应建在第三村庄使得最远村庄的人到医院看病所走的路最短，为6.相应的matlab程序如下：clear;clc;a=zeros(6);a(1,2)=2;a(1,3)=7;a(2,3)=4;a(2,4)=6;a(2,5)=8;a(3,4)=1;a(3,5)=3;a(4,5)=1;a(4,6)=6;a(5,6)=3;a=a;a=sparse(a);A=graphallshortestpaths(a,Directed,0)su=max(A)（2）A矩阵的第i行表示其他村庄到第i村庄的最短距离，设表示第i村庄的人数（i=1,2.6）,若在第i村上学，则所有小学生走的总路程为，则当时，即为所求。用matlab求解得：s= 0 2 6 7 8 11; 2 0 4 5 6 9; 6 4 0 1 2 5; 7 5 1 0 1 4; 8 6 2 1 0 3; 11 9 5 4 3 0; a=50 40 60 20 70 90; s=A*a结果如下列表：小学生走各村总路程S1S2S3S4S5S6213016701070104010501500从表中可以看出小学建在第4村总路程最短，为1040.第二题乙两个煤矿分别生产煤500万吨，供应A、B、C三个电厂发电需要，各电厂用量分别为300、300、400（单位:万吨）。已知煤矿之间、煤矿与电厂之间以及各电厂之间相互距离（单位：km）如表4.13表4.15所列。煤可以直接运达，也可以间接运达试确定从煤矿到各电厂间煤的最优调配方案。 表4.13 各煤矿之间的距离甲乙甲0120乙1000 表4.14 从两煤矿到三个发电厂之间的距离ABC甲15012080乙6016040 表4.15 三个电厂之间的距离ABCA070100B500120C1001500如图所示把甲，乙，A，B，C分别用1，2，3，4，5表示。由题意知这是一个多源多汇图。我们把源点聚为s点，收点聚为t点。这题就形成了求从s点到t点的最大流最小费用问题。s12345t1150可得到如下数学模型：model:sets:nodes/s,1,2,3,4,5,t/:d;arcs(nodes,nodes)/s 1,s 2,1 2,1 3,1 4,1 5,2 1,2 3,2 4,2 5,3 4,3 5,4 3,4 5,5 3,5 4,3 t,4 t,5 t/:b,c,f;endsetsdata:d=1000 0 0 0 0 0 -1000; b=0 0 120 150 120 80 100 60 160 40 70 100 50 120 100 150 0 0 0;c=500 500 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 300 300 400;enddatamin=sum(arcs:b*f);for(nodes(i):sum(arcs(i,j):f(i,j)-sum(arcs(j,i):f(j,i)=d(i);for(arcs:bnd(0,f,c);end求解结果：F（1,4）=300，F（1,5）=200，F（2,3）=300，F(2,5)=200，甲给B，C各运300,200万吨，乙给A,C各运300，200万吨。1000万吨媒运输的总路程为78000（Km*万吨）。三、微分方程问题评分：第一题 根据经验，当一种新商品投入市场后，随着人们对它的拥有量的增加，其销售量s(t)下降的速度与s(t)成正比。广告宣传可给销量添加一个增长速度，它与广告费a(t)成正比，但广告只能影响这种商品在市场上尚未饱和的部分（设饱和量为M）。建立一个销售s(t)的模型。若广告宣传只进行有限时间T，且广告费为常数a，问s(t)如何变化？解：设与s（t）是成正比的。所以在没有广告的影响下，又因为广告费可以给销量增添一个增长速度，与a（t）成正比，且只影响未饱和部分。所以建立如下模型：，（1）（1）当s(t)=M或a(t)=0时：（2）假设采用如下的广告策略： 当时，将a(t)=a/r带入（1）式得：可推出令=b, =c.原式可化为如下的常微分方程：若令s（t）=，求出解为：用matlab求y=dsolve(Ds+b*s=c,s(0)=s0,t)解得：y =(c - (c - b*s0)/exp(b*t)/b；简化为：第二题有高为4m的半椭球形容器，水从它的底部小孔流出。小孔横截面积为1cm2。开始时容器内盛满了水，求水从小孔流出过程中容器里水面的高度（水面与孔口中心的距离）随时间 变化的规律。 解：以底部中心为坐标原点，垂直向上为坐标的正向建立坐标系。由能量守恒原理得到：解得。设在微笑时间间隔内，水面高度由高度h降到h+dh（这里dh为负值）， 由物质守恒原理得：是 （1），其中：为底部小孔的面积。有：（2）把，和（2）代入（1）式，化简得：在考虑初始条件，得到如下微分方程模型：解得微分方程的解：t=dsolve(Dt=20000*pi*(4-h)2/16-1)/sqrt(2*g*h),t(4)=0,h);t=simple(t),pretty(t)解得：四、目标规划问题评分：一最近，某节能灯具厂接到了订购16000套A型和B型节能灯具的订购合同，合同中没有对这两种灯具各自的数量做要求，但合同要求工厂在一周内完成生产任务并交货。根据该厂的生产能力，一周内可以利用的生产时间为20000min，可利用的包装时间为36000min。生产完成和包装完成一套A型节能灯具各需要2min；生产完成和包装完成一套B型节能灯具分别需要1min和3min。每套A型节能灯具成本为7元，销售价格为15元，即利润为8元；每套B型节能灯具成本为14元，售价为20元，即利润为6元。厂长首先要求必须要按照合同完成订货任务，并且既不要有不足量，也不要有超过量。其次要求满意的销售额尽量达到或接近275000元。最后要求在生产总时间和包装总时间上可以增加，但超过量尽量地小。同时注意到增加生产时间要比增加包装时间困难得多。试为该节能灯具厂制定生产计划。 解：根据题意，确定目标和优先级：第一级目标：首先要求必须要按照合同完成订货任务，赋予优先因子为p1;第二级目标：其次要求满意的销售额尽量达到或接近275000元，赋予优先因子为p2；第三级目标：最后要求在生产总时间和包装总时间上可以增加，但超过量尽量地小，赋予优先因子为p3；建立相关目标约束，设生产A型和B型节能灯各x1，x2套。（1） 订单数量约束。设分别表示未达到和超过订单量的偏差量，又因为只要求完成16000套交货数量即可，则数量目标约束为（2） 销售利润目标。因为要求利润尽量完成和接近275000元，因此得（3） 生产总时间和包装总时间的限制。因为要求在生产总时间和包装总时间上可以增加，但超过量尽量地小，同时注意到增加生产时间要比增加包装时间困难得多，则可设二者的加权系数为0.4和0.6，则写出目标规划的数学模型，即写出相应的lingo程序如下：sets:a/1.3/:p,z,goal;b/1.2/:x;c/1.4/:g,dplus,dminus;d(c,b):e;obj(a,c)/1 1,2 2,3 3,3 4/:wplus,wminus;endsetsdata:ctr=?;goal=？0;g=16000 275000 20000 36000;e=1 1 15 20 2 1 2 3 ;wplus=1 0 0.4 0.6;wminus=1 1 0 0;enddatamin=sum(a:p*z);p(ctr)=1;for(a(i)|i#ne#ctr:p(i)=0);for(a(i):z(i)=sum(obj(i,j):wplus(i,j)*dplus(j)+wminus(i,j)*dminus(j);for(c(i):sum(b(j):e(i,j)*x(j)+dminus(i)-dplus(i)=g(i);for(a(i)|i#lt#size(a):bnd(0,z(i),goal);end运行结果如下： Objective value:3800.000 X(1) 9000.000 0.000000 X(2) 7000.000 0.000000 DPLUS(3) 5000.000 0.000000 DPLUS(4) 3000.000 0.000000求得的最优解为生产A型灯具9000套，B型灯具7000套。生产时间需要增加5000min，包装时间需要增加3000min，即可完成任务。五、评价与预测问题评分：一某商品的生产需要甲乙两种原料，产品利润，以及甲乙两种原料的市场供给等数据如表15-7所列，试预测2004年，甲的供应量为400KG，乙的供应量为500KG时的产品利润（要求建立灰色GM（1，1）。表15-7 原始数据表年度19902000200120022003i12345产品利润/元43837625105001131617818甲原料/KG83131180195306乙原料/KG146212233259404解：首先，为了保证建模方法的可行性，需要对已知数据做必要的检验处理。这里有三个变量，分别表示利润，甲原料，乙原料。对应三个原始数据，设参数数据分别为（1） 求级比（i=1,2,3），有 ，可以得=（0.5748,0.7262,0.9279,0.6351） =（0.6336,0.7279,0.9231,0.6373）（1）级比判断。由于（0.7165,1.3307），（i=1,2,3），所以需要对序列做必要的变换处理，使其落入可覆盖内。即取适当的常数c，作平移变换对作平移变换，k=1,2.5，则=（0.7621，0.8258，0.9529，0.7270）（0.7165,1.3307）对作平移变换，+100，k=1,2.5,则=（0.7922，0.8250，0.9492，0.7266）（0.7165,1.3307）对作平移变换，+150，k=1,2.5,=（0.8177，0.9452，0.9364，0.7383）（0.7165,1.3307）由于所有的（0.7165,1.3307），k=2,35,故可以用（i=1,2,3）作令人满意的GM（1,1）建模。2）GM(1,1)建模（1）对平移变换后数据（i=1,2,3）作一次累加，得到=（10383，13625，16500，1