基于双层共识控制的直流微电网优化调度（Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文章

---

💥1 概述

微电网 (mG) 是由分布式发电单元 (DGU)、存储设备和负载组成的电网。除了可再生能源整合、提高电能质量、降低传输损耗、能够在并网和孤岛模式下运行等诸多优势外，它们还与交流和直流运行标准兼容 [1]。特别是，直流微电网 (DCmGs) 最近获得了关注。它们的日益普及可归因于高效转换器的开发、与可再生能源（例如光伏模块）和电池的自然接口，以及本质上为直流的电子负载（各种电器、LED、电动汽车、计算机等）的可用性。

在孤岛 DCmG 中，电压稳定性至关重要，因为没有它，电压可能会突破临界水平并损坏连接的负载 [3]。因此，通常采用初级电压控制层来跟踪耦合点 (PC) 处的所需电压参考，从而将 DGU 连接到 DCmG。为此，文献中提出了几种方法，例如基于下垂控制[2]、[4]和即插即用控制[5]、[6]。除了电压稳定性，另一个理想的目标是确保电流共享，即 DGU 必须根据其额定电流共享 mG 负载。事实上，在没有它的情况下，未调节的电流可能会使发电机过载并最终导致 mG 故障。

本文提出了用于直流微电网中电流共享和电压平衡的新型二级调节器，由分布式发电机组、动态 RLC 线路和非线性 ZIP（恒定阻抗、恒定电流和恒定功率）负载组成。所提出的控制器位于初级电压控制层之上，并利用通过通信网络交换的信息来采取必要的控制措施。我们推导出平衡点存在和唯一性的充分条件，并证明在稳态下达到了预期目标。我们的控制设计只需要了解发电机组的本地参数，便于即插即用操作。我们提供电压稳定性分析，并通过仿真说明我们设计的性能和稳健性。所有结果都适用于任意（尽管是连接的）微电网和通信网络拓扑。

📚2 运行结果

复现结果图：

 

 

 

原文运行结果图：

可视化代码：

%% Plotting results
close all;
figure(1);
plot(t,y(:,1:N),'linewidth',2);
set(gcf, 'Position',[189, 611,560,310]);
set(gca,'fontsize',15, 'FontName', 'Times New Roman')
hold on
if expo==0 && det(P_l)~=0 && com==1
plot(t, ones(size(t))*Vmax,'b--','linewidth',2);
hold on
plot(t, ones(size(t))*Vmin,'b--','linewidth',2);
end
xlabel('Time (seconds)')
ylabel('DGU voltage')
grid on
legend('$V_1$','$V_2$','$V_3$','$V_4$','$V_5$','$V_6$','Location','north','Interpreter','latex','Orientation','horizontal')
hold off plots

figure(2);
plot(t,y(:,N+1:2*N),'linewidth',2);
set(gcf, 'Position',[189, 611,560,310]);
set(gca,'fontsize',15, 'FontName', 'Times New Roman')
legend('$I_{t1}$','$I_{t2}$','$I_{t3}$','$I_{t4}$','$I_{t5}$','$I_{t6}$','Location','southeast','Orientation','horizontal','Interpreter','latex')
xlabel('Time (seconds)')
ylabel('Filter currents')
grid on
hold off plots

if com==1
figure(3);
plot(t,y(:,N+1:2*N)*D,'linewidth',2);
set(gcf, 'Position',[189, 611,560,310]);
set(gca,'fontsize',15, 'FontName', 'Times New Roman')
leg=legend('$\frac{I_{t1}}{\bar{I}_{t1}}$','$\frac{I_{t2}}{\bar{I}_{t2}}$','$\frac{I_{t3}}{\bar{I}_{t3}}$','$\frac{I_{t4}}{\bar{I}_{t4}}$','$\frac{I_{t5}}{\bar{I}_{t5}}$','$\frac{I_{t6}}{\bar{I}_{t6}}$','Location','southeast','Orientation','horizontal','Interpreter','latex');
leg.FontSize = 20;
xlabel('Time (seconds)')
ylabel('Weighted filter currents')
grid on
hold off plots

figure(4);

plot(t,y(:,1:N)*D^(-1)*ones(N,1),'r','linewidth',2);
hold on
plot(t, (ones(1,N)*D^(-1)*Vref)*ones(size(t)),'b--','linewidth',2.5);
set(gcf, 'Position',[189, 611,560,310]);
set(gca,'fontsize',15, 'FontName', 'Times New Roman')
xlabel('Time (seconds)')
ylabel('Voltage regulation')
grid on
hold off plots
end

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

🌈4 Matlab代码、数据、文章