含电热联合系统的微电网运行优化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）热平衡
