untitled222bb.asv

这个是我的毕业设计做的我觉得对大家还是有写用的

clear all;
close all;
Size=80;   
G=1000;     
CodeL=21;
umax=10.;
umin=-10;
E=round(rand(Size,CodeL));    
for k=1:1:G
time(k)=k;
for s=1:1:Size
m=E(s,:);
y=0;
m1=m(1:1:CodeL);
for j=1:1:CodeL
   y=y+m1(j)*2^(j-1);
end
x=(umax-umin)*y/2097151+umin;
F(s)=(0.5+(sin(sum(x.^2)).^(1/2).*sin(sum(x.^2)).^(1/2)-0.5)/(1.0+0.01*sum(x.^2)).^2);;
end
fi=F;                          
[Oderfi,Indexfi]=sort(fi);     
Minfi=Oderfi(1);           
MinS=E(Indexfi(1),:);     
bfi(k)=Minfi;
   fi_sum=sum(fi);
   fi_Size=(Oderfi/fi_sum)*Size;
   fi_S=floor(fi_Size);       
   kk=1;
   for i=1:1:Size
      for j=1:1:fi_S(i)        
       TempE(kk,:)=E(Indexfi(i),:);  
         kk=kk+1;              
      end
   end
pc=0.60;
n=ceil(20*rand);
for i=1:2:(Size-1)
    temp=rand;
    if pc>temp                  
    for j=n:1:20
        TempE(i,j)=E(i+1,j);
        TempE(i+1,j)=E(i,j);
    end
    end
end
TempE(Size,:)=MinS;
E=TempE;
pm=0.1;   
 for i=1:1:Size
      for j=1:1:CodeL
         temp=rand;
         if pm>temp               
            if TempE(i,j)==0
               TempE(i,j)=1;
            else
               TempE(i,j)=0;
            end
        end
      end
   end
TempE(Size,:)=MinS;
E=TempE;
end

x

Min_Value=Minfi
MinS