%%基本粒子群算法
clear?all
clc
tic
%%參數設置
maxgen=100;???????%迭代次數
sizepop=20;???????%種群大小
c=1.49445;???????%速度更新參數，c1=c2
wmax=0.9;???????????%慣性權重
wmin=0.4
popmax=[1.06?8?5];?%發電機節點電壓、變壓器分接頭、電容器投入組數上限
popmin=[0.94?0?0];?%發電機節點電壓、變壓器分接頭、電容器投入組數下限
vmax=[0.2?2?2];????%速度上限
vmin=[-0.2?-2?-2];?%速度下限
dim=[6?4?2];???????%變量維數
shuju=case_ieee30;?%得到30節點的數據
%%初始化種群

for?i=1:sizepop
????pop（i1:dim（1））?????????????????????????????=rand（1dim（1））*（popmax（1）-popmin（1））+popmin（1）;?%初始化節點電壓
????pop（idim（1）+1:dim（1）+dim（2））???????????????=round（popmax（2）*rand（1dim（2）））;????????????????%初始化分接頭位置，取整
????pop（idim（1）+dim（2）+1:dim（1）+dim（2）+dim（3））?=round（popmax（3）*rand（1dim（3）））;????????????????%初始化投入組數，取整
????v（i1:dim（1））???????????????????????????????=0.2*rands（1dim（1））;
????v（idim（1）+1:dim（1）+dim（2））?????????????????=2*rands（1dim（2））;
????v（idim（1）+dim（2）+1:dim（1）+dim（2）+dim（3））???=2*rands（1dim（3））;
????%計算網損大小
????shuju.bus（[1?2?5?8?11?13]8）????????????????=pop（i1:dim（1））;????????????????????????????????%改30節點中的節點電壓
????shuju.gen（:6）??????????????????????????????=pop（i1:dim（1））;????????????????????????????????%改30節點中的節點電壓
????shuju.branch（[11?12?15?36]9）???????????????=pop（idim（1）+1:dim（1）+dim（2））*0.025+0.9;??????????????????%改變壓器變比
????shuju.bus（[10?24]6）????????????????????????=pop（idim（1）+dim（2）+1:dim（1）+dim（2）+dim（3））*10;????%改電容器容量
????[basemva?bus?gen?branch?success?et]?????????=runpf（shuju）;
????fitness（i）??????????????????????????????????=sum（branch（:14）+branch（:16））;????????????????%計算出適應度大小，即網損大小
end
[best?index]=min（fitness）;
gtsite=pop;?????????????%個體最佳位置
gtfit=fitness???????????%個體最佳適應度
popsite=pop（index:）;???%全局最佳位置
popfit=fitness（index）;??%全局最佳適應度
trace（1）=popfit;
for?j=1:maxgen
????w=wmax-（wmax-wmin）/maxgen*j;
????for?i=1:sizepop
????????%速度更新
????????v（i:）=w*v（i:）+c*rand*（gtsite（i:）-pop（i:））+c*rand*（popsite-pop（i:））;
????????%限制微粒速度
????????v1=v（i1:dim（1））;???????????????????????????????%節點電壓速度
????????v2=v（idim（1）+1:dim（1）+dim（2））;?????????????????%變比速度
????????v3=v（idim（1）+dim（2