24小时热门版块排行榜    

查看: 1005  |  回复: 6
【悬赏金币】回答本帖问题,作者s1y9shy将赠送您 5 个金币

s1y9shy

铁虫 (小有名气)

[求助] 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


这个程序老是运行不出来  这个是一个线性的  
还有一个问题也希望得到您的帮助    若是在目标函数中添入一个决策变量开根号  也就是非线性还能求解么  谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
原问题用文字描述下,包括各数据的意思。
2楼2014-04-15 16:20:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

s1y9shy

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2014-04-15 16:20:24
原问题用文字描述下,包括各数据的意思。

是这篇文献里面的,我的模型比它这个还要复杂一点,但是他这个我都求解不出来,麻烦了~~谢谢您呀
lingo 求解选址问题 本身金币太少 所以给不了多少悬赏  望见谅  真心求助
QQ图片20140415190813.jpg


lingo 求解选址问题 本身金币太少 所以给不了多少悬赏  望见谅  真心求助-1
QQ图片20140415190824.jpg


lingo 求解选址问题 本身金币太少 所以给不了多少悬赏  望见谅  真心求助-2
QQ图片20140415190827.jpg


lingo 求解选址问题 本身金币太少 所以给不了多少悬赏  望见谅  真心求助-3
QQ图片20140415190832.jpg


lingo 求解选址问题 本身金币太少 所以给不了多少悬赏  望见谅  真心求助-4
QQ图片20140415190836.jpg


lingo 求解选址问题 本身金币太少 所以给不了多少悬赏  望见谅  真心求助-5
QQ图片20140415190839.jpg


lingo 求解选址问题 本身金币太少 所以给不了多少悬赏  望见谅  真心求助-6
QQ图片20140415190842.jpg

3楼2014-04-15 19:10:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

s1y9shy

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2014-04-15 16:20:24
原问题用文字描述下,包括各数据的意思。

得从最下面一个图往上看,我传的时候弄反了,修改不了了,谢谢呀
4楼2014-04-15 19:12:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

s1y9shy

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by s1y9shy at 2014-04-15 19:10:12
是这篇文献里面的,我的模型比它这个还要复杂一点,但是他这个我都求解不出来,麻烦了~~谢谢您呀

QQ图片20140415190813.jpg

QQ图片20140415190824.jpg

QQ图片20140415190827.jpg

QQ图片20140415190832.jpg

QQ ...

得从最后一个图反着往上看~~上传图片的时候把顺序弄错了,无法再修改了,麻烦啦
5楼2014-04-15 19:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +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求解出相同结果。
Be water, my friend.
6楼2014-04-17 23:07:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

s1y9shy

铁虫 (小有名气)

引用回帖:
6楼: Originally posted by NicoleLas at 2014-04-17 23:07:45
我没有具体看模型,修改后至少有一个结果:
代码:
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
d ...

好的,我看看,我最近也没来上网,可不可以留一下联系方式,我最近写论文老得求解,有很多问题,希望您能帮助我,非常感谢,行么?我是真不认识人会算法的~~希望您可以抽空帮帮我,一定会很感恩的!!
7楼2014-04-26 14:32:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 s1y9shy 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 这个博士要读吗 +14 Sea Breeze 2024-06-16 24/1200 2024-06-20 09:25 by QHJ100100
[精细化工] 切削液交流群 +4 SZ吴彦祖 2024-06-17 4/200 2024-06-20 09:22 by 洗铁路的阿姨
[基金申请] 2024国社科通讯评审 +21 qsd10086 2024-06-13 44/2200 2024-06-20 09:22 by 全球漂
[基金申请] 刚刚收到科研之友邮件 +14 olivermiaoer 2024-06-19 18/900 2024-06-20 08:33 by olivermiaoer
[找工作] 药学硕士找不到工作,打算去做科研助理了 +9 pom戴墨镜 2024-06-14 22/1100 2024-06-20 05:47 by 投必得科研顾问
[考博] 2025年博士申请,可先做科研助理 +3 limit888 2024-06-18 6/300 2024-06-20 02:43 by limit888
[文学芳草园] 累并快乐着 +20 MYHLD521 2024-06-14 20/1000 2024-06-19 23:04 by tentoone
[基金申请] 江南大学到瑞士招聘,称取消非升即走,改预聘+长聘 +21 babu2015 2024-06-18 22/1100 2024-06-19 23:03 by feng6531
[基金申请] 青年基金会评专家到底是怎么会评的呀?主审专家是不是一般不会改动系统按函评给的顺序 5+3 他山攻玉之石 2024-06-18 14/700 2024-06-19 19:43 by liliwang215
[考博] 34岁读博士晚吗 +39 emitdne 2024-06-13 39/1950 2024-06-19 19:03 by Dqtjs
[论文投稿] 审稿 +5 香瓜木香 2024-06-19 6/300 2024-06-19 17:44 by xli1984
[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀? +7 65syn 2024-06-14 11/550 2024-06-19 10:19 by 65syn
[考博] 2025考博 +8 自强不息a?a 2024-06-15 13/650 2024-06-18 18:12 by 投必得科研顾问
[论文投稿] 投稿求助 10+3 2022_灵魂工程师 2024-06-16 8/400 2024-06-18 12:12 by 投必得科研顾问
[硕博家园] 博士毕业高校和就业的相关问题 +7 SCITOPPP 2024-06-14 11/550 2024-06-18 07:51 by yinxing1995
[论文投稿] 论文已接收,但发现修改稿传的是旧版该怎么办? +3 zsq54321 2024-06-15 6/300 2024-06-17 19:45 by wjykycg
[基金申请] 面青地会评时间 +8 tanjydd 2024-06-15 8/400 2024-06-17 17:08 by 小龙虾2008
[基金申请] 博士后创新人才支持计划公示 +9 aishida144 2024-06-14 15/750 2024-06-16 09:52 by msjy
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
信息提示
请填处理意见