2.845

2023影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

斐波那契树优化算法求解多峰函数全局最优解的可达性分析

董易 施心陵 王霞 王耀民 吕丹桔 张松海 李孙寸

高忠顺, 周炳海. 基于分枝搜索的带缓冲集束型设备调度与性能分析. 自动化学报, 2016, 42(1): 81-88. doi: 10.16383/j.aas.2016.c140142
引用本文: 董易, 施心陵, 王霞, 王耀民, 吕丹桔, 张松海, 李孙寸. 斐波那契树优化算法求解多峰函数全局最优解的可达性分析. 自动化学报, 2018, 44(9): 1679-1689. doi: 10.16383/j.aas.2017.c160703
GAO Zhong-Shun, ZHOU Bing-Hai. The Scheduling and Performance Analysis of Cluster Tools with Buffers Based on Branch Searching. ACTA AUTOMATICA SINICA, 2016, 42(1): 81-88. doi: 10.16383/j.aas.2016.c140142
Citation: DONG Yi, SHI Xin-Ling, WANG Xia, WANG Yao-Min, LV Dan-Ju, ZHANG Song-Hai, LI Sun-Cun. On Accessibility of Fibonacci Tree Optimization Algorithm for Global Optima of Multi-modal Functions. ACTA AUTOMATICA SINICA, 2018, 44(9): 1679-1689. doi: 10.16383/j.aas.2017.c160703

斐波那契树优化算法求解多峰函数全局最优解的可达性分析

doi: 10.16383/j.aas.2017.c160703
基金项目: 

国家自然科学基金 61661050

云南省教育厅科学研究基金项目 2017ZZX229

详细信息
    作者简介:

    施心陵 云南大学信息学院教授.主要研究方向为智能优化算法, 自适应信号处理与信息系统, 医学电子学.E-mail:xlshi@ynu.edu.cn

    王霞 云南大学信息学院博士研究生.2010年获得中南大学物理科学与技术学院硕士学位.主要研究方向为智能优化算法.E-mail:wangxiacsu@163.com

    王耀民 云南大学信息学院博士研究生.2013年获得南京邮电大学通信与信息工程学院硕士学位.主要研究方向为智能优化算法.E-mail:18988081898@189.cn

    吕丹桔 云南大学信息学院博士研究生.2013年获得云南大学信息学院硕士学位.主要研究方向为智能林业信息处理.E-mail:lvdanjv@hotmail.com

    张松海 云南大学信息学硕士研究生.主要研究方向为智能优化算法.E-mail:hai_zs@sina.com

    李孙寸 云南大学信息学院硕士研究生.主要研究方向为智能优化算法.E-mail:lisuncun@outlook.com

    通讯作者:

    董易 云南大学信息学院博士研究生.2009年获得云南大学信息学院硕士学位.主要研究方向为智能优化算法.本文通信作者.E-mail:yeo003@163.com

On Accessibility of Fibonacci Tree Optimization Algorithm for Global Optima of Multi-modal Functions

Funds: 

National Natural Science Foundation of China 61661050

Science Research Foundation of Yunnan Provincial Department of Education 2017ZZX229

