数学建模社区-数学中国

标题: 求大神帮我看看这个程序问题出在哪 [打印本页]
作者: lyfsky    时间: 2017-3-16 22:34
标题: 求大神帮我看看这个程序问题出在哪
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)));


作者: liwenhui    时间: 2017-3-17 09:28
这样就可以了：
  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

  6.   Model Class:                                        LP

  8.   Total variables:                     12
  9.   Nonlinear variables:                  0
  10.   Integer variables:                    0

  12.   Total constraints:                   14
  13.   Nonlinear constraints:                0

  15.   Total nonzeros:                      68
  16.   Nonlinear nonzeros:                   0

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

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

作者: liwenhui    时间: 2017-3-17 09:31
好像解都是0。你再考虑一下你的模型吧。




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5