主要内容
该程序实现一个综合能源系统的优化调度双层模型，上下层分别采用差分进化算法和规划算法进行求解。模型考虑了多种能源设备，包括燃气轮机、燃气锅炉、风电、光伏、储能设备等的协同运行，同时兼顾了能源供应商和用户的利益，需满足各种约束条件，包括设备出力约束、储能约束、负荷平衡约束等，以供应商和用户的收益、成本等为目标进行求解，同时考虑到负荷需求响应，通过分层模型将非线性求解模型转化为线性求解模型，实现了系统的经济和性能优化。程序采用matlab+cplex进行求解，注释清楚，有对应的文档说明，方便学习研究！
内容研究
1.模型简介
上层模型目标函数为：
为了便于理解，分块分析如下：
I3代表运营商售电收入，I4代表运营商售热收入，C_CCHP代表CCHP的运行成本，C_grid代表购电成本，FD代表储能的售电收益，FH代表售热能收益。
下层目标函数为：
逐项分析如下：
FU代表用户满意度，F5代表用户的购电，F6代表用户购热成本。
该目标的意思是用户满意度最高，用户购电和购热成本最低。
下层模型是采用规划算法结合cplex优化主体出力结果和目标值。
涉及到的约束主要有：
- 设备出力约束
- 设备爬坡约束
- 热电功率平衡约束
- 与电网交互约束
- 用户侧热电储能约束
- 用户侧负荷转移和削减约束
2 程序释义
程序包括的子程序较多，相互间的调用关系可参考下图：
具体每部分代码含义详见下载文件夹中说明文档。
部分代码
%基础热负荷
dh=[1000,1010,1043,1155,1172,1197,1158,1060,923,880,910,830,785,730,700,730,810,820,810,910,950,1010,1030,1020];
%基础电负荷
de=[540,528,504,516,524,520,612,652,696,856,956,968,992,864,720,676,680,668,708,976,968,960,572,532];
de1=[408,404,372,392,396,388,448,448,460,524,620,780,764,768,612,628,620,632,660,868,864,848,448,404];
%风电出力
Pwt=[203,277,264,331,137,81,72,141,43,12,20,12,5,48,86,346,287,530,491,448,603,601,403,380];
%光伏出力
Ppv=[0,0,0,0,0,0,97,220,336,410,486,444,453,445,442,325,202,140,29,0,0,0,0,0];
%电网分时电价、上网电价
grid_fs=[0.4,0.4,0.4,0.4,0.4,0.4,0.4,0.4,0.4,0.8,0.8,0.8,1.2,1.2,1.2,1.2,0.8,0.8,0.8,0.8,1.2,1.2,1.2,0.8];
grid_sw=[0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35];
ch_min=[0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15];%热价下限
ch_max=[0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5];%热价上限
Pice=sdpvar(1,24,'full');%电输出功率
Qgb=sdpvar(1,24,'full');%热输出功率
Pbuy=sdpvar(1,24,'full');%从电网购电电量
Psell=sdpvar(1,24,'full');%向电网售电电量
Pnet=sdpvar(1,24,'full');%交换功率
Temp_net=binvar(1,24,'full'); % 购|售电标志
le=sdpvar(1,24,'full');%可转移电负荷
lh=sdpvar(1,24,'full');%可转移热负荷
Pcharge=sdpvar(1,24,'full');UPcharge=binvar(1,24,'full');% 充电
Pdischarge=sdpvar(1,24,'full');UPdischarge=binvar(1,24,'full');% 放电
Hti=sdpvar(1,24,'full');UHti=binvar(1,24,'full');%充热
Hto=sdpvar(1,24,'full');UHto=binvar(1,24,'full');%放热
%用户偏好常系数
ve=1.5;
a_e=0.0009;
vh=1.1;
a_h=0.0011;
%燃气发电机、锅炉常数
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;%%电储能
运行结果
