CN111860754B - 基于蚁群以及遗传算法的agv调度方法 - Google Patents

基于蚁群以及遗传算法的agv调度方法 Download PDF

Info

Publication number
CN111860754B
CN111860754B CN202010678166.4A CN202010678166A CN111860754B CN 111860754 B CN111860754 B CN 111860754B CN 202010678166 A CN202010678166 A CN 202010678166A CN 111860754 B CN111860754 B CN 111860754B
Authority
CN
China
Prior art keywords
node
ant
vehicle
pheromone
path
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
CN202010678166.4A
Other languages
English (en)
Other versions
CN111860754A (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.)
Wuxi Hongyi Intelligent Technology Co ltd
Original Assignee
Wuxi Hongyi Intelligent Technology Co ltd
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 Wuxi Hongyi Intelligent Technology Co ltd filed Critical Wuxi Hongyi Intelligent Technology Co ltd
Priority to CN202010678166.4A priority Critical patent/CN111860754B/zh
Publication of CN111860754A publication Critical patent/CN111860754A/zh
Application granted granted Critical
Publication of CN111860754B publication Critical patent/CN111860754B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • 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/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • 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
    • G06Q10/00Administration; Management
    • G06Q10/08Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/087Inventory or stock management, e.g. order filling, procurement or balancing against orders

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Marketing (AREA)
  • Biomedical Technology (AREA)
  • Mathematical Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Educational Administration (AREA)
  • Computational Linguistics (AREA)
  • Game Theory and Decision Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Physiology (AREA)
  • Genetics & Genomics (AREA)
  • Accounting & Taxation (AREA)
  • Finance (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种AGV调度方法,尤其是一种基于蚁群以及遗传算法的AGV调度方法。在蚁群算法中引入车辆容量因素和时间窗因素改进蚂蚁状态转移概率,并且改进信息素挥发因子,使其能够随着计算进程自动调整,同时改进信息素更新策略,奖励超过全局最优解的精英蚂蚁。最后使用遗传算法中的选择、交叉、变异算子对蚁群算法得到的较优解进行局部优化,达到加快算法收敛速度,提高解的质量,能够明显解决传统优化算法进行路径规划时收敛速度慢、易陷入局部最优等缺点,能提高对实际问题的求解效率,减少迭代过程的盲目性。

Description

基于蚁群以及遗传算法的AGV调度方法
技术领域
本发明涉及一种AGV调度方法,尤其是一种基于蚁群以及遗传算法的AGV调度方法。
背景技术
随着现代工业和信息服务业的发展,人力资本愈发珍贵,人们越来越认识到仓储、货运等物流环节对提高产品利润的重要性。据资料显示,物流、运输等环节占整个制造业企业成本的50%以上,因此提高存储、运输的效率,充分发挥场地的性能优势,成为各大企业竞相发展的焦点。
自动化立体仓库是在现代物流系统中的一个重要组成部分,自动化立体仓库的智能化程度对整个物流业的发展有着重要的影响,智能化程度越高,则物流业越发达,合理的调度方案可以显著提高物流效率,从而节约物流成本,提升企业的效益。
蚁群算法(ant colony optimization,ACO)是一种启发式生物智能算法,以蚁群为研究对象,研究其在觅食过程中的群体性行为,并将这种群体智能行为抽象成数学算法,并应用到实际问题的求解过程中。遗传算法是一种模仿生物进化过程中,染色体的复制、交叉、变异等一系列变化过程,并将这种染色体的行为抽象成数学模型,应用到对实际问题的求解过程中。但是蚁群算法和遗传算法在解决AGV(Automated Guided Vehicle)调度的实际问题时容易出现无法收敛和陷入局部最优等不足,难以实现有效地调度。
发明内容
本发明的目的是克服现有技术中存在的不足,提供一种基于蚁群以及遗传算法的AGV调度方法,其能有效解决自动化仓储中多AGV协同的工作,提高AGV调度的效率,安全可靠。
按照本发明提供的技术方案,一种基于蚁群以及遗传算法的AGV调度方法,所述AGV调度方法包括如下步骤:
步骤1、对配送中心以及所述配送中心周围的n个客户节点,得到集合N={0,1,2,...,n},其中,节点0为配送中心,客户节点的集合为Nc,Nc={1,2,...,n},n为客户节点的数量;
对m个配送的车辆,得到车辆集合K,其中,K={1,2,...,m},在配送时,所有车辆必须从节点0开始并回到节点0;
对客户节点集合Nc中的任一节点i,存在所述节点i的时间窗[ai,bi],ai表示节点i最早开始接受货物的时间,bi为节点i最晚接受货物的时间,qi为节点i的需求量,车辆在节点i的服务完成时间为si;对车辆集合K中的任一车辆k,车辆k达到节点i的时刻为Tik,车辆k从节点i达到节点j的时刻为Tij
步骤2、根据步骤1中配送中心、客户节点集合Nc、车辆集合K,建立带时间窗的车辆路径的数学模型,具体为
Figure BDA0002584733290000021
Figure BDA0002584733290000022
Figure BDA0002584733290000023
Figure BDA0002584733290000024
Figure BDA0002584733290000025
Figure BDA0002584733290000026
Figure BDA0002584733290000027
Figure BDA00025847332900000213
Figure BDA0002584733290000028
其中,
Figure BDA0002584733290000029
Figure BDA00025847332900000210
步骤3、确定蚂蚁的数量R,为每个蚂蚁设置对应的禁忌表Tabu以及节点存储表Tau;初始化迭代参数NF、初始化信息素集Qr、初始化禁忌表Tabu以及初始化节点存储表Tau,并给定最大迭代次数NFmax、信息素重要度因子α、期望程度因子β、车辆容量因素的权衡因子θ、遗传算法中的交叉概率Pc以及变异概率Pm
步骤4、将R只蚂蚁初始化到不同的客户节点上,并将每个蚂蚁对应的客户节点写入每个蚂蚁的禁忌表Tabu表内;
步骤5、对于节点i上的蚂蚁r,将节点i加入蚂蚁r的节点存储表Tau中;计算得到蚂蚁r的转移概率,按轮盘赌的方式选择下一个客户节点,具体地,
Figure BDA00025847332900000211
其中,
Figure BDA00025847332900000212
为蚂蚁r在t时刻从节点i转移至节点j的概率,τij(t)为t时刻节点i与节点j之间的信息素强度;ηij(t)为t时刻节点i到节点j的期望启发函数,allow为蚂蚁r下一步允许访问的节点集合;ti为目前路径到达点i时经过的时间,(aj,bj)为节点j的时间窗,(az,bz)为节点z的时间窗,w1+w2=1;
步骤6、查询蚂蚁r的节点存储表Tau,当蚂蚁r遍历完成所有的节点后,转到步骤7,否则,根据步骤5将蚂蚁r成功转至节点j,将节点j加入禁忌表Tabu表;
在将蚂蚁r转移至节点j后,根据蚂蚁r的路径查询车辆的载重以及节点j的时间窗约束,当车辆在的载重大于车辆的额定载重Qe或节点j的时间不满足时间窗约束时,则将蚂蚁r返回到节点0,同时将节点0加入到节点存储表Tau表中,并转到步骤5;
步骤7:判断蚂蚁r是否回到节点0,若没有回到节点0,则将节点i加入到禁忌表Tabu中,并转至步骤8;否则,直接转至步骤8;
步骤8:判断r是否大于蚂蚁数R,若r>R,转到步骤9;否则,令r=r+1并转至步骤5;
步骤9:根据每个蚂蚁的节点存储表Tau,计算每个蚂蚁所走的路径长度,并将R个蚂蚁的路径按距离排列在集合Solution中,找出所有蚂蚁中行走距离最短的那只蚂蚁,将行走距离最短蚂蚁所行走的路径作为当前迭代过程中的最优解now_bestcost,若当前最优解now_bestcost小于全局最优解best_cost,将全局最优解best_cost更新为当前最优解now_bestcost;
步骤10、将当前迭代过程中的最优解now_bestcost复制W个,作为遗传算法中的部分初始种群,通过轮盘赌的方法每次从步骤9中所得的所有剩下的蚂蚁路径中选择两个行走路径进行比较,将路径长度较短的路线保留到初始种群中,直至种群数达到R个,将数量为R的初始种群中的所有解中全部去除配送中心0,得到R条客户节点的排列,每条这样的排列被称作染色体,以得到遗传算法的父代染色体;
步骤11、产生一个随机数rand1,0<rand1<1,当随机数rand1大于设定的交叉概率pc时,执行交叉步骤,否则,随机复制一个父代染色体作为子代染色体,重复执行上述操作,直至得到R个子代染色体;
步骤12、产生一个随机数rand2,0<rand2<1,当随机数rand2小于变异概率pm,执行变异操作,否则,随机复制一个父代染色体作为子代染色体,重复执行变异操作,直至再次得到R个子代染色体;
步骤13、对R个子代染色体进行解码,并计算各子代个体染色体解码后的路径长度;
步骤14:更新集合Solution中路径较长的解,更新节点存储表Tau,更新当前最优解now_bestcost、全局最优解best_cost以及信息素集Qr
步骤15、清空每只蚂蚁的禁忌表Tabu和节点存储表Tau,清空集合Solution,清空now_bestcost以及清空allow集合;
步骤16、判断迭代参数NF是否超过设置的最大迭代次数NFmax,若迭代参数NF大于最大迭代次数NCmax,则输出最优路径和最优路径长度;否则,NF=NF+1,跳到步骤4。
步骤14中,若当前最优解now_bestcost小于全局最优解best_cost,更新信息素集Qr为:
Figure BDA0002584733290000041
Figure BDA0002584733290000042
Figure BDA0002584733290000043
Figure BDA0002584733290000044
其中:τij(t+1)为更新后节点i和节点j之间的信息素强度;τij(t)为更新前的节点i和节点j之间的信息素强度;delta为惩罚值,能够奖励超过全局最优解的蚂蚁;ρ为信息素挥发系数,Δτij(t)为路径节点i到节点j上的信息素增加量,即所有蚂蚁在该段路径上释放的信息素之和,
Figure BDA0002584733290000046
为蚂蚁r在客户节点i到客户节点j之间所释放的信息素;R为蚁群的数量;Qr为信息素释放总量;Lr为蚂蚁r所走过路径(从节点i到节点j)的总长度。
信息素挥发系数ρ的更新为:
Figure BDA0002584733290000045
其中,:ρmin为信息素挥发系数ρ的最小值,x,y为系数。
本发明的优点:在蚁群算法中引入车辆容量因素和时间窗因素改进蚂蚁状态转移概率,并且改进信息素挥发因子,使其能够随着计算进程自动调整,同时改进信息素更新策略,奖励超过全局最优解的精英蚂蚁。最后使用遗传算法中的选择、交叉、变异算子对蚁群算法得到的较优解进行局部优化,达到加快算法收敛速度,提高解的质量,能够明显解决传统优化算法进行路径规划时收敛速度慢、易陷入局部最优等缺点,能提高对实际问题的求解效率,减少迭代过程的盲目性。
附图说明
图1为本发明进行蚁群算法以及遗传算法的流程图。
图2为本发明交叉得到子代染色体的示意图。
具体实施方式
下面结合具体附图和实施例对本发明作进一步说明。
如图1所示:为了能有效解决自动化仓储中多AGV协同的工作,提高AGV调度的效率,本发明的AGV调度方法包括如下步骤:
步骤1、对配送中心以及所述配送中心周围的n个客户节点,得到集合N={0,1,2,...,n},其中,节点0为配送中心,客户节点的集合为Nc,Nc={1,2,...,n},n为客户节点的数量;
对m个配送的车辆,得到车辆集合K,其中,K={1,2,...,m},在配送时,所有车辆必须从节点0开始并回到节点0;
对客户节点集合Nc中的任一节点i,存在所述节点i的时间窗[ai,bi],ai表示节点i最早开始接受货物的时间,bi为节点i最晚接受货物的时间,qi为节点i的需求量,车辆在节点i的服务完成时间为si;对车辆集合K中的任一车辆k,车辆k达到节点i的时刻为Tik,车辆k从节点i达到节点j的时刻为Tij
具体地,本发明所要解决的问题是带时间窗的车辆路径问题,车辆路径问题是一类经典的组合优化问题,也是物流领域研究的核心问题之一。如果在车辆路径问题中增加了对客户偏好服务时间的约束,即个性化地限定了客户被服务时间的范围,即延伸出了带时间窗的车辆路径问题。带时间窗的车辆路径问题可描述为:假设一个配送中心为周围若干个位于不同地理位置、且对货物送达时间有不同要求的客户节点提供配送服务。其中,配送中心用于运送的车辆都是同一型号的(即拥有相同的容量、速度);客户节点对车辆出入的时间有限制;车辆在所有客户节点有相同的停留服务时间;车辆需要从配送中心出发,经过所有的客户节点后回到配送中心。在这样的前提下,考虑如何合理安排车辆配送路线,既能在满足各项约束的条件下,将货物安全送到客户的手中,又能实现配送车辆行驶总路径最短。
因此,每个节点都包含三个信息,该节点的坐标、该节点的需求量qi以及该节点的时间窗。
步骤2、根据步骤1中配送中心、客户节点集合Nc、车辆集合K,建立带时间窗的车辆路径的数学模型,具体为
Figure BDA0002584733290000051
Figure BDA0002584733290000052
Figure BDA0002584733290000053
Figure BDA0002584733290000054
Figure BDA0002584733290000055
Figure BDA0002584733290000056
Figure BDA0002584733290000057
Figure BDA0002584733290000058
Figure BDA0002584733290000061
其中,
Figure BDA0002584733290000062
Figure BDA0002584733290000063
具体地,dijk为车辆k从节点i到节点j之间的欧式距离,公式1为车辆总行驶路径最小目标函数;公式2为配送车辆的载重约束;公式3表示每个客户节点只能由一辆车进行服务;公式4表示车辆服务某节点后必须从此节点出发去下一节点约束;公式5和公式6表示每个客户节点只允许被服务一次;公式7为消去子回路约束;公式8和公式9为时间窗约束。
步骤3、确定蚂蚁的数量R,为每个蚂蚁设置对应的禁忌表Tabu以及节点存储表Tau;初始化迭代参数NF、初始化信息素集Qr、初始化禁忌表Tabu以及初始化节点存储表Tau,并给定最大迭代次数NFmax、信息素重要度因子α、期望程度因子β、车辆容量因素的权衡因子θ、遗传算法中的交叉概率Pc以及变异概率Pm
具体地,初始迭代参数NF一般取1,最大迭代次数NFmax可取100,α是信息素在概率计算中的权重,表示信息素浓度对蚂蚁进行路径选择的影响程度。β为启发信息重要程度,表示路径上的期望度对蚂蚁进行路径选择的影响程度,θ是车辆容量因素的权衡因子,θ值越大,车辆选择货物需求量大的客户节点的可能性越大,车辆最终满载率可能最高。蚂蚁的数量R和客户节点数量可一致为100个。信息素集Qr为每两个客户节点间的信息素浓度集合,初始化信息素集Qr即把每两个节点的初始信息素浓度都设置为客户节点数(100)/所有节点(包括配送中心)间的总距离,总距离一般是指节点0到其他所有节点的总距离+节点1到其他节点的总距离+…+节点n到其他所有节点的总距离。每只蚂蚁都有一个禁忌表Tabu和一个节点存储表Tau,禁忌表Tabu存储蚂蚁经过的客户节点(不包含配送中心0),节点存储表Tau存储蚂蚁行走的路径中所有的节点(包括但并不止一个配送中心0),初始化禁忌表Tabu、初始化节点存储表Tau后,则禁忌表Tabu为空、节点存储表Tau内为0。此外,还包括集合Solution,所述集合Solution中储存的是算法运行后所得到的带时间窗的车辆路径规划问题的解。
遗传算法中的交叉概率Pc、变异概率Pm的具体情况可以根据需要进行选择设定,具体为本技术领域人员所熟知,此处不再赘述。
步骤4、将R只蚂蚁初始化到不同的客户节点上,并将每个蚂蚁对应的客户节点写入每个蚂蚁的禁忌表Tabu表内;
具体地,将R只蚂蚁随机化分布到客户节点上,各蚂蚁初始的化客户节点编号写入禁忌表Tabu表内。
步骤5、对于节点i上的蚂蚁r,将节点i加入蚂蚁r的节点存储表Tau中;计算得到蚂蚁r的转移概率,按轮盘赌的方式选择下一个客户节点,具体地,
Figure BDA0002584733290000071
其中,
Figure BDA0002584733290000072
为蚂蚁r在t时刻从节点i转移至节点j的概率,τij(t)为t时刻节点i与节点j之间的信息素强度;ηij(t)为t时刻节点i到节点j的期望启发函数,,一般地,ηij(t)为节点i到节点j之间距离的倒数,allow为蚂蚁r下一步允许访问的节点集合;ti为目前路径到达点i时经过的时间,(aj,bj)为节点j的时间窗,(az,bz)为节点z的时间窗,w1+w2=1;μij(t)为车辆容量因数,设车辆在节点i的累积载重量为Gi,qi为节点i的需求量,Qe为车辆的额定载重,则车辆容量因数μij(t)为/>
Figure BDA0002584733290000073
具体地,集合allow为蚂蚁r下一步允许访问的客户节点集合,即为NC集合中排除禁忌表Tabu中客户节点后所剩的允许访问的客户节点;从集合allow内依次选取客户节点,逐次累加节点i与各个客户节点间的转移概率,将这个累加值与(0,1)中的随机数比较,若累加值大于该随机数,则选择所累加的最后一个客户节点j为蚂蚁r行走的下一个节点(比如节点2到节点3的概率小于该随机数,然后就用这个概率加上从节点2到节点4的概率,累加后的概率和随机数比较,如果大于随机数,选择4节点作为下一个行走的节点),。一般地,w1=0.8,w2=0.2。
步骤6、查询蚂蚁r的节点存储表Tau,当蚂蚁r遍历完成所有的节点后,转到步骤7,否则,根据步骤5将蚂蚁r成功转至节点j,将节点j加入禁忌表Tabu表;
在将蚂蚁r转移至节点j后,根据蚂蚁r的路径查询车辆的载重以及节点j的时间窗约束,当车辆在的载重大于车辆的额定载重Qe或节点j的时间不满足时间窗约束时,则将蚂蚁r返回到节点0,同时将节点0加入到节点存储表Tau表中,并转到步骤5;
具体地,蚂蚁r遍历完成所有节点就是蚂蚁r遍历所有客户节点1次,也就是必须满足步骤2中的约束公式5和公式6,也相当于禁忌表Tabu中包含全部客户节点。所述车辆满足载重约束就是车辆目前所经过的客户节点的需求量qi的累加值不能超过车辆的额定载重Qe,也就是步骤2要求的约束公式2,所述车辆满足时间窗约束就是车辆到达客户节点的时间必须在客户节点的时间窗的范围内,即在客户节点的服务时间内,也就是约束公式8和公式9。
下面进行举例说明,具体地:如果有节点1,节点2,节点3,节点4,节点5,节点6,节点7,节点8,节点9,节点10这10个客户节点,那么解的形式就可能是0,1,2,3,0,4,5,6,0,7,8,9,0,10,0。这代表分4辆小车输送货物,第一辆小车路径0-1-2-3-0,第二辆小车路径0-4-5-6-0,第三辆小车路径0-7-8-9-0,第四辆小车路径0-10-0。判断是否满足载重,就是判断每辆小车是否超过其额定载重Qe,如第一辆载重就是节点1货物需求量q1+节点2货物需求量q2+节点3的货物需求量q3,判断其是否超过小车的额定载重Qe,小车的额定载重是固定的。
判断是否满足时间窗约束,以第一辆小车为例,在其走到节点1时,判断到达时间是否在节点1的服务时间窗的范围内。
步骤7:判断蚂蚁r是否回到节点0,若没有回到节点0,则将节点i加入到禁忌表Tabu中,并转至步骤8;否则,直接转至步骤8;
具体地,判断蚂蚁r是否回到节点0,即为判断节点存储表Tau中最后一个节点是否为节点0。
步骤8:判断r是否大于蚂蚁数R,若r>R,转到步骤9;否则,令r=r+1并转至步骤5;
具体地,r小于蚂蚁数R时,则表明存在蚂蚁为进行上述遍历。
步骤9:根据每个蚂蚁的节点存储表Tau,计算每个蚂蚁所走的路径长度,并将R个蚂蚁的路径按距离排列在集合Solution中,找出所有蚂蚁中行走距离最短的那只蚂蚁,将行走距离最短蚂蚁所行走的路径作为当前迭代过程中的最优解now_bestcost,若当前最优解now_bestcost小于全局最优解best_cost,将全局最优解best_cost更新为当前最优解now_bestcost;
具体地,经过上述步骤后,此时所有蚂蚁都已经周游所有客户节点。每只蚂蚁的节点存储表Tau即为每只蚂蚁所行走的所有节点的集合,计算各蚂蚁所走路径的长度,就是依次计算节点存储表Tau中的节点间的欧式距离并累加,然后按路径长短排列在集合Solution中。找出所有蚂蚁中行走距离最短的那只蚂蚁,也就是约束公式1,将其行走的路径作为当前迭代过程中的最优解now_bestcost,如果当前最优解now_bestcost小于全局最优解best_cost,将全局最优解best_cost更新为当前迭代过程中的当前最优解now_bestcost。一般地,全局最优解best_cost在初始状态下为无穷大。
接下来是利用遗传算法优化蚁群算法所求的路径集,遗传算法主要分为选择、交叉和遗传三个操作,即分别为步骤10、步骤11和步骤12。
步骤10、将当前迭代过程中的最优解now_bestcost复制W个,作为遗传算法中的部分初始种群,通过轮盘赌的方法每次从步骤9中所得的所有剩下的蚂蚁路径中选择两个行走路径进行比较,将路径长度较短的路线保留到初始种群中,直至种群数达到R个,将数量为R的初始种群中的所有解中全部去除配送中心0,得到R条客户节点的排列,每条这样的排列被称作染色体,以得到遗传算法的父代染色体;
具体地,遗传算法中选择就是按照设定的原则或方法,即选择种群中优秀的个体作为父代来繁衍子代的过程。将当前迭代过程中的最优解now_bestcost复制W个,W可以为R/5,或是其他所需不大于R的数量,具体可以根据需要进行选择,作为遗传算法中的部分初始种群(这里的初始种群就是将要被应用于遗传算法的带时间窗的车辆路径问题的若干解);然后通过轮盘赌的方法每次从步骤9中所得的所有剩下的蚂蚁路径中选择两个个体进行比较,将路径长度较短的路线保留到初始种群中,直至种群数达到R个。将数量为R的初始种群中的所有解中全部去除配送中心0,得到R条客户节点的排列,每条这样的排列被称作染色体,初始染色体的集合称之为遗传算法的父代染色体。
具体实施时,初始种群中的每一个解都相当于某一个蚂蚁的行驶路径,如节点0,节点1,节点2,节点3,节点0,节点4,节点5,节点6,节点0,节点7,节点8,节点9,节点0,去掉上述解中的配送中心编号,只保留客户节点的编号,并且客户节点编号的序列保持不变即为可行解对应的染色体结构,则上述解转换而来的染色体结构为:123456789,这样就可以在保证同一限定条件下,编码能够保证所有染色体的长度都是一样的,方便后续的交叉变异操作。
步骤11、产生一个随机数rand1,0<rand1<1,当随机数rand1大于设定的交叉概率pc时,执行交叉步骤,否则,随机复制一个父代染色体作为子代染色体,重复执行上述操作,直至得到R个子代染色体;
具体地,所述交叉操作即是从步骤10所得的父代染色体中随机选择两个染色体,首先在这两个染色体中随机选择两个点i,j,其中0≤i≤j≤n,即染色体长度为n。然后将亲代(从父代染色体中选择的两个染色体就是亲代1和亲代2)P1从i到j的基因填入子代相同的位置(染色体中的一个值代表一个基因)。接着将亲代P2的基因不重复地依次填入子代中;交换两个亲代染色体就可得到另一个子代染色体,如图2所示。
步骤12、产生一个随机数rand2,0<rand2<1,当随机数rand2小于变异概率pm,执行变异操作,否则,随机复制一个父代染色体作为子代染色体,重复执行变异操作,直至再次得到R个子代染色体;
具体地,变异操作时,首先随机选择两个父代染色体,在这两个父代染色体中随机确定两个基因位,然后把这两个基因位之间的染色体基因进行反序排列,即得到子代个体的染色体结构。
本发明实施例中,在遗传算法中,通过增加变异算子来模拟自然界中基因突变这一环节,变异算子的作用虽然小于交叉算子,但也是遗传算法所必须的,它能够有效地增强遗传算法的局部搜索能力,避免算法优化得到的解陷入局部最优,发生过早收敛现象。遗传算法中常用的变异算子有反转变异、2-opt变异等多种形式。),本发明实施例中,选用的变异算子是反转变异方式。
步骤13、对R个子代染色体进行解码,并计算各子代个体染色体解码后的路径长度;
具体地,解码时按照步骤2中的约束条件,依次将同时符合车辆载重约束和客户点时间窗约束的客户分配给配送车辆。以下面的染色体结构为例:(节点5,节点10,节点6,节点9,节点3,节点8,节点7,节点2,节点4,节点1)结合步骤2型中的约束,对其进行染色体解码:首先将客户5分配给车辆1,更新车辆1的时间和载重,再将其后的节点10分配给车辆1,检验载重约束和时间窗约束,如果载重和时间窗约束中有一个及以上不满足,则形成子路径(0,5,0),即车辆1只为客户点5配送货物;若都符合要求,则将客户10也车辆1配送,再依次检验后面的客户,只要车辆超载或是不满足时间窗约束,则启动下一辆车,依次检验下去,直到染色体结构中所有的客户点都被分配完。上述子代染色体解码出来的可行解为:(0,5,10,6,0,9,3,8,0,7,2,4,1,0)即共需要三辆车完成配送,其路径分别为(0,5,10,6,0),(0,9,3,8,0)和(0,7,2,4,10)。
步骤14:更新集合Solution中路径较长的解,更新节点存储表Tau,更新当前最优解now_bestcost、全局最优解best_cost以及信息素集Qr
具体地,运行到此步骤时本次迭代过程中的蚁群算法和遗传算法都已运行完毕,比较蚁群算法中每只蚂蚁行走的路径长度和各子代个体染色体解码后的路径长度,用子代个体染色体解码后路径长度较短的解替代集合Solution中蚂蚁行走的路径长度较长的解,同时替换相应的节点存储表Tau。如果遗传算法所得的最短路径长度小于当前迭代过程中的最优解now_bestcost,更新算法运行过程中的最优解best_cost;按照下式更新各节点间的信息素强度,即为更新信息素集Qr
Figure BDA0002584733290000101
上述公式的含义是:更新后客户节点i与客户节点j之间的信息素强度=原有信息素挥发后所剩的信息素强度+每只蚂蚁在这段距离内新产生的信息素强度+每次派出的最强蚂蚁的总路程的信息素强度更新。其中
Figure BDA0002584733290000102
Figure BDA0002584733290000103
Figure BDA0002584733290000104
其中:τij(t+1)为更新后节点i和j之间的信息素强度;τij(t)为更新前节点i和节点j之间的信息素强度;delta为惩罚值,能够奖励超过全局最优解的蚂蚁;ρ为信息素挥发系数;Δτij(t)为路径节点i到节点j上的信息素增加量,即所有蚂蚁在该段路径上释放的信息素之和,
Figure BDA0002584733290000112
为蚂蚁r在客户节点i到客户节点j之间所释放的信息素;R为蚁群的数量;Qm为信息素释放总量;Lr为蚂蚁r所走过路径(从节点i到节点j)的总长度。具体实施时,信息素释放总量Qm、蚂蚁r所走的总长度Lr均可以计算得到,从而可以计算得到每一只蚂蚁的
Figure BDA0002584733290000113
Δτij(t)是所有蚂蚁的/>
Figure BDA0002584733290000114
相加。
在经典蚁群算法中,信息素挥发系数是一个恒定不变的值,信息素挥发系数跟蚁群算法的全局搜索能力和算法的收敛速度有着紧密的联系;由于信息素挥发系数ρ的缘由,从来没有被搜寻的路径上的信息素会慢慢减为0,当信息素挥发系数ρ较大时,全局搜寻能力和随机搜索性能也会受到影响,容易陷入局部最优值;通过减小信息素挥发系数ρ能够增强算法的全局搜索能力和随机搜索性能,但同时会使得算法的收敛速度降低。
针对固定信息素挥发系数ρ值不变产生的问题,本发明信息素挥发系数ρ能够根据算法的进程自动地调整信息素挥发系数,意在通过自适应调整信息素挥发系数ρ来解决算法存在的全局性较差问题。由分析结果可知,ρ取值在0.3~0.8时效果比较好,设信息素挥发系数ρ在算法刚开始执行的时候初值为ρ=0.8,随着迭代的进行ρ按照下式调节:
Figure BDA0002584733290000111
其中,:ρmin为信息素挥发系数ρ的最小值,x,y为系数。一般地,ρmin=0.2;y可取25,x可取=1。
步骤15、清空每只蚂蚁的禁忌表Tabu和节点存储表Tau,清空集合Solution,清空now_bestcost以及清空allow集合;
具体地,清空每只蚂蚁的禁忌表Tabu和节点存储表Tau,清空集合Solution,清空now_bestcost以及清空allow集合,从而能便于进行后续的迭代步骤。
步骤16、判断迭代参数NF是否超过设置的最大迭代次数NFmax,若迭代参数NF大于最大迭代次数NCmax,则输出最优路径和最优路径长度;否则,NF=NF+1,跳到步骤4。
具体地,由上述可知,利用NF与最大迭代次数NFmax的关系,能对整个方法的过程进行控制。在迭代参数NF大于最大迭代次数NCmax,则整个方法结束,并输出最优路径和最优路径长度,所述最优路径即为步骤14中的全局最优解best_cost以及对应的路径长度。
综上,本发明在蚁群算法中引入车辆容量因素和时间窗因素改进蚂蚁状态转移概率,并且改进信息素挥发因子,使其能够随着计算进程自动调整,同时改进信息素更新策略,奖励超过全局最优解的精英蚂蚁。最后使用遗传算法中的选择、交叉、变异算子对蚁群算法得到的较优解进行局部优化,达到加快算法收敛速度,提高解的质量,能够明显解决传统优化算法进行路径规划时收敛速度慢、易陷入局部最优等缺点,能提高对实际问题的求解效率,减少迭代过程的盲目性。

