CN109886588A - 一种基于改进鲸鱼算法求解柔性作业车间调度的方法 - Google Patents

一种基于改进鲸鱼算法求解柔性作业车间调度的方法 Download PDF

Info

Publication number
CN109886588A
CN109886588A CN201910149784.7A CN201910149784A CN109886588A CN 109886588 A CN109886588 A CN 109886588A CN 201910149784 A CN201910149784 A CN 201910149784A CN 109886588 A CN109886588 A CN 109886588A
Authority
CN
China
Prior art keywords
whale
scheduling
algorithm
solution
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
Application number
CN201910149784.7A
Other languages
English (en)
Other versions
CN109886588B (zh
Inventor
栾飞
李富康
蔡宗琰
吴书强
杨�嘉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Rdy Technology Co ltd
Original Assignee
Changan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changan University filed Critical Changan University
Priority to CN201910149784.7A priority Critical patent/CN109886588B/zh
Publication of CN109886588A publication Critical patent/CN109886588A/zh
Application granted granted Critical
Publication of CN109886588B publication Critical patent/CN109886588B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于改进鲸鱼算法求解柔性作业车间调度的方法,包括以下步骤:1)建立柔性作业车间调度问题的数学模型;2)设置算法参数并生成初始种群;3)得到当前最优的调度解;4)判断当前迭代次数是否大于最大迭代次数;是,输出调度解;否,当前最优个体的计数器值是否不小于预设值;是,进行变邻域搜索操作,更新调度解;否,将调度解转换为鲸鱼个体位置向量,并保留调度解对应的鲸鱼个体;5)采用改进鲸鱼算法更新鲸鱼个体位置信息,将鲸鱼个体位置向量转换为调度解完成种群更新,迭代次数加1,返回3)。本方法可以很好的求解得柔性作业车间调度的所有最优解,并提高求解速度和精度。

Description

