《A fast and elitist multiobjective genetic algorithm: NSGA-II》阅读笔记
原文:K. Deb, A. Pratap, S. Agarwal and T. Meyarivan, "A fast and elitist multiobjective genetic algorithm: NSGA-II"
这篇论文对好的个体评价标准是:1)支配性强;2)在一个不拥挤的区域
一、论文中变量的含义:
M——目标函数的个数
N——种群中个体的数量
p——某个个体(解决方案)
np——支配p的解决方案数量(也就是有多少个解决方案能支配p)
Sp——被p支配的解决方案的集合
L——非支配个体的集合
L[i].m——集合L中第i个个体的第m个目标函数值
fmmax——第m个目标函数的最大值
fmmin——第m个目标函数的最小值
二、算法:
(一)快速非支配排序
上图解释如下:
-
第一段代码:求出个体之间的支配关系 并 找出属于第一支配前沿的个体
遍历群体中的每个个体p,遍历的时候做两个工作:1)计算支配当前个体的解决方案数量np;2)找出被当前个体支配的个体,并将它放入集合Sp中。
对当前个体来说,如果群体中的其他个体都不能支配它(即np=0),那么这个个体就属于第一支配前沿,将其放入集合F1。
//M——目标函数的个数
//N——种群中个体的数量
//p——某个个体(解决方案)
//np——支配p的解决方案数量(也就是有多少个解决方案能支配p)
//Sp——被p支配的解决方案的集合
//P——种群集
//个体p
for each p∈P
Sp=∅
np=0
for each q∈P
if (p≺q) then//如果p支配q
Sp=Sp∪{q}
else if (q≺p) then//如果q支配p
np=np+1
if np=0 then
prank=1
F1=F1∪{p}
-
第二段代码:求其他个体属于哪一支配前沿
这很像一个拓扑图,每个个体作为一个顶点,如果我支配你,那么就在你我之间画一个有向箭头(我—>你),当把所有的箭头画完之后就出现了一个有向无环图。
第一段代码从没有前驱的顶点开始拓扑排序,求出了第一支配前沿F1。求出之后将这些点以及连接它的箭头从图中删掉,这时又出现了没有前驱的顶点,这些顶点就属于F2…一直这样拓扑排序下去。
i=1
while Fi≠∅
Q=∅
for each p∈Fi
for each q∈Sp
nq=nq−1
if nq=0 then
qrank=i+1
Q=Q∪{q}
i=i+1
Fi=Q
(二)拥挤距离计算
一般调用此函数时,参数都为属于同一个支配前沿的个体集合Fi。因为有m个目标函数,所以对于每个函数该个体的拥挤距离都不同,每一个个体的拥挤距离是在m个目标函数中拥挤距离的和。