[討論]線性規劃建模問題

[ljuyentintho (小劉)]

各位大大好

我是用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可以自行給定

感謝好心的大大






--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.114.54.170
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1456381061.A.5BD.html
※ 編輯: ljuyentintho (140.114.54.170), 02/25/2016 14:19:08
推 sunev: ARJH要給定吧
→ ljuyentintho: 喔對了忘了說 ARJH可以自行給定
※ 編輯: ljuyentintho (140.114.54.170), 02/25/2016 14:24:55
推 sunev: 那就是單純的跑for-loop 把f Aineq Bineq Aeq Beq矩陣給建
→ sunev: 出來
→ ljuyentintho: 所以是包兩個for loop在前面這樣嗎?
推 sunev: 應該是三四個不定，重點是要跑過所有可能
→ sunev: 不確定的話，自己先手算幾個小的case看看
→ ljuyentintho: 可以請您有空的話示範Aeq跟beq的方法嗎?
→ ljuyentintho: 一條限制式就好.感謝您
推 sunev: x2      - 3.5x4  = 0
→ sunev:  Aeq   = [0  1  0  -3.5];
→ sunev:  beq   =  0;    ?
→ ljuyentintho: s大我的意思是如何跑for-loop把Aeq Beq矩陣建出來
推 sunev: 你有用過for-loop嗎？
→ ljuyentintho: 有 但是似乎不太正確 我是matlab新手
→ sunev: 把你的想法寫出來，還有可嘗試先從小的case試起
→ ljuyentintho: 以限制式(2)來說 我是這樣for a =1:A,for r =1:R,
→ ljuyentintho: Aeq   = [x(raj)]; beq   =  1;
推 sunev: 你沒有寫出raj的定義，而且Aeq裡要放的是x_raj的係數，不是
→ sunev: x_raj