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

P49 经济调度模型在确定节点电价中的应用

例3-9

cons=[];
% 行：机组阶段出力，列：时段
x=sdpvar(1,11);
%% 目标函数：购电成本最小：各机组的购电成本之和
% 目标函数:社会福利最大化
y=(41*x(1)+43*x(2)+44*x(3)+46*x(4)+47*x(5))-...
    (48*x(6)+47*x(7)+46*x(8)+45*x(9)+44*x(10)+41*x(11));
%% 约束条件

    cons=[cons,sum(x(1:5))==10+sum(x(6:11))];
P=[0 25
   0 20
   0 15
   0 20
   0 30
   0 17
   0 13
   0 10
   0 22
   0 18
   0 20];
   
 for i=1:11
    cons=[cons,P(i,1)<=x(i)<=P(i,2)];
 end
    
ops=sdpsettings('solver','gurobi','showprogress',1);
z=optimize(cons,y,ops);
x=double(x)
y=double(y)
if z.problem==0
    -value(y)
else
    disp('something wrong!');
end