%標準粒群優(yōu)化算法程序
% 2007.1.9 By 郟宣耀
%測試函數(shù)：f(x,y)=100(x^2-y)^2+(1-x)^2, -2.048<x,y<2.048
%求解函數(shù)最小值

global popsize;     %種群規(guī)模
%global popnum;      %種群數(shù)量
global pop;         %種群
%global c0;          %速度慣性系數(shù),為0—1的隨機數(shù)
global c1;          %個體最優(yōu)導向系數(shù)
global c2;          %全局最優(yōu)導向系數(shù)
global gbest_x;       %全局最優(yōu)解x軸坐標
global gbest_y;       %全局最優(yōu)解y軸坐標
global best_fitness;    %最優(yōu)解
global best_in_history; %最優(yōu)解變化軌跡
global x_min;           %x的下限
global x_max;           %x的上限
global y_min;           %y的下限
global y_max;           %y的上限
global gen;             %迭代次數(shù)
global exetime;         %當前迭代次數(shù)
global max_velocity;    %最大速度

initial;        %初始化

for exetime=1:gen
    outputdata;     %實時輸出結(jié)果
    adapting;       %計算適應值
    errorcompute(); %計算當前種群適值標準差
    updatepop;      %更新粒子位置
    pause(0.01);
end


clear i;
clear exetime;
clear x_max;
clear x_min;
clear y_min;
clear y_max;