pso.m

pso程序

%标准粒群优化算法程序
%测试函数：f(x,y)=100(x^2-y)^2+(1-x)^2, -2.048<x,y<2.048
%求解函数最小值

clc;
clf;
clear;

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坐标数组

initial;                % 初始化
date;                   % 加载原始点

for exetime=1:gen
    outputdata;         % 实时输出结果
    update_position;
    update_velocity;
    pause(0.01);
end
swarm(index_gbest,1,:)
swarm(index_gbest,4,1)