Claims (3)

1.一种基于蚁群以及遗传算法的AGV调度方法,其特征是,所述AGV调度方法包括如下步骤:
步骤1、对配送中心以及所述配送中心周围的n个客户节点,得到集合N={0,1,2,...,n},其中,节点0为配送中心,客户节点的集合为Nc,Nc={1,2,...,n},n为客户节点的数量;
对m个配送的车辆,得到车辆集合K,其中,K={1,2,...,m},在配送时,所有车辆必须从节点0开始并回到节点0;
对客户节点集合Nc中的任一节点i,存在所述节点i的时间窗[ai,bi],ai表示节点i最早开始接受货物的时间,bi为节点i最晚接受货物的时间,qi为节点i的需求量,车辆在节点i的服务完成时间为si;对车辆集合K中的任一车辆k,车辆k达到节点i的时刻为Tik,车辆k从节点i达到节点j的时刻为Tij
步骤2、根据步骤1中配送中心、客户节点集合Nc、车辆集合K,建立带时间窗的车辆路径的数学模型,具体为
Figure FDA0004180784110000011
Figure FDA0004180784110000012
Figure FDA0004180784110000013
Figure FDA0004180784110000014
Figure FDA0004180784110000015
Figure FDA0004180784110000016
Figure FDA0004180784110000017
Figure FDA0004180784110000018
Figure FDA0004180784110000019
其中,
Figure FDA00041807841100000110
Figure FDA00041807841100000111
dijk为车辆k从节点i到节点j之间的欧式距离;
步骤3、确定蚂蚁的数量R,为每个蚂蚁设置对应的禁忌表Tabu以及节点存储表Tau;初始化迭代参数NF、初始化信息素集Qr、初始化禁忌表Tabu以及初始化节点存储表Tau,并给定最大迭代次数NFmax、信息素重要度因子α、期望程度因子β、车辆容量因素的权衡因子θ、遗传算法中的交叉概率Pc以及变异概率Pm
步骤4、将R只蚂蚁初始化到不同的客户节点上,并将每个蚂蚁对应的客户节点写入每个蚂蚁的禁忌表Tabu表内;
步骤5、对于节点i上的蚂蚁r,将节点i加入蚂蚁r的节点存储表Tau中;计算得到蚂蚁r的转移概率,按轮盘赌的方式选择下一个客户节点,具体地,
Figure FDA0004180784110000021
其中,
Figure FDA0004180784110000022
为蚂蚁r在t时刻从节点i转移至节点j的概率,τij(t)为t时刻节点i与节点j之间的信息素强度;ηij(t)为t时刻节点i到节点j的期望启发函数,allow为蚂蚁r下一步允许访问的节点集合;ti为目前路径到达点i时经过的时间,(aj,bj)为节点j的时间窗,(az,bz)为节点z的时间窗,w1+w2=1;μij(t)为车辆容量因数;
步骤6、查询蚂蚁r的节点存储表Tau,当蚂蚁r遍历完成所有的节点后,转到步骤7,否则,根据步骤5将蚂蚁r成功转至节点j,将节点j加入禁忌表Tabu表;
在将蚂蚁r转移至节点j后,根据蚂蚁r的路径查询车辆的载重以及节点j的时间窗约束,当车辆在的载重大于车辆的额定载重Qe或节点j的时间不满足时间窗约束时,则将蚂蚁r返回到节点0,同时将节点0加入到节点存储表Tau表中,并转到步骤5;
步骤7:判断蚂蚁r是否回到节点0,若没有回到节点0,则将节点i加入到禁忌表Tabu中,并转至步骤8;否则,直接转至步骤8;
步骤8:判断r是否大于蚂蚁数R,若r>R,转到步骤9;否则,令r=r+1并转至步骤5;
步骤9:根据每个蚂蚁的节点存储表Tau,计算每个蚂蚁所走的路径长度,并将R个蚂蚁的路径按距离排列在集合Solution中,找出所有蚂蚁中行走距离最短的那只蚂蚁,将行走距离最短蚂蚁所行走的路径作为当前迭代过程中的最优解now_bestcost,若当前最优解now_bestcost小于全局最优解best_cost,将全局最优解best_cost更新为当前最优解now_bestcost;
步骤10、将当前迭代过程中的最优解now_bestcost复制W个,作为遗传算法中的部分初始种群,通过轮盘赌的方法每次从步骤9中所得的所有剩下的蚂蚁路径中选择两个行走路径进行比较,将路径长度较短的路线保留到初始种群中,直至种群数达到R个,将数量为R的初始种群中的所有解中全部去除配送中心0,得到R条客户节点的排列,每条这样的排列被称作染色体,以得到遗传算法的父代染色体;
步骤11、产生一个随机数rand1,0<rand1<1,当随机数rand1大于设定的交叉概率pc时,执行交叉步骤,否则,随机复制一个父代染色体作为子代染色体,重复执行上述操作,直至得到R个子代染色体;
步骤12、产生一个随机数rand2,0<rand2<1,当随机数rand2小于变异概率pm,执行变异操作,否则,随机复制一个父代染色体作为子代染色体,重复执行变异操作,直至再次得到R个子代染色体;
步骤13、对R个子代染色体进行解码,并计算各子代个体染色体解码后的路径长度;
步骤14:更新集合Solution中路径较长的解,更新节点存储表Tau,更新当前最优解now_bestcost、全局最优解best_cost以及信息素集Qr
步骤15、清空每只蚂蚁的禁忌表Tabu和节点存储表Tau,清空集合Solution,清空now_bestcost以及清空allow集合;
步骤16、判断迭代参数NF是否超过设置的最大迭代次数NFmax,若迭代参数NF大于最大迭代次数NCmax,则输出最优路径和最优路径长度;否则,NF=NF+1,跳到步骤4。
2.根据权利要求1所述基于蚁群以及遗传算法的AGV调度方法,其特征是,步骤14中,若当前最优解now_bestcost小于全局最优解best_cost,更新信息素集Qr为:
Figure FDA0004180784110000031
Figure FDA0004180784110000032
Figure FDA0004180784110000033
Figure FDA0004180784110000034
其中:τij(t+1)为更新后节点i和节点j之间的信息素强度;τij(t)为更新前的节点i和节点j之间的信息素强度;delta为惩罚值,能够奖励超过全局最优解的蚂蚁;ρ为信息素挥发系数,Δτij(t)为路径节点i到节点j上的信息素增加量,即所有蚂蚁在该段路径上释放的信息素之和,
Figure FDA0004180784110000035
为蚂蚁r在客户节点i到客户节点j之间所释放的信息素;R为蚁群的数量;Qr为信息素释放总量;Lr为蚂蚁r所走过路径即从节点i到节点j的总长度。
3.根据权利要求2所述基于蚁群以及遗传算法的AGV调度方法,其特征是,信息素挥发系数ρ的更新为:
Figure FDA0004180784110000041
其中,ρmin为信息素挥发系数ρ的最小值,x,y为系数。
CN202010678166.4A 2020-07-15 2020-07-15 基于蚁群以及遗传算法的agv调度方法 Active CN111860754B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010678166.4A CN111860754B (zh) 2020-07-15 2020-07-15 基于蚁群以及遗传算法的agv调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010678166.4A CN111860754B (zh) 2020-07-15 2020-07-15 基于蚁群以及遗传算法的agv调度方法

