function?[xmfv]?=?PSO（fitnessNc1c2wMD）
function?f=fitness（x）
f=x（1）.^2+x（2）.^2+x（3）.^2
end
format?long;
%100：N
%2:c1
%2：c2
%0.6:w
%1000:M
%3:D
%------初始化種群的個(gè)體------------

for?i=1:N

????for?j=1:D

????????x（ij）=randn;??%隨機(jī)初始化位置

????????v（ij）=randn;??%隨機(jī)初始化速度

????end

end

%------先計(jì)算各個(gè)粒子的適應(yīng)度，并初始化Pi和Pg----------------------

for?i=1:N

????p（i）=fitness（x（i:））;

????y（i:）=x（i:）;

end

pg?=?x（N:）;?????????????%Pg為全局最優(yōu)

for?i=1:（N-1）

????if?fitness（x（i