粒子群算法源代码

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
%------初始化种群的个体------------

for i=1:N

    for j=1:D

        x（ij）=randn;  %随机初始化位置

        v（ij）=randn;  %随机初始化速度

    end

end

%------先计算各个粒子的适应度，并初始化Pi和Pg----------------------

for i=1:N

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

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

end

pg = x（N:）;             %Pg为全局最优

for i=1:（N-1）

    if fitness（x（i