(完整版)萤火虫算法,matlab代码

1、rho=0.4;gamma=0.6;beta=0.08;nt=5;s=0.03;lo=5;rs=5.12;r0=5.12;%rs=2.048;%r0=2.048;%各参数初始化结束iter_max=200;%分配空间开始m=size(domx,1);globaln;n=50;address=zeros(n,m);value=zeros(n,1);li=zeros(n,1);rdi=zeros(n,1);%算法说明：荧火虫算法clc;%清屏clearall;%清除变量formatlong;%确定精度%各参数初始化开始domx=-5.12,5.12;-5.12,5.12;%domx=-2.048,

2、2.048;-2.048,2.048;%解空间%荧光素挥发因子%适应度提取比例%邻域变化率%邻域阀值(邻域荧火虫数)%步长%荧光素浓度%感知半径%感知半径%决策半径%最大迭代次数%函数空间维数%种群规模%分配荧火虫地址空间%分配适应度存放空间%分配荧光素存放空间%分配荧火虫决策半径存放空间uu=zeros(1,iter_max);%荧火虫常量初始化开始fori=1:maddress(:,i)=(domx(i,1)+(domx(i,2)-domx(i,1)*rand(n,1);%初始化地址%address(i,:)=5*rands(n,1);%随机产生初萤火虫所在位置endf=fun(addre

3、ss);x=-5.12:.1:5.12;%x=-2.048:.05:2.048;x,y=meshgrid(x);figure(1);plot3(address(:,1),address(:,2),f,k*)holdon;gridon;z=-(x.A2-10*cos(2*pi.*x)+10+y.A2-10*cos(2*pi.*y)+10);%Rastrigin函数%z=-(x.A2+y.A2);%目标函数J2%z=-(20+x.A2-10*cos(2*pi.*x)+y.A2-10*cos(2*pi.*y);%目标函数J1%目%z=-(-20*exp(-0.2*sqrt(x42+y.A2)/2)-

4、exp(cos(2*pi*x)+cos(2*pi*y)/2)+20+exp(1);标函数%z=-(0.5+sin(sqrt(x42+y.A2).A2-0.5)./(1+0.001*(x.A2+y.A2).A2);%目标函数J3mesh(x,y,z)xlabel(x轴);ylabel(y轴);zlabel(z轴);title(萤火虫初始分布图);li(:,1)=lo;%荧光素初始值rdi(:,1)=r0;%决策半径初始值t=1;%迭代累计量%荧火虫常量初始化结束%迭代开始while(t=iter_max)li=(1-rho)*li+gamma*fun(address);%.更新荧光素值li%各荧

5、火虫移动过程开始fori=1:nOrd_number=;%存放荧火虫序号forj=1:nif(norm(address(j,:)-address(i,:)rdi(i)&(li(i,1)li(j,1)%决策半径内找更优点Ord_number(numel(Ord_number)+1)=j;endend%计算Ord_number各元素被选择概率ifisempty(Ord_number)Ord_num_li=li(Ord_number,1);Sum_Ord_li=sum(Ord_num_li);Mol=Ord_num_li-li(i,1);Den=Sum_Ord_li-li(i,1);Pij=Mol.

6、/Den;Pij=cumsum(Pij);Pij=Pij./Pij(end);Pos=find(randPij);j=Ord_number(Pos(1);确定j位置%先判断Ord_number个数不为空%选出Ord_number荧光素%Ord_number荧光素和%分子(lj-li)%分母%计算各元素被选择概率%累计%归一化%确定位置%确定j的位置%荧火虫i向位置j移动address(i,:)=address(i,:)+s*(address(j,:)-address(i,:)/norm(address(j,:)-address(i,:);address(i,:)=range1(address(

7、i,:),domx);%P艮制范围end%更新决策半径rdi(i)=rdi(i)+beta*(nt-length(Ord_number);rdi(i)=min(rs,max(0,rdi(i);uu(t)=-max(fun(address);endt=t+1;%iter_max迭代结束g=;g=fun(address);x=-5.12:.1:5.12;%x=-2.048:.05:2.048;x,y=meshgrid(x);figure(2);plot3(address(:,1),address(:,2),g,k*);holdon;z=-(x.A2-10*cos(2*pi.*x)+10+y.A2-

8、10*cos(2*pi.*y)+10);%Rastrigin函数%z=-(x.A2+y.A2);%目标函数J2%目%z=-(20+x.A2-10*cos(2*pi*x)+y.A2-10*cos(2*pi*y);%目标函数J1%z=-(-20*exp(-0.2*sqrt(x42+y.A2)/2)-exp(cos(2*pi*x)+cos(2*pi*y)/2)+20+exp(1);标函数%z=-(0.5+sin(sqrt(x42+y.A2).A2-0.5)./(1+0.001*(x.A2+y.A2).A2);%目标函数J3gridon;mesh(x,y,z);title(萤火虫运动轨迹);xlabe

9、l(x轴);ylabel(y轴);zlabel(z轴);holdoff;endfigure(name,收敛曲线);plot(uu)title(萤火虫最优值曲线终止次数=num2str(iter_max);xlabel(迭代次数);ylabel(最优值);%输出最优结果value=fun(address);disp(最优值为：)num=find(value=max(value);%最大值序号MaxValue=-max(value)disp(最优解为：)BestAddress=address(num,:)%目标函数functiony=fun(x)%y=x(:,1).A2+x(:,2).A2;%y=0.5+sin(sqrt(x(:,1)42+x(:,2).A2).A2-0.5)./(1+0.001*(x(:,1).A2+x(:,2).A2).A2;%目标函数J3%y=20+x(:,1).A2-10*cos(2*pi*x(:,1)+x(:,2)42-10*cos(2*pi*x(:,2);%目标函数J1%y=-20*exp(-0.2*sqrt(x(:,1).A2+x(:,