【调峰】储能辅助电力系统调峰的容量需求研究（Matlab代码实现）

%% 连续变量
% 决策变量：火电出力, 风电出力, 调峰不足,u储能放电，u储能充电 ，储能容量,
P_fire = sdpvar(1,24);       % 火电出力出力
P_wind = sdpvar(1,24); %风电出力
P_loss = sdpvar(1,24);  %调峰不足
%储能
P_ch = sdpvar(1,24);   % 蓄电池充电功率
P_dis = sdpvar(1,24);   % 蓄电池放电功率
Uchr = binvar(1,24);   % 电池充电状态，1表示充电
U_dis = binvar(1,24);   % 电池放电状态，1表示放电
W_Bt= sdpvar(1,24);     % 蓄电池的荷电状态 
H_SS= sdpvar(1);      % 蓄电池容量 
%% 约束
constraint=[];
for t=1:24 % (1) 电能平衡约束    
    constraint = [constraint, P_fire(t)+P_wind(t)+P_loss(t)+P_dis(t)-P_ch(t)-Pload(t)==0 ]; 
end
constraint=[constraint, 0 <= H_SS <= E_max];
for t=1:24 %  蓄电池不等式约束
    constraint=[constraint, 0<=P_ch(t)<=Uchr(t)*E_max];
    constraint=[constraint, 0<=P_dis(t)<=U_dis(t)*E_max];
    constraint=[constraint, Uchr(t)+U_dis(t)<=1];
    constraint=[constraint, H_SS*soc_min<=W_Bt(t)<=H_SS*soc_max];
end
for t=1:24 % 蓄电池等式约束
     if t==1
        constraint=[constraint, W_Bt(t)==H_SS*soc0+P_ch(t)*eta-P_dis(t)/eta];
    else
        constraint=[constraint, W_Bt(t)==W_Bt(t-1)+P_ch(t)*eta-P_dis(t)/eta];
    end
end
 constraint = [constraint, H_SS*soc0==W_Bt(24) ]; 
%% 火电上下限约束
for t=1:24 
    constraint=[constraint, pf_min<=P_fire(t)<=pf_max] ; 
end
%% 火电爬坡约束
for t=2:24 
    constraint=[constraint, -delta<=P_fire(t)-P_fire(t-1)<=delta] ; 
end
%% 风电上下限约束
for t=1:24 
    constraint=[constraint, 0 <= P_wind(t) <= Pwind(t)] ; 
end
%% 调峰不足约束
for t=1:24
 constraint = [constraint, 0<= P_loss(t) <=Pload]; 
end