%適值計(jì)算
% 測試函數(shù)為f(x,y)=100(x^2-y)^2+(1-x)^2, -2.048<x,y<2.048

%計(jì)算適應(yīng)值并賦值
for i=1:popsize
    pop(i,8)=100.*(pop(i,1).^2-pop(i,2)).^2+(1-pop(i,1)).^2;
    if pop(i,7)>pop(i,8)    %若當(dāng)前適應(yīng)值優(yōu)于個(gè)體最優(yōu)值，則進(jìn)行個(gè)體最優(yōu)信息的更新
        pop(i,7)=pop(i,8);          %適值更新
        pop(i,5:6)=pop(i,1:2);      %位置坐標(biāo)更新
    end
end

%計(jì)算完適應(yīng)值后尋找當(dāng)前全局最優(yōu)位置并記錄其坐標(biāo)
best_fitness=min(pop(:,7));     %全局最優(yōu)值
gbest_x=pop(find(pop(:,7)==best_fitness),1);    %全局最優(yōu)粒子的位置   
gbest_y=pop(find(pop(:,7)==best_fitness),2);

%更新粒子速度與位置
for i