文章目录

  • 一、理论基础
  • 二、GSO算法
    • 1、萤火虫的部署(初始化)
    • 2、荧光素更新阶段
    • 3、移动概率计算阶段
    • 4、位置更新阶段
    • 5、邻域范围更新阶段
  • 三、MATLAB程序实现
    • 1、问题描述
    • 2、适应度函数
    • 3、主函数
    • 3、结果显示
    • 4、绘图
  • 四、参考文献
一、理论基础

2005年，印度学者K.N.Krishnanand和D.Ghose在IEEE群体智能会议上提出了一种新的群智能优化算法，人工萤火虫群优化(Glowworm Swarm Optimization, GSO)算法。2009年，剑桥学者Xin-She Yang根据自然界中萤火虫的发光行为提出萤火虫算法(Firefly Algorithm, FA)。自这两种萤火虫算法提出以来，各国学者对这两种算法进行了研究、改进和应用。经过几年的发展，在连续空间的寻优过程和一些生产调度方面萤火虫算法具有良好的应用前景。
GSO和FA有相似的方面，但在具体实现方面有一定差异。本文具体介绍GSO算法。

二、GSO算法

算法基本思想描述如下：在群体中，每个萤火虫个体被随机分布在目标函数定义的空间中，初始阶段，所有的萤火虫都具有相同的荧光素值和动态决策半径。其中，每个萤火虫个体根据来自动态决策半径内所有邻居萤火虫信号的强弱来决定其移动的方向。萤火虫的动态决策半径会随着在它范围内萤火虫个体的数目而变化，每个萤火虫的荧光素也会随着决策半径内萤火虫个体的数目而改变。萤火虫群优化算法是无记忆的，无需目标函数的全局信息和梯度信息，具有计算速度快，调节参数少，易于实现等特点。萤火虫进化过程中，每次迭代都由萤火虫的部署(初始化)、荧光素更新阶段、移动概率计算阶段、位置更新阶段、邻域范围更新阶段五个部分组成，现分别介绍如下：

1、萤火虫的部署(初始化)

人工萤火虫优化(GSO)算法_优化算法

2、荧光素更新阶段

人工萤火虫优化(GSO)算法_优化算法_02

3、移动概率计算阶段

人工萤火虫优化(GSO)算法_优化算法_03

4、位置更新阶段

人工萤火虫优化(GSO)算法_优化算法_04

5、邻域范围更新阶段

人工萤火虫优化(GSO)算法_优化算法_05

三、MATLAB程序实现

1、问题描述

Rosenbrock函数的三维立体图如图1所示。人工萤火虫优化(GSO)算法_优化算法_06

图1 Rosenbrock函数的三维立体图形

求解Rosenbrock函数在R = [ − 3 , 3 ] × [ − 3 , 3 ] \R=[-3,3]×[-3,3]R=[−3,3]×[−3,3]上的最大值。

2、适应度函数

适应度函数fun.m程序如下：

function z = fun(x)   
%% 适应度函数
z = 100.*(x(:, 1).^2-x(:, 2)).^2+(1-x(:, 1)).^2;
  • 1.
  • 2.
  • 3.

3、主函数

%% 清空环境变量
clc;
clear;

%% 初始化参数
domx = [-3, 3; -3, 3];       % 定义域
rho = 0.9;                   % 荧光素挥发因子
gamma = 0.1;                 % 适应度提取比例
beta = 0.58;                 % 邻域变化率
nt = 6;                      % 邻域阀值(邻域萤火虫数)
s = 0.03;                    % 步长
iot0 = 400;                  % 荧光素初始浓度
rs = 3;                      % 感知半径
r0 = 3;                      % 决策半径
m = size(domx, 1);           % 函数空间维数
n = 50;                      % 萤火虫数量
gaddress = zeros(n, m);      % 分配萤火虫地址空间
gvalue = zeros(n, 1);        % 分配适应度存放空间
ioti = zeros(n, 1);          % 分配荧光素存放空间
rdi = zeros(n, 1);           % 分配萤火虫决策半径存放空间

%% 萤火虫常量初始化
% 初始化地址
for i = 1:m
    gaddress(:, i) = domx(i, 1)+(domx(i, 2)-domx(i, 1))*rand(n, 1);
end
% 初始化荧光素浓度
ioti(:, 1) = iot0;
% 初始化决策半径
rdi(:, 1) = r0;
iter_max = 500;            % 最大迭代次数
t = 1;                     % 迭代计数器
yy = zeros(iter_max, 1);   % 各代最优解

