Strategy Self-adaptive Differential Evolution Algorithm Based on State Estimation Feedback
-
摘要: 借鉴闭环控制思想, 提出基于状态估计反馈的策略自适应差分进化(Differential evolution, DE)算法, 通过设计状态评价因子自适应判定种群个体所处于的阶段, 实现变异策略的反馈调节, 达到平衡算法全局探测和局部搜索的目的.首先, 基于抽象凸理论对种群个体建立进化状态估计模型, 提取下界估计信息并结合进化知识设计状态评价因子, 以判定当前种群的进化状态; 其次, 利用状态评价因子的反馈信息, 实现不同进化状态下策略的自适应调整以指导种群进化, 达到提高算法搜索效率的目的.另外, 20个典型测试函数与CEC2013测试集的实验结果表明, 所提算法在计算代价、收敛速度和解的质量方面优于主流改进差分进化算法和非差分进化算法.Abstract: Inspired by the idea of closed-loop control, a strategy self-adaptive differential evolution (DE) algorithm based on state estimation feedback is proposed, the stage of individual can be self-adaptively determined by designing the state judgment factor, and achieve the feedback adjustment of mutation strategies. Consequently, the algorithm can get a trade-off between the exploration and exploitation. Firstly, the estimation model of evolution state is established based on abstract convex theory, from which the underestimation information is extracted combining with the evolutionary information to design the state judgment factor, so that the evolution state of the current population is estimated. Secondly, according to the feedback information of the state judgment factor, the strategy in different evolution state is adaptively selected to guide the evolution of the population. Therefore, the searching efficiency of the algorithm can be improved. Additionally, experimental results of 20 benchmark functions and CEC2013 test set show that the proposed algorithm is superior to the main-stream differential evolution variants and non-differential evolution algorithms mentioned in this paper in terms of computational cost, convergence speed, and solution quality.
-
Key words:
- Differential evolution (DE) /
- state estimation /
- feedback /
- global optimization /
- abstract convex
-
工业过程运行反馈控制包括底层回路关键被控变量的反馈控制和上层运行指标的反馈控制, 也就是说, 工业过程运行反馈控制不仅包括保证过程控制系统关键被控变量的跟踪控制, 而且还要选择合适的关键被控变量设定值, 实现运行指标目标值的跟踪[1].如图 1所示, 传统的工业过程运行反馈控制过程中, 关键被控变量的设定值$ {Y^*} = y_j^*$, $j=1, 2$, $\cdots$, $n$由工序作业班的运行工程师根据运行指标目标值${R^*} = R_i^*$, $j=1, 2, \cdots, m$和多年积累的人工操作经验, 并结合各种运行工况信息人为给出.为实现工业运行过程的自动控制, 自上世纪80年代末以来, 很多学者开展了工业运行过程控制方法的研究.文献[2]基于分层递阶控制的架构和多层优化理论, 提出了反馈优化控制的思想.文献[3]通过离线选择与工业过程经济效益相关的被控变量的设定值, 提出了自优化控制的概念.文献[4]将底层回路控制与过程运行优化相结合, 提出了具有两层结构的实时优化(Real time optimization, RTO)控制方法, 上层采用静态模型优化经济性能指标, 产生底层控制回路的设定值, 通过底层控制器使被控变量跟踪设定值, 从而尽可能使过程运行在经济指标目标值附近.文献[5]将RTO与模型预测控制相结合, 提出了具有三层结构的运行反馈控制方法.此外, 还有一些基于神经网络、模糊推理、案例推理等智能技术的运行反馈控制方法, 例如文献[6]将案例推理、规则推理以及神经网络相结合, 提出了工业运行过程的混合智能控制方法; 文献[7]将神经网络与模糊推理相结合, 提出了一种设定值的混合监控方法.
上述运行反馈控制方法均假设底层过程控制可以跟踪运行控制给出的设定值, 没有考虑底层跟踪设定值的动态误差对整个运行过程优化和控制的影响.为解决这一问题, 文献[8-9]提出了使运行指标实际值与目标值偏差和控制回路输出与设定值跟踪误差的二次性能指标极小化的运行反馈控制方法.文献[10]提出了运行反馈解耦控制方法.上述方法均假设运行层的模型由底层关键被控变量与运行指标之间的静态模型精确描述.实际上, 运行指标反映产品在加工过程中的质量、效率、消耗等, 与底层控制回路的被控变量之间往往具有动态特性, 并且很难用精确的数学模型描述.
本文针对一类运行层为未知动态模型的工业运行过程, 提出一种新的多模型自适应控制方法.最早的多模型自适应控制方法通过线性模型和基于神经网络的非线性模型之间的切换不仅可以保证自适应系统有界输入和有界输出(Bound-input and bound-output, BIBO)稳定, 而且可以改善系统的跟踪性能[11], 但该方法只适用于单输入、单输出系统, 并且是在系统的未建模动态全局有界这一假设下实现的.文献[12]将上述方法推广到多变量系统, 提出了基于多模型与神经网络的多变量自适应控制方法, 放松了文献[11]对系统未建模动态全局有界的假设.文献[13-14]提出了多变量强耦合系统的多模型自适应解耦控制方法.文献[15]提出了参数跳变系统的多模型自适应控制方法.文献[16]提出了具有未知执行器非线性的多变量自适应控制方法.
上述多模型自适应控制方法都是针对底层被控对象设计的.这些方法采用带死区的投影算法对未知参数进行在线辨识.投影算法收敛速度慢, 对参数初值十分灵敏, 实际使用中只有当参数初值接近真值时, 才具有良好的收敛效果, 因此投影算法对过程的先验知识要求较高, 不适合应用于动态未知的工业运行过程.相比较, 递推最小二乘算法具有较快的收敛速度, 对参数初值不灵敏.本文提出的运行过程多模型自适应控制方法采用带死区的递推最小二乘方法对未知参数进行在线辨识.理论分析和仿真实验验证了所提方法的有效性.
1. 问题描述
工业运行过程动态模型由上层运行层的动态模型和底层被控对象的动态模型两部分组成.在本文中, 为了将问题简化, 底层被控对象由如下离散时间线性状态空间模型描述.
$ x(k+1)=\bar{A}x(k)+\bar{B}u(k) $
(1a) $ y(k)=\bar{C}x(k) $
(1b) 其中, $x\in{\bf R}^n$为被控对象状态, $u\in{\bf R}^m$为被控对象的控制输入, $y\in{\bf R}^m$是被控对象的测量输出, $\overline{A}$ $\in$ ${\bf R}^{n\times n}$, $\overline{B}\in{\bf R}^{n\times m}$, $\overline{C}\in{\bf R}^{m\times n}$为时不变矩阵.针对底层被控对象(1)设计极点配置控制器.
$ u(k)=-Kx(k)+L{{y}^{*}}(k) $
(2) 其中, $y^{*}(k)$为底层回路设定值, $K\in{\bf R}^{m\times n}$, $L$ $\in$ ${\bf R}^{m\times m}$为时不变矩阵.
为获得控制器参数矩阵$K$和$L$, 将式(2)代入式(1)得到闭环系统方程为
$ x(k+1)=(\bar{A}-\bar{B}K)x(k)+\bar{B}L{{y}^{*}}(k) $
(3a) $ y(k)=\bar{C}x(k) $
(3b) 为使闭环系统稳定, 并实现稳态跟踪, 应选择控制器参数矩阵和满足:
1) 矩阵$\overline{A}-\overline{B}K$稳定;
2) $\lim\nolimits_{z\rightarrow1}\overline{C}(zI_{n}-(\overline{A}-\overline{B}K))^{-1}\overline{B}L=I_{m}$, $L=$ $\lim\nolimits_{z\rightarrow1}(\overline{C}(zI_{n}- (\overline{A}-\overline{B}K))^{-1}\overline{B})^{-1}$.
由于上层运行层动态模型是底层关键被控变量与运行指标之间的函数, 它的输出与底层控制系统输出相关.在本文中, 考虑运行层模型为如下带有未建模动态的动态模型.
$ r(k+1)=Mr(k)+Ny(k)+\nu (k) $
(4) 其中, $r(k)$为运行过程输出, 即运行过程的工艺指标, $\nu(k)\in {\bf R}^m$为外部干扰或未建模动态, $M$, $N$ $\in$ ${\bf R}^{m\times m}$为时不变矩阵.工业过程运行控制系统涉及到底层关键被控变量的反馈控制和上层运行指标的反馈控制, 为充分考虑底层跟踪设定值的动态误差对整个运行过程控制的影响, 运行过程动态模型可看作是由底层基础反馈控制系统(3)和运行层动态模型(4)构成的广义模型.
$ x(k+1)=\widetilde{A}x(k)+\widetilde{B}{{y}^{*}}(k) $
(5a) $ r(k+1)=Mr(k)+\widetilde{C}x(k)+\nu (k) $
(5b) 其中, $\widetilde{A}=\bar{A}-\bar{B}K$, $\widetilde{B}=\bar{B}L, $ $\widetilde{C}=N\bar{C}, $满足${{\widetilde{C}}^{\text{T}}}\widetilde{C}$可逆.
假设 1. 未建模动态$\nu(k)$的差分项或变化率全局有界, 即对任意的$k > 0$, $\|\nu(k)-\nu(k-2)\|\leq\Gamma$, 其中, $\Gamma$为正常数.
本文的目标是将设定值$y^{*}(k)$看作控制输入, 确定一个多模型自适应控制器, 当其应用于不确定的运行过程(5)时, 闭环运行过程的输入、输出信号有界, 即闭环系统BIBO稳定, 并且运行过程输出$r(k)$尽可能跟踪事先指定的运行指标目标值$r^{*}(k)$的变化.由于未建模动态的存在, 单独使用线性控制器即使能保证闭环运行过程BIBO稳定, 也很难满足一定的跟踪性能.本文将基于带死区的递推最小二乘算法的线性鲁棒自适应控制器和具有未建模动态补偿的非线性自适应控制器与切换机制相结合, 提出的多模型自适应控制器不仅能够保证闭环运行过程BIBO稳定, 而且可使其具有良好的跟踪性能.
2. 基于带死区的递推最小二乘算法的线性鲁棒自适应控制
2.1 一步超前控制器设计
为进行控制器设计, 首先需要将广义模型(5)转化成差分方程形式, 为此引入后移算子$z^{-1}$, 于是式(5)可以重新整理为
$ A({{z}^{-1}})r(k+2)=B{{y}^{*}}(k)+C({{z}^{-1}})\nu (k+1) $
(6) 其中,
$ A({{z}^{-1}})=\widetilde{C}[{{I}_{n}}-\widetilde{A}{{z}^{-1}}]{{({{\widetilde{C}}^{\text{T}}}\widetilde{C})}^{-1}}{{\widetilde{C}}^{\text{T}}}({{I}_{m}}-M{{z}^{-1}}) $
$ B=\widetilde{C}\widetilde{B} $
$ C({z^{ - 1}}) = \widetilde C[{I_n} - \widetilde A{z^{ - 1}}]{({\widetilde C^{\rm{T}}}\widetilde C)^{ - 1}}{\widetilde C^{\rm{T}}} $
下面针对模型(6)设计一步超前控制器.引入如下一步超前最优性能指标:
$ J(k) = {\left\| {T({z^{ - 1}})r(k + 2) - R({z^{ - 1}}){r^*}(k + 2)} \right\|^2} $
(7) 其中, $r^{*}(k)=[r^{*}_{1}(k), r^{*}_{2}(k), \cdots, r^{*}_{m}(k)]^{\rm T}\in{\bf R}^m$为已知有界的运行指标目标值, $T(z^{-1})\in{\bf R}^{m\times m}$为稳定的对角加权多项式矩阵, 满足$T(0)$非奇异; $R(z^{-1})$ $\in$ ${\bf R}^{m\times m}$为对角加权多项式矩阵.引入方程
$ T({z^{ - 1}}) = H({z^{ - 1}})A({z^{ - 1}}) + {z^{ - 2}}G({z^{ - 1}}) $
(8) 为使$H(z^{-1})$和$G(z^{-1})$为唯一解或最小阶解, 由文献[17]可知, $H(z^{-1})$和$G(z^{-1})$都为1阶多项式矩阵, $T(z^{-1})$的阶次小于或等于3.易知, $H(0)=T(0)$.将式(6)两边乘$H(z^{-1})$并利用式(8), 得
$ \begin{array}{l} T({z^{ - 1}})r(k + 2) = G({z^{ - 1}})r(k) + \\ \;\;\;\;\;\;H({z^{ - 1}})B{y^*}(k) + H({z^{ - 1}})C({z^{ - 1}})\nu (k + 1) \end{array} $
(9) 定义时滞-差分算子$\Delta=1-z^{-2}$, 则式(9)转化为
$ \begin{array}{l} T({z^{ - 1}})r(k + 2) = G({z^{ - 1}})\Delta r(k) + \\ \;\;\;\;\;H({z^{ - 1}})B\Delta {y^*}(k) + T({z^{ - 1}})r(k) + \rho (k) \end{array} $
(10) 其中, $\rho(k)=H(z^{-1})C(z^{-1})[\nu(k+1)-\nu(k-1)]$.于是, 使性能指标(7)最小的一步超前最优控制$y^{*}(k)$通过下式计算.
$ \begin{array}{l} G({z^{ - 1}})\Delta r(k) + H({z^{ - 1}})B\Delta {y^*}(k) + \rho (k) = \\ \;\;\;\;\;\;R({z^{ - 1}}){r^*}(k + 2) - T({z^{ - 1}})r(k) \end{array} $
(11) 将式(11)代入模型(6), 得到运行过程闭环方程
$ T({z^{ - 1}})r(k + 2) = R({z^{ - 1}}){r^*}(k + 2) $
(12) 由式(12)可知, 若选择$R(z^{-1})=T(z^{-1})$, 则可消除运行过程的跟踪误差.
由于外部干扰或未建模动态往往是未知的, 当不考虑它对运行过程闭环系统的影响时, 可采用下面的线性控制器方程求取控制输入$y^{*}(k)$.
$ \begin{array}{l} G({z^{ - 1}})\Delta r(k) + H({z^{ - 1}})B\Delta {y^*}(k) = \\ \;\;\;\;\;\;R({z^{ - 1}}){r^*}(k + 2) - T({z^{ - 1}})r(k) \end{array} $
(13) 2.2 线性鲁棒自适应控制
运行过程的动态模型往往是未知的, 因此需要采用自适应方法在线获得控制器参数, 当组成$A(z^{-1})$, $B$, $C(z^{-1})$的参数阵未知时, 式(10)可看作控制器参数辨识方程, 为此记$\phi(k)=T(z^{-1})r(k)$, $G(z^{-1})$ $=G_0+G_1(z^{-1})$, $Q(z^{-1})=H(z^{-1})B=$ $Q_0 +Q_1(z^{-1})$, 并定义数据向量和参数矩阵分别为$\varphi(k)$ $=[\Delta r^{\rm T}(k), \Delta r^{\rm T}(k-1), \Delta {y^{*}}^{\rm T}(k)$, $\Delta {y^{*}}^{\rm T}(k-1)]^{\rm T}$和$\theta=[G_0, G_1, Q_0, Q_1]^{\rm T}$, 则控制器参数辨识方程(10)可以写为
$ \phi (k + 2) = {\theta ^{\rm{T}}}\varphi (k) + \phi (k) + \rho (k) $
(14) 线性控制器方程(13)可重新写为
$ \theta^{\rm T}\varphi(k)=R(z^{-1})r^{*}(k+2)-T(z^{-1})r(k) $
(15) 对于未知的参数矩阵$\theta$, 采用带死区的递推最小二乘方法进行在线辨识.
$ \hat \theta (k) = {\rm{proj}}\{ {\hat \theta ^ + }(k)\} $
(16a) $ \begin{array}{l} {{\hat \theta }^ + }(k) = \hat \theta (k - 2){\mkern 1mu} + \\ \;\;\;\;\;\;\;\;\;\;\;\frac{{\lambda (k)P(k - 2)\varphi (k - 2){e^{\rm{T}}}(k)}}{{1 + {\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)}} \end{array} $
(16b) $ \begin{array}{l} P(k) = P(k - 2) - \\ \;\;\;\;\;\;\;\;\;\;\frac{{\lambda (k)P(k - 2)\varphi (k - 2){\varphi ^{\rm{T}}}(k - 2)P(k - 2)}}{{1 + {\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)}} \end{array} $
(16c) $ e(k) = \phi (k) - \hat \phi (k) $
(16d) $ \hat \phi (k) = {\hat \theta ^{\rm{T}}}(k - 2)\varphi (k - 2) + \phi (k - 2) $
(16e) $ \lambda \left( k \right) = \left\{ \begin{array}{l} \frac{1}{2}, \;\;如果\left\| {e\left( k \right)} \right\|>2E\\ 0, \;\;否则 \end{array} \right. $
(16f) $ {\rm{proj}}\left\{ {{{\hat \theta }^ + }\left( k \right)} \right\} = \left\{ \begin{array}{l} {{\hat \theta }^ + }\left( k \right), \;\;\;\;\;\;\;\;\;\;\;\;\;\hat Q_0^ + \left( k \right)非奇异\\ {\left[ { \ldots , {Q_{\min }}, \ldots } \right]^{\rm{T}}}, \;\;\;否则 \end{array} \right. $
(16g) 其中, $[\varphi(0), \widehat{\theta}(0), P(0)]$为初始条件, $P(0)>0$为正定矩阵, $E$为$\rho(k)$的已知上界, $\widehat{\theta}(k)= [\widehat{G}_0(k)$, $\widehat{G}_1(k)$, $\widehat{Q}_0(k), \widehat{Q}_1(k)]^{\rm T}$为$k$时刻未知参数矩阵$\theta$的估计, $\widehat{\theta}^{+}(k)=[\widehat{G}_0(k), \widehat{G}_1(k), \widehat{Q}_0^{+}(k), \widehat{Q}_1(k)]^{\rm T}$, ${\rm proj}\{\cdot\}$为一投影算子, 满足式(16g).
由式(15)及确定性等价原则可知, 线性鲁棒自适应控制器为
$ {\hat \theta ^{\rm{T}}}(k)\varphi (k) = R({z^{ - 1}}){r^*}(k + 2) - T({z^{ - 1}})r(k) $
(17) 2.3 线性鲁棒自适应控制系统稳定性和性能
引理 1. 定义函数
$ V(k) = {\rm{tr}}\left[ {{{\widetilde \theta }^{\rm{T}}}(k){P^{ - 1}}(K)\widetilde \theta (k)} \right] $
则带死区的递推最小二乘辨识算法(16)具有如下性质:
1)
$ \begin{array}{l} V(k) - V(k - 2) \le \\ \;\;\;\; - \frac{{3\lambda (k){{\left\| {e(k)} \right\|}^2}}}{{8[1 + {\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)]}} - \\ \;\;\;\;\frac{{\lambda (k)[{{\left\| {e(k)} \right\|}^2} - 4{E^2}]}}{{4\{ 1 + [1 - \lambda (k)]{\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)\} }} \end{array} $
2)
$ \mathop {\lim }\limits_{k \to \infty } \left\| {\hat \theta (k) - \hat \theta (k - 2)} \right\| = 0 $
证明. 见附录A.
定理 1. 运行过程动态模型(5)或(6)满足假设1, 则当线性鲁棒自适应控制算法(16)应用于式(6)时, 闭环运行过程全局李雅普诺夫稳定, 并且广义跟踪误差满足${\lim _{k \to \infty }}\lambda (k)[{\left\| {\bar e(k)} \right\|^2} - 4{E^2}] = 0$, 其中, $\bar e(k): = T({z^{ - 1}})r(k) - R({z^{ - 1}}){r^*}(k)$.
证明. 由引理1的1)可知,
$ \mathop {\lim }\limits_{k \to \infty } \frac{{\lambda (k)[{{\left\| {e(k)} \right\|}^2} - 4{E^2}]}}{{4\{ 1 + [1 - \lambda (k)]{\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)\} }} = 0 $
(18) 由于$\overline{e}(k):=T(z^{-1})r(k)-R(z^{-1})r^*{(k)}$及$T(z^{-1})$的稳定性, 存在正常数$c_1$, $c_2$, $c_3$, $c_4$满足
$ \begin{array}{l} |{r_i}(k)| \le {c_1} + {c_2}\mathop {\max }\limits_{_{\scriptstyle0 \le \tau \le t\atop \scriptstyle1 \le i \le m}} |{{\bar e}_i}(\tau )|, \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;i = 1, 2, \cdots , m \end{array} $
(19) $ \begin{array}{l} |y_i^*(k - 2)| \le {c_3} + {c_4}\mathop {\max }\limits_{_{\scriptstyle0 \le \tau \le t\atop \scriptstyle1 \le i \le m}} |{r_i}(\tau )|, \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;i = 1, 2, \cdots , m \end{array} $
(20) 令
$ \begin{array}{l} X(k - 2) = \\ \;\;\;{[{r^{\rm{T}}}(k - 2), {r^{\rm{T}}}(k - 3), {y^*}^{\rm{T}}(k - 2), {y^*}^{\rm{T}}(k - 3)]^{\rm{T}}} \end{array} $
则存在正常数$c_5$, $c_6$满足
$ \left\| {X(k - 2)} \right\| \le {c_5} + {c_6}\mathop {\max }\limits_{0 \le \tau \le t} \left\| {\bar e(\tau )} \right\| $
(21) 由式(16d)可知,
$ e(k)=T(z^{-1})r(k)-R(z^{-1})r^{*}(k)=\overline{e}(k) $
(22) 因此, 由式(21)和式(22)可知
$ \left\| {X(k - 2)} \right\| \le {c_5} + {c_6}\mathop {\max }\limits_{0 \le \tau \le t} \left\| {e(\tau )} \right\| $
(23) 由式(23)可知, 单独采用线性鲁棒自适应控制算法时, 系统输入和输出信号的有界性由$e(k)$的有界性决定.下面假设$e(k)$无界.由式(16f)可知, 存在时刻$K_0>0$, 当$k>K_0$时, $\left\| {\mathit{e}(\mathit{k})} \right\| > {\rm{2}}E$并且$\lambda(k)$ $=$ $1/2$, 即式(18)的分子是一个正实序列.于是存在一单调递增序列$\left\| {\mathit{e}({\mathit{k}_n})} \right\|$, 使得
$ {\lim _{k \to \infty }}\left\| {e({k_n})} \right\| = \infty $
由式(23)可知
$ \begin{array}{l} \frac{{\lambda ({k_n})[{{\left\| {e({k_n})} \right\|}^2} - 4{E^2}]}}{{41 + [1 - \lambda ({k_n})]{\varphi ^{\rm{T}}}({k_n} - 2)P({k_n} - 2)\varphi ({k_n} - 2)}} \ge \\ \frac{{\lambda ({k_n})[{{\left\| {e({k_n})} \right\|}^2} - 4{E^2}]}}{{81 + [1 - \lambda ({k_n})][{{({c_5} + {c_6}\left\| {e({k_n})} \right\|)}^2}]\left\| {P({k_n} - 2)} \right\|}} \end{array} $
由于$\left\| {\mathit{P}({\mathit{k}_n} - {\rm{2}})} \right\|$为递减序列, 因此,
$ \left\| {P({k_n} - 2)} \right\| \le \left\| {P(0)} \right\| $
因此
$ \begin{array}{l} \mathop {\lim }\limits_{k \to \infty } \frac{{\lambda ({k_n})[{{\left\| {e({k_n})} \right\|}^2} - 4{E^2}]}}{{41 + [1 - \lambda ({k_n})]{\varphi ^{\rm{T}}}({k_n} - 2)P({k_n} - 2)\varphi ({k_n} - 2)}} \ge \\ \;\;\;\;\;\frac{1}{{8{c_6}\left\| {P(0)} \right\|}} > 0 \end{array} $
这与式(18)矛盾.故假设不成立, $e(k)$有界, 从而采用线性鲁棒自适应控制算法时, 闭环系统BIBO稳定.
注释 1. 单独使用线性鲁棒自适应控制器能够保证闭环运行过程全局李亚普洛夫稳定, 但是无法使运行过程具有良好的跟踪性能.为了改善运行过程的跟踪性能, 同时不影响其稳定性, 我们将线性鲁棒自适应控制器、基于神经网络的非线性控制器以及切换机制相结合, 提出一种新的多模型自适应控制方法.
3. 多模型自适应控制
3.1 多模型自适应控制算法
下面考虑多模型自适应控制问题, 为将问题简化, 采用两个模型, 多模型切换系统如图 2所示.
图 2中, 线性估计模型定义为
$ \widehat{\phi}_1(k)=\widehat{\theta}_1^{\rm T}(k-2)\varphi(k-2)+\phi(k-2) $
(24) 其中, $\widehat{\theta}_1(k)=[\widehat{G}_{1, 0}(k), \widehat{G}_{1, 1}(k), \widehat{Q}_{1, 0}(k), \widehat{Q}_{1, 1}(k)]^{\rm T}$为$k$时刻的基于线性模型(24)的估计, 采用式(16)在线辨识, 线性辨识误差为
$ {e_1}(k) = \phi (k) - {\hat \phi _1}(k) $
(25) 通过下式计算控制输入$y^*(k)$, 作为基于线性模型的控制器$y^*_1(k)$.
$ \hat \theta _1^{\rm{T}}(k)\varphi (k) = R({z^{ - 1}}){r^*}(k + 2) - T({z^{ - 1}})r(k) $
(26) 非线性估计模型定义为
$ {\hat \phi _2}(k) = \hat \theta _2^{\rm{T}}(k - 2)\varphi (k - 2) + \phi (k - 2) + \hat \rho (k - 2) $
(27) 其中,$\widehat{\theta}_2(k)=[\widehat{G}_{2, 0}(k), \widehat{G}_{2, 1}(k), \widehat{Q}_{2, 0}(k), \widehat{Q}_{2, 1}(k)]^{\rm T}$为$k$时刻$\theta$的基于非线性模型(27)的估计; $\widehat{\rho}(k)$为$\rho^*(k)$的神经网络估计, 其中, $\rho^*(k):=\Delta\phi(k+2)-\widehat{\theta}_2(k)^{\rm T}\varphi(k)$, 即
$ \widehat{\rho}(k)=NN[\widehat{W}(k), \varphi(k)] $
(28) 其中, $NN[\cdot]$表示神经网络结构; $\varphi(k)$为神经网络的输入向量; $\widehat{W}(k)$为$k$时刻理想权阵$W^*$的估计.与文献[12]类似, 除了要求参数阵的估计$\widehat{\theta}_2(k)$和权阵的估计$\widehat{W}(k)$有界, $\widehat{Q}_{2, 0}(k)$非奇异, 并未对$\widehat{\theta}_2(k)$的辨识算法和神经网络的结构以及权阵校正算法进行任何限制, 即
$ \widehat{\theta}_2(k), \widehat{W}(k);~~\widehat{Q}_{2, 0}(k)~\text{非奇异}, ~ \forall k $
(29) 非线性辨识误差为
$ e_2(k)=\phi(k)-\widehat{\phi}_2(k) $
(30) 因此, 根据式(10)和确定性等价原则, 通过下式计算控制输入$y^*(k)$, 作为基于非线性模型的控制器$y^*_2(k)$.
$ \begin{array}{l} {{\hat \theta }_2}(k)\varphi (k) + \hat \rho (k) = R({z^{ - 1}}){r^*}(k + 2){\mkern 1mu} - \\ \;\;\;\;\;\;\;\;T({z^{ - 1}})r(k) \end{array} $
(31) 切换准则为
$ \begin{array}{l} {J_j}\left( k \right) = \sum\limits_{i = 2}^k {\frac{{{\lambda _j}(k)[{{\left\| {{e_j}(k)} \right\|}^2} - 4{E^2}]}}{{4\{ 1 + [1 - {\lambda _j}(k)]{\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)\} }} + } \\ {c_0}\sum\limits_{l = k - N - 1}^k {\left( {\frac{1}{2} - {\lambda _j}(l){{\left\| {{e_j}(l)} \right\|}^2}} \right)} \end{array} $
(32) $ {\lambda _j}\left( k \right) = \left\{ \begin{array}{l} \frac{1}{2}, \;\;若\left\| {{e_j}(k)} \right\|>2E\\ 0, \;\;否则 \end{array} \right. $
(33) 其中, $N$是一个正整数, $c_0$是一个大于等于0的预先确定的常数.
每一时刻$k$, 比较$J_1(k)$和$J_2(k)$, 求出最小的$J^*(k)$, 选择与$J^*(k)$对应的自适应控制器$y_i^*(k)$, 并将其应用于运行过程.
3.2 多模型自适应控制系统稳定性和性能
定理 2. 运行过程动态模型(6)满足假设1, 则当基于多模型自适应控制算法(24)~ (33)用于运行过程(6)时, 闭环切换系统BIBO稳定.此外, 对于任意给定的正数$\varepsilon$, 存在时刻$K$, 当$k>K$时, 系统的广义跟踪误差满足$\left\| {\bar e(k)} \right\| \le 2E + \varepsilon $.
证明. 由引理1可知,
$ \mathop {\lim }\limits_{k \to \infty } \frac{{{\lambda _1}(k)[{{\left\| {{e_1}(k)} \right\|}^2} - 4{E^2}]}}{{4\{ 1 + [1 - {\lambda _1}(k)]{\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)\} }} = 0 $
(34) 由式(24)和式(25)可知,
$ \begin{array}{l} {e_1}(k) = \phi (k) - {{\hat \phi }_1}(k) = \\ \;\;\;\;\;\;\;\;\;\;\Delta \phi (k) - \hat \theta _1^{\rm{T}}(k - 2)\varphi (k - 2) = \\ \;\;\;\;\;\;\;\;\;T({z^{ - 1}})r(k) - T({z^{ - 1}})r(k - 2) - \\ \;\;\;\;\;\;\;\;\hat \theta _1^{\rm{T}}(k - 2)\varphi (k - 2) \end{array} $
(35) 由式(27)和式(30)可知
$ \begin{array}{l} {e_2}(k) = \phi (k) - {{\hat \phi }_2}(k) = \\ \;\;\;\;\;\;\Delta \phi (k) - \hat \theta _2^{\rm{T}}(k - 2)\varphi (k - 2) - \hat \rho (k - 2) = \\ \;\;\;\;\;\;T({z^{ - 1}})r(k) - T({z^{ - 1}})r(k - 2) - \\ \;\;\;\;\;\;\hat \theta _2^{\rm{T}}(k - 2)\varphi (k - 2) - \hat \rho (k - 2) \end{array} $
(36) 于是, 根据确定性等价原则, 每一时刻
$ \bar{e}(k)=e_1(k)~ \mbox{或}~e_2(k) $
(37) 由于每一时刻, 系统辨识误差$e(k)=e_1(k)$或$e_2(k)$, 故由式(21)可知, 存在正常数$c_7$, $c_8$满足
$ \left\| {X(k - 2)} \right\| \le {c_7} + {c_8}\mathop {\max }\limits_{0 \le \tau \le k} \left\| {e(\tau )} \right\| $
(38) 由式(33)可知, 切换函数$J_j(k)$ $(j=1, 2)$的第2项是有界的.因此由引理1可知, $J_1(k)$有界.对于$J_2(k)$有两种情况.
1) $J_2(k)$无界.由于$J_1(k)$有界, 因此存在时刻$K_1$使得当$k\geq K_1$时有$J_1(k)\leq J_2(k)$.故根据切换机制, 当$k\geq K_1+1$时, 系统辨识误差$e(k)=e_1(k)$满足
$ \mathop {\lim }\limits_{k \to \infty } \frac{{{\lambda _1}(k)[{{\left\| {e(k)} \right\|}^2} - 4{E^2}]}}{{4\{ 1 + [1 - {\lambda _1}(k)]{\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)\} }} = 0 $
(39) 其中,
$ \lambda \left( k \right) = \left\{ \begin{array}{l} \frac{1}{2}, \;\;若\left\| {e(k)} \right\|>2E\\ 0, \;\;否则 \end{array} \right. $
2) $J_2(k)$有界.由切换准则式(32)可知, $e_2(k)$满足
$ \mathop {\lim }\limits_{k \to \infty } \frac{{{\lambda _2}(k)[{{\left\| {{e_2}(k)} \right\|}^2} - 4{E^2}]}}{{4\{ 1 + [1 - {\lambda _2}(k)]{\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)\} }} \to 0 $
故系统辨识误差$e(k)=e_1(k)$或$e_2(k)$满足式(39).
其余部分的证明类似于定理1.
由式(39)和$X(k-2)$的有界性可知,
$ \mathop {\lim }\limits_{k \to \infty } \lambda (k)\left[ {{{\left\| {e(k)} \right\|}^2} - 4{E^2}} \right] $
即, 对任意小的正数$\varepsilon$, 存在时刻$K$, 当$k>K$时,
$ \left\| {e(k)} \right\| \le 2E + \varepsilon $
(40) 注释2. 由式(39)可知, 非线性辨识误差
$ {e_2}(k) = {\rho ^*}(k - 2) - \hat \rho (k - 2) $
(41) 适当选择神经网络结构和参数, 可以保证$\left\| {{\rho ^*}(k} \right.$ $-$ $\left. {2) - \hat \rho (k - 2)} \right\|<\varepsilon $.因此若运行过程选择非线性自适应控制器$y_2^*(k)$作为输入信号, 则由式(35)和式(36)可知, 广义跟踪误差$\left\| {\bar e(k)} \right\|<\varepsilon $满足.
4. 仿真实验
为验证本文所提方法的有效性, 首先考虑如下底层被控对象模型
$ \begin{array}{l} x(k + 1) = \left( {\begin{array}{*{20}{c}} {1.5}&6\\ 6&4 \end{array}} \right)x(k) + \\ \;\;\;\;\;\;\;\;\;\left( {\begin{array}{*{20}{c}} { - 4.2623}&{ - 3.8254}\\ {8.3534}&{6.1711} \end{array}} \right)u(k)\\ y(k) = \left( {\begin{array}{*{20}{c}} {0.1546}&{ - 0.012}\\ { - 0.0099}&{0.2281} \end{array}} \right)x(k) \end{array} $
(42) 其中, $x=[x_1, x_2]^{\rm T}$, $y=[y_1, y_2]^{\rm T}$, $u=[u_1, u_2]^{\rm T}$.为使底层闭环系统稳定, 并实现稳态跟踪, 选择如下极点配置控制器
$ u(t)=-Kx(t)+Ly^*(t) $
(43) 其中, $y^*(t)=[y_1^*, y_2^*]^{\rm T}$为底层回路设定值, 由后面的运行控制给出.
$ \begin{array}{l} K = \left( {\begin{array}{*{20}{c}} {7.1487}&{15.3085}\\ { - 8.8017}&{ - 19.0044} \end{array}} \right)\\ \;\;L = \left( {\begin{array}{*{20}{c}} { - 14.5}&{30.25}\\ { - 20.6}&{ - 35.6} \end{array}} \right) \end{array} $
(44) 假设运行层动态模型为
$ \begin{array}{l} r(k + 1) = \left( {\begin{array}{*{20}{c}} 1&0\\ 0&1 \end{array}} \right)r(k) + \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\left( {\begin{array}{*{20}{c}} { - 1.2893}&{ - 0.0678}\\ {0.2798}&{ - 4.3693} \end{array}} \right)y(k) + \nu (k) \end{array} $
(45) 其中,
$ r(k) = {[{r_1}, {r_2}]^{\rm{T}}} $
$ \begin{array}{l} \nu (k) = \left( {\begin{array}{*{20}{c}} {{\nu _1}}\\ {{\nu _2}} \end{array}} \right) = 0.01 \times \\ \;\;\;\;\;\;\;\;\;\;\left( {\begin{array}{*{20}{c}} {{\rm{sin}}({\rm{1 + }}y_{\rm{1}}^{{\rm{*2}}}(k{\rm{ - 1}}){\rm{ + }}r_{\rm{1}}^{\rm{2}}(k{\rm{ - 1}}){\rm{ + }}}\\ {r_2^2(k) - \frac{{{r_1}(k - 1) + {r_2}(k)}}{{1 + y_1^{*2}(k - 1) + r_1^2(k - 1) + r_2^2(k)}})}\\ {{\rm{sin}}({\rm{1 + }}y_{\rm{2}}^{{\rm{*2}}}(k{\rm{ - 1}}){\rm{ + }}r_{\rm{1}}^{\rm{2}}(k{\rm{ - 1}}){\mkern 1mu} {\rm{ + }}}\\ {r_2^2(k) - \frac{{{r_1}(k) + {r_2}(k - 1)}}{{1 + y_2^{*2}(k - 1) + r_1^2(k) + r_2^2(k - 1)}})} \end{array}} \right) \end{array} $
则由式(5)可知, 运行过程广义对象模型为
$ \begin{array}{l} x(k + 1) = \left( {\begin{array}{*{20}{c}} { - 1.7}&{ - 1.45}\\ {0.6}&{ - 6.6} \end{array}} \right)x(k) + \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\left( {\begin{array}{*{20}{c}} {17}&{7.25}\\ { - 6}&{33} \end{array}} \right){y^*}(k)\\ r(k + 1) = \left( {\begin{array}{*{20}{c}} 1&0\\ 0&1 \end{array}} \right)r(k) + \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\left( {\begin{array}{*{20}{c}} { - 0.2}&0\\ 0&{ - 1} \end{array}} \right)x(k) + \nu (k) \end{array} $
(46) 选择加权多项式矩阵
$ \begin{array}{l} T({z^{ - 1}}) = R({z^{ - 1}}) = \\ \;\;\;\;\;\;\;\;\;\;\;\;\left( \begin{array}{l} - 5 - 0.1{z^{ - 1}}\;\;\;\;\;\;\;\;\;0\\ \;\;\;\;\;\;0\;\;\;\;\;\;\;\;\; - 1 - 0.1{z^{ - 1}} \end{array} \right) \end{array} $
运行指标目标值为
$ {r^*}(k) = \left( \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;0.5\\ 0.5{\rm{sign}}\left( {\cos \left( {k \times \frac{\pi }{{50}}} \right)} \right) \end{array} \right) $
易知, 控制器真实参数阵为
$ \theta = \left( {{\theta _0}\;\;\;\;\;{\theta _1}} \right) $
其中,
$ \begin{array}{l} {\theta _0} = ( - 6.53, - 10.556, 1.43, 10.556, 17, 7.25, \\ \;\;\;\;\;\;\;\; - 2.86, - 52.78{)^{\rm{T}}} \end{array} $
$ \begin{array}{l} {\theta _1} = (21.6, - 36.53, - 21.6, 35.43, - 6, 33, \\ \;\;\;\;\;\;\;43.2, - 177.15{)^{\rm{T}}} \end{array} $
本仿真实验中, 我们假设它是未知的, 并根据先验知识选择待辨识控制器初始参数阵为
$ \begin{array}{l} \\ \begin{array}{*{20}{l}} {\hat \theta (0) = {{\left( {\begin{array}{*{20}{c}} { - 4}&{ - 7}&0&4&{19}&3&{ - 1}&{ - 35}\\ {16}&{ - 31}&{ - 15}&{30}&{ - 2}&{29}&{35}&{ - 160} \end{array}} \right)}^{\rm{T}}}} \end{array} \end{array} $
选择单隐层线性输出的静态BP网对$\rho^*(k)$进行估计, 其隐元数为20, 学习率为0.1;选择$c=1$, $N$ $=$ $2$.
图 3为单独采用线性鲁棒自适应控制方法时运行过程的运行指标目标值$r^*(k)$和运行过程输出$r(k)$.从图 3可以看出, 虽然该控制器可以使运行过程稳定, 但跟踪效果很差. 图 4为采用本文所提的多模型自适应控制方法时运行过程的运行指标目标值$r^*(k)$、运行过程输出$r(k)$和运行过程控制输入, 即底层设定值$y^*(k)$.与图 3相比, 图 4中的跟踪效果明显改善. 图 5为$\widehat{\theta}_1(k)$中16个参数的在线变化曲线. 图 6为底层极点配置控制系统的跟踪曲线.为进行比较, 仍以上述矩阵为控制器初始参数阵, 采用文献[12]提出的基于投影算法的多模型自适应控制方法对运行过程进行仿真, 运行过程跟踪结果如图 7所示.相应的$\widehat{\theta}_1(k)$中各参数的在线变化曲线如图 8所示.由图 4和图 7可知, 采用本文提出的基于最小二乘算法的多模型自适应控制方法时, 即使控制器初始参数阵离控制器真实参数阵较远, 仍具有有良好的跟踪效果.两相比较, 基于投影算法的多模型自适应控制方法对初始参数阵非常灵敏, 当初始参数阵远离控制器真实参数阵时, 控制效果较差.比较图 5和图 8可以看出, 最小二乘算法与投影算法相比具有更快的收敛速度.
5. 结论
工业运行过程应考虑底层跟踪设定值的动态误差对整个运行过程优化和控制的影响.现有的工业运行控制方法假设运行层为已知的线性静态模型.本文针对一类运行层为未知线性动态模型的工业运行过程, 提出了一种基于递推最小二乘算法的多模型自适应控制方法.通过理论分析和仿真比较, 验证了与现有的多模型自适应控制方法相比, 本文提出方法可以应用于工业运行过程, 并具有良好的跟踪效果.
工业过程运行控制是近年来控制领域比较热门的研究方向之一, 现有的方法针对的被控对象主要是线性的, 并且主要集中在理论研究上.在实际的工业过程控制中, 非线性动态无可避免, 当两层结构中的被控对象都是非线性时, 如何设计控制器, 如何将理论的研究成果进行实际应用具有一定的挑战.
附录A. 引理1的证明
证明. 当$\widehat{Q}_0^+(k)$非奇异时, $\widehat{\theta}(k)$与$\widehat{\theta}^+(k)$相等.由式(16d)和式(16e)可知,
$ e(k) = [{\theta ^{\rm{T}}} - {\hat \theta ^{\rm{T}}}(k - 1)]\varphi (k - 2) + \rho (k - 2) $
(A1) 令
$ L(k)=\frac{P(k-2)\varphi(k-2)}{1+\varphi^{\rm T}(k-2)P(k-2)\varphi(k-2)} $
(A2) 由式(16b)、式(16c)以及式(A1)和式(A2)可知,
$ \begin{array}{l} P(k) = P(k - 2) - \lambda (k)L(k){\varphi ^{\rm{T}}}(k - 2)P(k - 2) \Rightarrow \\ P(k) = [I - \lambda (k)L(k){\varphi ^{\rm{T}}}(k - 2)]P(k - 2) \Rightarrow \\ P(k){P^{ - 1}}(k - 2) = I - \lambda (k)L(k){\varphi ^{\rm{T}}}(k - 2) \end{array} $
(A3) 和
$ \begin{array}{l} \widetilde \theta (k) = \widetilde \theta (k - 2) + \lambda (k)L(k){e^{\rm{T}}}(k) = \\ \;\;\;\;\;\;\;\;\;\;[I - \lambda (k)L(k){\varphi ^{\rm{T}}}(k - 2)]\widetilde \theta (k - 2) + \\ \;\;\;\;\;\;\;\;\;\;\lambda (k)L(k){\rho ^{\rm{T}}}(k - 2) = \\ \;\;\;\;\;\;\;\;\;\;P(k){P^{ - 1}}(k - 2)\widetilde \theta (k - 2) + \\ \;\;\;\;\;\;\;\;\;\lambda (k)L(k){\rho ^{\rm{T}}}(k - 2) \end{array} $
(A4) 由式(A4)可知,
$ \begin{array}{l} {P^{ - 1}}(k)\widetilde \theta (k) - {P^{ - 1}}(k - 2)\widetilde \theta (k - 2) = \\ \;\;\;\;\;\;\lambda (k){P^{ - 1}}(k)L(k){\rho ^{\rm{T}}}(k - 2) \end{array} $
(A5) 由于$\varphi^{\rm T}(k-2)P(k-2)\varphi(k-2)\times I=\varphi(k-2)\varphi^{\rm T}(k-2)P(k-2)$, 因此由式(16c)可知
$ \begin{array}{l} \frac{{{P^{ - 1}}(k)P(k - 2)\varphi (k - 2)}}{{1 + {\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)}} = \\ \;\;\;\;\frac{{\varphi (k - 2)}}{{1 + [1 - \lambda (k)]{\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)}} \end{array} $
(A6) 令$Q(k):=\varphi^{\rm T}(k-2)P(k-2)\varphi(k-2)$, 则
$ \begin{array}{l} V(k) - V(k - 2) = \\ \;\;\;\;\;\;\;{\rm{tr}}\{ - \frac{{\lambda (k)[e(k){e^{\rm{T}}}(k) - 4\rho (k - 2){\rho ^{\rm{T}}}(k - 2)]}}{{4[1 + [1 - \lambda (k)]Q(k)]}} - \\ \;\;\;\;\;\;\;\left. {\frac{{3\lambda (k)e(k){e^{\rm{T}}}(k)[1 + Q(k)[1 - \frac{{4\lambda (k)}}{3}]]}}{{4[1 + [1 - \lambda (k)]Q(k)][1 + Q(k)]}}} \right\} \end{array} $
(A7) 由于
$ \frac{1+Q(k)\left[1-\frac{4\lambda(k)}{3}\right]}{1+[1-\lambda(k)]Q(k)}\geq \frac{1}{2} $
(A8) 因此
$ \begin{array}{l} V(k) - V(k - 2) \le \\ \;\;\;\;\;{\rm{tr}}\{ - \frac{{\lambda (k)[e(k){e^{\rm{T}}}(k){\rm{ - 4}}\rho (k{\rm{ - 2}}){\rho ^{\rm{T}}}(k{\rm{ - 2}})]}}{{{\rm{4}}[{\rm{1 + }}[{\rm{1 - }}\lambda (k)]Q(k)]}} - \\ \;\;\;\;\;\frac{{3\lambda (k)e(k){e^{\rm{T}}}(k)}}{{8[1 + Q(k)]}}\} \le \\ \;\;\;\;\; - \frac{{\lambda (k)[{{\left\| {e(k)} \right\|}^2} - 4{E^2}]}}{{4[1 + [1 - \lambda (k)]Q(k)]}} - \frac{{3\lambda (k){{\left\| {e(k)} \right\|}^2}}}{{8[1 + Q(k)]}} \end{array} $
(A9) 因此, 引理1中1)得证.由式(16b)可知,
$ \begin{array}{l} {\varphi ^{\rm{T}}}(k - 2)[\hat \theta (k) - \hat \theta (k - 2)] = \\ \;\;\;\;\;\;\frac{{\lambda (k){\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2){e^{\rm{T}}}(k)}}{{1 + {\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)}} \end{array} $
(A10) 因此
$ \begin{array}{l} {\varphi ^{\rm{T}}}(k - 2)[\hat \theta (k) - \hat \theta (k - 2)][\hat \theta (k) - \\ \;\;\;\;\;\;\;\hat \theta (k - 2){]^{\rm{T}}}\varphi (k - 2) \le \\ \;\;\;\;\;\;\;\frac{{\lambda (k){\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2){e^{\rm{T}}}(k)e(k)}}{{1 + {\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)}} \end{array} $
(A11) 故
$ {\left\| {\hat \theta (k) - \hat \theta (k - 2)} \right\|^2} \le \frac{{\lambda (k)\left\| {P(k - 2)} \right\|{e^{\rm{T}}}(k)e(k)}}{{1 + {\varphi ^{\rm{T}}}(k - 2)P(k - 2)\varphi (k - 2)}} $
(A12) 由引理1中的1)及$\| P(k-2)\|$的有界性, 可知2)成立.当$\widehat{Q}_0^+(k)$非奇异时, 记.由式(16g)可知${\left\| {{\theta ^ + }(k) - \theta } \right\|^2}$, 因此$V^+(k)\leq V(k)$, 故引理1\linebreak仍旧成立.
-
表 1 标准测试函数的参数
Table 1 The parameters of benchmark functions
函数名 表达式 维数($N$) 取值范围 全局最小值 Sphere $f_1({\bm x})=\sum\limits_{i=1}^{N}x_i^2$ 30, 50 $(-100, 100)^N$ 0 SumSquares $f_2({\bm x})=\sum\limits_{i=1}^{N}ix_i^2$ 30, 50 $(-10, 10)^N$ 0 Schwefel 2.22 $f_3({\bm x})=\sum\limits_{i=1}^{N}\left|x_i\right|+\prod\limits_{i=1}^{N}\left|x_i\right|$ 30, 50 $(-10, 10)^N$ 0 Exponential $f_4({\bm x})=-\exp(-0.5$$\sum\limits_{i=1}^{N}$$x_i^2)$ 30, 50 $(-1, 1)^N$ $-1$ Tablet $f_5({\bm x})=10^6x_1^2+\sum\limits_{i=2}^{N}x_i^2$ 30, 50 $(-100, 100)^N$ 0 Step $f_6({\bm x})=\sum\limits_{i=1}^{N}\left\lfloor x_i+0.5\right\rfloor^2$ 30, 50 $(-100, 100)^N$ 0 Zakharov $f_7({\bm x})=\sum\limits_{i=1}^{N}x_i^2+\sum\limits_{i=1}^{N}(0.5ix_i)^2+\sum\limits_{i=1}^{N}(0.5ix_i)^4$ 30, 50 $(-5, 10)^N$ 0 Rosenbrock $f_8({\bm x})=\sum\limits_{i=1}^{N-1}[100(x_{i+1}-x_i^2)^2+(x_i-1)^2]$ 30, 50 $(-30, 30)^N$ 0 Griewank $f_9({\bm x})=1+\frac{1}{4000}\sum\limits_{i=1}^{N}x_i^2-\prod\limits_{i=1}^{N}\cos(\frac {x_i}{\sqrt{i}})$ 30, 50 $(-600, 600)^N$ 0 Schaffer 2 $f_{10}({\bm x})=\sum\limits_{i=1}^{N-1}(x_i^2+x_{i+1}^2)^{0.25}(\sin^2(50(x_i^2+x_{i+1}^2)^{0.1})+1)$ 30, 50 $(-100, 100)^N$ 0 Schwefel 2.26 $f_{11}({\bm x})=-\sum\limits_{i=1}^{N}x_i\sin(\sqrt{\left\|x_i\right\|})$ 30, 50 $(-500, 500)^N$ $-12 569.18$ Himmelblau $f_{12}({\bm x})=N^{-1}\sum\limits_{i=1}^{N}(x_i^4-16x_i^2+5x_i)$ 30, 50 $(-5, 5)^N$ -78.3323 Levy and Montalvo 1 $f_{13}({\bm x})=\frac{\pi}{N}(10\sin^2(\pi y_1)+\sum\limits_{i=1}^{N-1}(y_i-1)^{2} [1+\\ 10\sin^2(\pi y_{i+1})]+(y_N-1)^2)$, $y_i=1+\frac{1}{4}(x_i+1)$ 30, 50 $(-10, 10)^N$ 0 Levy and Montalvo 2 $f_{14}({\bm x})=0.1(\sin^2(3\pi x_i)+\sum\limits_{i=1}^{N-1}(x_i-1)^{2} [1+\sin^2(3\pi x_{i+1})]+\\ (x_N-1)^2[1+\sin^2(2\pi x_N)])$ 30, 50 $(-5, 5)^N$ 0 Ackley $f_{15}({\bm x})=-20\exp(-0.02\sqrt{N^{-1}\sum\limits_{i=1}^{N}x_i^2}) -\exp(N^{-1}$ 30, 50 $(-30, 30)^N$ 0 $\sum\limits_{i=1}^{N}\cos(2\pi x_i))+20+e$ Rastrigin $f_{16}({\bm x})=10N+\sum\limits_{i=1}^{N}(x_i^2-10\cos(2\pi x_i))$ 30, 50 $(-5, 5)^N$ 0 Penalized 1 $f_{17}(x)=\frac{\pi}{N}\{\sum\nolimits_{i=1}^{N-1}(y_i-1)^2[1+ \sin(\pi y_{i+1})]+(y_N-1)^2+$ $(10\sin^2(\pi y_1))\}+\sum\nolimits_{i=1}^Nu(x_i, 10, 100, 4)$, $y_i=1+\frac{x_i+1}{4}$ $u(x_i, a, k, m) = \left\{ \begin{array}{ll} k(x_i-a)^m, ~~~~ x_i> a \\ 0, ~~~~-a\leq x_i\leq a \\ k(-x_i-a)^m, ~~~~x_i < -a \end{array} \right.$ 30, 50 $(-50, 50)^N$ 0 Penalized 2 $f_{18}(x)=0.1\{\sin^2(3\pi x_1)+\sum\nolimits_{i=1}^{N-1} (x_i-1)^2[1+\sin^2(3\pi x_{i+1})]+$ $(x_N-1)^2[1+\sin^2(2\pi x_N)]\}+\sum\nolimits_{i=1}^Nu(x_i, 5, 100, 4)$ 30, 50 $-(50, 50)^N$ 0 Neumaier $f_{19}({\bm x})=\sum\limits_{i=1}^{N}(x_i-1)^2- \sum\limits_{i=2}^{N}x_ix_{i-1}+\frac{N(N+4)(N-1)}{6}$ 30, 50 $(-900, 900)^N$ 0 Alpine $f_{20}({\bm x})=\sum\limits_{i=1}^{N}\left|x_i\sin x_i+0.1x_i\right|$ 30, 50 $(-10, 10)^N$ 0 表 2 SEFDE中参数$K$设置下的平均函数评价次数和成功率
Table 2 Average numbers of function evaluations and success rates of parameter settings $K$ in SEFDE
Fun $N$ $K=2$ $K=3$ $K=4$ $K=5$ $K=6$ FEs SR FEs SR FEs SR FEs SR FEs SR $f_{1}$ 30 1.51E+04 1.00 1.19E+04 1.00 1.26E+04 1.00 1.21E+04 1.00 1.18E+04 1.00 $f_{2}$ 30 1.17E+04 1.00 1.14E+04 1.00 1.19E+04 1.00 1.06E+04 1.00 1.14E+04 1.00 $f_{3}$ 30 1.85E+04 1.00 1.76E+04 1.00 1.73E+04 1.00 1.75E+04 1.00 1.83E+04 1.00 $f_{4}$ 30 7.88E+03 1.00 7.89E+03 1.00 7.85E+03 1.00 7.83E+03 1.00 7.87E+03 1.00 $f_{5}$ 30 2.94E+04 1.00 2.97E+04 1.00 2.93E+04 1.00 2.91E+04 1.00 2.90E+04 1.00 $f_{6}$ 30 8.94E+03 1.00 8.61E+03 1.00 8.56E+03 1.00 8.28E+03 1.00 8.45E+03 1.00 $f_{7}$ 30 1.46E+05 1.00 1.45E+05 1.00 1.46E+05 1.00 1.46E+05 1.00 1.44E+05 1.00 $f_{8}$ 30 1.24E+05 1.00 1.20E+05 1.00 1.21E+05 1.00 1.22E+05 1.00 1.21E+05 1.00 $f_{9}$ 30 1.28E+04 1.00 1.26E+04 1.00 1.25E+04 1.00 1.27E+04 1.00 1.24E+04 1.00 $f_{10}$ 30 1.00E+05 1.00 1.09E+05 1.00 1.10E+05 1.00 1.07E+05 1.00 1.09E+05 1.00 $f_{11}$ 30 4.43E+04 1.00 4.34E+04 1.00 4.57E+04 1.00 4.40E+04 1.00 4.48E+04 1.00 $f_{12}$ 30 1.24E+04 1.00 1.71E+04 1.00 1.74E+04 1.00 1.64E+04 1.00 1.70E+04 1.00 $f_{13}$ 30 1.26E+04 1.00 1.24E+04 1.00 1.26E+04 1.00 1.27E+04 1.00 1.26E+04 1.00 $f_{14}$ 30 1.18E+04 1.00 1.19E+04 1.00 1.20E+04 1.00 1.22E+04 1.00 1.19E+04 1.00 $f_{15}$ 30 2.02E+04 1.00 1.95E+04 1.00 2.01E+04 1.00 2.02E+04 1.00 1.98E+04 1.00 $f_{16}$ 30 5.60E+04 1.00 5.44E+04 1.00 5.62E+04 1.00 5.34E+04 1.00 5.66E+04 1.00 $f_{17}$ 30 1.90E+04 1.00 1.96E+04 1.00 1.91E+04 1.00 1.92E+04 1.00 1.98E+04 1.00 $f_{18}$ 30 2.32E+04 1.00 2.29E+04 1.00 2.32E+04 1.00 2.30E+04 1.00 2.29E+04 1.00 $f_{19}$ 30 1.72E+05 1.00 1.70E+05 1.00 1.67E+05 1.00 1.70E+05 1.00 1.68E+05 1.00 $f_{20}$ 30 3.62E+04 1.00 3.75E+04 1.00 3.13E+04 1.00 3.48E+04 1.00 3.47E+04 1.00 AVE 4.41E+04 1.000 4.42E+04 1.000 4.41E+04 1.000 4.40E+04 1.000 4.41E+04 1.000 表 3 函数评价次数和成功率对比数据
Table 3 Compared data on function evaluations and success rates
Fun $N$ DELU SHADE EPSDE CoDE SEFDE FEs SR FEs SR FEs SR FEs SR FEs SR $f_{1}$ 30 1.22E+04 1.00 2.35E+04 1.00 1.67E+04 1.00 4.02E+04 1.00 ${\bf1.21E+04}$ 1.00 $f_{2}$ 30 1.07E+04 1.00 2.16E+04 1.00 1.51E+04 1.00 3.63E+04 1.00 ${\bf1.06E+04}$ 1.00 $f_{3}$ 30 3.15E+04 1.00 3.42E+04 1.00 2.32E+04 1.00 5.64E+04 1.00 ${\bf1.75E+04}$ 1.00 $f_{4}$ 30 7.98E+03 1.00 1.67E+04 1.00 9.36E+03 1.00 2.24E+04 1.00 ${\bf7.83E+03}$ 1.00 $f_{5}$ 30 2.55E+04 1.00 2.67E+04 1.00 ${\bf1.83E+04}$ 1.00 4.13E+04 1.00 2.91E+04 1.00 $f_{6}$ 30 1.16E+04 1.00 1.19E+04 1.00 ${\bf8.33E+03}$ 1.00 2.01E+04 1.00 ${\bf8.28E+03}$ 1.00 $f_{7}$ 30 7.71E+04 1.00 ${\bf5.98E+04}$ 1.00 1.33E+05 1.00 7.80E+04 1.00 1.46E+05 1.00 $f_{8}$ 30 ${\bf7.09E+04}$ 1.00 1.09E+05 1.00 1.16E+05 0.87 2.37E+05 1.00 1.22E+05 1.00 $f_{9}$ 30 1.54E+04 1.00 2.56E+04 1.00 1.76E+04 0.83 4.39E+04 1.00 ${\bf1.27E+04}$ 1.00 $f_{10}$ 30 ${\bf9.78E+04}$ 1.00 1.08E+05 0.97 1.13E+05 1.00 1.73E+05 1.00 1.07E+05 1.00 $f_{11}$ 30 ${\bf1.34E+04}$ 1.00 9.67E+04 1.00 3.71E+04 1.00 6.00E+04 1.00 4.40E+04 1.00 $f_{12}$ 30 1.98E+04 1.00 3.35E+04 0.97 2.33E+04 1.00 3.56E+04 1.00 ${\bf1.64E+04}$ 1.00 $f_{13}$ 30 ${\bf1.05E+04}$ 1.00 1.67E+04 1.00 1.30E+04 1.00 2.63E+04 1.00 1.27E+04 1.00 $f_{14}$ 30 ${\bf8.71E+03}$ 1.00 1.71E+04 1.00 1.15E+04 1.00 2.70E+04 1.00 1.22E+04 1.00 $f_{15}$ 30 2.65E+04 1.00 3.21E+04 1.00 2.32E+04 1.00 5.68E+04 1.00 ${\bf2.02E+04}$ 1.00 $f_{16}$ 30 5.94E+04 1.00 1.24E+05 1.00 6.40E+04 1.00 1.30E+05 0.97 ${\bf5.34E+04}$ 1.00 $f_{17}$ 30 ${\bf9.17E+03}$ 1.00 1.92E+04 1.00 1.29E+04 1.00 3.03E+04 1.00 1.92E+04 1.00 $f_{18}$ 30 1.52E+04 1.00 ${\bf1.50E+04}$ 1.00 1.59E+04 0.97 3.54E+04 1.00 2.30E+04 1.00 $f_{19}$ 30 1.68E+05 0.97 ${\bf1.09E+05}$ 1.00 NA 0.00 2.69E+05 1.00 1.70E+05 1.00 $f_{20}$ 30 NA 0.00 NA 0.00 1.00E+05 1.00 NA 0.00 ${\bf3.48E+04}$ 1.00 AVE 4.96E+04 0.95 6.00E+04 0.947 5.36E+04 0.934 8.59E+04 0.949 $\textbf{4.40E+04}$ $\textbf{1.000}$ 表 4 DELU、SHADE、EPSDE、CoDE、SEFDE 对30维测试函数的优化结果, 平均值(标准差)
Table 4 Compared data of 30 D optimization results on DELU、SHADE、EPSDE、CoDE、SEFDE, Mean (Std)
Fun $N$ DELU SHADE EPSDE CoDE SEFDE Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) $f_{1}$ 30 9.25E-19(2.99E-36)$^+$ 2.87E-19(2.31E-19)$^+$ 3.87E-31(9.46E-31)$^+$ 2.16E-10(1.73E-10)$^+$ ${\bf4.01E-38}({\bf5.05E-38})$ $f_{2}$ 30 2.56E-32(4.00E-32)$^+$ 4.40E-20(3.25E-20)$^+$ 5.52E-31(1.94E-30)$^+$ 2.83E-11(2.61E-11)$^+$ ${\bf5.83E-42}({\bf9.48E-42})$ $f_{3}$ 30 1.46E-11(1.23E-11)$^+$ 4.98E-10(2.39E-10)$^+$ 1.29E-16(2.15E-16)$^+$ 3.86E-06(1.32E-06)$^+$ ${\bf2.03E-23}({\bf1.62E-23})$ $f_{4}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 4.07E-17(5.44E-17)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 1.30E-14(1.13E-14)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{5}$ 30 1.65E-20(4.70E-20)$^-$ 1.97E-18(2.27E-18)$^-$ ${\bf8.85E-30}({\bf1.75E-29})$$^-$ 4.04E-10(3.04E-10)$^+$ 1.06E-15(3.87E-16) $f_{6}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{7}$ 30 ${\bf1.82E-04}({\bf3.39E-04})$$^-$ 3.57E-04(1.00E-03)$^-$ 1.42E+01(2.25E+01)$^-$ 5.69E-04(7.53E-04)$^-$ 3.59E+01(1.12E+01) $f_{8}$ 30 ${\bf1.28E-04}({\bf1.36E-04})$$^-$ 1.81E+01(1.11E+00)$^-$ 8.61E+00(2.09E+00)$^-$ 1.97E+01(5.80E-01)$^-$ 2.54E+01(8.26E-01) $f_{9}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 2.47E-04(1.35E-03)$^+$ 6.57E-04(2.50E-03)$^+$ 2.47E-07(7.34E-07)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{10}$ 30 ${\bf2.03E-02}({\bf6.76E-03})$$^-$ 3.56E-01(5.96E-02)$^+$ 2.35E-01(1.31E-01)$^+$ 1.97E+00(4.24E-01)$^+$ 1.84E-01(7.45E-02) $f_{11}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 1.75E+02(5.85E+01)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 1.66E-02(3.13E-02)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{12}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 3.14E-05(1.02E-07)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{13}$ 30 2.19E-30(2.11E-30)$^+$ 5.53E-21(5.00E-21)$^+$ 2.91E-29(1.14E-28)$^+$ 1.80E-13(3.36E-13)$^+$ ${\bf1.57E-32}({\bf0.00E+00})$ $f_{14}$ 30 1.36E-32(3.13E-34)$^+$ 4.07E-21(6.03E-21)$^+$ 1.75E-32(9.04E-33)$^+$ 1.31E-13(1.51E-13)$^+$ ${\bf1.35E-32}({\bf0.00E+00})$ $f_{15}$ 30 2.42E-10(1.35E-10)$^+$ 1.19E-10(6.88E-11)$^+$ 6.04E-15(1.66E-15)$^+$ 3.75E-06(1.88E-06)$^+$ ${\bf4.00E-15}({\bf1.02E-15})$ $f_{16}$ 30 1.99E-06(2.50E-06)$^+$ 2.30E+01(2.38E+00)$^+$ 5.48E-02(1.39E-01)$^+$ 3.30E+01(5.81E+00)$^+$ ${\bf5.18E-07}({\bf1.12E-06})$ $f_{17}$ 30 4.57E-26(7.30E-26)$^+$ 2.85E-20(4.26E-20)$^+$ ${\bf6.02E-32}({\bf1.48E-31})$$^-$ 1.44E-12(1.26E-12)$^+$ 2.09E-26(1.95E-26) $f_{18}$ 30 3.80E-21(4.03E-21)$^+$ 3.51E-19(3.01E-19)$^+$ 3.66E-04(2.01E-03)$^+$ 2.45E-11(2.36E-11)$^+$ ${\bf2.48E-21}({\bf4.41E-21})$ $f_{19}$ 30 ${\bf3.67E+02}({\bf4.24E+02})$$^-$ 6.53E+02(6.50E+02)$^+$ 9.64E+02(4.81E+02)$^+$ 6.92E+02(8.00E+02)$^+$ 6.23E+02(5.05E+02) $f_{20}$ 30 2.20E-01(1.67E-05)$^+$ 1.60E-02(1.92E-03)$^+$ 1.38E-03(1.06E-03)$^+$ 1.75E+00(1.47E+00)$^+$ ${\bf2.07E-13}({\bf2.86E-13})$ $+/\approx/-$ 10/5/3 16/1/3 12/4/4 16/2/2 -/-/- 表 5 DELU、SHADE、EPSDE、CoDE、SEFDE对50维测试函数的优化结果, 平均值(标准差)
Table 5 Compared data of 50 D optimization results on DELU、SHADE、EPSDE、CoDE、SEFDE, Mean (Std)
Fun $N$ DELU SHADE EPSDE CoDE SEFDE Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) $f_{1}$ 50 4.24E-42(4.13E-42)$^+$ 2.05E-53(3.34E-53)$^+$ 2.52E-50(7.95E-50)$^+$ 4.94E-21(5.74E-21)$^+$ ${\bf2.30E-57}({\bf2.91E-57})$ $f_{2}$ 50 1.92E-32(4.01E-32)$^+$ ${\bf1.35E-53}({\bf3.43E-53})$$^-$ 1.80E-51(3.75E-51)$^-$ 4.46E-22(6.24E-22)$^+$ ${\bf1.59E-42}({\bf2.13E-42})$ $f_{3}$ 50 7.87E-29(9.66E-29)$^+$ 2.22E-27(1.67E-27)$^+$ 5.63E-29(1.68E-28)$^+$ 6.05E-12(2.77E-12)$^+$ ${\bf5.52E-29}({\bf4.63E-29})$ $f_{4}$ 50 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 9.99E-17(3.51E-17)$^+$ 1.33E-16(4.68E-17)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{5}$ 50 7.64E-21(8.52E-21)$^+$ ${\bf1.03E-52}({\bf1.07E-52})$$^-$ 3.04E-50(7.79E-50)$^-$ 1.01E-20(1.21E-20)$^+$ 1.50E-27(1.66E-27) $f_{6}$ 50 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 5.00E-01(7.07E-01)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{7}$ 50 1.77E-04(9.11E-05)$^-$ ${\bf8.22E-07}({\bf1.80E-06})$$^-$ 2.22E+02(6.11E+01)$^+$ 2.80E-04(3.97E-04)$^-$ 9.43E+01(9.34E+00) $f_{8}$ 50 ${\bf2.59E-06}({\bf1.29E-06})$$^-$ 1.34E+01(1.83E+00)$^-$ 2.72E+01(1.81E+01)$^-$ 5.16E+01(2.64E+01)$^+$ 3.73E+01(3.74E+00) $f_{9}$ 50 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 2.46E-03(5.69E-03)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{10}$ 50 2.28E-02(1.14E-02)$^-$ 5.45E-02(3.68E-02)$^-$ ${\bf3.12E-03}({\bf4.16E-03})$$^-$ 1.11E-01(7.80E-02)$^-$ 2.19E-01(2.67E-01) $f_{11}$ 50 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 1.19E+01(3.79E+01)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{12}$ 50 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 3.39E-01(3.95E-01)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{13}$ 50 9.42E-33(1.44E-48)$^+$ 9.42E-33(1.44E-48)$^+$ 9.46E-33(6.47E-35)$^+$ 4.76E-25(4.46E-25)$^+$ ${\bf1.44E-33}({\bf1.11E-33})$ $f_{14}$ 50 1.35E-32(2.88E-48)$^+$ 1.35E-32(2.88E-48)$^+$ 1.40E-32(8.62E-34)$^+$ 1.06E-24(8.10E-25)$^+$ ${\bf1.33E-32}({\bf0.00E+00})$ $f_{15}$ 50 6.13E-15(1.95E-15)$^+$ 7.11E-15(0.00E+00)$^+$ 1.14E-14(4.97E-15)$^+$ 8.59E-12(4.34E-12)$^+$ ${\bf1.04E-15}({\bf3.89E-15})$ $f_{16}$ 50 2.45E-07(1.58E-07)$^+$ 1.84E-01(5.31E-02)$^+$ 4.42E+00(1.09E+01)$^+$ 4.27E+01(7.39E+00)$^+$ ${\bf1.96E-07}({\bf1.10E-07})$ $f_{17}$ 50 3.25E-30(1.48E-30)$^-$ ${\bf9.42E-33}({\bf1.44E-48})$$^-$ 1.00E-32(1.96E-33)$^-$ 9.53E-24(4.52E-25)$^+$ 1.68E-26(3.74E-26) $f_{18}$ 50 3.08E-32(5.23E-33)$^-$ ${\bf1.36E-32}({\bf3.90E-34})$$^-$ 1.10E-03(3.47E-03)$^+$ 1.29E-22(1.22E-22)$^+$ 1.71E-24(3.73E-24) $f_{19}$ 50 8.90E+03(6.05E+03)$^-$ ${\bf1.04E+03}({\bf6.24E+02})$$^-$ 1.26E+04(1.65E+03)$^-$ 1.09E+04(1.44E+03)$^-$ 1.66E+04(2.24E+03) $f_{20}$ 50 2.70E-01(4.67E-02)$^+$ 2.60E-01(5.16E-02)$^+$ 2.80E-01(4.22E-02)$^+$ 3.20E-01(4.22E-02)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $+/\approx/-$ 9/5/6 10/2/7 12/2/6 12/5/3 -/-/- 表 6 CLPSO、CMA-ES、GL-25和SEFDE的优化结果性能对比数据, 平均值(标准差)
Table 6 Compared data of optimization results on CLPSO、CMA-ES、GL-25 and SEFDE, Mean (Std)
Fun $N$ CLPSO CMA-ES GL-25 SEFDE Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) $f_{1}$ 30 9.13E-14(3.33E-14)$^+$ 1.97E-29(2.07E-30)$^+$ 4.78E-87(1.69E-86)$^+$ ${\bf1.26E-103}({\bf2.39E-103})$ $f_{2}$ 30 9.79E-15(5.44E-15)$^+$ 3.75E-28(4.61E-29)$^+$ 2.66E-78(1.25E-77)$^+$ ${\bf2.40E-112}({\bf3.36E-112})$ $f_{3}$ 30 4.48E-09(1.48E-09)$^+$ 2.03E-14(9.76E-16)$^+$ 2.98E-28(1.12E-27)$^+$ ${\bf3.68E-60}({\bf3.36E-112})$ $f_{4}$ 30 3.37E-16(7.98E-17)$^+$ 4.81E-17(5.60E-17)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{5}$ 30 9.33E-14(6.40E-14)$^+$ 1.50E-24(2.15E-25)$^+$ ${\bf4.38E-85}({\bf2.40E-84})$$^-$ 2.42E-45(2.29E-45) $f_{6}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{7}$ 30 4.73E+00(1.20E+00)$^+$ ${\bf2.95E-27}({\bf4.82E-28})$$^-$ 3.14E-02(8.77E-02)$^-$ 2.66E-01(2.28E-01) $f_{8}$ 30 1.95E+01(2.62E+00)$^-$ ${\bf2.66E-01}({\bf1.01E+00})$$^-$ 2.21E+01(6.19E-01)$^+$ 2.19E+01(5.02E-01) $f_{9}$ 30 2.40E-09(3.67E-09)$^+$ 1.40E-03(3.70E-03)$^+$ 1.61E-15(4.55E-15)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{10}$ 30 1.54E-01(2.09E-02)$^-$ 2.50E+02(1.23E+01)$^+$ 2.95E+00(1.00E+00)$^+$ ${\bf1.30E-01}({\bf6.42E-02})$ $f_{11}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 5.19E+03(6.07E+02)$^+$ 4.46E+03(1.36E+03)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{12}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 1.30E+01(2.50E+00)$^+$ 3.14E-05(3.97E-14)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{13}$ 30 1.13E-17(6.31E-18)$^+$ 9.13E-01(1.06E+00)$^+$ 8.82E-31(4.10E-30)$^+$ ${\bf1.57E-32}({\bf0.00E+00})$ $f_{14}$ 30 5.94E-17(3.31E-17)$^+$ 1.10E-03(3.35E-03)$^+$ 1.28E-30(5.05E-30)$^+$ ${\bf1.35E-32}({\bf0.00E+00})$ $f_{15}$ 30 9.76E-08(2.38E-08)$^+$ 1.93E+01(1.97E-01)$^+$ 1.09E-13(1.91E-13)$^+$ ${\bf7.55E-15}({\bf0.00E+00})$ $f_{16}$ 30 9.41E-07(6.48E-07)$^+$ 2.20E+02(5.64E+01)$^+$ 3.07E+01(2.71E+01)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $f_{17}$ 30 3.49E-16(1.59E-16)$^+$ 1.73E-02(3.93E-02)$^+$ 1.26E+02(1.09E+01)$^+$ ${\bf1.57E-32}({\bf0.00E+00})$ $f_{18}$ 30 2.52E-14(1.22E-14)$^+$ 2.20E-03(4.47E-03)$^+$ 1.97E+03(1.93E+02)$^+$ ${\bf1.35E-32}({\bf0.00E+00})$ $f_{19}$ 30 6.26E+03(9.93E+02)$^+$ ${\bf5.59E-10}({\bf7.36E-11})$$^-$ 2.49E+03(4.07E+02)$^+$ 5.60E+02(1.36E+02) $f_{20}$ 30 2.33E-04(9.49E-05)$^+$ 8.96E-02(1.46E-01)$^+$ 3.55E-04(9.33E-04)$^+$ ${\bf1.23E-10}({\bf2.47E-10})$ $+/\approx/-$ 15/3/2 16/1/3 16/2/2 -/-/- 表 7 SHADE、IDE、ZPEDE、SinDE、SEFDE的优化结果性能对比数据, 平均值(标准差)
Table 7 Compared data of optimization results on SHADE、IDE、ZEPDE、SinDE、SEFDE, Mean (Std)
Fun N SHADE IDE ZEPDE SinDE SEFDE Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) Mean error (Std Dev) $F_{1}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 2.27E-13(1.53E-28)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $F_{2}$ 30 ${\bf2.66E+04}({\bf1.13E+04})$$^-$ 1.68E+06(4.23E+05)$^-$ 1.97E+05(7.53E+04)$^-$ 2.66E+06(8.33E+05)$^-$ 2.85E+07(4.18E+06) $F_{3}$ 30 8.80E+05(1.96E+06)$^+$ 1.38E+05(1.85E+05)$^+$ 1.50E+06(2.07E+06)$^+$ 1.01E+05(3.77E+05)$^+$ ${\bf5.80E+04}({\bf3.74E+04})$ $F_{4}$ 30 ${\bf1.61E-03}(1.41E-03)$$^-$ 6.85E+03(1.10E+03)$^+$ 7.66E-01(4.78E-01)$^-$ 8.28E+03(1.54E+03)$^+$ 3.13E+03(1.12E+03) $F_{5}$ 30 ${\bf0.00E+00}(0.00E+00)$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 1.14E-13(7.65E-29)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $F_{6}$ 30 4.28E+01(5.52E+00)$^+$ 4.34E+01(2.62E-04)$^+$ 4.34E+01(3.18E-13)$^+$ 4.34E+01(1.44E-14)$^+$ ${\bf1.52E+01}({\bf1.95E-01})$ $F_{7}$ 30 2.33E+01(9.32E+00)$^-$ 3.18E+00(1.55E+00)$^-$ 1.37E+01(4.88E+00)$^-$ ${\bf6.10E-01}(5.97E-01)$$^-$ 2.80E+01(7.96E+00) $F_{8}$ 30 2.09E+01(1.68E-01)$^+$ 2.11E+01(2.44E-02)$^+$ 2.11E+01(1.17E-01)$^+$ 2.11E+01(3.59E-02)$^+$ ${\bf2.09E+01}({\bf2.78E-02})$ $F_{9}$ 30 5.54E+01(1.98E+00)$^+$ 3.56E+01(5.54E+00)$^+$ 3.74E+01(5.85E+00)$^+$ 3.48E+01(4.34E+00)$^+$ ${\bf3.04E+01}({\bf6.01E-01})$ $F_{10}$ 30 7.37E-02(3.67E-02)$^+$ 4.38E-02(2.17E-02)$^+$ 1.37E-01(6.96E-02)$^+$ 7.93E-02(3.57E-02)$^+$ ${\bf2.71E-02}({\bf1.69E-03})$ $F_{11}$ 30 ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ ${\bf0.00E+00}({\bf0.00E+00})$$^\approx$ 3.65E-01(6.12E-01)$^+$ 5.92E+00(2.86E+00)$^+$ ${\bf0.00E+00}({\bf0.00E+00})$ $F_{12}$ 30 5.86E+01(1.11E+01)$^-$ 6.89E+01(8.82E+00)$^-$ 6.04E+01(1.76E+01)$^-$ ${\bf5.61E+01}(1.41E+01)$$^-$ 1.21E+02(9.62E+00) $F_{13}$ 30 1.45E+02(1.95E+01)$^+$ 1.34E+02(2.28E+01)$^+$ 1.32E+02(3.62E+01)$^+$ 1.39E+02(3.41E+01)$^+$ ${\bf1.24E+02}({\bf1.29E+00})$ $F_{14}$ 30 ${\bf3.45E-02}({\bf1.93E-02})$$^-$ 1.17E+02(8.38E+01)$^+$ 4.83E+00(2.70E+00)$^-$ 2.34E+02(9.23E+01)$^-$ 5.85E+00(2.05E+00) $F_{15}$ 30 6.82E+03(4.41E+02)$^+$ 6.54E+03(5.91E+02)$^+$ 6.59E+03(9.36E+03)$^+$ 6.80E+03(1.00E+03)$^+$ ${\bf5.95E+03}({\bf2.37E+02})$ $F_{16}$ 30 1.28E+00(2.07E-01)$^-$ 1.59E+00(2.36E-01)$^-$ ${\bf7.82E-01}({\bf6.74E-01})$$^-$ 2.08E+00(3.66E-01)$^-$ 1.02E+02(1.66E-01) $F_{17}$ 30 5.08E+01(4.27E-14)$^+$ 5.92E+01(1.41E+00)$^+$ 5.11E+01(1.60E-01)$^+$ 6.52E+01(3.47E+00)$^+$ ${\bf4.21E+01}({\bf6.06E+01})$ $F_{18}$ 30 1.37E+02(1.29E+01)$^-$ 1.68E+02(1.27E+01)$^-$ ${\bf1.03E+02}({\bf1.19E+01})$$^-$ 1.41E+02(2.27E+01)$^-$ 3.01E+02(9.05E+00) $F_{19}$ 30 2.64E+00(2.83E-01)$^-$ ${\bf2.24E+00}({\bf3.66E-01})$$^-$ 3.71E+00(7.55E-01)$^-$ 4.85E+00(8.82E-01)$^-$ 1.03E+02(1.58E-01) $F_{20}$ 30 1.93E+01(7.70E-01)$^-$ 1.93E+01(4.47E-01)$^-$ 1.97E+01(7.88E-01)$^-$ ${\bf1.92E+01}(7.52E-01)$$^-$ 1.12E+02(9.46E-02) $F_{21}$ 30 8.45E+02(3.63E+02)$^+$ 7.32E+02(3.82E+02)$^+$ 6.33E+02(4.48E+02)$^+$ 5.84E+02(4.22E+02)$^+$ ${\bf5.03E+02}({\bf4.03E+01})$ $F_{22}$ 30 ${\bf1.33E+01}({\bf7.12E+00})$$^-$ 6.88E+01(2.03E+01)$^-$ 4.23E+02(5.75E+02)$^+$ 3.51E+02(2.72E+02)$^+$ 2.41E+02(1.24E+01) $F_{23}$ 30 7.63E+03(6.58E+02)$^+$ 7.32E+03(6.92E+02)$^+$ 7.02E+03(8.73E+02)$^+$ 6.59E+03(8.47E+02)$^+$ ${\bf6.38E+03}({\bf2.25E+02})$ $F_{24}$ 30 2.34E+02(1.01E+01)$^+$ 2.02E+02(1.14E+00)$^-$ 2.35E+02(1.09E+01)$^+$ ${\bf2.00E+02}(1.34E-01)$$^-$ 2.32E+02(1.18E+00) $F_{25}$ 30 3.40E+02(3.09E+01)$^+$ 3.03E+02(1.09E+01)$^+$ 3.23E+02(1.31E+01)$^+$ 2.97E+02(1.33E+01)$^+$ ${\bf2.58E+02}({\bf1.05E+01})$ $F_{26}$ 30 2.58E+02(8.08E+01)$^-$ ${\bf2.23E+02}({\bf4.46E+01})$$^-$ 2.27E+02(6.20E+01)$^-$ 2.76E+02(5.96E+01)$^-$ 3.02E+02(1.90E-01) $F_{27}$ 30 9.36E+02(3.07E+02)$^+$ ${\bf3.58E+02}({\bf3.30E+01})$$^-$ 9.38E+02(1.40E+02)$^+$ 4.75E+02(1.55E+02)$^-$ 6.11E+02(1.62E+02) $F_{28}$ 30 4.58E+02(4.13E+02)$^+$ ${\bf4.00E+02}({\bf0.00E+00})$$^\approx$ ${\bf4.00E+02}({\bf0.00E+00})$$^\approx$ ${\bf4.00E+02}(0.00E+00)$$^\approx$ ${\bf4.00E+02}({\bf0.00E+00})$ $+/\approx/-$ 14/3/11 13/4/11 15/3/10 16/1/11 -/-/- -
[1] Storn R, Price K. Differential evolution - a simple and efficient heuristic for global optimization over continuous spaces. Journal of Global Optimization, 1997, 11(4): 341-359 doi: 10.1023/A:1008202821328 [2] Das S, Suganthan P N. Differential evolution: a survey of the state-of-the-art. IEEE Transactions on Evolutionary Computation, 2011, 15(1): 4-31 doi: 10.1109/TEVC.2010.2059031 [3] Raghunathan T, Ghose D. Differential evolution based 3-D guidance law for a realistic interceptor model. Applied Soft Computing, 2014, 16: 20-33 doi: 10.1016/j.asoc.2013.11.017 [4] 胡蓉, 钱斌.一种求解随机有限缓冲区流水线调度的混合差分进化算法.自动化学报, 2009, 35(12): 1580-1586 doi: 10.3724/SP.J.1004.2009.01580Hu Rong, Qian Bin. A hybrid differential evolution algorithm for stochastic flow shop scheduling with limited buffers. Acta Automatica Sinica, 2009, 35(12): 1580-1586 doi: 10.3724/SP.J.1004.2009.01580 [5] Pandit M, Srivastava L, Sharma M. Environmental economic dispatch in multi-area power system employing improved differential evolution with fuzzy selection. Applied Soft Computing, 2015, 28: 498-510 doi: 10.1016/j.asoc.2014.12.027 [6] Zhang G J, Zhou X G, Yu X F, Hao X H, Yu L. Enhancing protein conformational space sampling using distance profile-guided differential evolution. IEEE/ACM Transactions on Computational Biology and Bioinformatics, 2017, 14(6): 1288-1301 doi: 10.1109/TCBB.2016.2566617 [7] 朱腾, 王京春, 熊智华.基于改进DE-NMPC的酸碱中和反应pH值控制.自动化学报, 2010, 36(1): 159-163 doi: 10.3724/SP.J.1004.2010.00159Zhu Teng, Wang Jing-Chun, Xiong Zhi-Hua. DE-based nonlinear model predictive control of a pH neutralization process. Acta Automatica Sinica, 2010, 36(1): 159-163 doi: 10.3724/SP.J.1004.2010.00159 [8] Qiao D P, Pang G K H. A modified differential evolution with heuristic algorithm for nonconvex optimization on sensor network localization. IEEE Transactions on Vehicular Technology, 2016, 65(3): 1676-1689 doi: 10.1109/TVT.2015.2409319 [9] Das S, Mullick S S, Suganthan P N. Recent advances in differential evolution - an updated survey. Swarm and Evolutionary Computation, 2016, 27: 1-30 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=abbffecca110c67e0656d3a7220658fc [10] Dragoi E N, Dafinescu V. Parameter control and hybridization techniques in differential evolution: a survey. Artificial Intelligence Review, 2016, 45(4): 447-470 doi: 10.1007/s10462-015-9452-8 [11] Cheng M Y, Tran D H. Two-phase differential evolution for the multiobjective optimization of time - cost tradeoffs in resource-constrained construction projects. IEEE Transactions on Engineering Management, 2014, 61(3): 450-461 doi: 10.1109/TEM.2014.2327512 [12] Liu Z Z, Wang Y, Yang S X, Cai Z X. Differential evolution with a two-stage optimization mechanism for numerical optimization. In: Proceedings of the 2016 IEEE Congress on Evolutionary Computation (CEC). Vancouver, BC, Canada: IEEE, 2016. 3170-3177 [13] Tang L X, Dong Y, Liu J Y. Differential evolution with an individual-dependent mechanism. IEEE Transactions on Evolutionary Computation, 2015, 19(4): 560-574 doi: 10.1109/TEVC.2014.2360890 [14] Yu W J, Shen M E, Chen W N, Zhan Z H, Gong Y J, Lin Y, Liu Q, Zhang J. Differential evolution with two-level parameter adaptation. IEEE Transactions on Cybernetics, 2014, 44(7): 1080-1099 doi: 10.1109/TCYB.2013.2279211 [15] 张贵军, 陈铭, 周晓根.动态小生境半径两阶段多模态差分进化算法.控制与决策, 2016, 31(7): 1185-1191 http://d.old.wanfangdata.com.cn/Periodical/kzyjc201607006Zhang Gui-Jun, Chen Ming, Zhou Xiao-Gen. Two-stage differential evolution algorithm using dynamic niche radius for multimodal optimization. Control and Decision, 2016, 31(7): 1185-1191 http://d.old.wanfangdata.com.cn/Periodical/kzyjc201607006 [16] Price K, Storn R M, Lampinen J A. Differential Evolution - A Practical Approach to Global Optimization. Berlin Heidelberg: Springer, 2005: 1-24 [17] 张贵军, 何洋军, 郭海锋, 冯远静, 徐建明.基于广义凸下界估计的多模态差分进化算法.软件学报, 2013, 24(6): 1177-1195 http://d.old.wanfangdata.com.cn/Periodical/rjxb201306002Zhang Gui-Jun, He Yang-Jun, Guo Hai-Feng, Feng Yuan-Jing, Xu Jian-Ming. Differential evolution algorithm for multimodal optimization based on abstract convex underestimation. Journal of Software, 2013, 24(6): 1177-1195 http://d.old.wanfangdata.com.cn/Periodical/rjxb201306002 [18] 周晓根, 张贵军, 郝小虎.局部抽象凸区域剖分差分进化算法.自动化学报, 2015, 41(7): 1315-1327 doi: 10.16383/j.aas.2015.c140680Zhou Xiao-Gen, Zhang Gui-Jun, Hao Xiao-Hu. Differential evolution algorithm with local abstract convex region partition. Acta Automatica Sinica, 2015, 41(7): 1315-1327 doi: 10.16383/j.aas.2015.c140680 [19] Beliakov G. Extended cutting angle method of global optimization. Pacific Journal of Optimization, 2008, 4(1): 153-176 [20] Beliakov G, Lim K F. Challenges of continuous global optimization in molecular structure prediction. European Journal of Operational Research, 2007, 181(3): 1198-1213 doi: 10.1016/j.ejor.2005.08.033 [21] 周晓根, 张贵军, 郝小虎, 俞立.一种基于局部Lipschitz下界估计支撑面的差分进化算法.计算机学报, 2016, 39(12): 2631-2651 doi: 10.11897/SP.J.1016.2016.02631Zhou Xiao-Gen, Zhang Gui-Jun, Hao Xiao-Hu, Yu Li. Differential evolution algorithm based on local lipschitz underestimate supporting hyperplanes. Chinese Journal of Computers, 2016, 39(12): 2631-2651 doi: 10.11897/SP.J.1016.2016.02631 [22] Ali M M, Khompatraporn C, Zabinsky Z B. A numerical evaluation of several stochastic algorithms on selected continuous global optimization test problems. Journal of Global Optimization, 2005, 31(4): 635-672 doi: 10.1007/s10898-004-9972-2 [23] Sarker R A, Elsayed S M, Ray T. Differential evolution with dynamic parameters selection for optimization problems. IEEE Transactions on Evolutionary Computation, 2014, 18(5): 689-707 [24] Dragoi E N, Dafinescu V. Parameter control and hybridization techniques in differential evolution: a survey. Artificial Intelligence Review, 2015, 45(4): 447-470 [25] Zhou X G, Zhang G J. Abstract convex underestimation assisted multistage differential evolution. IEEE Transactions on Cybernetics, 2017, 47(9): 2730-2741 doi: 10.1109/TCYB.2017.2710626 [26] Zhou X G, Zhang G J, Hao X H, Yu L. A novel differential evolution algorithm using local abstract convex underestimate strategy for global optimization. Computers & Operations Research, 2016, 75: 132-149 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dd11fc41fe4a435f970157e38de1f775 [27] Tanabe R, Fukunaga A. Success-history based parameter adaptation for differential evolution. In: Proceedings of the 2013 IEEE Congress on Evolutionary Computation (CEC). Cancun, Mexico: IEEE, 2013. 71-78 [28] Mallipeddi R, Suganthan P N, Pan Q K, Tasgetiren M F. Differential evolution algorithm with ensemble of parameters and mutation strategies. Applied Soft Computing, 2011, 11(2): 1679-1696 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cface92266298dbf8ec3818bbe607897 [29] Wang Y, Cai Z X, Zhang Q F. Differential evolution with composite trial vector generation strategies and control parameters. IEEE Transactions on Evolutionary Computation, 2011, 15(1): 55-66 doi: 10.1109/TEVC.2010.2087271 [30] Corder G W, Foreman D I. Nonparametric Statistics for Non-Statisticians: A Step-by-Step Approach. Hoboken, NJ: Wiley 2009. 38-55 [31] Shang Y W, Qiu Y H. A note on the extended Rosenbrock function. Evolutionary Computation, 2006, 14(1): 119-126 [32] Liang J J, Qin A K, Suganthan P N, Baskar S. Comprehensive learning particle swarm optimizer for global optimization of multimodal functions. IEEE Transactions on Evolutionary Computation, 2006, 10(3): 281-295 [33] Hansen N, Ostermeier A. Completely derandomized self-adaptation in evolution strategies. Evolutionary Computation, 2001, 9(2): 159-195 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=58e04c72ec259f9411f7e06253975afb [34] García-Martínez C, Lozano M, Herrera F, Molina D, Sánchez A M. Global and local real-coded genetic algorithms based on parent-centric crossover operators. European Journal of Operational Research, 2008, 185(3): 1088-1113 doi: 10.1016/j.ejor.2006.06.043 [35] Fan Q Q, Yan X F. Self-adaptive differential evolution algorithm with zoning evolution of control parameters and adaptive mutation strategies. IEEE Transactions on Cybernetics, 2016, 46(1): 219-232 doi: 10.1109/TCYB.2015.2399478 [36] Draa A, Bouzoubia S, Boukhalfa I. A sinusoidal differential evolution algorithm for numerical optimisation. Applied Soft Computing, 2015, 27: 99-126 doi: 10.1016/j.asoc.2014.11.003 [37] Liang J J, Qu B Y, Suganthan P N, Hernández-Díaz A G. Problem definitions and evaluation criteria for the CEC 2013 special session on real-parameter optimization, Technical Report 201212, Computational Intelligence Laboratory, Zhengzhou University, Zhengzhou, China and Nanyang Technological University, Singapore, 2013. 期刊类型引用(10)
1. 宋红超,王昕,王振雷. 基于切换的非线性多模型二阶段广义预测控制. 控制理论与应用. 2024(11): 2147-2156 . 百度学术
2. 赵建国,杨春雨. 复杂工业过程非串级双速率组合分散运行优化控制. 自动化学报. 2023(01): 172-184 . 本站查看
3. 王建琴. 基于协同过滤算法的医院财务工作轮岗流程智能控制模型. 微型电脑应用. 2022(06): 54-56 . 百度学术
4. 陈尉. INFIT系统在火电机组辅助调频控制中的应用. 工业控制计算机. 2021(01): 26-27+30 . 百度学术
5. 陈铁锋,杜文. 基于距离度量的火电机组协调系统软切换预测控制. 热能动力工程. 2021(03): 87-92+121 . 百度学术
6. 王素珍,刘建锋. 改进型BP神经网络的非线性多模型自适应控制. 电光与控制. 2021(08): 1-5 . 百度学术
7. 李军,黄卫剑,万文军,刘哲. 一种新型反馈控制器的研究与应用. 控制理论与应用. 2020(02): 411-422 . 百度学术
8. 陈钊,乔侨. 分层结构多模型预测控制在630MW火电机组协调系统中的应用. 工业控制计算机. 2020(10): 4-6+10 . 百度学术
9. 王岩,王昕,王振雷. 一类时变系统的多模型切换动态调节控制算法. 控制理论与应用. 2020(12): 2501-2510 . 百度学术
10. 唐伟强,龙文堃,孙丽娟,黄小丽. 基于聚类方法和神经网络的非线性系统多模型自适应控制. 系统工程与电子技术. 2019(09): 2100-2106 . 百度学术
其他类型引用(7)
-