update_velocity.m

pso程序

function update_velocity

global swarm_size;      % 种群规模
global swarm;           % 种群,三维数组
global w;               % 速度惯性系数,为0—1的随机数
global c1;              % 个体最优导向系数
global c2;              % 全局最优导向系数
global x_min;           % x的下限
global x_max;           % x的上限
global y_min;           % y的下限
global y_max;           % y的上限
global gen;             % 迭代次数
global exetime;         % 当前迭代次数
global max_velocity;    % 最大速度
global best_in_history; % 最优解变化轨迹
global index_gbest;     % 全局最优解的粒子的序号
global dimension;       % 空间的维数
global points_size;     % points_size--行数，即点数
global x;               % 初始点的x坐标数组
global y;               % 初始点的y坐标数组

% 更新粒子的速度
for i = 1 : swarm_size
    for j = 1:dimension
        swarm(i, 2, j) = w*rand*swarm(i, 2, j) + c1*rand*(swarm(i, 3, j) - swarm(i, 1, j)) + c2*rand*(swarm(index_gbest, 3, j) - swarm(i, 1, j));   %x velocity component
      
%         % 限制粒子的运动速度
%         if abs(swarm(i, 2, j)) > max_velocity
%             if swarm(i, 2, j) > 0
%                 swarm(i, 2, j) = max_velocity;
%             else
%                 swarm(i, 2, j) = -max_velocity;
%             end
%         end
        
    end
end