gmjafsa.m

FSA算法的MATLAB实现

%人工鱼群算法
format long
Visual=1.5;
Step=0.1;
N=50;
Try_number=50;
a1=3;
b1=6;
a2=3;
b2=5;
d=[];
h=0.01;%允许的误差，循环结束条件
Friend_number=50;
k=0;
m=10;
X1=rand(N,1)*(b1-a1)+a1;
X2=rand(N,1)*(b2-a2)+a2;
X=[X1 X2];%人工鱼数量
for i=1:N
wwww=[X(i,1),X(i,2)];
d(i)=maxf(wwww);
end
[w,i]=min(d);
maxX=[X(i,1),X(i,2)];%初始公告板记录
maxY=w;%初始公告板记录
figurex=[];
figurey=[];
figurez=[];
figurex(numel(figurex)+1)=maxX(1);
figurey(numel(figurey)+1)=maxX(2);
figurez(numel(figurez)+1)=maxY;
kkk=0;
while(~(maxY<h))
    for i=1:N    
    XX=[X(i,1),X(i,2)];%人工鱼当前状态Xi
    nf=0;
    Xc=0;
    for j=1:N  %聚群行为开始
       XXX=[X(j,1),X(j,2)];
       if(norm(XXX-XX)<Visual)
       nf=nf+1;
       Xc=Xc+XXX;
       end
    end
    Xc=Xc/nf;
   if((maxf(Xc))<maxf(XX))
        XXnext1=Xc+rand*Step*(XX-Xc)/norm(XX-Xc);   
        if(XXnext1(1)>b1)
        XXnext1(1)=b1;   
        end
        if(XXnext1(1)<a1)
        XXnext1(1)=a1;   
        end
        if(XXnext1(2)>b2)
        XXnext1(2)=b2;   
        end
        if(XXnext1(2)<a2)
        XXnext1(2)=a2;   
        end
   else
       XXnext1=gmjprey(XX,Try_number,Visual,Step);
       if(XXnext1(1)>b1)
        XXnext1(1)=b1;   
        end
        if(XXnext1(1)<a1)
        XXnext1(1)=a1;   
        end
        if(XXnext1(2)>b2)
        XXnext1(2)=b2;   
        end
        if(XXnext1(2)<a2)
        XXnext1(2)=a2;   
        end
   end%聚群行为结束
   %maxX=XX;%追尾行为开始
   %maxY=maxf(XX);
   for j=1:Friend_number
     XXX=[X(j,1),X(j,2)];
     if(norm(XXX-XX)<Visual & maxf(XXX)<maxY)
         maxX=XXX;
         maxY=maxf(XXX);
     end
   end
   if((maxY)<maxf(XX))
       XXnext2=maxX+rand*Step*(XX-maxX)/norm(XX-maxX);  
       if(XXnext2(1)>b1)
        XXnext2(1)=b1;   
        end
        if(XXnext2(1)<a1)
        XXnext2(1)=a1;   
        end
        if(XXnext2(2)>b2)
        XXnext2(2)=b2;   
        end
        if(XXnext2(2)<a2)
        XXnext2(2)=a2;   
        end
   else
      XXnext2 =gmjprey(XX,Try_number,Visual,Step);
      if(XXnext2(1)>b1)
        XXnext2(1)=b1;   
        end
        if(XXnext2(1)<a1)
        XXnext2(1)=a1;   
        end
        if(XXnext2(2)>b2)
        XXnext2(2)=b2;   
        end
        if(XXnext2(2)<a2)
        XXnext2(2)=a2;   
        end
   end%追尾行为结束
   if(maxf(XXnext1)<maxf(XXnext2))
     X(i,1)=XXnext1(1);
     X(i,2)=XXnext1(2);
   else
     X(i,1)=XXnext2(1);
     X(i,2)=XXnext2(2);
   end
end %一次迭代结束  
for i=1:N
    XXXX=[X(i,1),X(i,2)];
       if maxf(XXXX)<maxY
            maxY=maxf(XXXX);
            maxX=XXXX;
            figurex(numel(figurex)+1)=maxX(1);
            figurey(numel(figurey)+1)=maxX(2);
            figurez(numel(figurez)+1)=maxY;
       end
end
    k=k+1; %进入下一次迭代
end   
maxX
maxY
plot3(figurex,figurey,figurez)