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