参考文献：

考虑用户舒适度的冷热电多能互补综合能源系统优化调度-邹云阳

主要内容：

针对综合能源系统优化调度的研究，具体包括：

1.建立含风光、P2G、燃气轮机、燃气锅炉等多能耦合元件的运行特性模型

2.电、热、冷、气多能稳态能流模型

3.冷、热负荷考虑热惯性得出，电、气负荷和风电光伏预测得出。

4.并网模式下考虑经济成本最优、碳排放最优的优化调度模型

部分程序：

clc;clear;close all;% 程序初始化

%% 读取数据

shuju=xlsread('shuju数据.xlsx'); %把一天划分为24小时

load_e=shuju(2,:); %电负荷

load_g=shuju(3,:); %气负荷

P_PV=shuju(4,:); %光电出力预期值

P_WT=shuju(5,:); %风电出力预期值

T_out=shuju(6,:); %室外温度

%% 各变量及常量定义

P_G3=sdpvar(1,24,'full');%微燃气轮机电功率出力

e_G3=0.26;%燃气轮机电效率

h_G3=0.68;%燃气轮机热效率

P_EH=sdpvar(1,24,'full');%余热锅炉输出热功率

EH=0.6;%余热回收效率

P_GH=sdpvar(1,24,'full');%燃气锅炉输出热功率

GH=0.85;%燃气锅炉效率

P_AC=sdpvar(1,24,'full');%吸收式制冷机输出冷功率

COP_AC=0.8;%吸收式制冷机制冷系数

P_EC=sdpvar(1,24,'full');%电制冷机输出冷功率

COP_EC=3;%电制冷机冷系数

P_EG=sdpvar(1,24,'full');%P2G设备输出气功率

EG=0.6;%P2G设备综合转换效率

Pbuy=sdpvar(1,24,'full');%从电网购电电量

Psell=sdpvar(1,24,'full');%向电网售电电量

Pnet=sdpvar(1,24,'full');%交换功率

Temp_net=binvar(1,24,'full'); % 购|售电标志

Gbuy=sdpvar(1,24,'full');%从气网购气量

R=0.93;cc=0.54; %房间热阻和热容，这里分析改变等效热阻Ｒ时对供冷系统的影响

load_h=sdpvar(1,24);

T_hui=sdpvar(1,24);

并网模式下，经济成本最优调度模型：

并网模式下，碳排放最优调度模型：