《量子粒子群算法》由会员分享，可在线阅读，更多相关《量子粒子群算法（3页珍藏版）》请在金锄头文库上搜索。
1、程序正文cleoasre aalll;%变量部分-po-p-si-z-e=50;%种群规模vartotal=2;%变量个数inertia=;%惯性因子selfw=;%自身因子globalw=;%全局因子mutatep=; %变异概率 maxgen=500; %限定代数%数组部分 varrange(1,1)=5社2第第变量变量最小值varrange(2,1)=-512; %第二变量最小值varrange(2,2)=512 % 第二变量最大值%粒子位置初始化 for for：1：popsizej=1:a1n:gvlaer(tio,tja)l=2*pi*rand;rom(i,1,j)二cos(ang
2、le(i,%第i个粒子(量子染色体)的幅角余弦 j)：%第i个粒子(量子染色体)的幅角正弦chrom(i,2,j)=sin(angle(i,S)l;chrom(i,1,j)=cos(angle(i,j)：%第 i个粒子自身最优位置的幅角余弦 sel fcnrAe(i，刃=2*p i *rand ie(i,j)；%第i个粒子自身最优位置的幅角正弦 dangle(i,j)=0：endend1%解空间变换for fio=r1:j1=:1p:o1p:s2izefor k=1:c1h:rvoamrxt(oit,ajl,k)=*(varrange(k,2)*(1+chrom(i,j,k)+varrange
3、(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%计算适应度 for fio=r1:j1=:1p:o1p:s2ize 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,
4、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)+fchefitness(1, 2e)nd maxfmiitn=ffiitt=nfeistsn(e1s,s2()1；,
5、b2e)s；tchain=2； if fitn ess(1,1)fi tness(1,2);maxfit=fitness(1,1)；bestchain= 1;emnidnfitindex=1;maxfitindex=1; avgfit=(fitness(1,1)+fitness(1,2)/(2*pop size);for fio=r2:j1=:1p:o1p:s2izeifminmfintfifti=tfnietsnse(sis,(ji),j);minfitind ex=i;einfd maxfmiatxfit=nfeistsn(eis,sj()i,j);maxfitindex=i;bestch
6、ai n=j;aevndgfit=avgfit+fitness(i,j)/(2*popsize);endend%初始化全局最优解g-l-o-an-g-l-e(1,:)=angle(maxfitindex,%获得全局最优相位:g l) ;ochrom (1 , :, :) =chrom (maxf i t i ndex, :, :%)初;始化全局最优位置 glochromx(1,:,:)=chromx(maxfi tindex,% 初始化全局最优位置 :);%获得全局最大适应度%-ofi-#-程序主循环开 ndT-%-_-妙三-粒m子位置移动for i=1:1:popsizefor j=1:t1
7、1:=vsaerltfoatnaglle(i,j)-ianf glt1epit1=t1-2*pi;endt2=gloangle(1,j)- ianf glt2epit2=t2-2*pi;end dangle(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+chro
8、m(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- for chrom：(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+c hromx(i,j,1)厂-(chromx(i,j,2)+1)*cos(abs(chromx(i,j,2)+1-,jromx)i;j计算)适应度(abs(chromx
9、(i,j,2)+1+chromx(i endend%粒子位置变异for i=1:1:popsizefor j=1:1:vartotalif randmutatep%变异后的相位%变异后的位置%变异后的位置angle(i,j)=pi/2- angle(i,j)； chrom(i,1,j)=cos(angle(i,j)；chrom(i,2,j)=sin(angle(i,chromx(i,1,j)=* chrom(i,1,j)；j解空間变换2)*(1+chrom(i,1,j)+varrange(j,1)*(1-chromx(i,2,j)=*(varrange(j,2)*(1+chrom(i,2,j)+varrange(j,1)*(1- chrom(i,2,j)；
