量子粒子群算法-程序

程序正文cleoasreaalll;;%變量部分-po--p-si-z--e=50;%種群規(guī)模vartotal=2;%變量個數(shù)inertia=;%慣性因子selfw=;%自身因子globalw=;%全局因子mutatep=;%變異概率maxgen=500;%限定代數(shù)%數(shù)組部分varrange(1,1)=5社2第第變量變量最小值varrange(2,1)=-512;%第二變量最小值varrange(2,2)=512%第二變量最大值%粒子位置初始化forfor：1：popsizej=1:a1n:gvlaer(tio,tja)l=2*pi*rand;rom(i,1,j)二cos(angle(i,%第i個粒子(量子染色體)的幅角余弦j))：%第i個粒子(量子染色體)的幅角正弦chrom(i,2,j)=sin(angle(i,S)l;chrom(i,1,j)=cos(angle(i,j))：%第i個粒子自身最優(yōu)位置的幅角余弦selfcnrAe(i，刃=2*pi*randie(i,j))；%第i個粒子自身最優(yōu)位置的幅角正弦dangle(i,j)=0：endend1%解空間變換forfio=r1:j1=:1p:o1p:s2izefork=1:c1h:rvoamrxt(oit,ajl,k)=*(varrange(k,2)*(1+chrom(i,j,k))+varrange(k,1)*(1-chrom(i,j,k)))iselfchromx(i,j,k)=*(varrange(k,2)*(1+selfchrom(i,j,k))+varrange(k,1)*(1-selfchernodm(i,j,k)))iendend1%計算適應度forfio=r1:j1=:1p:o1p:s2izefitness(i,j)=-chromx(i,j,1)*sin((abs(chromx(i,j,2)+1-chromx(i,j,1)))^*cos((abs(chromx(i,j,2)+1+chromx(i,j,1))廠-(chromx(i,j,2)+1)*cos((abs(chromx(i,j,2)+1-chromx(i,j,1)))^*sin((abs(chromx(i,j,2)+1+chromx(i,j,1))廠；selffitness(i,j)=-selfchromx(i,j,1)*sin((abs(selfchromx(i,j,2)+1-selfchromx(i,j,1)))^*cos((abs(selfchromx(i,j,2)+1+selfchromx(i,j,1))廠-(selfchromx(i,j,2)+1)*cos((abs(selfchromx(i,j,2)+1-2)+fche<Pfc(iomxti,；abs(selfchromx(i,j,endend%適應度統(tǒng)計miifnfit=fitness(1,1)；fitness(1,1)>fitness(1,2e)ndmaxfmiitn=ffiitt=nfeistsn(e1s,s2()1；,b2e)s；tchain=2；iffitness(1,1)>fitness(1,2);maxfit=fitness(1,1)；bestchain=1;emnidnfitindex=1;maxfitindex=1;avgfit=(fitness(1,1)+fitness(1,2))/(2*popsize);forfio=r2:j1=:1p:o1p:s2izeifminmfintf>ifti=tfnietsnse(sis,(ji),j);minfitindex=i;einfdmaxfmiatx<fit=nfeistsn(eis,sj()i,j);maxfitindex=i;bestchain=j;aevndgfit=avgfit+fitness(i,j)/(2*popsize);endend%初始化全局最優(yōu)解g--l-o-an-g-l--e(1,:)=angle(maxfitindex,%獲得全局最優(yōu)相位:gl);ochrom(1,:,:)=chrom(maxfitindex,:,:%)初;始化全局最優(yōu)位置glochromx(1,:,:)=chromx(maxfitindex,%初始化全局最優(yōu)位置:);%獲得全局最大適應度%-ofi-#--程序主循環(huán)開nd—T----------%--_-妙三--粒m子位置移動fori=1:1:popsize

forj=1:t11:=vsaerltfoatnaglle(i,j)-ianfglt1e<(-ip,ij);t1=t1+2*pi;endift1>pit1=t1-2*pi;endt2=gloangle(1,j)-ianfglt2e<(-ip,ij);t2=t2+2*pi;endift2>pit2=t2-2*pi;enddangle(i,j)=inertia*dangle(i,j)+selfw*rand*t1+globalw*rand*an；le(i,j)二angle(i,j)+dangl%移動后的相位(i,j);%移動后的位置chrom(i,1,j)二cos(angle(i,j))%移動后的位置chrom(i,2,j)=sin(angle(i,j));C解空間變換j)=*(varrange(j,2)*(1+chrom(i,1,j))+varrange(j,1)*(1-chrom(i,1,j)));endchromx(i,2,j)=*(varrange(j,2)*(1+chrom(i,2,j))+varrange(j,1)*(1-forchrom：(22,j)))；fitness(i,j)=-chromx(i,j,1)*sin((abs(chromx(i,j,2)+1-chromx(i,j,1)))^*cos((abs(chromx(i,j,2)+1+chromx(i,j,1))廠-(chromx(i,j,2)+1)*cos((abs(chromx(i,j,2)+1-,jromx)i;j計算))適應度((abs(chromx(i,j,2)+1+chromx(iendend%粒子位置變異fori=1:1:popsizeforj=1:1:vartotalifrand<mutatep%變異后的相位%變異后的位置%變異后的相位%變異后的位置%變異后的位置j))；chrom(i,2,j)=sin(angle(i,chromx(i,1,j)=*chrom(i,1,j)))；j解空間變換chromx(i,