| 查看: 1012 | 回复: 6 | |||
| 【悬赏金币】回答本帖问题，作者s1y9shy将赠送您 5 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
lingo 求解选址问题 本身金币太少 所以给不了多少悬赏 望见谅 真心求助已有2人参与
|
|||
|
model: sets： factory/1．．7/：a; warhouse/l．．4/：v,f,z,b； customer/1．．1 1/：c； linkl(factory,warhouse)：m,x; link2(warhouse,customer)：n,y； endsets data： a=130 24 82 60 1130 50 330; c=96.42 353.21 174.19 102.67 51.67 282.97 230.55 53.415 198.14 77.00 85.74; v=3 2 4 5; f=4000 3000 2000 1000； b=1000 1000 1000 1000； m=O 228 203 434 334 120 280 228 0 361 244 102 151 210 280 210 483 162 123 330 0 334 102 472 133 0 281 123； n=203 361 472 483 434 244 133 162 0 228 334 280 141 299 421 421 491 389 322 199 228 0 102 210 142 199 192 138 577 387 288 271 316 169 74 49 120 151 281 281 522 403 422 238； enddata min=@sum(warhouse(j)：@sum(factory(i)：m(i，j)*x(i,j)*z(j))) +@sum customer(k)：@sum(warhouse(j)：n(j，k)*y(i,j)*z(j)))+@sum (warhouse(j)：@sum(factory(i)：v(j)**x(i,j)*z(j)))+@sum(warhouse(j)：f(j)*z (j))； @for(factory(i)：@sum(warhouse(j)：x(i,j))<=a(i))； @for(warhous(j)：@sum(link2(j，k)：y(j，k)) =@sum(1ink1(i,j)：x(i,j)))； @for(customer(k)：@sum(warhouse(i)：y(j，k))>=c(k))； @for(warhouse：@bin(z))； @for(warhouse(j)：@sum(factory(i)：x(i,j))<=99999999*z(j)) @for(warhouse(j)：@sum(factory(i)：x(i,j))<=b(j))； end 这个程序老是运行不出来 这个是一个线性的 还有一个问题也希望得到您的帮助 若是在目标函数中添入一个决策变量开根号 也就是非线性还能求解么 谢谢 |
回复此楼» 猜你喜欢
大佬们，打的核磁氢谱与chemdraw预测的有偏差
已经有6人回复
-
请问有知道博士后面上专家评审什么时候结束的么？
已经有7人回复
-
为何全世界都用JCR分区，咱们要弄个中科院分区？意义何在？
已经有4人回复
-
offer选择
已经有10人回复
-
论文被硕导卡
已经有6人回复
-
三颈瓶加热
已经有5人回复
-
请问三个字的中文名字，在文中到底怎么表示？
已经有7人回复
-
高校辞职，要求赔偿，这到底合不合理
已经有36人回复
-
博后面上今天下午会公布吗？大家有无消息？
已经有14人回复
-
我在写关于多维空间的论文，希望能与大家交流，准备在science上发表
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有学lingo，会lingo的同学吗？
已经有3人回复
NicoleLas
木虫 (正式写手)
- 应助: 66 (初中生)
- 金币: 4288.7
- 红花: 26
- 帖子: 514
- 在线: 176.9小时
- 虫号: 2872321
- 注册: 2013-12-15
- 性别: GG
- 专业: 计算机软件
【答案】应助回帖
感谢参与，应助指数 +1
|
我没有具体看模型，修改后至少有一个结果： 代码： model: sets: factory/1..7/:a; warhouse/1..4/:v,f,z,b; customer/1..11/:c; link1(factory,warhouse):m,x; link2(warhouse,customer):n,y; endsets data: a=130 24 82 60 1130 50 330; c=96.42 353.21 174.19 102.67 51.67 282.97 230.55 53.415 198.14 77.00 85.74; v=3 2 4 5; f=4000 3000 2000 1000; b=1000 1000 1000 1000; m=0 228 203 434 334 120 280 228 0 361 244 102 151 210 280 210 483 162 123 330 0 334 102 472 133 0 281 123; n=203 361 472 483 434 244 133 162 0 228 334 280 141 299 421 421 491 389 322 199 228 0 102 210 142 199 192 138 577 387 288 271 316 169 74 49 120 151 281 281 522 403 422 238; enddata min=@sum(link1(i,j):m(i,j)*x(i,j)*z(j))+@sum(link2(j,k):n(j,k)*y(j,k)*z(j)) +@sum(link1(i,j):v(j)*x(i,j)*z(j))+@sum(warhouse(j):f(j)*z(j)); @for(factory(i):@sum(warhouse(j):x(i,j))<=a(i)); @sum(link2(j,k):y(j,k))=@sum(link1(i,j):x(i,j)); @for(customer(k):@sum(warhouse(j):y(j,k))>=c(k)); @for(warhouse:@bin(z)); @for(warhouse(j):@sum(factory(i):x(i,j))<=99999999*z(j)); @for(warhouse(j):@sum(factory(i):x(i,j))<=b(j)); end LINGO 14求解结果（部分）： Global optimal solution found. Objective value: 154529.9 Objective bound: 154529.9 Infeasibilities: 0.000000 Extended solver steps: 0 Total solver iterations: 41 Elapsed runtime seconds: 0.05 Model Class: MIQP Total variables: 76 Nonlinear variables: 74 Integer variables: 4 Total constraints: 28 Nonlinear constraints: 1 Total nonzeros: 278 Nonlinear nonzeros: 70 你可以使用网络流行的LINGO 9,10,11任一Crack Extended Version求解出相同结果。 |
6楼2014-04-17 23:07:45
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2014-04-15 16:20:24
3楼2014-04-15 19:10:12
4楼2014-04-15 19:12:30
5