Publications (2)

Publication Number Publication Date
CN111860754A CN111860754A (zh) 2020-10-30
CN111860754B true CN111860754B (zh) 2023-06-30

Family

ID=72983390

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010678166.4A Active CN111860754B (zh) 2020-07-15 2020-07-15 基于蚁群以及遗传算法的agv调度方法

Country Status (1)

Country Link
CN (1) CN111860754B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112327621B (zh) * 2020-11-02 2022-07-08 金航数码科技有限责任公司 一种基于蚁群算法的柔性生产线自适应控制系统及方法
CN112346839B (zh) * 2020-11-06 2023-06-20 东方红卫星移动通信有限公司 一种基于进化算法的关联任务调度方法
CN113009928A (zh) * 2020-12-31 2021-06-22 河北科技大学 基于改进蚁群算法的无人机冲突解脱方法
CN112926769A (zh) * 2021-02-08 2021-06-08 日日顺供应链科技股份有限公司 一种物流配送路径规划方法及装置
CN113009821B (zh) * 2021-02-10 2022-08-19 上海海事大学 一种基于信息素初始分配和动态更新的蚁群算法优化方法
CN113158580B (zh) * 2021-05-10 2021-12-10 南京林业大学 实木板材的一维排样方法
CN113222520A (zh) * 2021-06-16 2021-08-06 江苏佳利达国际物流股份有限公司 一种基于蚁群算法的货物优化配送方法及系统
CN113408823B (zh) * 2021-07-13 2022-12-13 重庆理工大学 一种城市突发事件紧急救援方法
CN113705879A (zh) * 2021-08-24 2021-11-26 武汉理工大学 一种多车场多车型车辆路径规划方法
CN113743608B (zh) * 2021-10-21 2024-01-30 安徽工程大学 一种基于基因库解决堆垛机路径规划的蚁群遗传混合算法
CN114722984B (zh) * 2022-04-20 2022-11-22 张家界富源仿真花有限公司 一种智能货物仓储优化方法
CN114780247B (zh) * 2022-05-17 2022-12-13 中国地质大学(北京) 一种流速和资源感知的流应用调度方法及系统
CN114971044A (zh) * 2022-06-02 2022-08-30 圆通速递有限公司 一种派件配送路径规划方法
CN115118724A (zh) * 2022-06-23 2022-09-27 福州大学 基于蚁群算法的多无人机辅助边缘计算系统部署优化方法
CN115617510A (zh) * 2022-09-28 2023-01-17 湖北公众信息产业有限责任公司 一种基于信息素遗传蚁群算法的服务器动态均衡方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104408589A (zh) * 2014-10-24 2015-03-11 陕西科技大学 基于混合粒子群算法的agv优化调度方法
CN105760954A (zh) * 2016-02-15 2016-07-13 南通大学 一种基于改进蚁群算法的泊车系统路径规划方法
CN106779252A (zh) * 2017-02-10 2017-05-31 泉州装备制造研究所 一种基于改进量子蚁群算法的agv实时路线规划方法
CN108460485A (zh) * 2018-03-05 2018-08-28 重庆邮电大学 一种基于改进蚁群算法的旅行商问题求解方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104408589A (zh) * 2014-10-24 2015-03-11 陕西科技大学 基于混合粒子群算法的agv优化调度方法
CN105760954A (zh) * 2016-02-15 2016-07-13 南通大学 一种基于改进蚁群算法的泊车系统路径规划方法
CN106779252A (zh) * 2017-02-10 2017-05-31 泉州装备制造研究所 一种基于改进量子蚁群算法的agv实时路线规划方法
CN108460485A (zh) * 2018-03-05 2018-08-28 重庆邮电大学 一种基于改进蚁群算法的旅行商问题求解方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
求解带时间窗车辆路径问题的动态混合蚁群优化算法;葛斌;韩江洪;魏臻;程磊;韩越;;模式识别与人工智能(第07期);全文 *

