2)).^2);%目标函数 J3Рgrid on;Рmesh(x,y,z);Рtitle('萤火虫运动轨迹');Рxlabel('x轴');Рylabel('y轴');Рzlabel('z轴');Рhold off;РendРfigure('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,:) Р%目标函数Рfunction y = fun(x)Р% y=x(:,1).^2+x(:,2).^2; Р% y=0.5+sin(sqrt(x(:,1).^2 + x(:,2).^2).^2-0.5)./(1+0.001*(x(:,1).^2 + x(:,2).^2)).^2;%目标函数 J3Р% y=20+x(:,1).^2-10*cos(2*pi*x(:,1))+x(:,2).^2-10*cos(2*pi*x(:,2)); %目标函数 J1 Р% y=-20*exp(-0.2*sqrt((x(:,1).^2+x(:,2).^2)/2))...Р% -exp((cos(2*pi*x(:,1))+cos(2*pi*x(:,2)))/2)+20+exp(1); %%目标函数Рy=x(:,1).^2-10*cos(2*pi.*x(:,1))+10+x(:,2).^2-10*cos(2*pi.*x(:,2))+10; %Rastrigin'函数Р y=-y;