算法源码-优化与控制：量子遗传算法代码.zip

function chrom=Qgate(chrom,fitness,best,binary) %% 量子旋转门调整策略 % 输入 chrom：更新前的量子比特编码 % fitness：适应度值 % best：当前种群中最优个体 % binary：二进制编码 % 输出 chrom：更新后的量子比特编码 sizepop=size(chrom,1)/2; lenchrom=size(binary,2); for i=1:sizepop for j=1:lenchrom A=chrom(2*i-1,j); % α B=chrom(2*i,j); % β x=binary(i,j); b=best.binary(j); if ((x==0)&(b==0))||((x==1)&(b==1)) delta=0; % delta为旋转角的大小 s=0; % s为旋转角的符号，即旋转方向 elseif (x==0)&(b==1)&(fitness(i)<best.fitness) delta=0.01*pi; if A*B>0 s=1; elseif A*B<0 s=-1; elseif A==0 s=0; elseif B==0 s=sign(randn); end elseif (x==0)&(b==1)&(fitness(i)>=best.fitness) delta=0.01*pi; if A*B>0 s=-1; elseif A*B<0 s=1; elseif A==0 s=sign(randn); elseif B==0 s=0; end elseif (x==1)&(b==0)&(fitness(i)<best.fitness) delta=0.01*pi; if A*B>0 s=-1; elseif A*B<0 s=1; elseif A==0 s=sign(randn); elseif B==0 s=0; end elseif (x==1)&(b==0)&(fitness(i)>=best.fitness) delta=0.01*pi; if A*B>0 s=1; elseif A*B<0 s=-1; elseif A==0 s=0; elseif B==0 s=sign(randn); end end e=s*delta; % e为旋转角 U=[cos(e) -sin(e);sin(e) cos(e)]; % 量子旋转门 y=U*[A B]'; % y为更新后的量子位 chrom(2*i-1,j)=y(1); chrom(2*i,j)=y(2); end end