QQ登录

只需要一步,快速开始

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

求大神帮我看看这个程序问题出在哪

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

11

主题

13

听众

118

积分

升级  9%

  • TA的每日心情
    开心
    2018-1-11 16:13
  • 签到天数: 20 天

    [LV.4]偶尔看看III

    国际赛参赛者

    自我介绍
    我叫刘焱芳,20岁,想好好学建模
    跳转到指定楼层
    1#
    发表于 2017-3-16 22:34 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    model:
    sets:
    house/1,2,3/:m,v;
    cargo/1,2,3,4/:p,q,c;
    link(house,cargo):x;
    endsets
    data:
    !每个货舱所能容纳的最大质量;
    m=10,16,8;
    !每个货舱所能容纳的最大体积;
    v=6800,8700,5300;
    !每种货物的质量;
    p=18,15,23,12;
    !每种货物的体积(m^3*t^-1);
    q=480,650,580,390;
    !每种货物的利润;
    c=3100,3800,3500,2850;
    enddata
    !目标;
    max=@sum(link(i,j):x(i,j)*c(j));
    !约束条件;
    @for(house(i):@sum(cargo(j):x(i,j)/m(i))=@sum(cargo(j):x(i+1,j)/m(i+1)));
    @for(house(i):@sum(cargo(j): x(i,j)<=m(i));
    @for(house(i):@sum(cargo(j):x(i,j)*q(j)<=v(i)));
    @for(cargo(j):@sum(house(i):x(i,j)<=p(j)));


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

    70

    主题

    65

    听众

    5195

    积分

    独孤求败

  • TA的每日心情
    擦汗
    2018-4-26 23:29
  • 签到天数: 1502 天

    [LV.Master]伴坛终老

    自我介绍
    紫薇软剑,三十岁前所用,误伤义士不祥,乃弃之深谷。 重剑无锋,大巧不工。四十岁前恃之横行天下。 四十岁后,不滞于物,草木竹石均可为剑。自此精修,渐进至无剑胜有剑之境。

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

    群组计量经济学之性

    群组LINGO

    这样就可以了:
    1. model:
    2. sets:
    3. house/1,2,3/:m,v;
    4. cargo/1,2,3,4/:p,q,c;
    5. link(house,cargo):x;
    6. endsets
    7. data:
    8. !每个货舱所能容纳的最大质量;
    9. m=10,16,8;
    10. !每个货舱所能容纳的最大体积;
    11. v=6800,8700,5300;
    12. !每种货物的质量;
    13. p=18,15,23,12;
    14. !每种货物的体积(m^3*t^-1);
    15. q=480,650,580,390;
    16. !每种货物的利润;
    17. c=3100,3800,3500,2850;
    18. enddata
    19. !目标;
    20. max=@sum(link(i,j):x(i,j)*c(j));
    21. !约束条件;
    22. @for(house(i):@sum(cargo(j):x(i,j)/m(i))=@sum(cargo(j)|i#le#2:x(i+1,j)/m(i+1)));
    23. @for(house(i):@sum(cargo(j): x(i,j))<=m(i));
    24. @for(house(i):@sum(cargo(j):x(i,j)*q(j))<=v(i));
    25. @for(cargo(j):@sum(house(i):x(i,j))<=p(j));
    26. end
    复制代码
    运行结果:
    1.   Global optimal solution found.
    2.   Objective value:                              0.000000
    3.   Infeasibilities:                              0.000000
    4.   Total solver iterations:                             0

    5.   Model Class:                                        LP

    6.   Total variables:                     12
    7.   Nonlinear variables:                  0
    8.   Integer variables:                    0

    9.   Total constraints:                   14
    10.   Nonlinear constraints:                0

    11.   Total nonzeros:                      68
    12.   Nonlinear nonzeros:                   0

    13.                                 Variable           Value        Reduced Cost
    14.                                    M( 1)        10.00000            0.000000
    15.                                    M( 2)        16.00000            0.000000
    16.                                    M( 3)        8.000000            0.000000
    17.                                    V( 1)        6800.000            0.000000
    18.                                    V( 2)        8700.000            0.000000
    19.                                    V( 3)        5300.000            0.000000
    20.                                    P( 1)        18.00000            0.000000
    21.                                    P( 2)        15.00000            0.000000
    22.                                    P( 3)        23.00000            0.000000
    23.                                    P( 4)        12.00000            0.000000
    24.                                    Q( 1)        480.0000            0.000000
    25.                                    Q( 2)        650.0000            0.000000
    26.                                    Q( 3)        580.0000            0.000000
    27.                                    Q( 4)        390.0000            0.000000
    28.                                    C( 1)        3100.000            0.000000
    29.                                    C( 2)        3800.000            0.000000
    30.                                    C( 3)        3500.000            0.000000
    31.                                    C( 4)        2850.000            0.000000
    32.                                 X( 1, 1)        0.000000            700.0000
    33.                                 X( 1, 2)        0.000000            0.000000
    34.                                 X( 1, 3)        0.000000            300.0000
    35.                                 X( 1, 4)        0.000000            950.0000
    36.                                 X( 2, 1)        0.000000            700.0000
    37.                                 X( 2, 2)        0.000000            0.000000
    38.                                 X( 2, 3)        0.000000            300.0000
    39.                                 X( 2, 4)        0.000000            950.0000
    40.                                 X( 3, 1)        0.000000            700.0000
    41.                                 X( 3, 2)        0.000000            0.000000
    42.                                 X( 3, 3)        0.000000            300.0000
    43.                                 X( 3, 4)        0.000000            950.0000

    44.                                      Row    Slack or Surplus      Dual Price
    45.                                        1        0.000000            1.000000
    46.                                        2        0.000000            38000.00
    47.                                        3        0.000000            98800.00
    48.                                        4        0.000000            129200.0
    49.                                        5        10.00000            0.000000
    50.                                        6        16.00000            0.000000
    51.                                        7        8.000000            0.000000
    52.                                        8        6800.000            0.000000
    53.                                        9        8700.000            0.000000
    54.                                       10        5300.000            0.000000
    55.                                       11        18.00000            0.000000
    56.                                       12        15.00000            0.000000
    57.                                       13        23.00000            0.000000
    58.                                       14        12.00000            0.000000
    复制代码

    四十岁后,不滞于物,草木竹石均可为剑。
    回复

    使用道具 举报

    liwenhui        

    70

    主题

    65

    听众

    5195

    积分

    独孤求败

  • TA的每日心情
    擦汗
    2018-4-26 23:29
  • 签到天数: 1502 天

    [LV.Master]伴坛终老

    自我介绍
    紫薇软剑,三十岁前所用,误伤义士不祥,乃弃之深谷。 重剑无锋,大巧不工。四十岁前恃之横行天下。 四十岁后,不滞于物,草木竹石均可为剑。自此精修,渐进至无剑胜有剑之境。

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

    群组计量经济学之性

    群组LINGO

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2024-5-29 04:18 , Processed in 0.319904 second(s), 61 queries .

    回顶部