你这个题目我试着做了一遍,很遗憾的是没有找到解。不知道你的模型具体是如何的。下面我把我的程序给你。不知道是否帮到你了。
model:
sets:
a/1..4/:w,t,bb;
b/1..5/:q;
r(a,b):x,c;
endsets
data:
t=1 0.7 1.5 0.5;
q=0.9 0.8 1.2 1.6 1;
w=2 2.5 1.3 3;
bb=50 40 40 50;
enddata
max=0.9*@sum(r(i,j):100*x(i,j)*c(i,j)-bb(i)*(1-x(i,j)));
@for(b(j):@sum(a(i):x(i,j))=3);
@for(b(j):@sum(r(i,j):c(i,j)*t(i))/(@sum(r(i,j):c(i,j)))=q(j));
@for(a(i):@for(b(j):c(i,j)
end
可以去matlab论坛询问大师
http://www.ilovematlab.cn/forum.php