简单的综合能源系统(IES)优化调度代码免费！含Matpower添加教程。

1.综合能源系统介绍

综合能源系统优化调度是一种综合考虑不同能源形式（电力、气、冷、热等）的能源系统，在实时运行和长期规划中进行优化决策的方法。系统的优化调度旨在提高能源利用效率、降低成本、减少碳排放，确保可靠性和可持续性。

以下是典型的综合能源系统优化调度框架及介绍：

能源输入：包括如风能、太阳能、传统电力、天然气等。

能源转换：系统中存在能源转换设备，如燃气轮机、燃气锅炉、制冷机等，它们将煤、天气等能源转换为电、热能。

储能技术：电池和蓄热槽等储能技术在综合能源系统中扮演重要角色。

需求侧管理：优化调度还需要考虑不同用户的能源需求，以便在不同用电峰谷期间实现最佳的能源供应。

环境和经济目标：综合能源系统优化调度通常还需要考虑环境可持续性和经济效益。它可以帮助降低碳排放，减少能源成本，并提高整体能源系统的效率。

介绍综合能源系统组成及优化目的后，本文推荐一份简单的综合能源系统优化调度代码，程序中主要包括电力平衡，供综合能源系统研究领域的伙伴参考学习！

%% 初始化
clc;
clear;
%% 各变量及常量定义
%------------------------变量定义-----------------------%

Pw=sdpvar(1,24,'full'); %风机出力
Ppv=sdpvar(1,24,'full');%光伏出力
Pbat=sdpvar(1,24,'full');%蓄电池出力
% Pde=Power(4,:);%柴油机组出力
% Lshift=zeros(1,24);
Pnet=sdpvar(1,24,'full');%交换功率
Pbuy=sdpvar(1,24,'full');%从电网购电电量
Psell=sdpvar(1,24,'full');%向电网售电电量
Temp_net=binvar(1,24,'full'); % 购/售电标志
Temp_cha=binvar(1,24,'full'); %充电标志
Temp_dis=binvar(1,24,'full'); %放电标志
Temp_static=binvar(1,24,'full'); %电池静置标志
Pcha=sdpvar(1,24);
Pdis=sdpvar(1,24);
Constraints = [];
%-------------------------常量定义-----------------------%

Load=[88.24 83.01 80.15 79.01 76.07 78.39  89.95 128.85 155.45 176.35 193.71 182.57 179.64 166.31 164.61 164.61 174.48 203.93 218.99 238.11 216.14 173.87 131.07 94.04];%风机预测出力

Pw=[66.9 68.2 71.9 72 78.8 94.8 114.3 145.1 155.5 142.1 115.9 127.1 141.8 145.6 145.3 150 206.9 225.5 236.1 210.8 198.6 177.9 147.2 58.7];%光伏预测出力

Ppv=[0 0 0 0 0.06 6.54 20.19 39.61 49.64 88.62 101.59 66.78 110.46 67.41 31.53 50.76 20.6 22.08 2.07 0 0 0 0 0];%分时电价

C_buy=[0.25  0.25 0.25 0.25 0.25 0.25 0.25 0.53 0.53 0.53 0.82 0.82 0.82 0.82 0.82 0.53 0.53 0.53 0.82 0.82 0.82 0.53 0.53 0.53];

C_sell=[0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.42 0.42 0.42  0.65  0.65 0.65 0.65 0.65 0.42 0.42 0.42 0.65 0.65 0.65 0.42 0.42 0.42];

%% 约束
for k = 1:24
  Constraints = [Constraints, -160<=Pnet(1,k)<=160,0<=Pbuy(1,k)<=160, -160<=Psell(1,k)<=0]; %主网功率交换约束
Constraints=[Constraints,Pnet(1,k)+Pw(1,k)+Ppv(1,k)==Load(1,k) +Pbat(1,k)];   %功率平衡约束
  Constraints = [Constraints, implies(Temp_net(1,k),[Pnet(1,k)>=0,Pbuy(1,k)==Pnet(1,k),Psell(1,k)==0])]; %购电情况约束
  Constraints = [Constraints, implies(1-Temp_net(1,k),[Pnet(1,k)<=0,Psell(1,k)==Pnet(1,k),Pbuy(1,k)==0])]; %售电情况约束
