基本粒子群算法的ieee30节点无功优化
%%基本粒子群算法clear all
clc
tic
%%参数设置
maxgen=100; %迭代次数
sizepop=20; %种群大小
c=1.49445; %速度更新参数，c1=c2
wmax=0.9; %惯性权重
wmin=0.4
popmax=; %发电机节点电压、变压器分接头、电容器投入组数上限
popmin=; %发电机节点电压、变压器分接头、电容器投入组数下限
vmax=; %速度上限
vmin=[-0.2 -2 -2]; %速度下限
dim=; %变量维数
shuju=case_ieee30; %得到30节点的数据
%%初始化种群
for i=1:sizepop
pop(i,1:dim(1)) =rand(1,dim(1))*(popmax(1)-popmin(1))+popmin(1); %初始化节点电压
pop(i,dim(1)+1:dim(1)+dim(2)) =round(popmax(2)*rand(1,dim(2))); %初始化分接头位置，取整
pop(i,dim(1)+dim(2)+1:dim(1)+dim(2)+dim(3)) =round(popmax(3)*rand(1,dim(3))); %初始化投入组数，取整
v(i,1:dim(1)) =0.2*rands(1,dim(1));
v(i,dim(1)+1:dim(1)+dim(2)) =2*rands(1,dim(2));
v(i,dim(1)+dim(2)+1:dim(1)+dim(2)+dim(3)) =2*rands(1,dim(3));
%计算网损大小
shuju.bus(,8) =pop(i,1:dim(1)); %改30节点中的节点电压
shuju.gen(:,6) =pop(i,1:dim(1)); %改30节点中的节点电压
shuju.branch(,9) =pop(i,dim(1)+1:dim(1)+dim(2))*0.025+0.9; %改变压器变比
shuju.bus(,6) =pop(i,dim(1)+dim(2)+1:dim(1)+dim(2)+dim(3))*10; %改电容器容量
=runpf(shuju);
fitness(i) =sum(branch(:,14)+branch(:,16)); %计算出适应度大小，即网损大小
end
=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(i,1:dim(1)); %节点电压速度
v2=v(i,dim(1)+1:dim(1)+dim(2)); %变比速度
v3=v(i,dim(1)+dim(2)+1:dim(1)+dim(2)+dim(3)); %电容器容量速度
v1(find(v1>vmax(1)))=vmax(1);
v1(find(v1<vmin(1)))=vmin(1);
v2(find(v2>vmax(2)))=vmax(2);
v2(find(v2<vmin(2)))=vmin(2);
v3(find(v3>vmax(3)))=vmax(3);
v3(find(v3<vmin(3)))=vmin(3);
v(i,:)=;
%位置更新
pop(i,:)=pop(i,:)+v(i,:);
%限制微粒位置
pop1=pop(i,1:dim(1));
pop2=pop(i,dim(1)+1:dim(1)+dim(2));
pop3=pop(i,dim(1)+dim(2)+1:dim(1)+dim(2)+dim(3));
pop1(find(pop1>popmax(1)))=popmax(1);
pop1(find(pop1<popmin(1)))=popmin(1);
pop2(find(pop2>popmax(2)))=popmax(2);
pop2(find(pop2<popmin(2)))=popmin(2);
pop3(find(pop3>popmax(3)))=popmax(3);
pop3(find(pop3<popmin(3)))=popmin(3);
pop2=round(pop2);
pop3=round(pop3);
pop(i,:)=;
end
for i=1:sizepop
%计算适应度
shuju.bus(,8) =pop(i,1:dim(1)); %改30节点中的节点电压
shuju.gen(:,6) =pop(i,1:dim(1)); %改30节点中的节点电压
shuju.branch(,9) =pop(i,dim(1)+1:dim(1)+dim(2))*0.025+0.9; %改变压器变比
shuju.bus(,6) =pop(i,dim(1)+dim(2)+1:dim(1)+dim(2)+dim(3))*10; %改电容器容量
=runpf(shuju);
fitness =sum(branch(:,14)+branch(:,16)); %计算出适应度大小，即网损大小
%个体适应度及位置更新
if fitness<gtfit(i)
gtfit(i)=fitness;
gtsite(i,:)=pop(i,:);
end
%全局适应度及位置更新
if fitness<popfit
popfit=fitness;
popsite=pop(i,:);
end
end
trace(j+1)=popfit; %记录每次迭代的最佳值
end
figure
plot(trace)
toc
页:
[1]