QQ登录

只需要一步,快速开始

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

求大神指出哪里错了!3Q!

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

1

主题

5

听众

151

积分

升级  25.5%

  • TA的每日心情

    2013-10-5 17:26
  • 签到天数: 47 天

    [LV.5]常住居民I

    自我介绍
    工科生
    跳转到指定楼层
    1#
    发表于 2012-8-30 16:52 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    1. model:
    2. sets:l/1..10/:y,k,h;
    3. m/1..5/:a;
    4. link(l,m):x,z,d;
    5. endsets

    6. data:
    7. d=5.26 5.19 4.21 4.00 2.95 2.74 2.46 1.90 0.64 1.27 1.90 0.99 1.90 1.13 1.27 2.25 1.48 2.04 3.09 3.51 5.89 5.61 5.61 4.56 3.51 3.65 2.46 2.46 1.06 0.57 0.64 1.76 1.27 1.83 2.74 2.60 4.21 3.72 5.05 6.10 4.42 3.86 3.72 3.16 2.25 2.81 0.78 1.62 1.27 0.50;
    8. a=78 85 124 85 85;
    9. y=81 79 87 68 74 87 68 74 87 87;
    10. k=61 68 64 68 71 81 68 84 87 81;
    11. enddata

    12. min=@sum(link:x*z*d);
    13. @for(m(j):@sum(link(i,j):z(i,j)*x(i,j))>=a(j));
    14. @for(l(i):@sum(link(i,j)|j#eq#1 #or# j#eq#2:x(i,j))<=y(i));
    15. @for(l(i):@sum(link(i,j)|j#ge#3 #and# j#le#5:x(i,j))<=k(i));
    16. @for(m(j):@sum(link(i,j):x(i,j))<=96);
    17. @for(m(j):@sum(link(i,j):z(i,j)*x(i,j)*h(i))/@sum(link(i,j):z(i,j)*x(i,j))>=0.285);
    18. @for(m(j):@sum(link(i,j):z(i,j)*x(i,j)*h(i))/@sum(link(i,j):z(i,j)*x(i,j))<=0.305);
    19. @for(link:@bin(z));
    20. end
    复制代码
    程序找不到最优解

    5d6034a85edf8db183b91ffa0923dd54574e749c.jpg (348.41 KB, 下载次数: 185)

    5d6034a85edf8db183b91ffa0923dd54574e749c.jpg

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

    32

    主题

    50

    听众

    2942

    积分

  • TA的每日心情
    慵懒
    2024-7-11 11:07
  • 签到天数: 543 天

    [LV.9]以坛为家II

    社区QQ达人 邮箱绑定达人 发帖功臣 新人进步奖 风雨历程奖 最具活力勋章

    群组小草的客厅

    群组数学建模认证项目实训

    群组我们一定会赢

    群组哈尔滨工业大学建模团

    群组聊天

    回复

    使用道具 举报

    且生        

    29

    主题

    9

    听众

    1500

    积分

    升级  50%

  • TA的每日心情
    慵懒
    2016-9-24 15:19
  • 签到天数: 412 天

    [LV.9]以坛为家II

    社区QQ达人

    群组学术交流A

    群组学术交流B

    群组2013认证赛B题讨论群组

    群组EXCEL

    回复

    使用道具 举报

    0

    主题

    5

    听众

    45

    积分

    升级  42.11%

  • TA的每日心情
    郁闷
    2012-9-6 20:11
  • 签到天数: 5 天

    [LV.2]偶尔看看I

    群组学术交流B

    群组数学建模

    回复

    使用道具 举报

    小卢        

    6

    主题

    6

    听众

    218

    积分

    升级  59%

  • TA的每日心情
    开心
    2016-5-22 22:09
  • 签到天数: 86 天

    [LV.6]常住居民II

    社区QQ达人

    群组学术交流B

    群组学术交流A

    本帖最后由 wujianjack2 于 2013-7-20 13:59 编辑

    好像在程序中把z当做y了,我改了下你的程序,不过最后没有灵敏度解,这个我也不清楚是什么原因了。。。
    程序:
    model:
    sets:l/1..10/:y,k,h;
    m/1..5/:a;
    link(l,m):x,z,d;
    endsets

    data:
    d=5.26 5.19 4.21 4.00 2.95 2.74 2.46 1.90 0.64 1.27 1.90 0.99 1.90 1.13 1.27 2.25 1.48 2.04 3.09 3.51 5.89 5.61 5.61 4.56 3.51 3.65 2.46 2.46 1.06 0.57 0.64 1.76 1.27 1.83 2.74 2.60 4.21 3.72 5.05 6.10 4.42 3.86 3.72 3.16 2.25 2.81 0.78 1.62 1.27 0.50;
    a=78 85 124 85 85;
    y=81 79 87 68 74 87 68 74 87 87;
    k=61 68 64 68 71 81 68 84 87 81;
    enddata

    min=@sum(link:x*z*d);
    @sum(link(i,j)|j#eq#1:z(i,j)*x(i,j))>=78;
    @sum(link(i,j)|j#ne#1#and#j#ne#3:z(i,j)*x(i,j))>=85;
    @sum(link(i,j)|j#eq#3:z(i,j)*x(i,j))>=124;
    @for(l(i):@sum(link(i,j)|j#eq#1 #and# j#eq#2:x(i,j))<=y(i));
    @for(l(i):@sum(link(i,j)|j#ge#3 #and# j#le#5:x(i,j))<=k(i));
    @for(m(j):@sum(l(i):x(i,j))<=96);
    @for(m(j)|j#ne#3#and#j#ne#4:(@sum(link(i,j):z(i,j)*x(i,j)*h(i))/@sum(link(i,j):z(i,j)*x(i,j)))>=0.285);
    @for(m(j)|j#ne#3#and#j#ne#4:(@sum(link(i,j):z(i,j)*x(i,j)*h(i))/@sum(link(i,j):z(i,j)*x(i,j)))<=0.305);
    @for(link:@bin(z));
    end
    回复

    使用道具 举报

    Smallvan        

    0

    主题

    3

    听众

    258

    积分

    升级  79%

  • TA的每日心情
    慵懒
    2013-4-7 12:39
  • 签到天数: 59 天

    [LV.5]常住居民I

    群组学术交流A

    回复

    使用道具 举报

    26

    主题

    10

    听众

    562

    积分

    升级  87.33%

  • TA的每日心情
    开心
    2015-5-2 15:22
  • 签到天数: 174 天

    [LV.7]常住居民III

    自我介绍
    自信

    社区QQ达人

    群组学术交流B

    群组学术交流A

    回复

    使用道具 举报

    26

    主题

    10

    听众

    562

    积分

    升级  87.33%

  • TA的每日心情
    开心
    2015-5-2 15:22
  • 签到天数: 174 天

    [LV.7]常住居民III

    自我介绍
    自信

    社区QQ达人

    群组学术交流B

    群组学术交流A

    还有那个x(i,j)上面有一横的那个是指,x的平均数吗?感觉你写代码的时候好想没有体现平均数吧
    回复

    使用道具 举报

    hzs2012        

    3

    主题

    6

    听众

    152

    积分

    升级  26%

  • TA的每日心情
    开心
    2013-2-4 14:45
  • 签到天数: 32 天

    [LV.5]常住居民I

    自我介绍
    hello world

    群组Matlab讨论组

    群组学术交流B

    回复

    使用道具 举报

    5

    主题

    3

    听众

    244

    积分

    升级  72%

  • TA的每日心情
    开心
    2013-12-30 15:23
  • 签到天数: 62 天

    [LV.6]常住居民II

    新人进步奖

    群组学术交流C

    群组数学之家

    群组学术交流D

    群组学术交流A

    群组学术交流B

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2025-7-7 05:43 , Processed in 0.760098 second(s), 106 queries .

    回顶部