个人主页:研学社的博客
欢迎来到本博客❤️❤️
博主优势:博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。
⛳️座右铭:行百里者,半于九十。
本文目录如下:
目录
1 概述
以含有光伏(PV)、风机 (W)、柴油发电机(de)、燃料电池(fc)、蓄电池(bat)的子微网1、子微网2 构成微网群进行算例分析, 微网群算例系统拓扑结构如下图。
微网群采用并联式结构,其中子网1中分布式电源包括:PV1, W1,fc1,bat1;子网2中分布式电源包括PV2,W2,de2,bat2。微网群与配电网之间的公共连接点群PCC保持闭合,网PCC1与网PCC2也均保持闭合。微网1,2和配电网之间可存在功率交换Pbuy1,,Pbuy2,,微网1与微网2之间可存在功率交换Pex1,Pex2。
本文建立的微网群以多目标协调优化为目的,其中各个子微网优化目标与总微网群的优化目标一致。微网群优化调度周期取一天电力调度论文,以1h作为一个调度时段,全天共分为24个时段电力调度论文,且算例中采用分时电价a1,a2,b1,b2的方式。
数学模型及详细文章讲解见第4部分。
2 运行结果
部分代码:
[best index]=min();
=pop; %个体最佳位置
gtfit=; %个体最佳适应度
=pop(index,:); %全局最佳位置
=(index); %全局最佳适应度
trace(1)=;
for k=1:
w=wmax-(wmax-wmin)/*j;
for i=1:
%速度更新
v(i,:)=w*v(i,:)+c*rand*((i,:)-pop(i,:))+c*rand*(-pop(i,:));
%限制微粒速度
v(find(v>vmax))=vmax;
v(find(v %位置更新
pop(i,:)=pop(i,:)+v(i,:);
%限制微粒位置
pop1=pop(i,1:dim);
pop2=pop(i,dim+1:2*dim);
pop3=pop(i,2*dim+1:3*dim);
pop4=pop(i,3*dim+1:4*dim);
pop5=pop(i,4*dim+1:5*dim);
pop1(find(pop1>(1)))=(1);
pop1(find(pop1
pop2(find(pop2>(2)))=(2);
pop2(find(pop2
pop3(find(pop3>(3)))=(3);
pop3(find(pop3
pop4(find(pop4>(4)))=(4);
pop4(find(pop4
pop5(find(pop5>(5)))=(5);
pop5(find(pop5
pop(i,:)=[pop1 pop2 pop3 pop4 pop5];
for j=1:14
pop(i,dim+j)=0;
pop(i,4*dim+j)=0;
end
pop(i,dim+23)=0;
pop(i,4*dim+23)=0;
pop(i,dim+24)=0;
pop(i,4*dim+24)=0;
for j=1:24
if abs(pop(i,j))>0&&abs(pop(i,j)) pop(i,j)=0;
end
%if abs(pop(i,3*dim+j))>0&&abs(pop(i,3*dim+j))
% pop(i,3*dim+j)=0;
% end
% if abs(pop(i,dim+j))>0&&abs(pop(i,dim+j)) %=+10000;
% pop(i,dim+j)=0;
%end
if abs(pop(i,4*dim+j))>0&&abs(pop(i,4*dim+j)) %=+10000;
pop(i,4*dim+j)=0;
end
if j=16
pop(i,2*dim+j)=0;
end
end
end
for i=1:
%计算适应度
%其他参数求解pbuy psell
pex2=-pop(i,2*dim+1:3*dim);
bs1=Load1-pop(i,1:dim)-Ppv1-pop(i,dim+1:2*dim)-Pw1-pop(i,2*dim+1:3*dim);
bs2=Load2-pop(i,3*dim+1:4*dim)-Ppv2-pex2-Pw2-pop(i,4*dim+1:5*dim);
3参考文献
文章来源:https://blog.csdn.net/weixin_46039719/article/details/128794190