%清空環境
clear
clc
%參數設置
w_max=0.9;
w_min=0.4;???????????????????????%慣性系數
c1=2;????????????????????????????%加速常數
c2=2;????????????????????????????%加速常數
Dim=3;???????????????????????????%維數
SwarmSize=100;???????????????????%粒子群規模
MaxIter=150;?????????????????????%最大迭代次數

Vmax=1.05;??????????????????????????%最大速度
Vmin=-1.05;?????????????????????????%最小速度

Ppv=[0?0?0?0?0?1?9?16?21?28.4?48?42?15?10?5?1?0?0?0?0?0?0?0?0];???????????????????????????????????????????????????%光伏上界
Pwt=[11?9.2?5.8?3.6?3?3.6?3.2?3.8?2?18?21?9?5.4?2.8?2.6?3.6?7.6?11?9.6?8?2.6?6.6?9.4?11.6];???????????????????????%風電上界
Pmt=[80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80?80];????????????????????????????????????%燃氣上界
Pba=[30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30?30];????????????????????????????????????%蓄電池充放電上界
Pba_1=[-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30?-30];??????????%蓄電池充放電下界
Pl=[53?53?50?50?51?57?63?70?75?77?80?78?74?72?72?77?80?85?85?88?82?78?71?65];

Pbat=60;
Pbatmax=60;
?y_fitness=zeros（124）;???????????????????????%預先生產5個空矩陣
?k_pv=zeros（124）;????????????????????????????%為繪圖做準備
?k_wt=zeros（124）;
?k_mt=zeros（124）;
?k_ba=zeros（124）;
?
for?t=1:24
?Ub=[Ppv（t）?Pwt（t）?Pmt（t）];???????????????????%風光燃參數上界
?Lb=[0?0?30];?????????????????????????????????%風光燃參數下界?
?iter=0;

?%%粒子群初始化
?Range=ones（SwarmSize1）*（Ub-Lb）;
?Swarm=rand（SwarmSizeDim）.*Range+ones（SwarmSize1）*Lb;???????????????????%初始化粒子群+ones（SwarmSize1）*Lb;?
?a4=zeros（SwarmSize1）;
?for?l=1:SwarmSize
?????a4（l）=Pl（t）-Swarm（l1）-Swarm（l2）-Swarm（l3）;
?end
?Swarm=[Swarm?a4];
?VStep=rand（SwarmSizeDim）*（Vmax-Vmin）+Vmin;????????%初始化速度
?fSwarm=zeros（SwarmSize1）;


?for?i=1:SwarmSize
??if?（（a4（i）>Pba（t））||（a4（i）????fSwarm（i:）=20000;
??else
????fSwarm（i:）=0.884*Swarm（i1）+0.573*Swarm（i2）+1.9352*Swarm（i3）+1.3822*Swarm（i4）;??%粒子群適應度計算
??end
?end

%%個體極值和群體極值
?[bestf?bestindex]=min（fSwarm）;
?zbest=Swarm（bestindex1:Dim）;?????????????????????????%全局最佳
?
?pbest=Swarm（:1:Dim