基于主从博弈的社区综合能源系统分布式协同优化运行策略(Matlab代码实现）

%燃气发电机、锅炉常数
ae=0.0013;
be=0.16;
ce=0;
ah=0.0005;
bh=0.11;
ch=0;
ce_ave=0.7;%平均电价约束
ch_ave=0.45;%平均热价约束
n_c=0.8;%热交换效率
n_ex=0.83;   %余热回收效率
n_ice=0.35;   %内燃机发电效率
%热储能
H_storage_max=1500; h_n=0.98;h_charge=0.98;h_discharge=1;%热储能容量/自损/充热/放热；
%电储能
E_storage_max=2000; e_n=1;e_charge=0.95;e_discharge=0.95;%电储能容量/自损/充电/放电;
bggin=1000;%%电储能
for i=1:24
    B(1,i)=bggin+Pcharge(1,i)*e_charge-Pdischarge(1,i); % 0.98为转换率
    bggin=B(1,i);
end
begin=1000;%%热储能
for i=1:24
    L(1,i)=begin*h_n+h_charge*Hti(1,i)-Hto(1,i);%%%热储能容量
    begin=L(1,i);
end
%约束条件
Constraints =[];
for i=1:24
    Constraints=[Constraints,200<=L(1,i)<=H_storage_max];
end
Constraints=[Constraints,L(1,24)>=800];
for i=1:24
    Constraints=[Constraints,0<=Hti(1,i)<=200*UHti(1,i)];
    Constraints=[Constraints,0<=Hto(1,i)<=150*UHto(1,i)];
end
for i=1:23
    Constraints=[Constraints,-300<=Hti(1,i+1)-Hto(1,i+1)-(Hti(1,i)-Hto(1,i))<=200];
end
for i=1:24
    Constraints=[Constraints,UHti(1,i)+UHto(1,i)<=1];
end
%% 电储能容量约束、充电约束、放电约束、状态约束、SOC约束
for i=1:24
    Constraints=[Constraints,0<=Pcharge(1,i)<=200*UPcharge(1,i)];
    Constraints=[Constraints,0<=Pdischarge(1,i)<=200*UPdischarge(1,i)];
end
%% 蓄电池爬坡约束
for i=1:24
    if  i>0&&i<24
        Constraints=[Constraints,-200<=Pcharge(1,i+1)-Pdischarge(1,i+1)-(Pcharge(1,i)-Pdischarge(1,i))<=200];
    elseif i==24
        Constraints=[Constraints,-200<=Pcharge(1,1)-Pdischarge(1,1)-(Pcharge(1,i)-Pdischarge(1,i))<=200];
    end
end
%% 蓄电池充放电约束
for i=1:24
    Constraints=[Constraints,UPcharge(1,i)+UPdischarge(1,i)<=1];
end
Constraints=[Constraints,sum(UPcharge(1,1:24)+UPdischarge(1,1:24))<=10];%考虑寿命
Constraints=[Constraints,B(1,24)==1000];
for i=1:24
    Constraints=[Constraints,400<=B(1,i)<=1600];
end