CN113408951B - 一种基于动态信息累计狮群最优柔性调度方法及系统 - Google Patents

一种基于动态信息累计狮群最优柔性调度方法及系统 Download PDF

Info

Publication number
CN113408951B
CN113408951B CN202110806542.8A CN202110806542A CN113408951B CN 113408951 B CN113408951 B CN 113408951B CN 202110806542 A CN202110806542 A CN 202110806542A CN 113408951 B CN113408951 B CN 113408951B
Authority
CN
China
Prior art keywords
optimal
population
individual
lion
formula
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.)
Active
Application number
CN202110806542.8A
Other languages
English (en)
Other versions
CN113408951A (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.)
Shandong University
Original Assignee
Shandong 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 Shandong University filed Critical Shandong University
Priority to CN202110806542.8A priority Critical patent/CN113408951B/zh
Publication of CN113408951A publication Critical patent/CN113408951A/zh
Application granted granted Critical
Publication of CN113408951B publication Critical patent/CN113408951B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06311Scheduling, planning or task assignment for a person or group
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial 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]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/04Manufacturing
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (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

本发明涉及一种基于动态信息累计狮群最优柔性调度方法及系统,包括:(1)以最高生产工时效益为目标,构建柔性车间调度的数学优化模型,确定优化目标为最大加工完成时间,确定加工工件和工序数量、加工机器数量、工序用时等先验参数;(2)使用动态信息累计狮群算法优化数学优化模型;初始化狮群算法种群,在种群个体和柔性调度方案之间建立映射关系;在种群更新迭代过程中,加入个体最优信息累计因子和种群最优信息累计因子;本发明可以得到用时更短的调度方案,节省生产时间,提高了生产效率,同时该所提算法对不同类型的调度问题具有很好的适应性和鲁棒性。

Description

一种基于动态信息累计狮群最优柔性调度方法及系统
技术领域
本发明涉及一种基于动态信息累计狮群最优柔性调度方法,属于群智能优化技术领域。
背景技术
生产调度问题广泛存在于工业生产领域,尤其是在各种制造企业中。工业生产的调度问题的研究主要集中在生产资源的使用和分配上。通常来说,针对不同类型和需求的生产加工任务,需要制定相应的生产调度方案和规划。制定方案通常有一个或多个优化目标作为依据,例如:生产时间,机器能耗,机器利用率等等。
随着科技水平和工业生产力的进步,调度问题变得越来越复杂化和规模化。从简单的置换流水线调度问题,到车间生产调度问题,再到柔性调度问题,乃至于各种特定限制条件和要求的调度问题都在不断地提出。而且生产车间调度问题已被证明是一种NP难题,找到它的最优解是很困难的。这对于这类问题的解决方法提出了很高的要求和严峻考验。
传统的数学规划算法由于计算复杂度高,鲁棒性较低,对调度问题的约束条件有高要求等原因,在解决大规模调度问题时常常很难获得令人满意的解决方案。
在这种情况下,狮群优化算法作为最新的群智能优化算法之一,也是一种元启发式算法,越来越受到工业生产调度领域的广泛关注。因为它具有高效性,智能性,通用性,鲁棒性好,适应性强的优点,所以渐渐成为解决调度问题的一种新的有效方法。它具有结构较为简单,收敛速度快,多样性好的优点,是一种比较优秀的优化算法。目前而言,狮群算法还存在一定的问题,如较易陷入局部极值等,还具有较大的改进空间和发展潜力。
发明内容
针对现有技术的不足,本发明提供了一种基于动态信息累计狮群最优柔性调度方法。通过在狮群算法的迭代更新机制中,增加了历史时序的累计信息,所采取措施有助于对于整个问题寻优空间进行彻底搜索,找到全局最优解,有效避开降低陷入局部最优的情况。
本发明还提供了一种基于动态信息累计狮群最优柔性调度系统。
术语解释:
1、群智能算法,是一种受自然生物群启发的新兴演化计算技术,灵感和理论来源是自然界的狮群相互协作,共同合作的捕猎过程。通过对狮子群体的组成成分和行为活动的观察分析,进而提出。
2、狮群种群,整个狮群在捕猎时分工明确,主要划分成狮王,母狮和幼狮三类群体。狮群优化算法的主要思想如下:从待寻优区域中的某一初始位置开始,其中具有最优目标函数值,也就是最佳适应度值的狮子确定为狮王,然后依据一定的比例,确定母狮,母狮相互配合寻优捕猎。如果发现比当前狮王位置更优的位置,该位置就会被狮王拥有。幼狮跟随母狮移动,学习捕猎,或在狮王附近活动、进食,成年后会被驱赶出狮群。整个狮群分工合作,按照自己的职责,不断重复搜寻,得出目标函数最优值。
3、种群个体,现假设一个优化问题f(x),它的解空间维度是D。狮群的种群数量为N。狮群的每一只狮子可以用一个D维的向量来表示。xi=(xi1,xi2,...,xiD),1≤i≤N。
本发明的技术方案为:
一种基于动态信息累计狮群最优柔性调度方法,包括步骤如下:
(1)以最高生产工时效益为目标,构建柔性车间调度的数学优化模型,确定优化目标为最大加工完成时间,确定先验参数,包括:工件数量n、机器数量m、加工工件编号和以及各自工序数量、加工机器编号和数量、每个工序的用时和对应机器编号;
(2)使用动态信息累计狮群算法优化数学优化模型;
初始化狮群算法种群,在种群个体和柔性调度方案之间建立映射关系;
在种群更新迭代过程中,加入个体最优信息累计因子和种群最优信息累计因子;在原来的基础上增加了两个自由度,进而增加了个体运动方向的多样性,可以彻底探索整个问题寻优空间。
最终输出狮子的最优位置和最优值,最优值便是优化后的最大加工完成时间,最优位置经过映射转换得到加工调度方案。
根据本发明优选的,步骤(1)中,柔性车间调度的数学优化模型为柔性车间调度模型,构建柔性车间调度模型,确定关键参数;
关键参数包括机器顺序矩阵M和加工时间矩阵T,机器顺序矩阵M中的任一元素Mi,j表示加工工件i的第j道工序的机器编号;加工时间矩阵T中的任一元素Ti,j为工件i的第j道工序需要的加工时间;机器顺序矩阵M和加工时间矩阵T两者均为n×m矩阵;
该柔性车间调度模型的优化目标函数如式(Ⅰ)所示:
Figure BDA0003166570230000021
式(Ⅰ)中,Cmax是加工完成时间,C(πi,k)代表编号为i的工件在编号为k的机器上的加工完毕时间;πi表示一个工件,i是它的编号,i=1,2,...,n;
C(πi,k)=max{C(πi-1,k),C(πi,k-1)}+p(πi,k),i=2,...,n;k=2,...,m,C(πi-1,k)是指编号为i-1的工件在编号为k的机器上的加工完毕时间,C(πi,k-1)是指编号为i的工件在编号为k-1的机器上的加工完毕时间;
该柔性车间调度模型的优化目为最小化加工完成时间:min{Cmax};
本方法中,将柔性车间调度中的机器设备选择和加工序列确定这两步合并为一步完成,这就需要以双层基因编码的方式确定种群个体。
柔性车间调度中,选择机器设备和确定加工序列:
以双层基因编码的方式确定种群个体,如式(Ⅱ)所示:
D=2×Np (Ⅱ)
式(Ⅱ)中,D是指种群个体维度,Np是指工件的工序总数;
在本发明中,加工工件和工序数量、加工机器数量、工序用时这些先验参数来自FJSP标准算例模型Brandimarte_MK系列。在实际生产中,则需要将待加工的工件依照工序时间先后进行步骤分解,一一确认。
根据本发明优选的,步骤(2)的具体实现过程包括:
A、初始化狮群种群,一个种群个体通过半活动解码算法映射为一种调度方案;
B、计算每个种群个体的适应度值,初始化种群个体最优和全局最优;
C、依照狮群算法的机制进行种群的更新迭代,在迭代计算中,加入了信息累计因子和动态调整学习方法。
根据本发明优选的,步骤A的具体实现过程包括:
定义一个种群个体Xi=[xi,1,xi,2,...,xi,n×m],为长度为n×m的实数串;
利用基于随机键编码的LOV(largest-order-value)规则实现从种群个体Xi=[xi,1,xi,2,...,xi,n×m]到离散加工排序πi={πi,1i,2,...,πi,n×m}的转化,具体是指:将Xi=[xi,1,xi,2,...,xi,n×m]按各个分量的数值大小升序排列进行标号,标号是1到n×m的整数,得到标号的中间序列
Figure BDA0003166570230000031
Figure BDA0003166570230000032
采用半活动解码算法映射确定调度方案,具体包括步骤如下:
步骤1:定义计数变量j=1,和计数向量k(ω)=1,ω=1,2,...n,g(ω)=1,ω=1,2,...m;
步骤2:依据加工工件πj的当前加工工序k(πj),确定对应的机器编号
Figure BDA0003166570230000041
步骤3:比较加工工件πj的上一道工序k(πj)-1的完成时间t1和机器
Figure BDA0003166570230000042
的上一个加工任务
Figure BDA0003166570230000043
完成时间t2;若k(πj)=1,则t1=0;若
Figure BDA0003166570230000044
则t2=0;
步骤4:
Figure BDA0003166570230000045
步骤5:令k(πj)=k(πj)+1,
Figure BDA0003166570230000046
步骤6:令j=j+1,若j≤n×m,转到步骤2,否则,结束。
根据本发明优选的,步骤B的具体实现过程包括:种群个体通过步骤A映射为一种调度方案,再根据式(Ⅰ)得到最大完成时间,也就是该种群个体的适应度值,将求取的该种群个体的适应度值即第一代狮群个体的适应度值初始化为每个种群个体最优,以所有种群个体最小的适应度值确定为全局最优。
根据本发明优选的,步骤C中,依照狮群算法的机制进行种群的更新迭代,更新迭代如式(Ⅲ)所示:
Figure BDA0003166570230000047
式(Ⅲ)中,
Figure BDA0003166570230000048
代表狮王自身更新获得的子代全局最优个体,gk为第k代种群的全局最优个体,
Figure BDA0003166570230000049
代表经历过k次迭代后第i个个体的历史最优位置,γ是依据正态分布产生的一个随机数;
动态调整学习方法是在狮子两两合作捕猎的过程中,通过权重因子的调整来加强较优个体的影响力,如式(Ⅳ)所示:
Figure BDA00031665702300000410
式(Ⅳ)中,K为(0,1)的削弱因子,起到削弱较差个体和增强优势个体的作用;
Figure BDA00031665702300000411
代表公式更新后的狮子个体,
Figure BDA00031665702300000412
是第i个狮子第k次迭代的历史最优位置,
Figure BDA00031665702300000413
是除了第i个狮子之外的另一个母狮子的历史最优位置,扰动因子αf,它是迭代次数t的动态函数,随着t的增加而减小,而逐渐趋于零;
通过式(Ⅳ)中的判断条件动态调整和选择,是指:先比较
Figure BDA0003166570230000051
Figure BDA0003166570230000052
的适应值,在适应值较大的位置(即离最优值距离较远的位置)上添加削弱因子K,以削弱较差位置个体的影响权重,増强较优位置个体的影响权重;
信息累计因子包括个体信息累计因子和全局信息累计因子,个体信息累计因子Cp和全局信息累计因子Cg分别如式(Ⅴ)、式(Ⅵ)所示:
Figure BDA0003166570230000053
Figure BDA0003166570230000054
式(Ⅴ)、式(Ⅵ)中,pi指的是第i个个体的历史最优位置,pg是整个种群的全局最优位置,k是当前迭代次数;
加入了信息累计因子和动态调整学习方法后的更新迭代如式(Ⅶ)所示:
Figure BDA0003166570230000055
根据本发明优选的,步骤C的具体实现过程包括:
1)计算新种群个体的适应度值,更新个体最优值和全局最优值;是指:
Figure BDA0003166570230000056
代表式(Ⅶ)更新后的狮子个体,
Figure BDA0003166570230000057
是第i个狮子第k次迭代的历史最优位置,若:
Figure BDA0003166570230000058
则令
Figure BDA0003166570230000059
此为新的个体最优;若存在:
Figure BDA00031665702300000510
pg k+1为新的全局最优。
2)根据新的个体最优值和全局最优值,更新个体信息累计因子和全局信息累计因子,是指:新的个体最优值pi和新的全局最优值pg,带入式(Ⅴ)、式(Ⅵ),计算新的Cp和Cg
3)判断是否达到迭代停止条件,若是,进行步骤4);否则,返回步骤1);迭代次数达到预设最大值T;
4)输出优化后的调度方案和目标函数值(加工完成用时);是指:输出最后一次迭代得到的pg,通过半活动解码算法映射为一种调度方案,即优化后的调度方案,该调度方案的加工完成时间就是目标函数值。
一种基于动态信息累计狮群最优柔性调度系统,用于实现基于动态信息累计狮群最优柔性调度方法,包括数学优化模型构建单元及数学优化模型优化单元;所述数学优化模型构建单元用于实现所述步骤(1),所述数学优化模型优化单元用于实现所述步骤(2)。
本发明的有益效果是:
在处理相关柔性调度问题时,本发明使用了一种新型动态信息累计狮群优化算法,可以得到用时更短的调度方案,节省生产时间,提高了生产效率,同时该所提算法对不同类型的调度问题具有很好的适应性和鲁棒性。
附图说明
图1为本发明基于动态信息累计狮群最优柔性调度方法的流程示意框图;
具体实施方式
下面结合附图和实施例对本发明作进一步说明,但不限于此。
实施例1
一种基于动态信息累计狮群最优柔性调度方法,如图1所示,包括步骤如下:
(1)以最高生产工时效益为目标,构建柔性车间调度的数学优化模型,确定优化目标为最大加工完成时间,确定先验参数,包括:工件数量n、机器数量m、加工工件编号和以及各自工序数量、加工机器编号和数量、每个工序的用时和对应机器编号;
(2)使用动态信息累计狮群算法优化数学优化模型;
初始化狮群算法种群,在种群个体和柔性调度方案之间建立映射关系;
在种群更新迭代过程中,加入个体最优信息累计因子和种群最优信息累计因子;在原来的基础上增加了两个自由度,进而增加了个体运动方向的多样性,可以彻底探索整个问题寻优空间。
最终输出狮子的最优位置和最优值,最优值便是优化后的最大加工完成时间,最优位置经过映射转换得到加工调度方案。
实施例2
根据实施例1所述的一种基于动态信息累计狮群最优柔性调度方法,其区别在于:
步骤(1)中,柔性车间调度的数学优化模型为柔性车间调度模型,构建柔性车间调度模型,确定关键参数;
关键参数包括机器顺序矩阵M和加工时间矩阵T,机器顺序矩阵M中的任一元素Mi,j表示加工工件i的第j道工序的机器编号;加工时间矩阵T中的任一元素Ti,j为工件i的第j道工序需要的加工时间;机器顺序矩阵M和加工时间矩阵T两者均为n×m矩阵;
该柔性车间调度模型的优化目标函数如式(Ⅰ)所示:
Figure BDA0003166570230000071
式(Ⅰ)中,Cmax是加工完成时间,C(πi,k)代表编号为i的工件在编号为k的机器上的加工完毕时间;πi表示一个工件,i是它的编号,i=1,2,...,n;
C(πi,k)=max{C(πi-1,k),C(πi,k-1)}+p(πi,k),i=2,...,n;k=2,...,m,C(πi-1,k)是指编号为i-1的工件在编号为k的机器上的加工完毕时间,C(πi,k-1)是指编号为i的工件在编号为k-1的机器上的加工完毕时间;
该柔性车间调度模型的优化目为最小化加工完成时间:min{Cmax};
本方法中,将柔性车间调度中的机器设备选择和加工序列确定这两步合并为一步完成,这就需要以双层基因编码的方式确定种群个体。
柔性车间调度中,选择机器设备和确定加工序列:
以双层基因编码的方式确定种群个体,如式(Ⅱ)所示:
D=2×Np (Ⅱ)
式(Ⅱ)中,D是指种群个体维度,Np是指工件的工序总数;
在本发明中,加工工件和工序数量、加工机器数量、工序用时这些先验参数来自FJSP标准算例模型Brandimarte_MK系列。在实际生产中,则需要将待加工的工件依照工序时间先后进行步骤分解,一一确认。
步骤(2)的具体实现过程包括:
A、初始化狮群种群,一个种群个体通过半活动解码算法映射为一种调度方案;
B、计算每个种群个体的适应度值,初始化种群个体最优和全局最优;
C、依照狮群算法的机制进行种群的更新迭代,在迭代计算中,加入了信息累计因子和动态调整学习方法。
步骤A的具体实现过程包括:
定义一个种群个体Xi=[xi,1,xi,2,...,xi,n×m],为长度为n×m的实数串;
利用基于随机键编码的LOV(largest-order-value)规则实现从种群个体Xi=[xi,1,xi,2,...,xi,n×m]到离散加工排序πi={πi,1i,2,...,πi,n×m}的转化,具体是指:将Xi=[xi,1,xi,2,...,xi,n×m]按各个分量的数值大小升序排列进行标号,标号是1到n×m的整数,得到标号的中间序列
Figure BDA0003166570230000072
Figure BDA0003166570230000081
采用半活动解码算法映射确定调度方案,具体包括步骤如下:
步骤1:定义计数变量j=1,和计数向量k(ω)=1,ω=1,2,...n,g(ω)=1,ω=1,2,...m;
步骤2:依据加工工件πj的当前加工工序k(πj),确定对应的机器编号
Figure BDA0003166570230000082
步骤3:比较加工工件πj的上一道工序k(πj)-1的完成时间t1和机器
Figure BDA0003166570230000083
的上一个加工任务
Figure BDA0003166570230000084
完成时间t2;若k(πj)=1,则t1=0;若
Figure BDA0003166570230000085
则t2=0;
步骤4:
Figure BDA0003166570230000086
步骤5:令k(πj)=k(πj)+1,
Figure BDA0003166570230000087
步骤6:令j=j+1,若j≤n×m,转到步骤2,否则,结束。
步骤B的具体实现过程包括:种群个体通过步骤A映射为一种调度方案,再根据式(Ⅰ)得到最大完成时间,也就是该种群个体的适应度值,将求取的该种群个体的适应度值即第一代狮群个体的适应度值初始化为每个种群个体最优,以所有种群个体最小的适应度值确定为全局最优。
步骤C中,依照狮群算法的机制进行种群的更新迭代,更新迭代如式(Ⅲ)所示:
Figure BDA0003166570230000088
式(Ⅲ)中,
Figure BDA0003166570230000089
代表狮王自身更新获得的子代全局最优个体,gk为第k代种群的全局最优个体,
Figure BDA00031665702300000810
代表经历过k次迭代后第i个个体的历史最优位置,γ是依据正态分布产生的一个随机数;
动态调整学习方法是在狮子两两合作捕猎的过程中,通过权重因子的调整来加强较优个体的影响力,如式(Ⅳ)所示:
Figure BDA00031665702300000811
式(Ⅳ)中,K为(0,1)的削弱因子,起到削弱较差个体和增强优势个体的作用;
Figure BDA00031665702300000812
代表公式更新后的狮子个体,
Figure BDA00031665702300000813
是第i个狮子第k次迭代的历史最优位置,
Figure BDA00031665702300000814
是除了第i个狮子之外的另一个母狮子的历史最优位置,扰动因子αf,它是迭代次数t的动态函数,随着t的增加而减小,而逐渐趋于零;
通过式(Ⅳ)中的判断条件动态调整和选择,是指:先比较
Figure BDA0003166570230000091
Figure BDA0003166570230000092
的适应值,在适应值较大的位置(即离最优值距离较远的位置)上添加削弱因子K,以削弱较差位置个体的影响权重,増强较优位置个体的影响权重;
信息累计因子包括个体信息累计因子和全局信息累计因子,个体信息累计因子Cp和全局信息累计因子Cg分别如式(Ⅴ)、式(Ⅵ)所示:
Figure BDA0003166570230000093
Figure BDA0003166570230000094
式(Ⅴ)、式(Ⅵ)中,pi指的是第i个个体的历史最优位置,pg是整个种群的全局最优位置,k是当前迭代次数;
加入了信息累计因子和动态调整学习方法后的更新迭代如式(Ⅶ)所示:
Figure BDA0003166570230000095
步骤C的具体实现过程包括:
1)计算新种群个体的适应度值,更新个体最优值和全局最优值;是指:
Figure BDA0003166570230000096
代表式(Ⅶ)更新后的狮子个体,
Figure BDA0003166570230000097
是第i个狮子第k次迭代的历史最优位置,若:
Figure BDA0003166570230000098
则令
Figure BDA0003166570230000099
此为新的个体最优;若存在:
Figure BDA00031665702300000910
pg k+1为新的全局最优。
2)根据新的个体最优值和全局最优值,更新个体信息累计因子和全局信息累计因子,是指:新的个体最优值pi和新的全局最优值pg,带入式(Ⅴ)、式(Ⅵ),计算新的Cp和Cg
3)判断是否达到迭代停止条件,若是,进行步骤4);否则,返回步骤1);迭代次数达到预设最大值T;
4)输出优化后的调度方案和目标函数值(加工完成用时);是指:输出最后一次迭代得到的pg,通过半活动解码算法映射为一种调度方案,即优化后的调度方案,该调度方案的加工完成时间就是目标函数值。
表1在mk系列的三个标准算例上对比了传统规划算法和基于动态信息累计狮群最优柔性调度方法的优化结果。表格统计了20次重复实验的测试数据,并进行分析。
表1
Figure BDA0003166570230000101
结果表明,在这些算例模型上,本发明的方法对比传统算法,在各项性能上均有所提高。说明基于动态信息累计狮群最优柔性调度方法在处理该类调度问题上具有优越性。
实施例3
一种基于动态信息累计狮群最优柔性调度系统,用于实现实施例1或2所述的基于动态信息累计狮群最优柔性调度方法,包括数学优化模型构建单元及数学优化模型优化单元;数学优化模型构建单元用于实现步骤(1),数学优化模型优化单元用于实现步骤(2)。

