粒子群Matlab

function [par par_best]=update_par（parpar_best）
    
    %Px=Px+Pv*t这里t=1Px为当前粒子的位置，Pv为当前粒子的速度
    par.x=par.x+par.vx;   
    par.y=par.x+par.vy;   
    par.fit=compute_fit（par）;    %计算当前粒子适应度
    %Pv=Pv+（c1*rand*（Gx-Px））+（c2*rand*（PBx-Px））
    %这里c1c2为加速因子
    %Gx为具有最佳适应度粒子的位置
    %PBx为当前粒子的最佳位置
    c1=1;
    c2=1;
    par.vx=par.vx+c1*rand（）*（par_best.x-par.x）+c2*rand（）*（par.bestx-par.x）;   
    par.vy=par.vy+c1*rand（）*（pa