布局优化基于萤火虫算法的无线传感器网(WSN)覆盖优化Matlab源码

Posted 2021-08-02 MatlabQQ1575304183

 一、WSN模型

二、萤火虫GSO算法

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

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

2、荧光素更新阶段

3、移动概率计算阶段

4、位置更新阶段

5、邻域范围更新阶段

三、代码

%% 初始化参数
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算法各代最优解变化');

 

四、参考文献