基于改进差分进化算法的微电网调度研究

%% Cuckoo Search Algorithm solving Micro-Grid problems code for Matlab % --------------------------Simple chord------------------ clc clear close all data=parameter; SearchDimension=96; % 问题维度 1-24 FC 25-48 MT 49-72 GS 73-96 GRID Popsize=50; % 种群数量 IterationMax=1000; % 迭代次数 data.N_numb=1;% 算法循环次数 [x_BT,SOC]=BT_dongtaiguihua(data); %%蓄电池的动态规划 Result_DE=zeros(data.N_numb,SearchDimension); BestIterValue_DE=zeros(data.N_numb,IterationMax); tic for i=1:data.N_numb data.i=i; [Result_DE(i,:),BestIterValue_DE(i,:)]=DE(Popsize,SearchDimension,data.SearchScope,IterationMax,data,x_BT); end toc figure hold on % x=[1:IterationMax]; plot(mean(BestIterValue_DE,1)) xlabel('迭代次数') ylabel('费用/元') legend('DE') disp('DE运算均值为') disp(mean(BestIterValue_DE(:,end))) % %% 绘图专用区 % %% 出力-负荷对比图 [~,column]=min(BestIterValue_DE(:,end)); P_vilolate=zeros(1,24); for i=1:24 tempHour=[i 24+i 48+i 72+i]; P_vilolate(1,i)=data.PV(i)+data.WT(i)+x_BT(i)+sum(Result_DE(column,tempHour));% 额外电量 end figure hold on plot(P_vilolate,'o'); plot(data.Demand','+'); plot(P_vilolate); plot(data.Demand'); xlabel('t') ylabel('P') title('各单元总出力-负荷对比图') legend('Micro-Grid system','总负荷','Location','NorthWest') %% 荷电状态 figure a=SOC(2:end); X=(1:24); i=1:24; A=a(1,i); plot(i,A,'o'); hold on plot(X,a); xlabel('时间') ylabel('荷电状态') legend('SOC') axis([0 25 0 1]) %xmin是x最小，xmax是x最大，ymin，ymax类似 % %% 出力图 figure FC=Result_DE(column,1:24); MT=Result_DE(column,25:48); GS=Result_DE(column,49:72); GRID=Result_DE(column,73:96); a=GRID; b=FC; c=MT; d=GS; e=x_BT'; f=data.PV; g=data.WT; h=data.Demand; i=1:24; plot(i,a,'-s',i,b,'-p',i,c,'-*',i,d,'-d',i,e,'-o',i,f,'-+',i,g,'-v',i,h,'-x'); xlabel('时间t/h') ylabel('功率P/KW') legend('GRID','FC','MT','GS','BT','PV','WT','总负荷','Location','NorthWest') figure y=[a' b' c' d' e' f' g']; bar(y,'stacked'); hold on; plot(h,'r','linewidth',2); legend('GRID','FC','MT','GS','BT','PV','WT','总负荷','Location','NorthWest','FontSize',5)