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

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

initial;        %初始化

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


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