%標準粒群優(yōu)化算法程序
% 2007.1.9 By 郟宣耀
%測試函數(shù)：f(x,y)=100(x^2-y)^2+(1-x)^2, -2.048<x,y<2.048 f=20+x.^2-10*cos(2*pi*x)+y.^2-10*cos(2*pi*y), -5.12<=x,y<=5.12
%求解函數(shù)最小值
%f(x)=x^2-2*x);
global d=
L=0.5ones(N,N)-diag(diag(d))
global N=100
global popsize;     %種群規(guī)模
global popnum;      %種群數(shù)量
global pop;         %種群
%global c0;          %速度慣性系數(shù),為0—1的隨機數(shù)
%global c1;          %個體最優(yōu)導(dǎo)向系數(shù)
global c;          %全局最優(yōu)導(dǎo)向系數(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;       %計算適應(yīng)值
 % errorcompute; %計算當前種群適值標準差
    updatepop;      %更新粒子位置
    pause(0.01);
end


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