gqa.m

遗传量子算法求解背包问题Matlab源码

%gqa
n=input('please input population size n=:');%群体规模
g=input('please input max-generation g=:');%进化代数
for number=1:30
     clc
     t=0;
     initialize;%初始化
     observe;%观测染色体，将量子态转化为二进制的问题解形式
     repair;%修复染色体，防止问题解超过重量限制
     evaluate;%评估问题解的适应度值，并从当代群体中选择出一个最优的，保留
      store;%当代最优跟保持的最优值比较，保留其中一个最优的
      while t<g  %进化循环
            t=t+1;
            observe;
            repair;
            evaluate;
            update;%量子门更新，产生下一代的量子态染色体
            store;
      end     
end
for t=1:g
 maxgqa10(t)=mean(MAX(:,t));
 avegqa10(t)=mean(ave(:,t)); 
end
 save data/data maxgqa10 avegqa10
 for number=1:30
     plot(MAX(number,:));
 end
    plot(maxgqa10,'r-');
hold on;
plot(avegqa10);