电力调度优化理论与应用例题代码（5）

P69 混合整数线性规划算法

例4-5

特别注意，表达式应该都是的格式，才能调用gurobi

cons=[];
% 行：机组阶段出力，列：时段
P=sdpvar(2,3);
F_G=sdpvar(2,3);
SU=sdpvar(2,3);

I=binvar(2,4);
open=binvar(2,3);
off=binvar(2,3);
%% 目标函数：购电成本最小：各机组的购电成本之和
for t=1:3
    F_G(1,t)=100*I(1,t+1)+8*P(1,t);
    F_G(2,t)=100*I(2,t+1)+8*P(2,t);
    SU(1,t)=200*open(1,t);
    SU(2,t)=200*open(2,t);
end
    
y=sum(F_G(:))+sum(SU(:));
%% 约束条件

% I(1,1)=0;
% I(2,1)=0;


Pd=[170
    240
    190];
   
 for t=1:3 

%  功率平衡约束
    cons=[cons,P(1,t)+P(2,t)==Pd(t)];
%  机组出力约束
    cons=[cons,50*I(1,t+1)<=P(1,t)<=200*I(1,t+1)];
    cons=[cons,20*I(2,t+1)<=P(2,t)<=100*I(2,t+1)];
  
%  启停状态的约束
    cons=[cons,