《数学实验实验报告.doc》由会员分享,可在线阅读,更多相关《数学实验实验报告.doc(10页珍藏版)》请在课桌文档上搜索。
1、-数学实验实验报告班级:学号:姓名;一、实验问题1指派问题考虑指定n个人完成n项任务每人单独承当一项任务,使所需的总完成时间本钱尽可能短. *指派问题的有关数据每人完成各任务所需的时间如下表所示,试建模并求解该指派问题。2二次指派问题*公司指派n个员工到n个城市工作每个城市单独一人,希望使所花费的总费用尽可能少。n个员工两两之间每个月通话的时间表示在下面的矩阵的上三角局部因为通话的时间矩阵是对称的,没有必要写出下三角局部,n个城市两两之间通话费率表示在下面的矩阵的下三角局部同样道理,因为通话的费率矩阵是对称的,没有必要写出上三角局部. 试求解该二次指派问题。3、金星第四章课后习题第1或3题任选
2、一题。二、问题的分析涉及的理论知识、数学建模与求解的方法等1)根据实际问题,建立数学优化模型 2)根据优化模型,利用LINGO 来求解模型。三、计算过程、结论和结果分析1.模型:model:sets:m/1.4/;n/1.4/;link(m,n):a,*;endsetsmin=sum(link(i,j):*(i,j)*a(i,j);for(m(i):sum(n(j):a(i,j)=1);for(n(j):sum(m(i):a(i,j)=1);data:*=15 18 21 24 19 23 22 18 26 18 16 19 19 21 23 17;enddataend结果:Global op
3、timal solution found. Objective value: 70.00000 Infeasibilities: 0.000000 Total solver iterations: 7 Variable Value Reduced Cost A( 1, 1) 0.000000 0.000000 A( 1, 2) 1.000000 0.000000 A( 1, 3) 0.000000 5.000000 A( 1, 4) 0.000000 10.00000 A( 2, 1) 1.000000 0.000000 A( 2, 2) 0.000000 1.000000 A( 2, 3)
4、0.000000 2.000000 A( 2, 4) 0.000000 0.000000 A( 3, 1) 0.000000 11.00000 A( 3, 2) 0.000000 0.000000 A( 3, 3) 1.000000 0.000000 A( 3, 4) 0.000000 5.000000 A( 4, 1) 0.000000 1.000000 A( 4, 2) 0.000000 0.000000 A( 4, 3) 0.000000 4.000000 A( 4, 4) 1.000000 0.000000 *( 1, 1) 15.00000 0.000000 *( 1, 2) 18.
5、00000 0.000000 *( 1, 3) 21.00000 0.000000 *( 1, 4) 24.00000 0.000000 *( 2, 1) 19.00000 0.000000 *( 2, 2) 23.00000 0.000000 *( 2, 3) 22.00000 0.000000 *( 2, 4) 18.00000 0.000000 *( 3, 1) 26.00000 0.000000 *( 3, 2) 18.00000 0.000000 *( 3, 3) 16.00000 0.000000 *( 3, 4) 19.00000 0.000000 *( 4, 1) 19.000
6、00 0.000000 *( 4, 2) 21.00000 0.000000 *( 4, 3) 23.00000 0.000000 *( 4, 4) 17.00000 0.000000 Row Slack or Surplus Dual Price1 70.00000 -1.000000 2 0.000000 -14.00000 3 0.000000 -18.00000 4 0.000000 -14.00000 5 0.000000 -17.00000 6 0.000000 -1.000000 7 0.000000 -4.000000 8 0.000000 -2.000000 9 0.0000
7、00 0.000000第1个人完成第2项,第2人完成第1项,第3人完成第3项,第4人完成第4项。最短时间为70.2.模型:model:sets:pe/1,2,3,4,5/;ci/1,2,3,4,5/;lin(pe,ci):*;lin1(pe,pe):d;lin2(ci,ci):c;endsetsmin=sum(lin1(i,j):d(i,j)*sum(lin2(a,b):c(a,b)*(a,i)*(b,j);for(pe(m):sum(ci(n):*(m,n)=1);for(ci(n):sum(pe(m):*(m,n)=1);for(lin:bin(*);data:d=0 1 1 2 3 1
8、0 2 1 2 1 2 0 1 2 2 1 1 0 1 3 2 2 1 0;c=0 5 2 4 1 5 0 3 0 2 2 3 0 0 0 4 0 0 0 5 1 2 0 5 0;enddataend结果:Local optimal solution found. Objective value: 50.00000 Objective bound: 50.00000 Infeasibilities: 0.000000 E*tended solver steps: 18 Total solver iterations: 802Variable Value Reduced Cost *( 1, 1
9、) 0.000000 14.00000 *( 1, 2) 0.000000 0.000000 *( 1, 3) 1.000000 5.999996 *( 1, 4) 0.000000 0.000000 *( 1, 5) 0.000000 6.000000 *( 2, 1) 0.000000 16.00000 *( 2, 2) 0.000000 9.999992 *( 2, 3) 0.000000 0.000000 *( 2, 4) 1.000000 1.999992 *( 2, 5) 0.000000 0.000000 *( 3, 1) 0.000000 14.00000*( 3, 2) 0.
10、000000 6.000000 *( 3, 3) 0.000000 0.8523365E-08 *( 3, 4) 0.000000 0.000000 *( 3, 5) 1.000000 0.000000 *( 4, 1) 1.000000 7.999996 *( 4, 2) 0.000000 0.000000 *( 4, 3) 0.000000 5.999996 *( 4, 4) 0.000000 6.000000 *( 4, 5) 0.000000 14.00000 *( 5, 1) 0.000000 0.000000 *( 5, 2) 1.000000 1.999981 *( 5, 3)
11、0.000000 0.000000 *( 5, 4) 0.000000 9.999981 *( 5, 5) 0.000000 15.99998 D( 1, 1) 0.000000 0.000000 D( 1, 2) 1.000000 0.000000 D( 1, 3) 1.000000 0.000000 D( 1, 4) 2.000000 0.000000 D( 1, 5) 3.000000 0.000000 D( 2, 1) 1.000000 0.000000 D( 2, 2) 0.000000 0.000000 D( 2, 3) 2.000000 0.000000 D( 2, 4) 1.0
12、00000 0.000000 D( 2, 5) 2.000000 0.000000 D( 3, 1) 1.000000 0.000000 D( 3, 2) 2.000000 0.000000 D( 3, 3) 0.000000 0.000000 D( 3, 4) 1.000000 0.000000 D( 3, 5) 2.000000 0.000000 D( 4, 1) 2.000000 0.000000 D( 4, 2) 1.000000 0.000000 D( 4, 3) 1.000000 0.000000 D( 4, 4) 0.000000 0.000000 D( 4, 5) 1.0000
13、00 0.000000 D( 5, 1) 3.000000 0.000000 D( 5, 2) 2.000000 0.000000 D( 5, 3) 2.000000 0.000000 D( 5, 4) 1.000000 0.000000 D( 5, 5) 0.000000 0.000000 C( 1, 1) 0.000000 0.000000 C( 1, 2) 5.000000 0.000000 C( 1, 3) 2.000000 0.000000 C( 1, 4) 4.000000 0.000000 C( 1, 5) 1.000000 0.000000 C( 2, 1) 5.000000
14、0.000000 C( 2, 2) 0.000000 0.000000 C( 2, 3) 3.000000 0.000000 C( 2, 4) 0.000000 0.000000 C( 2, 5) 2.000000 0.000000 C( 3, 1) 2.000000 0.000000 C( 3, 2) 3.000000 0.000000 C( 3, 3) 0.000000 0.000000 C( 3, 4) 0.000000 0.000000 C( 3, 5) 0.000000 0.000000 C( 4, 1) 4.000000 0.000000 C( 4, 2) 0.000000 0.0
15、00000 C( 4, 3) 0.000000 0.000000 C( 4, 4) 0.000000 0.000000 C( 4, 5) 5.000000 0.000000C( 5, 1) 1.000000 0.000000 C( 5, 2) 2.000000 0.000000 C( 5, 3) 0.000000 0.000000 C( 5, 4) 5.000000 0.000000 C( 5, 5) 0.000000 0.000000 Row Slack or Surplus Dual Price1 50.00000 -1.000000 2 0.000000 -10.00000 3 0.00
16、0000 -6.000000 4 0.000000 8.000000 5 0.000000 0.000000 6 0.000000 -0.8674480E-05 7 0.000000 -10.00000 8 0.000000 -16.00000 9 0.000000 -14.00000 10 0.000000 -12.00000 11 0.000000 -22.00000第一个人去第三个城市,第二个人去第四个城市,第三个人去第五个城市,第四个人去第一个城市,第五个人去第二个城市,最低费用50。3.模型:model: min=100*1+100*2+40*y1+40*y2+40*y3+40*y4
17、+40*y5;*1+*2+y14;*1+*2+y1+y23;*1+*2+y1+y2+y34;*2+y1+y2+y3+y46;*1+y2+y3+y4+y55;*1+*2+y3+y4+y56;*1+*2+y4+y58;*1+*2+y58;y1+y2+y3+y4+y53;gin(*1);gin(*2);gin(y1);gin(y2);gin(y3);gin(y4);gin(y5);endGlobal optimal solution found. Objective value: 820.0000 Objective bound: 820.0000 Infeasibilities: 0.000000
18、 E*tended solver steps: 0 Total solver iterations: 46 Variable Value Reduced Cost *1 3.000000 100.0000 *2 4.000000 100.0000 Y1 0.000000 40.00000 Y2 2.000000 40.00000 Y3 0.000000 40.00000 Y4 0.000000 40.00000 Y5 1.000000 40.00000 Row Slack or Surplus Dual Price 1 820.0000 -1.000000 2 3.000000 0.00000
19、0 3 6.000000 0.000000 4 5.000000 0.000000 5 0.000000 0.000000 6 1.000000 0.000000 7 2.000000 0.000000 8 0.000000 0.000000 9 0.000000 0.000000 10 0.000000 0.000000说明:该储蓄雇用7名全时效劳员,其中3名12:0013:00吃午餐,4名13:0014:00吃午餐;并在10:0014:00雇用2名、13:0017:00雇用1名半时效劳员,符合题目要求,且花费最少,为820元。第二问 y1+y2+y3+y4+y50;Global optimal solution found. Objective value: 560.0000 Objective bound: 560.0000 Infeasibilities: 0.000000 E*tended solver steps: 0 Total solver iterations: 2 Variable Value Reduced Cost *1 0.000000 100.0000 *2 0.000000 100.0000