%% 迭代寻优
while t <= iter_max
    % 更新荧光素浓度
    ioti = (1-rho)*ioti+gamma*fun(gaddress);
    % 各萤火虫移动过程开始
    for i = 1:n
        % 决策半径内找更优点
        Nit = [];                 % 存放萤火虫序号
        for j = 1:n
            if norm(gaddress(j, :)-gaddress(i, :)) < rdi(i) && ioti(i, 1) < ioti(j, 1)
                Nit(numel(Nit)+1) = j;
            end
        end
        % 找下一步移动的点开始
        if ~isempty(Nit)           
            Nitioti = ioti(Nit, 1);              % 选出Nit荧光素
            SumNitioti = sum(Nitioti);           % Nit荧光素和
            Molecular = Nitioti-ioti(i, 1);      % 分子
            Denominator = SumNitioti-ioti(i, 1); % 分母
            Pij = Molecular./Denominator;   % 计算Nit各元素被选择概率
            Pij = cumsum(Pij);        % 累计
            Pij = Pij./Pij(end);      % 归一化
            Pos = find(rand < Pij);   % 确定位置
            j = Nit(Pos(1));          % 确定j的位置
            % 萤火虫i向j移动一小步
            gaddress(i, :) = gaddress(i, :)+s*(gaddress(j, :)-gaddress(i, :))/norm(gaddress(j, :)-gaddress(i, :));
            % 边界处理(限制范围)
            gaddress(i, :) = min(gaddress(i, :), domx(1, 2));        
            gaddress(i, :) = max(gaddress(i, :), domx(1, 1)); 
            % 更新决策半径
            rdi(i) = rdi(i)+beta*(nt-length(Nit));
            if rdi(i, 1) < 0
                rdi(i, 1) = 0;
            end
            if rdi(i, 1) > rs
                rdi(i, 1) = rs;
            end
        end
    end
    % 每代最优解存入yy数组内
    yy(t) = max(fun(gaddress));
    % 迭代次数+1
    t = t+1;
end

%% 结果显示
gvalue = fun(gaddress);               % 求各个萤火虫的值
disp('最大值为：')
num = find(gvalue == max(gvalue));    % 最大值序号
MaxValue = max(gvalue)
disp('最优解为：')
BestAddress = gaddress(num, :)
figure;
plot(yy, 'r', 'linewidth', 2)
xlabel '迭代次数'; ylabel '函数值';
title 'GSO算法各代最优解变化';
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.
  • 7.
  • 8.
  • 9.
  • 10.
  • 11.
  • 12.
  • 13.
  • 14.
  • 15.
  • 16.
  • 17.
  • 18.
  • 19.
  • 20.
  • 21.
  • 22.
  • 23.
  • 24.
  • 25.
  • 26.
  • 27.
  • 28.
  • 29.
  • 30.
  • 31.
  • 32.
  • 33.
  • 34.
  • 35.
  • 36.
  • 37.
  • 38.
  • 39.
  • 40.
  • 41.
  • 42.
  • 43.
  • 44.
  • 45.
  • 46.
  • 47.
  • 48.
  • 49.
  • 50.
  • 51.
  • 52.
  • 53.
  • 54.
  • 55.
  • 56.
  • 57.
  • 58.
  • 59.
  • 60.
  • 61.
  • 62.
  • 63.
  • 64.
  • 65.
  • 66.
  • 67.
  • 68.
  • 69.
  • 70.
  • 71.
  • 72.
  • 73.
  • 74.
  • 75.
  • 76.
  • 77.
  • 78.
  • 79.
  • 80.
  • 81.
  • 82.
  • 83.
  • 84.
  • 85.
  • 86.
  • 87.
  • 88.
  • 89.
  • 90.

3、结果显示

Command Window中显示的结果为：

最大值为：

MaxValue =

       14416

最优解为：

BestAddress =

    -3    -3
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.
  • 7.
  • 8.
  • 9.
  • 10.
  • 11.

4、绘图

GSO算法各代最优解变化如图2所示。人工萤火虫优化(GSO)算法_优化算法_07

图2 GSO算法各代最优解变化

四、参考文献

代码下载或者仿真咨询添加QQ1575304183

[1] Krishnanand, K. N, Ghose, D. Detection of multiple source locations using a glowworm metaphor with applications to collective robotics[M]. 2005.
[2] 张谦. 基于群智能算法的无线传感器网络覆盖优化研究[D]. 2015.