【正文】 v2=chrom(index(2),pos)。 chrom(index(1),pos)=pick*v2+(1pick)*v1。 chrom(index(2),pos)=pick*v1+(1pick)*v2。 end ret=chrom。endfunction ret=Decode(lenchrom,bound,code,opts)% In this function ,it decode chromosome % lenchrom input : length of chromosome% opts input : tag of coding method% bound input : boundary of varibles% ret output: value of variblesswitch opts case 39。binary39。 % binary coding for i=length(lenchrom):1:1 data(i)=bitand(code,2^lenchrom(i)1)。 code=(codedata(i))/(2^lenchrom(i))。 end ret=bound(:,1)39。+data./(2.^lenchrom1).*(bound(:,2)bound(:,1))39。 case 39。grey39。 % grey coding for i=sum(lenchrom):1:2 code=bitset(code,i1,bitxor(bitget(code,i),bitget(code,i1)))。 end for i=length(lenchrom):1:1 data(i)=bitand(code,2^lenchrom(i)1)。 code=(codedata(i))/(2^lenchrom(i))。 end ret=bound(:,1)39。+data./(2.^lenchrom1).*(bound(:,2)bound(:,1))39。 case 39。float39。 % float coding ret=code。endfunction y=f1(t)y=*erfc((3000030*t)./sqrt(*t))。function ret=Mutation(pmutation,lenchrom,chrom,sizepop,opts,pop,bound)% In this function,it fulfils a mutation among chromosomes% pcorss input : probability of mutation% lenchrom input : length of a chromosome% chrom input : set of all chromosomes% sizepop input : size of population% opts input : tag for choosing method of crossover% pop input : current serial number of generation and maximum gemeration% ret output : new set of chromosomeswitch opts case 39。simple39。 % mutation at single position for i=1:sizepop % select child at random pick=rand。 while pick==0 pick=rand。 end index=ceil(pick*sizepop)。 pick=rand。 if pickpmutation continue。 end % mutation position pick=rand。 while pick==0 pick=rand。 end pos=ceil(pick*sum(lenchrom))。 v=bitget(chrom(index),pos)。 v=~v。 chrom1=bitset(chrom(index),pos,v)。 end ret=chrom。 case 39。float39。 % multiple position mutation for i=1:sizepop % select child at random pick=rand。 while pick==0 pick=rand。 end index=ceil(pick*sizepop)。 pick=rand。 if pickpmutation continue。 end % mutation position pick=rand。 while pick==0 pick=rand。 end pos=ceil(pick*sum(lenchrom))。 v=chrom(i,pos)。 v1=vbound(pos,1)。 v2=bound(pos,2)v。 pick=rand。 if pick delta=v2*(1pick^((1pop(1)/pop(2))^2))。 chrom(i,pos)=v+delta。 else delta=v1*(1pick^((1pop(1)/pop(2))^2))。 chrom(i,pos)=vdelta。 end end ret=chrom。endfunction ret=select(individuals,sizepop,opts)% In this function,it fulfils a selection among chromosomes% individuals input : set of individuals% sizepop input : size of population% opts input : tag for choosing method of selection% ret output : new set of individualsswitch opts case 39。roulette39。 % roulette wheel model sumfitness=sum(1./)。 sumf=(1./)./sumfitness。 index=[]。 for i=1:sizepop pick=rand。 while pick==0 pick=rand。 end for i=1:sizepop pick=picksumf(i)。 if pick0 index=[index i]。 break。 end end end =(index,:)。 =(index)。 ret=individuals。case 39。tournament39。 % tourment model allindex=[]。 for i=1:sizepop pick=rand(1,2)。 while prod(pick)==0 pick=rand(1,2)。 end index=ceil(pick.*sizepop)。 if (index(1))=(index(2)) allindex=[allindex index(1)]。 else allindex=[allindex index(2)]。 end end =(allindex,:)。 =(allindex)。 ret=individuals。 endz=zeros(1)。z=quadl(39。f139。,1000,9800)。disp