CN106373024A - 一种面向具有学习恶化效应无等待流水调度问题的方法 - Google Patents
一种面向具有学习恶化效应无等待流水调度问题的方法 Download PDFInfo
- Publication number
- CN106373024A CN106373024A CN201610785498.6A CN201610785498A CN106373024A CN 106373024 A CN106373024 A CN 106373024A CN 201610785498 A CN201610785498 A CN 201610785498A CN 106373024 A CN106373024 A CN 106373024A
- Authority
- CN
- China
- Prior art keywords
- scheduling
- task
- wait
- degradation effects
- flowshop
- 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
- 238000000034 method Methods 0.000 title claims abstract description 70
- 230000002542 deteriorative effect Effects 0.000 title abstract 6
- 238000012545 processing Methods 0.000 claims abstract description 24
- 238000004519 manufacturing process Methods 0.000 claims abstract description 23
- 230000000694 effects Effects 0.000 claims description 69
- 230000015556 catabolic process Effects 0.000 claims description 60
- 238000006731 degradation reaction Methods 0.000 claims description 60
- 238000003780 insertion Methods 0.000 claims description 14
- 230000037431 insertion Effects 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 14
- 238000007781 pre-processing Methods 0.000 claims description 6
- 230000037213 diet Effects 0.000 claims description 5
- 235000005911 diet Nutrition 0.000 claims description 5
- 238000004378 air conditioning Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000006866 deterioration Effects 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 230000006870 function Effects 0.000 abstract description 21
- 229910000831 Steel Inorganic materials 0.000 description 4
- 239000010959 steel Substances 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000032683 aging Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- 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/30—Computing systems specially adapted for manufacturing
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Manufacturing & Machinery (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- General Factory Administration (AREA)
Abstract
本发明公开了一种面向具有学习恶化效应无等待流水调度问题的解决方法,对加工工序具有无等待流水调度的特点和基于位置的学习恶化效应函数的特点,建立无等待流水调度学习恶化效应模型;基于迭代局部搜索方法最小化该模型的总延迟目标函数,得出任务在调度中处于不同位置时的实际处理时间,进而得到当前的最优加工工序。本发明通过结合问题自身的特点,高质高效地求解具有学习恶化效应无等待流水调度问题,最小化总延迟目标函数。
Description
技术领域
本发明涉及一种面向具有学习恶化效应无等待流水调度问题的方法,属于复杂调度计算技术领域。
背景技术
无等待流水调度问题是调度问题中的一类重要分支。该问题源自特定产品要求任务的各工序紧密连接的生产制造特性,也就是要求任务从开始加工到结束必须连续执行,任务在各机器上的加工均不能出现等待。这一要求有时会造成某些任务在特定机器上的开工时间推迟到前一个工序完成时。该问题具有重要的现实意义,广泛应用于如化工生产、塑模成型、食品生产和钢材轧制等生产制造领域。例如,在钢材轧制生产中,钢锭不能在各生产工序间等待,否则会造成钢锭冷却无法加工。
学习和恶化效应大多是由人类或具有学习恶化效应的机器参与生产制造而产生的,在许多实际生产制造系统中普遍存在。例如,工人们会由于多次重复操作而提高技能,从而更快地完成如机器的操作、原材料的处理等任务;与之相反,工人们可能因为工作的中断而忘记部分技能,或机器与其零部件由于老化而增加工序的处理时间。
关于具有学习恶化效应的无等待流水调度问题,目前少有研究。然而,这一效应又在现实应用场景中广泛存在,具有深入研究的意义。本发明中将提出一种面向具有学习恶化效应无等待流水调度问题的方法,以期有效地解决该问题,最小化调度的总延迟。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种面向具有学习恶化效应无等待流水调度问题的方法,结合问题本身无等待和具有学习恶化效应的特征,采用迭代局部搜索方法框架,最小化调度的总延迟目标函数。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种面向具有学习恶化效应无等待流水调度问题的方法,对加工工序利用无等待流水调度问题的特点和基于位置的学习恶化效应函数的特点,建立无等待流水调度学习恶化效应模型。基于迭代局部搜索方法最小化求解该无等待流水调度学习恶化效应模型的总延迟目标函数,得出任务在调度中处于不同位置时的实际处理时间,进而得到当前的最优加工工序。
具体包括以下步骤:
步骤1.获取待加工工序和基于位置的学习恶化效应函数,根据无等待流水调度的特点,任务的完工时间表示为该任务之前各任务的完工时间距离之和,而各工序的实际加工时间将根据给定的学习恶化效应得到,则根据给定的学习恶化效应得到的实际加工时间得到实际任务的完工时间,根据实际任务的完工时间建立无等待流水调度学习恶化效应模型的总延迟目标函数。
步骤2.初始解生成:将初始解π0初始化为一个空的调度,在空调度π0中插入位于起始处的虚任务π[0]。将任务集中的各任务按照截止期非降序排序,每次选取具有最小截止期的任务,将之移除出排序,并利用前向交换算子AFS(π0,|π0|)将该任务从当前调度的末尾开始,依次插入调度π0的各位置中,选取其中具有最小目标函数的调度作为新的已有调度π0。继续从排序中移除任务,直到所有任务都已在调度π0中,则初始解生成。当前调度π和全局最优调度π*初始化为初始解π0。
步骤3.局部搜索:对于当前调度π,每次随机选取一个位置k,k=1,2,…,n,对当前调度π应用前向交换算子AFS(π,k)和后向交换算子ABS(π,k),分别得到经过前向交换算子应用后的调度πAFS和经过后向交换算子应用后的调度πABS,分别将经过前向交换算子应用后的调度πAFS和经过后向交换算子应用后的调度πABS与当前调度π比较,从中选取最优的调度作为新的当前调度π,直到所有位置都已被选取并操作完成。若新的当前调度π优于全局最优调度π*,则用新的当前调度π更新全局最优调度π*,并重复以上过程,否则进入步骤4。
步骤4.停止条件判断:判断方法执行过程是否达到了方法停止条件,如达到,则方法停止执行,返回全局最优调度π*,否则进入步骤5。
步骤5.扰动:对于当前调度π,利用移除插入算子RI(π,i,j)随机选取当前调度π中的任一一个任务,移除该任务并插入新的随机位置j,反复进行该操作,得到一个临时调度π′。
步骤6.局部搜索:对临时调度π′重复步骤3中过程,得到新的临时调度π″。
步骤7.接受准则:若新的临时调度π″优于当前调度π,则用新的临时调度π″更新当前调度π,否则,以一定的概率用新的临时调度π″更新当前调度π,若当前调度π优于全局最优调度π*,则用当前调度π更新全局最优调度π*。转到步骤4。
优选的:所述步骤1中得到无等待流水调度学习恶化效应模型的总延迟目标函数的方法:
获取待加工工序和给定的学习恶化效应,所述加工工序加工工序包括任务集机器集每个任务预加工时间,任务集包括n个任务J1,J2,…,Jn,机器集包括m台机器M1,M2,…,Mn,每个任务包括m个工序,需要在m台机器上按顺序加工,且具有每个任务预加工时间,而各工序的实际加工时间将根据给定的学习恶化效应得到。n个任务在m台机器上的加工顺序相同,并且是确定的。一台机器一次只能执行一个任务的一个工序,一个任务不能同时在不同机器上加工。一个任务的相邻两工序必须连续加工,之间不允许出现中断或等待。不允许抢占。
根据无等待流水调度的特点,任务的完工时间表示为该任务之前各任务的完工时间距离之和,即:
其中,C[i](π)为任务的完工时间,π表示一个包含有一个初始虚任务π[0]和n个任务的调度,Ω表示所有调度的集合,Di,j,r是调度中在第r个位置的任务Ji和第r+1个位置的任务Jj的实际完工时间距离,其定义为:
其中,pi,j,r表示任务Ji在机器Mj上处于调度中第r个位置时的工序实际加工时间,pi,j,r表示任务Ji在机器Mj上处于调度中第r个位置时的工序实际加工时间,pi,j任务Ji在机器Mj上处于调度中第r个位置时的工序加工时间,α,β,γ,μ为学习恶化效应函数参数。
进而得出无等待流水调度学习恶化效应模型的总延迟目标函数:
其中,TT(π)表示总延迟目标,d[i]表示第i个任务的截止期,L[i](π)和T[i](π)分别是调度π中第i个任务的完工时间差以及延迟,其定义为:
L[i](π)=C[i](π)-d[i]。 (6)
T[i](π)=max{L[i](π),0}。 (7)
优选的:所述步骤7中若新的临时调度π″优于当前调度π,则用新的临时调度π″更新当前调度π,否则,以的概率用新的临时调度π″更新当前调度π,其中T是温度值,e自然对数的底数。
优选的:所述前向交换算子,记为AFS(π,k)。该算子将调度π的第k个任务π[k]从原位置逐个与其前方的任务交换,直到遇到虚任务π[0]为止,得出并返回所有交换产生的新调度中最优的调度π′。
后向交换算子,记为ABS(π,k)。该算子将调度π的第k个任务π[k]从原位置逐个与其后方的任务交换,直到与调度中最后一个任务π[n]完成交换为止,得出并返回所有交换产生的新调度中最优的调度π′。
移除插入算子,记为RI(π,i,j)。该算子将调度π的第i个任务π[i]从原位置移除并插入新位置j中,返回新生成的调度π′。
有益效果:本发明提供的一种面向具有学习恶化效应无等待流水调度问题的方法,相比现有技术,具有以下有益效果:
本发明通过结合问题自身的特点,高质高效地求解具有学习恶化效应无等待流水调度问题,最小化总延迟目标函数。
附图说明
图1是本发明实施例求解具有学习恶化效应无等待流水调度问题方法的总体架构图;
图2是本发明实施例中解的初始化和初步搜索流程图;
图3是本发明实施例中解的迭代求解过程流程图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示为一种
一种面向具有学习恶化效应无等待流水调度问题的方法,如图1-3所示,对加工工序利用无等待流水调度问题的特点和基于位置的学习恶化效应函数的特点,建立无等待流水调度学习恶化效应模型。基于迭代局部搜索方法最小化求解该无等待流水调度学习恶化效应模型的总延迟目标函数,得出任务在调度中处于不同位置时的实际处理时间,进而得到当前的最优加工工序。本发明利用无等待流水调度问题的特点和基于位置的学习恶化效应函数的特点,建立问题的学习恶化效应模型;基于迭代局部搜索方法框架,提出求解具有学习恶化效应的无等待流水调度问题,最小化调度的总延迟。针对一个需要求解的问题,首先结合问题特点,建立问题的学习恶化效应模型,通过初始解生成模块生成初始调度,并对该初始调度进行初步的局部搜索之后,迭代通过扰动模块、局部搜索模块和接受准则模块对调度进行提高,直到停止条件判断达到停止标准,方法结束并返回找到最优解。
具体包括以下步骤:
步骤1.获取待加工工序和给定的学习恶化效应,根据无等待流水调度的特点,任务的完工时间表示为该任务之前各任务的完工时间距离之和,而各工序的实际加工时间将根据给定的学习恶化效应得到,则根据给定的学习恶化效应得到的实际加工时间得到实际任务的完工时间,根据实际任务的完工时间建立无等待流水调度学习恶化效应模型的总延迟目标函数。
获取待加工工序和给定的学习恶化效应,所述加工工序加工工序包括任务集机器集每个任务预加工时间,任务集包括n个任务J1,J2,…,Jn,机器集包括m台机器M1,M2,…,Mn,每个任务包括m个工序,需要在m台机器上按顺序加工,且具有每个任务预加工时间,而各工序的实际加工时间将根据给定的学习恶化效应得到。n个任务在m台机器上的加工顺序相同,并且是确定的。一台机器一次只能执行一个任务的一个工序,一个任务不能同时在不同机器上加工。一个任务的相邻两工序必须连续加工,之间不允许出现中断或等待。不允许抢占。
学习恶化效应的特点是:由人类或具有学习恶化效应的机器参与生产制造而产生的任务的实际处理时间缩短或延长的现象,在许多实际生产制造系统中普遍存在。例如,工人们会由于多次重复操作而提高技能,从而更快地完成任务;同时,机器与其零部件的老化也可能增加工序的处理时间。
无等待流水调度的特点就是:任务的完工时间表示为该任务之前各任务的完工时间距离之和,即:
其中,C[i](π)为任务的完工时间,π表示一个包含有一个初始虚任务π[0]和n个任务的调度,Ω表示所有调度的集合,Di,j,r是调度中在第r个位置的任务Ji和第r+1个位置的任务Jj的实际完工时间距离,其定义为:
其中,pi,j,r表示任务Ji在机器Mj上处于调度中第r个位置时的工序实际加工时间,
进而得出无等待流水调度学习恶化效应模型的总延迟目标函数:
其中,TT(π)表示总延迟目标,d[i]表示第i个任务的截止期,L[i](π)和T[i](π)分别是调度π中第i个任务的完工时间差以及延迟,其定义为:
L[i](π)=C[i](π)-d[i]。 (6)
T[i](π)=max{L[i](π),0}。 (7)
而建立无等待流水调度学习恶化效应模型,即,任务的实际处理时间仅与任务在调度中所处的实际位置有关。
步骤2.初始解生成:将初始解π0初始化为一个空的调度,在空调度π0中插入位于起始处的虚任务π[0]。将任务集中的各任务按照截止期非降序排序,每次选取具有最小截止期的任务,将之移除出排序,并利用前向交换算子AFS(π0,|π0|)将该任务从当前调度的末尾开始,依次插入调度π0的各位置中,选取其中具有最小目标函数的调度作为新的已有调度π0。继续从排序中移除任务,直到所有任务都已在调度π0中,则初始解生成。当前调度π和全局最优调度π*初始化为初始解π0。
步骤3.局部搜索:对于当前调度π,每次随机选取一个位置k,k=1,2,…,n,对当前调度π应用前向交换算子AFS(π,k)和后向交换算子ABS(π,k),分别得到经过前向交换算子应用后的调度πAFS和经过后向交换算子应用后的调度πABS,分别将经过前向交换算子应用后的调度πAFS和经过后向交换算子应用后的调度πABS与当前调度π比较,从中选取最优的调度作为新的当前调度π,直到所有位置都已被选取并操作完成。若新的当前调度π优于全局最优调度π*,则用新的当前调度π更新全局最优调度π*,并重复以上过程,否则进入步骤4。
步骤4.停止条件判断:判断方法执行过程是否达到了方法停止条件,如达到,则方法停止执行,返回全局最优调度π*,否则进入步骤5。
步骤5.扰动:对于当前调度π,利用移除插入算子RI(π,i,j)随机选取当前调度π中的任一一个任务,移除该任务并插入新的随机位置j,反复进行该操作,得到一个临时调度π′。
步骤6.局部搜索:对临时调度π′重复步骤3中过程,得到新的临时调度π″。
步骤7.接受准则:若新的临时调度π″优于当前调度π,则用新的临时调度π″更新当前调度π,否则,以一定的概率用新的临时调度π″更新当前调度π,若当前调度π优于全局最优调度π*,则用当前调度π更新全局最优调度π*。转到步骤4。
优选的:若新的临时调度π″优于当前调度π,则用新的临时调度π″更新当前调度π,否则,以的概率用新的临时调度π″更新当前调度π,其中T是温度值,e自然对数的底数。
优选的:所述前向交换算子,记为AFS(π,k)。该算子将调度π的第k个任务π[k]从原位置逐个与其前方的任务交换,直到遇到虚任务π[0]为止,得出并返回所有交换产生的新调度中最优的调度π′。
后向交换算子,记为ABS(π,k)。该算子将调度π的第k个任务π[k]从原位置逐个与其后方的任务交换,直到与调度中最后一个任务π[n]完成交换为止,得出并返回所有交换产生的新调度中最优的调度π′。
移除插入算子,记为RI(π,i,j)。该算子将调度π的第i个任务π[i]从原位置移除并插入新位置j中,返回新生成的调度π′。
下面以一具体实施进行说明。
表1本发明实施例
包括学习恶化效应模型的建立11、初始解生成12、局部搜索13、停止条件14、扰动15和接受准则16。假设任务集合机器集合各任务在各机器上原有的基本加工时间和截止期见表1所示。建立的学习恶化效应函数如公式(8)所示:
pi,j,r=pi,j-pi,j[1-(r+1)-α]μ+γ[1-(βr+1)e-βr]×pi,j[1-(r+1)-α]μ (8)
其中,pi,j,r表示任务Ji在机器Mj上处于调度中第r个位置时的工序实际加工时间,pi,j任务Ji在机器Mj上处于调度中第r个位置时的工序加工时间,α,β,γ,μ为学习恶化效应函数参数,其中,α=0.65,β=0.01,μ=8。
图2是本发明实施例中解的初始化和初步搜索流程图。如图2所示,具体步骤如下:
步骤s201,令J0为一个在各机器上加工时间均为0且截止期为0的需任务,初始化初始解π0=(J0);
步骤s202,将任务集合中的任务按照截止期非降序排序,得到排序{J1,J2,J3,J4};
步骤s203,从排序中移除有最小截止期的任务,当前情况下为J1,利用AFS(π0,|π0|)算子将之插入初始解π0中,得到最优部分调度π0=(J0,J1);
步骤s204,判断排序中是否还有任务,若是,转步骤s203,否则,转步骤s205;
经过多次迭代后,由步骤s204转s205时,已生成初始解π0=(J0,J4,J3,J2,J1),其对应的总延迟为33.20。
步骤s205,用初始解π0对当前调度π和全局最优调度π*进行赋值,即得到:π=(J0,J4,J3,J2,J1),π*=(J0,J4,J3,J2,J1);
步骤s206,生成(1,2,3,4)的随机序列v=(3,2,1,4);
步骤s207,初始化循环计数器i=1;
步骤s208,选取随机序列v[1]对应的位置3,对当前调度π应用算子AFS(π,3)和ABS(π,3),分别得到新的调度πAFS=(J0,J2,J4,J3,J1)和πABS=(J0,J4,J3,J1,J2);
步骤s209,πAFS与π相比的目标增量是12.15,而πABS与π相比的目标增量是4.83,因此新的当前调度π=π;
步骤s210,循环计数器i=i+1;
步骤s211,判断i是否大于n=4,若是,转步骤s213,否则转步骤s208;
步骤s212,更新全局最优调度π*←π,转步骤s206;
经过多次循环后,i>n=4,由步骤s211转s213时,已生成当前调度π=(J0,J4,J3,J2,J1),其对应的总延迟为33.20。
步骤s213,判断当前调度π是否优于全局最优调度π*,若是,转步骤s212,否则转步骤s301。
在此实施例中,由于TT(π)=TT(π*),故由步骤s213转步骤s301。
图3是本发明实施例中解的迭代求解过程流程图。如图3所示,该过程具体步骤如下:
步骤s301,判断方法停止条件是否满足,若是,则转步骤s318,否则转步骤s302;
步骤s302,对于当前调度π,多次应用算子XX随机移除再插入任务,从而得到调度π′;假设此处得到的调度π′=(J0,J1,J4,J3,J2);
步骤s303,生成(1,2,3,4)的随机序列v=(1,2,4,3);
步骤s304,初始化循环计数器i=1;
步骤s305,选取随机序列v[1]对应的位置1,对当前调度π应用算子AFS(π,1)和ABS(π,1),分别得到新的调度πAFS=(J0,J1,J4,J3,J2)和πABS=(J0,J4,J3,J2,J1);
步骤s306,πAFS与π′相比的目标增量是0,而πABS与π′相比的目标增量是-14.99,因此新的当前调度π″=πABS;
步骤s307,循环计数器i=i+1;
步骤s308,判断i是否大于n=4,若是,转步骤s310,否则转步骤s305;
步骤s309,更新全局最优调度π*←π″,转步骤s303;
经过多次循环后,i>n=4,由步骤s308转s310时,已生成当前调度π″=(J0,J4,J3,J2,J1),其对应的总延迟为33.20。
步骤s310,判断当前调度π″是否优于全局最优调度π*,若是,转步骤s309,否则转步骤s311;
在此实施例中,由于TT(π″)=TT(π*),故由步骤s310转步骤s311。
步骤s311,判断调度π″是否优于调度π,若是,转步骤s312,否则转步骤s313;
在此实施例中,转步骤s313;
步骤s312,更新当前调度π←π″,转步骤s316;
步骤s313,生成一随机数μ;
步骤s314,判断是否有若是,转步骤s315,否则转步骤s316;
在此实施例中,假设判断条件不满足,转步骤s316;
步骤s315,更新当前调度π←π″;
步骤s316,判断当前调度π是否优于全局最优调度π*,若是,转步骤s317,否则转步骤s301;
在此实施例中,由于TT(π)=TT(π*),故由步骤s316转步骤s301。
步骤s317,更新全局最优调度π*←π,转步骤s301;
步骤s318,返回全局最优解π*;方法结束。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种面向具有学习恶化效应无等待流水调度问题的方法,其特征在于:对加工工序利用无等待流水调度问题的特点和基于位置的学习恶化效应函数的特点,建立无等待流水调度学习恶化效应模型;基于迭代局部搜索方法最小化求解该无等待流水调度学习恶化效应模型的总延迟目标函数,得出任务在调度中处于不同位置时的实际处理时间,进而得到当前的最优加工工序。
2.根据权利要求1所述的面向具有学习恶化效应无等待流水调度问题的方法,其特征在于:包括以下步骤:
步骤1.获取待加工工序和基于位置的学习恶化效应函数,根据无等待流水调度的特点,任务的完工时间表示为该任务之前各任务的完工时间距离之和,各工序的实际加工时间将根据给定的学习恶化效应得到,则根据给定的学习恶化效应得到的实际加工时间得到实际任务的完工时间,根据实际任务的完工时间建立无等待流水调度学习恶化效应模型的总延迟目标函数;
步骤2.初始解生成:将初始解π0初始化为一个空的调度,在空调度π0中插入位于起始处的虚任务π[0];将任务集中的各任务按照截止期非降序排序,每次选取具有最小截止期的任务,将之移除出排序,并利用前向交换算子AFS(π0,|π0|)将该任务从当前调度的末尾开始,依次插入调度π0的各位置中,选取其中具有最小目标函数的调度作为新的已有调度π0;继续从排序中移除任务,直到所有任务都已在调度π0中,则初始解生成;当前调度π和全局最优调度π*初始化为初始解π0;
步骤3.局部搜索:对于当前调度π,每次随机选取一个位置k,k=1,2,…,n,对当前调度π应用前向交换算子AFS(π,k)和后向交换算子ABS(π,k),分别得到经过前向交换算子应用后的调度πAFS和经过后向交换算子应用后的调度πABS,分别将经过前向交换算子应用后的调度πAFS和经过后向交换算子应用后的调度πABS与当前调度π比较,从中选取最优的调度作为新的当前调度π,直到所有位置都已被选取并操作完成;若新的当前调度π优于全局最优调度π*,则用新的当前调度π更新全局最优调度π*,并重复以上过程,否则进入步骤4;
步骤4.停止条件判断:判断方法执行过程是否达到了方法停止条件,如达到,则方法停止执行,返回全局最优调度π*,否则进入步骤5;
步骤5.扰动:对于当前调度π,利用移除插入算子RI(π,i,j)随机选取当前调度π中的任一一个任务,移除该任务并插入新的随机位置j,反复进行该操作,得到一个临时调度π′;
步骤6.局部搜索:对临时调度π′重复步骤3中过程,得到新的临时调度π″;
步骤7.接受准则:若新的临时调度π″优于当前调度π,则用新的临时调度π″更新当前调度π,否则,以一定的概率用新的临时调度π″更新当前调度π,若当前调度π优于全局最优调度π*,则用当前调度π更新全局最优调度π*;转到步骤4。
3.根据权利要求2所述的面向具有学习恶化效应无等待流水调度问题的方法,其特征在于:所述步骤1中得到无等待流水调度学习恶化效应模型的总延迟目标函数的方法:
获取待加工工序和给定的学习恶化效应,所述加工工序加工工序包括任务集机器集每个任务预加工时间,任务集包括n个任务J1,J2,…,Jn,机器集包括m台机器M1,M2,…,Mn,每个任务包括m个工序,需要在m台机器上按顺序加工,且具有每个任务预加工时间,而各工序的实际加工时间将根据给定的学习恶化效应得到;n个任务在m台机器上的加工顺序相同,并且是确定的;一台机器一次只能执行一个任务的一个工序,一个任务不能同时在不同机器上加工;一个任务的相邻两工序必须连续加工,之间不允许出现中断或等待;不允许抢占;
根据无等待流水调度的特点,任务的完工时间表示为该任务之前各任务的完工时间距离之和,即:
其中,C[i](π)为任务的完工时间,π表示一个包含有一个初始虚任务π[0]和n个任务的调度,Ω表示所有调度的集合,Di,j,r是调度中在第r个位置的任务Ji和第r+1个位置的任务Jj的实际完工时间距离,其定义为:
其中,pi,j,r表示任务Ji在机器Mj上处于调度中第r个位置时的工序实际加工时间,进而得出无等待流水调度学习恶化效应模型的总延迟目标函数:
其中,TT(π)表示总延迟目标,d[i]表示第i个任务的截止期,L[i](π)和T[i](π)分别是调度π中第i个任务的完工时间差以及延迟,其定义为:
L[i](π)=C[i](π)-d[i]; (6)
T[i](π)=max{L[i](π),0} (7)。
4.根据权利要求3所述的面向具有学习恶化效应无等待流水调度问题的方法,其特征在于:所述步骤7中若新的临时调度π″优于当前调度π,则用新的临时调度π″更新当前调度π,否则,以的概率用新的临时调度π″更新当前调度π,其中T是温度值,e自然对数的底数。
5.根据权利要求2所述的面向具有学习恶化效应无等待流水调度问题的方法,其特征在于:所述前向交换算子,记为AFS(π,k);该算子将调度π的第k个任务π[k]从原位置逐个与其前方的任务交换,直到遇到虚任务π[0]为止,得出并返回所有交换产生的新调度中最优的调度π′;
后向交换算子,记为ABS(π,k);该算子将调度π的第k个任务π[k]从原位置逐个与其后方的任务交换,直到与调度中最后一个任务π[n]完成交换为止,得出并返回所有交换产生的新调度中最优的调度π′;
移除插入算子,记为RI(π,i,j);该算子将调度π的第i个任务π[i]从原位置移除并插入新位置j中,返回新生成的调度π′。
6.根据权利要求2所述的面向具有学习恶化效应无等待流水调度问题的方法,其特征在于:所述无等待流水调度学习恶化效应模型为:
pi,j,r=pi,j-pi,j[1-(r+1)-α]μ+γ[1-(βr+1)e-βr]×pi,j[1-(r+1)-α]μ
其中,pi,j,r表示任务Ji在机器Mj上处于调度中第r个位置时的工序实际加工时间,pi,j任务Ji在机器Mj上处于调度中第r个位置时的工序加工时间,α,β,γ,μ为学习恶化效应函数参数。
7.根据权利要求6所述的面向具有学习恶化效应无等待流水调度问题的方法,其特征在于:所述无等待流水调度学习恶化效应模型的初步搜索方法:
步骤s201,令J0为一个在各机器上加工时间均为0且截止期为0的需任务,初始化初始解π0=(J0);
步骤s202,将任务集合中的任务按照截止期非降序排序,得到排序{J1,J2,J3,J4};
步骤s203,从排序中移除有最小截止期的任务,当前情况下为J1,利用前向交换算子AFS(π0,|π0|)将之插入初始解π0中,得到最优部分调度π0=(J0,J1);
步骤s204,判断排序中是否还有任务,若是,转步骤s203,否则,转步骤s205;
经过多次迭代后,由步骤s204转s205时,已生成初始解π0=(J0,J4,J3,J2,J1);
步骤s205,用初始解π0对当前调度π和全局最优调度π*进行赋值,即得到:π=(J0,J4,J3,J2,J1),π*=(J0,J4,J3,J2,J1);
步骤s206,生成(1,2,3,4)的随机序列v=(3,2,1,4);
步骤s207,初始化循环计数器i=1;
步骤s208,选取随机序列v[1]对应的位置3,对当前调度π应用前向交换算子AFS(π,3)和ABS(π,3),分别得到前向交换算子应用后的调度πAFS=(J0,J2,J4,J3,J1)和经过后向交换算子应用后的调度πABS=(J0,J4,J3,J1,J2);
步骤s209,πAFS与π相比的目标增量是12.15,而πABS与π相比的目标增量是4.83,因此新的当前调度π=π;
步骤s210,循环计数器i=i+1;
步骤s211,判断i是否大于n=4,若是,转步骤s213,否则转步骤s208;
步骤s212,更新全局最优调度π*←π,转步骤s206;
经过多次循环后,i>n=4,由步骤s211转s213时,已生成当前调度π=(J0,J4,J3,J2,J1);
步骤s213,判断当前调度π是否优于全局最优调度π*,若是,转步骤s212,否则转步骤s301。
8.根据权利要求7所述的面向具有学习恶化效应无等待流水调度问题的方法,其特征在于:所述无等待流水调度学习恶化效应模型的迭代求解方法:
步骤s301,判断方法停止条件是否满足,若是,则转步骤s318,否则转步骤s302;
步骤s302,对于当前调度π,多次应用移除插入算子RI(π,i,j)随机移除再插入任务,从而得到临时调度π′;假设此处得到的临时调度π′=(J0,J1,J4,J3,J2);
步骤s303,生成(1,2,3,4)的随机序列v=(1,2,4,3);
步骤s304,初始化循环计数器i=1;
步骤s305,选取随机序列v[1]对应的位置1,对当前调度π应用前向交换算子AFS(π,1)和后向交换算子ABS(π,1),分别得到新的前向交换算子应用调度πAFS=(J0,J1,J4,J3,J2)和后向交换算子应用调度πABS=(J0,J4,J3,J2,J1);
步骤s306,πAFS与π′相比的目标增量是0,而πABS与π′相比的目标增量是-14.99,因此新的当前调度π″=πABS;
步骤s307,循环计数器i=i+1;
步骤s308,判断i是否大于n=4,若是,转步骤s310,否则转步骤s305;
步骤s309,更新全局最优调度π*←π″,转步骤s303;
经过多次循环后,i>n=4,由步骤s308转s310时,已生成新的当前调度π″=(J0,J4,J3,J2,J1);
步骤s310,判断新的当前调度π″是否优于全局最优调度π*,若是,转步骤s309,否则转步骤s311;
步骤s311,判断新的当前调度π″是否优于当前调度π,若是,转步骤s312,否则转步骤s313;
步骤s312,更新当前调度π←π″,转步骤s316;
步骤s313,生成一随机数μ;
步骤s314,判断是否有若是,转步骤s315,否则转步骤s316;
步骤s315,更新当前调度π←π″;
步骤s316,判断当前调度π是否优于全局最优调度π*,若是,转步骤s317,否则转步骤s301;
步骤s317,更新全局最优调度π*←π,转步骤s301;
步骤s318,返回全局最优解π*。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610785498.6A CN106373024B (zh) | 2016-08-29 | 2016-08-29 | 一种面向具有学习恶化效应无等待流水调度问题的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610785498.6A CN106373024B (zh) | 2016-08-29 | 2016-08-29 | 一种面向具有学习恶化效应无等待流水调度问题的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106373024A true CN106373024A (zh) | 2017-02-01 |
CN106373024B CN106373024B (zh) | 2020-02-04 |
Family
ID=57899660
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610785498.6A Active CN106373024B (zh) | 2016-08-29 | 2016-08-29 | 一种面向具有学习恶化效应无等待流水调度问题的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106373024B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112364526A (zh) * | 2020-11-30 | 2021-02-12 | 安徽大学 | 一种基于果蝇算法的模糊批调度方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101334663A (zh) * | 2008-06-11 | 2008-12-31 | 清华大学 | 基于加工能力特征模型的机器组调度方法 |
US7512454B1 (en) * | 2002-05-31 | 2009-03-31 | Advanced Micro Devices, Inc. | Display unit with processor and communication controller |
-
2016
- 2016-08-29 CN CN201610785498.6A patent/CN106373024B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7512454B1 (en) * | 2002-05-31 | 2009-03-31 | Advanced Micro Devices, Inc. | Display unit with processor and communication controller |
CN101334663A (zh) * | 2008-06-11 | 2008-12-31 | 清华大学 | 基于加工能力特征模型的机器组调度方法 |
Non-Patent Citations (3)
Title |
---|
TIAN CHEN等: "Integrated Iterated Local Search for the Permutation Flowshop Problem with Tardiness Minimization", 《IEEE会议论文》 * |
徐海燕: "一类基于开工时间恶化效应和基于位置学习效应的单机调度问题的求解", 《湖南师范大学自然科学学报》 * |
徐海燕等: "基于学习和恶化效应模型的单机调度", 《东南大学学报(自然科学版)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112364526A (zh) * | 2020-11-30 | 2021-02-12 | 安徽大学 | 一种基于果蝇算法的模糊批调度方法及系统 |
CN112364526B (zh) * | 2020-11-30 | 2024-07-02 | 安徽大学 | 一种基于果蝇算法的模糊批调度方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106373024B (zh) | 2020-02-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104636813B (zh) | 一种求解车间作业调度问题的混合遗传模拟退火算法 | |
CN107357267B (zh) | 基于离散花朵授粉算法求解混合流水线调度问题的方法 | |
CN105974794A (zh) | 基于改进的猫群算法的离散制造车间调度方法 | |
CN109871270B (zh) | 调度方案生成方法及装置 | |
CN110046226A (zh) | 一种基于分布词向量cnn-rnn网络的图像描述方法 | |
CN112907603B (zh) | 一种基于Unet和分水岭算法的细胞实例分割方法 | |
CN105975701A (zh) | 一种基于混合模糊模型的并行调度拆卸路径生成方法 | |
CN112381343A (zh) | 一种基于遗传-骨干粒子群混合算法的柔性作业车间调度方法 | |
CN105224577B (zh) | 一种多标签文本分类方法及系统 | |
CN109344497A (zh) | 一种模具基础本体的设计实现方法 | |
CN110221907A (zh) | 一种基于edf算法和模糊集的实时任务调度方法 | |
CN106373024A (zh) | 一种面向具有学习恶化效应无等待流水调度问题的方法 | |
CN115049246A (zh) | 一种面向订单扰动下的离散制造生产线重调度方法 | |
CN115933568A (zh) | 一种多目标分布式混合流水车间调度方法 | |
CN107871034A (zh) | 基于变尺度教与学算法的公差分配多目标优化设计方法 | |
CN104537425A (zh) | 一种汽车空调出风口的生产装配过程的优化调度方法 | |
CN113051363A (zh) | 题库智能知识点标注系统 | |
CN113139720A (zh) | 一种具有学习效应的铸造车间生产调度优化方法 | |
CN110321799B (zh) | 一种基于sbr和平均类间距离的场景数选择方法 | |
CN104698838A (zh) | 基于论域动态划分和学习的模糊调度规则挖掘方法 | |
CN116339942A (zh) | 一种基于强化学习的分布式训练任务的自适应调度方法 | |
CN113485278B (zh) | 优化两生产指标的柔性作业车间调度多目标分布估计方法 | |
CN113393490A (zh) | 目标检测方法及装置、计算机可存储介质 | |
CN108107848B (zh) | 一种基于最小空闲时间的流水线车间调度方法 | |
CN112784844B (zh) | 一种制作语义分割网训练样本的方法、系统及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |