量子粒子群算法

glochromx(1,:,:)=chromx(maxfitindex,:,:); %初始化全局最优位置
glofit=fitness(maxfitindex,bestchain); %获得全局最大适应度
%---------程序主循环开始----------------------------
for gen=1:1:maxgen
end
end
%---------初始化全局最优解--------------------------
gloangle(1,:)=angle(maxfitindex,:); %获得全局最优相位
glochrom(1,:,:)=chrom(maxfitindex,:,:); %初始化全局最优位置
hromx(i,j,1)))^0.5)
-(chromx(i,j,2)+1)*cos((abs(chromx(i,j,2)+1-chromx(i,j,1)))^0.5)*sin((abs(chromx(i,j,2)+1+chromx(i
,j,1)))^0.5);
selffitness(i,j)=-selfchromx(i,j,1)*sin((abs(selfchromx(i,j,2)+1-selfchromx(i,j,1)))^0.5)*cos((abs(self
for i=2:1:popsize
for j=1:1:2
if minfit>fitness(i,j)
minfit=fitness(i,j);minfitindex=i;
end
if maxfit<fitness(i,j)
maxfit=fitness(i,j);maxfitindex=i;bestchain=j;
end
end
end
1%
---------计算适应度----------------------------
for i=1:1:popsize
for j=1:1:2
fitness(i,j)=-chromx(i,j,1)*sin((abs(chromx(i,j,2)+1-chromx(i,j,1)))^0.5)*cos((abs(chromx(i,j,2)+1+c
end
avgfit=avgfit+fitness(i,j)/(2*popsize);
end
end
%---------更新自身最优位置---------------------
for i=1:1:popsize
t1=fitness(i,1);
if t1<fitness(i,2)
chromx(i,j,2)+1+selfchromx(i,j,1)))^0.5)
-(selfchromx(i,j,2)+1)*cos((abs(selfchromx(i,j,2)+1-selfchromx(i,j,1)))^0.5)*sin((abs(selfchromx(i,j,
2)+1+selfchromx(i,j,1)))^0.5);
chromx(i,j,k)=0.5*(varrange(k,2)*(1+chrom(i,j,k))+varrange(k,1)*(1-chrom(i,j,k)));
selfchromx(i,j,k)=0.5*(varrange(k,2)*(1+selfchrom(i,j,k))+varrange(k,1)*(1-selfchrom(i,j,k)));
t1=t1-2*pi;
end
t2=gloangle(1,j)-angle(i,j);
if t2<-pi
t2=t2+2*pi;
end
if t2>pi
t2=t2-2*pi;
end
dangle(i,j)=inertia*dangle(i,j)+selfw*rand*t1+globalw*rand*t2;
for i=1:1:popsize
for j=1:1:vartotal
if rand<mutatep
angle(i,j)=pi/2-angle(i,j); %变异后的相位
chrom(i,1,j)=cos(angle(i,j)); %变异后的位置
chrom(i,2,j)=sin(angle(i,j)); %变异后的位置
angle(i,j)=angle(i,j)+dangle(i,j); %移动后的相位
chrom(i,1,j)=cos(angle(i,j)); %移动后的位置
chrom(i,2,j)=sin(angle(i,j)); %移动后的位置
%解空间变换
chromx(i,1,j)=0.5*(varrange(j,2)*(1+chrom(i,1,j))+varrange(j,1)*(1-chrom(i,1,j)));
,j,1)))^0.5);%计算适应度
end
end
%---------统计适应度----------------------------
minfit=fitness(1,1);
if fitness(1,1)>fitness(1,2)
minfit=fitness(1,2);
end
end
end
for j=1:1:2
fitness(i,j)=-chromx(i,j,1)*sin((abs(chromx(i,j,2)+1-chromx(i,j,1)))^0.5)*cos((abs(chromx(i,j,2)+1+c
hromx(i,j,1)))^0.5)
-(chromx(i,j,2)+1)*cos((abs(chromx(i,j,2)+1-chromx(i,j,1)))^0.5)*sin((abs(chromx(i,j,2)+1+chromx(i
%---------粒子位置移动------------------------------
for i=1:1:popsize
for j=1:1:vartotal
t1=selfangle(i,j)-angle(i,j);
if t1<-pi
t1=t1+2*pi;
end
if t1>pi
maxfit=fitness(1,2);bestchain=2;
if fitness(1,1)>fitness(1,2)
maxfit=fitness(1,1);bestchain=1;
end
minfitindex=1;maxfitindex=1;
avgfit=(fitness(1,1)+fitness(1,2))/(2*popsize);
chrom(i,1,j)=cos(angle(i,j)); %第 i 个粒子（量子染色体）的幅角余弦
chrom(i,2,j)=sin(angle(i,j)); %第 i 个粒子（量子染色体）的幅角正弦
selfangle(i,j)=2*pi*rand;
selfchrom(i,1,j)=cos(angle(i,j));%第 i 个粒子自身最优位置的幅角余弦
hromx(i,j,1)))^0.5)
-(chromx(i,j,2)+1)*cos((abs(chromx(i,j,2)+1-chromx(i,j,1)))^0.5)*sin((abs(chromx(i,j,2)+1+chromx(i
,j,1)))^0.5);%计算适应度
end
end
%---------粒子位置变异------------------------
selfchromx(i,:,:)=chromx(i,:,:);
selffitness(i,:)=fitness(i,:);
end
end
%---------更新全局最优位置--------------------
if glofit<fitness(maxfitindex,bestchain)
end
end
%---------适应度统计------------------------------
minfit=fitness(1,1);
if fitness(1,1)>fitness(1,2)
minfit=fitness(1,2);
end
maxfit=fitness(1,2);bestchain=2;
for j=1:1:2
if minfit>fitness(i,j)
minfit=fitness(i,j);minfitindex=i;
end
if maxfit<fitness(i,j)
maxfit=fitness(i,j);maxfitindex=i;bestchain=j;
end
avgfit=avgfit+fitness(i,j)/(2*popsize);
if fitness(1,1)>fitness(1,2);
maxfit=fitness(1,1);bestchain=1;
end
minfitindex=1;maxfitindex=1;
avgfit=(fitness(1,1)+fitness(1,2))/(2*popsize);