- 在线时间
- 72 小时
- 最后登录
- 2016-1-20
- 注册时间
- 2011-9-25
- 听众数
- 4
- 收听数
- 0
- 能力
- 0 分
- 体力
- 802 点
- 威望
- 0 点
- 阅读权限
- 40
- 积分
- 322
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 163
- 主题
- 8
- 精华
- 0
- 分享
- 1
- 好友
- 18
升级 7.33% TA的每日心情 | 开心 2015-9-17 11:12 |
---|
签到天数: 86 天 [LV.6]常住居民II
|
lingo程序总是出现unexpected jacobian overflow的错误,不知为什么?
我的程序为:
model:
title location problem;
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)*((a(i)-x(j))^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(demand(j):c(j,i)) <=e(i););
@for(supply: @bnd(0.5,X,8.75); @bnd(0.75,Y,7.75); );
end
答案标准程序为:
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
请问下各位大侠,我的程序哪儿出错了?和答案程序有什么不同? |
zan
|