一种最小化期望提前和拖期费用的多机调度方法技术

本发明专利技术公开了一种最小化期望提前和拖期费用的多机调度方法，属于机器调度技术领域；以工件顺序和交付期为决策变量，并考虑了多机器的情况。对于具有对称费用结构的单机和相同多机问题，给出了最小方差优先(SVF)的调度方法并根据该规则设计了贪心算法。对于一般费用结构的问题，给出了等价二阶锥规划模型；为了实现更高效的求解，针对其中的两机问题提出了基于模函数特性的启发式调度方法，针对多机问题提出了一种变邻域的方法。本发明专利技术所提出的方法考虑加工时间的不确定性，使决策更加符合实际情况，实用性更强，降低了决策成本；能够比较快速的求出精度很高的解，非常适合求解大规模的问题，可以为实际的调度问题提供决策。

【技术实现步骤摘要】
一种最小化期望提前和拖期费用的多机调度方法
本专利技术属于机器调度
，涉及一种最小化期望提前和拖期费用的多机调度方法，具体涉及一种以最小化期望E/T费用作为衡量指标的考虑随机加工时间的多机调度方法。
技术介绍
机器调度问题是一个被广泛研究的优化问题，同时在生产调度、机器制造等方面有广泛的实际应用。不确定性是生产系统中不可避免的因素，其对系统性能有着重要影响。机器的随机波动、启动时间的变化等都使得调度者无法确定准确的生产加工时间，因而考虑随机加工时间的机器调度问题更加符合实际情况。衡量指标包括与工件完成时间相关的指标和与交付期相关的指标，准时制(JIT)生产方式促进产生了一种新的衡量指标—提前和拖期(E/T)成本指标。提前到达可能会导致额外的库存成本、产品变质等，延期到达可能会带来顾客诚意损失和失销等惩罚。机器调度问题可以分为单机器和多机器，多机器调度问题研究工件在多个机器上加工的过程，其中每个工件只能在一个机器上加工。多机器调度问题在物流、制造业和服务业都有广泛应用，因此受到广泛关注；交付期分配和加工顺序是求解多机器调度问题重点关注的方面。与单机器交付期分配与加工顺序问题相比，多个相同机器的交付期分配和加工顺序方面的研究较少。与交付期分配和加工顺序有关的调度模型大致可以分为三类：1、交付期给定的情况下，求最优的工件加工顺序；2、工件加工顺序给定的情况下，求最优交付期；3、交付期和加工顺序都是内部给定。在工件加工顺序确定的交付期分配模型中，Sridharan等针对动态交付期设置问题提出了新的交付期分配规则(SL规则)，其中工件加工时间服从的分布已知；为获取SL规则的最优参数，作者分析了先到先服务规则下的M/M/1模型。在交付期给定的工件调度模型中，可以依据加工时间是否是随机变量进一步分类；交付期给定的单机排序模型的确定性版本是调度理论的基本范式。基于该调度模型，有许多变体。大部分研究都将加权的E/T成本之和作为刻画调度有效性的目标函数，也有一些研究将加权的二次E/T成本之和作为目标。相比于加工时间确定的情况，加工时间不确定的情况更具有挑战性，这是一个NP难问题；随机的调度问题采用了不同的衡量指标来衡量结果好坏，如加权E/T成本之和、最大延迟时间的期望、加权延迟时间之和的期望等。在多机器情况下，Liu等提出了一个两阶段的随机规划公式，以最大程度地减少期望的总加权E/T成本。加工顺序和交付期联合调度问题是指为各工件指定交付期并将工件排序以使得与交付期相关的性能指标最优。在这种情况下，交付期是决策变量，可以在与客户的谈判中使用或纳入进度表中，以适应各种工作的进度。此类模型通常会将交付期的惩罚纳入目标函数中，所有工件或有一个共同的交付期，或各自拥有各自的交付期。相比较而言，不同交付期的情况更难求解。该类模型也可以进一步分为加工时间确定的情况和加工时间不确定的情况，在加工时间不确定情况下，加工时间是随机变量，通常假定其概率分布已知。Baker在假设作业的加工时间服从正态分布的前提下研究了该问题，并开发了一种分支定界算法以找到单机E/T问题的最优解。其发现对于所有工件具有相同E/T成本参数的特殊情况，最小方差优先(SVF)规则是最佳的。在多机器情况下，Gu等提出了一种基于随机规划理论的随机期望值模型，并提出了一种量子遗传调度算法；与传统的遗传算法相比，该算法能够在没有过早收敛的情况下获得更好的调度方案；Otten等引入了E/T凹面的概念，并证明了在E/T凹面的假设下，最小方差优先(SVF)规则是最优的。通过对上述文献的回顾发现，关于加工时间不确定情况下的多机排序和交付期联合调度方面的研究较少，且难以实现高效求解。
技术实现思路
针对现有技术中存在的上述问题，本专利技术提供最小化期望提前和拖期费用的多机调度方法，其以最小化期望E/T成本为衡量指标、以加工顺序和交付期为决策变量；针对两机器和多机器的情况，分别提出了考虑模函数特性的启发式调度方法和一种变邻域的方法。本专利技术公开了一种最小化期望提前和拖期费用的多机调度方法，包括：对于一般费用结构的问题，采用等价二阶锥规划模型求解方法；对于具有对称费用结构的单机和相同多机问题，采用最小方差优先的调度方法。作为本专利技术的进一步改进，所述采用等价二阶锥规划模型的求解方法，包括：在原优化问题的基础上，基于各原优化问题转化成等价的混合0-1二阶锥规化模型，求解该模型获得精确解；其中，变量xk，ij∈{0，1}表示为工件i和j是否都在机器k上加工。作为本专利技术的进一步改进，所述采用等价二阶锥规划模型的求解方法，具体包括：设模型中有n个工件，工件集合为J＝{1，…，n}；有m个机器，机器集合为M＝{1，…，m}，工件i在机器k上的随机加工时间由pk，i表示，pk，i的方差由表示；变量xk，ij∈{0，1}表示为工件i和j是否都在机器k上加工，其中i≠j∈J，k∈M；变量xk，ii∈{0，1}表示为工件i是否在机器k上加工，其中i∈J，k∈M；用Jk，i表示在机器k上方差小于等于的工件集合；所考虑的优化问题的目标函数表述如下：(P1)模型的约束条件如下：将优化问题转化为如下等价的混合0-1二阶锥规化模型，目标函数具有以下形式：(P2)min∑k∈K∑i∈Jyk，i(16)(P2)模型的约束条件如下：xk，ij≥xk，ii+xk，jj-1(19)通过求解该模型获得精确解。作为本专利技术的进一步改进，还包括：对于具有一般费用结构的两机问题，采用基于模函数特性的启发式调度方法；对于具有一般费用结构的多机问题，采用变邻域的方法。作为本专利技术的进一步改进，所述采用基于模函数特性的启发式调度方法，包括：考虑工件集合N＝{1，…，n}的两机器问题，为分配到机器1上的工件，因此分配到机器2上的工件集合为N\S；fi(S)为当S被分配到机器i，与机器i相关的成本，f(S)＝f1(S)+f2(N\S)，所考虑的两机器问题可以被表述为用S*表示最优解；f(S)是一个超模函数，-f(S)是一个次模函数；次模函数-f(S)有如下命题：若且对于某个满足则即，工件k会被最优分配给机器2；若且对于某个k∈S满足则k∈S*；即，工件k会被最优分配给机器1。作为本专利技术的进一步改进，所述变邻域的方法包括：插入方法、两机器交换方法和全局交换方法。作为本专利技术的进一步改进，所述插入方法为：假设工件i当前在机器k上加工，计算当工件i在其他机器上加工时目标函数值是否减少，如果插入到某一机器使得目标函数值减小，则将工件i插入到该机器上；所述两机器交换方法为：假设工件i，j当前分别在机器p，q上加工，计算当工件i在机器q上加工，工件j在机器p上加工时目标函数值是否减少，如果函数值减小，则交换加工工件i和j的机器；所述全局交换方本文档来自技高网...