Claims (5)

1.一种基于动态信息累计狮群最优柔性调度方法,其特征在于,包括步骤如下:
(1)以最高生产工时效益为目标,构建柔性车间调度的数学优化模型,确定优化目标为最大加工完成时间,确定先验参数,包括:工件数量n、机器数量m、加工工件编号和以及各自工序数量、加工机器编号和数量、每个工序的用时和对应机器编号;
(2)使用动态信息累计狮群算法优化数学优化模型;
初始化狮群算法种群,在种群个体和柔性调度方案之间建立映射关系;
在种群更新迭代过程中,加入个体最优信息累计因子和种群最优信息累计因子;
最终输出狮子的最优位置和最优值,最优值便是优化后的最大加工完成时间,最优位置经过映射转换得到加工调度方案;
步骤(1)中,柔性车间调度的数学优化模型为柔性车间调度模型,构建柔性车间调度模型,确定关键参数;
关键参数包括机器顺序矩阵M和加工时间矩阵T,机器顺序矩阵M中的任一元素Mi,j表示加工工件i的第j道工序的机器编号;加工时间矩阵T中的任一元素Ti,j为工件i的第j道工序需要的加工时间;机器顺序矩阵M和加工时间矩阵T两者均为n×m矩阵;
该柔性车间调度模型的优化目标函数如式(Ⅰ)所示:
Figure FDA0003718382280000011
式(Ⅰ)中,Cmax是加工完成时间,C(πi,k)代表编号为i的工件在编号为k的机器上的加工完毕时间;πi表示一个工件,i是它的编号,i=1,2,...,n;
C(πi,k)=max{C(πi-1,k),C(πi,k-1)}+p(πi,k),i=2,...,n;k=2,...,m,C(πi-1,k)是指编号为i-1的工件在编号为k的机器上的加工完毕时间,C(πi,k-1)是指编号为i的工件在编号为k-1的机器上的加工完毕时间;
该柔性车间调度模型的优化目为最小化加工完成时间:min{Cmax};
柔性车间调度中,选择机器设备和确定加工序列:
以双层基因编码的方式确定种群个体,如式(II)所示:
D=2×Np (II)
式(II)中,D是指种群个体维度,Np是指工件的工序总数;
步骤(2)的具体实现过程包括:
A、初始化狮群种群,一个种群个体通过半活动解码算法映射为一种调度方案;
B、计算每个种群个体的适应度值,初始化种群个体最优和全局最优;
C、依照狮群算法的机制进行种群的更新迭代,在迭代计算中,加入了信息累计因子和动态调整学习方法;
步骤C中,依照狮群算法的机制进行种群的更新迭代,更新迭代如式(Ⅲ)所示:
Figure FDA0003718382280000021
式(Ⅲ)中,
Figure FDA0003718382280000022
代表狮王自身更新获得的子代全局最优个体,gk为第k代种群的全局最优个体,
Figure FDA0003718382280000023
代表经历过k次迭代后第i个个体的历史最优位置,γ是依据正态分布产生的一个随机数;
动态调整学习方法是在狮子两两合作捕猎的过程中,通过权重因子的调整来加强较优个体的影响力,如式(Ⅳ)所示:
Figure FDA0003718382280000024
式(Ⅳ)中,K为(0,1)的削弱因子,起到削弱较差个体和增强优势个体的作用;
Figure FDA0003718382280000025
代表公式更新后的狮子个体,
Figure FDA0003718382280000026
是第i个狮子第k次迭代的历史最优位置,
Figure FDA0003718382280000027
是除了第i个狮子之外的另一个母狮子的历史最优位置,扰动因子αf,它是迭代次数t的动态函数,随着t的增加而减小,而逐渐趋于零;
通过式(Ⅳ)中的判断条件动态调整和选择,是指:先比较
Figure FDA0003718382280000028
Figure FDA0003718382280000029
的适应值,在适应值较大的位置上添加削弱因子K,以削弱较差位置个体的影响权重,増强较优位置个体的影响权重;
信息累计因子包括个体信息累计因子和全局信息累计因子,个体信息累计因子Cp和全局信息累计因子Cg分别如式(V)、式(Ⅵ)所示:
Figure FDA00037183822800000210
Figure FDA0003718382280000031
式(V)、式(Ⅵ)中,pi指的是第i个个体的历史最优位置,pg是整个种群的全局最优位置,k是当前迭代次数;
加入了信息累计因子和动态调整学习方法后的更新迭代如式(Ⅶ)所示:
Figure FDA0003718382280000032
2.根据权利要求1所述的一种基于动态信息累计狮群最优柔性调度方法,其特征在于,步骤A的具体实现过程包括:
定义一个种群个体Xi=[xi,1,xi,2,...,xi,n×m],为长度为n×m的实数串;
利用基于随机键编码的LOV规则实现从种群个体
Figure FDA0003718382280000033
到离散加工排序πi={πi,1i,2,...,πi,n×m}的转化,具体是指:将Xi=[xi,1,xi,2,...,xi,n×m]按各个分量的数值大小升序排列进行标号,标号是1到n×m的整数,得到标号的中间序列
Figure FDA0003718382280000034
Figure FDA0003718382280000035
采用半活动解码算法映射确定调度方案,具体包括步骤如下:
步骤1:定义计数变量j=1,和计数向量k(ω)=1,ω=1,2,...n,g(ω)=1,ω=1,2,...m;
步骤2:依据加工工件πj的当前加工工序k(πj),确定对应的机器编号
Figure FDA0003718382280000036
步骤3:比较加工工件πj的上一道工序k(πj)-1的完成时间t1和机器
Figure FDA0003718382280000037
的上一个加工任务
Figure FDA0003718382280000038
完成时间t2;若k(πj)=1,则t1=0;若
Figure FDA0003718382280000039
则t2=0;
步骤4:
Figure FDA00037183822800000310
步骤5:令k(πj)=k(πj)+1,
Figure FDA00037183822800000311
步骤6:令j=j+1,若j≤n×m,转到步骤2,否则,结束。
3.根据权利要求1所述的一种基于动态信息累计狮群最优柔性调度方法,其特征在于,步骤B的具体实现过程包括:种群个体通过步骤A映射为一种调度方案,再根据式(Ⅰ)得到最大完成时间,也就是该种群个体的适应度值,将求取的该种群个体的适应度值即第一代狮群个体的适应度值初始化为每个种群个体最优,以所有种群个体最小的适应度值确定为全局最优。
4.根据权利要求1所述的一种基于动态信息累计狮群最优柔性调度方法,其特征在于,步骤C的具体实现过程包括:
1)计算新种群个体的适应度值,更新个体最优值和全局最优值;是指:
Figure FDA0003718382280000041
代表式(Ⅶ)更新后的狮子个体,
Figure FDA0003718382280000042
是第i个狮子第k次迭代的历史最优位置,若:
Figure FDA0003718382280000043
则令
Figure FDA0003718382280000044
此为新的个体最优;若存在:
Figure FDA0003718382280000045
pg k+1为新的全局最优;
2)根据新的个体最优值和全局最优值,更新个体信息累计因子和全局信息累计因子,是指:新的个体最优值pi和新的全局最优值pg,带入式(V)、式(Ⅵ),计算新的Cp和Cg
3)判断是否达到迭代停止条件,若是,进行步骤4);否则,返回步骤1);迭代次数达到预设最大值T;
4)输出优化后的调度方案和目标函数值即加工完成用时;是指:输出最后一次迭代得到的pg,通过半活动解码算法映射为一种调度方案,即优化后的调度方案,该调度方案的加工完成时间就是目标函数值。
5.一种基于动态信息累计狮群最优柔性调度系统,其特征在于,用于实现权利要求1-4任一所述基于动态信息累计狮群最优柔性调度方法,包括数学优化模型构建单元及数学优化模型优化单元;所述数学优化模型构建单元用于实现所述步骤(1),所述数学优化模型优化单元用于实现所述步骤(2)。
CN202110806542.8A 2021-07-16 2021-07-16 一种基于动态信息累计狮群最优柔性调度方法及系统 Active CN113408951B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110806542.8A CN113408951B (zh) 2021-07-16 2021-07-16 一种基于动态信息累计狮群最优柔性调度方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110806542.8A CN113408951B (zh) 2021-07-16 2021-07-16 一种基于动态信息累计狮群最优柔性调度方法及系统

