QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2009|回复: 2
打印 上一主题 下一主题

关于lingo代码,高手指点

[复制链接]
字体大小: 正常 放大
kenstorm        

5

主题

3

听众

281

积分

升级  90.5%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2009-5-12 03:31 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
MODEL:
Titel Location problem;
sets:
      demand/1..6/:a,b,d;
      supply/1..2/:x,y,e;
      Link(demand,supply):c;
endsets
data:
!Locations for the demand(需求点的位置);
a=1.25,8.75,0.5,5.75,3,7.25;
b=1.25,0.75,4.75,5,6.5,7.75;
!quantites of the demand and supply(供需量);
d=3,5,4,7,6,11;e=20,20;
enddata
init:
!initial Locations for the supply(初始点);
x,y=5 1 2 7;
endinit
!Objective function(目标);
[OBJ]min=@sum(Link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2));
!demand constraints(需求约束);
@for(demand(i):[DEMAND_CON]@sum(supply(j):c(i,j))=d(i););
!supply constraints(供应约束);
@for(supply(i):[SUPPLY_CON]@sum(demand(j):c(i,j))<=e(i););
@for(supplyfree(X);@free(Y););
END





model:
sets:
    demand/1..6/:a,b,d;
    supply/1..2/:x,y,e;
    link(demand,supply):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;e=20,20;
enddata
init:
x,y=5,1,2,7;
endinit
[OBJ]min=@sum(link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2));
@for(demand(i):[demand_con]@sum(supply(j):c(i,j))=d(i););
@for(supply(i):[supply_con]@sum(supply(j):c(i,j))<=e(i););
@for(supplyfree(x);@free(y););
end
以上两段代码,前者无法运算,但是我看不出有什么错误的地方,请高手指点,还有,后面一个如果加上Titel Location problem;也不能求解,请问加标题正确的代码是怎样的,谢谢啦
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
风口        

1

主题

4

听众

56

积分

升级  53.68%

该用户从未签到

新人进步奖

回复

使用道具 举报

kiddyse        

0

主题

3

听众

22

积分

升级  17.89%

该用户从未签到

新人进步奖

似乎是谢金星老师书上的源程序…我把程序粘出来吧…
MODEL:
Title Location Problem;
sets:
            demand/1..6/:a,b,d;
    supply/1..2/:x,y,e;
    link(demand,supply):c;
endsets
data:
!locations for the demand(需求点的位置);
a=1.25,8.75,0.5,5.75,3,7.25;
b=1.25,0.75,4.75,5,6.5,7.75;
!quantities of the demand and supply(供需量);
d=3,5,4,7,6,11; e=20,20;
enddata
init:
!initial locations for the supply(初始点);
x,y=5,1,2,7;
endinit
!Objective function(目标);
[OBJ] min=@sum(link(i,j): c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2) );
!demand constraints(需求约束);
@for(demand(i):[DEMAND_CON] @sum(supply(j):c(i,j)) =d(i););
!supply constraints(供应约束);
@for(supply(i):[SUPPLY_CON] @sum(demand(j):c(j,i)) <=e(i); );
@for(supply: @bnd(0.5,X,8.75); @bnd(0.75,Y,7.75); );
END
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2025-8-14 11:22 , Processed in 0.502667 second(s), 63 queries .

回顶部