function?PSOfirst（）
%%?清空環(huán)境
clear;
clc;

%%?參數(shù)設(shè)置
wcmax=0.8;
wcmin=0.7;%權(quán)值將影響PSO?的全局與局部搜優(yōu)能力，
%w的值較大，全局搜優(yōu)能力強，局部搜優(yōu)能力弱;反之，則局部搜優(yōu)能力增強，而全局搜優(yōu)能力減弱。
c1=0.01;%加速度，影響收斂速度
c2=0.01;
dim=24;%每一維表示在一個小時內(nèi)的交換功率
swarmsize=100;%粒子群規(guī)模，表示有100個解的空間
maxiter=80;%最大循環(huán)次數(shù)，影響時間
vmax=1;
vmin=-1;
PENUM=[192020201718171718191717201718171816181919192020];%由于電動車主的隨機出行，因此每一時刻能完成V2G功能的電動汽車數(shù)量不是固定的。

ub=8.*PENUM;%解向量的最大限制，就是充電功率上限
lb=-3.*PENUM;%解向量的最小限制，即放電功率上限

%%?種群初始化
range=ones（swarmsize1）*（ub-lb）;
swarm=rand（swarmsizedim）.*range+ones（swarmsize1）*lb;%粒子群位置矩陣
%對粒子進行更詳細的SOC約束，充放電不能使SOC超出預(yù)設(shè)范圍


%定義速度和位置???
vstep=rand（swarmsizedim）*（vmax-vmin）+vmin;%粒子群速度矩陣
fswarm=zeros（swarmsize1）;%預(yù)設(shè)空矩陣，存放適應(yīng)值
for?i=1:swarmsize
????X=swarm（i:）;
????[SUMEVPAf3f2]=f（X）;%目標函數(shù)，適應(yīng)度函數(shù)