Also Published As

Publication number Publication date
CN111860754A (zh) 2020-10-30

Similar Documents

Publication Publication Date Title
CN111860754B (zh) 基于蚁群以及遗传算法的agv调度方法
WO2022262469A1 (zh) 一种基于博弈论的工业园区物流调度方法及系统
CN110443397B (zh) 一种订单配送方法
CN112418497B (zh) 一种面向制造物联的物料配送路径优化方法
CN107145971A (zh) 一种动态调整的快递配送优化方法
CN112561194B (zh) 一种混合流水车间生产与物流集成调度方法及系统
CN110782086A (zh) 一种用于救援的带无人机车辆配送路径优化方法及系统
CN109345091A (zh) 基于蚁群算法的整车物流调度方法及装置、存储介质、终端
CN112378415B (zh) 一种工器具的调度规划方法、装置及设备
CN113359702B (zh) 一种基于水波优化-禁忌搜索的智能仓库agv作业优化调度方法
CN115186930A (zh) 一种不确定取货需求下带多时间窗的同时取送货车辆路径问题规划方法及系统
CN116187896A (zh) 绿色车辆路径问题求解方法、装置、计算机设备以及介质
CN113743860A (zh) 车辆物资装载方法、系统、电子设备及存储介质
CN115907254A (zh) 一种基于进化的多目标强化学习的车辆路线规划方法
CN116227773A (zh) 一种基于蚁群算法的配送路径优化方法
CN114048924A (zh) 基于混合遗传算法的多配送中心选址-配送路径规划方法
CN115759917A (zh) 一种基于改进混合蚁群算法的物流路径规划方法
CN116151497A (zh) 一种带有取送货的多隔间车辆路径规划方法及系统
CN113313451A (zh) 一种基于改进布谷鸟算法的多目标优化物流调度方法
CN112700190B (zh) 一种改进的扫描法与遗传模拟退火法配送托盘物资的方法
CN113741418B (zh) 异构车机编队协同路径的生成方法和装置
CN115310676A (zh) 一种时变路网下的路径优化方法、装置及存储介质
CN114757394A (zh) 基于工作量平衡的物流车辆路径优化方法、系统和介质
CN115470651A (zh) 一种基于蚁群算法的带道路与时间窗车辆路径优化方法
CN111445100B (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
CB02 Change of applicant information

Address after: 214174 21 Yanxin Road, Yanqiao street, Huishan District, Wuxi City, Jiangsu Province

Applicant after: Wuxi Hongyi Intelligent Technology Co.,Ltd.

Address before: 214174 21 Yanxin Road, Yanqiao street, Huishan District, Wuxi City, Jiangsu Province

Applicant before: WUXI HONGYI INTELLIGENT TECHNOLOGY CO.,LTD.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant