Many-objective Evolutionary Algorithm Driven by Indicator Under Adaptive Reference Point Adjustment
-
摘要: 在具有不同Pareto前沿形状的优化问题上, 基于参考点的高维多目标进化算法表现出较差的通用性. 为了解决这个问题, 提出参考点自适应调整下评价指标驱动的高维多目标进化算法(Many-objective evolutionary algorithm driven by evaluation indicator under adaptive reference point adjustment, MaOEA-IAR). MaOEA-IAR提出Pareto前沿形状监测基础上的参考点自适应策略, 利用该策略选择一组候选解作为初始参考点; 然后通过曲线参数对参考点位置进行调整; 将最终得到的能够适应不同Pareto前沿的参考点用于计算增强的反世代距离指标, 基于指标值设计适应度函数作为选择标准. 实验证明提出的算法在处理各种Pareto前沿形状的优化问题时能获得较好的性能, 算法通用性高.
-
关键词:
- 参考点自适应 /
- 评价指标 /
- 高维多目标 /
- Pareto前沿形状
Abstract: Many-objective evolutionary algorithms based on reference points show poor versatility on optimization problems with different shapes of Pareto fronts. To address this issue, this paper proposes a many-objective evolutionary algorithm driven by evaluation indicator under adaptive reference point adjustment (MaOEA-IAR). MaOEA-IAR proposes a reference point adaptation strategy based on the Pareto front shape monitoring, and uses this strategy to select a group of candidate solutions as the initial reference points, then adjusts their positions with the curve parameter. The final obtained reference points that can adapt to different Pareto fronts are used to calculate the enhanced inverted generational distance indicator, the fitness function is designed as selection criterion based on the indicator value. The experiment shows that the algorithm proposed in this paper can get good performance and high versatility when dealing with the optimization problems with various shapes of Pareto fronts. -
无线传感网(Wireless sensor networks, WSNs)由负责感知和传输任务的传感器节点组成, 具有功耗低、成本小、易部署的优势, 被广泛应用于环境监控、工业控制等领域[1-3]. 时间同步是无线传感网的关键技术之一. 无线传感网的许多基础功能, 例如数据融合、传输调度和目标定位等, 都需要节点间的精确时间同步[4-6]. 因此, 对无线传感网中的时间同步技术进行研究和探索具有重要意义.
由于时钟的启动时刻和振荡器的变化特性不同, 传感器节点时钟之间通常存在两方面的偏差. 一方面是时钟的初始相位偏差, 称之为时钟偏移; 另一方面是时钟运行速率的偏差, 称之为时钟漂移. 无线传感网中的节点之间要维持高精度同步, 就需要利用时间同步技术计算相对时钟漂移和偏移. 如果仅对时钟漂移进行估计, 那么只能校正节点之间的时钟速率偏差, 而相位偏差会一直存在, 导致同步无法实现. 反之, 若只估计时钟偏移, 虽然能够校正相位偏差, 但是振荡器的差异会导致同步时间较短, 从而需要频繁地进行重同步. 因此, 为了避免上述情况的出现, 就必须对时钟漂移和偏移参数进行联合估计.
无线传感网时间同步机制的设计所面临的一项重要挑战是: 如何尽可能地减少同步能耗. 由于网络中节点间的同步通常需要传输同步信息, 而同步信息传输所需能耗占据同步能耗的绝大部分. 因此, 一个有效的解决方案是设计以最小化同步信息传输数量为目标的低功耗同步协议. 隐含同步是其中的一个典型协议[7], 它基于无线媒介的广播特性, 利用监听策略隐式地获取同步信息, 可以显著减少同步信息传输数量. 在该机制中, 一个活跃节点(既发送同步信息又接收同步信息)和一个时钟源节点之间执行双向信息交换同步操作, 而一些位于这两个节点重叠通信范围内的节点(称之为隐含节点), 只需监听它们之间的信息交换过程就能实现与时钟源节点的同步. 在此过程中, 隐含节点仅接收了信息, 没有发送任何信息, 大幅度地降低了同步能耗. 文献[7-9]基于隐含同步机制, 在不同的网络场景下, 利用统计信号处理技术设计了多种同步算法, 实现节点之间的高精度低能耗同步. 文献[10-11]则将隐含同步机制与另一种同步机制(校正式同步)相结合, 节点仅以少量的能量就能在估计时钟参数的同时校正自己的本地时钟, 实现实时的同步. 此外, 文献[12]还将隐含同步机制应用到了水下无线传感网中, 提出适用于水下传感器节点的低功耗同步算法.
免时间戳同步是近年来提出的另一种低能耗的同步机制[13-16]. 由于其交互过程无需时间戳, 同步功能可以无缝嵌入现有网络数据流, 从而能够显著减少能耗. 在免时间戳同步中, 首先活跃节点发送不含时间戳的数据包给时钟源节点, 并记录此时的本地时间; 接着时钟源节点接收到数据包后在预定义的响应时间间隔返回不含时间戳的数据包; 最后活跃节点在接收到返回的数据包后记录自己的本地时间, 并通过预定义的响应时间间隔规则估计时钟参数, 与时钟源节点达到同步. 该机制利用接收方对发送方的预定义响应时间来传递同步信息, 避免了专用同步帧的传输. 文献[13]提出一种基于和响应策略的免时间戳同步协议, 其中时钟源节点响应时间和接收时间的和满足特定的规则. 通过在两个不同的时间间隔各返回一个响应数据包, 文献[14]提出一种跟随响应免时间戳同步协议. 但是, 由于跟随响应数据包的存在, 限制了该协议在实际网络中的应用. 为解决上述问题, 文献[15]提出一种基于动态响应的免时间戳同步协议, 利用特定的映射规则, 将连续的两个或多个同步周期的响应时间设置为不同值, 消除了跟随响应数据包的传输. 进一步, 文献[16]提出一种免时间戳交互与单向传输混合的同步机制, 实现了时钟漂移和偏移的联合估计.
将免时间戳同步与隐含同步相结合, 能够联合发挥两种同步机制的低能耗优势, 从而进一步降低无线传感网同步能耗开销. 文献[15]对其进行了初步的研究, 但是, 在基于动态响应的免时间戳同步方法中, 隐含节点只能估计时钟漂移, 不能估计时钟偏移, 从而无法实现隐含节点的完全同步. 因此, 针对此问题, 本文提出了一种在不传递时间戳的情况下, 隐含节点能够联合估计时钟漂移和偏移的低功耗同步协议. 主要贡献如下:
1)提出了一种基于免时间戳交互的隐含同步协议, 隐含节点记录监听到数据包的本地时间戳, 再结合预定义的响应规则就能完成对时钟漂移和偏移的联合估计, 以少量的能耗达到与时钟源节点的完全同步;
2)针对典型的高斯随机时延, 推导了隐含节点时钟漂移和偏移的联合最大似然估计器(Maximum likelihood estimator, MLE), 以及相应的克拉美罗下限(Cramer-Rao lower bound, CRLB);
3)仿真结果表明, 所提估计算法能够有效实现对隐含节点时钟漂移和偏移参数的免时间戳联合估计, 并具有达到CRLB的优良性能.
1. 同步协议描述
本节将对所提的基于免时间戳交互的隐含同步协议进行介绍, 并建立节点的时钟模型.
1.1 交互过程
考虑一个由多个节点构成的无线传感网, 其中节点
$ O $ 是时钟源节点, 为网络提供参考时间; 节点$ P $ 是活跃节点, 周期性地与时钟源节点$ O $ 进行双向信息交互; 节点$ Q $ 是位于活跃节点$ P $ 与时钟源节点$ O $ 公共广播领域内的隐含节点, 能够监听它们之间的信息交互过程, 如图1所示.当网络同步开始, 活跃节点
$ P $ 与时钟源节点$ O $ 之间进行免时间戳同步, 隐含节点$ Q $ 监听两节点的交互信息进行隐含同步. 具体的同步过程如图2所示, 以第$ j $ 轮通信过程为例, 详细步骤如下.
图 2 节点$Q$ 监听节点$O$ 与节点$P$ 之间的$N$ 轮免时间戳交互Fig. 2 The$N$ rounds of timestamp-free exchange between node$O$ and$P$ with node$Q$ overhearing them1)活跃节点
$ P $ 在本地时间为$ t_{1,j}^{\left( P \right)} $ 时刻发送不含时间戳信息的数据包给时钟源节点$ O $ .2)时钟源节点
$ O $ 接收到该数据包时, 记录本地时间为$ t_{2,j}^{\left( O \right)} $ . 接着等待一段时间, 在$ t_{3,j}^{\left( O \right)} $ 时刻返回一个不含时间戳的确认帧(Acknowledgement, ACK)给节点$ P $ . 这里$ t_{3,j}^{\left( O \right)} = \xi t_{2,j}^{\left( O \right)}-{{R}_{j}} $ , 其中$ \xi $ 是一个大于1且接近于1的已知常数,$ {{R}_{j}} $ 是一个等待间隔补偿.3)活跃节点
$ P $ 接收到来自时钟源节点$ O $ 的ACK时, 记录本地时间为$ t_{4,j}^{\left( P \right)} $ .4)隐含节点
$ Q $ 监听到来自活跃节点$ P $ 的数据包和时钟源节点$ O $ 的ACK时, 分别记录本地时间为$ t_{2,j}^{\left( Q \right)} $ 和$ t_{4,j}^{\left( Q \right)} $ .5)经过
$ N $ 轮的同步过程之后, 隐含节点$ Q $ 获得了一系列观测信息$ \left\{ t_{2,j}^{\left( Q \right)},t_{4,j}^{\left( Q \right)},{{R}_{j}} \right\}_{j = 1}^{N} $ . 至此, 隐含节点$ Q $ 可以利用该观测信息估计相对于时钟源节点$ O $ 的时钟漂移和偏移.在上述过程中, 活跃节点
$ P $ 与时钟源节点$ O $ 之间进行免时间戳同步, 隐含节点$ Q $ 监听它们之间的成对信息. 因此, 节点$ Q $ 只接收了信息, 节省了由发送报文产生的能耗.根据免时间戳同步的现有研究可知[16], 隐含节点
$ Q $ 若要估计相对于时钟源节点$ O $ 的时钟漂移与偏移, 除了记录接收数据包和ACK的本地时间戳之外, 还需要获取$ t_{2,j}^{\left( O \right)} $ (用于估计时钟偏移)和$ t_{3,j}^{\left( O \right)}- t_{2,j}^{\left( O \right)} $ (用于估计时钟漂移). 针对此, 本文通过预设时钟源节点返回ACK的时间$ t_{3,j}^{\left( O \right)} $ 同时获取$ t_{2,j}^{\left( O \right)} $ 和$ t_{3,j}^{\left( O \right)}-t_{2,j}^{\left( O \right)} $ . 变换步骤2)中的关系式可得:$$ \begin{equation} {R_j} = \left( {\xi - 1} \right)t_{2,j}^{\left( O \right)} - \left( {t_{3,j}^{\left( O \right)} - t_{2,j}^{\left( O \right)}} \right) \end{equation} $$ (1) 从式(1)中可以看出, 由于
$ \xi $ 大于1且接近于1, 因此$ t_{3,j}^{\left( O \right)}-t_{2,j}^{\left( O \right)} $ 的信息和部分$ t_{2,j}^{\left( O \right)} $ 的信息(即$ \left( \xi -1 \right)t_{2,j}^{\left( O \right)} $ )可以被嵌入到$ {{R}_{j}} $ 中.$ {{R}_{j}} $ 在本文中被定义为$ {{R}_{j}} = \left( \xi -1 \right)\left( j-1 \right)T $ , 其中$ T $ 是活跃节点$ P $ 与时钟源节点$ O $ 的固定交互周期. 隐含节点$ Q $ 在每轮监听过程中, 可以根据该定义式计算$ {{R}_{j}} $ , 从而隐含地获取到$ t_{3,j}^{\left( O \right)}-t_{2,j}^{\left( O \right)} $ 和部分$ t_{2,j}^{\left( O \right)} $ . 因此, 从获取同步信息的角度来看, 隐含节点$ Q $ 能够估计出两个同步参数(具体推导过程见第2节).$ {{R}_{j}} $ 的作用除了嵌入同步信息之外, 另一个作用是补偿时钟源节点$ O $ 的响应时间(因此称之为等待间隔补偿). 因为$ t_{2,j}^{\left( O \right)} $ 会随着交互次数的增多而不断增长, 如果没有$ {{R}_{j}} $ , 那么$ t_{3,j}^{\left( O \right)} = \xi t_{2,j}^{\left( O \right)} $ , 响应时间$ t_{3,j}^{\left( O \right)}- t_{2,j}^{\left( O \right)} = \left( \xi -1 \right)t_{2,j}^{\left( O \right)} $ 会逐渐增长, 最终会影响网络的正常运行. 此外, 设置$ \xi $ 为一个大于1且接近于1的常数有两个原因: 1)$ \xi $ 大于1可以保留部分$ t_{2,j}^{\left( O \right)} $ ; 2)$ \xi $ 接近于1可以确保响应时间合理, 因为$t_{3,j}^{\left( O \right)}- t_{2,j}^{\left( O \right)} = \left( \xi -1 \right)\left[ t_{2,j}^{\left( O \right)}-\left( j-1 \right)T \right]$ .1.2 时钟模型
网络中每个传感器节点的时钟相对于理想时间都存在时钟漂移和偏移. 其时钟模型可以表示为:
$$ \begin{equation} C\left( {{t_0}} \right) = \left( {1 + \alpha } \right){t_0} + \vartheta \end{equation} $$ (2) 其中,
$ C\left( {{t}_{0}} \right) $ 和$ {{t}_{0}} $ 分别表示节点的本地时间和理想时间,$ \alpha $ 和$ \vartheta $ 分别表示时钟漂移和偏移.同步信息在传输过程中会经历多种时延, 这些时延可分为固定时延和随机时延[17]. 其中固定时延包括传输时间、接收时间、传播时间, 随机时延包括发送时间、接收处理时间和信道访问时间. 一般情况下, 固定时延被假设为一个常数, 而随机时延则会根据不同的网络场景被建模为服从不同分布的随机变量. 例如, 延迟是由许多独立的随机过程相互叠加时, 根据中心极限定理, 随机时延可以被建模为高斯分布[18]. 此外, 该随机时延模型已被实验证明是合理的[19]. 因此, 在本文中, 假设固定时延是已知常数, 随机时延是独立同分布的高斯分布, 其均值为0, 方差为
$ {{\sigma }^{2}} $ .根据时钟模型式(2), 时钟源节点
$ O $ 在$ j $ 轮同步信息交换过程中, 本地时间戳$ t_{2,j}^{\left( O \right)} $ 可以表示为:$$ \begin{split} t_{2,j}^{\left( O \right)} = \;&\left( {1 + {\alpha ^{\left( {PO} \right)}}} \right)\times\\ &\left( {t_{1,j}^{\left( P \right)} + {d^{\left( {PO} \right)}} + \omega _j^{\left( {PO} \right)}} \right) + {\vartheta ^{\left( {PO} \right)}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \end{split} $$ (3) 其中,
$ {{\alpha }^{\left( PO \right)}} $ 和$ {{\vartheta }^{\left( PO \right)}} $ 分别表示节点$ P $ 相对节点$ O $ 的时钟漂移和偏移,$ {{d}^{\left( PO \right)}} $ 和$ \omega _{j}^{\left( PO \right)} $ 分别表示信息从节点$ P $ 传输到节点$ O $ 过程中的固定时延和随机时延.同样地, 时钟源节点
$ O $ 在$ j $ 轮交互过程中返回ACK的时间戳$ t_{3,j}^{\left( O \right)} $ 可以表示为:$$ \begin{split} t_{3,j}^{\left( O \right)} =\; & \left( {1 + {\alpha ^{\left( {QO} \right)}}} \right)\times\\ &\left( {t_{4,j}^{\left( Q \right)} - {d^{\left( {OQ} \right)}} - \omega _j^{\left( {OQ} \right)}} \right) + {\vartheta ^{\left( {QO} \right)}} \end{split} $$ (4) 其中,
$ {{\alpha }^{\left( QO \right)}} $ 和$ {{\vartheta }^{\left( QO \right)}} $ 分别表示节点$ Q $ 相对节点$ O $ 的时钟漂移和偏移,$ {{d}^{\left( OQ \right)}} $ 和$ \omega _{j}^{\left( OQ \right)} $ 分别表示信息从节点$ O $ 传输到节点$ Q $ 过程中的固定时延和随机时延.隐含节点
$ Q $ 在$ j $ 轮同步过程中监听到来自节点$ P $ 的数据包时的时间戳$ t_{2,j}^{\left( Q \right)} $ 可以表示为:$$ \begin{split} t_{2,j}^{\left( Q \right)} = \;& \left( {1 + {\alpha ^{\left( {PQ} \right)}}} \right)\times\\ &\left( {t_{1,j}^{\left( P \right)} + {d^{\left( {PQ} \right)}} + \omega _j^{\left( {PQ} \right)}} \right) + {\vartheta ^{\left( {PQ} \right)}} \end{split} $$ (5) 其中,
$ {{\alpha }^{\left( PQ \right)}} $ 和$ {{\vartheta }^{\left( PQ \right)}} $ 分别表示节点$ P $ 相对节点$ Q $ 的时钟漂移和偏移,$ {{d}^{\left( PQ \right)}} $ 和$ \omega _{j}^{\left( PQ \right)} $ 分别表示信息从节点$ P $ 传输到节点$ Q $ 过程中的固定时延和随机时延.2. 参数估计
基于上述的时间戳表达式, 本节将推导隐含节点时钟漂移和偏移的联合MLE以及相应的CRLB.
2.1 联合MLE
在实际无线传感网中,
$ {{\alpha }^{\left( PO \right)}} $ 、$ {{\alpha }^{\left( QO \right)}} $ 、$ {{\alpha }^{\left( PQ \right)}} $ 都非常小且接近于0, 导致累计的时钟偏移${{\alpha }^{\left( PO \right)}}\times \left( {{d}^{\left( PO \right)}} + \omega _{j}^{\left( PO \right)} \right)$ ,${{\alpha }^{\left( QO \right)}} \left( {{d}^{\left( OQ \right)}} + \omega _{j}^{\left( OQ\right)} \right)$ ,${{\alpha }^{\left( PQ \right)}} \Big( {{d}^{(PQ)}}+ \omega _{j}^{\left( PQ \right)} \Big)$ 相对较小. 因此, 为了简化推导, 忽略式(3)$ \sim $ 式(5)中的累计时钟偏移, 则第1节中的时间戳表达式可以分别重写为:$$ \begin{equation} t_{2,j}^{\left( O \right)} = {\alpha ^{\left( {PO} \right)}}t_{1,j}^{\left( P \right)} + t_{1,j}^{\left( P \right)} + {d^{\left( {PO} \right)}} + \omega _j^{\left( {PO} \right)} + {\vartheta ^{\left( {PO} \right)}} \end{equation} $$ (6) $$ \begin{equation} t_{3,j}^{\left( O \right)} = {\alpha ^{\left( {QO} \right)}}t_{4,j}^{\left( Q \right)} + t_{4,j}^{\left( Q \right)} - {d^{\left( {OQ} \right)}} - \omega _j^{\left( {OQ} \right)} + {\vartheta ^{\left( {QO} \right)}} \end{equation} $$ (7) $$ \begin{equation} t_{2,j}^{\left( Q \right)} = {\alpha ^{\left( {PQ} \right)}}t_{1,j}^{\left( P \right)} + t_{1,j}^{\left( P \right)} + {d^{\left( {PQ} \right)}} + \omega _j^{\left( {PQ} \right)} + {\vartheta ^{\left( {PQ} \right)}} \end{equation} $$ (8) 由于
$ t_{3,j}^{\left( O \right)} = \xi t_{2,j}^{\left( O \right)}-{{R}_{j}} $ , 将式(6)和式(7)代入可得:$$ \begin{split} {R_j} = \;&\xi t_{1,j}^{\left( P \right)}{\alpha ^{\left( {PO} \right)}} - t_{4,j}^{\left( Q \right)}{\alpha ^{\left( {QO} \right)}} + \xi t_{1,j}^{\left( P \right)}{\kern 1pt} -\\ &t_{4,j}^{\left( Q \right)} + {d^{\left( {OQ} \right)}} + \xi {d^{\left( {PO} \right)}} - {\vartheta ^{\left( {QO} \right)}}{\kern 1pt}+\\ &\xi {\vartheta ^{\left( {PO} \right)}} + \omega _j^{\left( {OQ} \right)} + \xi \omega _j^{\left( {PO} \right)} \end{split} $$ (9) 注意到, 式(9)中包含
$ {{\alpha }^{\left( PO \right)}} $ ,$ {{\alpha }^{\left( QO \right)}} $ ,$ {{\vartheta }^{\left( PO \right)}} $ ,$ {{\vartheta }^{\left( QO \right)}} $ 4个时钟参数, 而节点$ Q $ 要与节点$ O $ 达到时间同步, 只需估计$ {{\alpha }^{\left( QO \right)}} $ 和$ {{\vartheta }^{\left( QO \right)}} $ . 为了消除$ {{\alpha }^{\left( PO \right)}} $ 和$ {{\vartheta }^{\left( PO \right)}} $ 这2个无关参数, 将式(9)减去$ \xi $ 倍的式(8)可得:$$ \begin{split} {R_j} - \xi t_{2,j}^{\left( Q \right)} =\; &\xi t_{1,j}^{\left( P \right)}{\alpha ^{\left( {PO} \right)}} - \xi t_{1,j}^{\left( P \right)}{\alpha ^{\left( {PQ} \right)}} - t_{4,j}^{\left( Q \right)}{\alpha ^{\left( {QO} \right)}}-\\ &t_{4,j}^{\left( Q \right)} + {d^{\left( {OQ} \right)}} + \xi {d^{\left( {PO} \right)}} - \xi {d^{\left( {PQ} \right)}} -\\ &{\vartheta ^{\left( {QO} \right)}} + \xi {\vartheta ^{\left( {PO} \right)}} - \xi {\vartheta ^{\left( {PQ} \right)}} +\\ &\omega _j^{\left( {OQ} \right)} + \xi \omega _j^{\left( {PO} \right)} - \xi \omega _j^{\left( {PQ} \right)}\\[-12pt] \end{split} $$ (10) 根据3个节点之间时钟漂移和偏移的关系:
$ {{\alpha }^{\left( QO \right)}} = {{\alpha }^{\left( PO \right)}}-{{\alpha }^{\left( PQ \right)}} $ ,$ {{\vartheta }^{\left( QO \right)}} = {{\vartheta }^{\left( PO \right)}}-{{\vartheta }^{\left( PQ \right)}} $ , 可得系统的同步通式为:$$ \begin{split} {R_j} &- \xi t_{2,j}^{\left( Q \right)} + t_{4,j}^{\left( Q \right)} - {d^{\left( {OQ} \right)}} - \xi {d^{\left( {PO} \right)}} + \xi {d^{\left( {PQ} \right)}} = \\ &\left( {\xi t_{1,j}^{\left( P \right)} - t_{4,j}^{\left( Q \right)}} \right){\alpha ^{\left( {QO} \right)}} + \left( {\xi - 1} \right){\vartheta ^{\left( {QO} \right)}} +\\ &\;\,\omega _j^{\left( {OQ} \right)} + \xi \omega _j^{\left( {PO} \right)} - \xi \omega _j^{\left( {PQ} \right)} \\[-12pt]\end{split} $$ (11) 令
${{\Gamma }_{j}} ={{R}_{j}}\;-\;\xi t_{2,j}^{\left( Q \right)}\;+\;t_{4,j}^{\left( Q \right)}-{{d}^{\left( OQ \right)}}-\xi {{d}^{\left( PO \right)}}+ \xi {{d}^{\left( PQ \right)}}$ ,$ {{W}_{j}} = \omega _{j}^{\left( OQ \right)}+\xi \omega _{j}^{\left( PO \right)}-\xi \omega _{j}^{\left( PQ \right)} $ . 同步通式(11)可以简写为:$${\Gamma _j} = \left( {\xi t_{1,j}^{\left( P \right)} - t_{4,j}^{\left( Q \right)}} \right){\alpha ^{\left( {QO} \right)}} +\left( {\xi - 1} \right){\vartheta ^{\left( {QO} \right)}} + {W_j} $$ (12) 经观察发现, 式(12)还包含节点
$ P $ 的本地时间戳$ t_{1,j}^{\left( P \right)} $ , 节点$ Q $ 无法直接获取. 但是, 由于数据周期性发送,$ T = t_{1,j}^{\left( P \right)}-t_{1,j-1}^{\left( P \right)} $ , 隐含节点$ Q $ 可以根据$ t_{1,j}^{\left( P \right)} = \left( j-1 \right)T $ 计算得到$ t_{1,j}^{\left( P \right)} $ . 此外, 隐含节点$ Q $ 监听了$ N $ 轮的免时间戳交互, 可以获得$ N $ 个同步通式. 为了能够实现高精度的同步, 需利用$ N $ 个同步通式估计时钟漂移和偏移. 将所有的同步通式写成矩阵形式:$$ \begin{equation} \begin{split} \underbrace {\left[ {\begin{array}{*{20}{c}} {{\Gamma _1}}\\ \vdots \\ {{\Gamma _N}} \end{array}} \right]}_{ \buildrel \over {\boldsymbol{\Gamma}}} =\;& \underbrace {\left[ {\begin{array}{*{20}{c}} {\xi t_{1,1}^{\left( P \right)} - t_{4,1}^{\left( Q \right)}}&{\xi - 1}\\ \vdots & \vdots \\ {\xi t_{1,N}^{\left( P \right)} - t_{4,N}^{\left( Q \right)}}&{\xi - 1} \end{array}} \right]}_{ \buildrel \over {\boldsymbol{M}}}\times\\ &\underbrace {\left[ {\begin{array}{*{20}{c}} {{\alpha ^{\left( {QO} \right)}}}\\ {{\vartheta ^{\left( {QO} \right)}}} \end{array}} \right]}_{ \buildrel \over {\boldsymbol{\Theta }}} + \underbrace {\left[ {\begin{array}{*{20}{c}} {{W_1}}\\ \vdots \\ {{W_N}} \end{array}} \right]}_{ \buildrel \over {\boldsymbol{W}}} \end{split} \end{equation} $$ (13) 最大似然估计是一种能够简便地实现复杂估计问题求解的估计方法, 其是利用已知的观测数据, 反推最大概率产生该观测信息的参数值. 无线传感网中的时间同步就是通过观测的同步信息计算时钟参数. 此外, 在式(13)中,
$ \boldsymbol{\Gamma } $ 是由已知观测数据构成的矢量矩阵;$ \boldsymbol{M} $ 是秩为2的观测矩阵;$ \boldsymbol{\Theta } $ 是包含2个待估时钟参数的矢量, 并且与$ \boldsymbol{\Gamma } $ 呈线性关系;$ \boldsymbol{W} $ 是均值为0, 方差为$ \left( 1+2{{\xi }^{2}} \right){{\sigma }^{2}}\mathbf{1} $ 的高斯噪声矢量, 其中$ \mathbf{1} $ 是单位矩阵, 满足线性模型的最大似然估计定理. 因此, 根据文献[20]中的定理7.5, 时钟参数矢量$ \boldsymbol{\Theta } $ 的MLE为:$$ \begin{equation} {\boldsymbol{\Theta }} = {\left( {{{\boldsymbol{M}}^{\rm{H}}}{\boldsymbol{M}}} \right)^{ - 1}}{{\boldsymbol{M}}^{\rm{H}}}{\boldsymbol{\Gamma }} \end{equation} $$ (14) 将式(14)进行数学运算, 节点
$ Q $ 时钟漂移和偏移的最大似然估计器分别表示为:$$ \begin{equation} \begin{split} \left[ {\begin{array}{*{20}{c}} {{{\hat \alpha }^{\left( {QO} \right)}}}\\ {{{\hat \vartheta }^{\left( {QO} \right)}}} \end{array}} \right] =\;& \frac{1}{{\left( {\xi - 1} \right)\left( {N\displaystyle \sum\limits_{j = 1}^N {G_j^2 - {{\left( {\sum\limits_{j = 1}^N {{G_j}} } \right)}^2}} } \right)}}\times\\ &\left[ {\begin{array}{*{20}{c}} {\left( {\xi - 1} \right)\displaystyle \sum\limits_{j = 1}^N {{\Gamma _j}\left( {N{G_j} - \sum\limits_{j = 1}^N {{G_j}} } \right)} }\\ {\displaystyle \sum\limits_{j = 1}^N {{\Gamma _j}\left( {\displaystyle \sum\limits_{j = 1}^N {G_j^2} - {G_j}\displaystyle \sum\limits_{j = 1}^N {{G_j}} } \right)} } \end{array}} \right] \end{split} \end{equation} $$ (15) 其中,
$ {G_j} = \xi t_{1,j}^{\left( P \right)} - t_{4,j}^{\left( Q \right)} $ .所以, 利用估计器(15), 隐含节点
$ Q $ 可以估计自身的时钟漂移和偏移, 达到与时钟源节点$ O $ 的同步. 与现有的免时间戳同步与隐含同步相结合的同步机制相比, 本文所提的同步机制完成了对隐含节点完整时钟参数(时钟漂移和偏移)的估计, 实现了与时钟源节点的完全时间同步.2.2 克拉美罗下限(CRLB)
CRLB是衡量无偏估计器性能的一个理论标准[20], 无偏估计量的方差只能无限逼近或等于CRLB, 而不会小于CRLB, 当两者相等时, 表示估计器性能达到最优. 由于它容易被确定, 常被用来评估无偏估计器性能的好坏, 判断估计器性能是否达到了最优. 为了评估所提联合MLE的性能, 本节推导矢量参数
$\boldsymbol{\Theta } = {{\left[ {{\alpha }^{\left( QO \right)}} \;\; {{\vartheta }^{\left( QO \right)}} \ \right]}^{\rm{T}}}$ 的CRLB. 首先需要计算$ 2\times 2 $ 的费希尔信息矩阵, 然后再求其逆矩阵可分别得到$ {{\alpha }^{\left( QO \right)}} $ 和$ {{\vartheta }^{\left( QO \right)}} $ 的CRLB.由式(12)可得
$ \left( {{\alpha }^{\left( QO \right)}}, {{\vartheta }^{\left( QO \right)}} \right) $ 的对数似然函数如下:$$ \begin{split} &\ln L\left( {{\alpha ^{\left( {QO} \right)}},{\vartheta ^{\left( {QO} \right)}}} \right) = N\ln \frac{1}{\sigma{\sqrt {2\pi \left( {1 + 2{\xi ^2}} \right)} }}-\\ &\qquad\frac{{\displaystyle \sum\limits_{j = 1}^N {{{\left[ {{\Gamma _j} - {G_j}{\alpha ^{\left( {QO} \right)}} + \left( {\xi - 1} \right){\vartheta ^{\left( {QO} \right)}}} \right]}^2}} }}{{2\left( {1 + 2{\xi ^2}} \right){\sigma ^2}}} \end{split} $$ (16) 进一步, 似然函数分别对
$ {{\alpha }^{\left( QO \right)}} $ 和$ {{\vartheta }^{\left( QO \right)}} $ 求二阶导可得:$$ \begin{equation} \frac{{{\partial ^2}\ln L}}{{\partial {{\left( {{\alpha ^{\left( {QO} \right)}}} \right)}^2}}} = - \frac{1}{{\left( {1 + 2{\xi ^2}} \right){\sigma ^2}}}\sum\limits_{j = 1}^N {G_j^2} \end{equation} $$ (17) $$ \begin{equation} \frac{{{\partial ^2}\ln L}}{{\partial {{\left( {{\vartheta ^{\left( {QO} \right)}}} \right)}^2}}} = - \frac{{N{{\left( {\xi - 1} \right)}^2}}}{{\left( {1 + 2{\xi ^2}} \right){\sigma ^2}}} \end{equation} $$ (18) $$ \begin{equation} \frac{{{\partial ^2}\ln L}}{{\partial {\alpha ^{\left( {QO} \right)}}{\vartheta ^{\left( {QO} \right)}}}} = - \frac{{ {\xi - 1} }}{{\left( {1 + 2{\xi ^2}} \right){\sigma ^2}}}\sum\limits_{j = 1}^N {{G_j}} \end{equation} $$ (19) 再对式(17) ~ 式(19)求负期望可得费希尔信息矩阵如下:
$$ \begin{split} &{\bf{FIM}}\left( \lambda \right)=\\ & \left[ {\begin{array}{*{20}{c}} { - {\rm{E}} \left[ {\dfrac{{{\partial ^2}\ln L}}{{\partial {{\left( {{\alpha ^{\left( {QO} \right)}}} \right)}^2}}}} \right]}&{ - {\rm{E}} \left[ {\dfrac{{{\partial ^2}\ln L}}{{\partial {\alpha ^{\left( {QO} \right)}}{\vartheta ^{\left( {QO} \right)}}}}} \right]}\\ { - {\rm{E}} \left[ {\dfrac{{{\partial ^2}\ln L}}{{\partial {\alpha ^{\left( {QO} \right)}}{\vartheta ^{\left( {QO} \right)}}}}} \right]}&{ - {\rm{E}} \left[ {\dfrac{{{\partial ^2}\ln L}}{{\partial {{\left( {{\vartheta ^{\left( {QO} \right)}}} \right)}^2}}}} \right]} \end{array}} \right]=\;\;\;\;\\ & \dfrac{{\left[ {\begin{array}{*{20}{c}} {\displaystyle\sum\limits_{j = 1}^N {G_j^2} }&{\left( {\xi - 1} \right)\displaystyle\sum\limits_{j = 1}^N {{G_j}} }\\ {\left( {\xi - 1} \right)\displaystyle\sum\limits_{j = 1}^N {{G_j}} }&{N{{\left( {\xi - 1} \right)}^2}} \end{array}} \right]}}{{\left( {1 + 2{\xi ^2}} \right){\sigma ^2}}}\\[-25pt] \end{split} $$ (20) 最后, 求得费希尔信息矩阵的逆矩阵为:
$$ \begin{aligned} &{\bf{FI}}{{\bf{M}}^{ - 1}}\left( \lambda \right) = \\ &\frac{{\left( {1 + 2{\xi ^2}} \right){\sigma ^2}\left[ {\begin{array}{*{20}{c}} {N{{\left( {\xi - 1} \right)}^2}}&{ - \left( {\xi - 1} \right)\displaystyle\sum\limits_{j = 1}^N {{G_j}} }\\ { - \left( {\xi - 1} \right)\displaystyle\sum\limits_{j = 1}^N {{G_j}} }&{\displaystyle\sum\limits_{j = 1}^N {G_j^2} } \end{array}} \right]}}{{{{\left( {\xi - 1} \right)}^2}\left[ {N\displaystyle\sum\limits_{j = 1}^N {G_j^2} - {{\left( {\displaystyle\sum\limits_{j = 1}^N {{G_j}} } \right)}^2}} \right]}} \end{aligned} $$ (21) 费希尔信息矩阵的逆矩阵的对角元素即为相应参数的CRLB[20]. 因此, 时钟漂移
$ {{\alpha }^{\left( QO \right)}} $ 和时钟偏移$ {{\vartheta }^{\left( QO \right)}} $ 的CRLB分别表示为:$$ \begin{equation} {\mathop{\rm var}} \left( {{{\hat \alpha }^{\left( {QO} \right)}}} \right) \ge \frac{{N\left( {1 + 2{\xi ^2}} \right){\sigma ^2}}}{{N\displaystyle\sum\limits_{j = 1}^N {G_j^2} - {{\left( {\displaystyle\sum\limits_{j = 1}^N {{G_j}} } \right)}^2}}} \end{equation} $$ (22) $$ {\mathop{\rm var}} \left( {{{\hat \vartheta }^{\left( {QO} \right)}}} \right) \ge \\ \frac{{\left( {1 + 2{\xi ^2}} \right){\sigma ^2}\displaystyle\sum\limits_{j = 1}^N {G_j^2} }}{{{{\left( {\xi - 1} \right)}^2}\left[ {N\displaystyle\sum\limits_{j = 1}^N {G_j^2} - {{\left( {\displaystyle\sum\limits_{j = 1}^N {{G_j}} } \right)}^2}} \right]}} $$ (23) 3. 仿真验证与对比分析
本节基于MATLAB仿真平台, 验证所提隐含节点的联合最大似然估计器的有效性. 同时, 在估计性能、能量开销、计算数量三个方面, 将所提同步机制与现有的相关同步机制进行了对比.
3.1 仿真验证
在无线传感网中, 固定时延和同步周期通常是ms量级, 而对于典型的石英晶体振荡器, 时钟漂移的数量级通常为
$ {{10}^{-6}}\; $ s. 从而可以得出, 时钟漂移大约是时钟偏移和固定时延的$ {{10}^{-3}} $ 倍. 因此, 在仿真验证中, 将活跃节点$ P $ 相对于时钟源节点$ O $ 、活跃节点$ P $ 相对于隐含节点$ Q $ 的时钟漂移和偏移分别设置为:$ {{\alpha }^{( PO )}} = 0.003 $ ,$ {{\alpha }^{( PQ )}} = 0.001 $ ,$ {{\vartheta }^{( PO )}} = [ -5,5 ] $ ,$ {{\vartheta }^{( PQ )}} = [ -2.5,2.5 ] $ . 隐含节点$ Q $ 相对于时钟源节点$ O $ 的时钟漂移$ {{\alpha }^{( QO )}} $ 和偏移$ {{\vartheta }^{( QO )}} $ 可由它们之间的关系得出. 此外, 节点$ P $ 发送数据包的周期$ T $ 设置为80 ms, 高斯随机时延的标准差设置为$ \sigma = 0.2 $ , 信息传输链路中的固定时延取值范围分别设置为:$ {{d}^{( PO )}} = [ 3,13 ] $ ,${{d}^{( OP )}} = [ 0,10 ]$ ,${{d}^{( PQ )}} = [ 3,13 ]$ ,${{d}^{( OQ )}} = [ 0,10 ]$ , 节点$ O $ 在每轮同步过程中计算ACK响应时间$ t_{3,j}^{( O )} $ 的系数$ \xi = 1.4 $ . 在估计理论中, 通常选择均方误差(Mean square error, MSE)作为估计性能标准, 因此这里采用MSE来衡量所提时钟漂移$ {{\alpha }^{( QO )}} $ 和偏移$ {{\vartheta }^{( QO )}} $ 估计器的性能. 为了使仿真结果更加可靠, 在所有的仿真图中, 每个样本点的值都是10000次仿真运行后的平均值.图3为隐含节点
$ Q $ 相对于时钟源节点$ O $ 的时钟漂移最大似然估计器$ {{\hat{\alpha }}^{\left( QO \right)}} $ 的MSE和CRLB. 从图中可以看出, 时钟漂移最大似然估计器$ {{\hat{\alpha }}^{\left( QO \right)}} $ 的MSE曲线与CRLB曲线基本重合, 表明了其估计性能达到最优, 同时也验证了估计器的有效性. 而且随着观测次数的增加, 所提估计器的估计精度不断提高.
图 3 隐含节点$Q$ 时钟漂移估计${{\hat{\alpha }}^{\left(QO \right)}}$ 的MSE与CRLBFig. 3 MSE and CRLB of estimated clock skew${{\hat{\alpha }}^{\left(QO \right)}}$ for silent node$Q$ 图4所示为隐含节点
$ Q $ 相对于时钟源节点$ O $ 的时钟偏移最大似然估计器$ {{\hat{\vartheta }}^{\left( QO \right)}} $ 的仿真结果. 在图4中, 时钟偏移最大似然估计器$ {{\hat{\vartheta }}^{\left( QO \right)}} $ 的性能随着观测次数的增加不断提高, 而且MSE曲线与CRLB曲线基本重合, 表明该估计器的性能达到了最优, 同时验证了该估计器的有效性. 值得注意的是, 在仿真中, 系数$ \xi = 1.4 $ , 这表明隐含节点$ Q $ 在通信过程中只是获取到2/5的$ t_{2,j}^{\left( O \right)} $ , 就以此实现了相对于时钟源节点$ O $ 的时钟偏移估计. 此外, 时钟偏移的估计精度是影响同步精度的主要因素. 因为隐含节点$ Q $ 根据估计的时钟漂移与偏移调整本地时间, 减小与时钟源节点$ O $ 的时间偏差, 而时钟漂移的估计精度比时钟偏移的估计精度高出约$ {{10}^{6}} $ 个数量级, 其对同步精度的影响较小, 可以忽略不计. 所以可以利用时钟偏移的估计精度反映节点间的同步精度. 从图4中可以看出, 时钟偏移的估计精度为$ {{10}^{-1}} $ 数量级, 而时钟偏移被设置为数ms, 因此, 所提同步算法可达到数百μs的同步精度.
图 4 隐含节点$Q$ 时钟偏移估计${{\hat{\vartheta }}^{\left(QO \right)}}$ 的MSE与CRLBFig. 4 MSE and CRLB of estimated clock offset${{\hat{\vartheta }}^{\left(QO \right)}}$ for silent node$Q$ 从仿真结果可以看出, 隐含节点
$ Q $ 监听活跃节点$ P $ 和时钟源节点$ O $ 之间的免时间戳同步过程, 通过观测到的数据集$ \left\{ t_{2,j}^{\left( Q \right)},t_{4,j}^{\left( Q \right)}, {{R}_{j}} \right\}_{j = 1}^{N} $ 就可以同时估计相对于节点$ O $ 的时钟漂移和偏移, 达到与节点$ O $ 的同步. 与文献[15]中的基于动态响应的免时间戳的隐含同步方式相比, 所提的免时间戳隐含同步方法继承了现有免时间戳同步和隐含同步的优势, 能够大幅度地减少同步能量消耗, 同时也可以无缝地嵌入到网络数据流之中. 进一步, 隐含节点又能够联合估计相对于时钟源节点的时钟漂移和偏移, 实现隐含节点的完全同步, 弥补了目前两种机制结合使用的不足.3.2 对比分析
首先, 将所提的隐含节点
$ Q $ 的时钟漂移最大似然估计器$ {{\hat{\alpha }}^{\left( QO \right)}} $ 与类似的估计器(文献[15]中的式(26))进行估计性能比较. 为了确保比较的公平性, 在两种估计器的对比仿真中, 各个参数的初始化设置均相同, 文献[15]中两个连续周期内ACK的响应时间间隔设置为$ \left\{ 10,30 \right\} $ . 仿真结果如图3所示, 本文所提的隐含节点$ Q $ 的时钟漂移估计器的性能显著优于文献[15]中的时钟漂移估计器, 而且其估计精度大约高出101 ~ 103个数量级. 其原因是所提同步机制通过预设时钟源节点$ O $ 返回ACK的本地时间, 传递了部分时间戳信息$ \left( \xi -1 \right)t_{2,j}^{\left( O \right)} $ 给隐含节点$ Q $ , 时间戳信息的获取有利于时钟漂移估计精度的提高. 图5所示为不同系数$ \xi $ 下时钟漂移估计器性能对比结果, 其中设置$\xi = \{ 1.2,1.3,1.4, 1.5,1.6 \}$ . 从图中可以看出, 随着系数$ \xi $ 的增大, 所提时钟漂移估计器$ {{\hat{\alpha }}^{\left( QO \right)}} $ 的估计性能不断提高, 表明隐含节点$ Q $ 获取的时间戳信息越多, 时钟漂移的估计效果越好. 同时, 从仿真的角度说明, 相较于文献[15]中无法获取时间戳信息的隐含节点, 所提同步机制的隐含节点由于能够间接地获取部分时间戳信息, 因此其时钟漂移估计器性能更优. 但是, 需要注意的是, 不能为了提高时钟漂移的估计性能无限地增加$ \xi $ , 因为$ \xi $ 的取值必须确保响应时间间隔不会过大, 影响网络的正常运行.
图 5 不同系数$\xi $ 下时钟漂移估计器性能对比结果Fig. 5 The performance comparison results of clock skew estimator under different coefficient$\xi $ 其次, 将所提同步机制与隐含同步机制[7]、基于动态响应的免时间戳同步机制[15]进行信令流程及能耗开销对比. 无线传感网的时间同步依赖于信息交互, 交互流程所导致的能耗开销是影响同步能耗的重要因素. 在隐含同步机制中, 隐含节点监听的是传统的双向信息交互过程, 需要专用的同步帧传递时间戳信息, 同步帧的传输需要消耗额外的能量. 在基于动态响应的免时间戳同步机制中, 待同步节点无需发送时间戳信息, 在普通数据收发中实现同步, 无需消耗额外的能量. 在所提机制中, 隐含节点监听的是时钟源节点和活跃节点的免时间戳交互过程, 不需要专用同步帧, 没有额外的能量开销. 因此, 相较于隐含同步机制, 所提机制避免了网络中专用同步帧传输, 显著减少了能耗. 与免时间戳机制相比, 所提机制能够一次性实现公共广播领域内多个节点的同步, 而免时间戳机制只能实现一个节点的同步, 对于网络整体的能耗而言, 所提机制更具优势.
最后, 分析所提同步算法、隐含同步算法[7]、基于动态响应的免时间戳同步算法[15]以及免时间戳与隐含同步结合算法[15]的计算复杂度, 主要是对算法中加法和减法、乘法和除法计算数量进行比较, 对比结果如表1所示. 从表中可以看出, 4种同步算法时钟参数的计算复杂度均为
$ {\rm{O}}({{N}^{2}}) $ . 而本文算法时钟参数乘除法的计算数量高于隐含同步算法和免时间戳同步算法时钟参数计算数量, 这是因为时钟源节点$ O $ 返回ACK的本地时间中包含系数$ \xi $ , 导致算法中乘除法的数量相对较高. Wang等[16]研究表明, 传输1比特数据超过100米消耗的能量与执行300万条指令所需能量大致相等. 因此, 虽然本文算法的计算数量略高, 但能够在不传输专用同步帧的情况下实现多个节点的同步, 它仍然是一个高能效的同步方案. 相较于文献[15]中的免时间戳与隐含同步结合算法, 本文算法时钟漂移的计算复杂度略低于其时钟漂移计算复杂度, 而且本文算法还可以估计时钟偏移, 实现完整的同步.表 1 本文算法与隐含同步算法、免时间戳同步算法以及免时间戳和隐含同步结合算法的计算数量对比结果Table 1 The comparison results of the number of calculations among proposed algorithm, implicit synchronization algorithm, timestamp-free synchronization algorithm and combination algorithm of timestamp-free and implicit synchronization算法 参数 加减法数量 乘除法数量 本文算法 时钟漂移 $2{N^2} + 14N - 3$ ${N^2} + 11N + 3$ 时钟偏移 $5{N^2} + 13N - 2$ $4{N^2} + 11N + 4$ 隐含同步算法 时钟漂移 $2{N^2} + 12N - 3$ $3N + 3$ 时钟偏移 $5{N^2} + 11N - 3$ ${N^2} + 3N + 3$ 免时间戳同步算法 时钟漂移 ${N^2} + 6N - 2$ $4N + 3$ 免时间戳与隐含同步结合算法 时钟漂移 $4{N^2} + 10N - 1$ $4{N^2} + 12N + 1$ 4. 结论
本文提出了一种能够联合估计免时间戳同步场景中隐含节点的时钟漂移和偏移的同步协议. 在高斯随机时延下, 推导了隐含节点时钟漂移和偏移的最大似然估计器和相应的CRLB. 最后, 通过数值仿真验证了所提估计器的有效性, 并且与现有的相关同步方案在估计性能、能量开销和计算数量三个方面进行了对比, 论证了所提方案低功耗的优势.
-
图 6 5、10、15、25目标DTLZ4问题上的IGD进化轨迹
Fig. 6 Evolutionary trajectories of IGD on 5、10、15、25-objective DTLZ4
表 1 测试问题
Table 1 Test questions
问题 目标数目M 决策变量D PF DTLZ1 5, 10, 15, 25 M−1+5 线性 DTLZ2~DTLZ4 5, 10, 15, 25 M−1+10 凹型 DTLZ5~DTLZ6 5, 10, 15, 25 M−1+10 退化 DTLZ7 5, 10, 15, 25 M−1+20 断开 IDTLZ1 5, 10, 15, 25 M−1+5 倒置 IDTLZ2 5, 10, 15, 25 M−1+10 倒置 WFG1 5, 10, 15, 25 M−1+10 混合 WFG2 5, 10, 15, 25 M−1+10 断开 WFG3 5, 10, 15, 25 M−1+10 退化 WFG4~DTLZ9 5, 10, 15, 25 M−1+10 凹型 
下载: 导出CSV
表 2 不同目标对应的种群规模
Table 2 Population sizes corresponding to different objectives
目标数目M (H1, H2) 种群规模N 5 (6, 0) 126 10 (3, 2) 275 15 (2, 1) 135 25 (2, 1) 350
下载: 导出CSV
表 3 不同目标对应的终止条件
Table 3 Termination conditions corresponding to different objectives
目标数目M 种群规模N 进化代数 评价次数 5 126 800 100800 10 275 1000 275000 15 135 1300 175500 25 350 1500 525000
下载: 导出CSV
表 4 MOEA-DD、NSGA-III、RVEA、MOMBI-II、AR-MOEA以及MaOEA-IAR在DTLZ1~DTLZ4, WFG4~WFG9上获得IGD值的统一结果(均值和标准差). 最好的结果已被标记
Table 4 The statistical results (mean and standard deviation) of the IGD values obtained by MOEA-DD、NSGA-III、RVEA、MOMBI-II、AR-MOEA and MaOEA-IAR on DTLZ1~DTLZ4, WFG4~WFG9. The best results are highlighted
问题 M MOEA-DD NSGA-III RVEA MOMBI-II AR-MOEA MaOEA-IAR DTLZ1 5 6.3354×10−2
(8.24×10−5)+6.3588×10−2
(2.93×10−4)+6.3365×10−2
(1.44×10−4)+6.7530×10−2
(8.96×10−3)+1.3555×10−3
(7.87×10−6)+2.6870×10−1
(7.95×10−2)10 1.0784×10−1
(3.02×10−4)+1.2073×10−1
(2.91×10−2)+1.0804×10−1
(1.45×10−3)+2.1744×10−1
(2.32×10−2)+1.3680×10−3
(1.14×10−4)+4.7099×10+1
(1.33×10+1)15 1.4088×10−1
(5.94×10−3)+1.9693×10−1
(2.97×10−2)+1.6130×10−1
(7.54×10−3)+2.8789×10−1
(3.92×10−2)+3.1383×10−3
(1.58×10−3)+5.6383×10+1
(2.03×10+1)25 1.0354×10−1
(5.86×10−2)+2.7045×10−1
(1.96×10−1)=1.3372×10−1
(7.61×10−2)+2.0490×10−1
(1.19×10−1)+1.2392×10−1
(7.02×10−2)+2.4067×10−1
(1.44×10−1)DTLZ2 5 1.9506×10−1
(8.00×10−5)+1.9538×10−1
(1.18×10−4)+1.9508×10−1
(7.34×10−5)+2.0140×10−1
(1.93×10−3)=4.3609×10−3
(7.50×10−6)+2.0137×10−1
(1.61×10−3)10 4.2179×10−1
(8.12×10−4)−4.4544×10−1
(3.60×10−2)−4.2291×10−1
(1.61×10−3)−4.4619×10−1
(6.50×10−3)−4.9235×10−2
(1.43×10−4)−4.7588×10−2
(4.60×10−4)15 6.1801×10−1
(1.88×10−2)−6.4108×10−1
(2.01×10−2)−6.3958×10−1
(3.67×10−2)−8.5287×10−1
(6.41×10−2)−1.0058×10−2
(1.56×10−3)=1.1135×10−2
(1.23×10−3)25 4.3015×10−1
(2.45×10−1)−5.3577×10−1
(3.10×10−1)−4.3033×10−1
(2.45×10−1)−7.0528×10−1
(4.12×10−1)−6.1891×10−1
(3.45×10−1)−4.2923×10−1
(2.44×10−1)DTLZ3 5 1.9907×10−1
(3.45×10−3)+7.1945×10−1
(1.21×10+0)+2.0829×10−1
(2.11×10−2)+2.0571×10−1
(4.19×10−3)+4.4165×10−3
(5.43×10−5)+1.0190×10+0
(6.60×10−1)10 4.2392×10−1
(3.09×10−3)+1.1841×10+0
(1.08×10+0)+4.3390×10−1
(1.54×10−2)+6.9344×10−1
(1.66×10−1)+3.1091×10−2
(9.02×10−2)+4.8634×10+2
(9.91×10+1)15 6.6842×10−1
(2.26×10−1)+5.5656×10+0
(2.83×10+0)+8.3685×10−1
(3.72×10−1)+1.1033×10+0
(2.16×10−2)+3.8594×10−1
(5.94×10−1)+6.4283×10+2
(1.08×10+2)25 4.3041×10−1
(2.46×10−1)+8.3189×10+0
(6.16×10+0)−4.3076×10−1
(2.46×10−1)+8.7119×10−1
(5.08×10−1)+4.6567×10−1
(2.38×10−1)+9.8829×10−1
(5.81×10−1)DTLZ4 5 1.9508×10−1
(9.98×10−5)+2.4346×10−1
(9.84×10−2)−1.9514×10−1
(1.62×10−4)+2.5703×10−1
(9.71×10−2)−4.4014×10−3
(1.06×10−4)+2.0129×10−1
(1.72×10−3)10 4.2271×10−1
(1.73×10−3)−4.5035×10−1
(3.18×10−2)−4.3317×10−1
(2.13×10−3)−4.5814×10−1
(1.71×10−2)−5.7202×10−3
(5.61×10−4)−5.6898×10−3
(1.01×10−4)15 6.3495×10−1
(1.18×10−2)−6.4865×10−1
(2.13×10−2)−6.3283×10−1
(7.39×10−3)−6.5937×10−1
(2.31×10−2)−1.1470×10−2
(1.34×10−3)−9.9175×10−3
(1.41×10−3)25 4.3041×10−1
(2.45×10−1)−4.6255×10−1
(2.72×10−1)−4.3053×10−1
(2.46×10−1)−4.3098×10−1
(2.46×10−1)−5.8948×10−1
(2.86×10−1)−4.3021×10−1
(2.44×10−1)WFG4 5 1.2337×10+0
(4.54×10−3)−1.1663×10+0
(3.61×10−3)−1.1649×10+0
(2.28×10−3)−1.8356×10+0
(2.02×10−1)−1.1244×10+0
(1.47×10−2)−1.1192×10+0
(2.58×10−1)10 6.0999×10+0
(1.66×10−1)−4.5156×10+0
(4.39×10−2)−4.3816×10+0
(5.05×10−2)−5.6693×10+0
(4.27×10−1)−4.5463×10+0
(9.53×10−3)−4.0153×10+0
(2.71×10−2)15 9.4478×10+0
(7.47×10−1)−9.1431×10+0
(9.63×10−2)−9.4338×10+0
(3.83×10−1)−2.0211×10+1
(1.35×10+0)−9.3716×10+0
(3.68×10−2)−8.4000×10+0
(1.11×10−1)25 2.4032×10+1
(2.95×10−1)−1.3110×10+1
(7.97×10−1)−1.2394×10+1
(2.92×10−1)−3.7271×10+1
(3.18×10+0)−1.5196×10+1
(2.61×10−1)−1.1405×10+1
(1.89×10−2)WFG5 5 1.2116×10+0
(6.46×10−3)−1.1459×10+0
(4.71×10−3)−1.1554×10+0
(2.17×10−3)−1.5867×10+0
(1.16×10−1)−1.1648×10+0
(2.17×10−4)−1.1148×10+0
(1.31×10−2)10 6.2798×10+0
(1.39×10−1)−4.4654×10+0
(1.82×10−2)−4.3730×10+0
(4.27×10−2)−5.3610×10+0
(2.19×10−2)−4.5296×10+0
(9.52×10−3)−3.9834×10+0
(1.90×10−2)15 1.1246×10+1
(1.85×10−1)−8.9556×10+0
(7.70×10−2)−9.8645×10+0
(1.66×10−1)−2.4601×10+1
(1.34×10+0)−9.3564×10+0
(5.48×10−2)−8.2316×10+0
(8.70×10−2)25 2.2399×10+1
(1.96×10−1)−1.1842×10+1
(1.14×10+0)=1.1407×10+1
(5.45×10−2)=4.5754×10+1
(3.19×10+0)−1.4772×10+1
(2.79×10−1)−1.1370×10+1
(3.12×10−3)WFG6 5 1.2276×10+0
(1.14×10−2)−1.1617×10+0
(2.62×10−3)−1.1638×10+0
(2.18×10−3)−2.1021×10+0
(3.69×10−1)−1.1624×10+0
(1.24×10−3)−1.1443×10+0
(2.89×10−2)10 6.0366×10+0
(1.65×10−1)−4.5795×10+0
(1.88×10−2)−4.4020×10+0
(7.33×10−2)−5.2873×10+0
(2.86×10−2)−4.5333×10+0
(1.12×10−2)−4.0605×10+0
(4.78×10−2)15 1.0818×10+1
(1.01×10+0)−9.3642×10+0
(3.84×10−1)−1.0520×10+1
(4.62×10−1)−1.8732×10+1
(1.89×10+0)−9.3953×10+0
(5.07×10−2)−8.5397×10+0
(1.76×10−1)25 2.1950×10+1
(6.43×10−1)−1.5582×10+1
(9.19×10−1)−1.5281×10+1
(9.27×10−1)−3.5035×10+1
(2.81×10+0)−1.5206×10+1
(3.56×10−1)−1.1375×10+1
(1.29×10−2)WFG7 5 1.2454×10+0
(8.15×10−3)−1.1680×10+0
(3.67×10−3)−1.1693×10+0
(3.05×10−3)−1.8487×10+0
(2.46×10−1)−1.1787×10+0
(8.19×10−4)−1.1499×10+0
(3.17×10−2)10 5.1499×10+0
(2.73×10−1)−4.5094×10+0
(3.76×10−2)−4.3149×10+0
(5.44×10−2)−5.4063×10+0
(5.76×10−2)−4.5123×10+0
(1.07×10−2)−3.9423×10+0
(4.82×10−2)15 8.9281×10+0
(8.74×10−2)−9.0694×10+0
(9.98×10−2)−9.1250×10+0
(2.71×10−1)−1.7892×10+1
(1.65×10+0)−9.4000×10+0
(3.81×10−2)−8.7103×10+0
(1.27×10−1)25 1.6767×10+1
(1.55×10+0)−1.3522×10+1
(7.95×10−1)−1.1672×10+1
(3.01×10−1)=3.2655×10+1
(2.66×10+0)−1.5873×10+1
(1.94×10−1)−1.1572×10+1
(8.12×10−2)WFG8 5 1.2301×10+0
(8.57×10−3)−1.1867×10+0
(9.32×10−2)=1.1770×10+0
(4.82×10−3)−2.9462×10+0
(3.23×10−2)−1.1575×10+0
(2.95×10−3)−1.1675×10+0
(1.66×10−2)10 5.2132×10+0
(3.48×10−1)−4.5532×10+0
(3.00×10−1)−4.3661×10+0
(8.88×10−2)−6.0943×10+0
(8.02×10−1)−4.6273×10+0
(4.43×10−2)−4.1413×10+0
(7.20×10−2)15 9.2706×10+0
(3.66×10−1)−9.1735×10+0
(4.45×10−1)−9.9548×10+0
(7.95×10−1)−2.0437×10+1
(1.84×10+0)−9.3736×10+0
(7.14×10−2)−8.7680×10+0
(2.43×10−1)25 2.2579×10+1
(3.17×10+0)−1.6519×10+1
(6.99×10−1)−1.3760×10+1
(1.73×10+0)−3.9378×10+1
(2.20×10+0)−1.5213×10+1
(2.84×10−1)−1.1409×10+1
(2.03×10−2)WFG9 5 1.2087×10+0
(1.17×10−2)−1.1112×10+0
(1.08×10−2)−1.1315×10+0
(4.39×10−3)−2.6613×10+0
(1.67×10−1)−1.1476×10+0
(2.59×10−3)−1.0867×10+0
(1.42×10−2)10 5.2551×10+0
(3.81×10−1)−4.3117×10+0
(7.01×10−2)−4.3306×10+0
(6.41×10−2)−5.2773×10+0
(6.59×10−2)−4.4943×10+0
(1.01×10−2)−3.8995×10+0
(3.21×10−2)15 8.9598×10+0
(5.37×10−1)−8.6320×10+0
(1.41×10−1)−8.8404×10+0
(2.73×10−1)−2.5164×10+1
(1.97×10+0)−9.1233×10+0
(3.81×10−2)−8.1162×10+0
(1.37×10−1)25 1.5766×10+1
(3.46×10+0)−1.2211×10+1
(7.44×10−1)−1.2343×10+1
(1.28×10+0)−4.8512×10+1
(2.40×10+0)−1.5086×10+1
(2.89×10−1)−1.1166×10+1
(8.75×10−2)+/−/= 10/30/0 7/30/3 10/28/2 8/31/1 10/29/1 “+”表明该算法优于MaOEA-IAR, “−”劣于MaOEA-IAR, “=”则表示与MaOEA-IAR性能相似
下载: 导出CSV
表 5 MOEA-DD、NSGA-III、RVEA、MOMBI-II、AR-MOEA以及MaOEA-IAR在DTLZ1~DTLZ4, WFG4~WFG9上获得PD值的统一结果(均值和标准差). 最好的结果已被标记
Table 5 The statistical results (mean and standard deviation) of the PD values obtained by MOEA-DD、NSGA-III、RVEA、MOMBI-II、AR-MOEA and MaOEA-IAR on DTLZ1~DTLZ4, WFG4~WFG9. The best results are highlighted
问题 M MOEA-DD NSGA-III RVEA MOMBI-II AR-MOEA MaOEA-IAR DTLZ1 5 7.0347×10+5
(1.50×10+5)−6.6822×10+5
(2.89×10+5)−1.0144×10+6
(2.23×10+5)−6.5737×10+5
(5.56×10+5)−3.1124×10+5
(1.14×10+5)−1.1335×10+7
(2.48×10+7)10 3.9493×10+9
(4.51×10+8)−4.3920×10+9
(1.31×10+9)−4.4832×10+9
(5.54×10+8)−1.3231×10+8
(1.60×10+8)−3.9031×10+9
(3.47×10+8)−2.1838×10+12
(4.55×10+11)15 2.0771×10+11
(2.15×10+10)−7.2781×10+10
(4.17×10+10)−1.8665×10+11
(5.48×10+10)−2.8086×10+8
(4.77×10+8)−5.4706×10+10
(1.65×10+10)−2.4099×10+13
(8.65×10+12)25 4.2240×10+13
(1.16×10+12)+2.4959×10+13
(1.82×10+13)+1.5947×10+13
(3.82×10+12)+2.6381×10+9
(4.68×10+9)−1.2365×10+13
(6.44×10+12)+1.6674×10+11
(4.44×10+11)DTLZ2 5 4.5169×10+6
(3.69×10+5)−4.1624×10+6
(5.12×10+5)−4.4519×10+6
(4.97×10+5)−3.1301×10+6
(4.36×10+5)−2.2535×10+6
(1.23×10+5)−1.3815×10+7
(1.83×10+6)10 1.0700×10+10
(4.94×10+8)−1.5036×10+10
(4.96×10+9)−1.4787×10+10
(8.21×10+8)−6.3297×10+9
(9.76×10+8)−1.4371×10+10
(9.34×10+8)−4.8565×10+10
(8.50×10+9)15 2.3742×10+11
(5.12×10+10)−5.4986×10+10
(3.53×10+10)−1.5788×10+11
(3.98×10+10)−1.4169×10+9
(1.60×10+9)−2.2431×10+10
(1.21×10+10)−4.6071×10+11
(2.42×10+11)25 1.5327×10+13
(1.82×10+12)−1.9112×10+13
(9.88×10+12)−1.8880×10+13
(3.09×10+12)−1.4360×10+10
(2.20×10+10)−2.1209×10+12
(1.31×10+12)−1.0367×10+15
(3.68×10+13)DTLZ3 5 1.0522×10+7
(3.37×10+6 )+1.8368×10+7
(1.98×10+7 )+9.1974×10+6
(3.08×10+6 )+2.8101×10+6
(4.82×10+5 )+3.0005×10+6
(3.36×10+5 )+2.0182×10+6
(4.79×10+6)10 1.3601×10+10
(2.86×10+9)−2.2578×10+10
(1.68×10+10)−1.6599×10+10
(2.67×10+9)−1.1037×10+9
(1.66×10+9)−4.7006×10+9
(1.96×10+9)−2.9582×10+13
(4.66×10+12)15 3.6156×10+11
(8.59×10+11)−1.1462×10+12
(8.77×10+11)−1.2788×10+11
(2.39×10+11)−3.1557×10+8
(4.47×10+8)−4.9753×10+10
(4.57×10+10)−6.5391×10+14
(1.63×10+14)25 1.6252×10+13
(1.40×10+12)+3.2406×10+15
(2.26×10+15)+1.0735×10+13
(5.81×10+12)+7.3339×10+8
(1.67×10+9)−2.2266×10+13
(2.19×10+13)+9.6558×10+10
(2.77×10+11)DTLZ4 5 3.7620×10+6
(6.71×10+5)−3.3390×10+6
(8.43×10+5)−4.4024×10+6
(6.07×10+5)−2.1806×10+6
(1.14×10+6)−2.4008×10+6
(1.56×10+5)−1.3225×10+7
(1.43×10+6)10 1.3637×10+10
(1.08×10+9)+5.6960×10+9
(2.42×10+9)−1.0848×10+10
(2.29×10+9)+5.9777×10+7
(2.68×10+7)−4.0114×10+9
(1.10×10+9)−9.1937×10+9
(2.85×10+9)15 2.7750×10+9
(1.79×10+9)−5.2847×10+8
(5.23×10+8)−2.0648×10+9
(1.07×10+9)−4.0923×10+7
(1.33×10+8)−5.8144×10+9
(3.64×10+9)−3.9392×10+10
(1.98×10+10)25 3.7432×10+12
(9.47×10+11)−3.1839×10+11
(3.51×10+11)−3.8552×10+11
(1.51×10+11)−2.8053×10+8
(5.19×10+8)−1.2420×10+11
(5.99×10+10)−8.4925×10+14
(3.45×10+13)WFG4 5 1.2326×10+8
(6.16×10+6)−1.3502×10+8
(4.43×10+6)−1.2748×10+8
(7.49×10+6)−1.8889×10+7
(3.43×10+6)−6.4422×10+7
(4.60×10+6)−2.0686×10+8
(1.12×10+7)10 1.4338×10+11
(9.29×10+9)−2.9001×10+11
(1.68×10+10)−1.9392×10+11
(1.22×10+10)−7.4282×10+10
(2.70×10+10)−1.4488×10+11
(8.57×10+9)−6.6989×10+11
(2.52×10+10)15 3.4224×10+12
(5.49×10+11)−5.5478×10+12
(7.84×10+11)−4.0484×10+12
(1.04×10+12)−2.0839×10+10
(1.96×10+10)−2.0698×10+12
(6.60×10+11)−1.9906×10+13
(1.68×10+12)25 9.4061×10+13
(1.83×10+13)−2.2238×10+15
(6.40×10+14)−1.6309×10+14
(2.64×10+13)−2.5978×10+11
(6.99×10+11)−3.9968×10+14
(1.07×10+14)−1.0349×10+16
(5.36×10+14)WFG5 5 1.2123×10+8
(6.36×10+6)−1.4121×10+8
(8.96×10+6)−1.1886×10+8
(4.50×10+6)−2.3887×10+7
(4.70×10+6)−4.1612×10+7
(3.93×10+6)−1.9277×10+8
(1.04×10+7)10 1.9457×10+11
(1.24×10+10)−3.7692×10+11
(1.60×10+10)−2.2741×10+11
(1.20×10+10)−1.1103×10+11
(8.12×10+9)−1.2033×10+11
(6.96×10+9)−6.8603×10+11
(3.29×10+10)15 3.9863×10+12
(7.29×10+11)−8.8868×10+12
(7.95×10+11)−4.9690×10+12
(3.43×10+11)−2.0497×10+10
(3.24×10+10)−1.4482×10+12
(4.81×10+11)−2.0970×10+13
(1.56×10+12)25 2.8409×10+14
(2.86×10+13)−1.6088×10+15
(7.01×10+14)−2.7765×10+14
(3.00×10+13)−1.1499×10+12
(2.16×10+12)−5.7645×10+14
(6.09×10+13)−1.2220×10+16
(5.00×10+14)WFG6 5 1.1994×10+8
(8.37×10+6)−1.2908×10+8
(1.30×10+7)−1.1997×10+8
(7.99×10+6)−1.7038×10+7
(6.02×10+6)−5.4499×10+7
(6.80×10+6)−1.8108×10+8
(1.37×10+7)10 1.6619×10+11
(1.38×10+10)−2.7497×10+11
(1.56×10+10)−1.9869×10+11
(1.18×10+10)−9.2008×10+10
(6.73×10+9)−1.1142×10+11
(6.52×10+9)−6.2142×10+11
(2.98×10+10)15 3.5945×10+12
(7.09×10+11)−4.3937×10+12
(8.35×10+11)−2.6933×10+12
(5.32×10+11)−4.3260×10+10
(8.15×10+10)−5.8023×10+11
(5.03×10+11)−1.7544×10+13
(1.37×10+12)25 5.8142×10+13
(1.39×10+13)−2.2239×10+15
(4.43×10+14)−1.3831×10+14
(3.12×10+13)−1.0474×10+12
(4.06×10+12)−1.3621×10+14
(5.09×10+13)−1.0081×10+16
(4.97×10+14)WFG7 5 1.4307×10+8
(8.30×10+6)−1.5328×10+8
(1.08×10+7)−1.3594×10+8
(7.42×10+6)−2.0429×10+7
(7.25×10+6)−8.4724×10+7
(5.77×10+6)−2.0010×10+8
(1.05×10+7)10 2.5289×10+11
(1.42×10+10)−3.8153×10+11
(3.44×10+10)−2.6215×10+11
(1.24×10+10)−1.2614×10+11
(9.82×10+9)−1.8152×10+11
(7.90×10+9)−6.9190×10+11
(3.57×10+10)15 5.8545×10+12
(7.41×10+11)−9.3356×10+12
(1.86×10+12)−7.4133×10+12
(1.39×10+12)−6.2349×10+10
(7.43×10+10)−2.2599×10+12
(9.62×10+11)−1.4999×10+13
(1.54×10+12)25 6.8227×10+14
(1.72×10+14)−1.8449×10+15
(4.42×10+14)−7.9199×10+14
(1.16×10+14)−2.3362×10+12
(5.57×10+12)−8.9765×10+14
(1.59×10+14)−7.1165×10+15
(4.06×10+14)WFG8 5 1.5013×10+8
(8.11×10+6)−1.9783×10+8
(1.05×10+7)−1.7829×10+8
(7.32×10+6)−3.6308×10+7
(4.44×10+6)−1.4408×10+8
(7.57×10+6)−2.3335×10+8
(1.10×10+7)10 1.8219×10+11
(1.68×10+10)−3.8567×10+11
(5.74×10+10)−2.0046×10+11
(2.71×10+10)−8.2783×10+10
(3.03×10+10)−1.7126×10+11
(1.80×10+10)−6.0173×10+11
(7.82×10+10)15 4.0400×10+12
(9.81×10+11)−9.9234×10+12
(1.83×10+12)=3.1853×10+12
(9.93×10+11)−1.9227×10+11
(1.86×10+11)−1.6871×10+12
(1.26×10+12)−1.2895×10+13
(5.22×10+12)25 1.3971×10+14
(9.91×10+13)−4.8077×10+15
(6.01×10+14)−2.4789×10+14
(8.47×10+13)−4.7555×10+12
(8.90×10+12)−3.1451×10+14
(9.11×10+13)−1.0680×10+16
(8.26×10+14)WFG9 5 1.9141×10+8
(8.37×10+6)−2.3341×10+8
(1.10×10+7)−1.9389×10+8
(7.79×10+6)−3.6203×10+7
(5.92×10+6)−1.5846×10+8
(6.96×10+6)−2.6391×10+8
(1.35×10+7)10 3.7972×10+11
(3.13×10+10)−6.0898×10+11
(4.27×10+10)−4.2790×10+11
(1.62×10+10)−2.0148×10+11
(1.32×10+10)−3.5255×10+11
(1.94×10+10)−9.8787×10+11
(3.03×10+10)15 1.4912×10+13
(1.61×10+12)−2.0973×10+13
(2.10×10+12)−1.2833×10+13
(1.43×10+12)−5.5018×10+10
(8.57×10+10)−1.0318×10+13
(1.36×10+12)−3.2318×10+13
(1.52×10+12)25 4.2407×10+15
(8.53×10+14)−9.2528×10+15
(1.06×10+15)−3.7876×10+15
(6.04×10+14)−4.6745×10+12
(7.48×10+12)−5.3764×10+15
(3.32×10+14)−1.6297×10+16
(4.08×10+14)+/−/= 4/36/0 3/36/1 4/36/0 1/39/0 3/37/0 “+”表明该算法优于MaOEA-IAR, “−”劣于MaOEA-IAR, “=”则表示与MaOEA-IAR性能相似
下载: 导出CSV
表 6 MOEA-DD、NSGA-III、RVEA、MOMBI-II、AR-MOEA以及MaOEA-IAR在DTLZ5~DTLZ7, IDTLZ1~IDTLZ2, WFG1~WFG3上获得IGD值的统一结果(均值和标准差). 最好的结果已被标记
Table 6 The statistical results (mean and standard deviation) of the IGD values obtained by MOEA-DD、NSGA-III、RVEA、MOMBI-II、AR-MOEA and MaOEA-IAR on DTLZ5~DTLZ7, IDTLZ1~IDTLZ2, WFG1~WFG3. The best results are highlighted
问题 M MOEA-DD NSGA-III RVEA MOMBI-II AR-MOEA MaOEA-IAR DTLZ5 5 1.0286×10−1
(4.56×10−3)−1.5564×10−1
(6.36×10−2)−1.7052×10−1
(2.08×10−2)−2.6836×10−1
(1.37×10−2)−1.2293×10−1
(2.95×10−3)−9.2396×10−2
(1.58×10−2)10 1.3275×10−1
(1.77×10−2)+2.9646×10−1
(9.03×10−2)+3.7614×10−1
(6.44×10−2)+5.8454×10−1
(1.64×10−1)−9.4984×10−2
(3.11×10−3)+4.8043×10−1
(1.04×10−1)15 1.5268×10−1
(4.95×10−3)−2.5900×10−1
(4.20×10−2)−5.8957×10−1
(1.74×10−1)−6.4004×10−1
(6.56×10−2)−6.0343×10−1
(1.01×10−2)−1.3313×10−1
(1.25×10−2)25 1.4321×10−1
(4.38×10−2)+8.5001×10−1
(5.98×10−1)−3.6054×10−1
(2.03×10−1)−5.7672×10−1
(2.97×10−1)−9.6898×10−2
(1.91×10−2)−3.2482×10−2
(1.69×10−1)DTLZ6 5 1.0630×10−1
(9.13×10−3)+2.8728×10−1
(2.28×10−1)+1.5455×10−1
(3.80×10−2)+3.8061×10−1
(2.25×10−1)=3.3734×10−1
(3.50×10−2)−3.1149×10−1
(5.18×10−2)10 2.0616×10−1
(1.64×10−1)+1.5928×10+0
(5.87×10−1)+2.8394×10−1
(8.30×10−2)+5.4111×10−1
(2.36×10−1)+1.6497×10−1
(8.66×10−3)+2.4129×10+0
(3.39×10−1)15 2.0227×10−1
(1.34×10−1)+2.5854×10+0
(1.02×10+0)=3.2020×10−1
(1.73×10−1)+7.5239×10−1
(2.82×10−1)+1.2810×10−1
(1.01×10−1)+2.9171×10+0
(4.24×10−1)25 1.5709×10−1
(2.14×10−2)−6.1820×10+0
(3.10×10+0)−2.2612×10−1
(6.64×10−2)−5.4816×10−1
(2.83×10−1)−5.0812×10−1
(3.00×10−1)−1.0126×10−1
(2.94×10−2)DTLZ7 5 3.0002×10+0
(9.72×10−4)−3.3903×10−1
(1.25×10−2)+6.0021×10−1
(3.12×10−2)−4.8750×10−1
(8.74×10−2)−8.7917×10−3
(4.95×10−4)+4.0035×10−1
(8.59×10−2)10 2.0394×10+0
(3.09×10−1)−1.6480×10+0
(2.43×10−1)−1.8321×10+0
(4.74×10−1)−4.8283×10+0
(7.64×10−1)−2.8392×10−2
(3.28×10−3)+8.4705×10−1
(3.02×10−2)15 3.4630×10+0
(8.21×10−2)−8.8601×10+0
(7.49×10−1)−2.3863×10+0
(2.46×10−1)−1.1222×10+1
(1.43×10−1)−5.4456×10−1
(3.13×10−1)+1.6482×10+0
(3.88×10−2)25 4.9987×10+0
(4.01×10−2)−1.4671×10+1
(1.41×10+0)−3.5010×10+0
(4.87×10−1)+1.9012×10+1
(2.83×10−1)−3.7797×10+0
(1.17×10+0)+4.0533×10+0
(1.41×10−5)IDTLZ1 5 1.5305×10−1
(4.10×10−2)=1.3547×10−1
(1.54×10−2)=1.5234×10−1
(1.70×10−2)=1.1355×10−1
(3.26×10−4)=6.5841×10−2
(8.22×10−4)+2.5128×10−1
(1.65×10−1)10 2.2790×10−1
(1.05×10−2)−1.4175×10−1
(3.27×10−3)−2.5816×10−1
(2.39×10−2)−1.8633×10−1
(6.51×10−3)−1.2386×10−1
(3.42×10−2)=1.1027×10−1
(9.23×10−4)15 3.2984×10−1
(3.09×10−2)−2.1719×10−1
(4.73×10−3)−3.5590×10−1
(2.89×10−2)−2.1098×10−1
(1.14×10−2)−1.9277×10−1
(1.29×10−2)−1.7215×10−1
(7.39×10−2)25 3.2145×10−1
(1.81×10−2)−1.9936×10−1
(3.93×10−3)−3.8507×10−1
(3.13×10−2)−2.6360×10−1
(6.52×10−3)−4.3902×10−2
(1.27×10−1)−2.1252×10−2
(1.63×10−2)IDTLZ2 5 2.7867×10−1
(5.70×10−3)−2.4322×10−1
(6.88×10−3)−2.9833×10−1
(5.13×10−3)−3.1740×10−1
(1.15×10−3)−2.1362×10−1
(5.25×10−3)−2.0888×10−1
(7.30×10−3)10 7.3889×10−1
(6.96×10−3)−6.0395×10−1
(1.69×10−2)−6.7982×10−1
(7.75×10−3)−6.6946×10−1
(4.82×10−3)−4.4639×10−1
(6.90×10−3)−4.4152×10−1
(6.68×10−3)15 9.4524×10−1
(2.03×10−2)−7.6174×10−1
(1.42×10−2)−8.7051×10−1
(1.47×10−2)−8.5834×10−1
(7.35×10−3)−6.8015×10−1
(1.79×10−2)−5.7753×10−1
(9.19×10−3)25 1.1245×10+0
(1.07×10−2)−8.3706×10−1
(1.63×10−2)−1.0678×10+0
(1.60×10−2)−1.0151×10+0
(4.12×10−3)−8.9235×10−2
(2.87×10−2)-2.8098×10−2
(6.87×10−4)WFG1 5 8.4034×10−1
(1.15×10−1)−5.6876×10−1
(4.95×10−2)−5.3556×10−1
(4.24×10−2)−5.3968×10−1
(6.19×10−2)−5.1692×10−1
(1.26×10−2)−4.4951×10−1
(5.18×10−3)10 1.2623×10+0
(6.70×10−2)−1.2030×10+0
(7.30×10−2)−1.0736×10+0
(6.23×10−2)=1.2539×10+0
(5.56×10−2)−1.1573×10+0
(3.43×10−2)−1.0186×10+0
(3.21×10−2)15 1.9032×10+0
(6.34×10−2)−1.9456×10+0
(8.90×10−2)−1.8780×10+0
(7.27×10−2)−2.5378×10+0
(2.18×10−1)−1.7899×10+0
(4.09×10−2)=1.7788×10+0
(4.73×10−2)25 3.8022×10+0
(5.48×10−2)−3.0121×10+0
(4.19×10−1)−3.1303×10+0
(1.12×10−1)−3.6938×10+0
(9.71×10−2)−3.1942×10+0
(6.67×10−2)−2.6916×10+0
(1.82×10−1)WFG2 5 5.7666×10−1
(1.64×10−2)−4.6948×10−1
(3.19×10−3)+4.4930×10−1
(1.01×10−2)+5.1672×10−1
(6.95×10−2)+4.7596×10−1
(2.55×10−3)+5.4316×10−1
(2.97×10−2)10 1.4487×10+0
(2.06×10−2)−1.2010×10+0
(1.43×10−1)−1.1062×10+0
(3.99×10−2)−1.6478×10+0
(5.00×10−1)−1.4937×10+0
(4.68×10−2)−1.0054×10+0
(1.43×10−2)15 2.1614×10+0
(3.30×10−2)−1.7688×10+0
(7.55×10−2)=1.7815×10+0
(1.09×10−1)=7.5478×10+0
(2.36×10+0)−1.7084×10+0
(3.82×10−2)=1.7263×10+0
(1.03×10−1)25 4.0448×10+0
(1.16×10−2)−3.5543×10+0
(1.21×10−1)−2.8225×10+0
(1.52×10−1)−1.0109×10+1
(5.56×10+0)−3.1864×10+0
(1.60×10−1)−2.7976×10+0
(1.71×10−1)WFG3 5 7.4442×10−1
(3.28×10−2)−5.9020×10−1
(5.58×10−2)−6.8641×10−1
(1.11×10−1)−1.6953×10+0
(1.28×10−1)−6.8013×10−1
(1.54×10−1)−5.4840×10−1
(3.21×10−2)10 2.7654×10+0
(1.17×10−1)+1.4175×10+0
(4.13×10−1)+3.4602×10+0
(5.56×10−1)−8.5627×10+0
(1.89×10−1)−2.4166×10+0
(9.05×10−2)+3.0447×10+0
(5.68×10−1)15 6.4307×10+0
(6.67×10−1)+2.6293×10+0
(4.06×10−1)+6.6442×10+0
(1.19×10+0)+1.3900×10+1
(2.03×10−1)−5.5924×10+0
(2.08×10−1)+7.5080×10+0
(9.65×10−1)25 1.8160×10+1
(3.63×10−2)−1.2519×10+1
(1.40×10+0)−1.1530×10+1
(1.65×10+0)−2.7424×10+1
(3.33×10−2)−9.8933×10+0
(1.41×10−1)−6.8891×10+0
(5.11×10+0)+/−/= 7/24/1 7/22/3 7/22/3 3/27/2 11/18/3 “+”表明该算法优于MaOEA-IAR, “−”劣于MaOEA-IAR, “=”则表示与MaOEA-IAR性能相似
下载: 导出CSV
表 7 MOEA-DD、NSGA-III、RVEA、MOMBI-II、AR-MOEA以及MaOEA-IAR在DTLZ5~DTLZ7, IDTLZ1~IDTLZ2, WFG1~WFG3上获得PD值的统一结果(均值和标准差).最好的结果已被标记
Table 7 The statistical results (mean and standard deviation) of the PD values obtained by MOEA-DD、NSGA-III、RVEA、MOMBI-II、AR-MOEA and MaOEA-IAR on DTLZ5~DTLZ7, IDTLZ1~IDTLZ2, WFG1~WFG3. The best results are highlighted
问题 M MOEA-DD NSGA-III RVEA MOMBI-II AR-MOEA MaOEA-IAR DTLZ5 5 1.7925×10+7
(2.11×10+6)−3.5995×10+7
(4.18×10+6)+1.8932×10+7
(4.20×10+6)−1.1044×10+7
(1.59×10+6)−4.8118×10+7
(2.34×10+6)+2.5713×10+7
(2.69×10+6)10 1.9742×10+10
(2.54×10+9)−7.1635×10+10
(5.54×10+9)−2.9432×10+10
(3.78×10+9)−2.1484×10+9
(2.81×10+9)−8.3986×10+10
(3.99×10+9)−1.2087×10+11
(8.00×10+9)15 2.8397×10+11
(7.71×10+10)−1.1841×10+12
(1.38×10+11)−4.1423×10+10
(2.81×10+10)−6.5840×10+9
(4.48×10+9)−1.0895×10+12
(1.09×10+11)−2.9251×10+12
(2.88×10+11)25 4.2895×10+13
(9.57×10+12)+1.4562×10+11
(1.05×10+11)−1.7203×10+13
(1.37×10+13)+4.0293×10+10
(4.56×10+10)−1.0108×10+14
(1.23×10+13)−5.6358×10+14
(2.21×10+14)DTLZ6 5 2.8362×10+7
(3.50×10+6)−5.4971×10+7
(7.11×10+6)+2.6105×10+7
(7.26×10+6)−1.0490×10+7
(2.67×10+6)−5.2864×10+7
(3.81×10+6)+3.1205×10+7
(4.21×10+6)10 6.5293×10+10
(1.58×10+10)−1.8747×10+11
(4.95×10+10)−5.9078×10+10
(1.80×10+10)−2.1985×10+9
(2.83×10+9)−8.4598×10+10
(5.75×10+9)−3.9921×10+11
(2.70×10+10)15 9.1234×10+11
(3.10×10+11)−3.4359×10+12
(5.79×10+11)−7.4856×10+11
(3.95×10+11)−1.0372×10+10
(1.40×10+10)−4.2785×10+11
(3.48×10+11)−7.9758×10+12
(8.70×10+11)25 2.2080×10+14
(1.01×10+14)−2.2455×10+14
(1.56×10+13)−8.7904×10+13
(2.72×10+13)−1.5169×10+11
(1.15×10+11)−3.2540×10+13
(1.77×10+13)−2.3216×10+15
(1.08×10+15)DTLZ7 5 1.6766×10+3
(4.84×10+3)−2.1762×10+7
(4.10×10+6)=1.8563×10+7
(3.07×10+6)=5.1253×10+6
(7.63×10+5)−3.3243×10+7
(4.17×10+6)+2.0833×10+7
(7.02×10+6)10 2.1457×10+9
(1.38×10+9)−3.1794×10+10
(2.69×10+9)−1.6034×10+10
(4.07×10+9)−8.0445×10+9
(1.73×10+9)−2.5524×10+10
(3.65×10+9)−6.8683×10+10
(7.18×10+9)15 1.8990×10+10
(5.95×10+9)−4.7550×10+11
(9.43×10+10)−2.7698×10+11
(6.40×10+10)−4.9554×10+10
(2.76×10+10)−4.0328×10+11
(1.19×10+11)−2.5059×10+12
(3.10×10+11)25 4.2556×10+12
(1.02×10+12)−7.3834×10+13
(1.20×10+13)−4.7217×10+13
(1.00×10+13)−4.4333×10+13
(1.34×10+13)−1.3427×10+9
(6.62×10+8)−1.1453×10+14
(1.61×10+13)IDTLZ1 5 1.4745×10+6
(9.43×10+5)−6.3366×10+6
(3.40×10+6)−2.1610×10+6
(9.83×10+5)−8.4010×10+5
(2.39×10+5)−1.0739×10+7
(2.67×10+6)−1.9799×10+8
(2.32×10+8)10 1.0063×10+8
(1.13×10+8)−9.8146×10+9
(1.30×10+9)−1.0241×10+9
(9.00×10+8)−2.0726×10+9
(3.55×10+8)−7.1952×10+9
(1.15×10+9)−4.5218×10+10
(1.18×10+11)15 4.8749×10+10
(2.39×10+11)−2.7098×10+11
(4.55×10+10)−6.2342×10+10
(5.79×10+10)−3.9405×10+10
(8.72×10+9)−1.2277×10+11
(3.42×10+11)−3.2538×10+12
(1.06×10+13)25 3.8837×10+9
(8.47×10+9)−5.6426×10+13
(1.28×10+13)−7.6673×10+12
(5.54×10+12)−1.6889×10+12
(1.19×10+12)−1.6465×10+14
(4.61×10+14)−5.8266×10+14
(1.42×10+13)IDTLZ2 5 1.2432×10+7
(1.92×10+6)−3.5278×10+7
(7.22×10+6)−1.8123×10+7
(1.50×10+6)−4.8079×10+6
(6.93×10+5)−2.6241×10+7
(1.72×10+6)−5.5348×10+7
(2.41×10+6)10 3.7649×10+9
(3.38×10+8)−4.9622×10+10
(7.67×10+9)−1.1533×10+10
(1.30×10+9)−1.5701×10+10
(1.85×10+9)−9.1979×10+10
(3.71×10+9)−1.6089×10+11
(4.17×10+9)15 7.1923×10+10
(3.07×10+10)−1.0779×10+12
(1.50×10+11)−2.0643×10+11
(6.01×10+10)−2.9274×10+11
(4.07×10+10)−2.6072×10+12
(1.61×10+11)−5.4279×10+12
(1.82×10+11)25 4.4848×10+12
(1.99×10+12)−5.6171×10+14
(8.35×10+13)−1.4306×10+13
(8.54×10+12)−7.0067×10+13
(1.05×10+13)−8.5860×10+14
(4.94×10+13)−2.0550×10+15
(5.47×10+13)WFG1 5 7.7570×10+7
(4.76×10+6)+7.5733×10+7
(7.14×10+6)+6.6169×10+7
(5.02×10+6)+7.6250×10+6
(2.77×10+6)−3.0367×10+7
(2.17×10+6)+2.3345×10+7
(3.37×10+6)10 3.8144×10+10
(3.75×10+9)+8.4159×10+10
(8.64×10+9)+4.6067×10+10
(6.62×10+9)+6.4761×10+9
(1.55×10+9)−4.3604×10+10
(2.26×10+9)+1.8476×10+10
(2.80×10+9)15 1.3448×10+12
(2.24×10+11)+1.2676×10+12
(5.92×10+11)+6.2449×10+11
(1.04×10+11)=1.9369×10+10
(2.14×10+10)−6.0993×10+11
(4.26×10+10)=6.1772×10+11
(1.26×10+11)25 3.3402×10+13
(2.45×10+12)−1.6094×10+14
(4.13×10+13)+1.0889×10+14
(5.98×10+12)+4.8564×10+12
(1.42×10+12)−1.0664×10+14
(3.18×10+12)+8.4524×10+13
(1.67×10+13)WFG2 5 5.1237×10+7
(2.65×10+6)−5.4016×10+7
(3.16×10+6)−6.1875×10+7
(2.87×10+6)−1.1742×10+7
(3.17×10+6)−3.8312×10+7
(2.55×10+6)−7.6371×10+7
(4.16×10+6)10 3.9168×10+10
(1.43×10+9)−6.6429×10+10
(1.01×10+10)−5.3106×10+10
(2.41×10+9)−7.4067×10+9
(5.52×10+9)−3.9309×10+10
(2.97×10+9)−9.1141×10+10
(4.59×10+9)15 4.4512×10+11
(4.65×10+10)−1.5952×10+12
(2.02×10+11)−9.2321×10+11
(1.58×10+11)−1.2053×10+10
(1.52×10+10)−5.7677×10+11
(8.49×10+10)−2.2461×10+12
(1.65×10+11)25 4.1805×10+13
(3.61×10+12)−3.8842×10+14
(1.07×10+14)−1.7542×10+14
(2.40×10+13)−3.8888×10+11
(9.15×10+11)−1.6560×10+14
(1.30×10+13)−4.7295×10+14
(2.67×10+13)WFG3 5 1.2116×10+8
(9.07×10+6)−1.4756×10+8
(7.56×10+6)−1.3033×10+8
(1.40×10+7)−5.1396×10+7
(4.21×10+6)−1.5433×10+8
(7.28×10+6)=1.5501×10+8
(6.95×10+6)10 1.1096×10+11
(9.29×10+9)−2.4982×10+11
(2.39×10+10)−1.6581×10+11
(1.63×10+10)−9.4558×10+8
(1.51×10+8)−1.9220×10+11
(1.05×10+10)−3.1718×10+11
(2.00×10+10)15 3.8028×10+12
(7.77×10+11)−7.4227×10+12
(1.37×10+12)+4.1520×10+12
(7.83×10+11)−1.6200×10+10
(5.00×10+9)−2.7152×10+12
(9.23×10+11)−6.1382×10+12
(5.23×10+11)25 1.0941×10+14
(6.89×10+12)−1.9873×10+15
(6.87×10+14)=5.6632×10+14
(1.69×10+14)−5.1563×10+10
(1.66×10+10)−3.2575×10+14
(2.26×10+13)−1.5572×10+15
(1.06×10+14)+/−/= 4/28/0 7/23/2 4/26/2 0/32/0 6/24/2 “+”表明该算法优于MaOEA-IAR, “−”劣于MaOEA-IAR, “=”则表示与MaOEA-IAR性能相似
下载: 导出CSV
表 8 AR-MOEA和MaOEA-IAR在DTLZ1~DTLZ7, IDTLZ1~IDTLZ2, WFG1~WFG9上运行时间的统一结果(均值). 最好的结果已被标记
Table 8 The statistical results (mean) of the time obtained by AR-MOEA and MaOEA-IAR on DTLZ1~DTLZ7, IDTLZ1~IDTLZ2, WFG1~WFG9. The best results are highlighted
问题 M AR-MOEA MaOEA-IAR 问题 M AR-MOEA MaOEA-IAR DTLZ1 5 8.7284×10+1− 5.2877×10+1 WFG1 5 2.1497×10+2− 1.2453×10+2 10 3.2368×10+3− 1.4929×10+3 10 3.0802×10+3− 1.8838×10+3 15 1.5646×10+2+ 4.5986×10+2 15 3.3103×10+2− 2.9192×10+2 25 2.7908×10+3− 1.8621×10+3 25 2.4172×10+3+ 2.8807×10+3 DTLZ2 5 1.4045×10+2+ 3.1151×10+2 WFG2 5 2.5691×10+2− 1.9505×10+2 10 3.4692×10+3− 2.0823×10+3 10 3.6336×10+3− 3.2122×10+3 15 1.9603×10+2+ 4.9719×10+2 15 4.3556×10+2− 3.8969×10+2 25 3.1545×10+3+ 3.4599×10+3 25 3.4466×10+3= 3.3881×10+3 DTLZ3 5 7.6826×10+1− 3.3370×10+1 WFG3 5 3.4018×10+2− 3.2435×10+2 10 3.2959×10+3− 1.5038×10+3 10 4.3267×10+3− 4.1139×10+3 15 4.9418×10+2− 1.5822×10+2 15 5.1973×10+2− 4.3731×10+2 25 2.9858×10+3− 1.7848×10+3 25 3.8047×10+3= 3.7357×10+3 DTLZ4 5 1.4659×10+2+ 3.2766×10+2 WFG4 5 3.3620×10+2− 2.8286×10+2 10 3.8768×10+3− 2.1238×10+3 10 4.2819×10+3− 3.6326×10+3 15 1.9815×10+2+ 5.1050×10+2 15 3.8595×10+2+ 5.2267×10+2 25 3.1269×10+3+ 3.2684×10+3 25 3.7546×10+3− 3.3598×10+3 DTLZ5 5 2.2500×10+2− 1.3912×10+2 WFG5 5 3.3378×10+2− 2.9287×10+2 10 3.0904×10+3− 1.9358×10+3 10 4.2697×10+3− 3.7511×10+3 15 1.8227×10+2+ 4.5443×10+2 15 5.2633×10+2− 4.0438×10+2 25 3.0757×10+3− 1.9646×10+3 25 3.8177×10+3− 3.4843×10+3 DTLZ6 5 1.9591×10+2− 1.2208×10+2 WFG6 5 2.7412×10+2− 2.5095×10+2 10 3.5852×10+3− 1.8312×10+3 10 4.0843×10+3− 3.2581×10+3 15 1.5310×10+2+ 4.9933×10+2 15 3.1397×10+2+ 5.1200×10+2 25 2.6990×10+3− 2.0808×10+3 25 3.7598×10+3− 2.8335×10+3 DTLZ7 5 2.0799×10+2− 1.4207×10+2 WFG7 5 3.7091×10+2− 3.1761×10+2 10 3.8714×10+3− 2.3062×10+3 10 4.3950×10+3− 3.7082×10+3 15 4.8328×10+2− 2.0122×10+2 15 3.7417×10+2+ 5.2910×10+2 25 2.4555×10+3− 1.5033×10+3 25 3.8242×10+3− 3.2364×10+3 IDTLZ1 5 2.1931×10+2− 9.4000×10+1 WFG8 5 2.2219×10+2= 2.1770×10+2 10 3.1281×10+3− 1.9947×10+3 10 3.7129×10+3− 2.6444×10+3 15 3.4287×10+2− 2.5752×10+2 15 2.4552×10+2+ 4.5493×10+2 25 2.0253×10+3− 1.2549×10+3 25 3.7527×10+3− 2.5151×10+3 IDTLZ2 5 3.1594×10+2+ 3.3723×10+2 WFG9 5 3.8469×10+2− 3.2047×10+2 10 4.0353×10+3− 3.4738×10+3 10 4.3729×10+3− 3.9 448×10+3 15 4.4610×10+2− 4.2203×10+2 15 5.3114×10+2− 4.3118×10+2 25 2.5265×10+3− 1.7758×10+3 25 3.8142×10+3− 3.5823×10+3 +/−/= 10/26/0 +/−/= 5/28/3 “+”表明该算法优于MaOEA-IAR, “−”劣于MaOEA-IAR, “=”则表示与MaOEA-IAR性能相似
下载: 导出CSV
表 9 AR-MOEA和MaOEA-IAR在MaF1~MaF7上获得IGD值的统一结果(均值和标准差). 最好的结果已被标记
Table 9 The statistical results (mean and standard deviation) of the IGD values obtained by AR-MOEA and MaOEA-IAR on MaF1~MaF7. The best results are highlighted
问题 M AR-MOEA MaOEA-IAR MaF1 5 6.1773×10−2 (6.75×10−2)− 3.2287×10−3 (7.05×10−5) 10 1.0257×10−1 (1.11×10−1)− 6.5011×10−3 (1.44×10−4) 15 1.6987×10−1 (1.85×10−1)− 1.1253×10−2 (2.95×10−4) 25 1.7434×10−1 (1.89×10−1)− 9.2626×10−3 (1.97×10−4) MaF2 5 5.6263×10−2 (5.26×10−2)− 8.5071×10−3 (3.31×10−4) 10 8.3280×10−2 (8.77×10−2)− 6.4968×10−3 (9.61×10−5) 15 1.2184×10−1 (1.24×10−1)− 8.9506×10−3 (2.26×10−4) 25 1.0515×10−1 (1.21×10−1)− 6.2938×10−3 (4.48×10−4) MaF3 5 3.4660×10−2 (4.22×10−2)− 1.0505×10−3 (8.53×10−4) 10 3.4304×10−2 (4.06×10−2)− 4.8646×10−3 (8.72×10−3) 15 1.6992×10+0 (6.87×10+0)− 3.9708×10−3 (3.15×10−3) 25 8.5024×10−1 (2.47×10+0)− 2.0315×10−3 (2.00×10−3) MaF4 5 1.0234×10+0 (1.26×10+0)− 4.5541×10−2 (9.87×10−3) 10 3.6718×10+1 (4.35×10+1)− 1.7669×10+0 (1.34×10−2) 15 1.7572×10+3 (2.01×10+3)− 2.7158×10+2 (5.78×10−1) 25 1.8804×10+6 (2.17×10+6)− 1.7062×10+5 (3.15×10+2) MaF5 5 1.0214×10+0 (1.29×10+0)− 6.3123×10−2 (1.23×10−2) 10 3.7332×10+1 (4.78×10+1)− 6.3090×10−1 (9.64×10−2) 15 1.4852×10+3 (1.90×10+3)− 3.6171×10+0 (1.29×10+0) 25 1.7966×10+6 (2.32×10+6)− 7.1131×10+0 (8.87×10+0) MaF6 5 1.3381×10−3 (1.72×10−3)− 4.2532×10−6 (1.43×10−6) 10 3.1054×10−1 (1.30×10+0)+ 4.8255×10+0 (3.33×10+0) 15 2.1990×10−1 (7.89×10−1)− 4.6076×10−6 (1.37×10−6) 25 3.5425×10−1 (5.51×10−1)− 3.7808×10−6 (4.08×10−6) MaF7 5 1.3383×10−1 (1.61×10−1)= 1.8586×10−1 (2.35×10−1) 10 5.7112×10−1 (6.97×10−1)= 3.5443×10−1 (4.13×10−1) 15 2.2371×10+0 (2.25×10+0)− 7.1977×10−1 (7.75×10−1) 25 3.6581×10+0 (3.86×10+0)− 1.1867×10−1 (8.44×10−1) +/−/= 1/25/2 “+”表明该算法优于MaOEA-IAR, “−”劣于MaOEA-IAR, “=”则表示与MaOEA-IAR性能相似
下载: 导出CSV
表 10 不同参数下算法MaOEA-IAR的性能
Table 10 Performance of MaOEA-IAR under different parameters
问题 目标数 r初始值 5 7 9 11 13 DTLZ1 5 2.6870×10−1 2.6135×10−1 2.6051×10−1 2.6382×10−1 2.6351×10−1 10 4.8177×10+1 4.7009×10+1 4.6085×10+1 4.7519×10+1 4.7286×10+1 15 5.7026×10+1 5.6981×10+1 5.5765×10+1 5.6146×10+1 5.6083×10+1 25 4.5156×10−1 4.5013×10−1 3.6557×10−1 2.4067×10−1 2.4052×10−1 DTLZ5 5 9.2396×10−2 8.3615×10−2 8.0045×10−2 8.8463×10−2 8.8732×10−2 10 6.2597×10−1 4.7625×10−1 4.7044×10−1 4.7765×10−1 4.7829×10−1 15 1.3549×10−1 1.3158×10−1 1.2994×10−1 1.3036×10−1 1.3023×10−1 25 3.3464×10−2 3.3391×10−2 3.3301×10−2 3.3362×10−2 3.3485×10−2 IDTLZ1 5 2.5128×10−1 2.3743×10−1 2.1694×10−1 2.2748×10−1 2.2926×10−1 10 1.1236×10−1 1.1149×10−1 1.0017×10−1 1.0594×10−1 1.0542×10−1 15 1.8035×10−1 1.7891×10−1 1.7077×10−1 1.7901×10−1 1.7924×10−1 25 2.2693×10−2 2.2451×10−2 2.2279×10−2 2.2538×10−2 2.2523×10−2 WFG1 5 4.4951×10−1 4.3826×10−1 4.3007×10−1 4.3596×10−1 4.4001×10−1 10 1.1967×10+0 1.1051×10+0 1.0186×10+0 1.0598×10+0 1.0482×10+0 15 1.6956×10+0 1.5941×10+0 1.5076×10+0 1.6582×10+0 1.6677×10+0 25 2.8015×10+0 2.7693×10+0 2.6724×10+0 2.7795×10+0 2.7619×10+0 WFG2 5 5.4316×10−1 5.4029×10−1 5.4002×10−1 5.4174×10−1 5.4169×10−1 10 1.1086×10+0 1.0997×10+0 1.0051×10+0 1.1036×10+0 1.1097×10+0 15 1.8675×10+0 1.7649×10+0 1.7069×10+0 1.8007×10+0 1.7952×10+0 25 2.9598×10+0 2.9413×10+0 2.8039×10+0 2.8796×10+0 2.8976×10+0 WFG4 5 1.1192×10+0 1.1128×10+0 1.1065×10+0 1.1097×10+0 1.1106×10+0 10 4.1062×10+0 4.1057×10+0 4.0388×10+0 4.0263×10+0 4.0589×10+0 15 8.9645×10+0 8.8089×10+0 8.3986×10+0 8.4108×10+0 8.4004×10+0 25 1.9058×10+1 1.8769×10+1 1.6506×10+1 1.5375×10+1 1.6405×10+1
下载: 导出CSV
-
[1] Luo J P, Liu Q Q, Yang Y, Li X, Chen M R, Cao W M. An artificial bee colony algorithm for multi-objective optimisation. Applied Soft Computing, 2017, 50: 235-251 doi: 10.1016/j.asoc.2016.11.014 [2] Zhou Y L, Wang J H, Chen J, Gao S C, Teng L Y. Ensemble of many-objective evolutionary algorithms for many-objective problems. Soft Computing, 2017, 21(9): 2407-2419 doi: 10.1007/s00500-015-1955-3 [3] Jain H, Deb K. An evolutionary many-objective optimization algorithm using reference-point based nondominated sorting approach, part II: Handling constraints and extending to an adaptive approach. IEEE Transactions on Evolutionary Computation, 2014, 18(4): 602-622 doi: 10.1109/TEVC.2013.2281534 [4] Cinalli D, Martí L, Sanchez-Pi N, Garcia A C B. Collective intelligence approaches in interactive evolutionary multi-objective optimization. Logic Journal of the IGPL, 2020, 28(1): 95-108 doi: 10.1093/jigpal/jzz074 [5] Kudikala R, Mills A R, Fleming P J, Tanner G F, Holt J E. Aero engine health management system architecture design using multi-criteria optimization. In: Proceedings of the 15th Annual Conference Companion on Genetic and Evolutionary Computation. Amsterdam, The Netherlands: ACM, 2013. 185−186 [6] Saeidi A M, Hage J, Khadka R, Jansen S. A search-based approach to multi-view clustering of software systems. In: Proceedings of the 22nd IEEE International Conference on Software Analysis, Evolution, and Reengineering. Montreal, Canada: IEEE, 2015. 429−438 [7] 陈振兴, 严宣辉, 吴坤安, 白猛. 融合张角拥挤控制策略的高维多目标优化. 自动化学报, 2015, 41(6): 1145-1158Chen Zhen-Xing, Yan Xuan-Hui, Wu Kun-An, Bai Meng. Many-objective optimization integrating open angle based congestion control strategy. Acta Automatica Sinica, 2015, 41(6): 1145-1158 [8] Liu Y, Zhu N B, Li K L, Li M Q, Zheng J H, Li K Q. An angle dominance criterion for evolutionary many-objective optimization. Information Sciences, 2020, 509: 376-399 doi: 10.1016/j.ins.2018.12.078 [9] Yang S X, Li M Q, Liu X H, Zheng J H. A grid-based evolutionary algorithm for many-objective optimization. IEEE Transactions on Evolutionary Computation, 2013, 17(5): 721-736 doi: 10.1109/TEVC.2012.2227145 [10] 余伟伟, 谢承旺, 闭应洲, 夏学文, 李雄, 任柯燕, 等. 一种基于自适应模糊支配的高维多目标粒子群算法. 自动化学报, 2018, 44(12): 2278-2289Yu Wei-Wei, Xie Cheng-Wang, Bi Ying-Zhou, Xia Xue-Wen, Li Xiong, Ren Ke-Yan, et al. Many-objective particle swarm optimization based on adaptive fuzzy dominance. Acta Automatica Sinica, 2018, 44(12): 2278-2289 [11] Elarbi M, Bechikh S, Gupta A, Said L B, Ong Y S. A new decomposition-based NSGA-II for many-objective optimization. IEEE Transactions on Systems, Man, and Cybernetics: Systems, 2018, 48(7): 1191-1210 doi: 10.1109/TSMC.2017.2654301 [12] Zhang Q F, Li H. MOEA/D: A multiobjective evolutionary algorithm based on decomposition. IEEE Transactions on Evolutionary Computation, 2007, 11(6): 712-731 doi: 10.1109/TEVC.2007.892759 [13] Wang L P, Xu M N, Yu W, Qiu Q C, Wu F. Decomposition multi-objective evolutionary algorithm based on adaptive neighborhood adjustment strategy. IEEE Access, 2020, 8: 78639-78651. doi: 10.1109/ACCESS.2020.2990193 [14] Cheng R, Jin Y C, Olhofer M, Sendhoff B. A reference vector guided evolutionary algorithm for many-objective optimization. IEEE Transactions on Evolutionary Computation, 2016, 20(5): 773-791 doi: 10.1109/TEVC.2016.2519378 [15] Lu X, Tan Y Y, Zheng W, Meng L L. A decomposition method based on random objective division for MOEA/D in many-objective optimization. IEEE Access, 2020, 8: 103550-103564 doi: 10.1109/ACCESS.2020.2999417 [16] Deb K, Jain H. An evolutionary many-objective optimization algorithm using reference-point-based nondominated sorting approach, part I: Solving problems with box constraints. IEEE Transactions on Evolutionary Computation, 2014, 18(4): 577-601 doi: 10.1109/TEVC.2013.2281535 [17] 封文清, 巩敦卫. 基于在线感知Pareto前沿划分目标空间的多目标进化优化. 自动化学报, 2020, 46(8): 1628-1643Feng Wen-Qing, Gong Dun-Wei. Multi-objective evolutionary optimization with objective space partition based on online perception of Pareto front. Acta Automatica Sinica, 2020, 46(8): 1628-1643 [18] Beume N, Naujoks B, Emmerich M. SMS-EMOA: Multiobjective selection based on dominated hypervolume. European Journal of Operational Research, 2007, 181(3): 1653-1669 doi: 10.1016/j.ejor.2006.08.008 [19] Gómez R H, Coello C A C. Improved metaheuristic based on the R2 indicator for many-objective optimization. In: Proceedings of the 2015 Annual Conference on Genetic and Evolutionary Computation. Madrid, Spain: ACM, 2015. 679−686 [20] Liang Z P, Luo T T, Hu K F, Ma X L, Zhu Z X. An indicator-based many-objective evolutionary algorithm with boundary protection. IEEE Transactions on Cybernetics, 2021, 51(9): 4553-4566 doi: 10.1109/TCYB.2019.2960302 [21] Li F, Li T J, Zhang S N. R2 indicator and objective space partition based evolutionary algorithm for many-objective optimization. In: Proceedings of the 2019 IEEE Symposium Series on Computational Intelligence. Xiamen, China: IEEE, 2019. 1271−1278 [22] Sun Y A, Yen G G, Yi Z. IGD indicator-based evolutionary algorithm for many-objective optimization problems. IEEE Transactions on Evolutionary Computation, 2019, 23(2): 173-187 doi: 10.1109/TEVC.2018.2791283 [23] Hua Y C, Liu Q Q, Hao K R, Jin Y C. A survey of evolutionary algorithms for multi-objective optimization problems with irregular Pareto fronts. IEEE/CAA Journal of Automatica Sinica, 2021, 8(2): 303-318 doi: 10.1109/JAS.2021.1003817 [24] Ma X L, Yu Y A, Li X D, Qi Y T. A survey of weight vector adjustment methods for decomposition-based multiobjective evolutionary algorithms. IEEE Transactions on Evolutionary Computation, 2020, 24(4): 634-649 doi: 10.1109/TEVC.2020.2978158 [25] Yuan J W, Liu H L, Gu F Q, Zhang Q F, He Z S. Investigating the properties of indicators and an evolutionary many-objective algorithm using promising regions. IEEE Transactions on Evolutionary Computation, 2021, 25(1): 75-86 doi: 10.1109/TEVC.2020.2999100 [26] Liu S B, Lin Q Z, Tan K C, Gong M G, Coello C A C. A fuzzy decomposition-based multi/many-objective evolutionary algorithm. IEEE Transactions on Cybernetics, to be published [27] Van D A, Gary V, Lamont B. Multiobjective evolutionary algorithm research: A history and analysis. EvolutionaryComputation, 1998, 8(2): 125−147 [28] Zhou A M, Jin Y C, Zhang Q F, Sendhoff B, Tsang E. Combining model-based and genetics-based offspring generation for multi-objective optimization using a convergence criterion. In: Proceedings of the 2006 IEEE International Conference on Evolutionary Computation. Vancouver, Canada: IEEE, 2006. 892−899 [29] While L, Hingston P, Barone L, Huband S. A faster algorithm for calculating hypervolume. IEEE Transactions on Evolutionary Computation, 2006, 10(1): 29-38 doi: 10.1109/TEVC.2005.851275 [30] Liu C, Zhao Q, Yan B, Gao Y. A new hypervolume-based differential evolution algorithm for many-objective optimization. RAIRO-Operations Research, 2017, 51(4): 1301-1315 doi: 10.1051/ro/2017014 [31] Ishibuchi H, Imada R, Masuyama N, Nojima Y. Dynamic specification of a reference point for hypervolume calculation in SMS-EMOA. In: Proceedings of the 2018 IEEE Congress on Evolutionary Computation. Rio de Janeiro, Brazil: IEEE, 2018. 1−8 [32] Brockhoff D, Wagner T, Trautmann H. On the properties of the R2 indicator. In: Proceedings of the 14th Annual Conference on Genetic and Evolutionary Computation. Philadelphia, USA: ACM, 2012. 465−472 [33] Li F, Liu J C, Song Y X, Shang L L. An adaptive evolutionary multi-objective algorithm based on R2 indicator. In: Proceedings of the 2019 Chinese Control and Decision Conference. Nanchang, China: IEEE, 2019. 692−697 [34] Wang H D, Jin Y C, Yao X. Diversity assessment in many-objective optimization. IEEE Transactions on Cybernetics, 2017, 47(6): 1510-1522 doi: 10.1109/TCYB.2016.2550502 [35] Tian Y, Zhang X Y, Cheng R, Jin Y C. A multi-objective evolutionary algorithm based on an enhanced inverted generational distance metric. In: Proceedings of the 2016 IEEE Congress on Evolutionary Computation. Vancouver, Canada: IEEE, 2016. 5222−5229 [36] Tian Y, Cheng R, Zhang X Y, Cheng F, Jin Y C. An indicator-based multiobjective evolutionary algorithm with reference point adaptation for better versatility. IEEE Transactions on Evolutionary Computation, 2018, 22(4): 609-622 doi: 10.1109/TEVC.2017.2749619 [37] Deb K, Agrawal R B. Simulated binary crossover for continuous search space. Complex Systems, 1995, 9(3): 115-148 [38] Deb K, Goyal M. A combined genetic adaptive search (GeneAS) for engineering design. Computer Science and Informatics, 1996, 26(4): 30-45 [39] Das I, Dennis J E. Normal-boundary intersection: A new method for generating the Pareto surface in nonlinear multicriteria optimization problems. Siam Journal on Optimization, 1998, 8(3): 631-657 doi: 10.1137/S1052623496307510 [40] Li M Q, Zhen L L, Yao X. How to read many-objective solution sets in parallel coordinates. IEEE Computational Intelligence Magazine, 2017, 12(4): 88-100 doi: 10.1109/MCI.2017.2742869 期刊类型引用(3)
1. 井荣枝，王延堂，陈小乐，徐峰. 空间约束下多关节机械臂架末端柔顺控制方法. 机械设计与制造. 2024(11): 122-126 . 
百度学术2. 甘雨，郭鹏，林立栋. 基于变分贝叶斯推断的DPGMM风电机组异常数据识别研究. 动力工程学报. 2023(07): 885-892 . 百度学术3. 杨怡婷. 两轮自平衡可移动机器人能耗最优运动轨迹规划方法. 兰州文理学院学报(自然科学版). 2022(01): 59-63 . 百度学术其他类型引用(6)
-
下载:
计量
- 文章访问数: 880
- HTML全文浏览量: 219
- PDF下载量: 167
- 被引次数: 9