More Information
    Author Bio:

    Professor at the School of Information Science and Engineering, Yunnan University. His research interest covers intelligent optimization algorithm, adaptive signal processing, information system, and medical electronics

    Ph. D. candidate at the School of Information Science and Engineering, Yunnan University. She received her master degree from the College of Physical Science and Technology, Central South University in 2010. Her main research interest is intelligent optimization algorithm

    Ph. D. candidate at the School of Information Science and Engineering, Yunnan University. He received his master degree from the College of Telecommunications and Information Engineering, Nangjing University of Posts and Telecommunications in 2013. His main research interest is intelligent optimization algorithm

    Ph. D. candidate at the School of Information Science and Engineering, Yunnan University. She received her master degree from the School of Information Science and Engineering, Yunnan University in 2013. Her main research interest is intelligent forestry information processing

    Master student at the School of Information Science and Engineering, Yunnan University. His main research interest is intelligent optimization algorithm

    Master student at the Science and Engineering, Yunnan University. Her main research interest is intelligent optimization algorithm

    Corresponding author: DONG Yi Ph. D. candidate at the School of Information Science and Engineering, Yunnan University. He received his master degree from the School of Information Science and Engineering, Yunnan University in 2009. His main research interest is intelligent optimization algorithm. Corresponding author of this paper
  • 摘要: 为分析和验证斐波那契树优化算法(Fibonacci tree optimization algorithm,FTO)求解多峰函数全局最优解的算法性能,对算法的可达性问题进行研究.本文基于斐波那契法构造一个斐波那契树结构,在搜索空间中进行全局、局部交替搜索,不易陷入局部最优解.对斐波那契树优化算法基于该结构的可达性进行分析和证明.通过跟踪算法求解过程中坐标点的累积分布仿真实验和到达率的对比实验,分析和验证了算法求解多峰函数全局最优解的可达性.
  • 集束型设备已被广泛地应用于半导体制造的晶圆加工中[1-2].由于集束型设备加工晶圆时存在多重入流、驻留及资源约束等现象[3],使得集束型设备的调度变得非常复杂.近年来,为了提高设备的产能,集束型设备的结构随着不断改进也出现很多不同的变化,其中一些学者另辟蹊径,希望通过合理地设置缓冲模块来提高设备的产能[4].而对于增加缓冲模块的新型结构的集束型设备来说,现有的调度方法定然不能满足其生产要求,这就需要新的调度方法来支持其调度.

    集束型设备一般由多个晶圆加工模块(Processing module,PM)、一个晶圆搬运模块(Transport module,TM)和负责晶圆输入输出的卡匣模块(Cassette modules,CM) 组成.作为晶圆加工特有的工艺要求,驻留约束限制了晶圆在加工模块完成加工后所能停留的时间上限(Upper bound) ,即当晶圆在加工完成后,必须在规定时间范围内离开加工模块,否则残余气体或不适宜温度会导致晶圆出现质量问题,甚至不良品.驻留现象广泛存在于半导体制造过程中[5].针对一般的驻留约束问题,大部分研究都是通过调整晶圆的加工开始时间来满足晶圆的驻留约束[6],一些学者通过设置缓冲模块的方法,使晶圆在加工完成后能够尽早离开加工模块,释放资源,提前下一个晶圆的加工开始时间.

    目前,对于一般的集束型设备的周期调度研究已经比较成熟. Perkinson等[7] 与Venkatesh等[8]在不考虑驻留约束的情况下,分别采用推策略和交换策略来获得单臂与双臂集束型设备的最优的基本周期(Fundamental period,FP). Lim等[9] 为了有效提高整体调度周期,考虑晶圆搬运延迟,并提出了一个基于时间表技术的有效实时调度算法.

    然而,前面的文献没有考虑驻留约束对集束型设备的调度影响.目前已有相关文献对带驻留约束的集束型设备调度问题进行了研究,Wu等[10]为带驻留约束的集束型设备建立了Petri-网模型,提出了一个封闭式调度算法来获得最优周期的调度.Rostami等[11]运用线性规划和启发式算法来调度带驻留约束的单集束型设备.Zhou等[12]提出了基于侦误的调度启发式算法用于解决带驻留约束的集束型设备调度问题.Rostami等[13]对集束型设备带有加工模块和搬运模块的双驻留约束模型提出了一种调度算法,得到了问题的近优解.

    上述文献虽然考虑了驻留约束,却没有考虑缓冲模块对设备产能的影响.对于集束型设备的缓冲模块的研究目前正在起步,Ding等[14]在集束型设备群的研究中,对各个集束型之间的缓冲连接进行了研究.Dawande等[15]比较了制造单元里带输出缓冲和无缓冲的模型,证明带有输出缓冲的制造单元比双爪不带缓冲的制造单元的生产率提高了20%.Drobouchevitch等[16]表明在制造单元的每个机器已存在输出缓冲的基础上增加一个输入缓冲不会提高设备的产能.

    上述文献表明,目前对同时带驻留约束和缓冲模块的集束型设备调度问题尚无人研究.本文旨在上述文献的研究基础上,对带驻留约束的集束型设备进行了研究,通过在加工模块之间设置缓冲模块的方法,提高系统产能.该调度问题主要是确定优化的搬运顺序和搬运时间点以获得最小的周期调度.首先,对调度问题在不同搬运模式下建立了相应数学模型;再提出两种不同的调度算法;最后,评价算法的同时,分析缓冲对系统产能的影响.

    集束型设备里,晶圆由卡匣模块进入,按顺序到每个加工模块加工,完成所有加工后回到卡匣模块.本文的研究对象为加工模块之间带有缓冲模块的单臂集束型设备,如图 1所示.

    图 1  带缓冲的集束型设备示意图
    Fig. 1  The cluster tool with buffers

    为有效地描述带缓冲的集束型晶圆制造设备的调度问题,做如下的基本定义与假设:1) 加工模块,搬运模块和缓冲模块在同一时间均只能处理一个晶圆.2) 加工模块上存在驻留约束现象.3) 搬运模块搬运晶圆的时间和装卸时间是确定的.4) 每个晶圆的加工工艺流程都相同.5) 若晶圆进入缓冲模块,缓冲模块被视作虚拟的加工模块,其加工时间为零,且不存在驻留约束;否则,晶圆跳过缓冲模块直接到下一个加工模块进行加工.6) 晶圆在卡匣模块供应充足,晶圆完成最后一道加工工序后直接被搬运到卡匣模块.

    为了清晰地表述调度问题,现定义如下符号与变量:

    $m$ : 集束型设备加工模块的数量.

    $p_i$ : 晶圆 $j$ 在加工模块 $i$ 上的加工时间.

    $\delta_i$ : 搬运模块把晶圆从 $M_i$ 搬运到 $M_{i+1}$ 的时间,包括装载和卸载晶圆的时间 $\delta_i=\delta$ .

    $\theta_{ij}$ : 搬运模块从 $M_i$ 到 $M_j$ 的空移动的时间 $\theta_{ij}=\min(|i-j|,2m-|i-j|)\cdot\theta$ .

    $λ_k$ : 表示晶圆 $j$ 是否经过模块 $M_k$ ,如果是,那么 $λ_k=1$ ;否则, $λ_k=0$ .

    $\varphi(i)$ : 表示搬运模块从第 $\varphi(i)$ 模块搬运到 $(\varphi(i)+2-λ_{\varphi(i)+1})$ 模块.

    $s_j$ : 搬运模块的搬运模式.

    $T$ : 生产每个晶圆的周期时间.

    定义 1. 一个单元搬运模式(One-unit TM move style)是指一个晶圆进入加工模块到下一个晶圆开始进入加工模块之间的搬运模块的所有搬运动作.可以通过如下形式来表示: $s=[\varphi(0) ,\varphi(1) ,\varphi(2) ,\cdots,\varphi(m)]$ ,其中 $s$ 是由 $\{0,1,2,\cdots,m\}$ 的排列组成.

    由于搬运模式的重复性,假设 $\varphi(0)=0$ ,不失一般性, $\varphi[i\pm(\sum_{k=0}^{2m-1}\lambda_k+1)]=\varphi(i)$ . 两个加工模块的单集束型设备的搬运模式包括如下: $s_1=[0,1,3]$ , $s_2=[0,3,1]$ , $s_3=[0,1,2,3]$ , $s_4=[0,1,3,2]$ , $s_5=[0,2,1,3]$ , $s_6=[0,2,3,1]$ , $s_7=[0,3,1,2]$ , $s_8=[0,3,2,1]$ .

    定义 2. 给定一个搬运模式 $s$ ,其中 $\varphi(k)=i$ ,那么 $\varphi(k)$ 的反函数为 $\varphi^-(k)$ , $0\leq\varphi^-(k)\leq m$ , $k=0,1,\cdots,m$ , $\varphi^-(\varphi(i))=i$ ,且 $\varphi(\varphi^-(k))=k$ .

    从定义可知,对给定的搬运模式来说,如果搬运模块从 $M_k$ 卸载晶圆,该搬运动作处于单元搬运模式的 $\varphi^-(k)$ th位置,那么当存在 $\varphi^-(k)>\varphi^-(k-2+λ_{k-1})$ 时,从 $M_k$ 装载和卸载的同一片晶圆的搬运动作处于同一个单元搬运模式内;当存在 $\varphi^-(k)<\varphi^-(k-2+λ_{k-1})$ 时,从 $M_k$ 装载和卸载的同一片晶圆的搬运动作却分别处于两个先后的单元搬运模式内.所以 $u_k$ 被定义如下:

    \begin{align}u_k=\left\{\begin{array}{c}1,\quad \varphi^-(k)<\varphi^-(k-2+λ_{k-1})\\0,\quad \varphi^-(k)>\varphi^-(k-2+λ_{k-1})\end{array} \right. \end{align}

    (1)

    任何单元搬运模式内都存在一些不可避免的搬运动作,如晶圆卸载和装载、晶圆的搬运等,用 $IC(s)$ 来表示给定单元搬运模式 $s$ 中必然存在的搬运动作,其表达式如下:

    $IC(s)=\sum\limits_{i=0}^{2m-1}{(}{{\delta }_{\varphi (i)}}+{{\theta }_{\varphi (i)}}+2-{{\lambda }_{\varphi (i)+1}},\varphi (i+2-{{\lambda }_{i+1}}))\text{ }$

    (2)

    在单元搬运模式 $s$ 下,一片晶圆从装载至模块 $M_k$ 到加工完后被卸载这段时间内,搬运模块所耗费的所有的繁忙时间用 $z_k(s)$ 表示:

    ${{z}_{k}}(s)={{\theta }_{k,{{\tau }_{k}}(s)}}+\sum\limits_{i}^{{{b}_{k}}}{(}{{\delta }_{\varphi (i)}}+{{\theta }_{\varphi (i)}}+2-{{\lambda }_{\varphi (i)+1}},\varphi (i+2-{{\lambda }_{i+1}}))$

    (3)

    其中, $\tau_k(s)$ 为搬运动作 $k-1$ 的下一搬运动作,即 $\tau_k(s)=\varphi[\varphi^-(k-2+λ_{k-1})+1]$ ,而 $i$ 的初始值为: $\varphi^-(k-2+λ_{k-1})+1$ ,且 $b_k=\varphi^-(k)-1+u_k(\sum_{j=0}^{2m}λ_j+1) $ .

    若存在常数 $T$ 使得集束型设备生产状态经 $T$ 时间回到相同的状态,则该时间 $T$ 被称作周期时间.

    令 $w_{kj}$ 是指搬运模块在 $M_k$ 等待搬运晶圆 $j$ 的时间.由于周期生产的性质,假设 $\sigma(j\pm n)=\sigma(j)$ , $w_{k,\sigma(j\pm n)}=w_{k,\sigma(j)}$ ,从假设5) 和6) 可知,当 $k$ 是偶数时, $w_{kj}=0$ .

    若 $T(s)$ 是相应的搬运模式 $s$ 的周期时间,则目标函数 $T(s)$ 就有:

    \begin{align}T(s)=\min\left(IC(s)+\sum_{j=1}^n\sum_{k=0}^mw_{kj}\right) \end{align}

    (4)

    由加工时间需求的约束与假设1) 可知,晶圆在 $M_k$ 上加工需满足如下约束:

    $\sum\limits_{i}^{{{\varphi }^{-}}(k)}{{{w}_{\varphi (i)}}}\ge \max ({{p}_{k}}-{{z}_{k}},0),{{u}_{k}}=0$

    (5)

    $\sum\limits_{i=1}^{\Sigma {{\lambda }_{l}}}{+}\sum\limits_{i}^{{{\varphi }^{-}}(k)}{{{w}_{\varphi (i)}}}\ge \max ({{p}_{k}}-{{z}_{k}},0),{{u}_{k}}=1$

    (6)

    其中,式(5) 和式(6) 中 $i$ 的初始值为: $i=\varphi^-(k-2+λ_{k-1})+2-λ_{\varphi^-(k-2+λ_{k-1})+1}$ .

    由假设2) 的驻留约束有:

    \begin{align}z_k-p_k\leq a_k \end{align}

    (7)

    搬运模块的等待时需满足非负:

    \begin{align}w_k\geq 0 \end{align}

    (8)

    所以,对于给定的搬运模式 $s$ ,要优化目标函数(4) ,服从约束(5) $\sim$ (8) .其中的 $IC(s)$ 与决策变量 $w_{kj}$ 相互独立.

    本文是在搬运模式下建立相应的数学模型,在前面数学建模里,该模型的复杂度为 ${\rm O}(m× n)$ ,是多项式内可以求解. 而搬运模式的规模为 ${\rm O}((m+2) m!)$ ,那么整个问题解空间的规模为 ${\rm O}(n(m^2+2m)m!)$ .所以对于问题的求解分为两个阶段:加工等待时间优化和晶圆搬运顺序优化.然后分别构建全局搜索算法和分枝搜索算法对问题进行求解.

    在集束型设备中,为了方便构建算法,提出以下引理:

    引理 1. 生产单晶圆类型的集束型设备,使用超过两个单元搬运模式一般都不是最优的搬运策略.

    证明. 假设 $T_i^h$ 是在 $s_i$ 搬运模式下的周期时间, $T_i^h$ 表示在 $s_i$ 搬运模式下的周期时间,且有 $T_i^h=T_j^h$ .那么 $2T_i^h$ 就是在 $s_i$ 下生产两个晶圆的时间,而 $2T_j^h$ 为在 $s_j$ 下生产两个晶圆的时间.如果有 $T_i^h\geq T_j^h$ ,则 $2T_i^h\geq T_i^h+T_j^h\geq 2T_j^h$ ;否则,如果 $T_i^h\leq T_j^h$ ,则有 $2T_i^h\leq T_i^h+T_j^h\leq 2T_j^h$ .可知 $T_i^h+T_j^h$ 不可能小于 $\min(2T_i^h,2T_j^h)$ 的,即两种搬运模式的组合一般不是最优策略.对于超过两个周期的情况证法相似.

    定理 1.生产单晶圆类型的集束型设备,最优的搬运策略必定是单一的搬运模式.

    证明. 假设 $\Omega$ 为所有的搬运策略集合, $\Omega$ 由 $k$ 个集合即 $\{\Omega_1,\Omega_2,\cdots,\Omega_k\}$ 组成,其中 $\Omega_k$ 表示 $k$ 种不同搬运模式组成的搬运策略集合,那么必有最优策略 $X\in\Omega$ .假设定理1不成立,那么有 $X\in\Omega$ . 由引理1可知, $X\in\Omega|k\geq2$ ,则 $X\in {\bf C}_\Omega(\Omega_k)$ , ${\bf C}_\Omega(\Omega_k)$ ,即 $X\in\Omega_1$ ,与假设矛盾.所以定理成立.

    由定理1可知,最优搬运模式 $s_i$ 必然存在于有限搬运模式集合 $\{s\}$ .给定的晶圆搬运模式 $s_i$ 里,问题可被描述为前一部分提出的模型.分别用 $u_{ki}$ , $z_{ki}$ 和 $w_{ki}$ 表示在搬运模式 $s_i$ 下的 $u_k$ , $z_k$ 和 $w_k$ .

    单晶圆类型的调度搜索算法描述如下:

    输入. 搬运时间 $\delta$ 与 $\theta$ ,加工晶圆 $W_j$ 的加工时间 $p_k$ 和驻留约束时间 $a_k$ .

    输出. 晶圆 $W_j$ 的最优搬运策略周期时间 $T_j$ .

    步骤 1. 载入相关参数,比如 $m$ , $p_k$ , $a_k$ , $\delta$ 和 $\theta$ . 令 $i=0$ , $k=0$ 和 $T=+\infty$ .

    步骤 2.令 $i=i+1$ ,如果有 $italic>(m+2) m!||italic>G$ ,那么转到步骤7;否则,转到步骤3.

    步骤 3.令 $k=k+1$ ,如果 $k>2m+1$ ,那么转到步骤5;否则,搬运模式为 $s_i$ ,利用式(2) 计算出 $IC(s_i)$ ,利用式(3) 计算出 $z_{ki}$ .

    步骤 4.如果满足 $z_{ki}-p_k>a_k$ ,那么转到步骤2;否则,转到步骤3.

    步骤 5. 如果 $k$ 为偶数,就令 $w_{ki}=0$ .

    步骤 6. 优化的目标函数(4) ,服从约束(5) $\sim$ (8) ,求得所有 $w_{kt}$ 和 $T(s_i)$ .如果 $T(s_i)<T$ ,那么 $T=T(s_i)$ 且 $h=i$ ;否则, $T$ 保持不变.然后转回步骤2.

    步骤 7. $s_h$ 就是晶圆 $W_j$ 的最优搬运策略.

    问题的下界首先被用于提高分枝算法的速度,再被用来验证近似最优解的优劣.该部分松弛驻留约束(7) ,那么原问题就转换为以式(4) 为优化目标函数,并服从约束(5) 、(6) 和(8) 的松弛问题.按照文献 [7-8]的拉式与推式策略结合得到下界估算算法:

    步骤 1. 初始化数据,如果 $p_{\max} {\geq2m(\delta+\theta)}$ ,则 $T=p_{\max}+2(\delta+\theta)$ ; 否则,转到步骤2.

    步骤 2. 令 $i=1$ , $j=1$ , $K$ 取最大值下标.

    步骤 3. $n=i+j$ ,当 $n\leq m$ 时,转步骤4; 否则,当 $j\geq2$ 时转步骤8, $j=1$ 转步骤10.

    步骤 4. 若 $n=K$ ,则转步骤6; 否则,转步骤5.

    步骤 5. 若 $\sum_i^n{p_i}+(n-i)\delta\leq p_{\max}$ ,则 $j=j+1$ ,转到步骤3;否则,若 $\sum_i^{n-1}{p_i}+(n-1-i)\delta<p_{\max}$ ,转到步骤9;否则, $i=i+1$ , $j=1$ ,转步骤3.

    步骤 6. 如果 $\sum_i^{n-1}{p_i}+(n-1-i)\delta\leq p_{\max}$ ,转到步骤7.

    步骤 7. $p_i=\sum_i^{K-1}{p_i}+(K-1-i)\delta$ , $p_{i+k}=p_{K-1+k}$ , $m=m+i-K+1$ , $K=i+1$ , $i=K+1$ , $j=1$ ,转到步骤3.

    步骤 8.若 $\sum_i^{n-1}{p_i}+(n-1-i)\delta<p_{\max}$ ,转到步骤9,否则 $i=i+1$ , $j=1$ ,转步骤3.

    步骤 9. $p_i=\sum_i^{n-1}{p_i}+(n-1-i)\delta$ , $p_{i+k}=p_{n-1+k}$ , $m=m+i-n+1$ , $K=K+i-n+1$ ,若 $p_{\max}\geq2m(\delta+\theta)$ ,则 $T=p_{\max}+2(\delta+\theta)$ ; 否则, $i=i+1$ , $j=1$ ,转到步骤3.

    步骤 10. 如果当 $p_{\max}\geq 2m(\delta+\theta)$ ,那么有 $T=p_{\max}+2(\delta+\theta)$ ; 否则, $T=2(m+1) (\delta+\theta)$ .

    最后算法得到的最优解即为原问题的下界.

    本节对带缓冲的两加工模块的集束型设备详细分析,根据第一节的模型建立数学模型,得到每种搬运模式的周期时间的表达式,最后分析比较各自的周期时间.

    当 $s=s_5$ 时, $\varphi(k)$ 、 $\varphi^-(k)$ 、 $u_{k5}$ 和 $λ_k$ 的值如表 1所示.

    表 1  $s=s_5$ 时相应参数的值
    Table 1  The relevant parameter when $s=s_5$
    $k$ 0 1 2 3
    $\varphi(k)$ 0 2 1 3
    $\varphi^-(k)$ 0 2 1 3
    $u_k5$ $\setminus$ 0 1 0
    ${{\lambda }_{k}}$ 1 1 1 1
    下载: 导出CSV 
    | 显示表格

    那么问题可以转化为目标函数:

    $\min \sum\limits_{i=1}^{3}{({{w}_{i5}})}$

    约束条件为

    ${{w}_{15}}+{{w}_{25}}\ge \max ({{p}_{1}}-{{z}_{15}},0)\text{ }$

    (9)

    ${{w}_{25}}+{{w}_{35}}\ge \max ({{p}_{2}}-{{z}_{25}},0)\text{ }$

    (10)

    ${{w}_{15}}+{{w}_{35}}\ge \max ({{p}_{3}}-{{z}_{35}},0)\text{ }$

    (11)

    ${{z}_{15}}-{{p}_{1}}\le {{a}_{1}}$

    (12)

    ${{z}_{25}}-{{p}_{2}}\le {{a}_{2}}$

    (13)

    ${{z}_{35}}-{{p}_{3}}\le {{a}_{3}}$

    (14)

    解得:

    ${{w}_{15}}=\max ({{p}_{1}}-{{z}_{15}},0)\text{ }$

    (15)

    ${{w}_{25}}=0\text{ }$

    (16)

    ${{w}_{35}}=\max ({{p}_{3}}-{{z}_{35}}-{{w}_{15}},{{p}_{2}}-{{z}_{25}},0)\text{ }$

    (17)

    搬运模式 $s_5$ 的周期时间为 $T_5$ ,由式(4) 可知:

    \begin{align}T_5=IC(s_5) +\sum_{i=1}^3w_{i5} \end{align}

    (18)

    其中

    \begin{align}IC(s_5) =4\delta+4\theta \end{align}

    (19)

    把式(15) $\sim$ (17) 、(19) 带入式(18) 可得:

    \begin{align*}&T_5=4\delta+4\theta+w_{15}+w_{25}+w_{35}=\\&\quad 4\delta+4\theta+w_{15}+\max(p_3-\delta-3\theta-w_{15},0) =\\&\quad 4\delta+4\theta+\max(p_3-\delta-3\theta,w_{15})=\\&\quad 4\delta+4\theta+\max(p_3-\delta-3\theta,\max(p_3-\delta-3\theta,0) )=\\&\quad 3\delta+\theta+\max(p_3,\max(p_1,\delta+3\theta))=\\&\quad 3\delta+\theta+\max(p_3,p_1,\delta+3\theta)\end{align*}

    整理后:

    \begin{align}T_5=3\delta+\theta+\max(p_3,p_1,\delta+3\theta) \end{align}

    (20)

    相似地,能分别得到各种搬运模式 $s_1$ , $s_2$ , $s_3$ , $s_4$ , $s_5$ , $s_6$ , $s_7$ , $s_8$ 的周期时间 $T_1$ , $T_2$ , $T_3$ , $T_4$ , $T_5$ , $T_6$ , $T_7$ , $T_8$ 为

    ${{T}_{1}}=3\delta +{{p}_{1}}+{{p}_{3}}\text{ }$

    (21)

    ${{T}_{2}}=2\delta +\theta +\max ({{p}_{1}},{{p}_{3}},\delta +3\theta )$

    (22)

    ${{T}_{3}}=4\delta +{{p}_{1}}+{{p}_{3}}$

    (23)

    ${{T}_{4}}=2\delta +2\theta +\max ({{p}_{1}}+2\delta +2\theta ,{{p}_{3}})$

    (24)

    ${{T}_{6}}=2\delta +2\theta +\max ({{p}_{1}},{{p}_{3}}+2\delta +2\theta )$

    (25)

    ${{T}_{7}}=3\delta +\theta +\max ({{p}_{1}},{{p}_{3}},\delta +3\theta )$

    (26)

    ${{T}_{8}}=2\delta +2\theta +\max ({{p}_{1}},{{p}_{3}},2\delta +6\theta )$

    (27)

    定义 3. 若搬运模式 $s_i$ 下的周期时间 $T_i$ 大于或者等于搬运模式 $s_j$ 下的周期时间 $T_j$ ,则可说搬运模式 $s_i$ 受支配于搬运模式 $s_j$ ,即搬运模式 $s_j$ 相对于搬运模式 $s_i$ 具有一定优势.

    定理 2. 对于任意的晶圆,在满足驻留约束的前提下,搬运模式 $s_3$ 受支配于搬运模式 $s_1$ ,搬运模式 $s_7$ 受支配于搬运模式 $s_2$ .

    证明. 搬运模式 $s_3$ 与搬运模式 $s_1$ 的搬运动作相似,不同之处仅为搬运模式 $s_3$ 比搬运模式 $s_1$ 多了缓冲模块与加工模块之间的搬运.对周期时间进行分析, $T_1=3\delta+p_1+p_3$ 且 $T_3=4\delta+p_1+p_3$ ,有 $T_1<T_3$ ,则对于任意的晶圆,搬运模式 $s_3$ 受支配于搬运模式 $s_1$ .搬运模式 $s_7$ 里搬运模块到缓冲模块进行装卸,而搬运模式 $s_2$ 的晶圆直接跳过缓冲模块,对周期的表达式分析,有 $T_7=T_2+\delta>T_2$ ,则对于任意的晶圆,在满足时间驻留约束的前提下,搬运模式 $s_7$ 受支配于搬运模式 $s_2$ .

    定理 3. 对于任意晶圆,在满足驻留约束的前提下,搬运模式 $s_4$ , $s_5$ , $s_6$ , $s_8$ 受支配于搬运模式 $s_2$ .

    证明. 将各个搬运模式分别与搬运模式 $s_2$ 进行比较:

    1) 由式(22) 与(24) 可知,明显地,有不等式 $\max(p_1,p_3,\delta+3\theta)\leq\max(p_1+2\delta+3\theta,p_3+\theta)$ 成立,可推出 $T_2\leq T_4$ .

    2) 由式(22) 与(20) 可得, $T_5=T_2+\delta\geq T_2$ .

    3) 推理同1) 类似,由式(22) 与(25) 推得 $T_2\leq T_6$ .

    4) 由式(22) 与(27) 可知,明显地,有不等式 $\max(p_1,p_3,\delta+3\theta)\leq\max(p_1+\theta,p_3+\theta,2\delta+7\theta)$ 成立,可以推出 $T_2\leq T_8$ .

    总之,在满足驻留约束的前提下,搬运模式 $s_4$ , $s_5$ , $s_6$ , $s_8$ 受支配于搬运模式 $s_2$ .

    为更加方便地描述复杂条件,定义了如下的逻辑表达:

    \begin{align*}&A(x)=\{x\leq 4\theta\}\\[1mm]&B(x)=\{x\geq 4\theta\}\\[1mm]&C(x,y)=\{\max(x,y)\leq 2\delta+6\theta\}\\[1mm]&D(x,y)=\{\max(x,y)\geq 2\delta+6\theta\}\\[1mm]&E(x,y)=\{(x-y)\leq 2\delta+6\theta\}\\[1mm]&F(x,y)=\{(x-y)\geq 2\delta+6\theta\}\\[1mm]&g(x,y)=\{[A(y)\cap C(x,y)\cap E(x,y)]\cup\\[1mm]&\qquad [C(x,y)\cap F(x,y)]\}\\[1mm]&h(x,y)=\{D(x,y)\cup[B(y)\cap C(x,y)\cap E(x,y)]\}\end{align*}

    定理 4. 搬运模式之间的周期时间存在如下的大小关系:

    1) $T_1\leq T_2\Leftrightarrow p_1+p_3\leq 4\theta$ 而 $T_1\geq T_2\Leftrightarrow p_1+p_3\geq 4\theta$ .

    2) $T_1\leq T_4\Leftrightarrow p_1\leq 2\theta-\delta\Arrowvert p_3\leq 4\theta+\delta$ 而 $T_1\geq T_4\Leftrightarrow p_1\geq2\theta-\delta~\&~p_3\geq 4\theta+\delta$ .

    3) $T_1\leq T_5\Leftrightarrow p_1\leq\theta\Arrowvert p_3\leq\theta\Arrowvert p_1+p_3\leq\delta+3\theta$ 而 $T_1\geq T_5\Leftrightarrow p_1\geq\theta~\&~ p_3\geq\theta ~\&~p_1+p_3\geq\delta+3\theta$ .

    4) $T_1\leq T_6\Leftrightarrow p_1\leq4\theta+\delta\Arrowvert p_3\leq 2\theta-\delta$ 而 $T_1\geq T_6\Leftrightarrow p_1\geq4\theta+\delta~\&~ p_3\geq 2\theta-\delta$ .

    5) $T_1\leq T_8\Leftrightarrow \max(p_1,p_3) \leq2\theta-\delta\Arrowvert p_1+p_3\leq\delta+8\theta$ 而 $T_1\geq T_8\Leftrightarrow \max(p_1,p_3) \geq 2\theta-\delta~\&~p_1+p_3\geq\delta+8\theta$ .

    6) $T_4\leq T_5\Leftrightarrow p_3-p_1\geq3\theta+\delta$ 而 $T_4\geq T_5\Leftrightarrow p_3-p_1\leq3\theta+\delta$ .

    7) $T_4\leq T_6\Leftrightarrow p_1\leq p_3$ 而 $T_4\leq T_6\Leftrightarrow p_1\leq p_3$ .

    8) $T_4\leq T_7\Leftrightarrow p_3-p_1\geq3\theta+\delta$ 而 $T_4\geq T_7\Leftrightarrow p_3-p_1\leq3\theta+\delta$ .

    9) $T_4\leq T_8\Leftrightarrow g(p_3,p_1) $ 而 $T_4\geq T_8\Leftrightarrow h(p_3,p_1) $ .

    10) $T_5\leq T_6\Leftrightarrow p_1-p_3\leq2\theta+2\delta$ 而 $T_5\geq T_6\Leftrightarrow p_1-p_3\geq2\theta+2\delta$ .

    11) $T_5\leq T_8\Leftrightarrow \max(p_1,p_3) \leq2\theta+6\delta$ 而 $T_5\geq T_8\Leftrightarrow \max(p_1,p_3) \geq2\theta+6\delta$ .

    12) $T_6\leq T_8\Leftrightarrow g(p_1,p_3) $ 而 $T_6\leq T_8\Leftrightarrow h(p_1,p_3) $ .

    证明. 1) 从式(21) 与(22) 可知,若有 $T_1\leq T_2\Rightarrow\delta+p_1+p_3\leq\max(p_1+\theta,p_3+\theta,\delta+4\theta)$ ,由于 $\delta$ 表示晶圆移动时间而 $\theta$ 表示机械手空移动,存在 $\delta\geq\theta$ .由于 $\delta+p_1+p_3\leq\max(p_1+\theta,p_3+\theta)$ 始终成立,则需有 $\delta+p_1+p_3\leq \delta +4\theta\Rightarrow p_1+p_3\leq4\theta$ . 相似地, $T_1\geq T_2$ 需满足条件 $p_1+p_3\geq 4\theta$ .其他2) $\sim$ 12) 都可用类似的方法推理得到.

    根据定理2 $\sim$ 4,结合单晶圆类型调度的全局搜索算法,首先利用驻留约束确定 $s_2$ 的可行性得到两个分枝,其中的第一分枝利用定理2和定理3里的支配关系排除了 $s_4$ , $s_5$ , $s_6$ , $s_7$ , $s_8$ ,第二分枝利用定理4的条件分析再进行分枝.图 2为改进型搜索决策树示例.

    图 2  搜索决策树
    Fig. 2  The decision tree for searching

    为了对带缓冲的集束型设备进行评价,把不带缓冲的集束型设备用于对比分析. 运用量化的分析方法,为带有缓冲的单臂集束型设备相对于无缓冲的集束型设备的产能提高提供了较为可靠的数据.本节中假设 $\delta=4$ , $\theta=2$ , $a_i=5$ 为最基本的实验条件.由于增加缓冲的优势主要体现在最大加工时间的加工模块上,所以有 $p_{\max}=\max_{1\leq i\leq 2m}\{p_i\}$ .用C++分别对两种算法进行编程,统计50组算例数据,求得平均的结果,并对比两种算法的数值实验结果. 图 3为两种算法所需的时间,可以看出分枝搜索的速度会比全局搜索算法更加优越.但图 4表明在小规模时,两种算法都能得到较好的结果.

    图 3  两种算法所需的时间
    Fig. 3  Running time for the two algorithms
    图 4  两种算法的GAP值对比
    Fig. 4  GAP for the two algorithms

    将本文提出的算法与文献[13]的算法进行比较. 取加工时间为 $p_i\sim{\rm N}(60,2.5^2) $ ,如图 5所示为算法之间调度结果的比较.

    图 5  算法结果的比较
    Fig. 5  The comparison between the algorithms

    图 5可知,在小规模问题时,本文的算法与文献的算法得到的解基本相同,但随着规模的增加,本文的算法呈现出较为明显的优势,得到的解更优.

    本文定义了提高率 $R$ 来评价带缓冲模块相对于无缓冲模块的系统之间产量提高:

    $R=\frac{C{{T}_{G}}-C{{T}_{B}}}{C{{T}_{B}}}\times 100\%$

    上式表示的是带缓冲的集束型设备与无缓冲的设备的集束型设备产能的提高率.其中 $R$ 越大,表示提高率越大.其中 $CT_B$ 代表带缓冲的集束型设备的周期时间,而 $CT_G$ 为无缓冲的集束型设备的周期时间. $Bd={p_{\max}}/{(\delta+\theta)}$ 表示搬运模块的繁忙程度, $Bd$ 越大表示搬运块模块越空闲.

    令机械手繁忙程度 $Bd$ 取[0,10]的均匀分布,通过仿真分析,得到如图 6的结果.

    图 6  生产率提高与搬运模块繁忙程度的关系
    Fig. 6  The relationship of improvement rates with the busy degree of TM

    分析图 6可知,对于两集束型而言,仅当 $Bd\geq3$ 时,产能的提高率为正数,即仅当搬运模块繁忙程度 $Bd\geq3$ 时,缓冲模块能提高集束型设备的产能. 相反,当搬运模块繁忙程度 $Bd\leq3$ 时,产能提高率小于零,说明缓冲模块不提高产能,而对于三集束型而言, $Bd$ 的临界值为 $3.1$ .同时随着产能提高率的减少,最大加工时间 $p_{\max}$ 也相应减少,故TM在PM上等待时间也相应减少,直到晶圆加工完毕.若晶圆在PM上的等待搬运的时间大于驻留约束,那么搬运模块的搬运模式必须改变; 否则,晶圆会出现质量问题.无论搬运模块的繁忙程度如何的变化,最后TM都会处于搬运模式 $s_1$ ,说明在缓冲模块上的晶圆增加了额外的搬运.

    设置变量 $P={p_t}/{(\delta+\theta)}$ ,然后令 $Bd$ 分布取均匀分布[3, 9]. 经过仿真分析,得到结果如图 7所示.

    图 7  生产率提高与加工时间之间的关系
    Fig. 7  The relationship of improvement rates with processing time of PMs

    图 7可知,在所有的搬运模块繁忙程度中,产能提高率从不稳定改变到达稳定的状态.主要原因在于随着模块加工时间的增加,搬运模块空闲使得加工模块的驻留约束的影响被降低.

    本文针对带缓冲模块的集束型设备的调度问题建立了数学模型,构建全局搜索算法. 再通过具体分析两集束型设备,提出了分枝搜索算法.最后经仿真表明,分枝搜索算法比全局搜索算法更加迅速.同时通过和文献[13]的调度算法进行对比,表明本文的算法在加工模块比较少的时候调度结果提高不多,但随着规模增大,本文所得到的结果有明显的优势.通过对比带缓冲的集束型设备与无缓冲的集束型设备,证明了在加工模块之间增加缓冲,在一定程度上可以提高单臂集束型设备的产能.基于实验数据,仅当搬运模块TM的繁忙程度 $Bd$ 大于阈值时,设备产能才会提高.实验还表明驻留约束对于产能提高的影响是通过约束限制状态的改变产生的.随着加工模块的加工时间不断增加,产能提高越来越不显著,系统产能渐渐趋于平稳状态.


  • 本文责任编委 孙长银
  • 图  1  FTO算法基本结构

    Fig.  1  Basic structure of FTO

    图  2  斐波那契树生成过程示意图

    Fig.  2  Building procedure of Fibonacci tree

    图  3  FTO算法流程图

    Fig.  3  Algorithm flow chart of FTO

    图  4  函数图像和等高线示意图

    Fig.  4  Function graphs and contour graphs

    图  5  函数Damavandi求解坐标点的累积分布示意图

    Fig.  5  Cumulative distribution of points of function Damavandi

    图  6  函数Langermann求解坐标点的累积分布示意图

    Fig.  6  Cumulative distribution of points of function Langermann

    图  7  仿真实验收敛曲线

    Fig.  7  Convergence curves of simulation experiment

    图  8  对比实验的收敛曲线

    Fig.  8  Convergence curves of comparison experiment

    表  2  算法参数设置

    Table  2  Parameters setting of algorithms

    GA PSO CLPS FTO
    交叉概率 0.7 惯性权重 0.9~0.5 惯性权重 0.9~0.5 斐波那契树深度 6
    变异概率 0.01 加速常数$ c1 $ 2 加速常数$ c1 $ 1.49445
    精英保留概率 0.05 加速常数$ c2 $ 2 加速常数$ c2 $ 1.49445
    二进制位数 20 种群个数 20 学习概率 0.5~0.05
    种群个数 20 种群个数 20
    下载: 导出CSV

    表  1  函数局部极值点

    Table  1  Local optima of benchmark functions

    局部极值1 (全局最优解) 局部极值2
    Damavandi $ f(2, 2)=0 $ $ f(7, 7)=2 $
    Langermann $ f(2.003, 1.006)=-5.1621 $ $ f(7, 9)=-3 $
    下载: 导出CSV

    表  3  测试函数列表

    Table  3  Benchmark functions

    函数名称 表达式 定义域 全局最优解
    Damavandi $ [1-|\frac{\sin[\pi(x_1-2)]\sin[\pi(x_2-2)]}{\pi^2(x_1-2)(x_2-2)}|^5][2+(x_1-7)^2+(x_2-7)^2]- $ $ \sum\limits_{i=1}^5 c_i\cos(z_i\pi){\rm e}^{-\frac{z_i}{\pi}} $ $[0, 14]$ 0
    $ z_i=(x_1-a_i)^2+(x_2-b_i)^2 $
    Langermann $ a=[3, 5, 2, 1, 7] $ $[0, 10]$ $-5.1621$
    $ b=[5, 2, 1, 4, 9] $
    $ c=[1, 2, 5, 2, 3] $
    Michalewicz $ \sum\limits_{i=1}^d \sin(\frac{ix_i^2}{\pi})^{20}\sin(x_i) $ $[0, 5] $ $ -1.9679$
    YangStandingWave $ {\rm e}^{-\sum\limits_{i=1}^d(\frac{x_i}{15})^6}-2{\rm e}^{-\sum\limits_{i=1}^d x_i^2}\prod\limits_{i=1}^d \cos^2(x_i)$ $[-20, 20]$ $ -1 $
    cec2005 15 HCF1 $[-5, 5]$ 120
    cec2005 16 RHCF1 Hybrid composition functions $[-5, 5]$ 120
    cec2005 18 RHCF2 $[-5, 5]$ 10
    cec2005 21 RHCF3 $[-5, 5]$ 360
    下载: 导出CSV

    表  4  实验函数的最优解及可接受范围

    Table  4  The optimal solution and acceptable range of the functions in experiment

    函数 最优解 可接受范围
    Damavandi 0.0000E+00 $ < 1.0000$ E$-$ 02
    Langermann $-5.1621$ E+00 $ < -5.1000$ E+00
    Michalewicz $-1.9679$ E+00 $ < -1.9500$ E+00
    YangStandingWave $-1.0000$ E+00 $ < 0.0000$ E+00
    HCF1 1.2000E+02 $ < 1.2900$ E+02
    RHCF1 1.2000E+02 $ < 1.2900$ E+02
    RHCF2 1.0000E+01 $ < 1.1000$ E+01
    RHCF3 3.6000E+02 $ < 3.6100$ E+02
    下载: 导出CSV

    表  5  Damavandi实验结果

    Table  5  Result of Damavandi in experiment

    算法 到达均值 标准差1 未达均值 标准差2 $N$ $r$ (%)
    GA / / 2.0000E+00 /
    PSO -8.5487E-14 1.3213E-23 2.0000E+00 0.0000E+00 2 4
    CLPSO / / 2.0000E+00 /
    FTO 1.1168E-06 1.4894E-06 2.0000E+00 1.7888E-09 12 24
    下载: 导出CSV

    表  6  Langermann实验结果

    Table  6  Result of Langermann in experiment

    算法 到达均值 标准差1 未达均值 标准差2 $N$ $r$ (%)
    GA $-5.1621$ E+00 1.8078E-15 $-2.9938$ E+00 3.6444E-01 29 58
    PSO $-5.1621$ E+00 2.6968E-15 $-3.0000$ E+00 4.7475E-16 42 84
    CLPSO $-5.1615$ E+00 4.3819E-03 / / 50 100
    FTO $-5.1621$ E+00 5.0777E-08 / / 50 100
    下载: 导出CSV

    表  7  Michalewicz实验结果

    Table  7  Result of Michalewicz in experiment

    算法 到达均值 标准差1 未达均值 标准差2 $N$ $r$ (%)
    GA $-1.9679$ E+00 4.5325E-16 -1.7979E+00 9.6821E-02 25 50
    PSO $-1.9678$ E+00 3.9389E-05 -9.9196E-01 1.6075E-02 46 92
    CLPSO $-1.9679$ E+00 7.4886E-09 / / 50 100
    FTO $-1.9679$ E+00 7.0607E-07 / / 50 100
    下载: 导出CSV

    表  8  YangStandingWave实验结果

    Table  8  Result of YangStandingWave in experiment

    算法 到达均值 标准差1 未达均值 标准差2 $N$ $r$ (%)
    GA $-1.0000$ E+00 5.6277E-16 1.3173E-05 3.5265E-21 37 74
    PSO $-1.0000$ E+00 0.0000E+00 1.5571E-05 1.2516E-06 4 8
    CLPSO $-7.0208$ E-01 2.6415E-01 1.2852E-05 2.0573E-06 9 18
    FTO $-1.0000$ E+00 1.8775E-09 / / 50 100
    下载: 导出CSV

    表  9  HCF1实验结果

    Table  9  Result of HCF1 in experiment

    算法 到达均值 标准差1 未达均值 标准差2 $N$ $r$ (%)
    GA 1.2000E+02 3.6311E-06 2.6457E+02 1.2548E+02 27 54
    PSO 1.2000E+02 0.0000E+00 3.2000E+02 0.0000E+00 45 90
    CLPSO 1.2008E+02 2.7377E-01 1.4330E+02 1.4567E+01 48 96
    FTO 1.2000E+02 1.3480E-05 / / 50 100
    下载: 导出CSV

    表  10  RHCF1实验结果

    Table  10  Result of RHCF1 in experiment

    算法 到达均值 标准差1 未达均值 标准差2 $N$ $r$ (%)
    GA 1.2000E+02 0.0000E+00 2.1146E+02 5.0191E+01 7 14
    PSO 1.2000E+02 0.0000E+00 2.1597E+02 1.5598E+01 35 70
    CLPSO 1.2155E+02 2.6034E+00 1.9118E+02 $-3.3710$ E+01 41 82
    FTO 1.2000E+02 2.3951E-05 2.2000E+02 0.0000E+00 49 98
    下载: 导出CSV

    表  11  RHCF2实验结果

    Table  11  Result of RHCF2 in experiment

    算法 到达均值 标准差1 未达均值 标准差2 $N$ $r$ (%)
    GA 1.0006E+01 0.0000E+00 3.8551E+02 9.0210E+01 1 2
    PSO 1.0000E+01 0.0000E+00 3.8500E+02 6.0356E+01 14 28
    CLPSO / / 2.6217E+02 /
    FTO 1.0043E+01 3.1468E-02 1.1011E+02 8.2284E-02 44 88
    下载: 导出CSV

    表  12  RHCF3实验结果

    Table  12  Result of RHCF3 in experiment

    算法 到达均值 标准差1 未达均值 标准差2 $N$ $r$ (%)
    GA 3.6000E+02 0.0000E+00 7.8881E+02 1.6723E+02 1 2
    PSO / / 6.9800E+02 /
    CLPSO / / 6.0613E+02 /
    FTO 3.6000E+02 5.1019E-04 5.5474E+02 2.2942E+01 31 62
    下载: 导出CSV
  • [1] Boussad I, Lepagnot J, Siarry P. A survey on optimization metaheuristics. Information Sciences, 2013, 237:82-117 doi: 10.1016/j.ins.2013.02.041
    [2] 李建平, 宫耀华, 赵思远, 卢爱平, 李盼池.粒子群优化算法的改进及数值仿真.吉林大学学报(理学版), 2017, 55(2):322-332 http://d.old.wanfangdata.com.cn/Periodical/jldxzrkxxb201702024

    Li Jian-Ping, Gong Yao-Hua, Zhao Si-Yuan, Lu Ai-Ping, Li Pan-Chi. Improvement of particle swarm optimization algorithm and numerical simulation. Journal of Jilin University (Science Edition), 2017, 55 (2):322-332 http://d.old.wanfangdata.com.cn/Periodical/jldxzrkxxb201702024
    [3] 高雷阜, 佟盼.遗传-人工蜂群融合算法及其Markov收敛性分析.数学杂志, 2017, 37 (1):215-222 doi: 10.3969/j.issn.0255-7797.2017.01.023

    Gao Lei-Fu, Tong Pan. Combination algorithm of genetic-artificial bee colony and its Markov convergence analysis. Journal of Mathematics, 2017, 37 (1):215-222 doi: 10.3969/j.issn.0255-7797.2017.01.023
    [4] Whitley D. A genetic algorithm tutorial. Statistics and Computing, 1994, 4 (2):65-85 http://d.old.wanfangdata.com.cn/OAPaper/oai_doaj-articles_c227a5f223e4c1bc4e93b3863a92a817
    [5] Mareda T, Gaudard L, Romerio F. A parametric genetic algorithm approach to assess complementary options of large scale wind-solar coupling. IEEE/CAA Journal of Automatica Sinica, 2017, 4 (2), 260-272 doi: 10.1109/JAS.2017.7510523
    [6] Kennedy J, Eberhart R. Particle swarm optimization. In: Proceedings of the 1995 IEEE International Conference on Neural Networks. Perth, WA, Australia: IEEE, 1995, 4: 1942-1948
    [7] Poli R. The sampling distribution of particle swarm optimisers and their stability[Online], available: http://cswww.essex.ac.uk/technical-reports/2007/csm-465.pdf, 2007.
    [8] Zhang Y D, Wang S H, Ji G L. A comprehensive survey on particle swarm optimization algorithm and its applications. Mathematical Problems in Engineering, 2015, 2015:Article No. 931256 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=Doaj000004037005
    [9] Han Z, Zhao J, Wang W. An optimized oxygen system scheduling with electricity cost consideration in steel industry. IEEE/CAA Journal of Automatica Sinica, 2017, 4 (2), 216-222 doi: 10.1109/JAS.2017.7510439
    [10] 张军英, 许进, 保铮.遗传交叉运算的可达性研究.自动化学报, 2002, 28 (1):120-125 http://www.aas.net.cn/CN/abstract/abstract15509.shtml

    Zhang Jun-Ying, Xu Jin, Bao Zheng. Attainability of genetic crossover operator. Acta Automatica Sinica, 2002, 28 (1):120-125 http://www.aas.net.cn/CN/abstract/abstract15509.shtml
    [11] Zhang Y A, Ma Q L, Sakamoto M, Furutani H. Effect of mutation to distribution of optimum solution in genetic algorithm. Natural Computing. Tokyo, Japan:Springer, 2010. 380-387
    [12] Liu Q F, Wei W H, Yuan H Q, Zhan Z H, Li Y. Topology selection for particle swarm optimization. Information Sciences, 2016, 363:154-173 doi: 10.1016/j.ins.2016.04.050
    [13] Clerc M, Kennedy J. The particle swarm-explosion, stability, and convergence in a multidimensional complex space. IEEE Transactions on Evolutionary Computation, 2002, 6 (1):58-73 doi: 10.1109/4235.985692
    [14] 刘建华, 刘国买, 杨荣华, 胡文瑜.粒子群算法的交互性与随机性分析.自动化学报, 2012, 38 (9):1471-1484 http://www.aas.net.cn/CN/abstract/abstract17757.shtml

    Liu Jian-Hua, Liu Guo-Mai, Yang Rong-Hua, Hu Wen-Yu. Analysis of interactivity and randomness in particle swarm optimization. Acta Automatica Sinica, 2012, 38 (9):1471-1484 http://www.aas.net.cn/CN/abstract/abstract17757.shtml
    [15] Leboucher C, Shin H S, Siarry P, Le Ménec S, Chelouah R, Tsourdos A. Convergence proof of an enhanced particle swarm optimisation method integrated with evolutionary game theory. Information Sciences, 2016, 346-347:389-411 doi: 10.1016/j.ins.2016.01.011
    [16] Subasi M, Yildirim N, Yildiz B. An improvement on Fibonacci search method in optimization theory. Applied Mathematics and Computation, 2004, 147 (3):893-901 doi: 10.1016/S0096-3003(02)00828-7
    [17] 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 doi: 10.1109/TEVC.2005.857610
    [18] Suganthan P N, Hansen N, Liang J J, Deb K, Chen Y P, Auger A, Tiwari S. Problem definitions and evaluation criteria for the CEC 2005 special session on real-parameter optimization. Singapore: Nanyang Technological University[Online], available: http://web.mysites.ntu.edu.sg/epn-sugan/PublicSite/Shared
  • 期刊类型引用(6)

    1. 黄年昌,杨阳,张强,韩军功. 基于深度学习的RGB-D图像显著性目标检测前沿进展. 计算机学报. 2025(02): 284-316 . 百度学术
    2. 闫梦凯,钱建军,杨健. 弱对齐的跨光谱人脸检测. 自动化学报. 2023(01): 135-147 . 本站查看
    3. 王利锋,辛丽平,刘家硕,鞠莲. 基于热红外视频图像监测的海面溢油识别技术研究. 海洋学报. 2022(05): 148-160 . 百度学术
    4. 蔺素珍,张海松,禄晓飞,李大威,李毅. RBNSM:一种复杂背景下红外弱小目标检测新方法. 红外技术. 2022(07): 667-675 . 百度学术
    5. 刘晓玲,牛海春,宋海燕,秦富贞. 复杂环境下弱信号中的红外小目标自动检测. 激光杂志. 2020(10): 82-86 . 百度学术
    6. 刘松涛,姜康辉,刘振兴. 基于区域协方差和目标度的航空侦察图像舰船目标检测. 系统工程与电子技术. 2019(05): 972-980 . 百度学术

    其他类型引用(2)

  • 加载中
  • 图(8) / 表(12)
    计量
    • 文章访问数:  2372
    • HTML全文浏览量:  519
    • PDF下载量:  669
    • 被引次数: 8
    出版历程
    • 收稿日期:  2016-10-09
    • 录用日期:  2017-08-08
    • 刊出日期:  2018-09-20

    目录

    /

    返回文章
    返回