%gqa
n=input('please input population size n=:');%群體規(guī)模
g=input('please input max-generation g=:');%進(jìn)化代數(shù)
for number=1:30
     clc
     t=0;
     initialize;%初始化
     observe;%觀測(cè)染色體，將量子態(tài)轉(zhuǎn)化為二進(jìn)制的問(wèn)題解形式
     repair;%修復(fù)染色體，防止問(wèn)題解超過(guò)重量限制
     evaluate;%評(píng)估問(wèn)題解的適應(yīng)度值，并從當(dāng)代群體中選擇出一個(gè)最優(yōu)的，保留
      store;%當(dāng)代最優(yōu)跟保持的最優(yōu)值比較，保留其中一個(gè)最優(yōu)的
      while t<g  %進(jìn)化循環(huán)
            t=t+1;
            observe;
            repair;
            evaluate;
            update;%量子門(mén)更新，產(chǎn)生下一代的量子態(tài)染色體
            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);