基于绿证-阶梯式碳交易交互的源荷互补调度优化

23年新鲜代码，基本完成四个场景的复现。

针对多能精合的区域综合能源系统的低经济运行问题，提出基于绿证-阶梯式碳交易交与的源荷互补优化调度模型。首先，通过引入绿证-阶梯式碳交易交互机制来提高源侧可再生能源的消纳水平和降低系统碳排放量，其次，在负荷侧引入考用户满意度的激励型需求响应和调峰收益来实现热电负荷的"峰填公”。最后以日运行成本最小化为标。

部分代码

clear;clc;close all;
tic
%%
%% 各类负荷与风光数据
PL=[1231,1152,1126,1146,1179,1388,1632,2165,2296,2357,2235,2004,1969,2001,2152,2484,2676,2889,2976,2898,2549,2104,1807,1528];%电负荷kW
HL=[533,551,560,569,640,1066,2069,2351,2250,2202,2167,1847,1616,1430,1359,1439,1563,1654,2297,2410,2235,963,933,684];%热负荷
WT=[407,423,449,355,300,348,530,630,611,654,952,1007,865,754,741,757,945,1075,1192,1147,1082,1024,848,657];%风
PV=[0,0,0,0,0,0,46,122,238,375,433,463,420,376,330,203,35,0,0,0,0,0,0,0];%光
%碳交易
%电量、热量碳排放配额
a1=1.08;%kg/kwh，单位发电量产生的co2。
a2=0.23;%单位发热产生的co2
a3=0.043;%kg/kwh,单位热量无偿排放量
a4=0.72;%无偿获得的单位排放量
%电价气价
pwt=0.3;
ppv=0.35;
pe=[0.32 0.32 0.32 0.32 0.32 0.61 0.61 0.61 0.92 0.92 0.92 0.61 0.61 0.61 0.61 0.61 0.92 0.92 1 1 1 0.61 0.61 0.61];%电价
pch4=[2.3 2.3 2.3 2.3 2.3 2.3 2.3 2.3 3 3 3 3 3 3 2.8 2.8 2.8 2.8 3.3 3.3 3.3 2.8 2.3 2.3,];%元/m3
%能源设备效率
K=9.5;%天然气热值
K0=0.58;%GF发电效率
K1=0.05;%GF余热回收 
K2=0.95;%GB效率
K3=0.45;%CHP气热
K4=0.3;%CHP气电
%% CHP热电联产
CHP=sdpvar(1,24);%燃气轮机气功率
PCHP=0.4*CHP;   %%电功率
HCHP=0.6*CHP;  %%热功率