%程序初始化

gen=50;     %設置進化代數
popsize=50;     %設置種群規模大小
best_in_history(gen)=inf;   %初始化全局歷史最優解
best_in_history(:)=inf;   %初始化全局歷史最優解
max_velocity=0.5;
%popnum=1;       %設置種群數量

pop(popsize,8)=0;   %初始化種群,創建popsize行5列的0矩陣
%種群數組第1列為x軸坐標，第2列為y軸坐標，第3列為x軸速度分量，第4列為y軸速度分量
%第5列為個體最優位置的x軸坐標，第6列為個體最優位置的y軸坐標
%第7列為個體最優適值，第8列為當前個體適應值

for i=1:popsize
    pop(i,1)=4*rand()-2;     %初始化種群中的粒子位置，值為-2—2，步長為其速度
    pop(i,2)=4*rand()-2;     %初始化種群中的粒子位置，值為-2—2，步長為其速度
    pop(i,5:6)=pop(:,1:2);  %初始狀態下個體最優值等于初始位置
    pop(:,3:4)=rand(popsize,2)*0.02-0.01;    %初始化種群微粒速度，值為-0.01—0.01，間隔為0.0001
    pop(:,7:8)=inf;

c1=2;
c2=2;
x_min=-2;
y_min=-2;
x_max=2;
y_max=2;

gbest_x=pop(1,1);   %全局最優初始值為種群第一個粒子的位置
gbest_y=pop(1,2);