粒子群算法pso优化源代码

粒子群算法（PSO）是一种基于群体协同的优化算法。该算法最初是由Kennedy和Eberhart在1995年提出的，它模拟鸟群或鱼群等生物群体的行为，通过在搜索空间中寻找最优解来进行优化。

PSO算法中的优化问题通常由目标函数确定。算法通过初始化粒子群中每个粒子的位置和速度，然后使用目标函数来评价每个粒子的适应度值。接着，每个粒子通过和其他粒子的交互来更新自己的位置和速度，并且更新最优解。最后，纪录最好的粒子状态，返回最优解。

以下是一个基于MATLAB编写的PSO优化源代码，代码实现了在目标函数中最小化的问题。假设目标函数为f(x)=-cos(x[1])*cos(x[2])*exp(-(x[1]-pi)^2-(x[2]-pi)^2)，其中x[1]和x[2]为粒子群中的每个粒子的向量：

function [Gbest,Gfit] = pso(n,dim,lb,ub,max_iter)

% n代表粒子群的数量，dim代表每个粒子向量的维度数，
% lb和ub分别代表向量的下限和上限，max_iter代表最大迭代次数。

% 初始化粒子和速度
V = zeros(n,dim);
X = repmat(lb,n,1) + repmat((ub-lb),n,1).*rand(n,dim);

% 计算适应度
fit = zeros(n,1);
for i=1:n
fit(i) = fitness(X(i,:)); % 调用fitness函数计算每个粒子的适应度
end

% 纪录最好的粒子和适应度
[Gfit,idx] = min(fit);
Gbest = X(idx,:);

% 开始迭代
for iter=1:max_iter
% 更新速度和位置
for i=1:n
r1 = rand(1,dim);
r2 = rand(1,dim);
V(i,:) = V(i,:) + r1.*(Gbest-X(i,:)) + r2.*(X(idx,:) - X(i,:));
X(i,:) = X(i,:) + V(i,:);
end

% 边界约束
X(X<lb) = lb(X<lb);
X(X>ub) = ub(X>ub);

% 更新适应度
for i=1:n
fit(i) = fitness(X(i,:)); % 调用fitness函数计算每个粒子的适应度
end

% 更新最好的粒子和适应度
[fit_i,idx_i] = min(fit);
if fit_i < Gfit
Gfit = fit_i;
Gbest = X(idx_i,:);
end
end

% fitness函数
function [f] = fitness(x)
f = -cos(x(1))*cos(x(2))*exp(-(x(1)-pi)^2-(x(2)-pi)^2);
end

end

需要注意的是，PSO算法的性能是与参数的选择密切相关的。在实际应用中，很多赋值都是基于经验值。需要根据不同的问题进行讨论和选择。