PlatEMO：基于MATLAB的多目标优化工具箱完全指南

初做进化计算或者多目标优化研究时大部分人第一周都会做同一件事把论文里的算法伪代码翻译成 MATLAB 代码。这个过程听起来不难真做起来却全是坑测试函数写得不一致、指标计算结果对不上、对比算法复现不出论文效果、画出来的 Pareto 前沿还歪歪扭扭。等你把这些基础设施全部搞定两周时间已经过去了。PlatEMO 解决的就是这个问题。它是一个基于 MATLAB 平台的进化多目标优化工具箱内置了百余种优化算法、数十个标准测试问题和多种性能评价指标。它不只是一个算法集合更像是一套标准化实验框架算法、测试问题、评价指标、可视化都被统一封装你既可以直接调用现成算法做实验也可以把自己的研究成果以固定格式“插”进去和内置算法公平对比。这篇文章会从实际使用角度讲清楚 PlatEMO 是什么、怎么安装、GUI 和命令行两种实验方式怎么用、如何自定义算法和测试问题最后给出常见问题和工程建议。如果你正在做多目标优化相关的研究或开发建议先收藏再阅读。1. 为什么说 PlatEMO 是“一个搞定所有优化问题的工具”很多人第一次听到 PlatEMO 时会以为它只是“又一个大而全的算法库”。这种理解不能说错但会低估它的价值。它的意义不在于算法数量多而在于把进化多目标优化实验的整个流程标准化了。在没有 PlatEMO 之前做一组对比实验大致要经历以下步骤从论文中复现算法 A、B、C 的源代码每个作者的代码风格、变量命名、输入输出格式完全不同。在相同的测试函数上运行这些算法但测试函数的主体实现可能完全不同比如 ZDT1 在某些代码里是 30 维在另一些代码里是 10 维。计算 IGD、HV 等评价指标每个算法代码里可能都自带一套指标实现结果互相之间没有可比性。画图展示 Pareto 前沿有些算法输出的是目标值有些输出的是决策变量值需要手动转换。PlatEMO 把上述所有环节统一成了一个框架。它做的事情可以概括为三点算法部分内置了大量已经实现的进化算法从经典的 NSGA-II、MOEA/D到近年的各种改进算法都可以直接调用。问题部分内置了 ZDT、DTLZ、WFG、LZ 等多个系列的测试问题并统一了问题定义和维度设置。指标部分内置了 IGD、HV、GD、SP 等常用性能指标算法跑完直接计算不需要自己重新实现。最关键的在于这三个部分都遵循同一套接口规范。当你想对比自己的新算法和 NSGA-II 时不需要去网上找 NSGA-II 的代码也不需要担心对方代码里测试函数写得对不对。你只需要把自己的算法按照 PlatEMO 的模板写好然后由 PlatEMO 统一调度所有算法跑在同一个测试问题上指标计算方式完全一致。这才是它被称为“一个搞定所有优化问题”的真正原因。2. PlatEMO 的核心概念与设计逻辑在进入安装和使用之前有必要先理解 PlatEMO 的几个核心概念。因为它的所有扩展都围绕这几个概念展开。2.1 Algorithm算法在 PlatEMO 中每个算法都是一个继承自ALGORITHM基类的类。算法类需要实现run方法或者更准确地说实现进化迭代的主循环。PlatEMO 已经提供了大量的算法实现包括经典算法NSGA-II、NSGA-III、MOEA/D、SPEA2、PESA-II 等。基于分解的算法MOEA/D 及其多种变体。基于指标的算法IBEA、SMS-EMOA、HypE 等。基于种群的算法RVEA、LSMOF、KNEA 等。一些新兴的基于深度学习的算法也在陆续加入。对于使用者来说如果只是想跑对比实验内置算法基本够用。2.2 Problem测试问题测试问题在 PlatEMO 中也以类的形式存在继承自PROBLEM基类。每个问题类需要定义M目标个数。D决策变量个数。lower和upper决策变量的上下界。CalObj方法根据决策变量计算目标函数值。内置的测试问题覆盖了连续优化和离散优化中的很多经典测试集包括 ZDT、DTLZ、WFG、LZ、UF、MOP 等。在设置实验时可以直接通过ZDT1、DTLZ2这种函数句柄的方式引用。2.3 Metric评价指标PlatEMO 内置了多种性能指标最常用的是 IGD 和 HV。IGDInverted Generational Distance反转世代距离衡量算法求得的解集与真实 Pareto 前沿之间的逼近程度值越小越好。HVHypervolume超体积衡量解集在目标空间中覆盖的体积值越大越好。GD、SP 等也都可以通过配置文件或命令行参数指定。指标的一致性对于实验公平性非常重要这也是 PlatEMO 的重要价值之一。2.4 GUI 与命令行两种模式PlatEMO 提供两种操作方式GUI 模式通过图形界面选择算法、测试问题、种群大小和终止条件点击运行即可看到实时更新。适合演示和快速验证。命令行模式通过main(-algorithm, ...)函数批量运行实验适合需要跑多组对比实验的研究场景。这两者共用同一套底层框架实验结果格式也完全一致。下面会分别演示。3. 环境准备与安装PlatEMO 是一个 MATLAB 工具箱安装前需要确认本机已经安装 MATLAB。关于 MATLAB 版本的要求不同版本的 PlatEMO 略有差异建议参考官方 README本文重点演示通用使用思路。3.1 下载 PlatEMO从 GitHub 上搜索 BIMK/PlatEMO 即可找到官方仓库也可以从 MATLAB File Exchange 或相关发布页面下载压缩包。下载完成后解压到一个路径中建议不要放在 MATLAB 的安装目录下而是放在自己的工作目录中比如D:\PlatEMO或~/PlatEMO。3.2 添加路径到 MATLAB打开 MATLAB在主页点击“设置路径”将解压后的 PlatEMO 根目录以及所有子目录添加到 MATLAB 搜索路径中。也可以用命令完成% 将路径替换为你自己的解压路径 addpath(genpath(D:\PlatEMO)); savepath;这里真正容易踩坑的地方是如果使用addpath(genpath(...))但没有savepath重启 MATLAB 后路径会丢失。而使用了savepath之后后续每次打开 MATLAB 都能直接使用 PlatEMO。如果不想永久修改路径也可以在使用前手动执行一次addpath(genpath(...))。3.3 验证安装在 MATLAB 命令行输入PlatEMO如果安装成功会弹出 PlatEMO 的图形界面窗口。如果没有弹出说明路径没有添加成功或者下载的压缩包不完整。4. GUI 模式快速上手对于第一次接触 PlatEMO 的用户建议先通过 GUI 模式跑通一个最小实验。这样既能验证安装是否正确也能直观感受算法的收敛过程。4.1 启动 GUI 并选择算法在 MATLAB 命令行输入PlatEMO后主界面会分为几个区域算法选择区域从下拉列表中选择要运行的算法比如 NSGA-II。问题选择区域从下拉列表中选择测试问题比如 ZDT1。参数设置区域设置种群大小、目标函数个数、决策变量个数、最大评价次数等。运行区域点击运行按钮开始实验。通常第一次实验可以设置为算法 NSGA-II问题 ZDT1种群大小 100目标个数 2决策变量个数 30最大评价次数 10000。4.2 观察实验过程点击运行后GUI 会实时显示当前种群的 Pareto 前沿图像随着迭代次数增加前沿点会逐渐向真实 Pareto 前沿靠近。这个实时可视化过程对于理解算法的收敛行为和参数影响非常有帮助。4.3 批量比较多个算法在 GUI 中也可以选择多个算法进行对比PlatEMO 会以不同颜色展示各算法获得的 Pareto 前沿。这种直观对比在论文的初步实验阶段非常实用。GUI 模式适合快速验证但如果需要跑大量实验比如 10 个算法在 20 个测试问题上各跑 30 次再用 GUI 一个个点就非常低效这时必须使用命令行模式。5. 命令行模式与批量实验命令行模式是 PlatEMO 真正发挥效率的地方。通过一个main函数可以在一条命令中指定算法、问题、参数完成一次完整的实验。5.1 基本命令示例main(-algorithm, NSGAII, -problem, ZDT1, -N, 100, -M, 2, -D, 30, -maxFE, 10000);这条命令执行的功能是-algorithm指定使用 NSGA-II 算法。-problem指定测试问题为 ZDT1。-N种群大小为 100。-M目标函数个数为 2。-D决策变量个数为 30。-maxFE最大函数评价次数为 10000。运行结束后PlatEMO 会保存实验结果到 Data 文件夹包括算法获得的最终种群、运行过程数据等。5.2 批量实验脚本实际研究中通常需要在多个算法和多个问题上做对比实验。这时可以写一个 MATLAB 脚本用循环调用main函数% 文件路径run_comparison.m algorithms {NSGAII, MOEAD, SPEA2}; problems {ZDT1, ZDT2, DTLZ2}; for i 1 : length(algorithms) for j 1 : length(problems) main(-algorithm, algorithms{i}, ... -problem, problems{j}, ... -N, 100, ... -M, 2, ... -D, 30, ... -maxFE, 10000, ... -save, 1); end end一次性跑完所有组合每个组合的结果都保存在独立的文件夹中。命令行模式支持一个非常有用的参数-save设置为 1 时保存结果设置为 0 时只运行不保存。做大规模参数扫描时建议先设置为 0 做小规模试运行确认没有报错再正式运行。5.3 理解输出结构PlatEMO 运行完成后会在Data目录下生成以运行时间戳和实验参数命名的文件夹。其中包含Population最终获得的种群文件。Indicator计算出的性能指标结果。一些运行日志文件。这种统一输出格式的最大优点是无论你用了什么算法、什么问题最终拿到的数据接口是完全一致的。后续做统计分析和绘图就变得非常简单。6. 自定义算法把论文里的算法写成 PlatEMO 格式如果只是调用现成算法PlatEMO 和普通算法库没有本质区别。真正让 PlatEMO 有价值的是它的可扩展性你可以把自己的算法按照标准模板写成 PlatEMO 的算法类然后直接和内置算法对比。6.1 创建算法类在 PlatEMO 的根目录下的Algorithms文件夹中新建一个文件比如MyAlgorithm.m内容模板如下% 文件路径Algorithms/MyAlgorithm.m classdef MyAlgorithm ALGORITHM methods function run(Algorithm, Problem) % 初始化种群 Population Problem.Initialization(); % 进化主循环 while Algorithm.NotTerminated(Population) % 选择父代 MatingPool TournamentSelection(2, Problem.N, Population.objs); % 交叉变异产生子代 Offspring OperatorGA(Population(MatingPool)); % 环境选择选择下一代种群 Population EnvironmentalSelection([Population, Offspring], Problem.N); end end end end这个模板展示了 PlatEMO 算法类的基本结构算法类继承自ALGORITHM。run方法接收两个参数Algorithm和Problem。通过Problem.Initialization()初始化种群。Algorithm.NotTerminated(Population)判断是否达到终止条件。EnvironmentalSelection是 PlatEMO 提供的通用环境选择函数根据支配关系保留优秀个体。需要注意的是不同算法的主要区别在于“如何选择父代”和“如何更新种群”PlatEMO 并不限制你必须用什么算子。你可以完全自定义选择和更新过程例如引入自适应算子、局部搜索、分解策略等。6.2 调用自己的算法保存算法类之后在命令行中就可以和调用内置算法一样使用main(-algorithm, MyAlgorithm, -problem, ZDT1, -N, 100, -M, 2, -D, 30, -maxFE, 10000);这样就完成了从自定义算法到对比实验的整个流程。相比从零写实验代码这个模板明显缩短了开发时间。6.3 新手容易犯的三个错误自定义算法时有三个典型的错误经常遇到第一run方法中忘记调用NotTerminated导致死循环。PlatEMO 的终止控制是通过NotTerminated实现的它会在达到最大评价次数时返回 false。第二初始化种群时直接使用随机数生成而不是使用Problem.Initialization()。这会导致决策变量的上下界没有被正确读取算法求解的问题维度可能和设定不符。第三环境选择后没有更新种群就进入下一轮迭代导致种群数量不断增加或减少最终报错。7. 自定义测试问题除了自定义算法PlatEMO 还支持自定义测试问题。这在验证算法在特定工程问题上的表现时尤其有用。7.1 创建问题类与算法类类似测试问题也需要继承自PROBLEM基类。在 PlatEMO 根目录的Problems文件夹下新建文件MyProblem.m% 文件路径Problems/MyProblem.m classdef MyProblem PROBLEM methods function Setting(Problem) % 设置目标个数、决策变量个数和边界 Problem.M 2; Problem.D 30; Problem.lower zeros(1, Problem.D); Problem.upper ones(1, Problem.D); end function f CalObj(Problem, X) % 自定义目标函数这里以 ZDT1 为例 n size(X, 2); g 1 9 * mean(X(:, 2:end), 2); f(:, 1) X(:, 1); f(:, 2) g .* (1 - sqrt(X(:, 1) ./ g)); end end end关键方法说明Setting(Problem)问题设置阶段必须设置目标个数M、决策变量个数D、以及决策变量边界lower和upper。CalObj(Problem, X)目标函数计算函数输入变量矩阵X输出目标函数值矩阵f。7.2 调用自定义问题调用时只需替换-problem参数即可main(-algorithm, NSGAII, -problem, MyProblem, -N, 100, -M, 2, -D, 30, -maxFE, 10000);将实际问题建模成 PlatemO 的PROBLEM类后所有内置算法都可以直接用于求解该问题。这是 PlatEMO 从学术工具走向实际工程应用的关键环节。8. 运行结果与性能指标验证实验跑完不等于工作结束还需要量化评估算法性能。PlatEMO 内置的指标计算模块让这一过程变得标准化。8.1 查看指标结果在命令行模式下PlatEMO 运行完成后会计算默认指标并保存。如果想直接查看某个算法的 IGD 和 HV可以手动调用指标计算函数。PlatEMO 提供了一些工具函数可以通过结果文件夹中的数据计算% 示例加载种群结果并计算 IGD % 假设 Data 目录中存在运行结果路径以实际为准 resultPath .\Data\...; load(fullfile(resultPath, Population.mat)); % 如果测试问题已知真实 Pareto 前沿可以调用 % IGD(Population.objs, PF)PlatEMO 对不同测试问题内置了真实 Pareto 前沿数据运行结果中会保存这些信息计算指标时可以直接使用。8.2 IGD 与 HV 怎么看IGD 衡量的是算法获得的解集与真实 Pareto 前沿的平均距离IGD 越小说明解集越逼近真实前沿。HV 衡量的是解集在目标空间中的覆盖体积HV 越大说明解集的分布性和收敛性越好。在很多论文中IGD 和 HV 是通用的性能报告指标。如果你在自己的算法中观察到 IGD 下降但 HV 没有明显提升说明算法虽然更逼近前沿但在覆盖度上还有提升空间。两个指标结合才能对算法性能有一个更全面的判断。8.3 如何判断实验是否成功判断实验成功有几条基本原则算法正常结束没有报错。获得的 Pareto 前沿和目标空间中的真实前沿形状相近。多次独立运行的结果稳定指标波动在可接受范围内。同组对比实验中不同算法在相同条件下运行评价次数一致。如果运行失败先看 MATLAB 命令行窗口中的红色报错信息这是最直接的排查入口。9. 常见问题与排查方法PlatEMO 在安装和使用过程中有一些高频问题整理成下表问题现象可能原因排查方式解决方案打开 MATLAB 后输入 PlatEMO 报错“未定义函数或变量”PlatEMO 路径未添加或没有执行 savepath检查路径设置执行which PlatEMO重新addpath(genpath(...))并savepathGUI 能打开但运行实验时进度条不动种群数量和评价次数设置过大或算法本身效率较低查看命令行是否输出日志检查内存占用减小种群规模或最大评价次数先跑小规模测试报错“类 ALGORITHM 未找到”工具箱根目录没有被完整添加到路径中检查addpath(genpath(...))是否覆盖了Algorithms和Problems等目录确保使用genpath添加所有子目录自定义算法运行时报错“未定义函数 OperatorGA”调用了 PlatEMO 中不存在的算子名称或版本中函数名不同查看当前版本支持的算子函数列表替换为当前版本中实际存在的算子检查算子参数格式对比实验中各算法评价次数不一致不同算法使用了不同的最大评价次数或不同种群大小检查每个运行的参数配置统一-N和-maxFE参数结果指标出现 NaN 或 Inf测试问题定义有误或目标值出现除零等情况检查 CalObj 中是否有除零或 log 负数等操作在 CalObj 中对边界情况增加保护逻辑这里特别提醒一下PlatEMO 不同版本的函数名和目录结构可能会有调整代码示例中的类名和函数名以你实际下载版本中的官方模板为准。如果发现某个函数不存在优先查看该版本的官方文档或模板文件。10. 最佳实践与工程建议经过一段时间的使用我总结了一些 PlatEMO 的实践建议对准备把它用于正式研究或项目开发的同学会比较有用。10.1 先跑通最小示例再做大实验不要一上来就配置 10 个算法、20 个问题、每个跑 30 次。建议先用一个算法、一个问题、最大评价次数 1000跑一个只需要几秒钟的实验确认整套流程没有问题。再逐步扩大实验规模。10.2 自定义算法从复制现有算法开始如果你想实现一个新算法最快的路径不是从空模板开始写而是先复制一个与你的算法思路最接近的内置算法文件在它的基础上修改。比如你想做一个基于 MOEA/D 的改进算法就把 MOEAD.m 复制一份保持它的框架不变只修改核心部分。这样既不会遗漏 PlatEMO 接口也更容易定位问题。10.3 实验参数需要统一记录多目标优化实验结果受种群大小、评价次数、交叉变异概率等参数影响很大。建议在实验前把这些参数写入一个配置文件每次运行都从配置中读取。PlatEMO 的命令行参数虽然直观但如果你手动修改了几次参数后忘记记录实验结果的可信度就会打折扣。10.4 合理设置最大评价次数最大评价次数maxFE是影响实验结果稳定性的关键参数。设置过小算法没有收敛IGD 指标虚高设置过大实验时间过长。一般建议根据测试问题的复杂度确定比如在 ZDT 系列问题上 10000 到 50000 次评价可以观察到明显收敛更复杂的 DTLZ 问题则需要更多。10.5 运行结果尽量独立保存PlatEMO 默认会把每次运行的结果保存到 Data 文件夹。如果做大规模对比实验建议在每次运行后把重要的指标结果同步到统一的表格或文件中避免实验结束后要从几十个子文件夹里手工整理数据。10.6 注意 MATLAB 版本兼容性PlatEMO 使用了较新的 MATLAB 面向对象语法老版本 MATLAB 可能出现语法不支持的问题。如果在老旧环境中使用建议先检查官方说明中标注的最低 MATLAB 版本。此外在虚拟机中运行 MATLAB 本身性能会下降如果实验规模较大优先在物理机或计算集群上运行。10.7 安全与合规提示如果你要在生产环境或团队共享环境中使用 PlatEMO注意确认工具箱的许可证类型和引用要求。学术使用场景下通常在论文中引用 PlatEMO 的论文即可。如果涉及商业项目建议先与原作者确认授权范围。11. 关于优化算法学习的延伸思考使用 PlatEMO 的价值不只是“跑实验方便”。它提供了一个非常好的学习路径你可以在同一个框架下比较不同算法的行为差异通过可视化观察算法是如何一步步逼近 Pareto 前沿的。这种直观理解比单纯读论文来得深刻得多。如果你正在学习进化多目标优化建议按照以下路径实践用 PlatEMO 跑通 NSGA-II 和 MOEA/D对比两者的 Pareto 前沿差异。修改种群大小、评价次数、交叉变异参数观察对结果的影响。在 DTLZ 系列问题上观察不同算法的表现差异理解算法在多目标场景下的瓶颈。用 PlatEMO 自定义一个简单的改进算法比如在 NSGA-II 中加入新的环境选择策略跑实验验证指标变化。逐步扩展到真实工程问题把实际问题建模成测试问题类用 PlatEMO 内置算法求解。只要把前两步做完你对 PlatEMO 的理解就已经超过大多数只跑过 demo 的初学者。而当你开始尝试自定义算法时才真正进入了进化计算研究的大门。