clear?all;
clc;
format?short;?????
tic

XDCRL=10;??%蓄電池容量
pgridmax=7;%交互功率
m=100;?????%種群規模數
dd=500;????%迭代次數
c1=1;
c2=1;
vmax=1;?
vmin=-1;
wmax=0.9;?
wmin=0.4;?
d=25;%24時刻蓄電池SOC0
CFX?=[];%懲罰項
cfx?=?0;%懲罰項

PV=load（‘Ppv2.txt‘）;%光伏發電量
%?Pev=zeros（124）;
Pev=load（‘PEV.txt‘）;%電動汽車
LOAD=?load（‘Load2.txt‘）;%負荷
sub=load（‘C_sub2.txt‘）;
buy=load（‘C_buy2.txt‘）;
sell=load（‘C_sell2.txt‘）;
xdccl=zeros（m24）;%蓄電池出力
v=zeros（md）;%初始化速度
SOC=rand（md）;%蓄電池的SOC
grid=zeros（m24）;%微網
GRID=zeros（124）;
fx=zeros（1m）;

%%初始化粒子更新
for?n=1:m;?%每一次更新粒子群（100個粒子）都有判斷是否滿足條件
for?i=2:d??%保證蓄電池的SOC大于0.3小于0.95????
????if?SOC（ni）<0.3
????????SOC（ni）=0.3;
????end
????if?SOC（ni）>0.95
????????SOC（ni）=0.95;
????end
????SOC（n1）=0.4;%蓄電池的SOC不變，都為0.4
????SOC（n25）=0.4;%蓄電池的SOC不變，都為0.4
end
?%儲能系統的功率約束
for?i=1:24?
????%蓄電池
????j?=SOC（ni）-SOC（ni+1）;
????if?j?>?0.3
????????j?=?0.3;%放電
????????SOC（ni+1）=?SOC（ni）?-?j;
????elseif?j?<=?-0.3;
????????j?=?-0.3;%充電
?????????SOC（ni+1）=?SOC（ni）?-?j;
????end
????xdccl（ni）?=j*XDCRL;%蓄電池出力?
??
end???
%確定大電網出力
??for?i?=?1:24
???
???grid（ni）=LOAD（i）+Pev（i）?-?PV（i）-?xdccl（ni）;?
???if?grid（ni）>pgridmax??%買電
???????grid（ni）=pgridmax;
??????????
???end
???if?grid（ni）<-pgridmax??%賣電
???????grid（ni）=-pgridmax;
%??????????grid（ni）=LOAD（i）+Pev（i）?-?PV（i）-?xdccl（ni）;
???end
??
??end
cfx?=?（SOC（n1）-SOC（n25））*1000000;
CFX=[CFX;cfx];
end?%初始粒子群的更新結束

%需要帶入函數計算運行成本的單元：大電網出力，蓄電池出力
for?i=1:m
??fx（i）=fitnessEcoVir（CFX（i:）grid（i:）xdccl（i:）subPVbuysell）;
end
[fgbest?best_x]=min（fx）;%fitgbest全局最佳適應度值
gbest=SOC（best_x:）;%全局最優位置?
pbest=SOC;%個體最優位置的初始值，粒子的初始位置
fpbest=fx;%個體最佳適應度值?


for?t=1:dd
????
??for?f=1:m??%m=100
????%量子粒子更新
????alpha=（-0.8*（t-1））/（dd-1）+1.2;
????mbest=sum（pbest）/m;?????
????fi1=rand（1d）;
????fi2=rand（1d）;
????p=（2*fi2.*pbest（f:）+2.1*fi1.*gbest）./（2*fi2+2.1*fi1）;
????u=rand（1d）;
????b=alpha*abs（mbest-SOC（f:））;
????v=-log（u）;
????SOC（f:）=p+（（-1）.^ceil（0.5+rand（1d）））.*b.*v;??
%?????%普通粒子更新
%?????v（f:）?=v（f:）+c1*rand*（pbest（f:）-SOC（f:））+c2*rand*（gbest-SOC（f:））;%%SOC速度更新
%?????SOC（f:）=SOC（f:）+v（f:）;%%SOC位置更新
%?????if（v（f:）>vmax）
%??????????????v（f:）=vmax;
%?????elseif（v（f:）%??????????????v（f:）=vmin;?
%?????end
????for?h=2:d
???????if?SOC（fh）<0.3
????????SOC（fh）=0.3;
???????end
???????if?SOC（fh）>0.95
????????SOC（fh）=0.95;
???????end
????????SOC（f1）=0.4;
????end
????for?i=1:24
????????????j?=SOC（fi）-SOC（fi+1）;
????????????if?j?>?0.3
????????????????j?=?0.3;
????????????????SOC（fi+1）=?SOC（fi）?-?j;
????????????elseif?j<-0.3
????????????????j=-0.3;
????????????????SOC（fi+1）=?SOC（fi）?-?j;
????????????end
????????????xdccl（fi）?=j*XDCRL;??%蓄電池出力
????end
????for?i=1:24
????
?????????grid（fi）=LOAD（i）+Pev（i）?-?PV（i）-?xdccl（fi）;??
??????if?grid（f:）>pgridmax??%買電
?????????grid（f:）=pgridmax;
??????????
??????elseif?grid（f:）<-pgridmax??%賣電
?????????grid（f:）=-pgri

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4629??2019-04-02?17:30??singleobject.m