CN103116324A - 基于指标预测和在线学习的微电子生产线调度方法 - Google Patents
基于指标预测和在线学习的微电子生产线调度方法 Download PDFInfo
- Publication number
- CN103116324A CN103116324A CN2012105437877A CN201210543787A CN103116324A CN 103116324 A CN103116324 A CN 103116324A CN 2012105437877 A CN2012105437877 A CN 2012105437877A CN 201210543787 A CN201210543787 A CN 201210543787A CN 103116324 A CN103116324 A CN 103116324A
- Authority
- CN
- China
- Prior art keywords
- scheduling
- fuzzy
- current generation
- gene
- individual
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
基于指标预测和在线学***均流经时间和最大化瓶颈机器利用率为调度目标的微电子生产线可产生较好的调度效果。
Description
技术领域
本发明属于先进制造、自动化和信息领域,具体涉及具有可重入特征的微电子生产线的调度方法。
背景技术
针对较大规模具有可重入特征的微电子生产线调度问题,目前有如下几类工件排序方法:基于模型的方法、智能优化方法、基于启发式规则的方法。基于模型的排序方法需要对调度问题附加一些较强且不切实际的假设条件,因而该类方法主要用来对一些典型且规模较小的微电子生产线调度问题进行性能分析;由于研究的问题一般具有较大的规模,基于计算智能的方法也难以直接应用;基于启发式规则的排序方法实施简单,但该类方法考虑的因素少,对规模较大的微电子生产线调度问题一般很难获得令人满意的解;
发明内容
为解决以最小化平均流经时间和最大化瓶颈机器利用率为调度目标的较大规模微电子生产线调度问题,本发明提出一种基于指标预测和在线学***均流经时间和最大化瓶颈机器利用率为调度目标的微电子生产线可产生较好的调度效果。
1.基于指标预测和在线学***均流经时间和最大化瓶颈机器利用率为调度目标的较大规模微电子生产线的一种 调度方法,且所述方法是在计算机上依次按如下步骤实现的:
步骤(1):初始化,设定如下基本变量
设定问题变量:
N:微电子生产线中的工序总数
Cn:第n个工序的机器组总数,n=1,2,…,N
μnl:第n个工序第l个机器组包含的机器总数,n=1,2,…,N;l=1,2,…,Cn
L:产品类型总数
Vl:属于第l类产品的lot总数,l=1,2,…,L
Jq:第q个lot,q=1,2,…,Q
Gq:lot Jq所包含的操作总数
Uqg:操作Oqg的可加工机器组集合
Tqg:操作Oqg的可加工机器组总数,即Tqg=||Uqg||
bnc:第n个工序第c个机器组缓存
设定算法相关参数:
W:用于迭代式分解的时间窗口长度
κ:每个当前阶段调度子问题需固定加工开始时间和加工机器的操作与当前阶段调度子问题包含操作的比例
M:用于多模糊规则MFRs学习的训练数据规模
M1:每次进行多模糊规则在线学习时需增加和移除的数据规模
minFS:模糊规则的最小支撑度
步骤(2):采集包括上述工序数、品种数、各Lot的工艺路径和所包含的lot数、各工序的机器组数、各机器组中的机器个数、各操作在各机器组上的加工时间在内的调度相关信息并存储至调度数据库中,并形成待求解的微电子生产线调度问题实例。
步骤(3):微电子生产线调度问题迭代分解
在当前调度时刻,基于时间窗口长度W及各lot的工艺路径按如下方法将整个待调度问题的操作分解为已调度操作集B、当前阶段调度子问题的操作集H和当前阶段剩余调度子问题对应的操作集R。
步骤(3.1):确定已调度操作集B
所有已被确定加工机器及加工开始时间的操作构成操作集B。
步骤(3.2):确定当前阶段调度子问题的操作集H
设hq为lotJq在操作集B中的最大操作号,即
hq=max{g|Oqg∈B}q=1,2,…,Q
操作集H根据下式确定
其中,
步骤(3.3):确定用于预测的当前阶段剩余调度子问题对应的操作集R
操作集R由调度问题中不属于操作集B或H的操作组成,即
步骤(4):采用基于资源冲突程度特征值的自适应差分进化方法SDEH求解当前阶段多目标调度子问题
针对上述微电子生产线调度问题,采用如下SDEH算法求解当前阶段多目标调度子问题,形成当前阶段全局调度解;
步骤(4.1):个体编码和种群初始化
个体由在负载较大的前几个机器上加工的若干个操作串组成,每个操作串对应一个基因串,即个体由若干个基因串组成。
采用启发式规则:最短剩余加工时间优先(Shortest Remaining Processing Time,SRPT)、最先进入队列的优先(First Come First Served,FCFS),最短加工时间优先(Shortest ProcessingTime,SPT)、最长剩余加工时间优先(Maximum Remaining Processing Time,MRPT)及上述规则的随机加权规则,求解当前阶段调度子问题,根据下述方法获得负载较大的前几个机器,并形成由Z个个体组成的初始种群P:
设Ψ为预先给定的个体基因长度的最大值,为当前迭代阶段第n个工序第c个机器组缓存bnc中等待加工操作的总数,是当前调度子问题中在第n个工序第c个机器组上加工操作的总数;根据的大小对bnc进行排序,,设排列结果为bη1,bη2,…,bηΓ,其中,Γ为缓存的总数;
然后,令
则个体由当前调度子问题中在与缓存bη1,bη2,…,bηΓ对应的机器组上加工的操作串组成;每个基因串中的基因表示对应操作在相应机器组上的加工优先级,越靠前其加工优先级越高。
步骤(4.2)初始化
设总迭代步数为Λ,令θ=1,初始复制概率为Rep。
步骤(4.3):交叉
采用顺序交叉方法对随机配对的Z/2个体对进行交叉操作,形成子代种群P1,其中,对配对的两个个体的对应基因串单独进行顺序交叉。
步骤(4.4):变异
基于变异概率λ,对需变异的父代个体基因串S,从当前种群随机选择3个不同个体,与S对应的基因串分别设为S1、S2和S3,产生随机整数r∈{1,2,...,m},m为基因串S的基因长度;然后,根据下述流程生成子代个体对应的基因串S′:
步骤(4.4.1):生成取值在区间[0,1]内的m个随机数。
步骤(4.4.2):若上述第j(j=1,2,...,m)个随机数不小于基因复制概率Repij,且j≠r,则复制基因串S中的第j个基因作为子代个体相应基因串S′的第j个基因,设被复制基因的总数为Υ,将所有Υ个被复制基因从S删除,同时删除S1、S2和S3中对应位置的基因,那么基因串S、S1、S2和S3的长度均为m-Υ;Repij为对应基因的复制概率,根据下式自适应调整:
Repij=Rep×βij
其中,βij为相应基因通过预测获得的资源冲突程度值,其定义为:设与基因xij对应的操作为Ouv,Δ为通过求解前一迭代阶段剩余调度子问题获得的调度问题解的总数,其中在δuv个解中,操作Ouv与其它操作间存在资源冲突,则
βij=δuv/Δ
步骤(4.4.3):设迭代次数i=1。
步骤(4.4.6):若i<m-Υ,则令i=i+1,执行步骤(4.4.4);否则,执行步骤(4.5)。
步骤(4.5):轮盘赌选择,其中,对当前调度子问题中的每个个体,基于预测机制,针对当前阶段剩余调度子问题,采用所提出的MFRs对当前阶段剩余调度子问题求解,获得每个个体对应的全局调度目标函数值。
步骤(4.6):迭代终止条件判别
若θ≤Λ,则执行第5.3步,否则执行步骤(4.7);
步骤(4.7):构成当前阶段的全局调度解
由第四步所得到的当前阶段调度子问题的优化解F1和当前阶段剩余调度子问题相应的优化解F2共同构成当前阶段全局调度解F=F1∪F2。
第五步:固定当前阶段调度子问题中部分操作的加工机器和加工开始时间
步骤(5.1):若当前阶段全局调度解F优于迄今已获得的全局优化解Fo,则令Fo=F。
步骤(5.2):根据迄今已获得的全局优化解Fo和预先给定的待固定操作的比例κ,基于时间轴依次固定当前阶段调度子问题中部分操作的加工机器和加工开始时间。
步骤(6):采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习
在当前阶段调度子问题求解完成后,按下述方法利用相关调度数据采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习。
步骤(6.1):生成训练数据
按下述方法,根据当前阶段调度子问题具有最好全局调度性能的解F1生成用于MFRs自适应调整的M1个数据对,若当前阶段为第一阶段,则生成M个数据对;第m个数据对具有如下形式:
步骤(6.1.2):对上述M个决策时刻两个操作的每个个属性值按如下方法进行归一化处理,获得用于MFRs自适应调整的数据对:
m=1,2,…,M;k=1,2,…,K;j=1,2
步骤(6.2):若为第一阶段,则根据生成的M个数据对生成有效的MFRs;否则根据下述提出的调整方法TILM对用于求解剩余调度子问题的MFRs进行在线调整,其中,MFRs中每条规则具有如下形式:
…,
上述模糊规则具有多个前件和一个后件,是k个条件属性,k=1,2,…,K,为的第ikp个语言变量值,Yq是属性的语言变量值总数;当p=1,对应于需要排序的两个操作中的第一个操作;当p=2,对应于第二个;是由同一属性两个操作的不同取值对应的2维模糊网格。D是结论属性,表示两个操作那个具有更高的调度优先级,从而,分类标签只有两个取值,iR=1表示第一个操作具有更高的调度优先级,iR=2表示第二个操作具有更高的调度优先级;CF(R)为规则R的确信度。
步骤(6.3):初始化
计算min FSa
min FSa=(min FS×M-M1)/M
步骤(6.4):数据移除
按下述方法更新保留的模糊网格。
步骤(6.4.1):从M个训练数据中移除最先生成的M1个训练数据。
步骤(6.4.2):对由M个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格, 计算其在剩余M-M1个训练数据下的模糊支撑度FSr。
步骤(6.4.3):若FSr≥minFSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.5):数据新增
按下述方法更新保留的模糊网格。
步骤(6.5.1):将求解当前阶段调度子问题过程中新生成的M1个训练数据增加到剩余的训练数据中。
步骤(6.5.2):对由剩余的M-M1个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格,计算其在加入新生成数据后的M个训练数据下的模糊支撑度FSa。
步骤(6.5.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.6):增加新的模糊网格
根据新生成的M1个训练数据,增加新的模糊网格。
步骤(6.6.1):根据新生成的M1个训练数据,生成模糊支撑度大于或等于min FS的模糊网格。
步骤(6.6.2):对于生成的模糊网格,基于M个训练数据,按步骤(6.5.2)计算其模糊支撑度FSa。
步骤(6.6.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.7):获得有效的MFRs
将所保留的模糊网格生成对应的多模糊规则MFRs。
步骤(6.7.1):根据上述保留的模糊网格获得有效的MFRs;其中,每条模糊分类规则R的模糊支撑度FS(R)作为其确信度,即:
CF(R)=FS(R)
步骤(6.7.2):考虑到两个操作在输入顺序改变的情形下,其调度优先级应该是不变的,从而,在多模糊规则中,每条规则R均生成一条逆规则R′与其对应,逆规则具有如下形式:
…,
其中,CF(R′)为规则R′的确信度,本发明中令CF(R′)=CF(R)
步骤(7):判定算法是否结束
若待调度问题的所有操作均已调度完成,则结束算法;否则,执行步骤(3)。
根据上述基于指标预测和在线学***均流经时间。
附图说明
图1:基于指标预测和在线学习的微电子生产线调度硬件***结构示意图。
图2:基于指标预测和在线学习的微电子生产线调度方法流程示意图。
图3:多模糊规则在线学习方法结构示意图。
图4:传统差分进化方法(DE)和DPL算法的平均流经时间优化曲线。
具体实施方式
本发明调度方法依赖于相关数据采集***,有调度***客户端和调度服务器实现。在实际微电子生产线调度中应用本发明的软硬件架构示意图如图1所示,本发明的实施方式如下。步骤(1):采集上述微电子生产线工序数、品种数、各品种的工艺路径和所包含的lot数、各工序的机器组数、各机器组中的机器个数、各操作在各机器组上的加工时间在内的调度相关信息并存储至调度数据库中,并形成待求解的微电子生产线调度问题实例。
步骤(2):设定算法相关参数
W:用于迭代式分解的时间窗口长度;
κ:每个当前阶段调度子问题需固定加工开始时间和加工机器的操作与当前阶段调度子问题包含操作的比例;
M:用于多模糊规则MFRs学习的训练数据规模;
M1:每次进行多模糊规则在线学习时需增加和移除的数据规模;
min FS:模糊规则的最小支撑度。
步骤(3):微电子生产线调度问题迭代分解
在当前调度时刻,基于时间窗口长度W及各lot的工艺路径按如下方法将整个待调度问题的操作分解为已调度操作集B、当前阶段调度子问题的操作集H和当前阶段剩余调度子问题对应的操作集R,具体过程可参加“发明内容”相关部分。
步骤(4):采用基于资源冲突程度特征值的自适应差分进化方法SDEH求解当前阶段多目标调度子问题
针对上述微电子生产线调度问题,采用如下SDEH算法求解当前阶段多目标调度子问题,形成当前阶段全局调度解。
步骤(4.1):个体编码和种群初始化
个体由在负载较大的前几个机器上加工的若干个操作串组成,每个操作串对应一个基因串,即个体由若干个基因串组成。
采用SRPT、FCFS、SPT和MRPT启发式规则及上述规则的随机加权规则,求解当前阶段调度子问题,形成由Z个个体组成的初始种群P,其中,其负载最大的机器由所提出的相 应方法确定。上述规则的含义见“发明内容”个体由在负载较大的前几个机器上加工的若干个操作串组成,每个操作串对应一个基因串,即个体由若干个基因串组成。
步骤(4.2)初始化
设总迭代步数为Λ,令迭代次数θ=1,初始复制概率为Rep。
步骤(4.3):交叉
采用顺序交叉方法对随机配对的Z/2个体对进行交叉操作,形成子代种群P1,其中,对配对的两个个体的对应基因串单独进行顺序交叉。
步骤(4.4):变异
基于变异概率λ,对需变异的父代个体基因串S,从当前种群随机选择3个不同个体,与S对应的基因串分别设为S1、S2和S3,产生随机整数r∈{1,2,...,m},m为基因串S的基因长度;然后,根据所提出方法生成子代个体对应的基因串S′。
步骤(4.5):轮盘赌选择,其中,对当前调度子问题中的每个个体,基于预测机制,针对当前阶段剩余调度子问题,采用所提出的MFRs对其对应的全局调度目标函数值进行预测。
步骤(4.6):迭代终止条件判别
若θ≤Λ,则执行第5.3步,否则执行步骤(4.7)。
步骤(4.7):构成当前阶段的全局调度解
由第四步所得到的当前阶段调度子问题的优化解F1和当前阶段剩余调度子问题相应的优化解F2共同构成当前阶段全局调度解F=F1∪F2;
第五步:固定当前阶段调度子问题中部分操作的加工机器和加工开始时间
步骤(5.1):若当前阶段全局调度解F优于迄今已获得的全局优化解Fo,则令Fo=F。
步骤(5.2):根据迄今已获得的全局优化解Fo和预先给定的待固定操作的比例κ,基于时间轴依次固定当前阶段调度子问题中部分操作的加工机器和加工开始时间。
步骤(6):采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习
在当前阶段调度子问题求解完成后,按下述方法利用相关调度数据采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习。
步骤(6.1):生成训练数据
按“发明内容”中方法,根据当前阶段调度子问题具有最好全局调度性能的解F1生成用于MFRs自适应调整的M1个数据对,若当前阶段为第一阶段,则生成M个数据对。
步骤(6.2):若为第一阶段,则根据生成的M个数据对生成有效的MFRs;否则根据”发明 内容”提出的调整方法TILM对用于求解剩余调度子问题的MFRs进行在线调整,其中,MFRs中每条规则具有“发明内容”中给出的形式。
步骤(6.3):初始化
计算min FSa
min FSa=(minFS×M-M1)/M
步骤(6.4):数据移除
按下述方法更新保留的模糊网格。
步骤(6.4.1):从M个训练数据中移除最先生成的M1个训练数据。
步骤(6.4.2):对由M个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格,
计算其在剩余M-M1个训练数据下的模糊支撑度FSr,详见“发明内容”。
步骤(6.4.3):若FSr≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.5):数据新增
按下述方法更新保留的模糊网格。
步骤(6.5.1):将求解当前阶段调度子问题过程中新生成的M1个训练数据增加到剩余的训练数据中。
步骤(6.5.2):对由剩余的M-M1个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格,计算其在加入新生成数据后的M个训练数据下的模糊支撑度FSa,详见“发明内容”。
步骤(6.5.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.6):增加新的模糊网格
根据新生成的M1个训练数据,增加新的模糊网格。
步骤(6.6.1):根据新生成的M1个训练数据,生成模糊支撑度大于或等于minFS的模糊网格。
步骤(6.6.2):对于生成的模糊网格,基于M个训练数据,按步骤(6.5.2)计算其模糊支撑度FSa。
步骤(6.6.3):若FSa≥minFSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.7):获得有效的MFRs
将所保留的模糊网格生成对应的多模糊规则MFRs,详见“发明内容”.
步骤(7):判定算法是否结束
若待调度问题的所有操作均已调度完成,则结束算法;否则,执行步骤(3);
基于指标预测和在线学习的微电子生产线调度方法的参数选择如下:
●用于迭代分解过程中的算法参数:
时间窗口长度W:设为所加工产品中所有操作平均加工时间的10倍与最大加工时间的最大值;
每个调度子问题需要固定加工开始时间和加工机器的操作的比例κ=80%。
●用于多模糊规则在线学习的算法参数:
用于多模糊规则MFRs学习的训练数据规模M=200;
每次训练数据增加和移除的规模M1=20;
模糊规则的最小支撑度min FS=0.4;
考虑到上述多模糊规则的学习过程需要较长的计算时间,本算法在每隔5个迭代阶段学习一次。
●用于求解当前阶段调度子问题的自适应差分进化算法参数:
种群规模Z=20;
总迭代步数为Λ=20;
初始复制概率Rep=0.5;
变异概率λ=0.2。
根据上述所提出的基于指标预测和在线学习的微电子生产线调度方法,本发明做了大量的仿真试验,运行的硬件环境为:P4 2.8GHz CPU,1024M RAM,操作***为Windows、UNIX。
由于篇幅所限,仅列出部分实验结果,表1给出了3个产品在相应机器组上的加工时间,其中,∞表示相应机器组具有足够的加工能力,“-”表示对应产品不需要在相应机器组上加工;表2-4分别给出了3个产品的工艺路径:
表1工艺路径参数
表2产品1的工艺路径
表3产品2的工艺路径
表4产品3的工艺路径
实验一:传统差分进化算法(DE)与DPL算法的比较:表5给出了DE与DPL的数值比较结果。其中,P表示四个不同的问题实例集,每个集合包括包括25个问题实例,4个集合内的实例分别具有40,50,60和70个lot,但不同的实例为上述3种产品的不同组合。Ct表示相应算法的平均计算时间,Obj1和Obj2分别表示平均流经时间和瓶颈机器利用率。Ip表示DPL相对DE的最大平均改进百分比。
针对实验一给出的问题实例,图4给出了DE和DPL算法的平均流经时间优化曲线。从中可以看出,DPL优于DE。
表5DE与DPL性能比较
实验二:DPL与启发式规则的比较;针对与实验二相同的问题实例集,表6给出了DPL与各种启发式规则的性能比较,从中可看出,在平均流经时间上DPL均优于启发式规则。
表6DPL与各种启发式规则的性能比较
实验三:基因复制概率自适应调整对DPL算法的影响:表7给出了DPL与无自适应调整的DPL(简称为DPLO)的性能比较。从中可以看出基因复制概率自适应调整对DPL算法有一定的影响。
表7DPL与DPLO的性能比较
Claims (1)
1.基于指标预测和在线学***均流经时间和最大化瓶颈机器利用率为调度目标的较大规模微电子生产线的一种调度方法,且所述方法是在计算机上依次按如下步骤实现的:
步骤(1):初始化,设定如下基本变量
设定问题变量:
N:微电子生产线中的工序总数
Cn:第n个工序的机器组总数,n=1,2,…,N
μnl:第n个工序第l个机器组包含的机器总数,n=1,2,…,N;l=1,2,…,Cn
L:产品类型总数
Vl:属于第l类产品的lot总数,l=1,2,…,L
Jq:第q个lot,q=1,2,…,Q
Gq:lot Jq所包含的操作总数
Oqg:lot Jq的第g个操作,g=1,2,…,Gq
Uqg:操作Oqg的可加工机器组集合
Tqg:操作Oqg的可加工机器组总数,即Tqg=||Uqg||
bnc:第n个工序第c个机器组缓存
设定算法相关参数:
W:用于迭代式分解的时间窗口长度
κ:每个当前阶段调度子问题需固定加工开始时间和加工机器的操作与当前阶段调度子问题包含操作的比例
M:用于多模糊规则MFRs学习的训练数据规模
M1:每次进行多模糊规则在线学习时需增加和移除的数据规模
min FS:模糊规则的最小支撑度
步骤(2):采集包括上述工序数、品种数、各Lot的工艺路径和所包含的lot数、各工序的机器组数、各机器组中的机器个数、各操作在各机器组上的加工时间在内的调度相关信息并存储至调度数据库中,并形成待求解的微电子生产线调度问题实例;
步骤(3):微电子生产线调度问题迭代分解
在当前调度时刻,基于时间窗口长度W及各lot的工艺路径按如下方法将整个待调度问题 的操作分解为已调度操作集B、当前阶段调度子问题的操作集H和当前阶段剩余调度子问题对应的操作集R;
步骤(3.1):确定已调度操作集B
所有已被确定加工机器及加工开始时间的操作构成操作集B;
步骤(3.2):确定当前阶段调度子问题的操作集H
设hq为lot Jq在操作集B中的最大操作号,即
hq=max{g|Oqg∈B}q=1,2,…,Q
操作集H根据下式确定
其中,
步骤(3.3):确定用于预测的当前阶段剩余调度子问题对应的操作集R
操作集R由调度问题中不属于操作集B或H的操作组成,即
步骤(4):采用基于资源冲突程度特征值的自适应差分进化方法SDEH求解当前阶段多目标调度子问题
针对上述微电子生产线调度问题,采用如下SDEH算法求解当前阶段多目标调度子问题,形成当前阶段全局调度解;
步骤(4.1):个体编码和种群初始化
个体由在负载较大的前几个机器上加工的若干个操作串组成,每个操作串对应一个基因串,即个体由若干个基因串组成。
采用启发式规则:最短剩余加工时间优先(Shortest Remaining Processing Time,SRPT)、最先进入队列的优先(First Come First Served,FCFS),最短加工时间优先(Shortest ProcessingTime,SPT)、最长剩余加工时间优先(Maximum Remaining Processing Time,MRPT)及上述规则的随机加权规则,求解当前阶段调度子问题,根据下述方法获得负载较大的前几个机器,并形成由Z个个体组成的初始种群P:
设Ψ为预先给定的个体基因长度的最大值,为当前迭代阶段第n个工序第c个机器组缓存bnc中等待加工操作的总数,是当前调度子问题中在第n个工序第c个机器组上加工操作的总数;根据的大小对bnc进行排序,,设排列结果为bη1,bη2,…,bηΓ,其中,Γ为缓存的 总数;
然后,令
则个体由当前调度子问题中在与缓存bη1,bη2,…,bηγ对应的机器组上加工的操作串组成;每个基因串中的基因表示对应操作在相应机器组上的加工优先级,越靠前其加工优先级越高;
步骤(4.2)初始化
设总迭代步数为Λ,令θ=1,初始复制概率为Rep;
步骤(4.3):交叉
采用顺序交叉方法对随机配对的Z/2个体对进行交叉操作,形成子代种群P1,其中,对配对的两个个体的对应基因串单独进行顺序交叉;
步骤(4.4):变异
基于变异概率λ,对需变异的父代个体基因串S,从当前种群随机选择3个不同个体,与S对应的基因串分别设为S1、S2和S3,产生随机整数r∈{1,2,...,m},m为基因串S的基因长度;然后,根据下述流程生成子代个体对应的基因串S′:
步骤(4.4.1):生成取值在区间[0,1]内的m个随机数;
步骤(4.4.2):若上述第j(j=1,2,...,m)个随机数不小于基因复制概率Repij,且j≠r,则复制基因串S中的第j个基因作为子代个体相应基因串S′的第j个基因,设被复制基因的总数为Υ,将所有Υ个被复制基因从S删除,同时删除S1、S2和S3中对应位置的基因,那么基因串S、S1、S2和S3的长度均为m-Υ;Repij为对应基因的复制概率,根据下式自适应调整:
Repij=Rep ×βij
其中,βij为相应基因通过预测获得的资源冲突程度值,其定义为:设与基因xij对应的操作为Ouv,Δ为通过求解前一迭代阶段剩余调度子问题获得的调度问题解的总数,其中在δuv个解中,操作Ouv与其它操作间存在资源冲突,则
βij=δuv/Δ
步骤(4.4.3):设迭代次数i=1;
步骤(4.4.5):将填入到S′最左边的空白位置;
步骤(4.4.6):若i<m-Υ,则令i=i+1,执行步骤(4.4.4);否则,执行步骤(4.5);
步骤(4.5):轮盘赌选择,其中,对当前调度子问题中的每个个体,基于预测机制,针对当前阶段剩余调度子问题,采用所提出的MFRs对当前阶段剩余调度子问题求解,获得每个个体对应的全局调度目标函数值;
步骤(4.6):迭代终止条件判别
若θ≤Λ,则执行第5.3步,否则执行步骤(4.7);
步骤(4.7):构成当前阶段的全局调度解
由第四步所得到的当前阶段调度子问题的优化解F1和当前阶段剩余调度子问题相应的优化解F2共同构成当前阶段全局调度解F=F1∪F2;
第五步:固定当前阶段调度子问题中部分操作的加工机器和加工开始时间
步骤(5.1):若当前阶段全局调度解F优于迄今已获得的全局优化解Fo,则令Fo=F;
步骤(5.2):根据迄今已获得的全局优化解Fo和预先给定的待固定操作的比例κ,基于时间轴依次固定当前阶段调度子问题中部分操作的加工机器和加工开始时间;
步骤(6):采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习
在当前阶段调度子问题求解完成后,按下述方法利用相关调度数据采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习;
步骤(6.1):生成训练数据
按下述方法,根据当前阶段调度子问题具有最好全局调度性能的解F1生成用于MFRs自适应调整的M1个数据对,若当前阶段为第一阶段,则生成M个数据对;第m个数据对具有如下形式:
步骤(6.1.2):对上述M个决策时刻两个操作的每个个属性值按如下方法进行归一化处理,获得用于MFRs自适应调整的数据对:
m=1,2,…,M;k=1,2,…,K;j=1,2
步骤(6.2):若为第一阶段,则根据生成的M个数据对生成有效的MFRs;否则根据下述提出的调整方法TILM对用于求解剩余调度子问题的MFRs进行在线调整,其中,MFRs中每条规则具有如下形式:
…,
上述模糊规则具有多个前件和一个后件,是k个条件属性,k=1,2,…,K,为的第ikp个语言变量值,Yq是属性的语言变量值总数;当p=1,对应于需要排序的两个操作中的第一个操作;当p=2,对应于第二个;是由同一属性两个操作的不同取值对应的2维模糊网格。D是结论属性,表示两个操作那个具有更高的调度优先级,从而,分类标签只有两个取值,iR=1表示第一个操作具有更高的调度优
先级,iR=2表示第二个操作具有更高的调度优先级;CF(R)为规则R的确信度;
步骤(6.3):初始化
计算min FSa
min FSa=(minFS×M-M1)/M
步骤(6.4):数据移除
按下述方法更新保留的模糊网格;
步骤(6.4.1):从M个训练数据中移除最先生成的M1个训练数据:
步骤(6.4.2):对由M个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格, 计算其在剩余M-M1个训练数据下的模糊支撑度FSr;
步骤(6.4.3):若FSr≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格;
步骤(6.5):数据新增
按下述方法更新保留的模糊网格;
步骤(6.5.1):将求解当前阶段调度子问题过程中新生成的M1个训练数据增加到剩余的训练数据中:
步骤(6.5.2):对由剩余的M-M1个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格,计算其在加入新生成数据后的M个训练数据下的模糊支撑度FSa;
步骤(6.5.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格;
步骤(6.6):增加新的模糊网格
根据新生成的M1个训练数据,增加新的模糊网格:
步骤(6.6.1):根据新生成的M1个训练数据,生成模糊支撑度大于或等于min FS的模糊网格:
步骤(6.6.2):对于生成的模糊网格,基于M个训练数据,按步骤(6.5.2)计算其模糊支撑度FSa;
步骤(6.6.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格;
步骤(6.7):获得有效的MFRs
将所保留的模糊网格生成对应的多模糊规则MFRs;
步骤(6.7.1):根据上述保留的模糊网格获得有效的MFRs;其中,每条模糊分类规则R的模糊支撑度FS(R)作为其确信度,即:
CF(R)=FS(R)
步骤(6.7.2):考虑到两个操作在输入顺序改变的情形下,其调度优先级应该是不变的,从而,在多模糊规则中,每条规则R均生成一条逆规则R′与其对应,逆规则具有如下形式:
…,
其中,CF(R′)为规则R′的确信度,本发明中令CF(R′)=CF(R)
步骤(7):判定算法是否结束
若待调度问题的所有操作均已调度完成,则结束算法;否则,执行步骤(3)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210543787.7A CN103116324B (zh) | 2012-12-17 | 2012-12-17 | 基于指标预测和在线学习的微电子生产线调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210543787.7A CN103116324B (zh) | 2012-12-17 | 2012-12-17 | 基于指标预测和在线学习的微电子生产线调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103116324A true CN103116324A (zh) | 2013-05-22 |
CN103116324B CN103116324B (zh) | 2014-06-04 |
Family
ID=48414735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210543787.7A Expired - Fee Related CN103116324B (zh) | 2012-12-17 | 2012-12-17 | 基于指标预测和在线学习的微电子生产线调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103116324B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103324175A (zh) * | 2013-06-09 | 2013-09-25 | 广东豪美铝业股份有限公司 | 一种铝型材生产线的管控*** |
CN104536412A (zh) * | 2014-12-23 | 2015-04-22 | 清华大学 | 基于指标预报和解相似度分析的光刻工序动态调度方法 |
CN104698838A (zh) * | 2014-12-23 | 2015-06-10 | 清华大学 | 基于论域动态划分和学习的模糊调度规则挖掘方法 |
CN108255141A (zh) * | 2018-01-17 | 2018-07-06 | 北京理工大学 | 一种装配调度信息生成方法及*** |
CN110955197A (zh) * | 2018-09-27 | 2020-04-03 | 财团法人资讯工业策进会 | 产线自动配置装置及其方法 |
WO2021043007A1 (zh) * | 2019-09-06 | 2021-03-11 | 北京北方华创微电子装备有限公司 | 物料加工路径选择方法及装置 |
WO2021068274A1 (zh) * | 2019-10-10 | 2021-04-15 | 清华大学 | 基于操作完工时间快速预测的集成电路生产线调度方法 |
CN117215249A (zh) * | 2023-11-01 | 2023-12-12 | 北京航空航天大学 | 可重入制造***的动态调控方法及相关产品 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6456996B1 (en) * | 1998-06-05 | 2002-09-24 | I2 Technologies Us, Inc. | Computer implemented scheduling system and process using abstract local search technique |
US20030195794A1 (en) * | 2002-04-15 | 2003-10-16 | Yasuda To-Oru | Manufacture managing method |
CN101216710A (zh) * | 2007-12-28 | 2008-07-09 | 东南大学 | 一种由计算机实现的自适应选择动态生产调度控制*** |
CN101782769A (zh) * | 2010-03-08 | 2010-07-21 | 清华大学 | 一种基于指标补偿的平均流经时间快速预测方法 |
CN101788819A (zh) * | 2010-03-08 | 2010-07-28 | 清华大学 | 大规模生产过程一种基于迭代式分解和流松驰的调度方法 |
CN101794115A (zh) * | 2010-03-08 | 2010-08-04 | 清华大学 | 一种基于规则参数全局协调优化的调度规则智能挖掘方法 |
JP2011034432A (ja) * | 2009-08-04 | 2011-02-17 | Kobe Steel Ltd | スケジューリングシステム及びスケジューリング方法、並びにスケジューリングプログラム |
-
2012
- 2012-12-17 CN CN201210543787.7A patent/CN103116324B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6456996B1 (en) * | 1998-06-05 | 2002-09-24 | I2 Technologies Us, Inc. | Computer implemented scheduling system and process using abstract local search technique |
US20030195794A1 (en) * | 2002-04-15 | 2003-10-16 | Yasuda To-Oru | Manufacture managing method |
CN101216710A (zh) * | 2007-12-28 | 2008-07-09 | 东南大学 | 一种由计算机实现的自适应选择动态生产调度控制*** |
JP2011034432A (ja) * | 2009-08-04 | 2011-02-17 | Kobe Steel Ltd | スケジューリングシステム及びスケジューリング方法、並びにスケジューリングプログラム |
CN101782769A (zh) * | 2010-03-08 | 2010-07-21 | 清华大学 | 一种基于指标补偿的平均流经时间快速预测方法 |
CN101788819A (zh) * | 2010-03-08 | 2010-07-28 | 清华大学 | 大规模生产过程一种基于迭代式分解和流松驰的调度方法 |
CN101794115A (zh) * | 2010-03-08 | 2010-08-04 | 清华大学 | 一种基于规则参数全局协调优化的调度规则智能挖掘方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103324175A (zh) * | 2013-06-09 | 2013-09-25 | 广东豪美铝业股份有限公司 | 一种铝型材生产线的管控*** |
CN104536412A (zh) * | 2014-12-23 | 2015-04-22 | 清华大学 | 基于指标预报和解相似度分析的光刻工序动态调度方法 |
CN104698838A (zh) * | 2014-12-23 | 2015-06-10 | 清华大学 | 基于论域动态划分和学习的模糊调度规则挖掘方法 |
CN104698838B (zh) * | 2014-12-23 | 2017-03-29 | 清华大学 | 基于论域动态划分和学习的模糊调度规则挖掘方法 |
CN108255141A (zh) * | 2018-01-17 | 2018-07-06 | 北京理工大学 | 一种装配调度信息生成方法及*** |
CN108255141B (zh) * | 2018-01-17 | 2019-11-26 | 北京理工大学 | 一种装配调度信息生成方法及*** |
CN110955197A (zh) * | 2018-09-27 | 2020-04-03 | 财团法人资讯工业策进会 | 产线自动配置装置及其方法 |
WO2021043007A1 (zh) * | 2019-09-06 | 2021-03-11 | 北京北方华创微电子装备有限公司 | 物料加工路径选择方法及装置 |
US11703839B2 (en) | 2019-09-06 | 2023-07-18 | Beijing Naura Microelectronics Equipment Co., Ltd. | Material processing path selection method and device |
WO2021068274A1 (zh) * | 2019-10-10 | 2021-04-15 | 清华大学 | 基于操作完工时间快速预测的集成电路生产线调度方法 |
CN117215249A (zh) * | 2023-11-01 | 2023-12-12 | 北京航空航天大学 | 可重入制造***的动态调控方法及相关产品 |
CN117215249B (zh) * | 2023-11-01 | 2024-02-13 | 北京航空航天大学 | 可重入制造***的动态调控方法及相关产品 |
Also Published As
Publication number | Publication date |
---|---|
CN103116324B (zh) | 2014-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103116324A (zh) | 基于指标预测和在线学习的微电子生产线调度方法 | |
Young et al. | Optimizing deep learning hyper-parameters through an evolutionary algorithm | |
Ergezer et al. | Oppositional biogeography-based optimization | |
Sarkheyli et al. | Robust optimization of ANFIS based on a new modified GA | |
Hassan et al. | A hybrid of multiobjective Evolutionary Algorithm and HMM-Fuzzy model for time series prediction | |
US11783195B2 (en) | Process and system including an optimization engine with evolutionary surrogate-assisted prescriptions | |
CN104536412B (zh) | 基于指标预报和解相似度分析的光刻工序动态调度方法 | |
Joy et al. | Batch Bayesian optimization using multi-scale search | |
Shao et al. | A self-guided differential evolution with neighborhood search for permutation flow shop scheduling | |
CN106411896A (zh) | 基于apde‑rbf神经网络的网络安全态势预测方法 | |
Shrivastava et al. | GLAD: Learning sparse graph recovery | |
Bhatt et al. | Genetic algorithm applications on job shop scheduling problem: A review | |
CN103324954A (zh) | 一种基于树结构的图像分类方法及其*** | |
Donate et al. | Evolutionary optimization of sparsely connected and time-lagged neural networks for time series forecasting | |
Defourny et al. | Scenario trees and policy selection for multistage stochastic programming using machine learning | |
Buruzs et al. | Advanced learning of fuzzy cognitive maps of waste management by bacterial algorithm | |
CN109446420A (zh) | 一种跨域协同过滤方法和*** | |
Diwekar et al. | Efficient ant colony optimization (EACO) algorithm for deterministic optimization | |
Rad et al. | GP-RVM: Genetic programing-based symbolic regression using relevance vector machine | |
Jaddi et al. | Taguchi-based parameter designing of genetic algorithm for artificial neural network training | |
Pham et al. | Unsupervised training of Bayesian networks for data clustering | |
Bezáková et al. | Graph model selection using maximum likelihood | |
CN104698838A (zh) | 基于论域动态划分和学习的模糊调度规则挖掘方法 | |
Noekhah et al. | Software reliability prediction model based on ICA algorithm and MLP neural network | |
Liu et al. | Contrastive divergence learning for the restricted Boltzmann machine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140604 Termination date: 20181217 |
|
CF01 | Termination of patent right due to non-payment of annual fee |