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

%計算適應值并賦值
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)    %若當前適應值優于個體最優值，則進行個體最優信息的更新
        pop(i,7)=pop(i,8);          %適值更新
        pop(i,5:6)=pop(i,1:2);      %位置坐標更新
    end
end

%計算完適應值后尋找當前全局最優位置并記錄其坐標
best_fitness=min(pop(:,7));     %全局最優值
gbest_x=pop(find(pop(:,7)==best_fitness),1);    %全局最優粒子的位置   
gbest_y=pop(find(pop(:,7)==best_fitness),2);

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