✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。
🍎个人主页：Matlab科研工作室
🍊个人信条：格物致知。
更多Matlab仿真内容点击👇
⛄ 内容介绍
复现论文
受多能源整合优势的启发，本文建立了一个基于交互控制的双层两阶段框架，以实现互联多能源系统（MES）之间的最优能源供应。在较低层次上，每个 MES 通过解决成本最小化问题自主确定每个可控资产的最佳设定点，其中采用滚动水平优化来处理负载和可再生能源的随机特性。进一步实现了放宽存储互补约束优化模型凸化的技术，其数学证明验证了放宽的正确性。在上层，建立了一个协调器，以最小化协作互连 MES 的总成本，同时防止变压器过载。这个协作问题在基于市场清算机制的两阶段过程中进一步分解和迭代解决。该方法的一个显着特点是它兼容操作时间要求，同时保留了每个MES的可扩展性、信息隐私和操作权限。通过对自主协作优化机制进行详细分析的模拟案例验证了所提出框架的有效性。
⛄ 部分代码
%%文章复现《Collaborative Autonomous Optimization of Interconnected Multi-Energy
%%Systems with Two-Stage Transactive Control Framework》
%%绘图部分xticks函数可能2017以下版本无法使用，如果报错把它去掉
%%绘图部分比较简单，我不做过多注释
%%数据读取函数，依次读取绘图需要的数据以及变量值
%%绘图函数以及数据读取函数
function [] = drawMES(t, result_Ele, ymin, ymax)
figure;
H1 = bar(t, result_Ele / 1000); hold on;
H1(1).EdgeColor = 'none';
H1(2).EdgeColor = 'none';
H1(3).EdgeColor = 'none';
color_mes1 = ColorHex('4083ff') / 255;
color_mes2 = ColorHex('005aff') / 255;
color_mes3 = ColorHex('3200ff') / 255;
H1(1).FaceColor = color_mes1;
H1(2).FaceColor = color_mes2;
H1(3).FaceColor = color_mes3;
le = legend([H1(1), H1(2), H1(3)],'MES_1', 'MES_2', 'MES_3','Orientation','horizontal');
ylabel('power(MW)');
xlim([0, 25]);
ylim([ymin, ymax]/1000);
xticks(0: 6 : 24);
xticklabels({ '0:00','6:00','12:00','18:00','24:00' });
set(le,'Box','off');
set(gcf,'Position',[0 0 500 200]);
end
⛄ 运行结果
⛄ 参考文献
[1] Cheng, Y. , P. Zhang , and X. Liu . "Collaborative Autonomous Optimization of Interconnected Multi-Energy Systems with Two-Stage Transactive Control Framework." MDPI, Open Access Journal, 10.3390/en13010171. 2019.
