%% 决策变量
P_GT_e=sdpvar(1,24); %GT的输出电功率
P_GT_h=sdpvar(1,24); %GT的输出热功率
P_GT_g=sdpvar(1,24); %GT消耗天然气功率
P_h_WHB=sdpvar(1,24); %GT供给WHB的热功率
P_h_ORC=sdpvar(1,24); %GT供给ORC的热功率
P_WHB_h=sdpvar(1,24); %WHB输出热功率
P_ORC_e=sdpvar(1,24); %ORC输出电功率
P_CHP_e=sdpvar(1,24); %chp的输出电功率
P_CHP_h=sdpvar(1,24); %chp的输出热功率
P_DG=sdpvar(1,24); %风电消纳功率
P_g_GB=sdpvar(1,24); %输入GB设备的天然气功率
P_GB_h=sdpvar(1,24); %GB设备输出的热功率
P_e_P2G=sdpvar(1,24); %P2G耗电功率
P_P2G_g=sdpvar(1,24); %P2G输出的天然气功率
%电热气储能
P_ES1_cha=sdpvar(1,24);P_ES2_cha=sdpvar(1,24);P_ES3_cha=sdpvar(1,24); %充放功率
P_ES1_dis=sdpvar(1,24);P_ES2_dis=sdpvar(1,24);P_ES3_dis=sdpvar(1,24);
S_1=sdpvar(1,24);S_2=sdpvar(1,24);S_3=sdpvar(1,24); %各储能的实时容量状态
B_ES1_cha=binvar(1,24);B_ES2_cha=binvar(1,24);B_ES3_cha=binvar(1,24); %充标志
B_ES1_dis=binvar(1,24);B_ES2_dis=binvar(1,24);B_ES3_dis=binvar(1,24); %放标志
P_e_buy=sdpvar(1,24); %购电功率
P_g_buy=sdpvar(1,24); %购气功率
%需求响应负荷
P_e_load=sdpvar(1,24); %经过需求响应后实际的电负荷
P_e_tran=sdpvar(1,24);  %参与DR的可转移电负荷
P_e_tranin=sdpvar(1,24); %可转移电负荷的转入功率
P_e_tranout=sdpvar(1,24); %可转移电负荷的转出功率
Zeta_etranin=binvar(1,24);%可转移电负荷的转入功率二进制变量
Zeta_etranout=binvar(1,24);%可转移电负荷的转出功率二进制变量
P_e_re=sdpvar(1,24);  %参与DR的可替代电负荷
P_e_rein=sdpvar(1,24); %可替代电负荷的转入功率
P_e_reout=sdpvar(1,24); %可替代电负荷的转出功率
Zeta_erein=binvar(1,24);%可替代电负荷的转入功率二进制变量
Zeta_ereout=binvar(1,24);%可替代电负荷的转出功率二进制变量
P_h_load=sdpvar(1,24); %经过需求响应后实际的热负荷
P_h_tran=sdpvar(1,24);  %参与DR的可转移热负荷
P_h_tranin=sdpvar(1,24); %可转移热负荷的转入功率
P_h_tranout=sdpvar(1,24); %可转移热负荷的转出功率
Zeta_htranin=binvar(1,24);%可转移热负荷的转入功率二进制变量
Zeta_htranout=binvar(1,24);%可转移热负荷的转出功率二进制变量
P_h_re=sdpvar(1,24);  %参与DR的可替代热负荷
P_h_rein=sdpvar(1,24); %可替代热负荷的转入功率
P_h_reout=sdpvar(1,24); %可替代热负荷的转出功率
Zeta_hrein=binvar(1,24);%可替代热负荷的转入功率二进制变量
Zeta_hreout=binvar(1,24);%可替代热负荷的转出功率二进制变量
P_g_load=sdpvar(1,24); %经过需求响应后实际的热负荷
P_g_tran=sdpvar(1,24);  %参与DR的可转移热负荷
P_g_tranin=sdpvar(1,24); %可转移热负荷的转入功率
P_g_tranout=sdpvar(1,24); %可转移热负荷的转出功率
Zeta_gtranin=binvar(1,24);%可转移热负荷的转入功率二进制变量
Zeta_gtranout=binvar(1,24);%可转移热负荷的转出功率二进制变量
P_g_re=sdpvar(1,24);  %参与DR的可替代热负荷
P_g_rein=sdpvar(1,24); %可替代热负荷的转入功率
P_g_reout=sdpvar(1,24); %可替代热负荷的转出功率
Zeta_grein=binvar(1,24);%可替代热负荷的转入功率二进制变量
Zeta_greout=binvar(1,24);%可替代热负荷的转出功率二进制变量