Publications (2)

Publication Number Publication Date
CN113408951A CN113408951A (zh) 2021-09-17
CN113408951B true CN113408951B (zh) 2022-08-12

Family

ID=77686628

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110806542.8A Active CN113408951B (zh) 2021-07-16 2021-07-16 一种基于动态信息累计狮群最优柔性调度方法及系统

Country Status (1)

Country Link
CN (1) CN113408951B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114565245A (zh) * 2022-02-16 2022-05-31 石河子大学 一种计及电冷热的综合能源系统优化调度方法及装置
CN115564146B (zh) * 2022-11-18 2023-04-07 吉林大学 一种以改进的粒子群算法为内核的分组车间作业调度方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107817772A (zh) * 2017-10-17 2018-03-20 西南交通大学 一种柔性作业车间调度优化方法
CN109784603A (zh) * 2018-11-15 2019-05-21 长安大学 一种基于混合鲸鱼群算法求解柔性作业车间调度的方法
CN112668864A (zh) * 2020-12-24 2021-04-16 山东大学 一种基于狮群算法的车间生产排产方法及系统
CN112836949A (zh) * 2021-01-26 2021-05-25 西安理工大学 一种基于狮群算法的水库群双目标优化调度方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113011612B (zh) * 2021-03-23 2022-09-20 合肥工业大学 基于改进灰狼算法的生产与维修调度方法及系统
CN112862380B (zh) * 2021-04-12 2022-03-18 山东大学 一种基于混合算法的项目型产品装配车间人员调度方法、设备及存储介质

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107817772A (zh) * 2017-10-17 2018-03-20 西南交通大学 一种柔性作业车间调度优化方法
CN109784603A (zh) * 2018-11-15 2019-05-21 长安大学 一种基于混合鲸鱼群算法求解柔性作业车间调度的方法
CN112668864A (zh) * 2020-12-24 2021-04-16 山东大学 一种基于狮群算法的车间生产排产方法及系统
CN112836949A (zh) * 2021-01-26 2021-05-25 西安理工大学 一种基于狮群算法的水库群双目标优化调度方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"Joint Task Planning of UAV Groups using Improved Multi-objective Lion Swarm Optimization";Shengwei Dong;《2020 39th Chinese Control Conference》;20200729;全文 *
"Power System Optimization Based on Multi-agent Structure and Lion Swarm Optimization";Ying Guo 等;《Journal of Physics: Conference Series》;20201217;全文 *
"Prediction of equipment performance index based on improved chaotic lion swarm optimization–LSTM";Zhe Yang等;《Soft Computing》;20191101;全文 *
"基于改进NSGA-II算法的水库双目标优化调度";刘东 等;《西安理工大学学报》;20200515;全文 *

