约束多目标优化问题[1] (constrained multi-objective optimization problem, CMOP)广泛存在于科学研究和工程应用中,不仅涉及多个冲突目标函数的优化,还需满足多种约束条件,是一类复杂的优化问题,对现有求解方法提出了严峻挑战。
为了求解CMOP,学术界通过结合多目标进化算法和约束处理技术设计了多种约束多目标进化算法(constrained multi-objective evolutionary algorithm,CMOEA)[2-3],其中主流的算法包括基于多阶段的CMOEA[4-5]、基于多种群的CMOEA[6-7]和基于多任务的CMOEA[8-9]。基于多任务的CMOEA通过设计简化的辅助任务,降低了约束或目标的求解难度[8,10]。进一步,算法从辅助任务向主任务(即原CMOP)迁移有效的知识,帮助求解主任务。例如,Qiao等[8]提出了第一个基于多任务的CMOEA,创建了任意数量的辅助种群,并通过知识迁移提升每个种群的性能。进一步,Qiao等[9]提出了基于动态辅助任务的方法,该辅助种群动态缩放约束边界,提升了和主种群的相似性,从而提升了知识迁移的有效性。近期,多位学者提出了更多类型的辅助任务和知识定义方法[11-12],丰富了该类方法的理论,提升了该类方法的有效性。
然而,现有基于多任务的方法忽视了调节种群间的竞争关系,采用均质化计算资源分配策略,导致进化后期低效种群仍产生冗余个体,显著降低了收敛效率。此外,种群间信息交互局限于子代分享,子代生成过程缺乏跨种群协作,导致高质量个体的信息未能被充分利用。
为了克服以上缺陷,本文提出了一种基于竞争与合作多任务的约束多目标进化算法(competition and cooperation multitasking based constrained multi-objective evolutionary algorithm, CCMT),其包含3种任务及其对应的进化种群。主任务即原CMOP,其种群优先保存可行帕累托非支配解;辅助任务1为删除约束的无约束多目标优化问题,其种群忽视约束以加强目标信息的利用,从而加速算法的收敛速度;辅助任务2为约束边界动态缩放的CMOP,其种群动态利用有效的不可行解信息,以维持算法的多样性。此外,本文设计了一种基于竞争的资源分配策略,通过测评各任务种群的历史表现,动态调整其在随后进化过程中所使用的计算资源,以适合当前地形的任务发挥更大的价值。进一步,本文还设计了基于父代聚合和子代扩散的协同优化策略,3个种群在子代产生和环境选择过程中通过合作来提高信息利用效率。在38个标准测试函数上的实验结果表明,所提算法明显优于5种现有的较先进的算法。
不失一般性,一个CMOP可以用以下公式表示[13]:
(1)
式中:f(x) 代表单个优化目标;F(x)代表一个M维的优化目标向量;x代表一个D维决策向量;SD代表决策空间;g(x)和h(x)分别代表不等式约束和等式约束;p和(n-p)分别代表不等式约束和等式约束的数量。
为了计算等式约束的违反度,首先,一般通过松弛参数δ将等式约束转换为不等式约束[14]:
|hi(x)|-δ≤0,i=p+1,p+2,…,n。
(2)
其次,每个约束的违反度可以由下式计算:
(3)
最后,x的总约束违反度[14]可以由下式计算:
(4)
当G(x)等于0时,x满足所有约束,被称为可行解,否则被称为不可行解。给定两个可行解x1和x2,对任意j∈(1,2,…,M), fj(x1)≤fj(x2),并且存在k∈(1,2,…,M),使得fk(x1)<fk(x2),则称解x1帕累托支配解x2。若x*不被任何其他解支配,则称其为帕累托最优解[15]。求解CMOP的目的是获得一组可行的帕累托最优解集(constrained Pareto set, CPS),其在目标空间中的映射向量被称为约束帕累托前沿(constrained Pareto front, CPF)。若不考虑约束,CPF将变为无约束帕累托前沿(unconstrained PF,UPF)[16]。
图1为算法CCMT的框架。如图1所示,算法CCMT框架中包含一个主任务(表示为T1)和两个辅助任务(分别表示为T2和T3),并创建3个种群分别求解3个任务(分别表示为P1、P2和P3)。第一个任务T1为用式(1)表示的原CMOP,该任务被称为主任务,需要算法找到可行帕累托最优解。因此,种群P1使用约束支配准则(constrained dominance principle, CDP)[15]更新种群,优先保存约束值小的解,使得P1可以快速定位到可行域并保存可行解。
图1 算法CCMT的框架
Figure 1 Framework of the algorithm CCMT
第二个任务T2为无约束辅助任务[6],其公式为
min F(x)=(f1(x),f2(x),…,fM(x))。
(5)
该辅助任务删除了所有约束,最优解为UPF。通过舍弃约束,该辅助任务能够帮助主任务种群快速探索整个目标空间,从而快速找到有潜力的搜索区域。因为T2为无约束多目标任务,种群P2使用无约束非支配排序方法[15]更新种群。
第3个任务T3为约束缩放任务[11],公式为
(6)
式中:ε代表约束缩放边界,是一个动态变化的值,在本文中,ε等于后代种群中所有不可行解的平均约束值。为了更有效地探索约束边界内的区域,种群P3使用epsilon方法[17]更新种群,其中约束值小于ε的解使用基于多目标的约束处理方法[18]进行排序;约束值大于ε的解使用CDP方法进行排序。
这3个任务具有不同的最优解,使得3个种群有不同的进化轨迹,因此,执行知识迁移可以提高种群多样性。此外,两个辅助任务分别利用全局解和缩放区域中的解,对主任务具有互补的辅助效果。
CCMT算法提供了详细的伪代码,如算法1所示。首先,为每个任务初始化一个大小为NP的特定种群。此外,初始化奖励矩阵R和每个任务的选择概率s1、s2和s3。随后,算法进入主循环。第5行,执行轮盘赌操作确定选择第k∈{1,2,3}个任务作为进化任务。第6行,第k个任务执行父代聚合操作产生一个杂交父代种群HP。第7行,HP使用差分进化操作[9] (DE/current-to-rand/1和DE/current-to-pbest/1)产生后代种群O。第8行,执行后代扩散策略完成环境选择过程。最后,更新3个任务的选择概率。
算法1: CCMT算法。
输入: 种群规模NP、父代聚合概率ρ;
输出: P1。
① 初始化种群P1、P2、P3,每个种群大小为NP;
② 初始化奖励矩阵R和选择概率(s1,s2,s3);
③ t=0;
④ While没有达到最大迭代次数
⑤ k←轮盘赌(s1,s2,s3);
⑥ HP←父代聚合(P1,P2,P3,k,ρ)→算法2;
⑦ O←HP使用差分进化操作产生NP个后代个体;
⑧ P1,P2,P3←后代扩散(P1,P2,P3,O) →算法3;
⑨ t=t + 1;
⑩ For k=1 : 3
和
根据式(8)、式(9)计算第k个任务的全局奖励和种群奖励;
Rk,t←根据式(7)计算第k个任务的奖励;
End
根据式(10)更新(s1,s2,s3);
End While
现有基于多任务的算法中,每个种群被分配相同的计算资源,忽略了种群在不同进化阶段适应性的动态变化,因此,有必要调整种群间的资源分配[11]。图2提供了一个小可行域问题[19]上种群在进化后期的分布。P2通过忽视约束定位在无约束帕累托前沿上,P3通过约束缩放定位在约束搜索边界内部,而P1通过优化原问题定位在局部约束帕累托前沿。如图2所示,无约束帕累托前沿和约束帕累托前沿相距较远,导致P2的更新无法继续为寻找约束帕累托前沿提供有效的信息。因此,有必要实时调整种群间的计算资源。基于此,本节提出了一种基于竞争的资源分配策略。
图2 小可行域问题上3个种群在进化后期的分布
Figure 2 Distribution of three populations at the late evolutionary phase on a problem with small feasible regions
如算法1所示,CCMT基于任务的选择概率使用轮盘赌方法选择有效的进化种群,其中选择概率由任务种群在过去LR代的累计奖励计算得来。
第t代中第k个任务的奖励Rk,t为
(7)
式中:
为全局奖励,反映任务对全局最优的贡献程度;
为种群奖励,反映第k个种群的多样性和收敛性的改进速度;参数α控制两种奖励的权重。如果一个任务有助于更频繁地更新全局最优值,或者促进其种群内个体的更快进化,它将会获得更高的奖励。
和
的具体计算过程为
(8)
(9)
式中:
表示直到第t代所有可行解在每个目标上的最小值;
表示前者支配后者。如果sign函数中的布尔值为真,则返回值为1,否则返回值为0;num_O代表在环境选择中成功进入下一代的后代个体数量,num_O越大,代表对应种群进化得越快。
基于上述定义,每个任务的选择概率可用式(10)更新:
(10)
式中:T为最大进化代数;参数β控制选择概率的更新方式。在进化前期,每个任务有相同的选择概率,使得每个任务都能被选中以评估其表现;在进化后期,种群在历史过程中的表现奖励被用于更新选择概率,使得适合当前环境的任务能被更频繁地选中,从而发挥此任务种群的效果,避免低效种群浪费计算资源。
种群协同优化包含两个步骤:基于父代聚合的交配父代选择和基于后代扩散的环境选择。
2.3.1 基于父代聚合的交配父代选择
现有方法仅在环境选择过程中分享子代个体完成信息分享,然而,这种方式未能充分利用不同种群的信息产生高质量的杂交子代[20]。如图2所示,小可行域导致P1收敛至局部约束帕累托前沿上。虚线箭头是由P2的个体和P3的个体产生的搜索方向,由该搜索方向产生的子代个体可以定位于小可行域上。因此,基于种群合作的子代产生方式能提供更多的搜索多样性,有助于种群跳出局部最优。基于此,本节设计了基于父代聚合的交配父代选择策略。
算法2提供了交配父代选择的伪代码,其目的是从3个种群中选择包含NP个父代个体的交配池HP。首先,当前进化种群Pk由算法1的第5行决定。对于种群Pk的第i个个体xi,k,若随机数rand小于ρ,则它与其他种群的个体进行比较,从而确定能够进入交配池的个体。具体步骤为从另外两个种群中随机选择一个种群Ph,并从中随机选择一个个体xr,h。然后使用CDP方法比较xi,k和xr,h,其中较好的个体被添加进集合HP中。若随机数rand大于等于ρ,直接将xi,k添加进HP中。
此方法中,基于概率的方法能从另外两个种群中选择有效的个体参与后代产生,这有助于信息交流。此外,基于CDP的比较方法使得另外两个种群中约束违反度较低的解参与后代产生,能够帮助种群更快收敛至可行域。因此,基于父代聚合的交配父代选择方法能够维持种群间的合作并且加快收敛。
算法2: 父代聚合算法。
输入: 种群规模NP、父代聚合概率ρ、种群P1、种群P2、种群P3、被选中种群的索引k、聚合概率β;
输出: 交配父代种群HP。
① 设置HP等于空集;
② for i=1 : NP
③ 选择Pk中第i个个体,记作xi,k;
④ if rand<ρ
⑤ 随机选择一个任务Th;
⑥ 从Ph中随机选择一个个体xr,h;
⑦ if xr,h
CDP xi,k
⑧ x*←xr,h;
⑨ elseif xi,k
CDP xr,h
⑩ x*←xi,k;
else
x*←从xr,h和xi,k中随机选择一个个体;
End if
else
x*←xi,k;
End if
HP←x*;
End for
2.3.2 基于后代扩散的环境选择
算法3提供了该方法的伪代码,主要目的是使用杂交后代种群O更新3个任务种群。3个任务种群分别通过CDP方法、无约束非支配排序方法和epsilon方法完成环境选择过程。
算法3: 后代扩散算法。
输入: 种群规模NP、种群P1、种群P2、种群P3、后代种群O;
输出: P1、P2、P3。
① TP1←P1∪O;
② P1←使用CDP方法从TP1中选择NP个个体;
③ TP2←P2∪O;
④ P2←使用无约束非支配排序方法从TP2中选择NP个个体;
⑤ TP3←P3∪O;
⑥ P3←使用epsilon方法从TP3中选择NP个个体。
CCMT的时间复杂度主要来源于父代聚合操作、子代产生操作和后代扩散操作。前两种操作的时间复杂度分别为O(NP)和O(NP·D)。后代扩散操作中,CDP、无约束非支配选择和改进epsilon方法的时间复杂度均为O(M·NP2)。因此,CCMT的时间复杂度为O(NP)+O(NP·D)+O(M·NP2)=O(M·NP2),与大多数约束多目标进化算法的时间复杂度相同[6-7]。
对比方法。本文选择了几个最新的CMOEA作为对比方法,包括CMOEA_MS[5]、cDPEA[7]、EMCMO[8]、MTCMO[9]和CMOEMT[11]。CMOEA_MS为基于多阶段策略的方法,学习约束的优先级,并逐渐合并有效的约束条件以利用多样化的不可行信息。cDPEA是基于多种群策略的算法,同时使用CDP和惩罚函数方法平衡算法的收敛性和多样性。最后3个为基于多任务的算法,但他们设计了不同的辅助任务和知识迁移策略,因此适用于求解不同类型的问题。这5个对比算法从单约束信息利用角度、约束处理技术角度、知识迁移策略角度和辅助任务角度设计了不同的策略,因此具有代表性并且被选为对比算法。
测试集。本文选择DASCMOP[21]、LIRCMOP[19]和SDC[22]作为测试集,它们分别包含9,14,15个测试函数,这3个测试集共包含38个测试函数。前两个测试集在目标空间包含了离散可行域和小可行域。SDC在决策空间包含了大量局部不可行域。这3个测试集具有不同的特性,以评估算法的综合性能。
实验参数设置。对于两目标和三目标问题,所有算法种群大小分别设置为100和120,最大函数评价次数设置为20万,每个算法的运行次数设置为30[23]。对于使用的差分进化算法操作,比例因子F和杂交率CR分别从集合{0.6, 0.8, 1.0}和集合{0.1, 0.2, 1.0}中随机选择。此外,参数β设置为0.1,参数α设置为0.5,参数ρ设置为0.3,参数LR设置为10。对比算法的参数设置和原文保持一致。本文所有实验均在服务器上运行,其硬件环境配置为Inter (R) Xeon (R) Gold 6230 CPU@2.10 GHz 640 GB RAM。
评价指标。使用的评价指标为反向世代距离IGD和超体积HV。IGD计算最优解和算法输出可行非支配解集的平均距离,而HV计算参考点和算法输出可行非支配解集围成的目标空间中区域的体积。IGD越小或者HV越大,代表算法的收敛性和多样性越好。
表1提供了6种算法在LIRCMOP测试集上的IGD结果,其中每个问题上的最优结果用灰色表示。表2提供了所有算法在DASCMOP和SDC测试集上的统计IGD结果。图3提供了所有算法在3个测试集上的HV的对比情况。
表1 6种算法在LIRCMOP测试集上的IGD
Table 1 IGD of six algorithms on LIRCMOP benchmark
问题IGD 均值 (IGD标准差)CMOEA_MScDPEAEMCMOMTCMOCMOEMTCCMTLIRCMOP13.0385e-1(1.26e-1)-2.6626e-2(8.99e-3)+1.2591e-1(7.62e-2)-2.4149e-2(9.16e-3)+2.5455e-2(2.68e-2)+4.4518e-2(1.98e-2)LIRCMOP21.9266e-1(9.40e-2)-6.3464e-3(4.91e-4)+1.1472e-1(6.42e-2)-1.2825e-2(5.81e-3)=3.6101e-2(4.79e-2)=1.6459e-2(1.14e-2)LIRCMOP33.2548e-1(1.68e-1)-7.9127e-3(5.33e-3)+1.3209e-1(7.42e-2)-2.7255e-2(2.05e-2)+4.6272e-2(6.69e-2)-3.7920e-2(1.54e-2)LIRCMOP42.7296e-1(7.33e-2)-6.1471e-3(3.78e-3)+1.5981e-1(5.45e-2)-2.7872e-2(2.07e-2)+6.5701e-2(8.80e-2)=3.4114e-2(1.53e-2)LIRCMOP54.2658e-1(5.74e-1)-1.2335e-2(5.50e-3)-8.6057e-3(8.96e-4)-7.6624e-2(2.58e-1)-1.3204e-2(8.77e-3)-6.8212e-3(3.97e-4)LIRCMOP63.9039e-1(5.84e-1)=7.3729e-3(5.45e-4)-7.7869e-3(1.03e-3)-8.8336e-2(3.02e-1)-1.1301e-2(4.94e-3)-6.1947e-3(2.36e-4)LIRCMOP73.3850e-2(5.63e-2)=8.1932e-3(1.88e-3)-7.5581e-3(2.17e-4)-1.3309e-2(1.39e-2)-1.0820e-2(5.08e-3)-7.3008e-3(3.64e-4)LIRCMOP81.6450e-2(2.80e-2)=1.3456e-2(2.62e-2)=7.5481e-3(2.54e-4)=7.6846e-3(4.73e-4)=8.6779e-3(7.31e-4)-7.4540e-3(4.60e-4)LIRCMOP94.9078e-1(8.74e-2)-1.9781e-1(4.67e-2)-2.3729e-2(2.63e-2)-3.2604e-1(8.23e-2)-8.4322e-2(5.59e-2)-8.6783e-3(2.96e-3)LIRCMOP101.7763e-1(9.98e-2)-1.3828e-2(2.02e-2)-7.9222e-3(7.62e-4)-8.7950e-3(1.46e-3)-1.1255e-2(3.28e-3)-5.6205e-3(2.64e-4)LIRCMOP111.6427e-1(1.11e-1)-5.2423e-2(3.90e-2)-4.3364e-3(2.97e-3)-1.6867e-1(1.43e-1)-5.2631e-3(1.85e-3)-2.4632e-3(4.92e-5)LIRCMOP122.2497e-1(8.64e-2)-1.3314e-1(3.04e-2)-4.8508e-3(1.14e-3)-1.8514e-1(4.58e-2)-3.2491e-2(3.17e-2)-3.1923e-3(5.53e-5)LIRCMOP131.0913e-1(1.91e-3)-1.0982e-1(2.00e-3)-1.0438e-1(1.84e-3)-8.7360e-1(3.36e-1)-1.0222e-1(6.49e-3)-9.3323e-2(1.70e-3)LIRCMOP149.2574e-2(9.16e-4)-9.4087e-2(1.36e-3)-9.1582e-2(1.16e-3)-8.6899e-1(3.80e-1)-9.1673e-2(1.56e-3)-8.7767e-2(5.98e-4)
注:括号后的符号“+”、“-”及“ =”分别表示对比算法显著好于、明显差于和相似于CCMT,下同。
表2 在DASCMOP和SDC测试集上,基于IGD指标对比算法优于、差于或等于CCMT的函数数量
Table 2 On the DASCMOP and SDC benchmarks, the number of functions with algorithm better than, worse than, or equal to CCMT based on the IGD metric
对比算法函数数量测试集DASCMOP测试集SDC+-=+-=CMOEA_MS081294cDPEA261294EMCMO180294MTCMO171474CMOEMT180591
图3 CCMT与其他算法在3个测试集上的HV的对比情况
Figure 3 Comparison of HV between CCMT and other algorithms on three test sets
LIRCMOP和DASCMOP包含单模态特性的目标空间约束,这需要算法首先利用UPF信息快速逼近可行域,然后利用UPF和CPF中间的不可行解维持多样性并找到更多的可行域。表1显示CCMT在大部分LIRCMOP函数上取得了最优结果,这得益于CCMT中两个特性互补的辅助任务。对比算法中,只有cDPEA和MTCMO分别在3个和1个函数上取得了最优结果。LIRCMOP1~LIRCMOP4包含了非常小的可行域,需要算法利用多样化的不可行解,从而推动种群从不可行域逼近可行域。EMCMO仅包含了无约束辅助任务,因此在大多数问题上取得了较差的表现。MTCMO仅包含了使用epsilon方法的辅助任务以维持种群多样性,因此在部分问题上取得了较好的表现。虽然CMOEMT同时包含了无约束辅助任务和约束缩放辅助任务,但是种群独立产生后代,不充分的信息迁移使得算法产生了较差的结果。表2显示CCMT在DASCMOP测试集中的至少6个问题上显著优于5个对比算法,表明了所提算法的优异性。
SDC的决策空间包含了大量局部不可行域,导致算法容易陷入局部最优。所有方法中,CMOEA_MS、EMCMO和CMOEMT在3个、1个和1个函数上未能取得100%的可行率。CCMT和MTCMO在8个和5个函数上取得了最优结果,且表2的统计结果显示CCMT至少在7个问题上显著优于5个对比算法,表明了所提方法的优异性。
经3个测试集共38个测试函数上的Wilcoxon 秩和检验,在IGD指标上,CCMT优于对比算法CMOEA_MS、cDPEA、EMCMO、MTCMO、CMOEMT的函数数量分别为28,24,30,23,28个,而差于5个对比算法的函数数量分别为2,7,3,8,7个。此外, 图3显示CCMT在HV指标上依然取得了优异的表现,它优于5个对比算法的函数数量分别为32,28,30,24,28个,而差于5个对比算法的函数数量分别为2,5,4,6,4个。总之,IGD和HV结果显示了所提方法在3个测试集上显著优于对比算法。
图4提供了所有算法在所有测试函数上的统计对比结果,其中R+和R-分别代表CCMT和对比算法的表现,符号“*”代表两个算法之间存在显著差异。对于每个对比算法,R+明显大于R-,并且除了MTCMO在IGD指标上与CCMT未存在显著差异,其余所有算法在两个指标上均与CCMT存在显著差异,证明CCMT显著优于对比算法。图5提供了每个算法基于Friedman检验的平均排名,CCMT取得了最优的排名。
图4 统计性检验结果
Figure 4 Statistical test results
图5 不同算法的IGD和HV排名
Figure 5 IGD and HV rankings of different algorithms
本节进行了参数分析实验,对于4个参数,每个参数选取了5个不同的值,参数敏感性分析如表3所示。
表3 参数敏感性分析
Table 3 Parameter sensitive analysis
LR排名α排名β排名ρ排名12.602.603.203.252.80.32.50.12.20.12.9102.20.52.10.22.60.32.4503.10.83.50.33.40.52.71003.41.03.60.53.60.72.8
参数LR决定了利用历史信息的代数,LR越大,种群利用历史信息的程度越深,从而使得种群偏向于全局探索。参数α控制全局奖励和种群奖励的比例,且α越小使得种群越偏向于全局探索。因此,LR和α共同控制了算法全局探索和局部开发的平衡。表3显示两个参数的最优值均为5个参数的中间值,表明算法赋予了全局探索和局部开发相似的权重,使得种群既能寻找更优的解,同时也能维持多样性。
参数β控制何时开始竞争阶段。太小的β值使得奖励矩阵未能记录有效信息;而较大的β值使得资源分配未能被充分执行。因此,由表3可知,最优β值是一个较小的值。
参数ρ控制父代聚合概率。较大的值使得其他种群深度参与子代产生,提供的多样性可能会影响种群的进化方向;而太小的值则会抑制任务间的信息交流,导致有效信息无法被利用。因此,由表3可知,最优ρ值为较小的0.3。
CCMT中包括两个主要创新点:资源分配和种群协同进化。为了验证资源分配的有效性,设计了无资源分配的算法变种CCMT-noRA,其中每个任务的选择概率在进化过程中保持相同。为了验证种群协同进化的有效性,设计了两个算法变种:第一个是无父代聚合的算法变种CCMT-noPA,第二个是无子代聚合的算法变种CCMT-noOD。算法变种和CCMT的统计对比IGD结果如表4所示。在IGD指标上,CCMT优于对比变种算法CCMT-noRA、CCMT-noPA、CCMT-noOD的函数数量为21,25,33个,证明了所提策略的有效性。CCMT-noRA优于CCMT的函数数量为5个,这表明资源分配策略能够在不同的问题和不同的进化阶段为更有效的任务分配更多的计算资源。CCMT-noPA优于CCMT的函数数量为3个,CCMT-noOD优于CCMT的函数数量为0个。上述结果表明在子代产生和环境选择过程中的知识迁移均是有效的。
表4 基于IGD指标,变种算法优于、差于或等于CCMT的函数数量
Table 4 The number of functions with the variant algorithm better than, worse than, or equal to CCMT based on the IGD metric
变种算法函数数量+-=CCMT-noRA52112CCMT-noPA32510CCMT-noOD0335
本文提出了一种基于竞争与合作多任务的约束多目标进化算法CCMT。该算法中设置了一个主任务和两个互补的辅助任务,同时设计了基于竞争的资源分配策略和基于合作的种群协同进化策略。所提方法在3个不同特性的测试集上和5个最新的算法进行了对比,实验结果表明所提算法寻找可行解的能力较强,且多样性和收敛性明显优于对比算法,能够为决策者提供多样化的可行解。
虽然所提方法相比于5个最新的对比方法取得了显著优势,但该方法在具有非常小可行域的函数上未能取得最佳结果,因此未来将继续研究如何提升算法在该类问题上的效果,例如可以通过适应度地形分析技术提取小可行域的地形特征,从而指导种群在小可行域进行更为细致的搜寻。此外,基于多任务的约束多目标优化仍然是一个比较新的研究方向,目前研究主要关注于辅助任务的设计,在知识定义、知识迁移、任务内与任务间信息利用平衡等方面还需要进行深入研究。
[1] 李二超, 李进. 两阶段三存档集约束优化算法(TSDA)[J]. 郑州大学学报(工学版), 2018, 39(6): 23-29.LI E C, LI J. Constraint optimization algorithm with two-stage and three-archive (TSDA)[J]. Journal of Zhengzhou University (Engineering Science), 2018, 39(6): 23-29.
[2] 陈少淼, 陈瑞, 梁伟, 等. 面向复杂约束优化问题的进化算法综述 [J]. 软件学报, 2023, 34(2): 565-581.CHEN S M, CHEN R, LIANG W, et al. Overview of evolutionary algorithms for complex constrained optimization problems [J]. Journal of Software, 2023, 34(2): 565-581.
[3] 王林锋, 揭丽琳, 黎明, 等. 基于自适应双阶段分级均衡的约束高维多目标进化算法 [J]. 控制与决策, 2024, 40(5): 1512-1522.WANG L F, JIE L L, LI M, et al. Constrained many-objective evolutionary algorithm based on adaptive two-stage hierarchical equilibrium [J]. Control and Decision, 2024, 40(5): 1512-1522.
[4] FAN Z, LI W, CAI X, et al. Push and pull search for solving constrained multi-objective optimization problems [J]. Swarm and Evolutionary Computation, 2019, 44: 665-679.
[5] TIAN Y, ZHANG Y, SU Y, et al. Balancing objective optimization and constraint satisfaction in constrained evolutionary multiobjective optimization [J]. IEEE Transactions on Cybernetics, 2021, 52(9): 9559-9572.
[6] TIAN Y, ZHANG T, XIAO J, et al. A coevolutionary framework for constrained multiobjective optimization problems [J]. IEEE Transactions on Evolutionary Computation, 2020, 25(1): 102-116.
[7] MING M, TRIVEDI A, WANG R, et al. A dual-population-based evolutionary algorithm for constrained multiobjective optimization [J]. IEEE Transactions on Evolutionary Computation, 2021, 25(4): 739-753.
[8] QIAO K J, YU K J, QU B Y, et al. An evolutionary multitasking optimization framework for constrained multiobjective optimization problems [J]. IEEE Transactions on Evolutionary Computation, 2022, 26(2): 263-277.
[9] QIAO K J, YU K J, QU B Y, et al. Dynamic auxiliary task-based evolutionary multitasking for constrained multiobjective optimization [J]. IEEE Transactions on Evolutionary Computation, 2022, 27(3): 642-656.
[10] LI G, WANG Z, GAO W, et al. Decoupling constraint: task clone-based multi-tasking optimization for constrained multi-objective optimization [J]. IEEE Transactions on Evolutionary Computation, 2025, 29(2): 404-417.
[11] MING F, GONG W, WANG L, et al. Constrained multiobjective optimization via multitasking and knowledge transfer [J]. IEEE Transactions on Evolutionary Computation, 2022, 28(1): 77-89.
[12] YE Q, WANG W, LI G, et al. A self-organizing assisted multi-task algorithm for constrained multi-objective optimization problems [J]. Information Sciences, 2024, 664: 120339-120360.
[13] WANG J, LI Y, ZHANG Q, et al. Cooperative multiobjective evolutionary algorithm with propulsive population for constrained multiobjective optimization [J]. IEEE Transactions on Systems, Man, and Cybernetics: Systems, 2021, 52(6): 3476-3491.
[14] GARCíA J L L, MONROY R, HERN
NDEZ V A S, et al. COARSE-EMOA: An indicator-based evolutionary algorithm for solving equality constrained multi-objective optimization problems [J]. Swarm and Evolutionary Computation, 2021, 67: 100983-100997.
[15] DEB K, PRATAP A, AGARWAL S, et al. A fast and elitist multiobjective genetic algorithm: NSGA-Ⅱ[J]. IEEE Transactions on Evolutionary Computation, 2002, 6(2): 182-197.
[16] WANG C, LIU Z, QIU J, et al. Adaptive constraint handling technique selection for constrained multi-objective optimization [J]. Swarm and Evolutionary Computation, 2024, 86: 101488-101503.
[17] SONG S, ZHANG K, ZHANG L, et al. A dual-population algorithm based on self-adaptive epsilon method for constrained multi-objective optimization [J]. Information Sciences, 2024, 655: 119906-119923.
[18] LIU Z Z, QIN Y, SONG W, et al. Multiobjective-based constraint-handling technique for evolutionary constrained multiobjective optimization: a new perspective [J]. IEEE Transactions on Evolutionary Computation, 2022, 27(5): 1370-1384.
[19] FAN Z, LI W, CAI X, et al. An improved epsilon constraint-handling method in MOEA/D for CMOPs with large infeasible regions [J]. Soft Computing, 2019, 23: 12491-12510.
[20] MING M, WANG R, ISHIBUCHI H, et al. A novel dual-stage dual-population evolutionary algorithm for constrained multiobjective optimization [J]. IEEE Transactions on evolutionary Computation, 2021, 26(5): 1129-1143.
[21] FAN Z, LI W, CAI X, et al. Difficulty adjustable and scalable constrained multiobjective test problem toolkit [J]. Evolutionary Computation, 2020, 28(3): 339-378.
[22] QIAO K, LIANG J, YU K, et al. Evolutionary constrained multiobjective optimization: Scalable high-dimensional constraint benchmarks and algorithm [J]. IEEE Transactions on Evolutionary Computation, 2023, 28(4): 965-979.
[23] ZHONG X, YAO X, QIAO K, et al. A multitasking-based constrained multi-objective evolutionary algorithm with forward and backward stages [J]. IEEE Transactions on Emerging Topics in Computational Intelligence, 2024, 8(5): 3474-3488.