量子粒子群算法程序.doc

《量子粒子群算法程序.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、%---------程序正文-------------------clearall;closeall;%---------变量部分-------------------popsize=50;vartotal=2;%种群规模%变量个数%惯性因子%自身因子%全局因子%变异概率%限定代数inertia=0.5;selfw=2.0;globalw=2.0;mutatep=0.05;maxgen=500;%---------数组部分-----------------------varrange(1,1)=-512;%第一变量最小值varrange(1,2

2、)=512;%第一变量最大值varrange(2,1)=-512;%第二变量最小值%第二变量最大值varrange(2,2)=512;%---------粒子位置初始化-----------------fori=1:1:popsizeforj=1:1:vartotalangle(i,j)=2*pi*rand;chrom(i,1,j)=cos(angle(i,j));chrom(i,2,j)=sin(angle(i,j));selfangle(i,j)=2*pi*rand;%第i个粒子（量子染色体）的幅角余弦%第i个粒子（量子染色体）的幅角正弦s

3、elfchrom(i,1,j)=cos(angle(i,j));%第i个粒子自身最优位置的幅角余弦selfchrom(i,2,j)=sin(angle(i,j));%第i个粒子自身最优位置的幅角正弦dangle(i,j)=0;endend1%---------解空间变换------------------------fori=1:1:popsizeforj=1:1:2fork=1:1:vartotalchromx(i,j,k)=0.5*(varrange(k,2)*(1+chrom(i,j,k))+varrange(k,1)*(1-chrom(

4、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)));endendend1%---------计算适应度----------------------------fori=1:1:popsizeforj=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,

5、j,2)+1+chromx(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(selfchromx(i,j,2)+1+selfchromx(i,j,1))

6、)^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);endend%---------适应度统计------------------------------minfit=fitness(1,1);iffitness(1,1)>fitness(1,2)minfit=fitness(1,2);endmaxfit=fitness(1,

7、2);bestchain=2;iffitness(1,1)>fitness(1,2);maxfit=fitness(1,1);bestchain=1;endminfitindex=1;maxfitindex=1;avgfit=(fitness(1,1)+fitness(1,2))/(2*popsize);fori=2:1:popsizeforj=1:1:2ifminfit>fitness(i,j)minfit=fitness(i,j);minfitindex=i;endifmaxfit

8、j);maxfitindex=i;bestchain=j;endavgfit=avgfit+fitness(i,j)/(2*popsize);ende