人工鱼群算法源代码

人工鱼群算法源代码%人工鱼群算法；formatlong；Visual=2.5;；Step=0.3;；N=50;；Try_number=50;；a1=-10;；b1=10;；a2=-10;；b2=10;；d=;；h=1e-1;；Friend_number=50;；k=0;；m=50;；X1=rand(N,1)*(b1-a1)+a1;；X2=rand(N,1)*(b2-a2)%人工鱼群算法format longVisual=2.5;Step=0.3;N=50;Try_number=50;a1=-10;b1=10;a2=-10;b2=10;d=;h=1e-1;Friend_number=50;k=0;m=50;X1=rand(N,1)*(b1-a1)+a1;X2=rand(N,1)*(b2-a2)+a2;X=X1 X2;%人工鱼数量for i=1:Nwwww=X(i,1),X(i,2);d(i)=maxf(wwww);endw,i=max(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(km)for i=1:NXX=X(i,1),X(i,2);%人工鱼当前状态Xinf=0;Xc=0;for j=1:N %聚群行为开始XXX=X(j,1),X(j,2);if(norm(XXX-XX)maxf(XX)XXnext1=XX+rand*Step*(Xc-XX)/norm(Xc-XX); if(XXnext1(1)b1)XXnext1(1)=b1;endif(XXnext1(1)b2)XXnext1(2)=b2;endif(XXnext1(2)b1)XXnext1(1)=b1;endif(XXnext1(1)b2)XXnext1(2)=b2;endif(XXnext1(2)a2)XXnext1(2)=a2;endend%聚群行为结束%maxX=XX;%追尾行为开始%maxY=maxf(XX);for j=1:Friend_numberXXX=X(j,1),X(j,2);if(norm(XXX-XX)maxY) maxX=XXX;maxY=maxf(XXX);endendif(maxY)maxf(XX)XXnext2=XX+rand*Step*(maxX-XX)/norm(maxX-XX); if(XXnext2(1)b1)XXnext2(1)=b1;endif(XXnext2(1)b2)XXnext2(2)=b2;endif(XXnext2(2)b1)XXnext2(1)=b1;endif(XXnext2(1)b2)XXnext2(2)=b2;endif(XXnext2(2)maxf(XXnext2)X(i,1)=XXnext1(1);X(i,2)=XXnext1(2);elseX(i,1)=XXnext2(1);X(i,2)=XXnext2(2);endend %一次迭代结束for i=1:NXXXX=X(i,1),X(i,2);if maxf(XXXX)maxYmaxY=maxf(XXXX);maxX=XXXX;figurex(numel(figurex)+1)=maxX(1); figurey(numel(figurey)+1)=maxX(2); figurez(numel(figurez)+1)=maxY;endendk=k+1; %进入下一次迭代endmaxXmaxYplot3(figurex,figurey,figurez)function XXnext = gmjprey(XX,Try_number,Visual,Step) pp=0;for j=1:Try_numberXXj=XX+rand*Step*Visual;if(maxf(XX)maxf(XXj)XXnext=XX+rand*Step*(XXj-XX)/norm(XXj-XX); pp=1;breakendendif(pp)XXnext=XX+rand*Step;endfunction y = maxf(QQ)y=(sin(QQ(1)/QQ(1)*(sin(QQ(2)/QQ(2);人工鱼群算法的仿真程序-matlab23tic；figure(1);holdon；ezplot(x*sin(10*pi*；%参数设置；fishnum=50;%生成50只人工鱼；MAXGEN=50;%最多迭代次数；try_number=100;%最多试探次数；visual=1;%感知距离；delta=0.618;%拥挤度因子；step=0.1;%步长；%初始化鱼群；lb_ub=-1,ticfigure(1);hold onezplot(x*sin(10*pi*x)+2,-1,2);% 参数设置fishnum=50; %生成50只人工鱼MAXGEN=50; %最多迭代次数try_number=100;%最多试探次数visual=1; %感知距离delta=0.618; %拥挤度因子step=0.1; %步长% 初始化鱼群lb_ub=-1,2,1;X=AF_init(fishnum,lb_ub);LBUB=;fori=1:size(lb_ub,1)LBUB=LBUB;repmat(lb_ub(i,1:2),lb_ub(i,3),1);endgen=1;BestY=-1*ones(1,MAXGEN); %每步中最优的函数值BestX=-1*ones(1,MAXGEN); %每步中最优的自变量besty=-100; %最优函数值Y=AF_foodconsistence(X);while genYi2X(:,i)=Xi1;Y(1,i)=Yi1;elseX(:,i)=Xi2;Y(1,i)=Yi2;endendYmax,index=max(Y);figure(1);plot(X(1,index),Ymax,.,color,gen/MAXGEN,0,0)ifYmaxbestybesty=Ymax;bestx=X(:,index);BestY(gen)=Ymax;BestX(:,gen)=X(:,index); elseBestY(gen)=BestY(gen-1);BestX(:,gen)=BestX(:,gen-1); endgen=gen+1;endplot(bestx(1),besty,ro,MarkerSize,100) xlabel(x)ylabel(y)title(鱼群算法迭代