拓冰网站定制 
简介这份资源面向具备一定Python基础、希望系统掌握多目标优化算法的学习者与工程研究人员围绕非支配排序遗传算法NSGA-II展开帮助解决多个相互制约目标下的帕累托最优求解问题。压缩包共10个文件约518KB包含4个ipynb交互式笔记本、1个py脚本、1个pdf文档、1个txt说明及若干zbak备份文件覆盖算法实现、Pareto前沿最优解选取与工具库对比等模块。其中笔记本完整呈现从初始种群生成、分层排序、密度评估到精英保留与遗传算子应用的六步流程并涉及DEAP、Platypus等框架的编程实践pdf则补充了前沿解筛选思路。已有94人学习下载适合作为课程设计、科研入门或工程优化的实操参考便于读者对照代码理解算法内在机理并快速迁移到自身问题建模中。1. 从一组 Pareto 前沿说起NSGA-II 到底解决了什么工程问题如果你做过参数调优大概率遇到过这种场景模型准确率上去了推理延迟也炸了成本压下来了精度又没法看。单目标优化只能给你一个“最优解”但工程里真正难的是在多个互相拉扯的目标之间找一组“谁也不比谁差”的候选方案。NSGA-II非支配排序遗传算法第二代就是干这个的——它不返回一个解而是返回一条 Pareto 前沿让你根据业务偏好去挑。这篇笔记围绕 Python 实现 NSGA-II 多目标优化算法把 Jupyter 代码详解与下载这件事拆开讲透算法流程怎么走、代码每一块在干什么、参数怎么调、Jupyter 里怎么跑通、哪些坑我踩过。适合已经会写 Python、想把这套算法落到自己优化问题上的工程师也适合刚接触多目标优化、想找一个能直接复现的 Jupyter 实现的人。2. NSGA-II 的算法骨架非支配排序、拥挤度与精英保留2.1 为什么是 NSGA-II而不是加权求和或 MOEA/D很多人第一反应是“多目标不就是把几个目标加权加起来当单目标做吗”。加权求和的问题在于权重极难定而且它只能找到 Pareto 前沿上的凸部分遇到凹前沿直接失效。NSGA-II 的核心优势有三个一是用非支配排序给种群分层保证收敛到前沿二是用拥挤度距离维持解的多样性避免全挤在一个点三是精英保留策略父代和子代合并后再选防止已经找到的好解丢失。这三点加起来让它在两到三个目标的优化问题上非常稳。常见做法是目标数不超过 3 时优先用 NSGA-II超过 3 个目标再考虑 NSGA-III 或 MOEA/D。2.2 非支配排序与拥挤度的计算逻辑非支配排序的本质是对每个个体统计有多少个体支配它 domination count 以及它支配了哪些个体。第一层是 domination count 为 0 的个体把它们移除后被它们支配的个体 count 减一再找出新的 0 层依次类推。拥挤度则是把同一层内的个体按每个目标排序取相邻个体在各自目标上的距离归一化后求和。边界个体的拥挤度设为无穷大保证它们一定被保留。def fast_non_dominated_sort(population): # population: list of individuals, each has .objectives (list of float) fronts [[]] for p in population: p.domination_count 0 p.dominated_solutions [] for q in population: if p.dominates(q): p.dominated_solutions.append(q) elif q.dominates(p): p.domination_count 1 if p.domination_count 0: p.rank 0 fronts[0].append(p) i 0 while fronts[i]: next_front [] for p in fronts[i]: for q in p.dominated_solutions: q.domination_count - 1 if q.domination_count 0: q.rank i 1 next_front.append(q) i 1 fronts.append(next_front) return fronts[:-1] # 最后一层为空这段代码里dominates方法需要你自己在个体类里实现判断条件是对所有目标p 不差于 q且至少有一个目标严格优于 q。domination_count和dominated_solutions是挂在个体对象上的临时属性每次排序前要重置。参数上唯一要注意的是 fronts 的初始化fronts [[]]后面 append 空列表会导致多一层所以返回时用fronts[:-1]去掉。2.3 选择、交叉、变异在 NSGA-II 里的具体接法NSGA-II 的选择用二元锦标赛但比较规则是先比 rankrank 小的赢rank 相同比拥挤度拥挤度大的赢。交叉和变异就是常规的模拟二进制交叉SBX和多项式变异。SBX 的分布指数eta_c一般取 15 到 20多项式变异的eta_m取 20 左右变异概率p_m通常取1/n_vars。这些参数不是玄学eta越大子代越靠近父代探索能力越弱太小又退化成随机搜索。def tournament_selection(pop, k2): candidates random.sample(pop, k) best candidates[0] for c in candidates[1:]: if c.rank best.rank: best c elif c.rank best.rank and c.crowding_distance best.crowding_distance: best c return best锦标赛大小k默认 2 就够k 越大选择压力越大收敛快但容易早熟。拥挤度在边界个体上是float(inf)比较时不会出错但如果你用 numpy 的 inf 要注意别参与后续的归一化计算。3. 在 Jupyter 里把 NSGA-II 跑起来从环境到第一个 Pareto 前沿3.1 Jupyter 环境准备与依赖安装Jupyter Notebook 的安装方式现在主流是走 Anaconda 或者 pip。如果你用 Anacondaconda install jupyter notebook就行如果走 pippip install notebook也可以。Python 版本建议 3.9 以上numpy 和 matplotlib 是必装的。我一般会在项目目录下建一个虚拟环境避免和系统 Python 打架。python -m venv nsga2_env source nsga2_env/bin/activate # Windows 用 nsga2_env\Scripts\activate pip install numpy matplotlib notebook jupyter notebook启动后浏览器会自动打开如果没打开终端里会有一行带 token 的地址复制到浏览器即可。Jupyter 默认存放地址是启动时所在的目录如果你想让它固定在某个文件夹可以在启动前cd过去或者改配置文件里的notebook_dir。常见坑是 Anaconda 的 Jupyter 突然打不开多半是端口被占或者配置文件被改坏了用jupyter notebook --port 8889换个端口试试。3.2 用 Jupyter 组织 NSGA-II 代码的推荐结构一个能复现的 NSGA-II 实现我建议在 Jupyter 里按 cell 分块第一个 cell 放导入和全局参数第二个 cell 定义个体类和支配判断第三个 cell 写非支配排序和拥挤度第四个 cell 写交叉变异第五个 cell 写主循环第六个 cell 做可视化。这样调试的时候可以单独跑某个 cell不用每次从头执行。Jupyter 一个 cell 只输出最后一个表达式的结果如果你要同时看多个变量用print或者display。import numpy as np import matplotlib.pyplot as plt import random # 全局参数 POP_SIZE 100 N_GEN 200 N_VARS 2 ETA_C 20 ETA_M 20 P_M 1.0 / N_VARS BOUNDS [(-5, 5), (-5, 5)]POP_SIZE太小前沿覆盖不全太大跑得慢100 到 200 是常见起点。N_GEN看问题复杂度简单问题 100 代就收敛复杂问题可能要 500 代以上。BOUNDS是决策变量的上下界SBX 和多项式变异都要用。3.3 主循环与结果可视化主循环的逻辑是初始化种群评估目标非支配排序加拥挤度然后循环选择、交叉、变异生成子代父子合并再排序选前 N 个。可视化用 matplotlib 画散点图横纵轴分别是两个目标值。def nsga2_main(): pop [Individual(random.uniform(*BOUNDS[i]) for i in range(N_VARS)) for _ in range(POP_SIZE)] for ind in pop: ind.evaluate() fronts fast_non_dominated_sort(pop) assign_crowding_distance(fronts) for gen in range(N_GEN): offspring [] while len(offspring) POP_SIZE: p1 tournament_selection(pop) p2 tournament_selection(pop) c1, c2 sbx_crossover(p1, p2, ETA_C, BOUNDS) mutate(c1, ETA_M, P_M, BOUNDS) mutate(c2, ETA_M, P_M, BOUNDS) c1.evaluate(); c2.evaluate() offspring.extend([c1, c2]) combined pop offspring[:POP_SIZE] fronts fast_non_dominated_sort(combined) assign_crowding_distance(fronts) pop [] for front in fronts: if len(pop) len(front) POP_SIZE: pop.extend(front) else: front.sort(keylambda x: x.crowding_distance, reverseTrue) pop.extend(front[:POP_SIZE - len(pop)]) break return popcombined的长度是 2 倍POP_SIZE选前POP_SIZE个的时候如果某一层放不下就按拥挤度从大到小截断。这个截断逻辑是 NSGA-II 精英保留的关键写错会导致种群数量不稳定。跑完后用plt.scatter把第一层前沿画出来横纵轴分别是f1和f2。4. 参数调优与常见翻车现场排查4.1 交叉变异参数怎么设才不跑偏SBX 的eta_c控制子代和父代的接近程度值越大子代越像父代。我一般从 15 开始试如果前沿收敛太慢就降到 10如果种群多样性不够就升到 25。多项式变异的eta_m类似20 是默认值。变异概率p_m用1/n_vars是经验公式变量多的时候自动降低变异强度。注意 SBX 和多项式变异都要处理边界超出BOUNDS的变量要截断或者反射回来否则目标函数可能报错。4.2 非支配排序结果不对的三种典型原因第一种是dominates方法写反了把“不差于”写成了“严格优于”导致支配关系全错。第二种是排序前没重置domination_count第二次调用时 count 还是上一轮的值。第三种是 fronts 初始化多了一层空列表导致 rank 偏移。排查方法很简单拿两个明显有支配关系的个体手动跑一遍dominates看返回值对不对然后在排序函数里打印每层的个体数量正常应该是递减的。4.3 Jupyter 里跑大规模种群的性能问题Python 原生循环做非支配排序是 O(MN²)种群 500 以上、代数 500 以上就会明显卡。Jupyter 里可以用%%time魔法命令测每个 cell 的耗时。优化方向有两个一是用 numpy 向量化支配判断把目标值存成矩阵用广播比较二是换用 pymoo 这类成熟库它底层用 numpy 优化过。如果只是学习算法流程原生实现够用如果要上生产建议直接看 pymoo 的 NSGA-II 实现。5. 避坑与常见问题那些让我重跑一晚上的细节5.1 现象前沿点全挤在一起多样性极差原因拥挤度计算时没有归一化或者边界个体的拥挤度没设成无穷大。解决在assign_crowding_distance里对每个目标先取 min 和 max归一化后再算相邻距离边界个体直接赋float(inf)。5.2 现象种群数量每代都在变最后只剩几个解原因精英保留的截断逻辑写错了比如pop.extend(front)之后没有 break或者截断时用了front[:POP_SIZE]而不是front[:POP_SIZE - len(pop)]。解决严格按“能放整层就放整层放不下就按拥挤度截断”的逻辑写截断后立即 break。5.3 现象目标函数报数组维度错误原因决策变量是 list目标函数里直接当 numpy 数组用或者evaluate返回的不是标量。解决在evaluate里把变量转成 numpy 数组再算返回值确保是 float。Jupyter 里可以用%debug进事后调试看具体是哪一维对不上。5.4 现象跑了几十代前沿不再改善原因变异概率太低或者eta_m太大种群失去探索能力。解决把p_m临时调大或者给eta_m加个退火策略前期小后期大。也可以每隔几代随机重置一小部分个体。5.5 现象Jupyter 内核频繁挂掉原因种群太大内存爆了或者某个 cell 里有死循环。解决先用小种群跑通流程确认逻辑无误再放大。Jupyter 里可以用top或者任务管理器看内存占用超过 80% 就要考虑减种群或减代数。6. 进阶技巧用 pymoo 交叉验证你的实现以及一个收敛判断的小习惯自己手写 NSGA-II 最大的价值是理解流程但验证结果对不对我一般会拿 pymoo 跑同一个问题做对比。pymoo 的 API 很简洁定义问题、算法、终止条件三行就能跑。如果你的前沿和 pymoo 的前沿基本重合说明实现没问题如果差很多优先检查支配判断和拥挤度。from pymoo.algorithms.moo.nsga2 import NSGA2 from pymoo.problems import get_problem from pymoo.optimize import minimize problem get_problem(zdt1) algorithm NSGA2(pop_size100) res minimize(problem, algorithm, (n_gen, 200), seed1, verboseFalse) print(res.F.shape) # Pareto 前沿的目标值zdt1是标准测试问题前沿是凸的适合入门验证。res.F是前沿上的目标值矩阵行数是解的数量。跑完后可以和你的实现画在同一张图上对比。另一个习惯是每跑完一次把第一层前沿的f1最小值、f2最小值记下来看随代数怎么变。如果连续 20 代这两个值都不再下降基本可以停。这个判断比固定代数更省时间也避免跑过头。Jupyter 里可以写个简单的回调函数每 10 代打印一次当前前沿的极值。最后说个血泪教训别在 Jupyter 里直接跑 1000 代不存中间结果内核一挂全没了。我现在的习惯是每 50 代把种群的目标值存成 npy 文件文件名带代数这样即使后面崩了也能从最近一次恢复。希望帮到你。本文还有配套的精品资源点击获取 Python实现NSGA-II多目标优化：Jupyter代码详解与下载
简介这份资源面向具备一定Python基础、希望系统掌握多目标优化算法的学习者与工程研究人员围绕非支配排序遗传算法NSGA-II展开帮助解决多个相互制约目标下的帕累托最优求解问题。压缩包共10个文件约518KB包含4个ipynb交互式笔记本、1个py脚本、1个pdf文档、1个txt说明及若干zbak备份文件覆盖算法实现、Pareto前沿最优解选取与工具库对比等模块。其中笔记本完整呈现从初始种群生成、分层排序、密度评估到精英保留与遗传算子应用的六步流程并涉及DEAP、Platypus等框架的编程实践pdf则补充了前沿解筛选思路。已有94人学习下载适合作为课程设计、科研入门或工程优化的实操参考便于读者对照代码理解算法内在机理并快速迁移到自身问题建模中。1. 从一组 Pareto 前沿说起NSGA-II 到底解决了什么工程问题如果你做过参数调优大概率遇到过这种场景模型准确率上去了推理延迟也炸了成本压下来了精度又没法看。单目标优化只能给你一个“最优解”但工程里真正难的是在多个互相拉扯的目标之间找一组“谁也不比谁差”的候选方案。NSGA-II非支配排序遗传算法第二代就是干这个的——它不返回一个解而是返回一条 Pareto 前沿让你根据业务偏好去挑。这篇笔记围绕 Python 实现 NSGA-II 多目标优化算法把 Jupyter 代码详解与下载这件事拆开讲透算法流程怎么走、代码每一块在干什么、参数怎么调、Jupyter 里怎么跑通、哪些坑我踩过。适合已经会写 Python、想把这套算法落到自己优化问题上的工程师也适合刚接触多目标优化、想找一个能直接复现的 Jupyter 实现的人。2. NSGA-II 的算法骨架非支配排序、拥挤度与精英保留2.1 为什么是 NSGA-II而不是加权求和或 MOEA/D很多人第一反应是“多目标不就是把几个目标加权加起来当单目标做吗”。加权求和的问题在于权重极难定而且它只能找到 Pareto 前沿上的凸部分遇到凹前沿直接失效。NSGA-II 的核心优势有三个一是用非支配排序给种群分层保证收敛到前沿二是用拥挤度距离维持解的多样性避免全挤在一个点三是精英保留策略父代和子代合并后再选防止已经找到的好解丢失。这三点加起来让它在两到三个目标的优化问题上非常稳。常见做法是目标数不超过 3 时优先用 NSGA-II超过 3 个目标再考虑 NSGA-III 或 MOEA/D。2.2 非支配排序与拥挤度的计算逻辑非支配排序的本质是对每个个体统计有多少个体支配它 domination count 以及它支配了哪些个体。第一层是 domination count 为 0 的个体把它们移除后被它们支配的个体 count 减一再找出新的 0 层依次类推。拥挤度则是把同一层内的个体按每个目标排序取相邻个体在各自目标上的距离归一化后求和。边界个体的拥挤度设为无穷大保证它们一定被保留。def fast_non_dominated_sort(population): # population: list of individuals, each has .objectives (list of float) fronts [[]] for p in population: p.domination_count 0 p.dominated_solutions [] for q in population: if p.dominates(q): p.dominated_solutions.append(q) elif q.dominates(p): p.domination_count 1 if p.domination_count 0: p.rank 0 fronts[0].append(p) i 0 while fronts[i]: next_front [] for p in fronts[i]: for q in p.dominated_solutions: q.domination_count - 1 if q.domination_count 0: q.rank i 1 next_front.append(q) i 1 fronts.append(next_front) return fronts[:-1] # 最后一层为空这段代码里dominates方法需要你自己在个体类里实现判断条件是对所有目标p 不差于 q且至少有一个目标严格优于 q。domination_count和dominated_solutions是挂在个体对象上的临时属性每次排序前要重置。参数上唯一要注意的是 fronts 的初始化fronts [[]]后面 append 空列表会导致多一层所以返回时用fronts[:-1]去掉。2.3 选择、交叉、变异在 NSGA-II 里的具体接法NSGA-II 的选择用二元锦标赛但比较规则是先比 rankrank 小的赢rank 相同比拥挤度拥挤度大的赢。交叉和变异就是常规的模拟二进制交叉SBX和多项式变异。SBX 的分布指数eta_c一般取 15 到 20多项式变异的eta_m取 20 左右变异概率p_m通常取1/n_vars。这些参数不是玄学eta越大子代越靠近父代探索能力越弱太小又退化成随机搜索。def tournament_selection(pop, k2): candidates random.sample(pop, k) best candidates[0] for c in candidates[1:]: if c.rank best.rank: best c elif c.rank best.rank and c.crowding_distance best.crowding_distance: best c return best锦标赛大小k默认 2 就够k 越大选择压力越大收敛快但容易早熟。拥挤度在边界个体上是float(inf)比较时不会出错但如果你用 numpy 的 inf 要注意别参与后续的归一化计算。3. 在 Jupyter 里把 NSGA-II 跑起来从环境到第一个 Pareto 前沿3.1 Jupyter 环境准备与依赖安装Jupyter Notebook 的安装方式现在主流是走 Anaconda 或者 pip。如果你用 Anacondaconda install jupyter notebook就行如果走 pippip install notebook也可以。Python 版本建议 3.9 以上numpy 和 matplotlib 是必装的。我一般会在项目目录下建一个虚拟环境避免和系统 Python 打架。python -m venv nsga2_env source nsga2_env/bin/activate # Windows 用 nsga2_env\Scripts\activate pip install numpy matplotlib notebook jupyter notebook启动后浏览器会自动打开如果没打开终端里会有一行带 token 的地址复制到浏览器即可。Jupyter 默认存放地址是启动时所在的目录如果你想让它固定在某个文件夹可以在启动前cd过去或者改配置文件里的notebook_dir。常见坑是 Anaconda 的 Jupyter 突然打不开多半是端口被占或者配置文件被改坏了用jupyter notebook --port 8889换个端口试试。3.2 用 Jupyter 组织 NSGA-II 代码的推荐结构一个能复现的 NSGA-II 实现我建议在 Jupyter 里按 cell 分块第一个 cell 放导入和全局参数第二个 cell 定义个体类和支配判断第三个 cell 写非支配排序和拥挤度第四个 cell 写交叉变异第五个 cell 写主循环第六个 cell 做可视化。这样调试的时候可以单独跑某个 cell不用每次从头执行。Jupyter 一个 cell 只输出最后一个表达式的结果如果你要同时看多个变量用print或者display。import numpy as np import matplotlib.pyplot as plt import random # 全局参数 POP_SIZE 100 N_GEN 200 N_VARS 2 ETA_C 20 ETA_M 20 P_M 1.0 / N_VARS BOUNDS [(-5, 5), (-5, 5)]POP_SIZE太小前沿覆盖不全太大跑得慢100 到 200 是常见起点。N_GEN看问题复杂度简单问题 100 代就收敛复杂问题可能要 500 代以上。BOUNDS是决策变量的上下界SBX 和多项式变异都要用。3.3 主循环与结果可视化主循环的逻辑是初始化种群评估目标非支配排序加拥挤度然后循环选择、交叉、变异生成子代父子合并再排序选前 N 个。可视化用 matplotlib 画散点图横纵轴分别是两个目标值。def nsga2_main(): pop [Individual(random.uniform(*BOUNDS[i]) for i in range(N_VARS)) for _ in range(POP_SIZE)] for ind in pop: ind.evaluate() fronts fast_non_dominated_sort(pop) assign_crowding_distance(fronts) for gen in range(N_GEN): offspring [] while len(offspring) POP_SIZE: p1 tournament_selection(pop) p2 tournament_selection(pop) c1, c2 sbx_crossover(p1, p2, ETA_C, BOUNDS) mutate(c1, ETA_M, P_M, BOUNDS) mutate(c2, ETA_M, P_M, BOUNDS) c1.evaluate(); c2.evaluate() offspring.extend([c1, c2]) combined pop offspring[:POP_SIZE] fronts fast_non_dominated_sort(combined) assign_crowding_distance(fronts) pop [] for front in fronts: if len(pop) len(front) POP_SIZE: pop.extend(front) else: front.sort(keylambda x: x.crowding_distance, reverseTrue) pop.extend(front[:POP_SIZE - len(pop)]) break return popcombined的长度是 2 倍POP_SIZE选前POP_SIZE个的时候如果某一层放不下就按拥挤度从大到小截断。这个截断逻辑是 NSGA-II 精英保留的关键写错会导致种群数量不稳定。跑完后用plt.scatter把第一层前沿画出来横纵轴分别是f1和f2。4. 参数调优与常见翻车现场排查4.1 交叉变异参数怎么设才不跑偏SBX 的eta_c控制子代和父代的接近程度值越大子代越像父代。我一般从 15 开始试如果前沿收敛太慢就降到 10如果种群多样性不够就升到 25。多项式变异的eta_m类似20 是默认值。变异概率p_m用1/n_vars是经验公式变量多的时候自动降低变异强度。注意 SBX 和多项式变异都要处理边界超出BOUNDS的变量要截断或者反射回来否则目标函数可能报错。4.2 非支配排序结果不对的三种典型原因第一种是dominates方法写反了把“不差于”写成了“严格优于”导致支配关系全错。第二种是排序前没重置domination_count第二次调用时 count 还是上一轮的值。第三种是 fronts 初始化多了一层空列表导致 rank 偏移。排查方法很简单拿两个明显有支配关系的个体手动跑一遍dominates看返回值对不对然后在排序函数里打印每层的个体数量正常应该是递减的。4.3 Jupyter 里跑大规模种群的性能问题Python 原生循环做非支配排序是 O(MN²)种群 500 以上、代数 500 以上就会明显卡。Jupyter 里可以用%%time魔法命令测每个 cell 的耗时。优化方向有两个一是用 numpy 向量化支配判断把目标值存成矩阵用广播比较二是换用 pymoo 这类成熟库它底层用 numpy 优化过。如果只是学习算法流程原生实现够用如果要上生产建议直接看 pymoo 的 NSGA-II 实现。5. 避坑与常见问题那些让我重跑一晚上的细节5.1 现象前沿点全挤在一起多样性极差原因拥挤度计算时没有归一化或者边界个体的拥挤度没设成无穷大。解决在assign_crowding_distance里对每个目标先取 min 和 max归一化后再算相邻距离边界个体直接赋float(inf)。5.2 现象种群数量每代都在变最后只剩几个解原因精英保留的截断逻辑写错了比如pop.extend(front)之后没有 break或者截断时用了front[:POP_SIZE]而不是front[:POP_SIZE - len(pop)]。解决严格按“能放整层就放整层放不下就按拥挤度截断”的逻辑写截断后立即 break。5.3 现象目标函数报数组维度错误原因决策变量是 list目标函数里直接当 numpy 数组用或者evaluate返回的不是标量。解决在evaluate里把变量转成 numpy 数组再算返回值确保是 float。Jupyter 里可以用%debug进事后调试看具体是哪一维对不上。5.4 现象跑了几十代前沿不再改善原因变异概率太低或者eta_m太大种群失去探索能力。解决把p_m临时调大或者给eta_m加个退火策略前期小后期大。也可以每隔几代随机重置一小部分个体。5.5 现象Jupyter 内核频繁挂掉原因种群太大内存爆了或者某个 cell 里有死循环。解决先用小种群跑通流程确认逻辑无误再放大。Jupyter 里可以用top或者任务管理器看内存占用超过 80% 就要考虑减种群或减代数。6. 进阶技巧用 pymoo 交叉验证你的实现以及一个收敛判断的小习惯自己手写 NSGA-II 最大的价值是理解流程但验证结果对不对我一般会拿 pymoo 跑同一个问题做对比。pymoo 的 API 很简洁定义问题、算法、终止条件三行就能跑。如果你的前沿和 pymoo 的前沿基本重合说明实现没问题如果差很多优先检查支配判断和拥挤度。from pymoo.algorithms.moo.nsga2 import NSGA2 from pymoo.problems import get_problem from pymoo.optimize import minimize problem get_problem(zdt1) algorithm NSGA2(pop_size100) res minimize(problem, algorithm, (n_gen, 200), seed1, verboseFalse) print(res.F.shape) # Pareto 前沿的目标值zdt1是标准测试问题前沿是凸的适合入门验证。res.F是前沿上的目标值矩阵行数是解的数量。跑完后可以和你的实现画在同一张图上对比。另一个习惯是每跑完一次把第一层前沿的f1最小值、f2最小值记下来看随代数怎么变。如果连续 20 代这两个值都不再下降基本可以停。这个判断比固定代数更省时间也避免跑过头。Jupyter 里可以写个简单的回调函数每 10 代打印一次当前前沿的极值。最后说个血泪教训别在 Jupyter 里直接跑 1000 代不存中间结果内核一挂全没了。我现在的习惯是每 50 代把种群的目标值存成 npy 文件文件名带代数这样即使后面崩了也能从最近一次恢复。希望帮到你。本文还有配套的精品资源点击获取