CN107169594B - 一种车辆路径问题的优化方法及装置 - Google Patents

一种车辆路径问题的优化方法及装置 Download PDF

Info

Publication number
CN107169594B
CN107169594B CN201710276184.8A CN201710276184A CN107169594B CN 107169594 B CN107169594 B CN 107169594B CN 201710276184 A CN201710276184 A CN 201710276184A CN 107169594 B CN107169594 B CN 107169594B
Authority
CN
China
Prior art keywords
individuals
population
individual
optimal
initial population
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
CN201710276184.8A
Other languages
English (en)
Other versions
CN107169594A (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.)
Shenzhen University
Original Assignee
Shenzhen 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 Shenzhen University filed Critical Shenzhen University
Priority to CN201710276184.8A priority Critical patent/CN107169594B/zh
Publication of CN107169594A publication Critical patent/CN107169594A/zh
Application granted granted Critical
Publication of CN107169594B publication Critical patent/CN107169594B/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/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/047Optimisation of routes or paths, e.g. travelling salesman problem
    • 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Operations Research (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (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)
  • Feedback Control In General (AREA)

Abstract

本发明适用计算机技术领域,提供了一种车辆路径问题的优化方法及装置,该方法包括:生成车辆路径问题优化目标的初始种群,计算该种群中个体的适应度值,通过每个个体邻居范围内的最优个体引导雇佣蜂进行邻域搜索,生成第一种群,通过第一种群中部分优秀个体以及第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群,通过侦察蜂随机生成新个体替代第二种群中连续更新失败次数达到更新阈值的失败个体,当个体评价次数超过最大评价次数时,将第二种群中最优个体作为最优路径输出,否则将第二种群设置为初始种群,重新执行通过雇佣蜂进行的邻域搜索操作,从而加强了人工蜂群中蜜蜂的信息交流,提高了车辆路径问题的优化效率和优化效果。

Description

一种车辆路径问题的优化方法及装置
技术领域
本发明属于计算机技术领域,尤其涉及一种车辆路径问题的优化方法及装置。
背景技术
车辆路径问题是NP-hard组合优化问题,目前用于求解该类问题的算法主要有精确算法、传统启发式算法以及现代启发式算法。其中,精确算法可以找到最优解,但只适用于规模较小的问题,随着问题规模的增大以及计算机能力的限制,精确算法的求解效率降低,传统启发算法(例如扫描法、节约法和两阶段算法等)可以求解较大规模问题,提高问题的求解效率,但容易陷入局部最优,现代启发式算法(例如粒子群算法、模拟退火算法、禁忌搜索算法和遗传算法等)是现在使用最多的,可以求解大规模问题,通过规则提高了搜索效率,解的质量也优于传统启发式算法,但这些算法或多或存在收敛速度慢、容易陷入局部最优以及参数难调整的问题。
人工蜂群算法因为其结果简单易实现、参数少以及性能好的特点,一经提出便受到了各界学者的关注,但现有的人工蜂群算法存在大部分进化算法都存在的缺点——收敛速度慢,通过现有的人工蜂群算法解决车辆路径问题效率较低。
发明内容
本发明的目的在于提供一种车辆路径问题的优化方法及装置,旨在解决由于现有技术无法提供一种车辆路径问题优化的有效方法,导致车辆路径问题优化效率较低的问题。
一方面,本发明提供了一种车辆路径问题的优化方法,所述方法包括下述步骤:
根据车辆路径问题的优化目标,生成所述优化目标的初始种群,计算所述初始种群中每个个体的适应度值;
确定所述初始种群中个体的邻居范围并搜索每个邻居范围内的最优个体,通过所述每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索,生成第一种群;
在所述第一种群中选取部分优秀个体,同归根据预设的搜索公式以及公式使用概率,通过观察蜂对所述部分优秀个体进行邻域搜索,生成第二种群;
当所述第二种群中存在连续更新失败次数达到预设更新阈值的失败个体时,通过侦察蜂随机生成用来替代所述失败个体的新个体;
当个体评价次数超过预设最大评价次数时,将所述第二种群中的最优个体输出为所述车辆路径问题的最优路径,否则将所述第二种群设置为所述初始种群,重新执行确定所述初始种群中个体的邻居范围的操作。
另一方面,本发明提供了一种车辆路径问题的优化装置,所述装置包括:
初始化模块,用于根据车辆路径问题的优化目标,生成所述优化目标的初始种群,计算所述初始种群中每个个体的适应度值;
第一种群生成模块,用于确定所述初始种群中个体的邻居范围并搜索每个邻居范围内的最优个体,通过每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索,生成第一种群;
第二种群生成模块,用于在所述第一种群中选取部分优秀个体,通过根所述部分优秀个体和所述第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群;
个体替换模块,用于当所述第二种群中存在连续更新失败次数达到预设更新阈值的失败个体时,通过侦察蜂随机生成用来替代所述失败个体的新个体;以及
循环输出模块,用于当个体评价次数超过预设最大评价次数时,将所述第二种群中的最优个体输出为所述车辆路径问题的最优路径,否则将所述第二种群设置为所述初始种群,重新执行确定所述初始种群中个体的邻居范围的操作。
本发明在生成车辆优化目标的初始种群后,建立初始种群中每个个体的邻居范围,通过每个邻居范围内的最优个体引导雇佣蜂对初始种群中的个体进行邻域搜索,生成第一种群,从第一种群中选取部分优秀个体,通过部分优秀个体和第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群,当第二种群中存在连续更新失败次数达到预设更新阈值的失败个体时,通过侦察蜂随机生成用来替代失败个体的新个体,当个体评价次数超过预设最大次数时,将第二种群中的最优个体输出为车辆路径问题的最优路径,否则将第二种群设置为下一次迭代的初始种群,并重新执行确定初始种群中个体的邻居范围的操作,从而通过加强蜜蜂间的信息交流,使得雇佣蜂和观察蜂的搜索不是盲目独立的,有效地提高了人工蜂群算法的收敛速度和优化性能,进而有效地提高了车辆路径问题的优化效率和优化效果。
附图说明
图1是本发明实施例一提供的车辆路径问题的优化方法的实现流程图;
图2是本发明实施例一提供的通过雇佣蜂进行邻域搜索的示例图;
图3是本发明实施例二提供的车辆路径问题的优化装置的结构示意图;以及
图4是本发明实施例二提供的车辆路径问题的优化装置的优化结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
以下结合具体实施例对本发明的具体实现进行详细描述:
实施例一:
图1示出了本发明实施例一提供的车辆路径问题的优化方法的实现流程,为了便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在步骤S101中,根据车辆路径问题的优化目标,生成优化目标的初始种群,计算初始种群中个体的适应度值。
本发明实施例适用于求解车辆路径问题(VRP)的系统或平台。根据车辆路径问题的优化目标,例如路程最短、费用最少或者时间最短等优化目标,建立相应的数学模型,即优化模型。在建立好优化模型后,初始化相关参数和初始种群,其中,相关参数可包括种群个体数量、最大评价次数、个体连续更新失败的最大次数(更新阈值)、个体评价次数以及记录所有个体连续更新失败次数的集合等。
在本发明实施例中,可根据人工蜂群算法中的种群初始化公式,生成一组多维随机值,第i个多维随机值表示为Xi=(Xi,1,Xi,2,…),i=1,2,...,SN,SN为种群个体数量,由于人工蜂群初始化得到的每个多维随机值是一组连续的实数,而车辆行走的一条路径为一列整数序列,因此需要对每个多维随机值进行离散化,离散化后的多维随机值为不同的整数序列。将离散化后的多维随机值设置为初始种群中的个体,并根据构建的优化模型计算每个个体的适应度值(或目标函数值)。
具体地,人工蜂群中的种群初始化公式可表示为:
Figure BDA0001278325090000041
其中,Xij为第i个多维随机值的第j维分量,
Figure BDA0001278325090000042
为预设的第j维分量的最小值,
Figure BDA0001278325090000043
为预设的第j分量的最大值,rand0,1为0-1之间的随机值。
具体地,在对每个多维随机值进行离散化时,可分别将每个多维随机值中的分量进行排序,得到每个多维随机值中每个分量的序列号,因此每个多维随机值可对应一组序列号。例如,当一个多维随机值为(0.3,0.5,0.8,0.1)时,离散化后该多维随机值为(2,3,4,1),这些分量可依次表示路径节点2、3、4、1。
在步骤S102中,确定初始种群中个体的邻居范围并搜索每个邻居范围内的最优个体,通过每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索,生成第一种群。
在本发明实施例中,为初始种群中每个个体建立一个邻居范围,一个邻居范围可认为是每只雇佣蜂的局部视野,根据每个邻居范围中个体的适应度值搜索出该邻居范围内的最优个体,再通过这些最优个体对雇佣蜂的邻域搜索进行引导,使得雇佣蜂具有良好的方向性。邻域搜索后得到初始种群中每个个体对应的新个体,为了便于区分,将这些新个体称为第一新个体。将初始种群中每个个体与自身对应的第一新个体进行适应度值的比较,以对每个个体与对应的第一新个体间的优劣性进行评价,从而将每个个体与对应的第一新个体间的优秀个体保留下来,构成第一种群。其中,在雇佣蜂对一个个体进行一次邻域搜索后,可对个体评价次数进行加一操作。
具体地,可根据初始种群中当前的个体与初始种群中其它个体间的距离平均值,建立该个体的邻居范围,具体地,当个体Xi与个体Xj的距离dij<r·mdi时,可认为个体Xj是个体Xi的邻居,其中,r为预设参数,mdi为个体Xi与初始种群中其它个体的距离平均值。
具体地,通过每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索的公式为:
Vi,j=Xi,j+rand0,1·(Xnbest-Xi,j),其中,Vi,j为初始种群中第i个个体在邻域搜索后得到第一新个体Vi的第j维,Xi,j为初始种群中第i个个体Xi的第j维,Xnbest为第i个个体的邻居范围内的最优个体。在邻域搜索时,先确定当前个体变异的维度j,再通过邻域搜索公式对当前个体维度j的分量(例如图2中个体X的分量3)进行变异,得到变异后的分量(例如图2中第一新个体Xnew的分量6),查找变异后的分量在当前个体中的维度(例如图2中个体X的分量6所在维度),将该维度的分量设置为当前个体维度j的分量。
在步骤S103中,在第一种群中选取部分优秀个体,通过部分优秀个体和第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群。
在本发明实施例中,根据第一种群种个体的适应度值,在第一种群种按照预设的选取比例q选择B=q·SN个优秀个体,为便于区分,将这些优秀个体统称为部分优秀个体。将部分优秀个体的顺序打乱,并根据预设公式使用概率Pstr在部分优秀个体中选择前floor(Pstr·B)个个体作为第一部分个体,剩余的B-floor(Pstr·B)个个体作为第二部分个体,并通过该部分优秀个体以及第一种群中的最优个体,引导雇佣蜂采用预设第一公式对第一部分个体进行邻域搜索、预设第二公式对第二部分个体进行邻域搜索,生成部分优秀个体中每个个体对应的新个体,为了便于区分,将这些新个体称为第二新个体。根据部分优秀个体每个个体的适应度值和对应第二新个体的适应度值,保留部分优秀个体和对应的第二新个体之间的优秀个体,构成第二种群。从而通过在第一种群种选取部分优秀个体进行观察蜂的邻域搜索,有效地降低了通过观察蜂的搜索、选择压力,提高了观察蜂邻域搜索的效率。其中,观察蜂每对一个个体进行邻域搜索后,可对个体评价次数进行加一操作。
具体地,第一公式为Vi,j'=XPbest,j+φ·(Xi,j-Xk,j),第二公式为Vi,j'=Xbest,j+φ·(Xpbest,j-Xi,j),其中,XPbest为从部分优秀个体中随机选取的一个个体,Xk为从第二种群中随机选取的不同于个体Xi和个体XPbest的个体,Xbest为第二种群中的最优个体,Xi为部分优秀个体中的第i个个体,Vi'为部分优秀个体中的第i个个体Xi邻域搜索后得到的第二新个体,φ为[-1,1]中的一个随机实数。
优选地,当部分优秀个体中的个体比对应的第二新个体差时,计算该个体与对应的第二新个体间的适应度值差值(适应度值间的差值),当该个体的邻域搜索是采用第一公式完成时,将该适应度值差值累加到第一公式的提升度上,否则将该适应度值差值累加到第二公式的提升度上,继而根据第一公式的提升度和第二公式的提升度更新公式使用概率,从而有效提高通过观察蜂进行的邻域搜索的效果。具体地,根据第一公式的提升度和第二公式的提升度更新公式使用概率时,先计算第一公式和第二公式的平均提升度,当公式的平均提升度之和为0时,公式使用概率为0.5,否则,公式使用概率的公式为:
Pstr=0.1+0.8·aID1/(aID1+aID2),
其中,aID1为第一公式的平均提升度,aID1=ID1/floor(Pstr·B),aID2为第二公式的平均提升度,aID2=ID2/(B-floor(Pstr·B)),ID1为第一公式的提升度,ID2为第二公式的提升度。
在本发明实施例中,可将选取比例的倒数值设置为在第一种群中选取部分优秀个体的次数阈值,每完成一次上述的部分优秀个体选取、观察蜂邻域搜索和公式使用概率更新,就对当前的选取次数进行加一操作,并判断当前的选取次数是否超过次数阈值,当未超过时,将第二种群设置为第一种群,并跳转至在第一种群中选取部分优秀个体的操作。
在步骤S104中,第二种群中是否存在连续更新失败次数达到预设更新阈值的失败个体。
在本发明实施例中,雇佣蜂的邻域搜索和观察蜂的邻域搜索即分别对初始种群中个体进行更新和对第一种群中个体进行更新,当更新前的个体优于更新后的个体,可认为该个体此次更新失败。当第二种群中存在连续更新失败次数达到预设更新阈值的失败个体时执行步骤S105,否则执行步骤S106。
在步骤S105中,通过侦察蜂随机生成用来替代失败个体的新个体。
在本发明实施例中,在侦察蜂阶段,通过人工蜂群中的种群初始化公式生成一个多维随机值,并对该多维随机值进行离散化,得到一个新个体,用该新个体代替连续更新失败次数达到预设更新阈值的失败个体。其中,侦察蜂每更新一个个体后,可对个体评价次数进行加一操作。
在步骤S106中,个体评价次数是否超过预设最大评价次数。
在本发明实施例中,在初始化时,个体评价次数为0,每完成一次雇佣蜂阶段、观察蜂阶段的邻域搜索以及侦察蜂的个体更新,可对个体评价次数进行更新。当个体评价次数超过预设的最大评价次数时,认为车辆路径问题的优化完成,执行步骤S108,否则执行步骤S107。
在步骤S107中,将第二种群设置为初始种群,重新执行确定初始种群中个体的邻居范围的操作。
在本发明实施例中,当个体评价次数未超过最大评价次数时,将第二种群设置为下次迭代过程的初始种群,循环执行雇佣蜂、观察蜂的邻域搜索以及侦察蜂的个体更新等操作。
在步骤S108中,将第二种群中的最优个体输出为车辆路径问题的最优路径。
在本发明实施例中,当个体评价次数超过最大评价次数时,将第二种群中的最优个体作为车辆路径问题的最优路径输出,完成对车辆路径问题的优化。
在本发明实施例中,在雇佣蜂阶段通过邻居范围内的最优个体引导雇佣蜂的邻域搜索,加强了蜂群局部信息的交流,在观察蜂阶段通过部分优秀个体、整个种群中的最优个体引导观察蜂的邻域搜索,降低观察蜂搜索负担的同时加强了蜂群全局信息的交流,从而有效地提高了人工蜂群算法的收敛速度和优化性能,进一步有效地提高了车辆路径问题的优化效率。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,所述的存储介质,如ROM/RAM、磁盘、光盘等。
实施例二:
图3示出了本发明实施例二提供的车辆路径问题的优化装置的结构,为了便于说明,仅示出了与本发明实施例相关的部分,其中包括:
初始化模块31,用于根据车辆路径问题的优化目标,生成优化目标的初始种群,计算初始种群中个体的适应度值;
第一种群生成模块32,用于确定初始种群中个体的邻居范围并搜索每个邻居范围内的最优个体,通过每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索,生成第一种群;
第二种群生成模块33,用于在第一种群中选取部分优秀个体,通过部分优秀个体和第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群;
个体替换模块34,用于当第二种群中存在连续更新失败次数达到预设更新阈值的失败个体时,通过侦察蜂随机生成用来替代失败个体的新个体;
循环输出模块35,用于当个体评价次数超过预设最大评价次数时,将第二种群中的最优个体输出为车辆路径问题的最优路径,否则将第二种群设置为初始种群,重新执行确定初始种群中个体的邻居范围的操作。
优选地,如图4所示,初始化模块31包括模型建立模块411和种群生成模块412,其中:
模型建立模块411,用于根据车辆路径问题的优化目标,建立相应的优化模型,初始化优化模型的相关参数;以及
种群初始化模块412,用于种群生成模块,用于根据人工蜂群的种群初始化公式生成一组多维的随机值,离散化所有随机值并将离散化后的随机值设置为初始种群中的个体,根据优化模型计算初始种群中每个个体的适应度值。
优选地,第一种群生成模块32包括邻居确定模块421邻居搜索模块422第一邻域搜索模块423和第一更新模块424,其中:
邻居确定模块421,用于根据初始种群中每个个体间的距离,确定初始种群中每个个体的邻居范围;
邻居搜索模块422,用于根据每个邻居范围中所有个体的适应度值,搜索每个邻居范围的最优个体;
第一邻域搜索模块423,用于根据每个邻居范围内的最优个体引导雇佣蜂对初始种群中对应的个体进行邻域搜索,生成初始种群中每个个体对应的第一新个体;以及
第一更新模块424,用于根据初始种群每个个体的适应度值和对应的第一新个体的适应度值,将初始种群中每个个体与对应的第一新个体中的优秀个体设置为第一种群中的个体。
优选地,第二种群生成模块33包括部分选取模块431、第二邻域搜索模块432和第二更新模块433,其中:
部分选取模块431,用于根据第一种群中个体的适应度值和预设选取比例选取部分优秀个体,打乱部分优秀个体的顺序,将顺序打乱的部分优秀个体根据预设的公式使用概率划分为第一部分个体和第二部分个体;
第二邻域搜索模块432,用于通过部分优秀个体和第一种群中的最优个体,引导观察蜂采用预设的第一公式和第二公式分别对对第一部分个体、第二部分个体进行邻域搜索,生成部分优秀个体对应的第二新个体;以及
第二更新模块433,用于根据部分优秀个体与对应的第二新个体的适应度值,将部分优秀个体与对应的第二新个体中的优秀个体设置为第二种群的个体。
优选地,第二种群生成模块33还包括差值计算模块434、概率更新模块435和选取循环模块436,其中:
差值计算模块434,用于当部分优秀个体对应的第二新个体优于部分优秀个体时,计算部分优秀个体与对应的第二新个体的适应度值差值;
概率更新模块435,用于根据适应度值差值分别计算第一公式的提升度和第二公式的提升度,根据第一公式第二公式的提升度更新公式使用概率;以及
选取循环模块436,用于判断部分优秀个体的选取次数是否超过预设次数阈值,当选取次数未超过次数阈值时,将第二种群设置为第一种群并重新执行选取部分优秀个体的操作。
在本发明实施例中,在雇佣蜂阶段通过邻居范围内的最优个体引导雇佣蜂的邻域搜索,加强了蜂群局部信息的交流,在观察蜂阶段通过部分优秀个体、整个种群中的最优个体引导观察蜂的邻域搜索,降低观察蜂搜索负担的同时加强了蜂群全局信息的交流,从而有效地提高了人工蜂群算法的收敛速度和优化性能,进一步有效地提高了车辆路径问题的优化效率。
在本发明实施例中,车辆路径问题的优化装置的各模块可由相应的硬件或软件模块实现,各模块可以为独立的软、硬件模块,也可以集成为一个软、硬件模块,在此不用以限制本发明。本发明实施例中各模块的具体实施方式可参考前述实施例一中各步骤的描述,在此不再赘述。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种车辆路径问题的优化方法,其特征在于,所述方法包括下述步骤:
根据车辆路径问题的优化目标,生成所述优化目标的初始种群,计算所述初始种群中个体的适应度值;
确定所述初始种群中个体的邻居范围并搜索每个邻居范围内的最优个体,通过所述每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索,生成第一种群,通过每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索的公式为Vi,j=Xi,j+rand0,1·(Xnbest-Xi,j),其中,Vi,j为初始种群中第i个个体在邻域搜索后得到第一新个体Vi的第j维,Xi,j为初始种群中第i个个体Xi的第j维,Xnbest为第i个个体的邻居范围内的最优个体;
在所述第一种群中选取部分优秀个体,通过所述部分优秀个体和所述第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群;
当所述第二种群中存在连续更新失败次数达到预设更新阈值的失败个体时,通过侦察蜂随机生成用来替代所述失败个体的新个体;
当个体评价次数超过预设最大评价次数时,将所述第二种群中的最优个体输出为所述车辆路径问题的最优路径,否则将所述第二种群设置为所述初始种群,重新执行确定所述初始种群中个体的邻居范围的操作;
在所述第一种群中选取部分优秀个体,通过所述部分优秀个体和所述第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群的步骤,包括:
根据所述第一种群中个体的适应度值和预设选取比例选取所述部分优秀个体,打乱所述部分优秀个体的顺序,将所述顺序打乱后的部分优秀个体根据预设的公式使用概率划分为第一部分个体和第二部分个体;
通过所述部分优秀个体和所述第一种群中的最优个体,引导所述观察蜂采用预设第一公式Vi,j'=XPbest,j+φ·(Xi,j-Xk,j)和第二公式Vi,j'=Xbest,j+φ·(Xpbest,j-Xi,j)分别对所述第一部分个体、第二部分个体进行邻域搜索,生成所述部分优秀个体对应的第二新个体,其中,XPbest为从部分优秀个体中随机选取的一个个体,Xk为从第二种群中随机选取的不同于个体Xi和个体XPbest的个体,Xbest为第二种群中的最优个体,Xi为部分优秀个体中的第i个个体,Vi'为部分优秀个体中的第i个个体Xi邻域搜索后得到的第二新个体,φ为[-1,1]中的一个随机实数;
根据所述部分优秀个体与对应的第二新个体的适应度值,将所述部分优秀个体与对应的第二新个体中的优秀个体设置为所述第二种群的个体。
2.如权利要求1所述的方法,其特征在于,根据车辆路径问题的优化目标,生成所述优化目标的初始种群,计算所述初始种群中每个个体的适应度值的步骤,包括:
根据所述车辆路径问题的优化目标,建立相应的优化模型,初始化所述优化模型的相关参数;
根据人工蜂群的种群初始化公式生成一组多维的随机值,离散化所有随机值并将所述离散化后的随机值设置为所述初始种群中的个体,根据所述优化模型计算所述初始种群中每个个体的适应度值。
3.如权利要求1所述的方法,其特征在于,确定所述初始种群中个体的邻居范围并搜索每个邻居范围内的最优个体,通过所述每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索,生成第一种群的步骤,包括:
根据所述初始种群中每个个体间的距离,确定所述初始种群中每个个体的邻居范围;
根据所述每个邻居范围中所有个体的适应度值,搜索所述每个邻居范围的最优个体;
通过所述每个邻居范围内的最优个体引导所述雇佣蜂对所述初始种群中对应的个体进行邻域搜索,生成所述初始种群中每个个体对应的第一新个体;
根据所述初始种群每个个体的适应度值与对应的第一新个体的适应度值,将所述初始种群中每个个体与对应的第一新个体中的优秀个体设置为所述第一种群的个体。
4.如权利要求1所述的方法,其特征在于,在所述第一种群中选取部分优秀个体,通过所述部分优秀个体和所述第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群的步骤,还包括:
当所述部分优秀个体对应的第二新个体优于所述部分优秀个体时,计算所述部分优秀个体与对应的第二新个体的适应度值差值;
根据所述适应度值差值分别计算所述第一公式的提升度和所述第二公式的提升度,根据所述第一公式和第二公式的提升度更新所述公式使用概率;
判断所述部分优秀个体的选取次数是否超过预设次数阈值,当所述选取次数未超过所述次数阈值时,将所述第二种群设置为所述第一种群并重新执行选取所述部分优秀个体的操作。
5.一种车辆路径问题的优化装置,其特征在于,所述装置包括:
初始化模块,用于根据车辆路径问题的优化目标,生成所述优化目标的初始种群,计算所述初始种群中每个个体的适应度值;
第一种群生成模块,用于确定所述初始种群中个体的邻居范围并搜索每个邻居范围内的最优个体,通过每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索,生成第一种群,通过每个邻居范围内的最优个体引导雇佣蜂进行邻域搜索的公式为Vi,j=Xi,j+rand0,1·(Xnbest-Xi,j),其中,Vi,j为初始种群中第i个个体在邻域搜索后得到第一新个体Vi的第j维,Xi,j为初始种群中第i个个体Xi的第j维,Xnbest为第i个个体的邻居范围内的最优个体;
第二种群生成模块,用于在所述第一种群中选取部分优秀个体,通过所述部分优秀个体和所述第一种群中的最优个体引导观察蜂进行邻域搜索,生成第二种群;
个体替换模块,用于当所述第二种群中存在连续更新失败次数达到预设更新阈值的失败个体时,通过侦察蜂随机生成用来替代所述失败个体的新个体;以及
循环输出模块,用于当个体评价次数超过预设最大评价次数时,将所述第二种群中的最优个体输出为所述车辆路径问题的最优路径,否则将所述第二种群设置为所述初始种群,重新执行确定所述初始种群中个体的邻居范围的操作;
所述第二种群生成模块包括:
部分选取模块,用于根据所述第一种群中个体的适应度值和预设选取比例选取所述部分优秀个体,打乱所述部分优秀个体的顺序,将所述顺序打乱的部分优秀个体根据预设的公式使用概率划分为第一部分个体和第二部分个体;
第二邻域搜索模块,用于通过所述部分优秀个体和所述第一种群中的最优个体,引导所述观察蜂采用预设第一公式Vi,j'=XPbest,j+φ·(Xi,j-Xk,j)和第二公式Vi,j'=Xbest,j+φ·(Xpbest,j-Xi,j)分别对所述第一部分个体、第二部分个体进行邻域搜索,生成所述部分优秀个体对应的第二新个体,其中,XPbest为从部分优秀个体中随机选取的一个个体,Xk为从第二种群中随机选取的不同于个体Xi和个体XPbest的个体,Xbest为第二种群中的最优个体,Xi为部分优秀个体中的第i个个体,Vi'为部分优秀个体中的第i个个体Xi邻域搜索后得到的第二新个体,φ为[-1,1]中的一个随机实数;以及
第二更新模块,用于根据所述部分优秀个体与对应的第二新个体的适应度值,将所述部分优秀个体与对应的第二新个体中的优秀个体设置为所述第二种群的个体。
6.如权利要求5所述的装置,其特征在于,所述初始化模块包括:
模型建立模块,用于根据所述车辆路径问题的优化目标,建立相应的优化模型,初始化所述优化模型的相关参数;以及
种群初始化模块,用于根据人工蜂群的种群初始化公式生成一组多维的随机值,离散化所有随机值并将离散化后的随机值设置为所述初始种群中的个体,根据所述优化模型计算所述初始种群中每个个体的适应度值。
7.如权利要求5所述的装置,其特征在于,所述第一种群生成模块包括:
邻居确定模块,用于根据所述初始种群中每个个体间的距离,确定所述初始种群中每个个体的邻居范围;
邻居搜索模块,用于根据每个邻居范围中所有个体的适应度值,搜索所述每个邻居范围的最优个体;
第一邻域搜索模块,用于根据所述每个邻居范围内的最优个体引导所述雇佣蜂对所述初始种群中对应的个体进行邻域搜索,生成所述初始种群中每个个体对应的第一新个体;以及
第一更新模块,用于根据所述初始种群每个个体的适应度值和对应的第一新个体的适应度值,将所述初始种群中每个个体与对应的第一新个体中的优秀个体设置为所述第一种群中的个体。
8.如权利要求5所述的装置,其特征在于,所述第二种群生成模块还包括:
差值计算模块,用于当所述部分优秀个体对应的第二新个体优于所述部分优秀个体时,计算所述部分优秀个体与对应的第二新个体的适应度值差值;
概率更新模块,用于根据所述适应度值差值分别计算所述第一公式的提升度和所述第二公式的提升度,根据所述第一公式第二公式的提升度更新所述公式使用概率;以及
选取循环模块,用于判断所述部分优秀个体的选取次数是否超过预设次数阈值,当所述选取次数未超过所述次数阈值时,将所述第二种群设置为所述第一种群并重新执行选取所述部分优秀个体的操作。
CN201710276184.8A 2017-04-25 2017-04-25 一种车辆路径问题的优化方法及装置 Active CN107169594B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710276184.8A CN107169594B (zh) 2017-04-25 2017-04-25 一种车辆路径问题的优化方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710276184.8A CN107169594B (zh) 2017-04-25 2017-04-25 一种车辆路径问题的优化方法及装置

Publications (2)

Publication Number Publication Date
CN107169594A CN107169594A (zh) 2017-09-15
CN107169594B true CN107169594B (zh) 2020-08-04

Family

ID=59812316

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710276184.8A Active CN107169594B (zh) 2017-04-25 2017-04-25 一种车辆路径问题的优化方法及装置

Country Status (1)

Country Link
CN (1) CN107169594B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108538372A (zh) * 2018-03-07 2018-09-14 深圳市南山区人民医院 医疗服务价格调整方法、装置、服务器及存储介质
CN112232602A (zh) * 2020-11-19 2021-01-15 湘潭大学 一种大规模节点的物流配送路径优化方法及系统
CN116989819B (zh) * 2023-09-27 2024-01-02 湘江实验室 一种基于模型解的路径确定方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104299052A (zh) * 2014-10-12 2015-01-21 刘岩 一种基于人工蜂群和神经网络的建筑能耗预测方法
CN104808665A (zh) * 2015-04-16 2015-07-29 上海大学 基于多目标人工蜂群算法的多机器人路径规划方法
CN105704729A (zh) * 2016-01-22 2016-06-22 南京大学 一种采用改进的人工蜂群算法的无线传感器部署方法
CN105956706A (zh) * 2016-05-09 2016-09-21 苏州泛能电力科技有限公司 一种基于改进型人工蜂群算法的配电网无功优化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8660322B2 (en) * 2011-08-25 2014-02-25 King Saud University Passive continuous authentication method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104299052A (zh) * 2014-10-12 2015-01-21 刘岩 一种基于人工蜂群和神经网络的建筑能耗预测方法
CN104808665A (zh) * 2015-04-16 2015-07-29 上海大学 基于多目标人工蜂群算法的多机器人路径规划方法
CN105704729A (zh) * 2016-01-22 2016-06-22 南京大学 一种采用改进的人工蜂群算法的无线传感器部署方法
CN105956706A (zh) * 2016-05-09 2016-09-21 苏州泛能电力科技有限公司 一种基于改进型人工蜂群算法的配电网无功优化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A quick artificial bee colony (q ABC) algorithm and its performance on;D. Karaboga 等;《Applied Soft Computing》;20141231;全文 *
人工蜂群算法的改进与应用;董琴;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160225;全文 *
基于蜂群算法的随机需求车辆路径优化问题研究;王连稳;《电子世界》;20131231;全文 *

Also Published As

Publication number Publication date
CN107169594A (zh) 2017-09-15

Similar Documents

Publication Publication Date Title
CN111310999B (zh) 一种基于改进蚁群算法的仓库移动机器人路径规划方法
CN112036512B (zh) 基于网络裁剪的图像分类神经网络架构搜索方法和装置
CN107169594B (zh) 一种车辆路径问题的优化方法及装置
CN108882256B (zh) 一种无线传感器网络节点覆盖优化方法和装置
WO2019154215A1 (zh) 机器人运行路径生成方法、计算设备及存储介质
CN113422695B (zh) 一种提高物联网拓扑结构鲁棒性能的优化方法
US9674083B2 (en) Path calculation order deciding method, program and calculating apparatus
CN113111185A (zh) 一种用于知识图谱补全任务的关键路径构建方法和系统
CN110046774A (zh) 一种配电网故障定位方法、装置及相关设备
Huang et al. Ponas: Progressive one-shot neural architecture search for very efficient deployment
CN115169517A (zh) 一种改进的灰狼优化算法的方法
CN108133426B (zh) 一种社交网络链路推荐方法
CN110852435A (zh) 一种基于神经进化计算模型
CN104657901A (zh) 一种基于随机游走的标签传播社区发现方法
KR101565694B1 (ko) 무선 센서 네트워크에서 인공벌군집 클러스터링 설계 방법 및 시스템
CN111412795B (zh) 测试点设置方案生成方法及装置
CN116151384B (zh) 量子电路处理方法、装置及电子设备
CN111292062A (zh) 基于网络嵌入的众包垃圾工人检测方法、系统及存储介质
CN113590912B (zh) 融合节点相对位置和绝对度分布的跨社交网络对齐方法
CN108924196A (zh) 工业互联网绿色能源管理系统
CN109039698A (zh) 工业互联网智能服务处理方法、可读存储介质、终端
CN110048945B (zh) 一种节点移动性分簇方法及系统
Han et al. Research on limited buffer scheduling problems in flexible flow shops with setup times
CN116611527B (zh) 量子电路处理方法、装置及电子设备
CN110766126B (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