北交大LINGO考试1 转运问题等.docx

物流软件实验报告 实验编号： 学号： 序号： 姓名： 班级： 2015年10月一，用lingo解决运输问题（1） 建立数学模型设Cij为从工厂Ai销售到销售点Bj运费单价，Xij为从Ai到B觉得运输量，因此总运费为第i个产地的运出量应小于等于该地的需求量第j个销售地的运入量应等于该地的需求量因此得到数学模型如下：Min；s.t. Xij=0, i=1,2,m, j=1,2,n.（2） 编写相关的模型如下MODEL:SETS:warehouse/1,2,3,4/:a;customer/1,2,3,4,5/:b;ROUTES(WAREHOUSE,CUSTOMER):c,x;ENDSETSDATA:a=20,15,25,30;b=25,20,15,20,10;c=3,11,3,10,4, 1,9,2,8,3, 7,4,10,5,6, 5,8,2,6,2;enddataOBJmin=sum(ROUTES:c*x);FOR(WAREHOUSE(I):SUPSUM(CUSTOMER(J):X(I,J)=0 X2=0（2）用lingo编程MODEL:SETS:WAREHOUSE/W1,W2/:PRODUCE;CENTER/P1,P2,P3/;CUSTOMER/D1,D2,D3,D4/:REQUIRE;LINKI(WAREHOUSE,CENTER):CI,XI;LINKII(CENTER,CUSTOMER):CII,XII;ENDSETSDATA:PRODUCE=11,9;REQUIRE=3,6,5,6;CI=13,12,15, 11,13,22;CII=3,11,3,10, 1,9,2,8, 7,4,10,5;ENDDATAobjmin=sum(linki:ci*xi)+sum(linkii:cii*xii);for(WAREHOUSE(i):supsum(center(j):xi(i,j)=produce(i);for(center(j):midsum(WAREHOUSE(i):xi(i,j) = sum(customer(k):xii(j,k);for(customer(k):demsum(center(j):xii(j,k)=require(k);End！LINKI和LINKII都是派生集合名，CI,CII是属性，表对应的运费，XI,XII也是属性，表对应的决策运量。!CUSTOMER是顾客需求集合名，D1，D2，D3，D4是元素，require是属性！CENTER是物流中心集合名，P1，P2，P3是元素。！WAREHOUSE是仓库集合名W1，W2是元素，PRODUCE是属性（3） 运行结果如下Global optimal solution found. Objective value: 347.0000 Total solver iterations: 0 Variable Value Reduced Cost PRODUCE( W1) 11.00000 0.000000 PRODUCE( W2) 9.000000 0.000000 REQUIRE( D1) 3.000000 0.000000 REQUIRE( D2) 6.000000 0.000000 REQUIRE( D3) 5.000000 0.000000 REQUIRE( D4) 6.000000 0.000000 CI( W1, P1) 13.00000 0.000000 CI( W1, P2) 12.00000 0.000000 CI( W1, P3) 15.00000 0.000000 CI( W2, P1) 11.00000 0.000000 CI( W2, P2) 13.00000 0.000000 CI( W2, P3) 22.00000 0.000000 XI( W1, P1) 0.000000 3.000000 XI( W1, P2) 5.000000 0.000000 XI( W1, P3) 6.000000 0.000000 XI( W2, P1) 9.000000 0.000000 XI( W2, P2) 0.000000 0.000000 XI( W2, P3) 0.000000 6.000000 CII( P1, D1) 3.000000 0.000000 CII( P1, D2) 11.00000 0.000000 CII( P1, D3) 3.000000 0.000000 CII( P1, D4) 10.00000 0.000000 CII( P2, D1) 1.000000 0.000000 CII( P2, D2) 9.000000 0.000000 CII( P2, D3) 2.000000 0.000000 CII( P2, D4) 8.000000 0.000000 CII( P3, D1) 7.000000 0.000000 CII( P3, D2) 4.000000 0.000000 CII( P3, D3) 10.00000 0.000000 CII( P3, D4) 5.000000 0.000000 XII( P1, D1) 3.000000 0.000000 XII( P1, D2) 0.000000 2.000000 XII( P1, D3) 5.000000 0.000000 XII( P1, D4) 1.000000 0.000000 XII( P2, D1) 0.000000 0.000000 XII( P2, D2) 0.000000 2.000000 XII( P2, D3) 0.000000 1.000000 XII( P2, D4) 5.000000 0.000000 XII( P3, D1) 0.000000 9.000000 XII( P3, D2) 6.000000 0.000000 XII( P3, D3) 0.000000 12.00000 XII( P3, D4) 0.000000 0.000000 Row Slack or Surplus Dual Price OBJ 347.0000 -1.000000 SUP( W1) 0.000000 1.000000 SUP( W2) 0.000000 0.000000 MID( P1) 0.000000 -11.00000 MID( P2) 0.000000 -13.00000 MID( P3) 0.000000 -16.00000 DEM( D1) 0.000000 -14.00000 DEM( D2) 0.000000 -20.00000 DEM( D3) 0.000000 -14.00000 DEM( D4) 0.000000 -21.00000从结果可以看出，W1向P1，P2，P3分别运送0，5，6； W2向P1，P2，P3分别运送9，0，0。 P1 向D1，D2，D3，D4,分别配送3，0，5，1；P2 向D1，D2，D3，D4分别配送0，0，0，5；P3 向D1，D2，D3，D4分别配送0，6，0，0。最少运费为347元/单位（4） 改变以上的扩建规模，经过分析，每个仓库最多可以扩建2单位。有9种可能性1， W1+1，W2+0（无解）2， W1+1，W2+1（无解）3， W1+1，W2+2（388）4， W1+2，W2+0（无解）5， W1+2，W2+1（387）6， W1+2，W2+2（387）7， W1+0，W2+0（无解）8， W1+0，W2+1（无解）9， W1+0，W2+2（无解）应该选择方案5，即扩建W1两单位，扩建W2一单位。相应程序如下MODEL:SETS:WAREHOUSE/W1,W2/:PRODUCE;CENTER/P1,P2,P3/;CUSTOMER/D1,D2,D3,D4/:REQUIRE;LINKI(WAREHOUSE,CENTER):CI,XI;LINKII(CENTER,CUSTOMER):CII,XII;ENDSETSDATA:PRODUCE=13,11;REQUIRE=6,6,5,6;CI=13,12,15, 11,13,22;CII=3,11,3,10, 1,9,2,8, 7,4,10,5;ENDDATAobjmin=sum(linki:ci*xi)+sum(linkii:cii*xii);for(WAREHOUSE(i):supsum(center(j):xi(i,j)=produce(i);for(center(j):midsum(WAREHOUSE(i):xi(i,j) = sum(customer(k):xii(j,k);for(customer(k):demsum(center(j):xii(j,k)=require(k);end结果运行Global optimal solution found. Objective value: 387.0000 Total solver iterations: 11 Variable Value Reduced Cost PRODUCE( W1) 13.00000 0.000000 PRODUCE( W2) 11.00000 0.000000 REQUIRE( D1) 6.000000 0.000000 REQUIRE( D2) 6.000000 0.000000 REQUIRE( D3) 5.000000 0.000000 REQUIRE( D4) 6.000000 0.000000 CI( W1, P1) 13.00000 0.000000 CI( W1, P2) 12.00000 0.000000 CI( W1, P3) 15.00000 0.000000 CI( W2, P1) 11.00000 0.000000 CI( W2, P2) 13.00000 0.000000 CI( W2, P3) 22.00000 0.000000 XI( W1, P1) 0.000000 3.000000 XI( W1, P2) 7.000000 0.000000 XI( W1, P3) 6.000000 0.000000 XI( W2, P1) 10.00000 0.000000 XI( W2, P2) 0.000000 0.000000 XI( W2, P3) 0.000000 6.000000 CII( P1, D1) 3.000000 0.000000 CII( P1, D2) 11.00000 0.000000 CII( P1, D3) 3.000000 0.000000 CII( P1, D4) 10.00000 0.000000 CII( P2, D1) 1.000000 0.000000 CII( P2, D2) 9.000000 0.000000 CII( P2, D3) 2.000000 0.000000 CII( P2, D4) 8.000000 0.000000 CII( P3, D1) 7.000000 0.000000 CII( P3, D2) 4.000000 0.000000 CII( P3, D3) 10.00000 0.000000 CII( P3, D4) 5.000000 0.000000 XII( P1, D1) 5.000000 0.000000 XII( P1, D2) 0.000000 2.000000 XII( P1, D3) 5.000000 0.000000 XII( P1, D4) 0.000000 0.000000 XII( P2, D1) 1.000000 0.000000 XII( P2, D2) 0.000000 2.000000 XII( P2, D3) 0.000000 1.000000 XII( P2, D4) 6.000000 0.000000 XII( P3, D1) 0.000000 9.000000 XII( P3, D2) 6.000000 0.000000 XII( P3, D3) 0.000000 12.00000 XII( P3, D4) 0.000000 0.000000 Row Slack or Surplus Dual Price OBJ