%----------------------蓄电池约束--------------------%

 sum_bat=zeros(1,24);
  Constraints = [Constraints, -40<=Pbat(1,k)<=40,0<=Pcha(1,k)<=40,-40<=Pdis(1,k)<=0];%电池充放电约束
  Constraints = [Constraints, implies(Temp_cha(1,k),[Pbat(1,k)>=0,Pcha(1,k)==Pbat(1,k),Pdis(1,k)==0])];%充电情况约束
  Constraints = [Constraints, implies(Temp_dis(1,k),[Pbat(1,k)<=0,Pdis(1,k)==Pbat(1,k),Pcha(1,k)==0])];%放电情况约束
  Constraints = [Constraints, implies(Temp_static(1,k),[Pbat(1,k)==0,Pdis(1,k)==0,Pcha(1,k)==0])];%静置情况约束
  Constraints = [Constraints,Temp_cha(1,k)+Temp_dis(1,k) +Temp_static(1,k)==1];
%    sum_bat(1,k+1)=sum_bat(1,k)+Pcha(1,k)+Pdis(1,k);%计算SOC
%   Constraints=[Constraints,-30<=sum(Pdis(1,1:k)+Pcha(1,1:k))<=165] ;%SOC约束，电池容量300kwh，初始S0C为0.4，0.3<=SOC<=0.95
end
  Constraints=[Constraints,sum(Pdis+Pcha)==0] ;%ST=S0，始末SOC相等约束

%% 目标函数
 F=0;
%-----------------总费用--------------------%
for k = 1:24
  F = F + 0.52*Pw(:,k)+0.72*Ppv(:,k)+C_buy(:,k)*Pbuy(:,k) +C_sell(:,k)*Psell(1,k)+0.2*Pdis(1,k);
end

ops=sdpsettings('solver', 'cplex');%参数指定程序用cplex求解器
optimize(Constraints,F,ops)
value(F)%费用
value(Pcha)
value(Pdis)
%% 画图
x=1:24;
PP=[Pbuy;-Pdis;Pw;Ppv];
PP_neg=[Psell;-Pcha];
figure
bar(PP','stack');
h=legend('交换功率','蓄电池','风机出力','光伏出力','Location','NorthWest');
set(h,'Orientation','horizon')
hold on
bar(PP_neg','stack');
plot(x,value(Load),'r','linewidth',2);
xlabel('时段');ylabel('功率/kW');
hold off
for k=1:24
s(k)=value(sum(Pdis(1,1:k)+Pcha(1,1:k)))/300+0.4;
soc(k+1)=s(k);
end
soc(1)=0.4;
xx=0:24;
figure
plot(xx,soc,'r');
xlabel('时段');ylabel('SOC值');
title('蓄电池SOC状态');

2. Matpower 工具包添加方法

接上条Yalmip+Cplex工具包的添加方法，Matpower步骤类似Yalmip，供您练手！

MATPOWER（MATLAB Power System Optimization Toolbox）用于进行电力系统的建模、仿真、分析和优化。旨在帮助电力系统工程师和研究人员分析电力系统的各种运行和规划问题，以下是MATPOWER的主要特点和功能：

电力系统建模：MATPOWER提供了用于建模电力系统的工具，包括各种元件（发电机、变压器、线路、负荷等）的建模方法，以及电力系统拓扑的表示。

潮流分析：MATPOWER支持直流（DC）和交流（AC）潮流分析，用于计算电力系统中各个节点的电压和相角，以评估系统的稳定性和性能。

稳定性分析：MATPOWER可以用于评估电力系统的稳定性，包括暂态稳定性和静态稳定性。这对于分析系统在各种故障条件下的行为非常有用。

市场运营：MATPOWER可以用于建立电力市场模型，其中发电机、负荷和交易商之间的市场参与者可以进行竞争性交易，以确定电力价格和分配资源。

电力系统优化：MATPOWER包括各种电力系统优化问题的求解器，如潮流调整、机组组合优化、经济调度等，以帮助电力系统规划者和运营商做出最佳决策。

MATPOWER下载：https://matpower.org/download/

安装步骤：

步骤一：将下载好的Matpower压缩包解压，将解压文件夹放在MATLAB安装路径的toolbox文件中。

步骤二：打开MATLAB软件，点击设置路径—添加并包含子文件夹-选中toolbox中Matpower文件—保存—关闭即可，如下图所示：

以上内容对您有帮助吗？如果有的话麻烦给小编一个关注吧，继续推送更有干货的内容！

后面将推送大量电力系统优化调度及预测的代码，复现TOP期刊论文，欢迎关注哦 ~