小生境pso matlab

以下是一个简单的PSO算法的Matlab实现：

% PSO算法Matlab实现
% 问题：求f(x) = x^2的最小值

clear all;
clc;

% 参数设置
w = 0.6; % 惯性权重
c1 = 1.8; % 学习因子
c2 = 1.8; % 学习因子
maxgen = 100; % 最大迭代次数
sizepop = 20; % 群体大小
Vmax = 5; % 最大速度限制
Xmax = 5; % X的范围
Xmin = -5;

% 初始化粒子位置和速度
for i = 1:sizepop
    X(i) = Xmax - (Xmax - Xmin) * rand(1); % 粒子位置
    V(i) = Vmax - 2 * Vmax * rand(1); % 粒子速度
end

% 计算每个粒子的适应值
for i = 1:sizepop
    fitness(i) = X(i)^2;
end

% 记录全局最优解
[bestfitness, bestindex] = min(fitness);
bestX = X(bestindex);

% 开始迭代
for i = 1:maxgen
    % 更新每个粒子的速度和位置
    for j = 1:sizepop
        V(j) = w * V(j) + c1 * rand(1) * (bestX - X(j)) + c2 * rand(1) * (X(bestindex) - X(j));
        if V(j) > Vmax
            V(j) = Vmax;
        elseif V(j) < -Vmax
            V(j) = -Vmax;
        end
        X(j) = X(j) + V(j);
        if X(j) > Xmax
            X(j) = Xmax;
        elseif X(j) < Xmin
            X(j) = Xmin;
        end
    end
    
    % 计算每个粒子的适应值
    for j = 1:sizepop
        fitness(j) = X(j)^2;
    end
    
    % 更新全局最优解
    [tempfitness, tempindex] = min(fitness);
    if tempfitness < bestfitness
        bestfitness = tempfitness;
        bestindex = tempindex;
        bestX = X(bestindex);
    end
    
    % 显示当前迭代次数和最优解
    disp(['迭代次数：', num2str(i), '，最优解：', num2str(bestfitness)]);
end

% 输出最优解
disp(['PSO算法求解f(x) = x^2的最小值为：', num2str(bestfitness)]);

这里实现的是求解$f(x) = x^2$的最小值。你可以根据需要修改目标函数以及相应的参数设置。