技术领域
[0001] 本
发明涉及一种解决柔性作业车间调度问题的方法，特别是一种采用改良免疫
遗传算法的柔性作业车间调度方法。
背景技术
[0002] 柔性作业车间调度问题是经典作业车间调度问题的延伸，调度的目标是为每项作业的每道工序分配适当的机器，并确定各机器上工序的加工顺序，使作业的完工时间最短或其他指标达到最优。由于柔性作业车间的灵活性较高，使得柔性作业车间调度问题的求解需要非常大的组合搜索空间，且问题的复杂度较高。
计算机编程采用免疫遗传算法(IGA)解决柔性作业车间调度问题时，存在易陷入局部最优的缺点。
发明内容
[0003] 本发明的目的在于，提供一种采用改良免疫遗传算法的柔性作业车间调度方法。本发明将免疫遗传算法(IGA)与模拟
退火算法(SA)进行了结合，提出一种应用于车间作业调度的混合免疫遗传算法，效率更高，偏差更小。
[0004] 本发明的技术方案：采用改良免疫遗传算法的柔性作业车间调度方法，包括以下步骤，
[0005] 步骤1：建立车间作业调度问题数学模型，基于作业车间调度问题的基本信息，利用直接编码方式，将车间调度问题转化为求解遗传算法中个体优良
染色体问题，从而确定种群中个体的染色体基因序列；设置种群规模sizepop，退火初始
温度T0，退火终止温度Tend，退火系数γ，以及最大
迭代次数MAXGEN的大小；
[0006] 步骤2：初始化种群进化次数Gen，计算种群中各个个体的适应度Agv、平均适应度浓度Cv值以及期望繁殖概率ev，再对优良的个体染色体进行
破碎处理，提取其基因
片段作为
疫苗；
[0007] 步骤3：计算种群进
化成熟度M以及精英保留个数Ne，并从种群中选择适应度最高的Ne个个体进行精英保留，并对种群中一部分个体进行疫苗接种；
[0008] 步骤4：对父代种群个体进行交叉操作；
[0009] 步骤5：对临时子代种群中的个体适应度与平均适应度 比较，若个体适应度大于平均适应度，则采用换位变异方式变异；反之，则采取反转变异方式进行变异，根据自适应变异公式计算变异概率，开始变异操作；
[0010] 步骤6：模拟退火选择，通过退火初始温度T0以及退火系数γ计算当前退火温度Tk，利用Tk计算接受概率p，检测注射疫苗的
抗体的适应度并将其适应度值与旧抗体的适应度值进行比较，若新抗体适应度更高，则以抗体的繁殖概率ev选取新抗体进入下一代；若旧抗体适应度更高，接受新的抗体的概率为p*ev，更新种群；
[0011] 步骤7：判断迭代次数是否大于最大迭代次数MAXGEN，若是，则进入下一步；若不是，则返回步骤3；
[0012] 步骤8：判断当前退火温度Tk与终止温度Tend的大小，若大于终止温度，则返回步骤3；反之，则输出最优结果。
[0013] 上述的采用改良免疫遗传算法的柔性作业车间调度方法，所述步骤1的直接编码方式为：首先按照染色体上加工工序的从左到右依次降低工序的优先权重，然后安排染色体上的工序序列的第1道工序进行加工，再然后依次安排各
工件剩余工序，并将其插入到对应机器上可行的最佳加工时刻加工，最终达到所有的加工工序都安排在对应机器上最佳可行的时刻加工，得到加工时间。
[0014] 上述的采用改良免疫遗传算法的柔性作业车间调度方法，步骤3中，所述种群进化成熟度M以及精英保留个数Ne计算方法为：
[0015] 单个抗体v的适应度和浓度分别为Agv，Cv,则成熟度为：
[0016] M＝α∑Agv+(1-α)∑Cv；
[0017] 保留精英个体数的公式为：
[0018]
[0019] 式中：Ne—精英保留个体数；N—抗体总数；b—预设正整数；ceil—向上取整。
[0020] 上述的采用改良免疫遗传算法的柔性作业车间调度方法，所述步骤5中，自适应变异公式为：
[0021]
[0022] 式中：pm—变异概率；farg—群体的平均适应度值；fmax—群体最大适应度值；f—变异个体适应度值；取pc1＝0.1，pc2＝0.001。
[0023] 上述的采用改良免疫遗传算法的柔性作业车间调度方法，所述步骤6中，繁殖概率ev计算公式为：
[0024] 式中：χ—常数，Agv—抗体与
抗原的亲和
力即结合程度，Agv的计算公式为：
[0025]
[0026] 式中：Gv—抗体v对应的目标函数，Agv值越大，表明抗原与抗体的结合程度越高；Cv—抗体v浓度，Cv的计算公式为：
[0027]
[0028]
[0029] 式中：N—抗体总数，Acxy—抗体x与抗体y的相似度；T(一般取值0.9～1)—预先设定的阙值；Abxy—抗体x与抗体y之间的亲和力，其计算公式为：
[0030]
[0031] 式中：L—染色体基因的长度，l—染色体x与y的同一
位置相同编码基因的个数。
[0032] 上述的采用改良免疫遗传算法的柔性作业车间调度方法，所述步骤6过程中，概率p的公式为：
[0033]
[0034] Tk＝γk-1T0
[0035] 式中：T0—退火初始温度；γ∈(0，1)—降温系数；Tk—当前退火温度；p—Tk温度下的接受概率；Δfv—新个体适应度与旧体适应度之差；fvs—新个体适应度值；fv(s-1)—上一次迭代中旧个体适应度值。
[0036] 有益效果：本
专利算法将免疫遗传算法(IGA)与模拟退火算法(SA)进行了结合，提出一种应用于车间作业调度的混合免疫遗传算法。
[0037] 此算法设计了一种基于适应度和浓度的自适应精英保留策略且重新设置了变异算子，即将变尺度变异和自适应变异算子进行了融合。将模拟退火算法的局部搜索理论与免疫遗传算法的全局性、多样性和自适应性结合起来，互相取长补短，形成性能优良的全局寻优的混合免疫遗传算法，并设置了自适应精英选择操作，最后将其应用在作业车间调度问题的求解中。
[0038] 为验证算法的有效性与高效性，下面用Muth和Thompson提出的著名的MT06作为该算法的测试基准。
[0039] 利用该算法不同参数对MT06问题仿真试验20次，取最优的算法参数，参数设置如下：种群规模为47；交叉概率0.6～0.9，初始值为0.9；变异概率为0.01；阙值T为0.95；初始温度为10000；衰减系数为0.90；迭代次数为110。
[0040] 从表3中可见，就MT06问题而言本算法比一般的IGA，SA算法得到最优解的概率要高很多，而且平均偏差仅为2.03％；由表4可知，对于MT20问题，本专利算法比其他算法得到的结果更优，对MT06和MT10问题，虽然各个算法都能得到最优值，但本专利算法的出错率更低。综上所述，本专利算法前期能够快速收敛，后期能够跳出局部收敛，防止“早熟”，显示了其良好的全局寻优能力。
[0041] 表3基于MT06问题的不同算法的比较
[0042]
[0043] 表4不同算法性能比较
[0044]
[0045] 由上可知：
[0046] (1)本发明算法通过对MT06问题的仿真，生成了最优的调度方案。其中，最大完成时间达到了最优的55s，说明了算法的可行性；算法在第四十次就达到了收敛，如图2所示，说明了算法的高效性，图3为本发明算法对MT06问题得到的工件加工甘特图。
[0047] (2)就MT06问题而言，混合IGA算法比普通的SA算法和IGA算法击中最优解的概率更高，而且偏差更小。
[0048] (3)在求解相对复杂的MT10问题时，混合IGA算法比普通的SA算法和IGA算法离最优解更近，离最优解相差只有1秒。
附图说明
[0049] 附图1为本发明解码方法的调度方案图；
[0050] 附图2为对MT06问题的免疫遗传算法收敛曲线图；
[0051] 附图3为MT06问题得到的工件加工甘特图。
具体实施方式
[0052] 下面结合
实施例对本发明作进一步的说明，但并不作为对本发明限制的依据。
[0053] 实施例：
[0054] JSP数学描述
[0055] 车间作业调度问题(JSP)主要研究：n个工件，在m台机器上加工完成，每个工件i的工序数为oi，每个工件的不同工序在不同的机器上加工完成，其相应的加工时间也各不相同。假设：A.同一机器，可以加工多道工序，但是同一时刻最多只有一个工件在上面加工。B.每个工件的每道工序加工机器和时间已知。C.工序开始加工后，不能中断加工。D.不同工件的工序加工没有先后顺序的要求。E.同一工件在加工过程中，有预设的加工顺序。
[0056] 表1数学符号定义
[0057]
[0058] JSP问题数学模型通常可以做如下描述：
[0059] (1)目标函数：最小化最大完成时间。
[0060]
[0061] 相关约束：
[0062] (2)工艺约束下的各工件的加工顺序。同一工件必须在前一道工序加工完成后才能加工后一道工序。
[0063] cik-pik+ξ(1-aihk)≥cih (2)
[0064] (3)工件加工顺序。保证同一时刻，不同工件不会出现在同一台机器上。
[0065] cjk-pjk+ξ(1-xijk)≥cik (3)
[0066] (4)工件各工序完成时间需大于等于零，即：
[0067] cik≥0 (4)
[0068] (5)指示系数数值的分配情况。
[0069]
[0070] (6)指示变量数值分配情况。
[0071]
[0072] 求解JSP问题的改进IGA算法
[0073] 编码与解码
[0074] 考虑到直接编码方式编码和解码方式简单、柔性高和防死
锁等优点，采用一种基于工序的直接编码方式。直接编码过程就是：在求解有m台机器加工n个工件的JSP问题时，假设每个工件有m道工序，这些工序分别在m台机器上加工，所以编码的染色体组成的基因个数为n*m，每条染色体都可视为一个调度方案，染色体中同一工件的加工工序序号相同。在一条染色体中，若同一工件的加工工序序号第k次出现，则这个第k次出现的基因表示这一工件的第k道工序。例如：一个3×3的JSP问题如表2所示，可将它的一个工序序列设为[3
1 1 2 3 2 2 3 1]。其中，1，2，3分别表示工件wp1,wp2和wp3，工序序列中从左到右三个1分别表示工件wp1的第一道工序，第二道工序和第三道工序，2和3同1类似。
[0075] 表2一个3*3的JSP问题
[0076]
[0077]
[0078] 解码就是将计算机语言转化到实际问题的一种映射过程，在车间作业调度问题中，解码就是将染色体转化成一个调度方案。本专利采用的解码方法是一种插入式贪婪解码算法。贪婪解码过程为：首先按照染色体上加工工序的从左到右依次降低工序的优先权重，然后安排染色体上的工序序列的第1道工序进行加工，再然后依次安排各工件剩余工序，并将其插入到对应机器上可行的最佳加工时刻加工，最终达到所有的加工工序都安排在对应机器上最佳可行的时刻加工。例如：上述的3×3的调度问题，其中一条染色体为[3 1 1 2 3 2 2 3 1]，对应加工时间序列为[1 3 2 2 3 2 3 3 4]，机器
序列表述为[2 1 3 2 3
1 3 1 2]，该染色体通过贪婪式解码后可得到：机器1上的工件加工顺序为1-3-2，机器2上的工件加工顺序为2-1-3，机器3上的工件加工顺序为2-3-1。使用贪婪解码方法为图1所示。
由图可知使用贪婪解码的时间为10。
[0079] 精英保留策略
[0080] 在抗体评价过程中，一方面提高了适应度高的个体的选择概率，而另一方面又降低了浓度高的个体的选择概率。而在抑制高浓度的个体过程中，很容易丢失最优解。所以利用精英保留策略来保留已获得的优质的进化成果。传统的精英保留策略只保留某一固定数量适应度最高的个体，这种保留策略不利于算法的收敛，且容易导致算法陷入局部最优。基于此问题本专利提出了一种自适应精英保留策略。
[0081] 为了确定种群进化程度，引入了成熟度的概念。若单个抗体v的适应度和浓度分别为Agv，Cv,则成熟度为：
[0082] M＝α∑Agv+(1-α)∑Cv (7)
[0083] 式中：M—种群进化的成熟度；α—常数，为了加权适应度与浓度的响。
[0084] 精英保留在种群的成熟度较低的时候应该保留较多的精英个体，加快收敛速度，成熟度较高时，选择较少精英个体，提高种群的多样性。所以，其保留精英个体数的公式为：
[0085]
[0086] 式中：Ne—精英保留个体数；N—抗体总数；b—预设正整数，防止初期保留太多精英个体，导致算法陷入局部最优；ceil—向上取整。
[0087] 变异算子
[0088] 为了有效的提高算法的收敛速度和避免算法陷入局部最优，该算法将变尺度变异与自适应概率变异进行了融合，提出了混合变异算子。其中，自适应变异概率为：
[0089]
[0090] 式中：pm—变异概率；farg—群体的平均适应度值；fmax—群体最大适应度值；f—变异个体适应度值；取pc1＝0.1，pc2＝0.001。
[0091] 在算法早期，种群大多数是由随机产生的，离最优解相差甚远，种群中的大多数抗体的适应度低，大概率的采用反转变异可以提高算法在早期的查找效率。在算法后期，种群抗体的适应度高，如果也采用反转法变异将会产生过多的冗余信息。为了避免盲目搜索，提高算法执行效率，又不至于陷入局部最优，该算法采用基于自适应变异概率的小步成熟机制的换位变异法进行变异操作。
[0092] 混合变异算子的流程为：
[0093] I.计算当代种群平均适应度与需要变异抗体的适应度。
[0094] II.比较I中两适应度的大小，若变异抗体适应度大于平均适应度，选择小步成熟机制的换位法变异；若变异抗体适应度小于等于平均适应度，则采用反转变异法[0095] III.根据自适应变异公式计算变异概率，开始变异操作。
[0096] 免疫算子
[0097] 疫苗的提取与接种
[0098] 免疫操作分两步进行，即提取疫苗和接种疫苗。适当的进行疫苗注射可以有效的避免种群退化和提高种群的适应度，从而大大加快了算法的收敛速度和提高算法的寻优能力。
[0099] 本专利算法所采用的疫苗提取方法为基于加工机器的基因片段提取疫苗。
[0100] 提取疫苗的流程为：
[0101] ①将上代抗体中亲和力最大的若干个体
抽取出来。
[0102] ②将抽取出来的抗体按加工机器打碎成基因片段，将所得到的基因片段保存在疫苗库中。
[0103] 接种疫苗是利用前面提取的疫苗来
修改所选的部分个体染色体的某些位置的某些基因，从而有效的避免种群退化和提高种群的适应能力。
[0104] 疫苗接种操作如下：
[0105] 1)将父代种群的部分个体选取出来。
[0106] 2)随机选取疫苗库中的疫苗作为注射疫苗。
[0107] 3)取第1)步操作中的未注射疫苗的一个个体作为疫苗注射对象。
[0108] 4)将疫苗的基因顺序代替注射对象同一机器所对应位置基因的顺序。
[0109] 5)重复2)到4)步操作，直到第1)步操作所选取出来的个体全部注射疫苗。
[0110] 抗体评价及退火免疫选择
[0111] 本专利以抗体v的繁殖概率ev即免疫选择率作为个体的评价标准。繁殖概率是由抗体的亲和度与该抗体的浓度共同决定的。利用繁殖概率这种选择机制选择抗体既有效的鼓励了适应度高的个体，又适当的抑制了浓度高的抗体，其计算公式为：
[0112]
[0113] 式中：χ—常数，Agv—抗体与抗原的亲和力即结合程度，Agv的计算公式为：
[0114]
[0115] 式中：Gv—抗体v对应的目标函数，Agv值越大，表明抗原与抗体的结合程度越高；Cv—抗体v浓度，Cv的计算公式为：
[0116]
[0117]
[0118] 式中：N—抗体总数，Acxy—抗体x与抗体y的相似度；T(一般取值0.9～1)—预先设定的阙值；Abxy—抗体x与抗体y之间的亲和力，其计算公式为：
[0119]
[0120] 式中：L—染色体基因的长度，l—染色体x与y的同一位置相同编码基因的个数。
[0121] 为了更好地接受适应度高的抗体，并且适当的抑制浓度高的抗体，提出了基于抗体评价的模拟退火选择。
[0122] 模拟退火选择可以分为检测抗体和选择抗体两步来进行。第一步，检测抗体，检测注射疫苗的抗体的适应度并将其适应度值与旧抗体的适应度值进行比较，若新抗体适应度更高，则以抗体的繁殖概率选取新抗体进入下一代；若旧抗体适应度更高，则说明出现了种群退化现象，需要进入第二步操作。第二步是利用模拟退火的Metropolis接受准则来决定接受新抗体的概率。即①若Δfv＝(fvs-fv(s-1))≥0，则接受新抗体的概率为ev；②若Δfv＝(fvs-fv(s-1))＜0，则接受新的抗体的概率为p*ev。其中，概率p的公式为：
[0123]
[0124] Tk＝γk-1T0 (16)
[0125] 式中：T0—退火初始温度；γ∈(0，1)—降温系数；Tk—当前退火温度；p—Tk温度下的接受概率；Δfv—新个体适应度与旧体适应度之差；fvs—新个体适应度值；fv(s-1)—上一次迭代中旧个体适应度值。
[0126] 由上述分析得出：采用改良免疫遗传算法的柔性作业车间调度方法，其具体实施步骤如下：采用改良免疫遗传算法的柔性作业车间调度方法，包括以下步骤，[0127] 步骤1：建立车间作业调度问题数学模型，基于作业车间调度问题的基本信息，利用直接编码方式，将车间调度问题转化为求解遗传算法中个体优良染色体问题，从而确定种群中个体的染色体基因序列；设置种群规模sizepop，退火初始温度T0，退火终止温度Tend，退火系数γ，以及最大迭代次数MAXGEN的大小；
[0128] 步骤2：初始化种群进化次数Gen，计算种群中各个个体的适应度Agv、平均适应度浓度Cv值以及期望繁殖概率ev，再对优良的个体染色体进行破碎处理，提取其基因片段作为疫苗；
[0129] 步骤3：计算种群进化成熟度M以及精英保留个数Ne，并从种群中选择适应度最高的Ne个个体进行精英保留，并对种群中一部分个体进行疫苗接种；
[0130] 步骤4：对父代种群个体进行交叉操作；
[0131] 步骤5：对临时子代种群中的个体适应度与平均适应度 比较，若个体适应度大于平均适应度，则采用换位变异方式变异；反之，则采取反转变异方式进行变异，根据自适应变异公式计算变异概率，开始变异操作；
[0132] 步骤6：模拟退火选择，通过退火初始温度T0以及退火系数γ计算当前退火温度Tk，利用Tk计算接受概率p，检测注射疫苗的抗体的适应度并将其适应度值与旧抗体的适应度值进行比较，若新抗体适应度更高，则以抗体的繁殖概率ev选取新抗体进入下一代；若旧抗体适应度更高，接受新的抗体的概率为p*ev，更新种群；
[0133] 步骤7：判断迭代次数是否大于最大迭代次数MAXGEN，若是，则进入下一步；若不是，则返回步骤3；
[0134] 步骤8：判断当前退火温度Tk与终止温度Tend的大小，若大于终止温度，则返回步骤3；反之，则输出最优结果。
[0135] 所述步骤1的直接编码方式为：首先按照染色体上加工工序的从左到右依次降低工序的优先权重，然后安排染色体上的工序序列的第1道工序进行加工，再然后依次安排各工件剩余工序，并将其插入到对应机器上可行的最佳加工时刻加工，最终达到所有的加工工序都安排在对应机器上最佳可行的时刻加工，得到加工时间。
[0136] 步骤3中，所述种群进化成熟度M以及精英保留个数Ne计算方法为：
[0137] 单个抗体v的适应度和浓度分别为Agv，Cv,则成熟度为：
[0138] M＝α∑Agv+(1-α)∑Cv；
[0139] 保留精英个体数的公式为：
[0140]
[0141] 式中：Ne—精英保留个体数；N—抗体总数；b—预设正整数；ceil—向上取整。
[0142] 所述步骤5中，自适应变异公式为：
[0143]
[0144] 式中：pm—变异概率；farg—群体的平均适应度值；fmax—群体最大适应度值；f—变异个体适应度值；取pc1＝0.1，pc2＝0.001。
[0145] 所述步骤6中，繁殖概率ev计算公式为：
[0146] 式中：χ—常数，Agv—抗体与抗原的亲和力即结合程度，Agv的计算公式为：
[0147]
[0148] 式中：Gv—抗体v对应的目标函数，Agv值越大，表明抗原与抗体的结合程度越高；Cv—抗体v浓度，Cv的计算公式为：
[0149]
[0150]
[0151] 式中：N—抗体总数，Acxy—抗体x与抗体y的相似度；T(一般取值0.9～1)—预先设定的阙值；Abxy—抗体x与抗体y之间的亲和力，其计算公式为：
[0152]
[0153] 式中：L—染色体基因的长度，l—染色体x与y的同一位置相同编码基因的个数。
[0154] 所述步骤6过程中，概率p的公式为：
[0155]
[0156] Tk＝γk-1T0
[0157] 式中：T0—退火初始温度；γ∈(0，1)—降温系数；Tk—当前退火温度；p—Tk温度下的接受概率；Δfv—新个体适应度与旧体适应度之差；fvs—新个体适应度值；fv(s-1)—上一次迭代中旧个体适应度值。
