Lingo

計算機學(xué)院張亞玲數(shù)學(xué)規(guī)劃軟件某公司有6個建筑工地，位置坐標為(ai,bi)(單位：公里)，水泥日用量di

(單位：噸）假設(shè)：料場和工地之間有直線道路（1）現(xiàn)有2料場，位于A(5,1)，B(2,7)，記(xj,yj)，j=1,2，日儲量ej各有20噸。制定每天的供應(yīng)計劃：即從A,B兩料場分別向各工地運送多少噸水泥，使總的噸公里數(shù)最小。i１２３４５６a２8.750.55.7537.25b３0.754.7556.57.75d1.25547611例選址問題決策目標約束解：決策變量：料場j到工地i的運量——

cij

——12維線性規(guī)劃模型？lingo——表達式？目標：噸公里約束：需求供應(yīng)建筑工地位置坐標(ai,bi)、水泥日用量di

：對每個建筑工地（6個）都有一個對應(yīng)的值都是一個由6個元素組成的數(shù)組——是已知的料場位置坐標(xj,yj)、日儲量ej對每個料場（2個）都有一個對應(yīng)的值都是一個由2個元素組成的數(shù)組——目前是已知的料場到建筑工地的供應(yīng)計劃cij

對每個料場與建筑工地之間（6×2）都有一個對應(yīng)的值是一個6×2個元素組成的矩陣——是未知數(shù)LINDO無數(shù)組，每個變量輸入——麻煩（1）Lingo的集合Set——下標集合——100個工地？1、Lingo的集合Set及其屬性（2）集合Set及其屬性定義數(shù)組下標集合demand/1..6/———表示6個建筑工地a,b,d稱為該集合的屬性———表示坐標(ai,bi)、水泥日用量di定義數(shù)組下標集合supply/1..2/———表示6個建筑工地該集合的屬性x,y,e———表示坐標(xj,yj)、日儲量ej

定義數(shù)組下標集合link(demand,supply)———表示6×2個料場到建筑工地的連接該集合的屬性c———表示每個料場與建筑工地之間供應(yīng)計劃cij1到6的整數(shù)建立下標集合（3）Lingo建模語言——集合段數(shù)據(jù)段需求點的位置供需量sets:demand/1..6/:a,b,d;supply/1..2/:x,y,e;link(demand,supply):c;endsetsdata:a=1.25,8.75,0.5,5.75,3,7.25;b=1.25,0.75,4.75,5,6.5,7.75;d=3,5,4,7,6,11;e=20,20;x,y=5,1,2,7;enddata賦值需求供應(yīng)連接LINGO建模語言也稱為矩陣生成器（MATRIXGENERATOR）。類似DEMAND和SUPPLY直接把元素列舉出來的集合，稱為基本集合(primaryset),而把LINK這種基于其它集合而派生出來的二維或多維集合稱為派生集合(derivedset)。由于是DEMAND和SUPPLY生成了派生集合LINK，所以DEMAND和SUPPLY稱為LINK的父集合。建立下標集合例3

選址問題需求點的位置供需量sets:demand/1..6/:a,b,d;supply/1..2/:x,y,e;link(demand,supply):c;endsetsdata:a=1.25,8.75,0.5,5.75,3,7.25;b=1.25,0.75,4.75,5,6.5,7.75;d=3,5,4,7,6,11;e=20,20;x,y=5,1,2,7;enddata賦值需求供應(yīng)連接基本集合派生集合目標：噸公里例3選址問題min=@sum(link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2));sets:demand/1..6/:a,b,d;supply/1..2/:x,y,e;link(demand,supply):c;endsets約束：非負@for(demand(i):@sum(supply(j):c(i,j))=d(i););@for(supply(j):@sum(demand(i):c(i,j))<=e(j););約束：

需求

供應(yīng)OK目標與約束段MODEL:TitleLocationProblem;sets:demand/1..6/:a,b,d;supply/1..2/:x,y,e;link(demand,supply):c;endsetsdata:a=1.25,8.75,0.5,5.75,3,7.25;b=1.25,0.75,4.75,5,6.5,7.75;d=3,5,4,7,6,11;e=20,20;x,y=5,1,2,7;enddatainit:endinitmin=@sum(link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2));@for(demand(i):@sum(supply(j):c(i,j))=d(i););@for(supply(j):@sum(demand(i):c(i,j))<=e(j););END集合段數(shù)據(jù)段初始段供應(yīng)約束需求需求點的位置供需量供應(yīng)初始點目標需求約束連接結(jié)果：總噸公