Also Published As

Publication number Publication date
CN113408951A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
CN113408951B (zh) 一种基于动态信息累计狮群最优柔性调度方法及系统
CN106959675B (zh) 一种面向流水车间的多目标调度优化方法
CN108389001B (zh) 一种基于step-nc的智能非线性工艺规划方法
CN109946965B (zh) 一种基于改进的多目标Jaya算法的离散制造车间排产方法
CN106875094A (zh) 一种基于多色集合遗传算法的多目标车间调度方法
CN108460463B (zh) 基于改进遗传算法的高端装备流水线生产调度方法
CN110618668A (zh) 柔性生产的绿色动态调度方法
CN111382915B (zh) 一种共融agv的柔性作业车间调度方法
CN110543151A (zh) 基于改进nsga-ⅱ求解车间节能调度问题的方法
CN112907150B (zh) 一种基于遗传算法的生产排程方法
CN109669423B (zh) 基于改进多目标灰狼算法获得零件加工最优调度方案的方法
CN112163808B (zh) 一种基于对立学习的自适应鲸鱼算法求解物流中心选址问题的方法
CN115130749A (zh) 融合nsga-ⅲ和topsis的数据驱动多目标优化方法
CN109917754B (zh) 一种基于多种群分布估计算法的机器人装配单元多目标布局优化方法
CN105550825A (zh) 云计算环境下的基于MapReduce并行化的柔性工厂作业调度方法
Wang et al. Genetic algorithm in flexible work shop scheduling based on multi-objective optimization
CN113281993B (zh) 一种贪心k-均值自组织神经网络多机器人路径规划方法
CN113326970A (zh) 一种混流装配线排序优化方法
CN113112121A (zh) 一种基于多目标非支配排序的车间布局调度优化方法
CN116880424A (zh) 一种基于多目标优化的多机器人调度方法以及装置
CN103473614A (zh) 基于碳排放评估模型的低碳工艺规划方法
CN116680958A (zh) 一种面向风机叶片双边协作的机器人喷涂轨迹优化方法
CN112487723B (zh) 一种基于多参数的中水热源泵系统优化设计方法
CN114676987A (zh) 一种基于超启发式算法的智能柔性作业车间主动调度方法
CN113935610A (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