基于粒子群优化算法的微电网调度(光伏、储能、电动车、电网交互)（Matlab代码实现）

%SOC约束
     fx(f)=fitnessEcoVir(CFX,grid(f,:),xdccl(f,:),sub,PV,buy,sell);%适应度更新
%个体最优
    if fx(f)<fpbest(f)
        pbest(f,:)=SOC(f,:);
        fpbest(f)=fx(f);
    end
    %群体最优
    if fx(f)<fgbest
        gbest=SOC(f,:);
        fgbest=fx(f);   
      GRID=grid(f,:);
    end
  end
     MINIMUM(t)=fgbest;  
end
db=gbest;
for j=1:24
    xdc(j)=(gbest(j)-gbest(j+1))*XDCRL;  
end
% for j=1:24
%     GRID(j)=LOAD(j)+Pev(i) - PV(j)- xdc(j);
% end
GRID  % 电网
xdc    %蓄电池
a=sum(xdc)
figure(1);
fitgbest=min(MINIMUM)
plot(MINIMUM,'Color','b','Marker','.','LineWidth',1.5)
xlabel('迭代次数')
ylabel('最优解')
figure(2);
l=1:24;
plot(l,GRID,'o-',l,xdc,'*-',l,PV,'--',l,Pev,'-','LineWidth',1.5)
legend('主网出力','蓄电池出力','光伏出力','电动车出力')
xlabel('小时/h')
ylabel('功率/kW')
% 
% set(gca,'xtick',0:2:24) 
% hold off
SOC=db(1:25)
figure(2)
ll=0:24;
figure(3);
plot(ll,SOC,'Color','r','Marker','.','LineWidth',1.5)
xlabel('小时/h')
ylabel('SOC')
xlim([0 24])
title('蓄电池SOC')
set(gca,'xtick',0:2:24)