CN107918806A - 一种智能的优化调度方法 - Google Patents
一种智能的优化调度方法 Download PDFInfo
- Publication number
- CN107918806A CN107918806A CN201711117188.8A CN201711117188A CN107918806A CN 107918806 A CN107918806 A CN 107918806A CN 201711117188 A CN201711117188 A CN 201711117188A CN 107918806 A CN107918806 A CN 107918806A
- Authority
- CN
- China
- Prior art keywords
- size
- sequence
- individual
- max
- stage
- 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
- 238000005457 optimization Methods 0.000 title abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000009826 distribution Methods 0.000 claims description 12
- 241000254158 Lampyridae Species 0.000 claims description 6
- 238000010521 absorption reaction Methods 0.000 claims description 6
- 238000012163 sequencing technique Methods 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 6
- 241000131894 Lampyris noctiluca Species 0.000 abstract 2
- 230000010006 flight Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 230000007547 defect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 239000012847 fine chemical Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 239000004065 semiconductor 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06316—Sequencing of tasks or work
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Primary Health Care (AREA)
- Manufacturing & Machinery (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种智能的优化调度方法。本发明利用改进的群智能优化方法对多处理机任务混合流水车间调度问题进行求解,生成最优调度方案。提出了一种基于Levy飞行的萤火虫算法,通过引入Levy飞行机制来避免算法陷入局部最优解,提升全局搜索的广度。此外,本发明在萤火虫算法中引入一系列工序调整规则,提升算法的收敛速度和精度。本发明能够有效缩短处理机的空闲等待时间,从而提高生产效率,提升生产车间的经济效益。
Description
技术领域
本发明涉及生产调度领域,具体地,涉及一种智能的优化调度方法。
背景技术
在工业界的实际生产过程中,许多工作的一个或多个阶段都由多台处理机共同完成。鉴于流水作业特征以及某些工序上存在并行机的特点,多处理机任务混合流水车间调度问题(hybrid flow-shop scheduling problem with multiprocessor tasks,HFSPMT)有着很强的工业背景,广泛存在于机械制造、冶金、精细化工、半导体制造,微处理器系统,机器视觉,物流等多个领域。
HFSPMT具有较大的复杂性,属于NP-hard问题。因此,传统的精确求解算法与启发式算法分别在求解时间与解的精度上无法满足生产要求,不适合用于HFSPMT的求解。而智能算法能够兼顾求解时间与解的精度,常用于各种类型的调度问题求解。目前,对HFSPMT的求解大多局限于遗传算法、粒子群算法等相对成熟的智能算法,而较少涉及新型的智能算法。另外,在求解复杂问题时,已有的智能算法通常容易陷入局部最优,并且存在收敛速度缓慢、精度不高等缺陷。
发明内容
为了克服目前求解HFSPMT智能算法的不足,丰富新型智能算法在HFSPMT中的应用,本发明目的在于提供一种能有效求解HFSPMT,生成最优调度方案的新型智能方法。
本发明解决其技术问题所采用的技术方案是:一种智能的优化调度方法,该方法应用基于Levy飞行的萤火虫算法以及新的解码方式。萤火虫算法作为一种新的群智能方法,在求解复杂优化问题时有着更强的寻优能力。引入Levy飞行可在局部最优值持续不变时,扩大搜索范围,跳出局部最优,避免发生过早收敛。新的解码算法则能使算法生成最优的调度方案。具体过程包括以下几个步骤:
1)已知一个包含n个工作的集合J={1,2,…,n},在有k个阶段的流水线上被处理,每个阶段i有mi个平行处理机,i=1,2,…,k,将每个工作视作k个任务的一个序列,每个阶段的任务必须在前一阶段的任务完成后才能够被处理。一件工作中的每个任务都需要对应阶段的一个或多个处理机同时连续地处理一段时间。用sizeij与pij表示工作j在阶段i所需的处理机数量与花费的时间;i=1,2,…,k,j∈J。将求解问题需要的size和p矩阵输入系统。
2)参数设置,种群个体数N、最大迭代次数tmax、随机参数α、个体吸引力β0、介质吸收率γ,Levy随机参数c;其中令N=20,tmax=500,α=0.5,β0=0.2,γ=1,c=1.5。
3)种群个体初始化。
生成种群X=(x1,x2,…,xN),种群中的第s个个体xs=(xs1,…,xsn),xsj为0~n之间的实数,s∈{1,2,…,N},j∈{1,2,…,n}。由于个体xs的坐标是连续的实数,而工作序列是离散的整数序列,用最小排序方法将连续坐标转化为工作序列,即将个体xs=(xs1,…,xsn)的各个维度从小到大排序,排序的序号构成的整数序列作为初始工作序列π1。
4)计算每个个体对应的最大完成时间Cmax。
萤火虫算法的目标函数为序列对应的最大完成时间Cmax。本发明基于先到先得的原则,根据前一阶段各工作的完成时间顺序构造下一阶段的工作序列,然后根据一定规则对生成的工作序列进行适当调整,灵活地进行工作排序,减少加工过程的空闲时间,最终得到最优的调度方案及最大完成时间Cmax。
4.1)i=1时,根据构造出第1个阶段的调度方案。其中h∈J,π1(h)表示序列π1中第h个元素的值,表示工作π1(h)在第1个阶段所需的处理机数量。
4.2)将第i阶段的各工作加工完成时间进行非递减顺序排序,得到第i+1个阶段的调度序列πi+1。对于任意的πi(h)和πi(l),h,l∈J,当且h<l时,进行πi+1排序时假定
4.3)令j=1,根据适当规则调整序列πi+1中第j和j+1个工作的次序。
对于阶段i中相邻的三个工作A,B,C,即A=πi,p,B=πi,p+1,C=πi,p+2。根据STiA与STiB的关系,可分别提出以下规则来得到更优的处理序列。
a.STiA>STiB
规则1-i,若sizeiB+sizeiC>mi,sizeiA+sizeiC>mi,则交换A与B的处理顺序。
规则1-ii,若sizeiB+sizeiC≤mi,sizeiA+sizeiC>mi,sizeiA+sizeiB>mi,且max{{STiB+piB,STiA}+piA,STiC}<max{STiA+piA+piB,STiC+piC},则交换A与B的处理顺序。
规则1-iii,若sizeiB+sizeiC>mi,sizeiA+sizeiC≤mi,则交换A与B的处理顺序。
规则1-iv,若sizeiB+sizeiC≤mi,sizeiA+sizeiC≤mi,则交换A与B的处理顺序。
b.STiA=STiB
规则2,若sizeiA+sizeiB>mi,sizeiB+sizeiC>mi,sizeiA+sizeiC≤mi,则交换A与B的处理顺序。
c.STiA<STiB
规则3,若STiA+piA>STiB,sizeiA+sizeiB>mi,sizeiB+sizeiC>mi,sizeiA+sizeiC≤mi,则交换A与B的处理顺序。
其中,STiA、STiB、STiC分别表示工作A、B、C在阶段i的最早开始时间,sizeiA、sizeiB、sizeiC分别表示工作A、B、C在阶段i所需的处理机数量,piA、piB、piC分别表示工作A、B、C在阶段i的处理时间,mi表示阶段i的处理机总量
4.4)若j=n-2计算πi+1中各工作的完成时间,继续;否则j=j+1,转至步骤4.3)。
4.5)若i=k-1,πi+1各工作的完成时间中最大值即为目标函数最大完成时间Cmax,继续;否则i=i+1,转步骤4.2)。
4.6)将该个体Cmax与全局最优值Gbest进行比较,若Cmax<Gbest,则令Gbest=Cmax,同时i=1,2,…,k,用记录该个体在阶段i的调度序列πi。
5)对于每个个体,如在种群中有Cmax更小的个体,则按照式(4)向该个体移动。
其中,v表示个体移动的速度,xL和yL是符合正态分布的随机数,它们的标准差为和c是Levy随机参数,和满足:
式(4)表示个体xs向更亮个体xb移动。其中Γ表示欧拉第二积分伽马函数,α为[0,1]间的值,β0为距离为0时的吸引力,常数γ为介质的吸收率,rsb为个体xs和个体xb的之间的欧几里德距离;c∈[0.3,1.99],决定了Levy飞行分布的性质,c越小则Levy飞行的特性越显著,c>2时基本等同于正态分布;v∈[-1,1]。随机搜索的方向由v决定,步长由v与随机参数α决定。Levy飞行行走步长满足一个重尾的稳定分布,在这种形式的行走中,短距离的探索与偶尔较长距离的行走相间。因此,引入Levy飞行能扩大算法搜索范围,增加种群多样性,并且更容易跳出局部最优点。
6)迭代次数达到tmax输出Gbest和调度方案i=1,2,…,k,否则转至步骤4)。
本发明的技术构思为:利用改进的群智能优化方法对多处理机任务混合流水车间调度问题进行求解,在提高算法随机搜索性能的同时,通过Levy飞行的自调节步长避免陷入局部最优。
本发明的有益效果主要表现在:提出的工序调整规则能加快算法的收敛速度,缩减算法的运行时间,通过引入Levy飞行机制来避免算法陷入局部最优解,提升全局搜索的广度,所述的智能的优化调度方法具有更强的随机搜索性能,同时能够避免陷入局部最优,生成最优调度方案。
附图说明
图1是本发明的流程图。
图2是根据技术方案4.3)中规则1-i交换工作A,B的示意图。
图3是根据技术方案4.3)中规则1-ii交换工作A,B的示意图。
图4是根据技术方案4.3)中规则2交换工作A,B的示意图。
图5是根据技术方案4.3)中规则2交换工作A,B的示意图。
图6是根据技术方案4.3)中规则3交换工作A,B的示意图。
图7是根据技术方案4.3)中规则3调整调度方案的示意图。
图8是示例的已知数据图。
图9是示例的调度结果甘特图。
具体实施方式
下面根据附图具体说明本发明。
假设2阶段混合流水车间要安排9个工作,每个阶段有5个处理机。操作时间和每个工作需要的处理机参照图8。
参照图1。本方法流程包括以下几个步骤:
1)矩阵size和p的值根据图8输入。
2)参数设置,种群个体数N、最大迭代次数tmax、随机参数α、个体吸引力β0、介质吸收率γ,Levy随机参数c;其中令N=20,tmax=500,α=0.5,β0=0.2,γ=1,c=1.5。
3)种群个体初始化。
生成种群X=(x1,x2,…,xN),种群中的第s个个体xs=(xs1,…,xs9),xsj为0~9之间的实数,s∈{1,2,…,N},j∈{1,2,…,9}。由于个体xs的坐标是连续的实数,而工作序列是离散的整数序列,用最小排序方法将连续坐标转化为工作序列,即将个体xs=(xs1,…,xs9)的各个维度从小到大排序,排序的序号构成的整数序列作为初始工作序列π1。以表1为例,来说明π1:
表1:最小排序方法示例
j | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
xsj | 1.89 | 0.65 | 1.24 | 2.43 | 5.89 | 5.12 | 3.95 | 7.53 | 8.27 |
π1中位置 | 3 | 1 | 2 | 4 | 7 | 6 | 5 | 8 | 9 |
根据表1所示的最小排序方法,得到第1阶段的工作序列π1=(231476589)。
4)计算每个个体对应的最大完成时间Cmax,整个过程参照图9。
萤火虫算法的目标函数为序列对应的最大完成时间Cmax。本发明基于先到先得的原则,根据前一阶段各工作的完成时间顺序构造下一阶段的工作序列,然后根据一定规则对生成的工作序列进行适当调整,灵活地进行工作排序,减少加工过程的空闲时间,最终得到最优的调度方案及最大完成时间Cmax。
4.1)i=1时,根据构造出第1个阶段的调度方案,h∈{1,2,…,9},π1(h)表示序列π1中第h个元素的值,表示工作π1(h)在第1个阶段所需的处理机数量。根据π1=(231476589)和相应的处理机要求从0时刻开始分配任务,整个过程保证调度方案不违背序列π1的先后关系,最后得到阶段1所有工作的调度方案,参照图9。
4.2)将第1阶段的各工作加工完成时间进行非递减顺序排序,得到第2个阶段的调度序列π2=(213764859)。
4.3)令j=1,根据规则开始调整序列π2中第j和j+1个工作的次序。
调整过程中发现,对于任意相邻的三个工作A,B,C,都有STiA<STiB。对照规则3,只有j=5时满足条件,交换A与B的处理顺序可得到π2=(213746859)。具体过程如图3-7所示。
4.4)若j=7计算π2中各工作的完成时间,继续;否则j=j+1,转至步骤4.3)。
4.5)i=k-1,因此π2=(213746859)各工作的完成时间中最大值即为目标函数最大完成时间Cmax,Cmax=20,调度方案参照图9。
4.6)将该个体Cmax=20与全局最优值Gbest进行比较,若Cmax<Gbest,则令Gbest=Cmax,同时记录该个体在两个阶段的调度序列π1和π2。
5)对于每个个体,如在种群中有Cmax更小的个体,则按照式(4)向该个体移动。
其中,v表示个体移动的速度,xL和yL是符合正态分布的随机数,它们的标准差为和c是Levy随机参数,和满足:
式(4)表示个体xs向更亮个体xb移动。其中Γ表示欧拉第二积分伽马函数,α为[0,1]间的值,β0为距离为0时的吸引力,常数γ为介质的吸收率,rsb为个体xs和个体xb的之间的欧几里德距离;c∈[0.3,1.99],决定了Levy飞行分布的性质,c越小则Levy飞行的特性越显著,c>2时基本等同于正态分布;v∈[-1,1]。随机搜索的方向由v决定,步长由v与随机参数α决定。Levy飞行行走步长满足一个重尾的稳定分布,在这种形式的行走中,短距离的探索与偶尔较长距离的行走相间。因此,引入Levy飞行能扩大算法搜索范围,增加种群多样性,并且更容易跳出局部最优点。
6)迭代次数达到tmax输出Gbest和调度方案否则转至步骤4)。
上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
Claims (1)
1.一种智能的优化调度方法。其特征是:应用基于Levy飞行的萤火虫算法以及新的解码方式,能有效求解多处理机任务混合流水车间调度问题,生成最优调度方案。具体过程包括以下几个步骤:
1)已知一个包含n个工作的集合J={1,2,…,n},在有k个阶段的流水线上被处理,每个阶段i有mi个平行处理机,i=1,2,…,k,将每个工作视作k个任务的一个序列,每个阶段的任务必须在前一阶段的任务完成后才能够被处理。一件工作中的每个任务都需要对应阶段的一个或多个处理机同时连续地处理一段时间。用sizeij与pij表示工作j在阶段i所需的处理机数量与花费的时间;i=1,2,…,k,j∈J。将求解问题需要的size和p矩阵输入系统。
2)参数设置,种群个体数N、最大迭代次数tmax、随机参数α、个体吸引力β0、介质吸收率γ,Levy随机参数c;其中令N=20,tmax=500,α=0.5,β0=0.2,γ=1,c=1.5。
3)种群个体初始化。
生成种群X=(x1,x2,…,xN),种群中的第s个个体xs=(xs1,…,xsn),xsj为0~n之间的实数,表示个体在空间中的位置,s∈{1,2,…,N},j∈{1,2,…,n}。由于个体xs的坐标是连续的实数,而工作序列是离散的整数序列,用最小排序方法将连续坐标转化为工作序列,即将个体xs=(xs1,…,xsn)的各个维度从小到大排序,排序的序号构成的整数序列作为初始工作序列π1。
4)计算每个个体对应的最大完成时间Cmax。
萤火虫算法的目标函数为序列对应的最大完成时间Cmax。本发明基于先到先得的原则,根据前一阶段各工作的完成时间顺序构造下一阶段的工作序列,然后根据一定规则对生成的工作序列进行适当调整,灵活地进行工作排序,减少加工过程的空闲时间,最终得到最优的调度方案及最大完成时间Cmax。
4.1)i=1时,根据构造出第1个阶段的调度方案。其中h∈J,π1(h)表示序列π1中第h个元素的值,该值代表的是一个工作,表示工作π1(h)在第1个阶段所需的处理机数量。
4.2)将第i阶段的各工作πi(h)的加工完成时间进行非递减顺序排序,得到第i+1个阶段的调度序列πi+1。对于任意的πi(h)和πi(l),h,l∈J,当且h<l时,进行πi+1排序时假定
4.3)令j=1,根据适当规则调整序列πi+1中第j和j+1个工作的次序。
对于阶段i中相邻的三个工作A,B,C,即A=πi,p,B=πi,p+1,C=πi,p+2,p=1,2,…,n-2。用STiA、STiB、STiC分别表示工作A、B、C在阶段i的最早开始时间,sizeiA、sizeiB、sizeiC分别表示工作A、B、C在阶段i所需的处理机数量,piA、piB、piC分别表示工作A、B、C在阶段i的处理时间,mi表示阶段i的处理机总量。
根据STiA与STiB的关系,可分别提出以下规则来得到更优的处理序列。
a.STiA>STiB
规则1-i,若sizeiB+sizeiC>mi,sizeiA+sizeiC>mi,则交换A与B的处理顺序。
规则1-ii,若sizeiB+sizeiC≤mi,sizeiA+sizeiC>mi,sizeiA+sizeiB>mi,且max{{STiB+piB,STiA}+piA,STiC}<max{STiA+piA+piB,STiC+piC},则交换A与B的处理顺序。
规则1-iii,若sizeiB+sizeiC>mi,sizeiA+sizeiC≤mi,则交换A与B的处理顺序。
规则1-iv,若sizeiB+sizeiC≤mi,sizeiA+sizeiC≤mi,则交换A与B的处理顺序。
b.STiA=STiB
规则2,若sizeiA+sizeiB>mi,sizeiB+sizeiC>mi,sizeiA+sizeiC≤mi,则交换A与B的处理顺序。
c.STiA<STiB
规则3,若STiA+piA>STiB,sizeiA+sizeiB>mi,sizeiB+sizeiC>mi,sizeiA+sizeiC≤mi,则交换A与B的处理顺序。
4.4)若j=n-2计算πi+1中各工作的完成时间,继续;否则j=j+1,转至步骤4.3)。
4.5)若i=k-1,πi+1各工作的完成时间中最大值即为目标函数最大完成时间Cmax,继续;否则i=i+1,转步骤4.2)。
4.6)将该个体Cmax与全局最优值Gbest进行比较,若Cmax<Gbest,则令Gbest=Cmax,同时用记录该个体在阶段i的调度序列πi。
5)对于每个个体,如在种群中有Cmax更小的个体,则按照式(4)向该个体移动。
<mrow>
<mi>v</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>x</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>|</mo>
<msub>
<mi>y</mi>
<mi>L</mi>
</msub>
<msup>
<mo>|</mo>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,v表示个体移动的速度,xL和yL是符合正态分布的随机数,它们的标准差为和c是Levy随机参数,和满足:
<mrow>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>x</mi>
<mi>L</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>c</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>c</mi>
</mrow>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msup>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>c</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>y</mi>
<mi>L</mi>
</msub>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>x</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msup>
<msub>
<mi>&gamma;r</mi>
<mrow>
<mi>s</mi>
<mi>b</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&alpha;</mi>
<mi>v</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式(4)表示个体xs向更亮个体xb移动。其中Γ表示欧拉第二积分伽马函数,α为[0,1]间的值,β0为距离为0时的吸引力,常数γ为介质的吸收率,rsb为个体xs和个体xb的之间的欧几里德距离;c∈[0.3,1.99],决定了Levy飞行分布的性质,c越小则Levy飞行的特性越显著,c>2时基本等同于正态分布;v∈[-1,1]。随机搜索的方向由v决定,步长由v与随机参数α决定。Levy飞行行走步长满足一个重尾的稳定分布,在这种形式的行走中,短距离的探索与偶尔较长距离的行走相间。因此,引入Levy飞行能扩大算法搜索范围,增加种群多样性,并且更容易跳出局部最优点。
6)迭代次数达到tmax输出Gbest和调度方案否则转至步骤4)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711117188.8A CN107918806B (zh) | 2017-11-13 | 2017-11-13 | 一种智能的优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711117188.8A CN107918806B (zh) | 2017-11-13 | 2017-11-13 | 一种智能的优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107918806A true CN107918806A (zh) | 2018-04-17 |
CN107918806B CN107918806B (zh) | 2021-01-26 |
Family
ID=61895472
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711117188.8A Expired - Fee Related CN107918806B (zh) | 2017-11-13 | 2017-11-13 | 一种智能的优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107918806B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109491344A (zh) * | 2018-12-12 | 2019-03-19 | 合肥工业大学 | 面向航天发动机研制过程的智能协同调度方法和系统 |
CN111506413A (zh) * | 2020-07-02 | 2020-08-07 | 上海有孚智数云创数字科技有限公司 | 一种基于业务效率优化的智能任务调度方法及系统 |
CN112766615A (zh) * | 2021-02-22 | 2021-05-07 | 首钢京唐钢铁联合有限责任公司 | 一种平行工序间量平衡的产线分工方法及装置 |
CN112836846A (zh) * | 2020-12-02 | 2021-05-25 | 红云红河烟草(集团)有限责任公司 | 一种卷烟发货多库点多方向联运调度双层优化算法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106611215A (zh) * | 2016-05-14 | 2017-05-03 | 四川用联信息技术有限公司 | 一种新的布谷鸟搜索算法解决作业车间调度问题 |
CN106651280A (zh) * | 2017-03-23 | 2017-05-10 | 广东工业大学 | 一种集装箱轮船物流运输调度方法及系统 |
-
2017
- 2017-11-13 CN CN201711117188.8A patent/CN107918806B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106611215A (zh) * | 2016-05-14 | 2017-05-03 | 四川用联信息技术有限公司 | 一种新的布谷鸟搜索算法解决作业车间调度问题 |
CN106651280A (zh) * | 2017-03-23 | 2017-05-10 | 广东工业大学 | 一种集装箱轮船物流运输调度方法及系统 |
Non-Patent Citations (3)
Title |
---|
DONNY SUTANTYO ETAL.: "Collective-adaptive Lévy flight for underwater multi-robot exploration", 《2013 IEEE INTERNATIONAL CONFERENCE ON MECHATRONICS AND AUTOMATION》 * |
包晓晓 等: "改进的萤火虫算法求解具有学习退化效应的JSP问题", 《数学理论与应用》 * |
彭郎军: "基于萤火虫算法的柔性作业车间调度问题研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109491344A (zh) * | 2018-12-12 | 2019-03-19 | 合肥工业大学 | 面向航天发动机研制过程的智能协同调度方法和系统 |
CN109491344B (zh) * | 2018-12-12 | 2020-07-07 | 合肥工业大学 | 面向航天发动机研制过程的智能协同调度方法和系统 |
CN111506413A (zh) * | 2020-07-02 | 2020-08-07 | 上海有孚智数云创数字科技有限公司 | 一种基于业务效率优化的智能任务调度方法及系统 |
CN112836846A (zh) * | 2020-12-02 | 2021-05-25 | 红云红河烟草(集团)有限责任公司 | 一种卷烟发货多库点多方向联运调度双层优化算法 |
CN112836846B (zh) * | 2020-12-02 | 2022-07-08 | 红云红河烟草(集团)有限责任公司 | 一种卷烟发货多库点多方向联运调度双层优化算法 |
CN112766615A (zh) * | 2021-02-22 | 2021-05-07 | 首钢京唐钢铁联合有限责任公司 | 一种平行工序间量平衡的产线分工方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107918806B (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107918806B (zh) | 一种智能的优化调度方法 | |
CN107703900A (zh) | 一种高效的优化调度方法 | |
Aghajani et al. | Balancing of robotic mixed-model two-sided assembly line with robot setup times | |
CN112783172B (zh) | 基于离散鲸鱼优化算法的agv与机器集成调度方法 | |
Fan et al. | Solving the two-stage hybrid flow shop scheduling problem based on mutant firefly algorithm | |
CN113516429B (zh) | 一种基于网络拥堵模型的多agv全局规划方法 | |
CN106611275A (zh) | 针对作业车间生产问题的排产算法 | |
CN114186791A (zh) | 一种面向多型号小批量的复杂装备产品装调生产动态调度方法 | |
CN107831740A (zh) | 一种应用于笔记本零件的分布式生产制造过程中的优化调度方法 | |
CN109445386A (zh) | 一种基于onba的云制造任务最短生产时间调度方法 | |
Kim et al. | Scheduling of dual-gripper robotic cells with reinforcement learning | |
CN115034444B (zh) | 基于学习效应的多目标双柔性作业车间调度方法及系统 | |
Chamnanlor et al. | Hybrid genetic algorithms for solving reentrant flow-shop scheduling with time windows | |
CN107437121B (zh) | 适用于多机器同时处理单工件的生产工艺过程控制方法 | |
CN113326970B (zh) | 一种混流装配线排序优化方法 | |
CN108181810B (zh) | 一种基于区间不确定性的鲁棒单机调度方法 | |
CN107831746A (zh) | 一种高效的航空发动机装配车间调度系统 | |
CN107808222A (zh) | 一种智能的集装箱码头装卸调度装置 | |
CN111596622A (zh) | Ecm规则分布估计算法的柔性作业车间调度方法 | |
CN110689216B (zh) | 一种基于混合果蝇算法的机器人装配单元混合产品jit节能生产优化方法 | |
CN108873850A (zh) | 一种机械制造自动化生产过程的优化调度方法 | |
CN108107848B (zh) | 一种基于最小空闲时间的流水线车间调度方法 | |
CN107918805A (zh) | 一种智能的航空发动机装配车间调度系统 | |
CN112631214B (zh) | 基于改进入侵杂草优化算法的柔性作业车间分批调度方法 | |
CN108215202A (zh) | 一种考虑打印质量的3d打印分批控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210126 Termination date: 20211113 |