Question

LINGO问题

model:
sets:
jumin/1,2,3,4/;
shuiku/1..3/:S;
var/1,2/;
gongshui(shuiku,jumin):X,P,A;
yongshui(var,jumin):R;
endsets

data:
S=50,60,50;
P=160,130,220,170,
140,130,190,150,
190,200,230,0;
A=1,1,1,1,
1,1,1,1,
1,1,1,0;
R=30 70 10 10
50 70 20 40;
enddata

max=450*(@sum(gongshui:X))-@sum(gongshui:P*X);
@for(shuiku(i):@sum(jumin(j):X(i,j)*A(i,j))<=S(i));
@for(jumin(j):@sum(shuiku(i):X(i,j)*A(i,j))>=R(1,j);
@sum(shuiku(i):X(i,j)*A(i,j))<=R(1,j)+R(2,j););
end

Answer 1

你把矩阵相乘的原理搞错了，一个3*4的矩阵无法乘以一个3*4的矩阵，只能乘以一个4*3的矩阵，你把矩阵x定义为4*3的就行了。

model:
sets:
jumin/1..4/;
shuiku/1..3/:S;
var/1,2/;
gongshui(shuiku,jumin):A,P;
gg(jumin,shuiku):X;
kk(shuiku,shuiku):Q;
yongshui(var,jumin):R;
endsets
data:
S=50,60,50;
P=160,130,220,170,
140,130,190,150,
190,200,230,0;
A=1,1,1,1,
1,1,1,1,
1,1,1,0;
R=30 70 10 10
50 70 20 40;
@text()=@status();
enddata
max=450*(M-N);
M=@sum(gg:x);
@for(kk(i,j):q(i,j)=@sum(jumin(k):P(i,k)*X(k,j)));
N=@sum(kk:q);
@for(shuiku(i):@sum(jumin(j):X(j,i)*A(i,j))<=S(i));
@for(jumin(j):@sum(shuiku(i):X(j,i)*A(i,j))>=R(1,j);
@sum(shuiku(i):X(j,i)*A(i,j))<=R(1,j)+R(2,j););
end

Answer 2

意思是该问题为无界解，你看看是不是啥地方符号搞错了或者漏了约束条件

Answer 3

无界解 肯定是你的约束有问题 你的约束X(4,4)可以取无穷大 自己看看题目跟你写的程序吧