LINGO應(yīng)用實(shí)例

1、LINGO應(yīng)用實(shí)例班級(jí)：2011級(jí)軟件工程1班；姓名：林保京；學(xué)號(hào):P111713222某公司有6 個(gè)建筑工地要開工，每個(gè)工地的位置（用平面坐標(biāo)a, b 表示，距離單位：公里）及水泥日用量d（噸）由下表給出。目前有兩個(gè)臨時(shí)料場(chǎng)位于P (5, 1), Q (2, 7) ，日儲(chǔ)量各有20 噸。假設(shè)從料場(chǎng)到工地之間均有直線道路相連，試制定每天的供應(yīng)計(jì)劃，即從A, B 兩料場(chǎng)分別向各工地運(yùn)送多少噸水泥，使總的噸公里數(shù)最小。為了進(jìn)一步減少噸公里數(shù)，打算舍棄兩個(gè)臨時(shí)料場(chǎng)，改建兩個(gè)新的，日儲(chǔ)量仍各為20 噸，問應(yīng)建在何處，節(jié)省的噸公里數(shù)有多大。記工地的位置為(ai,bi)，水泥的日用量為di ,i=16;料

2、場(chǎng)的位置為(xj,yj),日儲(chǔ)量為ej,j=1,2;從料場(chǎng)j到工地i的運(yùn)送量為cij,這個(gè)優(yōu)化問題的數(shù)學(xué)規(guī)劃模型是：minf=j-12i-16cijxj-ai2+yj-bi2st j-12cij=di,i=16i-16cijej,j=1,2將其轉(zhuǎn)化成LINGO代碼：MODEL:Title Location Problem;sets:demand/1.6/:a,b,d;supply/1.2/:x,y,e;link(demand,supply):c;endsetsdata:!需求點(diǎn)設(shè)置;a=1.25,8.75,0.5,5.75,3,7.25;b=1.25,0.75,5.75,5,6.5,7.75;

3、!需求量;d=3,5,4,7,6,11;e=20,20;enddatainit:!初始點(diǎn);x,y=5,1,2,7;endinitOBJ min=sum(link(i,j):c(i,j)*(x(j)-a(i)2+(y(j)-b(i)2)(1/2);for(demand(i):DEMAND_CONsum(supply(j):c(i,j)=d(i););for(supply(i):SUPPLY_CONsum(demand(j):c(j,i)=e(i););for(supply:free(X); free(Y););END運(yùn)行結(jié)果如下：Local optimal solution found. Obj

4、ective value: 84.41406 Infeasibilities: 0.000000 Total solver iterations: 132 Model Title: Location Problem Variable Value Reduced Cost A( 1) 1.250000 0.000000 A( 2) 8.750000 0.000000 A( 3) 0.5000000 0.000000 A( 4) 5.750000 0.000000 A( 5) 3.000000 0.000000 A( 6) 7.250000 0.000000 B( 1) 1.250000 0.00

5、0000 B( 2) 0.7500000 0.000000 B( 3) 5.750000 0.000000 B( 4) 5.000000 0.000000 B( 5) 6.500000 0.000000 B( 6) 7.750000 0.000000 D( 1) 3.000000 0.000000 D( 2) 5.000000 0.000000 D( 3) 4.000000 0.000000 D( 4) 7.000000 0.000000 D( 5) 6.000000 0.000000 D( 6) 11.00000 0.000000 X( 1) 3.113090 0.000000 X( 2)

6、7.250000 -0.1954955E-06 Y( 1) 6.044738 0.000000 Y( 2) 7.750000 -0.8808641E-08 E( 1) 20.00000 0.000000 E( 2) 20.00000 0.000000 C( 1, 1) 3.000000 0.000000 C( 1, 2) 0.000000 3.701915 C( 2, 1) 0.000000 0.5747176 C( 2, 2) 5.000000 0.000000 C( 3, 1) 4.000000 0.000000 C( 3, 2) 0.000000 4.410404 C( 4, 1) 7.

7、000000 0.000000 C( 4, 2) 0.000000 0.2961603 C( 5, 1) 6.000000 0.000000 C( 5, 2) 0.000000 3.960913 C( 6, 1) 0.000000 4.474589 C( 6, 2) 11.00000 0.000000 Row Slack or Surplus Dual Price OBJ 84.41406 -1.000000 DEMAND_CON( 1) 0.000000 -5.143988 DEMAND_CON( 2) 0.000000 -7.158911 DEMAND_CON( 3) 0.000000 -2.629659 DEMAND_CON( 4) 0.000000 -2.836331 DEMAND_CON( 5) 0.000000 -0.4690982 DE