(完整版)螢火蟲算法,matlab代碼

1、% 算法說(shuō)明：熒火蟲算法clc;%清屏clear all;%清除變量format long;%確定精度% 各參數(shù)初始化開始domx=-5.12,5.12;-5.12,5.12;% domx=-2.048,2.048;-2.048,2.048;%解空間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;%各參數(shù)初始化結(jié)束iter_max=200;% 分配空間開始%熒光素?fù)]發(fā)因子%適應(yīng)度提取比例%鄰域變化率%鄰域閥值(鄰域熒火蟲數(shù))%步長(zhǎng)%熒光素濃度%感知半徑%感知半徑 %決策半徑%最

2、大迭代次數(shù)m=size(domx,1);global n;n=50;address=zeros(n,m);value=zeros(n,1);li=zeros(n,1);rdi=zeros(n,1);%函數(shù)空間維數(shù)%種群規(guī)模%分配熒火蟲地址空間%分配適應(yīng)度存放空間%分配熒光素存放空間%分配熒火蟲決策半徑存放空間uu=zeros(1,iter_max);% 熒火蟲常量初始化開始for i=1:maddress(:,i)=(domx(i,1)+(domx(i,2)-domx(i,1)*rand(n,1); % 初始化地址%address(i,:)=5*rands(n,1);%隨機(jī)產(chǎn)生初螢火蟲所在位置

3、endf=fun(address);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*')hold on;grid on;z=-(x.A2-10*cos(2*pi.*x)+10+y.A2-10*cos(2*pi.*y)+10); %Rastrigin'函數(shù)% z=-(x.A2+y.A2); % 目標(biāo)函數(shù) J2% z=-(20+xA2-10*cos(2*pi.*x)+yA2-10*cos(2*pi.*y); % 目標(biāo)函數(shù)

4、 J1% z=-(-20*exp(-0.2*sqrt(x.A2+y.A2)/2)-exp(cos(2*pi*x)+cos(2*pi*y)/2)+20+exp(1);%目標(biāo)函數(shù)% z=-(0.5+sin(sqrt(x.A2 + y.A2).A2-0.5)./(1+0.001*(x.A2 + y.A2).A2);% 目標(biāo)函數(shù)J3mesh(x,y,z) xlabel('x 軸 ');ylabel('y 軸 ');zlabel('z 軸 ');title(' 螢火蟲初始分布圖');li(:,1)=lo;%熒光素初始值rdi(:,1)=r0

5、;%決策半徑初始值t=1;%迭代累計(jì)量%熒火蟲常量初始化結(jié)束% 迭代開始while(t<=iter_max)li=(1-rho)*li+ gamma * fun(address);%.更新熒光素值li%各熒火蟲移動(dòng)過(guò)程開始for i=1:nOrd_number=;%存放熒火蟲序號(hào)for j=1:nif (norm(address(j,:)-address(i,:)<rdi(i) && (li(i,1)<li(j,1) % 決策半徑內(nèi)找更優(yōu)點(diǎn)Ord_number(numel(Ord_number)+1)=j;endendj 位置%計(jì)算Ord_number 各元素

6、被選擇概率，確定if isempty(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./Den;Pij=cumsum(Pij);Pij=Pij./Pij(end);Pos=find(rand<Pij);j=Ord_number(Pos(1);%先判斷Ord_number 個(gè)數(shù)不為空%選出Ord_number 熒光素%Ord_number 熒光素和%分子(lj-li)%分母%計(jì)算各元素被選擇概率%累計(jì)

7、%歸一化%確定位置%確定j 的位置%熒火蟲i 向位置 j 移動(dòng)address(i,:)=address(i,:)+ s*(address(j,:)-address(i,:)/norm(address(j,:)-address(i,:); address(i,:)=range1(address(i,:),domx);% 限制范圍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 迭代結(jié)束g=;g=fun

8、(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*');hold on;z=-(x.A2-10*cos(2*pi.*x)+10+y.A2-10*cos(2*pi.*y)+10); %Rastrigin'函數(shù)% z=-(x.A2+y.A2);% 目標(biāo)函數(shù) J2% z=-(20+x.A2-10*cos(2*pi*x)+y.A2-10*cos(2*pi*y); %目標(biāo)函數(shù) J1%目% z=-(-20*ex

9、p(-0.2*sqrt(xA2+yA2)/2)-exp(cos(2*pi*x)+cos(2*pi*y)/2)+20+exp(1);標(biāo)函數(shù)% z=-(0.5+sin(sqrt(x.A2 + y.A2).A2-0.5)./(1+0.001*(x.A2 + y.A2).A2);% 目標(biāo)函數(shù)J3grid on;mesh(x,y,z);title(' 螢火蟲運(yùn)動(dòng)軌跡');xlabel('x 軸 ');ylabel('y 軸 ');zlabel('z 軸 '); hold off;endfigure('name','

10、收斂曲線');plot(uu)title('螢火蟲最優(yōu)值曲線''終止次數(shù)='num2str(iter_max);xlabel(' 迭代次數(shù)');ylabel('最優(yōu)值')； % 輸出最優(yōu)結(jié)果value=fun(address);disp('最優(yōu)值為：)num=find(value=max(value);% 最大值序號(hào) MaxValue=-max(value)disp('最優(yōu)解為：)BestAddress = address(num,:)%目標(biāo)函數(shù)function y = fun(x)% y=x(:,1)A2+x(:,2)A2;% y=0.5+sin(sqrt(x(:,1).A2 + x(:,2).A2).A2-0.5)./(1+0.001*(x(:,1).A2 + x(:,2).A2).A2;%目標(biāo)函數(shù) J3% y=20+x(:,1)A2-10*cos(2*pi*x(:,1)+x(:,2)A2-10*cos(2*pi*x(:,2); % 目標(biāo)函數(shù) J1% y=-20*exp(-0.2*sqrt(x(: