之三(Lingo中的基本集合与派生集合).doc

三．Lingo中的基本集合与派生集合 例4：料场选址问题 六个建筑工地的位置（用平面坐标a、b表示，距离单位：km）及其对水泥的日用量（用d表示，单位：t）由下表给出： 工地的位置（a，b）及 水泥日用量d -------------------------------------------------------------------------------------------------- 工地编号 1 2 3 4 5 6 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 -------------------------------------------------------------------------------------------------- 现有两个临时料场位于P(5,1), Q(2,7)，每日提供水泥的最大能力分别为20t. 假设从料场到各工地均有直线道路连接，运输费用与距离、重量成正比。 （1）请制定运输计划，使总运费尽量低。 （2）进一步调整这两个临时料场的位置，使总运费最低。 解：第i号工地：位置，水泥日用量， i=1,2,3,4,5,6. 第j号料场：位置，水泥日供应能力， j=1,2. 从j号料场向i号工地的日运输水泥量记为 . 注意：在问题（1）中，为已知数据，故决策变量为，共12个； 在问题（2）中，待定，故决策变量为，共16个。 从j号料场到i号工地，距离为，送去重量为的水泥，二者乘积即为运输费。 目标函数： 约束条件： 满足需求： 供应能力： ，j=1,2 非负性： 这就是本题的优化模型。 尝试用Lingo求解该模型时，六个建筑工地作为一个集合gdjh，两个料场作为一个集合lcjh。接下来就会遇到困难：决策变量不仅是依赖于集合gdjh的属性，而且是依赖于集合lcjh的属性，这样的属性应该如何定义呢？ 根据两个基本集合gdjh与lcjh构造一个派生集合gdlcjh，再把定义为这个集合的属性。 先看本题第(1)问的Lingo程序，再看注解： model: sets: gdjh/1..6/:a,b,d; lcjh/1,2/:x,y,e; gdlcjh(gdjh,lcjh):c; endsets data: 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; x,y=5,1,2,7;e=20,20; enddata min=@sum(gdlcjh(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^0.5); @for(gdjh(i):@sum(lcjh(j):c(i,j))=d(i)); @for(lcjh(j):@sum(gdjh(i):c(i,j))=e(j)); end 注解：（1）定义派生集合及其属性的命令格式为 派生集合名（基本集合1，基本集合2）：属性 （2）赋值语句“x,y=5,1,2,7”的赋值顺序是“x(1)=5,y(1)=1,x(2)=2,y(2)=7”，而不是“x(1),x(2),y(1),y(2)”；在程序中，该语句可换成“x=5,2;y=1,7”，功能相同。 （3）当表达式中出现的下标符号多于1个时，必须指明针对哪个符号做运算。 执行结果： Global optimal solution found at iteration: 4 Objective value: 136.2275 Variable Value Reduced Cost C( 1, 1) 3.000000 0.000000 C( 1, 2) 0.000000 3.85220