【技术保护点】
1.一种最小化期望提前和拖期费用的多机调度方法，其特征在于，包括：/n对于一般费用结构的问题，采用等价二阶锥规划模型求解方法；/n对于具有对称费用结构的单机和相同多机问题，采用最小方差优先的调度方法。/n

【技术特征摘要】
1.一种最小化期望提前和拖期费用的多机调度方法，其特征在于，包括：
对于一般费用结构的问题，采用等价二阶锥规划模型求解方法；
对于具有对称费用结构的单机和相同多机问题，采用最小方差优先的调度方法。


2.如权利要求1的多机调度方法，其特征在于，所述采用等价二阶锥规划模型的求解方法，包括：
在原优化问题的基础上，基于将原优化问题转化成等价的混合0-1二阶锥规化模型，求解该模型获得精确解；其中，变量xk,ij∈{0,1}表示为工件i和j是否都在机器k上加工。


3.如权利要求2所述的多机调度方法，其特征在于，所述采用等价二阶锥规划模型的求解方法，具体包括：
设模型中有n个工件，工件集合为J＝{1,…,n}；有m个机器，机器集合为M＝{1,…,m}，工件i在机器k上的随机加工时间由pk,i表示，pk,i的方差由表示；变量xk,ij∈{0,1}表示为工件i和j是否都在机器k上加工，其中i≠j∈J,k∈M；变量xk,ii∈{0,1}表示为工件i是否在机器k上加工，其中i∈J,k∈M；用Jk,i表示在机器k上方差小于等于的工件集合；
所考虑的优化问题的目标函数表述如下：

(P1)模型的约束条件如下：












将优化问题转化为如下等价的混合0-1二阶锥规化模型，目标函数如下：
(P2)min∑k∈K∑i∈Jyk,i(6)
(P2)模型的约束条件如下：






xk,ij≥xk,ii+xk,jj-1(9)



通过求解该模型获得精确解。


4.如权利要求1～3中任一项所述的多机调度方法，其特...

【专利技术属性】
技术研发人员：张玉利，周超越，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京;11