一种基于改进鲸鱼算法求解柔性作业车间调度的方法
技术领域
本发明属于柔性作业车间调度领域,涉及一种基于改进鲸鱼算法求解柔性作业车间调度的方法。
背景技术
柔性作业车间调度问题(Flexible Job-shop Scheduling Problem,FJSP)是经典调度问题的扩展,由于其考虑了工件柔性加工路径的特性,更符合车间生产的实际情况,但机器约束的减少,也增加了问题的求解难度,因此,该问题的求解算法也成为近些年车间调度领域研究的热点之一。目前,智能优化算法是解决这一问题的主要方法,它为FJSP的解决提供了更多思路和途径,而引起了国内外学者的广泛关注。
虽然近年来学者们提出了各种元启发式算法用于求解FJSP,包括基于云计算的改进细菌觅食算法、混合遗传算法以及粒子群算法等等,但目前仍没有那一种算法可以求解得FJSP的所有最优解,因此学者们仍在不断积极探索,以期找到一种能够更好得求解得最优解的求解方法。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供一种基于改进鲸鱼算法求解柔性作业车间调度的方法。
为达到上述目的,本发明采用以下技术方案予以实现:
一种基于改进鲸鱼算法求解柔性作业车间调度的方法,包括以下步骤:
步骤1:对柔性作业车间调度的机器选择和工序排序进行编码,建立柔性作业车间调度问题的数学模型;
步骤2:设置改进鲸鱼算法种群规模、最大迭代次数和位置元素,并利用混沌反向学习策略和搜索法生成初始种群;
步骤3:定义适应度函数:
其中,Fmax表示工件的最大完工时间,表示工件i的完工时间,n表示工件数;计算种群中各调度解的适应度值,得到当前最优的调度解;
步骤4:判断当前迭代次数是否大于最大迭代次数;否,则执行步骤5-8;是,则执行步骤9;
步骤5:为最优鲸鱼个体设置计数器,若改进鲸鱼算法迭代一次后最优鲸鱼个体未发生改变,则计数器值加1,判断当前最优鲸鱼个体的计数器值是否不小于预设值;是,则执行步骤6;否,则执行步骤7;
步骤6:对调度解执行变邻域搜索操作,更新调度解;
步骤7:将调度解转换为鲸鱼个体位置向量,并保留调度解对应的鲸鱼个体;
步骤8:采用改进鲸鱼算法更新鲸鱼个体位置信息;将鲸鱼个体位置向量转换为调度解S完成种群更新,令t=t+1;其中,t为当前迭代次数,转至步骤3;
步骤9:输出调度解。
本发明进一步的改进在于:
改进鲸鱼算法包括环绕捕食、发泡网攻击和随机搜索;其中:
环绕捕食的数学模型为:
其中:t表示当前迭代次数,是系数向量,是当前最优鲸群个体的位置向量,是当前鲸群个体的位置向量,| |是取绝对值,·是逐元素相乘,在迭代过程中其内部元素非线性收敛因子a在发泡网攻击和随机搜索中从2线性递减到0,表达式为tmax为最大迭代次数;为随机向量,r取值范围为[0,1];
发泡网攻击包括收缩包围和螺旋式位置更新;收缩包围时的变动范围随着的减少而缩小,的取值范围为[-1,1];螺旋式位置更新的数学模型为:
其中,为鲸鱼最佳位置与猎物之间的距离;b为定义的控制对数螺线形状的常数,l为[-1,1]之间的随机数;同时引入自适应权重ω:
则发泡网攻击的数学模型为:
随机搜索的数学模型为:
其中, 为从当前种群中随机选择的个体位置矢量,p是[0,1]中的一个随机数。
步骤1中采用等长的两段式编码方式对柔性作业车间调度的机器选择和工序排序进行编码。
步骤6的具体方法为:
步骤6-1:构建邻域结构N1、N2和N3
步骤6-2:进行変邻域搜索,具体包括以下步骤:
步骤6-2-1:将当前最优个体位置X*作为变邻域初始解X,设置λ←1,qmax←3最大迭代次数ηmax
步骤6-2-2:设置q←1;
步骤6-2-3:若q=1,则X'←N1(X);若q=2,则X'←N2(X);若q=3,则X'←N3(X);
步骤6-2-4:将X'作为初始解进行局部搜索,得到最优解X”;
步骤6-2-5:若X”优于X,则X←X”,并设置q←1;否则,设置q←q+1;
步骤6-2-6:判断q>qmax是否满足,若满足,则设置q←q+1,转到步骤6-2-7;否则,转到步骤6-2-3;
步骤6-2-7:判断η>ηmax是否满足,若满足,结束;否则,转到步骤6-2-2。
步骤6-1包括以下步骤:
步骤6-1-1:构建邻域结构N1:在工序排序部分的码段中任选两个元素,将选中的两个元素进行位置交换操作;
步骤6-1-2:构建邻域结构N2:在工序排序部分的码段中任选两个元素,将选中的两个元素中位置靠后的元素插入到位置靠前的元素前面的位置;
步骤6-1-3:构建邻域结构N3:在机器分配部分的码段中任选一个元素,该元素的可加工设备数量大于1台,然后将该元素的对应的加工设备变为其他可选设备,进而通过方案转换个体位置的方式生成新的排序方案。
步骤6-2-4中局部搜索采用阈值接受法,包括以下步骤:
步骤6-2-4-1:获取初始解X',并设置阈值δ>0,γ←1,ρ←1和终止条件γmax
步骤6-2-4-2:若ρ=1,则X”←N1(X')∪N3(X');若ρ=0,则X”←N2(X')∪N3(X');
步骤6-2-4-3:判断Cmax(X”)-Cmax(X')≤δ是否满足;如果满足,则X'←X”;否则,设置ρ←|ρ-1|;
步骤6-2-4-4:令γ←γ+1,判断γ>γmax是否满足;如果满足,则X”←X',转到步骤6-2-4-5;否则,转到步骤6-2-4-2;
步骤6-2-4-5:局部搜索结束。
将调度解S转换为鲸鱼个体位置向量的具体方法为:
S1:机器选择:采用下式将工序可选机器集中已选机器的序号转换为鲸鱼个体位置向量元素值,
x(i)=[2m/(s(i)-1)](n(i)-1)-m,s(i)≠1
其中:x(i)表示个体位置向量的第i个元素;s(i)表示元素i对应工序可以选择的机器个数;m表示机器总数;n(i)∈[1,s(i)]表示选定机器在可选机器集里面的序号;
如果s(i)=1,则x(i)在[-m,m]内任意取值;
S2:工序排序:首先生成[-m,m]内的一组与工序排序相对应的随机数,并按升序排列规则为每个随机数赋予一个唯一的ROV值,使得每个ROV值对应一个工序,然后根据工序的编码顺序对ROV值进行重排,重排后的ROV值对应的随机数顺序即为鲸鱼个体位置向量中各元素的值,组合为鲸鱼个体位置向量。
将鲸鱼个体位置向量转换为调度解S的具体方法为:
R1:机器选择:按照下式
得到选定机器在机器集里面的序号;其中:xij表示i工件的第j道工序所对应的个体位置元素;sij表示i工件的第j道工序可以选择的机器个数;nij∈[1,sij]表示选定机器在机器集里面的序号;如果sij=1,则xij在[-m,m]内任取数值,2m表示个体位置向量的长度;
R2:工序排序:首先给每一个鲸鱼个体位置元素按升序的顺序赋予对应的ROV值,然后以ROV值作为元素编号,将ROV值与工序相对应构造得出相应的调度解S。
与现有技术相比,本发明具有以下有益效果:
通过采用基于混沌反向学习策略和搜索法的种群初始化方法来提高初始种群的质量;其次,针对当前最优个体执行变邻域搜索操作,增强了算法的局部搜索能力;采用两段式编码方式,有效的实现了柔性作业车间调度最优解的离散值和改进鲸鱼算法的连续值之间的转换;引入非线性收敛因子和自适应权重平衡了算法的全局搜索和局部开发的能力,使算法可以有效的求解得到柔性作业车间调度的所有最优解。通过对15个标准算例的仿真实验,验证了算法各种改善机制的有效性,同时也测验了算法在求解速度和精度方面的优越性,为后续的深入工作提供了一种新途径。
附图说明
图1为本发明的求解柔性作业车间调度的方法流程图;
图2为本发明的FJSP编码方案示意框图;
图3为本发明的调度方案转换为个体位置向量转换过程框图;
图4为本发明的个体位置向量转换为调度方案转换框图;
图5为本发明实施例3的迭代收敛曲线图;
图6为本发明实施例3的调度结果甘特图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面结合附图对本发明做进一步详细描述:
参见图1,本发明基于改进鲸鱼算法求解柔性作业车间调度的方法,包括以下步骤:
步骤1:对柔性作业车间调度的机器选择和工序排序进行编码,建立柔性作业车间调度问题的数学模型;
步骤2:设置改进鲸鱼算法种群规模、迭代次数和位置元素,并利用混沌反向学习策略和搜索法生成初始种群;
步骤3:定义适应度函数:
其中,Fmax表示工件的最大完工时间,表示工件i的完工时间,n表示工件数;计算初始种群中各调度解的适应度值,得到当前最优的调度解;
步骤4:判断是否满足算法的终止条件;若不满足,则执行步骤5-8;如若满足,则执行步骤9;
步骤5:为最优鲸鱼个体设置计数器,若改进鲸鱼算法迭代一次后未发生改变,则计数器值加1,判断当前最优个体的计数器值是否满足≥预设值,如果满足,则执行步骤6;否则,执行步骤7;
步骤6:对调度解执行变邻域搜索操作,更新调度解;
步骤7:将调度解S转换为鲸鱼个体位置向量,并保留调度解所对应的鲸鱼个体;
步骤8:采用改进鲸鱼算法更新鲸鱼个体位置信息;将鲸鱼个体位置向量转换为调度解完成种群更新,令t=t+1;其中,t为当前迭代次数,转至步骤3;
步骤9:输出最优结果和对应调度解。
下面详细介绍本发明方法的内容:
一、柔性作业车间调度的问题模型
FJSP指的是n个工件如何在m台机器上安排加工的问题,其中,每个工件i的加工由Ji道工序所组成,每道工序可以在一台以上的设备上完成加工,每台机器能够加工的工件种类和加工所需要的时间已知,在此基础上,寻找最合理的加工分配方案,力求最小化最大完工时间,该问题满足的假设条件如下:
所有工件在初始时刻均处于可加工状态;(2)一道工序只能完整的在同一台机器上加工完成,中间不允许移动;(3)一台机器同时只能加工一道工序;(4)同一个工件的工序之间加工顺序是固定的,不同工件的工序之间加工顺序相互独立。
令Fmax表示工件的最大完工时间,表示工件i的完工时间,据此可以将目标函数表示为:
二、基本鲸鱼算法
鲸鱼优化算法(WOA)是模仿座头鲸觅食行为的新型智能优化算法。座头鲸发现猎物后,以螺旋方式收缩包围猎物,同时发出气泡网进行觅食。具体觅食方式包括“环绕捕食”“发泡网攻击”“随机捕食”等三种捕食方式。其中“发泡网攻击”又包括“收缩包围”和“螺旋式位置更新”两个步骤,因此座头鲸的觅食方式有以下数学描述:
2.1环绕捕食
由于猎物的位置(最佳位置)在搜索空间中是未知的,所以WOA假设当前鲸群最佳个体位置为最接近目标猎物的个体位置,鲸群其他个体将通过朝最佳个体位置方向移动进行更新。该行为表示为:
其中,t表示当前迭代次数,是系数向量,是当前最优个体的位置向量,是当前鲸群个体的位置向量,| |是取绝对值,·是逐元素相乘,在迭代过程中其内部元素(发泡网攻击和搜索捕食阶段)从2线性递减到0,表达式为其中M为最大迭代次数;为随机向量,r取值范围为[0,1]。
2.2发泡网攻击
觅食过程中,座头鲸沿着螺旋路径逐渐收缩到猎物周围,为了对这种行为进行建模,假设在优化过程中,更新鲸鱼的位置时,以50%的概率作为阈值来选择收缩包围机制或螺旋式位置更新:
2.2.1收缩包围机制
这种行为是通过减少公式(4)中的a值来实现的;其中的变动范围也随着的减少而缩小,此时的取值在[-1,1]之间。
2.2.2螺旋式位置更新
首先计算鲸群个体和位于最优位置处的猎物之间的距离,然后在鲸群个体和猎物的位置之间创建一个螺旋数学模型,以模仿座头鲸的螺旋形运动。数学模型如下所示:
则座头鲸的发泡网觅食行为的数学模型如下:
p是[0,1]中的一个随机数。
2.3随机搜索
随机搜索是指鲸群个体根据随机选择的鲸鱼个体更新位置,而不是目前找到的最佳鲸鱼个体。在值大于1或者小于-1时,迫使鲸鱼个体偏离当前猎物位置,以此搜索其他更合适的猎物。数学模型如下:
其中,是从当前种群中随机选择的个体位置矢量。
三、改进鲸鱼算法求解柔性作业车间调度问题
3.1FJSP的编码机制
FJSP包含机器选择和工序排序两个子问题,针对这一特点,采用等长的两段式编码方式来对问题加以描述。如对于3个工件每个工件都有2道工序的FJSP,其个体位置向量总长度为12,个体位置元素取值范围为[-3,3],其中前6位表示6道工序的机器选择,后6位表示加工顺序的安排,编码方案如图2所示,OPij表示工件i的第j道工序。
3.2转换机制
FJSP的解是离散值,WOA的个体位置向量是连续值。为了解决这一差异,有效地利用WOA来求解FJSP,采用下述方法实现二者的转换。
1)调度方案转换为个体位置向量
i)机器选择:按照式(11)将工序可选机器集中已选机器的序号转换为对应个体位置向量元素值,如下式所示:
x(i)=[2m/(s(i)-1)](n(i)-1)-m,s(i)≠1 (11)
其中:x(i)表示个体位置向量的第i个元素;s(i)表示元素i对应工序可以选择的机器个数;n(i)∈[1,s(i)]表示选定机器在可选机器集里面的序号。如果s(i)=1,则x(i)在[-m,m]内任意取值。
ii)工序排序:首先生成[-m,m]内的一组与工序排序方案相对应的随机数,并按升序排列ROV(Rankedordervalue)规则为每个随机数赋予一个唯一的ROV值,使得每个ROV值对应一个工序,然后根据工序的编码顺序对ROV值进行重排,重排后的ROV值所对应的随机数顺序即为个体位置向量中各元素的值。转换过程如图3所示。
2)个体位置向量转换为调度方案
i)机器选择:按照式(11)的反向推导可得机器的编号如下式(12)所示。
其中:xij表示i工件的第j道工序所对应的个体位置元素;sij表示i工件的第j道工序可以选择的机器个数;nij∈[1,sij]表示选定机器在机器集里面的序号。如果sij=1,则xij在[-m,m]内任取数值,2m表示个体位置向量的长度。
ii)工序排序:首先给每一个个体位置元素按升序的顺序赋予对应的ROV值,然后以ROV值作为元素编号,将其与工序相对应即可构造得出相应的新的排序方案。转换过程如图4所示。
3.3种群初始化
种群初始化将针对机器分配和工序排序两个子问题分层展开,然后利用公式(11)将其转换为个体位置向量。首先対机器分配部分采用混沌反向学习策略生成一组较优的分配方案;然后,采用搜索的方式随机生成若干组工序排序方案,并将其依次与已生成的机器分配方案组合之后求对应的适应度值,选择适应度值最好的方案对为初始解,以此来完成初始种群的生成。
3.4非线性收敛因子与自适应权重更新策略
在WOA中,其全局搜索和局部开发能力主要依赖于收敛因子a,其对于算法的求解能力有很大的的影响,其值大,则算法搜索能力较强,可以在较大区域内展开搜索,其值小,则只能在小范围内进行搜索,因此,为了提高WOA的搜索效率,在算法前期搜索阶段可以将a设置大一些,有利于在较大范围展开寻优,后期设置小一些则便于在最优解周围进行针对性的精细搜索。对收敛因子a进行改进,其公式为:
通过对a的改进可以一定程度上提高算法的寻优能力,但并不能达到有效平衡算法局部搜索和全局探索能力的目的,为此,引入自适应权重ω与非线性收敛因子a相配合来协调算法全局搜索和局部开发的能力,通过自适应权重ω来提高算法的的寻优性能,其公式如(14)所示,改进后的更新公式如式(15)所示。
则随机搜索的数学模型为:
其中, 为从当前种群中随机选择的个体位置矢量,p是[0,1]中的一个随机数。
3.5邻域搜寻模式
在基本WOA搜索过程中,种群个体根据当前最优个体X*来更新位置信息,因此,X*的多样性也决定了WOA算法跳出局部最优解的能力,影响算法的计算性能,本文将变领域搜索嵌入到WOA中,并为X*设置“计数器”,若WOA迭代一次后X*未发生改变,则“计数器”加1,否则不变。当计数器超过预设的上限(本文取15),则对执行变邻域搜索,使其跳出局部最优值。变邻域搜索的算法步骤如下:
构建邻域结构:1)邻域结构N1:在工序排序部分的码段中任选两个元素,将选中的两个元素进行位置交换操作。
2)邻域结构N2:在工序排序部分的码段中任选两个元素,将选中的两个元素中位置靠后的元素插入到位置靠前的元素前面的的位置。
3)邻域结构N3:在及其分配部分的码段中任选一个元素,该元素的可加工设备数量大于1台,然后将该元素的对应的加工设备变为其他可选设备,进而通过方案转换个体位置的方式生成新元素值。
进行変邻域搜索:
步骤1:将当前最优个体位置X*作为变邻域初始解X,设置λ←1,qmax←3最大迭代次数ηmax
步骤2:设置q←1;
步骤3:若q=1,则X'←N1(X);若q=2,则X'←N2(X);若q=3,则X'←N3(X);
步骤4:将X'作为初始解进行局部搜索,得到最优解X”;
步骤5:若X”优于X,则X←X”,并设置q←1;否则,设置q←q+1;
步骤6:判断q>qmax是否满足,若满足,则设置q←q+1,转到步骤7;否则,转到步骤3;
步骤7:判断η>ηmax是否满足,若满足,结束;否则,转到步骤2。
邻域搜索中的局部搜索采用阈值接受法,具体如下:
步骤1:获取初始解X',并设置阈值δ>0,γ←1,ρ←1和终止条件γmax
步骤2:若ρ=1,则X”←N1(X')∪N3(X');若ρ=0,则X”←N2(X')∪N3(X');
步骤3:判断Cmax(X”)-Cmax(X')≤δ是否满足;如果满足,则X'←X”;否则,设置ρ←|ρ-1|;
步骤4:令γ←γ+1,判断γ>γmax是否满足;如果满足,则X”←X',转到步骤5;否则,转到步骤2;
步骤5:局部搜索结束。
四、算例分析
柔性作业车间调度问题的局部化和多目标优化算法(算法1)、基于禁忌搜索的柔性作业车间调度算法(算法2)中的FJSP基准Brandimarte算例(实施例1)和Kacem算例(实施例2)进行仿真实验。仿真环境为:采用MATILAB2016a程序语言,在Windows10操作系统下,配置为8G内存;CPUIntel(R)Core(TM)i5-8250U@1.60GHz主频1.80GHz的计算机上进行。算法的参数设置为:种群规模为100,螺旋常数b取1,最大迭代次数为400,领域搜索和局部搜索的最大迭代次数和设置为10。
首先分别对引入非线性收敛因子和自适应权重以及变邻域方法后算法的有效性进行对比分析,将各算法在MK01~MK10算例下独立运行20次进行比较,如表1所示,其中Best表示20次仿真得到的最优值,Avg表示20次仿真结果的平均值。WOA表示传统的鲸鱼算法,WOA-1表示在WOA中引入非线性收敛因子和自适应权重,WOA-2表示在WOA中引入变领域搜索策略,IWOA表示本文提出的算法。从表1中数据可以看出,IWOA算法的计算结果优于其他两种算法,但计算时间较其他两种算法有所延长,这是由于引入变邻域搜索机制和非线性收敛因子等,一定程度上增加了算法的时间复杂度。
为了进一步验证IWOA求解FJSP问题的有效性,将其与现有文献中的不同算法进行比较,对比结果如表2和表3所示。在表2中,Time表示对应算例的计算时间,Avg-T表示5个标准算例的平均计算时间,LB表示对应算例的已知最优解,从表2中的数据可以明显看出,IWOA在求解Kacem算例时,大多数情况下均能获得最优值。同时,基于蚁群优化算法在柔性作业车间调度问题中的应用(算法3)在配置为2.4GHzCPUand1GBRAM的计算机上对5个算例的平均计算时间214.8s,针对柔性作业车间调度问题提出的一种高效领域结构的混合禁忌搜索算法(算法4)在配置为1.6GHzCPUand512MBRAM的计算机上4978.8s的平均计算时间,混合灰狼优化算法(HGWO)求解柔性作业车间调度问题(算法5)在配置为PentiumCPUG2030@3.00GHz,2.99GHz的计算机上19s的平均计算时间,本文算法在配置为Intel(R)Core(TM)i5-8250U@1.60GHz1.80GHz的计算机上的平均计算时间仅为4s,虽然各算法的仿真平台不完全一样,但也可以看出IWOA在计算速度方面有明显的提高。
由表3可知,在求解Brandimarte的十个基准算例过程中,IWOA可以求得其中4个算例的最优解,优于针对柔性作业车间调度问题的“高尔夫社区”算法(算法8)(1个),混合灰狼优化算法(HGWO)求解柔性作业车间调度问题(算法5)(2个),一种求解柔性作业车间调度问题的有效分布式粒子群算法(算法7)(1个),基于蚁群优化算法在柔性作业车间调度问题中的应用(算法3)(3个),和基于柔性作业车间调度问题的人工免疫算法(算法6)(4个)和针对柔性作业车间调度问题提出的一种高效领域结构的混合禁忌搜索算法(算法4)(4个)相同。综合表2和表3的实验结果可知IWOA在求解以最小化最大完工时间为优化目标的FJSP方面具有一定的可行性和优越性。
表1三种算法的对比结果
表2 Kacem算例对比结果
表3 Brandimarte算例对比结果
实施例3
采用改进鲸鱼算法以一个3×6的FJSP调度问题为例进行仿真求解,其加工信息表如表4所示。其对应迭代收敛曲线如图5所示,得到的调度结果甘特图如图6所示。
表4加工信息表
本发明对柔性作业车间调度问题的特点,提出一种用于求解FJSP的改进的鲸鱼优化算法:1)算法采用两段式编码方式,设计了一种基于混沌反向学习策略和搜索法的种群初始化方法来提高初始种群的质量;其次针对当前最优个体执行变邻域搜索操作,以增强算法的局部搜索能力;此外,引入非线性收敛因子和自适应权重以平衡算法的全局搜索和局部开发的能力。2)通过对15个标准算例的仿真实验,验证了算法各种改善机制的有效性,同时也测验了算法在求解速度和精度方面的优越性。3)为下一步将把WOA进行更深层次的改进以用来求解更复杂的车间调度优化问题提供了一种新途径。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。

