各位大大好
我是用CPLEX内建的函式库去写
如果要处理的问题是二元整数规划
% Maximize x1 + 2 x2 + 3 x3 + x4
% Subject to
% - x1 + x2 + x3 + 10 x4 <= 20
% x1 - 3 x2 + x3 <= 30
% x2 - 3.5x4 = 0
% Binary Integer
% x1 x2 x3 x4
那么只要这样打以下
f = [-1 -2 -3 -1]';
Aineq = [-1 1 1 10;
1 -3 1 0];
bineq = [20 30]';
Aeq = [0 1 0 -3.5];
beq = 0;
options = cplexoptimset;
options.Display = 'on';
[x, fval, exitflag, output] = cplexbilp (f, Aineq, bineq, Aeq, beq, ...
[ ], options);
fprintf ('\nSolution status = %s\n', output.cplexstatusstring);
fprintf ('Solution value = %d\n', fval);
disp ('Values = ');
disp (x');
这样就能正确跑出最佳解了
问题来了
但是如果今天模型的限制式长这样
http://imgur.com/Pz1X1VN
该怎么处理变量前面包了两个sigma呢?
另外ARJH可以自行给定
感谢好心的大大

ARJH要给定吧

喔对了忘了说 ARJH可以自行给定

那就是单纯的跑for-loop 把f Aineq Bineq Aeq Beq矩阵给建出来

所以是包两个for loop在前面这样吗?

应该是三四个不定，重点是要跑过所有可能不确定的话，自己先手算几个小的case看看

可以请您有空的话示范Aeq跟beq的方法吗?一条限制式就好.感谢您

x2 - 3.5x4 = 0Aeq = [0 1 0 -3.5];beq = 0; ?

s大我的意思是如何跑for-loop把Aeq Beq矩阵建出来

你有用过for-loop吗？

有 但是似乎不太正确 我是matlab新手

把你的想法写出来，还有可尝试先从小的case试起

以限制式(2)来说 我是这样for a =1:A,for r =1:R,Aeq = [x(raj)]; beq = 1;

你没有写出raj的定义，而且Aeq里要放的是x_raj的系数，不是x_raj