多园区综合能源系统优化调度+碳排放和绿色证书联合交易

clc;
clear;
mydata=xlsread('data7.xlsx');
data1=mydata(1:24,3:8);
data2=mydata(27:50,3:8); 
data3=mydata(53:76,3:8);
data4=mydata(79:102,3:8);
data(:,:,1)=data1;
data(:,:,2)=data2;
data(:,:,3)=data3;

paras(:,1,1:4)=mydata(105:134,2:5);
clear data1 data2 data3 data4;
%系统参数

%% 燃气轮机参数
%CCHP1
MT_max1=paras(2,1,1);
MT_min1=paras(3,1,1);
eta_MT1=paras(4,1,1);
alpha_MT1=paras(5,1,1);

%CCHP2
MT_max2=paras(2,1,2);
MT_min2=paras(3,1,2);
eta_MT2=paras(4,1,2);
alpha_MT2=paras(5,1,2);

%CCHP3
MT_max3=paras(2,1,3);
MT_min3=paras(3,1,3);
eta_MT3=paras(4,1,3);
alpha_MT3=paras(5,1,3);

%CCHP4
MT_max4=paras(2,1,4);
MT_min4=paras(3,1,4);
eta_MT4=paras(4,1,4);
alpha_MT4=paras(5,1,4);
%% 冷热电、光伏，购售电价，天然气价格
price_E_sell=0.606; 
%CCHP1
price_E_buy1=data(:,5,1);
price_G1=paras(1,1,1);
load_E1=data(:,1,1);
load_H1=data(:,2,1);
load_C1=data(:,3,1);
P_PV1=data(:,4,1);
%CCHP2
price_E_buy2=data(:,5,2);
price_G2=paras(1,1,2);
load_E2=data(:,1,2);
load_H2=data(:,2,2);
load_C2=data(:,3,2);
P_PV2=data(:,4,2);
%CCHP3
price_E_buy3=data(:,5,3);
price_G3=paras(1,1,3);
load_E3=data(:,1,3);
load_H3=data(:,2,3);
load_C3=data(:,3,3);
P_PV3=data(:,4,3);

%% 燃气锅炉参数
%CCHP1
GB_max1=paras(9,1,1);
GB_min1=paras(10,1,1);
eta_GB1=paras(11,1,1);
%CCHP2
GB_max2=paras(9,1,2);
GB_min2=paras(10,1,2);
eta_GB2=paras(11,1,2);
%CCHP3
GB_max3=paras(9,1,3);
GB_min3=paras(10,1,3);
eta_GB3=paras(11,1,3);

%% 吸附式制冷机、电制冷机运行参数
%CCHP1
EC_max1=paras(12,1,1);
EC_min1=paras(13,1,1);
eta_EC1=paras(14,1,1);
AC_max1=paras(15,1,1);
AC_min1=paras(16,1,1);
eta_AC1=paras(17,1,1);
%CCHP2
EC_max2=paras(12,1,2);
EC_min2=paras(13,1,2);
eta_EC2=paras(14,1,2);
AC_max2=paras(15,1,2);
AC_min2=paras(16,1,2);
eta_AC2=paras(17,1,2);
%CCHP3
EC_max3=paras(12,1,3);
EC_min3=paras(13,1,3);
eta_EC3=paras(14,1,3);
AC_max3=paras(15,1,3);
AC_min3=paras(16,1,3);
eta_AC3=paras(17,1,3);

%% 余热锅炉、余热回收装置参数
%CCHP1
Rcv_max1=paras(6,1,1);
Rcv_min1=paras(7,1,1);
eta_Rcv1=paras(8,1,1);
HE_max1=paras(18,1,1);
HE_min1=paras(19,1,1);
eta_HE1=paras(20,1,1);
%CCHP2
Rcv_max2=paras(6,1,2);
Rcv_min2=paras(7,1,2);
eta_Rcv2=paras(8,1,2);
HE_max2=paras(18,1,2);
HE_min2=paras(19,1,2);
eta_HE2=paras(20,1,2);
%CCHP3
Rcv_max3=paras(6,1,3);
Rcv_min3=paras(8,1,3);
eta_Rcv3=paras(8,1,3);
HE_max3=paras(18,1,3);
HE_min3=paras(19,1,3);
eta_HE3=paras(20,1,3);

%% 蓄电池参数
%CCHP1
BT_charge_max1=paras(21,1,1);
BT_discharge_max1=paras(22,1,1);
eta_charge1=paras(23,1,1);
eta_discharge1=paras(24,1,1);
Cap_BT1=paras(25,1,1);
BT_miu1=paras(26,1,1);
BT_loss1=paras(27,1,1);
BT_W0_1=paras(28,1,1);
BT_Wmax1=paras(29,1,1);
BT_Wmin1=paras(30,1,1);
%CCHP2
BT_charge_max2=paras(21,1,2);
BT_discharge_max2=paras(22,1,2);
eta_charge2=paras(23,1,2);
eta_discharge2=paras(24,1,2);
Cap_BT2=paras(25,1,2);
BT_miu2=paras(26,1,2);
BT_loss2=paras(27,1,2);
BT_W0_2=paras(28,1,2);
BT_Wmax2=paras(29,1,2);
BT_Wmin2=paras(30,1,2);
%CCHP3
BT_charge_max3=paras(21,1,3);
BT_discharge_max3=paras(22,1,3);
eta_charge3=paras(23,1,3);
eta_discharge3=paras(24,1,3);
Cap_BT3=paras(25,1,3);
BT_miu3=paras(26,1,3);
BT_loss3=paras(27,1,3);
BT_W0_3=paras(28,1,3);
BT_Wmax3=paras(29,1,3);
BT_Wmin3=paras(30,1,3);
%% 风机参数
P_wind1=data(:,6,1);
P_wind2=data(:,6,2);
P_wind3=data(:,6,3);
 

MATLAB代码实现-多园区综合能源系统优化调度+碳排放和绿色证书联合交易