粒子群优化算法matlab代码请写出来

以下是粒子群优化算法的 MATLAB 代码示例：

function [best_pos, best_val] = PSO(fitness_func, dim, lb, ub, max_iter, swarm_size, w, c1, c2)

% 初始化粒子群
pos = rand(swarm_size, dim) .* (ub - lb) + lb; % 粒子位置
vel = zeros(swarm_size, dim); % 粒子速度
pbest_pos = pos; % 个体最优位置
pbest_val = inf(swarm_size, 1); % 个体最优值
gbest_pos = zeros(1, dim); % 全局最优位置
gbest_val = inf; % 全局最优值

% 迭代优化
for iter = 1:max_iter
% 更新粒子速度和位置
vel = w * vel + c1 * rand(swarm_size, dim) .* (pbest_pos - pos) + c2 * rand(swarm_size, dim) .* (repmat(gbest_pos, swarm_size, 1) - pos);
pos = pos + vel;
% 边界处理
pos(pos < lb) = lb(pos < lb);
pos(pos > ub) = ub(pos > ub);
% 计算适应度值
val = fitness_func(pos);
% 更新个体最优位置和全局最优位置
update = val < pbest_val;
pbest_val(update) = val(update);
pbest_pos(update, :) = pos(update, :);
[min_val, min_idx] = min(pbest_val);
if min_val < gbest_val
gbest_val = min_val;
gbest_pos = pbest_pos(min_idx, :);
end
end

% 返回最优解
best_pos = gbest_pos;
best_val = gbest_val;

end

注意：这只是一个示例代码，实际使用时需要根据具体问题进行修改和调整。