遺傳算法源程序

遺傳算法源程序functionGenetic(StrAimFunc)%Thisissimplegeneticalgorithm(SGA)%Inthisfunction,itfulfilsgeneticalgorithm%Thesecanbemodifiedasyoulikemaxgen=200;sizepop=100;AimFunc=StrAimFunc;countingfitnessfselect='roulette';'tournament';'roulette'fcode='grey';'float';'grey';'binary'pcross=[0.6];crossover,between0and1fcross='uniform';'float';'simple';'uniform'pmutation=[0.2];mutation,between0and1fmutation='simple';'float';'simple';lenchrom=[1010101010];variblebound=[01;...01;...01;...01;...01];%%maximumgeneration%sizeofpopulation%thisisfunctionof%methodofselect%youcanchoose%methodofcoding%youcanchoose%probablityof%methodofcrossover%youcanchoose%probabilityof%methodofmutation%youcanchoose%lengthofbitofeveryindividuals=struct('fitness',zeros(1,sizepop),...%'value',zeros(1,sizepo%structure%structureofpopulation%averagefitnessof'chrom',[]);avgfitness=[];population%bestfitnessofpopulation%chromosomeof%bestfitnessofpopulation%chromosomeofbest%producenewpopulationatrandomindividuals.chrom(i,:)=Code(lenchrom,fcode,bound);x=Decode(lenchrom,bound,individuals.chrom(i,:),fcode);%caculatefitness%修改處%individuals.fitness(i)=AimFunc(x,heli,id_data,0);individuals.fitness(i)=feval(AimFunc,x);end%findminimumvaluewhichisbest[bestfitnessbestindex]=min(individuals.fitness);bestchrom=individuals.chrom(bestindex,:);avgfitness=sum(individuals.fitness)/sizepop;%recordaverageandbestfitnessofeverygenerationtrace=[avgfitnessbestfitness];%evolutionbeginfori=1:maxgen%selectionindividuals=Select(individuals,sizepop,fselect);avgfitness=sum(individuals.fitness)/sizepop;%crossoverindividuals.chrom=Cross(pcross,lenchrom,individuals.chrom,...sizepop,fcross,[imaxgen]);%mutationindividuals.chrom=Mutation(pmutation,lenchrom,individuals.chrom,.sizepop,fmutation,[imaxgen],bound);%calculatefitnessforj=1:sizepopx=Decode(lenchrom,bound,individuals.chrom(j,:),fcode);%修改處%individuals.fitness(j)=AimFunc(x,heli,id_data,0);individuals.fitness(j)=feval(AimFunc,x);end%substitutechromosomeofworestfitness%findminimumvaluewhichisbest[newbestfitness,newbestindex]=min(individuals.fitness);[worestfitness,worestindex]=max(individuals.fitness);%substitutechromosomeofworestfitnessifbestfitness>newbestfitnessbestfitness=newbestfitness;bestchrom=individuals.chrom(newbestindex,:);endindividuals.chrom(worestindex,:)=bestchrom;individuals.fitness(worestindex)=bestfitness;avgfitness=sum(individuals.fitness)/sizepop;trace=[trace;avgfitnessbestfitness];endend%drawfitnessofeverygenerationhfig=findobj('Tag','trace');%Seeifitisopenifishandle(hfig)figure(hfig);elsehfig=figure('Tag','trace');endfigure(hfig);[rc]=size(trace);plot([1:r]',trace(:,1),'r-',[1:r]',trace(:,2),'b--');title(['適應(yīng)度曲線’'終止代數(shù)='num2str(maxgen)]);xlabel('進(jìn)化代數(shù)');ylabel(適應(yīng)度');legend('平均適應(yīng)度','最佳適應(yīng)度');disp(適應(yīng)度變量');x=Decode(lenchrom,bound,bestchrom,fcode);%showincommandwindowdisp([bestfitnessx]);%functionret=Code(lenchrom,opts,bound)%Inthisfunction,itconvertsasetvariblesintoachromosome%lenchrominput:lengthofchromosome%optsinput:tagofcodingmethod%boundinput:boundaryofvaribles%retoutput:chromosomeswitchoptscase'binary'%binarycodingpick=rand(1,sum(lenchrom));bits=ceil(pick-0.5);temp=sum(lenchrom)-1:-1:0;ret=sum(bits.*(2.Atemp));case'grey'%greycodingpick=rand(1,sum(lenchrom));bits=ceil(pick-0.5);greybits=bits;fori=2:length(greybits)greybits(i)=bitxor(bits(i-1),bits(i));endtemp=sum(lenchrom)-1:-1:0;ret=sum(greybits.*(2.Atemp));case'float'%floatcodingpick=rand(1,length(lenchrom));ret=bound(:,1)'+(bound(:,2)-bound(:,1))'.*pick;end%functionret=Decode(lenchrom,bound,code,opts)%Inthisfunction,itdecodechromosome%lenchrominput:lengthofchromosome%optsinput:tagofcodingmethod%boundinput:boundaryofvaribles

%retoutput:valueofvariblesswitchoptscase'binary'%binarycodingfori=length(lenchrom):-1:1data(i)=bitand(code,2Alenchrom(i)-1);code=(code-data(i))/(2Alenchrom(i));endret=bound(:,1)'+data./(2.Alenchrom-1).*(bound(:,2)-bound(:,1))';case'grey'%greycodingfori=sum(lenchrom):-1:2code=bitset(code,i-1,bitxor(bitget(code,i),bitget(code,i-1)));endfori=length(lenchrom):-1:1data(i)=bitand(code,2Alenchrom(i)-1);code=(code-data(i))/(2Alenchrom(i));endret=bound(:,1)'+data./(2.Alenchrom-1).*(bound(:,2)-bound(:,1))';case'float'%floatcodingret=code;end%functioninputinputinputinputinput:probabilityofmutation:lengthofachromosome:setofallchromosomes:sizeofpopulation:tagforchoosingmethodofcrossoverinputinputinputinputinput:probabilityofmutation:lengthofachromosome:setofallchromosomes:sizeofpopulation:tagforchoosingmethodofcrossoverinput:currentserialnumberofgenerationandgemerationoutput:newsetofchromosomecase'simple'%mutationatsinglepositionfori=1:sizepop%selectchildatrandompick=rand;whilepick==0pick=rand;endindex=ceil(pick*sizepop);pick=rand;ifpick>pmutationcontinue;end%mutationpositionpick=rand;whilepick==0pick=rand;endpos=ceil(pick*sum(lenchrom));v=bitget(chrom(index),pos);v=~v;chrom1=bitset(chrom(index),pos,v);endret=chrom;case'float'%multiplepositionmutationfori=1:sizepop%selectchildatrandompick=rand;whilepick==0pick=rand;endindex=ceil(pick*sizepop);pick=rand;ifpick>pmutationcontinue;end%mutationpositionpick=rand;whilepick==0pick=rand;endpos=ceil(pick*sum(lenchrom));v=chrom(i,pos);v1=v-bound(pos,1);v2=bound(pos,2)-v;pick=rand;ifpick>0.5delta=v2*(1-pickA((1-pop(1)/pop(2))A2));chrom(i,pos)=v+delta;elsedelta=v1*(1-pickA((1-pop(1)/pop(2))A2));chrom(i,pos)=v-delta;endendret=chrom;end%functionret=select(individuals,sizepop,opts)%Inthisfunction,itfulfilsaselectionamongchromosomes%individualsinput:setofindividuals%sizepopinput:sizeofpopulation%optsinput:tagforchoosingmethodofselection%retoutput:newsetofindividualsswitchoptscase'roulette'%roulettewheelmodelsumfitness=sum(1./individuals.fitness);sumf=(1./individuals.fitness)./sumfitness;index=[];fori=1:sizepoppick=rand;whilepick==0pick=r