智能算法(GA、DBO等)求解零空闲流水车间调度问题(NIFSP)-编程知识

先做一个声明：文章是由我的个人公众号中的推送直接复制粘贴而来，因此对智能优化算法感兴趣的朋友，可关注我的个人公众号：启发式算法讨论。我会不定期在公众号里分享不同的智能优化算法，经典的，或者是近几年提出的新型智能优化算法，并附MATLAB代码。

本期的主要参考资料：

[1] 潘全科, 高亮, 李新宇. 流水车间调度及其优化算法[M]. 武汉: 华中科技大学出版社, 2013.

不允许机器停止运转，这就是所谓的零空闲调度问题。图1(a)所示为一个3×3的零空闲流水车间调度问题的甘特图，由于机器连续运转的要求，当第一个工件到达机器2时，该机器并不能立即开工。与图1(b)所示的置换流水车间调度相比，机器2的开工时间有所滞后。因此，零空闲流水车间调度问题不同于一般的置换流水车间调度问题。

图1

01
问题描述

这里主要考虑以最大完成时间为优化目标的零空闲流水车间调度问题(no-idle flow-shop scheduling problem, NIFSP)，记为Fm|perm, no-idle|Cmax。当m≥3时，这类调度问题就是NP-Hard 问题。

Fm|perm, no-idle|Cmax问题可描述为：有n个工件按照相同的工艺路线在m台机器上加工，约定所有机器上工件的加工次序都相同,要求在同一机器上加工的相邻两工件之间没有空闲时间。假设机器之间存在无限大的缓冲区，一个工件不能同时由多台机器加工，一台机器也不能同时加工多个工件。已知工件在各机器上的加工时间。问题是如何安排各工件的生产次序，使得最大完成时间取值最小。

02
数学模型

（以下内容截自推文开头提到的参考书籍，潘老师的那本书。）

03
加工性能指标计算

最大完成时间(Cmax)是研究零空闲流水车间调度问题最常用的加工性能指标。NIFSP的最大完工时间有五种计算方法，其时间复杂度均为O(nm)。

方法一：计算机器之间的开工时间差

方法二：Kalezynski和Kamburowski方法(转化为F2|perm|Cmax)

方法三：前向计算法

方法四：反向计算法

方法五：双向计算法

这里主要介绍前向计算法。（以下内容截自推文开头提到的参考书籍，潘老师的那本书。）其他计算方法也可以在这本书籍里查阅。选择前向计算法是为了方便画甘特图。

04
智能算法(GA、DBO等)编码方法

对于遗传算法(GA)，因为其算法本身是离散的，通过选择、交叉、变异产生下一代。因此，一条染色体就代表一种调度方案。即工件的排序即是它的个体编码。例如，10个工件的排序方案，用MATLAB初始化GA的一个个体(一条染色体)就是：

x=randperm(10);

效果如下所示：

但是对于粒子群优化(PSO)、麻雀搜索算法(SSA)、蜣螂优化(DBO)等，它们本身是针对连续优化问题提出的，所以在编码时需要经过进一步的处理。与GA一样，一个调度方案(工件排序)表示一个个体，可以采用SPV规则，将实数编码转成整数编码。例如，10个工件的排序方案，用MATLAB初始化DBO的一个个体(一条染色体)就是：

jobNum=10; % 工件数x=unifrnd(0,1,[1 jobNum]); % 产生10个[0,1]之间随机数os = 1:1:jobNum; % 产生从1到10的数列[~, up_index] = sort(x); % 对x进行降序排序, 得到位置序列x = os(up_index); % 按照位置序列排序工件, 得到一个调度方案

效果如下：

此外，与SPV规则相反，Li等提出最大排序值法(Largest rank value, LRV)，也是将连续值映射成离散排列常用的方法之一。如图2所示，LRV将代表种群个体的一组连续值按降序排列生成一组工件排序。(参考文献：[2] LI X, YIN M. An opposition-based differential evolution algorithm for permutation flow shop scheduling based on diversity measure [J]. Advances in Engineering Software, 2013, 55(8): 10-31.)

图2 最大排序值法的表示方法

05
数值实验

这里对DBO求解NIFSP的效果进行简单测试，调度问题算例选用Rec(21个)。最大迭代次数T设置为2000，种群规模NP设为60。下面展示的结果都是算法随机运行一次得到的结果。

首先，以Rec05(20工件×5机器为例)，展示DBO随机运行一次的求解结果。图3绘制了种群每代的最优适宜度收敛曲线和平均适宜度收敛曲线：

图3 DBO-NIFSP对于Rec05的收敛曲线

图4绘制了调度结果的甘特图：

图4 DBO-NIFSP对于Rec05的甘特图

其次，以Rec11(20工件×10机器为例)，展示DBO随机运行一次的求解结果，如图5和图6所示。

图5 DBO-NIFSP对于Rec11的收敛曲线

图6 DBO-NIFSP对于Rec11的甘特图

最后，以Rec41(75工件×20机器为例)，展示DBO随机运行一次的求解结果，如图7和图8所示。

图7 DBO-NIFSP对于Rec41的收敛曲线

图8 DBO-NIFSP对于Rec41的甘特图

06
MATLAB代码

智能算法(GA、PSO、DE、GWO、SSA、DBO等)求解零空闲流水车间调度问题(no-idle flow-shop scheduling problem, NIFSP)的MATLAB代码，其中：main.m是主函数，直接运行即可；以算法简称命名的.m算法代码；gantt_chart.m用来绘制甘特图；objective.m是目标函数，即计算Makespan；method.pdf用来说明Makespan的计算方法，代码采用的是前向计算方法；Car.xlsx和Rec.xlsx是流水车间调度的两个经典测试集。

输出结果包括Makespan、工件排序、计算时间、最优适宜度收敛曲线、平均适宜度收敛曲线、甘特图。

博主选择了十种算法来求解NIFSP。主要是几种经典算法和几个近几年的高引算法。对应的MATLAB代码链接如下：

遗传算法(GA)求解NIFSP
差分进化(DE)求解NIFSP	关注公众号，里面有链接
粒子群优化(PSO)求解NIFSP	关注公众号，里面有链接
灰狼优化(GWO)求解NIFSP	关注公众号，里面有链接
鲸鱼优化算法(WOA)求解NIFSP	关注公众号，里面有链接
哈里斯鹰优化(HHO)求解NIFSP	关注公众号，里面有链接
麻雀搜索算法(SSA)求解NIFSP	关注公众号，里面有链接
非洲秃鹫优化算法(AVOA)求解NIFSP	关注公众号，里面有链接
蜣螂优化(DBO)求解NIFSP	关注公众号，里面有链接
星鸦优化算法(NOA)求解NIFSP	关注公众号，里面有链接
以上十种智能优化算法(GA、DE、PSO、GWO、WOA、HHO、SSA、AVOA、DBO、NOA)求解NIFSP的全家桶	关注公众号，里面有链接