参考文献：

考虑柔性负荷的综合能源系统低碳经济优化调度_薛开阳

考虑用户侧柔性负荷的社区综合能源系统日前优化调度_刘蓉晖

主要内容：

基础模型参考刘蓉晖的论文，主要做了场景1、2、3；碳交易模型采用薛开阳论文中的。

采用CPIEX求解某微网的运行优化情况，下层优化得出的微网向配电网购电或售电功率，以及各机组的出力.基于能源集线器概念，结合需求侧柔性负荷的可平移、可转移、可削减特性，构建了含风光储、燃气轮机、柔性负荷等在内的 IES 模型。 综合考虑了系统运行成本和碳交易成本，建立了以总成本最低为优化目标的 IES 低碳经济调度模型，采用cplex求解器对算例进行求解。

部分代码：

clc;clear;close all;

%读取数据

%电负荷、热负荷、光伏、风机、购电价、售电价

e_load=[160150140140130135150180215250275320335290260275270280320360345310220160];%电负荷

h_load=[135140 150 135 140 120 115 100 115 115 160 180 190 170 140 130 145 200 220 230 160 150 140 130];%热负荷

ppv=[0 0 00 010 15 25 45 75 90 100 80 100 50 40 30 15 10 0 0 0 0 0 ];%光伏预测数据

pwt=[60 65 70 75 80 85 90 100 125 150 130 110 100 120 125 130 140 160 180 200 175 160 155 150];%风机预测数据

buy_price=[0.250.25 0.25 0.25 0.25 0.25 0.25 0.53 0.53 0.53 0.82 0.82 0.82 0.82 0.8 0.53 0.53 0.53 0.82 0.82 0.82 0.53 0.53 0.53];%购电价

sell_price=[0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.42 0.42 0.42 0.65 0.65 0.65 0.65 0.65 0.42 0.42 0.42 0.65 0.65 0.65 0.42 0.42 0.42];%售电价

%需求响应数据

Pcut=[10 10 10 10 10 10 15 15 25 50 50 50 50 50 50 50 50 50 50 50 40 40 15 10];%可削减电负荷

Temp_Pcut=binvar(1,24,'full'); % 电负荷削减标志

PPcut=sdpvar(1,24,'full');%电负荷消减量

n1=zeros(1,1);%消减连续

Hcut=[25 25 25 25 25 25 25 25 30 40 40 40 40 40 40 40 40 40 50 50 30 30 20 15];%可削减热负荷

Temp_Hcut=binvar(1,24,'full'); % 热负荷削减标志

HHcut=sdpvar(1,24,'full');%热负荷消减量

n2=zeros(1,1);%消减连续

Ptran=[0 0 0 0 0 0 0 0 0 0 0 0 25 25 25 25 0 0 0 0 0 0 0 0 ];%可转移电负荷

Temp_Ptran=binvar(1,24,'full'); % 可转移电负荷 转移标志

PPtran=sdpvar(1,24,'full');%电负荷转移量

Pshift1=[0 0 0 0 0 0 0 0 0 0 0 25 25 0 0 0 0 0 0 0 0 0 0 0 ];%可平移电负荷1

Temp_Pshift1=binvar(1,24,'full'); % 可平移电负荷1 平移标志

PPshift1=sdpvar(1,24,'full');%可平移电负荷1量

Pshift2=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 25 25 25 0 0 ];%可平移电负荷2

Temp_Pshift2=binvar(1,24,'full'); % 可平移电负荷2 平移标志

PPshift2=sdpvar(1,24,'full');%可平移电负荷2量

Hshift=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 45 45 0 0 0 0 ];%可平移热负荷

Temp_Hshift=binvar(1,24,'full'); % 可平移热负荷 平移标志

HHshift=sdpvar(1,24,'full');%可平移热负荷量

场景1：

场景2：

场景3：