数学建模社区-数学中国

标题: 新手求高手帮助 [打印本页]
作者: fxw1062    时间: 2014-7-18 12:02
标题: 新手求高手帮助
MODEL:
SETS:
Shijian/1..7/:a,b;!a表示小时数，b表示需求;
Jiqi/1..4/:c,d,e,f,g,h;!c表示机器可用数，d表示每台机器最小输出功率，
e表示每台机器最大输出功率，f表示固定成本，g表示边际成本，h表示启动成本;
LINKS(Shijian,Jiqi):m,n,k;  !m表示运行机器数，n表示机器使用功率，k表示新开启机器数;
ENDSETS

DATA:     !初始化;
a=6,3,3,2,4,4,2;
b=12000,32000,25000,36000,25000,30000,18000;
c=10,4,8,3;
d=750,1000,1200,1800;
e=1750,1500,2000,3500;
f=2250,1800,3750,4800;
g=2.7,2.2,1.8,3.8;
h=5000,1600,2400,1200;
ENDDATA

p1=@sum(Shijian(i):@sum(Jiqi(j):a(i)*m(i,j)*f(j)));!固定成本和;
p2=@sum(Shijian(i):@sum(Jiqi(j):(n(i,j)-d(j))*m(i,j)*a(i)*g(j)));!边际成本和;
p3=@sum(Jiqi(j):m(1,j)*h(j))+@sum(Jiqi(j):@sum(Shijian(i)|i#ge#2:
@if(m(i,j)#gt#m(i-1,j),(m(i,j)-m(i-1,j)*h(j)),0));!启动成本;
MIN=p1+p2+p3;!目标函数;

!约束条件;
@for(Shijian(i):@sum(Jiqi(j):0.8*m(i,j)*n(i,j))=b(i));!需求约束;
@for(links(i,j):@gin(m(i,j)));
@for(links(i,j):@gin(k(i,j)));
@for(Shijian(i):@for(Jiqi(j):m(i,j),=c(j)));
@for(links(i,j):@bnd(d(j),n(i,j),e(j)));
END
解决的是一个电力生产问题，为什么约束条件的for提示错误
作者: madio    时间: 2014-7-18 17:17
  1. MODEL:
  2. SETS:
  3. Shijian/1..7/:a,b;
  4. Jiqi/1..4/:c,d,e,f,g,h;
  5. LINKS(Shijian,Jiqi):m,n,k;  
  6. ENDSETS

  8. DATA:     
  9. a=6,3,3,2,4,4,2;
  10. b=12000,32000,25000,36000,25000,30000,18000;
  11. c=10,4,8,3;
  12. d=750,1000,1200,1800;
  13. e=1750,1500,2000,3500;
  14. f=2250,1800,3750,4800;
  15. g=2.7,2.2,1.8,3.8;
  16. h=5000,1600,2400,1200;
  17. ENDDATA
  18. MIN=p1+p2+p3;
  19. p1=@sum(Shijian(i):@sum(Jiqi(j):a(i)*m(i,j)*f(j)));
  20. p2=@sum(Shijian(i):@sum(Jiqi(j):(n(i,j)-d(j))*m(i,j)*a(i)*g(j)));
  21. p3=@sum(Jiqi(j):m(1,j)*h(j))+@sum(Jiqi(j):@sum(Shijian(i)|i#ge#2:
  22. @if(m(i,j)#gt#m(i-1,j),(m(i,j)-m(i-1,j)*h(j)),0)));


  25. @for(Shijian(i):@sum(Jiqi(j):0.8*m(i,j)*n(i,j))=b(i));
  26. @for(links(i,j):@gin(m(i,j)));
  27. @for(links(i,j):@gin(k(i,j)));
  28. @for(Shijian(i):@for(Jiqi(j):m(i,j)=c(j)));
  29. @for(links(i,j):@bnd(d(j),n(i,j),e(j)));
  30. END
复制代码
改了一下，代码是没有问题了，但是没有可行解，是不是哪个约束条件写的和模型不一样。
作者: fxw1062    时间: 2014-7-19 13:02
madio 发表于 2014-7-18 17:17
改了一下，代码是没有问题了，但是没有可行解，是不是哪个约束条件写的和模型不一样。

谢谢你啦  可能是约束出了点问题吧！我很想知道哪里出了问题，避免再犯错了。
作者: madio    时间: 2014-7-19 13:21
fxw1062 发表于 2014-7-19 13:02
谢谢你啦  可能是约束出了点问题吧！我很想知道哪里出了问题，避免再犯错了。

你对一下代码就知道了，是少了一个右括号
作者: fxw1062    时间: 2014-7-19 15:20
madio 发表于 2014-7-19 13:21
你对一下代码就知道了，是少了一个右括号

ok   谢啦  真的是新手，好多都不会弄，以后有问题还来求教！



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