MATLAB PSO算法结果去重：生成唯一解矩阵

本代码使用MATLAB实现粒子群优化算法（PSO），并对算法运行结果进行去重，生成包含唯一解的矩阵。

目标:

循环执行1000次PSO算法，每次获取三个参数 [xopt, yopt, fopt]。
将这三个参数存储到一个矩阵中，但需要确保矩阵中每个元素都是唯一的。
设置一个阈值 eps，如果两个解之间的差异小于eps，则认为它们是相同的。

代码实现:

mat = []; % 初始化矩阵
eps = 1e-6; % 设置eps值
for i = 1:1000
    [xopt, yopt, fopt] = PSO(@Shubert, np, maxiter, bounds); % 调用PSO函数
    if isempty(mat) % 如果矩阵为空，直接将结果存入矩阵中
        mat = [xopt, yopt, fopt];
    else
        % 判断当前结果是否和矩阵中已有的结果相同
        flag = false;
        for k = 1:size(mat, 1)
            if abs(mat(k,1) - xopt) < eps && abs(mat(k,2) - yopt) < eps && abs(mat(k,3) - fopt) < eps
                flag = true;
                break;
            end
        end
        % 如果不相同，则将结果存入矩阵中
        if ~flag
            mat = [mat; xopt, yopt, fopt];
        end
    end
end
% 输出矩阵
disp(mat);

代码说明:

mat 为存储唯一解的矩阵，初始化为空矩阵。
eps 设置为1e-6，用于判断两个解是否相同。
循环1000次执行PSO算法，每次获取结果 [xopt, yopt, fopt]。
如果 mat 为空，则直接将结果存入矩阵中。
如果 mat 不为空，则遍历矩阵，判断当前结果是否与矩阵中已有结果相同。
如果结果与矩阵中已有结果相同，则标记 flag 为 true，并跳出循环。
如果结果与矩阵中已有结果不同，则将结果存入矩阵中。
最后输出 mat，即包含所有唯一解的矩阵。

注意事项:

@Shubert、np、maxiter、bounds 等参数需根据具体应用场景进行设置。
eps 的大小会影响解的判定，需要根据实际情况进行调整。

应用场景:

使用PSO算法进行优化问题求解时，可以利用本代码进行结果去重，得到唯一解集合。
可以通过对唯一解的分析，深入理解PSO算法的运行机制和优化效果。

其他:

本代码仅供参考，可以根据实际需求进行修改和调整。
为了更清晰地展示代码，省略了一些注释，实际应用中建议添加详细的注释。
欢迎讨论和交流！