Claims (8)

1.一种基于改进鲸鱼算法求解柔性作业车间调度的方法,其特征在于,包括以下步骤:
步骤1:对柔性作业车间调度的机器选择和工序排序进行编码,建立柔性作业车间调度问题的数学模型;
步骤2:设置改进鲸鱼算法种群规模、最大迭代次数和位置元素,并利用混沌反向学习策略和搜索法生成初始种群;
步骤3:定义适应度函数:
其中,Fmax表示工件的最大完工时间,表示工件i的完工时间,n表示工件数;计算种群中各调度解的适应度值,得到当前最优的调度解;
步骤4:判断当前迭代次数是否大于最大迭代次数;否,则执行步骤5-8;是,则执行步骤9;
步骤5:为最优鲸鱼个体设置计数器,若改进鲸鱼算法迭代一次后最优鲸鱼个体未发生改变,则计数器值加1,判断当前最优鲸鱼个体的计数器值是否不小于预设值;是,则执行步骤6;否,则执行步骤7;
步骤6:对调度解执行变邻域搜索操作,更新调度解;
步骤7:将调度解转换为鲸鱼个体位置向量,并保留调度解对应的鲸鱼个体;
步骤8:采用改进鲸鱼算法更新鲸鱼个体位置信息;将鲸鱼个体位置向量转换为调度解S完成种群更新,令t=t+1;其中,t为当前迭代次数,转至步骤3;
步骤9:输出调度解。
2.根据权利要求1所述的基于改进鲸鱼算法求解柔性作业车间调度的方法,其特征在于,所述改进鲸鱼算法包括环绕捕食、发泡网攻击和随机搜索;其中:
环绕捕食的数学模型为:
其中:t表示当前迭代次数,是系数向量,是当前最优鲸群个体的位置向量,是当前鲸群个体的位置向量,||是取绝对值,·是逐元素相乘,在迭代过程中其内部元素非线性收敛因子a在发泡网攻击和随机搜索中从2线性递减到0,表达式为tmax为最大迭代次数;为随机向量,r取值范围为[0,1];
发泡网攻击包括收缩包围和螺旋式位置更新;收缩包围时的变动范围随着的减少而缩小,的取值范围为[-1,1];螺旋式位置更新的数学模型为:
其中,为鲸鱼最佳位置与猎物之间的距离;b为定义的控制对数螺线形状的常数,l为[-1,1]之间的随机数;同时引入自适应权重ω:
则发泡网攻击的数学模型为:
随机搜索的数学模型为:
其中, 为从当前种群中随机选择的个体位置矢量,p是[0,1]中的一个随机数。
3.根据权利要求1所述的基于改进鲸鱼算法求解柔性作业车间调度的方法,其特征在于,所述步骤1中采用等长的两段式编码方式对柔性作业车间调度的机器选择和工序排序进行编码。
4.根据权利要求1所述的基于改进鲸鱼算法求解柔性作业车间调度的方法,其特征在于,所述步骤6的具体方法为:
步骤6-1:构建邻域结构N1、N2和N3
步骤6-2:进行変邻域搜索,具体包括以下步骤:
步骤6-2-1:将当前最优个体位置X*作为变邻域初始解X,设置λ←1,qmax←3最大迭代次数ηmax
步骤6-2-2:设置q←1;
步骤6-2-3:若q=1,则X'←N1(X);若q=2,则X'←N2(X);若q=3,则X'←N3(X);
步骤6-2-4:将X'作为初始解进行局部搜索,得到最优解X”;
步骤6-2-5:若X”优于X,则X←X”,并设置q←1;否则,设置q←q+1;
步骤6-2-6:判断q>qmax是否满足,若满足,则设置q←q+1,转到步骤6-2-7;否则,转到步骤6-2-3;
步骤6-2-7:判断η>ηmax是否满足,若满足,结束;否则,转到步骤6-2-2。
5.根据权利要求4所述的基于改进鲸鱼算法求解柔性作业车间调度的方法,其特征在于,所述步骤6-1包括以下步骤:
步骤6-1-1:构建邻域结构N1:在工序排序部分的码段中任选两个元素,将选中的两个元素进行位置交换操作;
步骤6-1-2:构建邻域结构N2:在工序排序部分的码段中任选两个元素,将选中的两个元素中位置靠后的元素插入到位置靠前的元素前面的位置;
步骤6-1-3:构建邻域结构N3:在机器分配部分的码段中任选一个元素,该元素的可加工设备数量大于1台,然后将该元素的对应的加工设备变为其他可选设备,进而通过方案转换个体位置的方式生成新的排序方案。
6.根据权利要求4所述的基于改进鲸鱼算法求解柔性作业车间调度的方法,其特征在于,所述步骤6-2-4中局部搜索采用阈值接受法,包括以下步骤:
步骤6-2-4-1:获取初始解X',并设置阈值δ>0,γ←1,ρ←1和终止条件γmax
步骤6-2-4-2:若ρ=1,则X”←N1(X')∪N3(X');若ρ=0,则X”←N2(X')∪N3(X');
步骤6-2-4-3:判断Cmax(X”)-Cmax(X')≤δ是否满足;如果满足,则X'←X”;否则,设置ρ←|ρ-1|;
步骤6-2-4-4:令γ←γ+1,判断γ>γmax是否满足;如果满足,则X”←X',转到步骤6-2-4-5;否则,转到步骤6-2-4-2;
步骤6-2-4-5:局部搜索结束。
7.根据权利要求1所述的基于改进鲸鱼算法求解柔性作业车间调度的方法,其特征在于,所述将调度解S转换为鲸鱼个体位置向量的具体方法为:
S1:机器选择:采用下式将工序可选机器集中已选机器的序号转换为鲸鱼个体位置向量元素值,
x(i)=[2m/(s(i)-1)](n(i)-1)-m,s(i)≠1
其中:x(i)表示个体位置向量的第i个元素;s(i)表示元素i对应工序可以选择的机器个数;m表示机器总数;n(i)∈[1,s(i)]表示选定机器在可选机器集里面的序号;
如果s(i)=1,则x(i)在[-m,m]内任意取值;
S2:工序排序:首先生成[-m,m]内的一组与工序排序相对应的随机数,并按升序排列规则为每个随机数赋予一个唯一的ROV值,使得每个ROV值对应一个工序,然后根据工序的编码顺序对ROV值进行重排,重排后的ROV值对应的随机数顺序即为鲸鱼个体位置向量中各元素的值,组合为鲸鱼个体位置向量。
8.根据权利要求1所述的基于改进鲸鱼算法求解柔性作业车间调度的方法,其特征在于,所述将鲸鱼个体位置向量转换为调度解S的具体方法为:
R1:机器选择:按照下式
得到选定机器在机器集里面的序号;其中:xij表示i工件的第j道工序所对应的个体位置元素;sij表示i工件的第j道工序可以选择的机器个数;nij∈[1,sij]表示选定机器在机器集里面的序号;如果sij=1,则xij在[-m,m]内任取数值,2m表示个体位置向量的长度;
R2:工序排序:首先给每一个鲸鱼个体位置元素按升序的顺序赋予对应的ROV值,然后以ROV值作为元素编号,将ROV值与工序相对应构造得出相应的调度解S。
CN201910149784.7A 2019-02-28 2019-02-28 一种基于改进鲸鱼算法求解柔性作业车间调度的方法 Active CN109886588B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910149784.7A CN109886588B (zh) 2019-02-28 2019-02-28 一种基于改进鲸鱼算法求解柔性作业车间调度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910149784.7A CN109886588B (zh) 2019-02-28 2019-02-28 一种基于改进鲸鱼算法求解柔性作业车间调度的方法

Publications (2)

Publication Number Publication Date
CN109886588A true CN109886588A (zh) 2019-06-14
CN109886588B CN109886588B (zh) 2024-01-02

Family

ID=66929894

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910149784.7A Active CN109886588B (zh) 2019-02-28 2019-02-28 一种基于改进鲸鱼算法求解柔性作业车间调度的方法

Country Status (1)

Country Link
CN (1) CN109886588B (zh)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110443433A (zh) * 2019-08-15 2019-11-12 哈尔滨工业大学 基于改进鲸鱼优化算法的数字微流控芯片的调度优化方法
CN110531716A (zh) * 2019-08-12 2019-12-03 陕西科技大学 基于离散鲸鱼算法求解低碳车间调度问题的方法
CN110543151A (zh) * 2019-08-12 2019-12-06 陕西科技大学 基于改进nsga-ⅱ求解车间节能调度问题的方法
CN110782085A (zh) * 2019-10-23 2020-02-11 武汉晨曦芸峰科技有限公司 一种铸造生产调度方法及系统
CN111079987A (zh) * 2019-11-28 2020-04-28 电子科技大学 基于遗传算法的半导体车间生产调度方法
CN111967654A (zh) * 2020-07-27 2020-11-20 西安工程大学 一种基于混合遗传算法求解柔性作业车间调度的方法
CN112163808A (zh) * 2020-09-23 2021-01-01 贵州工程应用技术学院 一种基于对立学习的自适应鲸鱼算法求解物流中心选址问题的方法
CN112257922A (zh) * 2020-10-21 2021-01-22 福州大学 一种柔性作业车间调度优化方法
CN112766343A (zh) * 2021-01-12 2021-05-07 郑州轻工业大学 一种基于改进woa-svm的网络安全态势评估方法
CN112783172A (zh) * 2020-12-31 2021-05-11 重庆大学 基于离散鲸鱼优化算法的agv与机器集成调度方法
CN113960964A (zh) * 2021-09-22 2022-01-21 哈尔滨工业大学 基于仿真优化的柔性流水车间生产调度系统
CN114021745A (zh) * 2021-10-22 2022-02-08 三峡大学 一种面向水电站设备的虚拟检修流程优化方法
CN114489061A (zh) * 2022-01-17 2022-05-13 湖北工业大学 一种基于鲸鱼优化算法的平滑路径规划方法
CN114692345A (zh) * 2022-06-02 2022-07-01 河南科技学院 一种起重机臂架生产智能调度方法
CN114839929A (zh) * 2022-03-17 2022-08-02 兰州理工大学 一种用于电解铝的节能流水线调度鲸鱼优化方法及集成系统
CN114841611A (zh) * 2022-05-27 2022-08-02 盐城工学院 一种基于改进海洋捕食者算法求解作业车间调度的方法
CN114925915A (zh) * 2022-05-23 2022-08-19 广东工业大学 一种针对大规模下料的切割路径规划方法及系统
CN114936081A (zh) * 2022-05-30 2022-08-23 南京信息工程大学 一种基于woa算法的雾计算应用程序任务调度方法
CN115081754A (zh) * 2022-08-19 2022-09-20 合肥工业大学 基于混合鲸鱼-变邻域搜索的生产与维修调度方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05225203A (ja) * 1992-02-17 1993-09-03 Nippon Telegr & Teleph Corp <Ntt> ジョブショップスケジューリング問題解決方式
US20080300705A1 (en) * 2007-05-31 2008-12-04 International Business Machines Corporation Integration of job shop scheduling with discrete event simulation for manufacturing facilities
CN108021658A (zh) * 2017-12-01 2018-05-11 湖北工业大学 一种基于鲸鱼优化算法的大数据智能搜索方法及系统
CN108873835A (zh) * 2018-06-12 2018-11-23 昆明理工大学 一种半导体集成电路制造中光刻工艺的优化调度方法
CN109255426A (zh) * 2018-08-17 2019-01-22 长安大学 一种权重动态调整的地质断层参数粒子群优化搜索算法
CN109345005A (zh) * 2018-09-12 2019-02-15 中国电力科学研究院有限公司 一种基于改进鲸鱼算法的综合能源系统多维寻优方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05225203A (ja) * 1992-02-17 1993-09-03 Nippon Telegr & Teleph Corp <Ntt> ジョブショップスケジューリング問題解決方式
US20080300705A1 (en) * 2007-05-31 2008-12-04 International Business Machines Corporation Integration of job shop scheduling with discrete event simulation for manufacturing facilities
CN108021658A (zh) * 2017-12-01 2018-05-11 湖北工业大学 一种基于鲸鱼优化算法的大数据智能搜索方法及系统
CN108873835A (zh) * 2018-06-12 2018-11-23 昆明理工大学 一种半导体集成电路制造中光刻工艺的优化调度方法
CN109255426A (zh) * 2018-08-17 2019-01-22 长安大学 一种权重动态调整的地质断层参数粒子群优化搜索算法
CN109345005A (zh) * 2018-09-12 2019-02-15 中国电力科学研究院有限公司 一种基于改进鲸鱼算法的综合能源系统多维寻优方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GUOHUI ZHANG ETAL.: "Solving flexible job shop scheduling problems with transportation time based on improved genetic algorithm", 《MATHEMATICAL BIOSCIENCES AND ENGINEERING》 *
王家海等: "基于改进鱼群算法的柔性作业车间调度问题研究", 《机电一体化》 *
龙文等: "求解大规模优化问题的改进鲸鱼优化算法", 《系统工程理论与实践》 *

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110531716A (zh) * 2019-08-12 2019-12-03 陕西科技大学 基于离散鲸鱼算法求解低碳车间调度问题的方法
CN110543151A (zh) * 2019-08-12 2019-12-06 陕西科技大学 基于改进nsga-ⅱ求解车间节能调度问题的方法
CN110443433A (zh) * 2019-08-15 2019-11-12 哈尔滨工业大学 基于改进鲸鱼优化算法的数字微流控芯片的调度优化方法
CN110443433B (zh) * 2019-08-15 2023-01-17 哈尔滨工业大学 基于改进鲸鱼优化算法的数字微流控芯片的调度优化方法
CN110782085B (zh) * 2019-10-23 2022-03-29 武汉晨曦芸峰科技有限公司 一种铸造生产调度方法及系统
CN110782085A (zh) * 2019-10-23 2020-02-11 武汉晨曦芸峰科技有限公司 一种铸造生产调度方法及系统
CN111079987A (zh) * 2019-11-28 2020-04-28 电子科技大学 基于遗传算法的半导体车间生产调度方法
CN111967654A (zh) * 2020-07-27 2020-11-20 西安工程大学 一种基于混合遗传算法求解柔性作业车间调度的方法
CN112163808A (zh) * 2020-09-23 2021-01-01 贵州工程应用技术学院 一种基于对立学习的自适应鲸鱼算法求解物流中心选址问题的方法
CN112163808B (zh) * 2020-09-23 2022-04-01 贵州工程应用技术学院 一种基于对立学习的自适应鲸鱼算法求解物流中心选址问题的方法
CN112257922A (zh) * 2020-10-21 2021-01-22 福州大学 一种柔性作业车间调度优化方法
CN112257922B (zh) * 2020-10-21 2022-06-14 福州大学 一种柔性作业车间调度优化方法
CN112783172A (zh) * 2020-12-31 2021-05-11 重庆大学 基于离散鲸鱼优化算法的agv与机器集成调度方法
CN112783172B (zh) * 2020-12-31 2023-03-31 重庆大学 基于离散鲸鱼优化算法的agv与机器集成调度方法
CN112766343A (zh) * 2021-01-12 2021-05-07 郑州轻工业大学 一种基于改进woa-svm的网络安全态势评估方法
CN112766343B (zh) * 2021-01-12 2022-03-11 郑州轻工业大学 一种基于改进woa-svm的网络安全态势评估方法
CN113960964A (zh) * 2021-09-22 2022-01-21 哈尔滨工业大学 基于仿真优化的柔性流水车间生产调度系统
CN114021745A (zh) * 2021-10-22 2022-02-08 三峡大学 一种面向水电站设备的虚拟检修流程优化方法
CN114021745B (zh) * 2021-10-22 2024-07-05 三峡大学 一种面向水电站设备的虚拟检修流程优化方法
CN114489061A (zh) * 2022-01-17 2022-05-13 湖北工业大学 一种基于鲸鱼优化算法的平滑路径规划方法
CN114489061B (zh) * 2022-01-17 2023-05-16 湖北工业大学 一种基于鲸鱼优化算法的平滑路径规划方法
CN114839929A (zh) * 2022-03-17 2022-08-02 兰州理工大学 一种用于电解铝的节能流水线调度鲸鱼优化方法及集成系统
CN114839929B (zh) * 2022-03-17 2024-03-08 兰州理工大学 一种用于电解铝的节能流水线调度鲸鱼优化方法及集成系统
CN114925915A (zh) * 2022-05-23 2022-08-19 广东工业大学 一种针对大规模下料的切割路径规划方法及系统
CN114925915B (zh) * 2022-05-23 2023-03-21 广东工业大学 一种针对大规模下料的切割路径规划方法及系统
CN114841611A (zh) * 2022-05-27 2022-08-02 盐城工学院 一种基于改进海洋捕食者算法求解作业车间调度的方法
CN114936081A (zh) * 2022-05-30 2022-08-23 南京信息工程大学 一种基于woa算法的雾计算应用程序任务调度方法
CN114692345A (zh) * 2022-06-02 2022-07-01 河南科技学院 一种起重机臂架生产智能调度方法
CN115081754A (zh) * 2022-08-19 2022-09-20 合肥工业大学 基于混合鲸鱼-变邻域搜索的生产与维修调度方法

Also Published As

Publication number Publication date
CN109886588B (zh) 2024-01-02

Similar Documents

Publication Publication Date Title
CN109886588A (zh) 一种基于改进鲸鱼算法求解柔性作业车间调度的方法
CN109886589B (zh) 一种基于改进鲸鱼优化算法求解低碳车间调度的方法
CN109784603B (zh) 一种基于混合鲸鱼群算法求解柔性作业车间调度的方法
Huang et al. Survey on multi-objective evolutionary algorithms
CN109902873A (zh) 一种基于改进型鲸鱼算法的云制造资源配置的方法
CN110471419B (zh) 一种基于差分进化烟花算法的多机器人协同路径规划方法
CN104834809B (zh) 基于人工蜂群搜索的七自由度机械臂逆向运动学求解方法
CN112865195A (zh) 一种基于改进的鲸鱼算法的机组组合出力分配方法及系统
CN112577507A (zh) 基于哈里斯鹰优化算法的电动汽车路径规划方法
CN114493337A (zh) 基于改进粒子群遗传混合算法的柔性作业车间调度方法
CN116859903A (zh) 基于改进哈里斯鹰优化算法的机器人平滑路径规划方法
CN108399105B (zh) 一种基于改进头脑风暴算法的软硬件划分方法
CN113848889A (zh) 一种基于鲸鱼优化算法和人工势场法结合的路径规划方法
Thanh et al. A survey on hybridizing genetic algorithm with dynamic programming for solving the traveling salesman problem
CN114662638A (zh) 基于改进人工蜂群算法的移动机器人路径规划方法
CN110210072B (zh) 基于近似模型及差分进化算法求解高维优化问题的方法
CN113743566A (zh) 产品拆卸序列优化方法、装置、计算机设备和存储介质
CN114358414A (zh) 一种自适应并行的社区养老服务点选址优化方法及系统
CN114089699B (zh) 一种基于细胞型膜计算模型求解混合流水车间调度的方法
Ding et al. Multi-working mode product color planning using evolutionary algorithm and swarm intelligence
CN108921354A (zh) 一种基于粒子群优化的蚁群算法求解tsp问题的方法
CN115655279A (zh) 一种基于改进鲸鱼算法的海上无人救援飞艇路径规划方法
CN114565242A (zh) 一种综合低碳能源基地选址智能优化方法及装置
CN113191486A (zh) 基于参数服务器架构的图数据与参数数据的混合划分方法
CN115169754B (zh) 能源调度方法、装置、电子设备及存储介质

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
TR01 Transfer of patent right

Effective date of registration: 20240409

Address after: 710000, Chang'an Shuangchuang Center 11D-1-16, Wenyuan Middle Road, Guodu Street Office, Chang'an District, Xi'an City, Shaanxi Province

Patentee after: Xi'an Gangfeng Energy Technology Co.,Ltd.

Country or region after: China

Address before: 710064 middle section of South Second Ring Road, Beilin District, Xi'an City, Shaanxi Province

Patentee before: CHANG'AN University

Country or region before: China

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240415

Address after: Building B0806, Innovation Plaza, 2007 Pingshan Avenue, Pingshan Street, Pingshan District, Shenzhen City, Guangdong Province, 518118

Patentee after: SHENZHEN RDY TECHNOLOGY CO.,LTD.

Country or region after: China

Address before: 710000, Chang'an Shuangchuang Center 11D-1-16, Wenyuan Middle Road, Guodu Street Office, Chang'an District, Xi'an City, Shaanxi Province

Patentee before: Xi'an Gangfeng Energy Technology Co.,Ltd.

Country or region before: China

TR01 Transfer of patent right