%標準粒群優化算法程序
% 2007.1.9 By 郟宣耀
%測試函數：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
%求解函數最小值
%f(x)=x^2-2*x);
global N;
N=10;
global d;
d=0.2*ones(N,N);
L=0.5*ones(N,N)-diag(diag(d))

global popsize;     %種群規模
global popnum;      %種群數量
global pop;         %種群
%global c0;          %速度慣性系數,為0—1的隨機數
%global c1;          %個體最優導向系數
global c;          %全局最優導向系數
c=0.2*ones(N,N);
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.01);
end


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