含电热联合系统的微电网运行优化matlab程序（yalmip+cplex）（yalmip+gurobi）
参考文献：含电热联合系统的微电网运行优化

在当前能源互联网迅速发展及电热联系日渐紧密的环境下，提出基于电热联合调度的区域并网型微电网运行优化模型。综合网内储能特性、分时电价、电热负荷与分布式电源的时序特征，以包含风机、光伏电池、热电联产系统、电锅炉、燃料电池和储能系统的并网型微电网为例，采用Cplex
优化软件求得调度周期内各微电源最佳出力及总运行成本，并与两种常见电热调度方式进行比较。仿真算例表明：联合调度模型能实现电热统一协调调度并降低微电网运行成本。该模型可为电热之间能源互联及规划运营提供参考。

1 微电网内常见单元的数学模型
1.1 典型微电网结构

1.2 热电联产系统

1.3 电锅炉

1.4 储能系统

1.5 燃料电池

2 微电网经济性模型及求解方法
2.1 目标函数

2.2 约束条件
1）能量平衡约束。

2）联络线功率约束。

3）储能的相关约束。

4）可控机组运行条件。

2.3 求解算法

3 matlab程序

// 含电热联合系统的微电网运行优化_李正茂
%% 含电热联合系统的微电网运行优化_李正茂
% 采用Cplex优化软件求得调度周期内各微电源最佳出力及总运行成本clc; 
clear; 
close all; 
warning off;%% 定义决策变量
% 输入参数
parameter_cplex;
Ppv_1 = sdpvar(1,24);
Pwt_1 = sdpvar(1,24);
% 电设备
Pmt = sdpvar(1,24);            % 燃气轮机电功率
Peb = sdpvar(1,24);            % 电锅炉电功率 
Pfc = sdpvar(1,24);              % 燃料电池电功率
Pgrid = sdpvar(1,24);             % 电网功率                   >0 买电，<0 卖电
Eb= sdpvar(1,24);               % 电储能的容量
Pbch = sdpvar(1,24);          % 电储能充电功率
Pbdis = sdpvar(1,24);         % 电储能放电功率
% 热设备
Pmth = sdpvar(1,24);          % 燃气轮机热功率
Pheb = sdpvar(1,24);          % 电锅炉热功率
Eh= sdpvar(1,24);               % 热储能的容量
Phch = sdpvar(1,24);          % 热储能储热功率
Phdis = sdpvar(1,24);         % 热储能放热功率
% 辅助变量
Ubch = binvar(1,24);   % 电池充电状态，1表示充电
Ubdis = binvar(1,24);   % 电池放电状态，1表示放电
Uhch = binvar(1,24);   % 热储能储热状态，1表示储热
Uhdis = binvar(1,24);   % 热储能放热状态，1表示放热
% 初始化
objective=0;
constraint=[];
.............略

4 程序结果
1）电平衡

2）热平衡