说明书
免费：https://download.csdn.net/download/qq_50594161/87625438
MATLAB代码：考虑用户舒适度的冷热电多能互补综合能源系统优化调度
关键词：用户舒适度 综合能源 PMV 优化调度
参考文档：《冷热电气多能互补的微能源网鲁棒优化调度》基础模型加舒适度部分模型；
仿真平台：MATLAB+yalmip+cplex
主要内容：代码主要做的是考虑用户舒适度的冷热电多能互补综合能源系统优化调度模型，在传统的冷热电联供型综合能源系统的基础上，进一步考虑了热惯性以及用户的舒适度，并用预测平均投票数PMV对用户的舒适度进行衡量，且通过改变PMV的数值，可以对比不同舒适度要求对于综合能源系统调度结果的影响。同时，代码还补充性的考虑了碳排放交易机制，并设置经济性最优以及碳排放最优两种对比场景，从而丰富算例，效果非常明显。
%并网模式下，经济成本最优调度模型
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);
%% 约束条件
C =[];
%% 热冷负荷
%对于热负荷，利用PMV值得到供热时室内温度
T_out_hot=sdpvar(1,26); %供热时室外温度
for i=3:26
C=[C,T_out_hot(1)==0;T_out_hot(2)==0];%-1、0时刻的供热室外温度
C=[C,T_out_hot(i)==T_out(i-2)];
end
T_in_hot=sdpvar(1,26); %供热时室内温度
T_gong=sdpvar(1,26); %供水温度
PMV=2;%PMV值，公式4,这里分析改变PMV指标时对供热系统的影响
for i=3:26
C=[C,T_in_hot(1)==25;T_in_hot(2)==25];%-1、0时刻的供热室内温度
C=[C,-PMV<=(0.303*exp(-0.036*70)+0.028)*(70-0-3.05*(10^(-3))*(5733-6.99*70-2000)-0.42*(70-58.15)-1.7*(10^(-5))*70*(5867-2000)-1.4*(10^(-3))*70*(34-T_in_hot(i)-3.96*(10^(-8))*1.15*((32+273)^4-(29.7+273)^4))-1.15*4.7*(32-T_in_hot(i)))<=PMV];%公式4
end
