CN112884229B - 基于差分进化算法的大型公共场所人流引导路径规划方法 - Google Patents

基于差分进化算法的大型公共场所人流引导路径规划方法 Download PDF

Info

Publication number
CN112884229B
CN112884229B CN202110219913.2A CN202110219913A CN112884229B CN 112884229 B CN112884229 B CN 112884229B CN 202110219913 A CN202110219913 A CN 202110219913A CN 112884229 B CN112884229 B CN 112884229B
Authority
CN
China
Prior art keywords
path
population
vector
target
test vector
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
CN202110219913.2A
Other languages
English (en)
Other versions
CN112884229A (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.)
Sino Singapore International Joint Research Institute
Original Assignee
Sino Singapore International Joint Research Institute
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 Sino Singapore International Joint Research Institute filed Critical Sino Singapore International Joint Research Institute
Priority to CN202110219913.2A priority Critical patent/CN112884229B/zh
Publication of CN112884229A publication Critical patent/CN112884229A/zh
Application granted granted Critical
Publication of CN112884229B publication Critical patent/CN112884229B/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
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • 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
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Economics (AREA)
  • General Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Artificial Intelligence (AREA)
  • Health & Medical Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Marketing (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • General Business, Economics & Management (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • Mathematical Physics (AREA)
  • Development Economics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于差分进化算法的大型公共场所人流引导路径规划方法。本发明将路径规划问题表示为一个连续优化问题,利用差分进化算法结合多目标优化技术解决大型公共场所的人流引导路径规划问题,将路径方案编码为一组前后相关互联的基因序列,在差分进化过程中,基于人群行为仿真的结果对路径方案进行适应度评估,引导种群进行个体选择,将安全性和效率作为多目标选择操作的两个目标,应用经典的NSGA‑II算法进行种群个体选择,得到的路径方案既有着较好的安全性,又能引导人群高效通过场景,有效地克服了现有技术优化时评估准则过于单一的不足,减少了路径交叉,兼顾了安全性和效率。

Description

基于差分进化算法的大型公共场所人流引导路径规划方法
技术领域
本发明涉及路径规划和智能计算技术领域,具体涉及一种基于差分进化算法的大型公共场所人流引导路径规划方法。
背景技术
大型公共场所的人流引导问题是公共安全的重要问题。它要求在给定场景下提供路牌或路线指示,从而使得行人能快速、安全地到达各自的目的地,同时尽可能地避免不同路线之间的交叉导致的拥堵。
路径规划方法主要包括根据场景信息进行的全局路径规划和根据实时信息进行的局部路径规划。传统的路径规划方法包括Dijkstra算法和A*算法等,近年来也有仿生学方法被应用于路径规划,如蚁群算法、遗传算法等。这些方法往往针对给定的静态场景提供一组最短或较短的路径,但在像车站和商场这些人流密集、行人路线复杂的场景下,仅优化距离得到的路径往往不够安全,行走方向不同的人流之间的交叉可能会导致人群拥堵,而严重的拥堵则可能引发安全事故。在这类场景下,路径的距离或人群通过场景的效率往往不是唯一应该考虑的的优化目标,而现有的路径规划方法大多仅适用于优化单一目标,故而,在路径复杂的场景下,考虑路径的距离提高效率的同时减少路径交叉,从而提高路径安全性成为一个需要解决的问题。
此外,对于多出入口的场景,还需要将出口的选择作为规划内容的一部分,如何为走向各目的地的行人合理分配出口,既提高出口利用率,也控制出口的负载,同样是一个既关乎安全性也关乎效率的需要解决的问题。
发明内容
本发明的目的是为了解决现有技术中的上述缺陷,提供一种基于差分进化算法的大型公共场所人流引导路径规划方法。差分进化算法是一种基于种群的启发式搜索算法,种群中的每个种群个体对应一个解向量,称之为目标向量,其进化流程包括变异、杂交和选择等操作,生成测试向量,通过对目标向量和测试向量进行适应度评估,更新种群个体,经过不断地进化迭代,保留优良个体,淘汰劣质个体,引导搜索向最优解逼近。本发明提出的路径规划方法利用差分进化算法结合多目标优化技术解决路径规划问题,有效地克服了现有技术优化时评估准则过于单一的不足,减少了路径交叉,兼顾了安全性和效率,提高了其实用性。
本发明的目的可以通过采取如下技术方案达到:
一种基于差分进化算法的大型公共场所人流引导路径规划方法,其路径规划方法包括以下步骤:
S1、执行初始化操作:根据场景中的可行区域,随机生成大小为N的种群Q,所述种群包含N个种群个体;用前后基因互相关联的编码方式编码路径方案,为每个种群个体构建目标向量,路径方案由多条路径组成,对路径方案进行编码即对路径用前后基因互相关联的编码方式进行编码;基于人群行为的仿真结果计算所述目标向量的目标函数值,在仿真过程中,利用智能体模拟人群行为。
S2、执行个体交叉操作:为每个所述种群个体构建测试向量;基于人群行为仿真的结果计算测试向量的目标函数值;进行适应度评估,执行多目标选择操作更新种群个体的目标向量。
S3、执行局部搜索操作:更新种群个体测试向量;基于人群行为仿真的结果计算测试向量的目标函数值;进行适应度评估,执行多目标选择操作更新种群个体的目标向量。
S4、重复执行S2、S3步骤,直到达到预设迭代次数或路径方案的安全性和效率都满足预设需求。
S5、将最后一次迭代得到的路径方案作为最终路径方案,并对路径规划方案进行后期处理。
进一步地,步骤S1中目标向量包含规划方案的路线信息和终点信息,按如下方式构建:
Pi,k=[θi,k,1,θi,k,2,…,θi,k,j,…,θi,k,D,ei,k]
Figure BDA0002954352820000031
其中,Xi表示种群的第i个个体,θi,k,j表示第i个个体对应的路径方案中的第k条路径的第j个关键点与前一个点的连线与水平方向的夹角,Pi,k表示第i个个体对应的路径方案中的第k条路径,ei,k表示第i个个体第k条路径的终点,i=1,2,......,N,N为种群个体的数量、k=1,2,……,M,M为一个路径方案中的路径数量、j=1,2,......,D,D为每条路径的关键点个数,关键点不包括起点和终点。
每个种群个体被编码成一个M×(D+1)维的实数向量,每个个体Xi包含M条路径,每条路径用一个(D+1)维的向量表示,按前后基因互相关联的编码方式对路径方案中的每一条路径进行编码:路径在编码中简单地用关键点序列表示,每条路径中包含D个关键点,前D个实数是角的弧度值,表示该数对应的关键点与前一个关键点连线与水平线之间的夹角,第D+1个实数ei,k映射为具体目的地;θi,k1对应的关键点的前一个点是路径的起点,θi,kD对应的关键点的后一个点是ei,k,除了第一个关键点以外,其余关键点的前一个点均为关键点。后续通过小规模仿真,可以在关键点间得到更多间距更小的点;关键点为路径上依次按预设距离排列的点,即关键点在路径上均匀排列,但每条路径的关键点序列位于后面的关键点可能是冗余的,并不在路径上,故而路径的长度各不相同。
进一步地,目标函数值按如下方式计算:
f1=∑H(C)
f2=L/T
其中,C是以预设边长对场景对进行网格划分后得到的某个方形区域的中心点;H(C)是不同路径的行人经过点C的所在网格的频率的熵;f1则是所有网格中心点熵的和,是整个路径方案的熵;f1反映路径方案的安全性,f1越小,表明场景中同一位置上的人流方向不一样的情况越少,即人流交叉带来的安全隐患越小;L是在进行评估的仿真结束时仍然留在场景内的人数,T是进行评估的仿真过程中进入场景的总人数,f2反映在给定路径方案下人群通过场景的效率,f2越小,在给定路径方案下人群通过场景的效率就越高;两个目标函数值又被合称为适应度。
对于每个路径方案,进行预定时长的人群行为仿真,该仿真用多个智能体模拟行人产生人群动态,以计算两个目标函数值。在仿真过程中,各个智能体在不同时刻进入场景,并沿着路径方案中给出的路径走向各自的出口。
进一步地,所述步骤S2中的测试向量按如下方式构建:
Figure BDA0002954352820000041
其中,ui,m是所述测试向量第i个向量中的第m个元素,xi,m是种群目标向量的第i个向量中的第m个元素,F和CR是范围在(0,1)之间的参数,F是缩放比例因子,CR是个体交叉发生的概率,Xb和Xd是在种群中随机选择的除Xi外的两个不同的种群个体,xb,m和xd,m则是这两个种群个体的目标向量中的第m个元素。
进一步地,所述步骤S3的局部搜索中更新测试向量按如下方式操作:
对当前种群中的每个个体,根据目标向量的路径方案计算测试向量;选定对测试向量路径方案熵值贡献最大的路径,应用该路径的起点和终点,按前后基因互相关联的编码方式,重新随机生成一条路径;用生成的这条路径代替测试向量路径方案中对熵值贡献最大的那条路径,即淘汰掉带来最多安全隐患的路径。路径方案中的一条路径对路径方案熵值的贡献可以通过计算该路径方案中不含此条路径的其他路径的熵值来得到,若此计算得到的熵值越小,说明这条路径对整个路径方案熵值的贡献越大。
进一步地,所述步骤S5中后期处理为平滑处理,所述平滑处理为利用贝塞尔曲线对路径方案中的路径进行平滑,此外,还可以采用B-样条曲线来平滑路径方案中的路径。
进一步地,所述步骤S2和步骤S3中适应度评估,按如下方式操作:
比较种群个体目标向量和测试向量的适应度,适应度包括两个目标函数值,如果测试向量的两个目标函数值都小,即测试向量支配目标向量,则用测试向量取代目标向量;如果目标向量的两个目标函数值都小,即目标向量支配测试向量,则舍弃测试向量;如果测试向量有目标函数值比目标向量的目标函数值小,目标向量也有目标函数值比测试向量的目标函数值小,即目标向量和测试向量互不支配,此时将测试向量记录下来。
进一步地,所述步骤S2和步骤S3中多目标选择操作,按如下方式操作:
将记录下来的测试向量加入当前种群,执行NSGA-II算法的基于非支配排序和拥挤距离的多目标选择操作,逐个选择优秀的个体,使得当前种群的大小保持为N。除了NSGA-II算法,多目标选择还可以采用NSGA-III算法,NSGA-III为NSGA-II的改进版本,专门针对高维多目标优化。
进一步地,所述步骤S2和步骤S3中基于非支配排序和拥挤距离操作具体过程如下:
将所有种群个体根据支配关系分配为指定个前沿,前一个前沿支配后一个前沿,同一个前沿内的种群个体互相不支配;如果一个种群个体的两个目标函数值都优于另一个种群个体的目标函数值,则称前者支配后者;
计算前沿中每个种群个体的拥挤距离,并将种群个体按照拥挤距离从大到小排序,拥挤距离反映了每个种群个体周围的种群个体密度,为了保持种群中个体的多样性,选择操作将在互不支配的个体中优先选择周围种群个体密度稀疏的种群个体;
优先选择靠前的前沿中的种群个体加入新种群,在同一个前沿中则优先将拥挤距离大的种群个体加入新种群,直到新种群大小为N。
进一步地,所述步骤S3的局部搜索中根据目标向量的路径方案计算测试向量是将目标向量的路径方案直接赋值给测试向量。
进一步地,所述步骤S1中的人群行为仿真基于RVO2和A*算法实现。RVO2算法用于实现碰撞避免,也可以选择社会力模型来实现碰撞避免。A*算法用于实现点到点的路径导航,还可以采用Dijkstra算法、D*算法等进行路径导航。仿真分为两步:第一步是利用单个智能体根据关键点进行仿真,得到点间距离比关键点间距离更近的路径;第二步是用多个智能体模拟现实中的行人,智能体从各个入口持续不断地进入场景,沿着第一步仿真得到的路径,走向出口。
进一步地,所述步骤S2和步骤S3中的应用多目标选择操作更新种群个体的目标向量时,还可以采用线性加权的方式,即考虑对两个目标函数值线性加权,将其转化为单目标问题,可以根据实际应用场景有侧重性地考虑路径规划方案的安全性或效率。相较而言,多目标优化无需设计权值,具有参数较少的优点。
本发明相对于现有技术具有如下的优点及效果:
本发明将路径规划问题表示为一个连续优化问题,利用差分进化算法结合多目标优化技术解决大型公共场所的人流引导路径规划问题,采用一种新的前后基因相互关联的编码方式将路径方案编码为一组前后相关联的基因序列,在差分进化过程中,基于人群行为仿真的结果对路径方案进行适应度评估,引导种群更新进行种群个体选择,将安全性和效率作为多目标选择操作的两个目标,应用经典的NSGA-II基于非支配排序和拥挤距离的多目标选择操作进行种群个体选择,得到的路径方案既有着较好的安全性,又能引导人群高效通过场景,有效地克服了现有技术优化时评估准则过于单一的不足,减少了路径交叉,兼顾了安全性和效率。
附图说明
图1是本发明实施例中公开的基于差分进化算法的大型公共场所人流引导路径规划方法的流程图;
图2是本发明实施例中广州体育西路地铁站示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图对本发明的方法作更详细的描述。
实施例
对于大型公共场景,假设有M条不同的路径需要规划,M取决于人流有多少个主要的流向,将出、入口分组,每组入口与每组出口之间都有一条路径,人群可以根据自身所在的位置以及要抵达的出口,确定自己要跟随的路线指引。由于进行分组操作,一条路径可能对应多个入口和出口,每条规划路径均包含出口的选择以及中间的路线。我们的目标是在有效引导人群快速到达终点的同时,尽可能减少场景内不同路线上的人群之间的交汇,以提高安全性,因此,效率和安全性是本发明使用的多目标差分进化算法中的两个目标,对应了两个目标函数值,也即适应度。目标函数值按如下方式计算:
f1=∑H(C)
f2=L/T
其中,C是以预设边长对场景进行网格划分后得到的某个方形区域的中心点,网格划分将场景划分为多个大小相同的方形,便于计算熵,每个经过网格的模拟行人的智能体都会被用于计算该网格的熵;H(C)是在进行评估的仿真过程中,不同路线的模拟行人的智能体经过点C所在的网格的频率的熵,f1则是所有网格的熵的和,是整个路径方案的熵,f1反映路径方案的安全性,f1越小,表明场景中同一位置上的人流方向不一样的情况越少,即人流交叉带来的安全隐患越小;L是在进行评估的仿真结束时,仍然留在场景内的人数,T是进行评估的仿真过程中进入场景的总人数,f2反映在给定路径规划方案下人群通过场景的效率。
对于每个路径方案,将基于用于碰撞避免的RVO2算法和用于路径导航的A*算法,进行持续2分钟的仿真,该仿真用多个智能体模拟行人产生人群动态,以计算两个目标函数值,在仿真过程中,各个智能体在不同时刻进入场景,并沿着该算法给出的路径方案走向各自的出口。
本实施例中,依据两个目标函数值也即适应度对路径方案进行适应度评估,适应度评估方式如下:
比较种群个体目标向量和测试向量的适应度,适应度包括两个目标函数值,如果测试向量的两个目标函数值都小,即测试向量支配目标向量,则用测试向量取代目标向量;如果目标向量的两个目标函数值都小,即目标向量支配测试向量,则舍弃测试向量;如果测试向量有目标函数值比目标向量的目标函数值小,目标向量也有目标函数值比测试向量的目标函数值小,即目标向量和测试向量互不支配,此时将测试向量记录下来。
本实施例中,应用用多目标优化技术更新种群,按如下方式操作:
在每个个体都生成过一个测试向量以后,将记录下来的测试向量加入到种群中,用基于非支配排序和拥挤距离的选择操作淘汰差的个体,使得新种群大小保持为N,选择的过程如下:
首先,将所有种群个体根据支配关系分配为多个前沿,前一个前沿内的种群个体支配后一个前沿内的种群个体,同一个前沿内的种群个体互相不支配。如果一个种群个体的两个目标函数值都优于另一个种群个体,则前者支配后者;
计算前沿中每个种群个体的拥挤距离,并将种群个体按照拥挤距离从大到小排序,拥挤距离反映了每个种群个体周围的种群个体密度,为了保持种群中个体的多样性,选择操作将在互不支配的个体中优先选择周围种群个体密度稀疏的种群个体;
优先选择靠前的前沿中的种群个体加入新种群,在同一个前沿中则优先将拥挤距离大的种群个体加入新种群,直到新种群大小为N。
为了减小搜索空间,提高搜索效率,需要对不同路线的人群进行分组,并为种群个体的路线规划设置关键点间距离。关键点由进化算法搜索得到,关键点间的距离是一个预设的参数,参照路径整体的曲折程度和算力设定,关键点之间的距离越短,关键点的连线越接近方案对应的实际路径,如果方向要频繁发生改变,则距离要适当设置得小一些,但是为了提高搜索效率,距离不应当设置得过小,因为路径上的关键点如果比较多,搜索会变慢,适当减少关键点的个数,还可以简化目标向量,提高进化算法的搜索效率。
进行分组时,将出入口所在分组相同的人群分为一组。在执行进化算法的过程中,为了提高搜索效率,路径由数量较少、间隔距离较大的多个关键点表示,不足以在实际情形中引导行人,故在进行适应度评估的仿真之前,需要先用单个智能体来模拟单个行人进行简单仿真以得到间距较小的更多的点来表示路径。该仿真与用于适应度评估的仿真不同,单个智能体不会反映人群动态的安全性和效率,仅用于在路径上生成的更多的点,这些点的点间距比关键点之间的距离更小,该智能体根据路径方案中的关键点,结合传统的路径规划算法,如A*、Dijkstra和D*算法,对每条路径走出一条由间距比关键点之间的距离更小的更多的点组成的新路径,这条新路径在实际应用中比关键点的引导能力更强,可用于直接在地面绘制路线。在随后用于正式的适应度评估的仿真过程中,模拟多个行人的多个智能体将依据自己的起点和终点,跟随先前走出的路径点序列走向目的地。
图1给出了本发明实施例中基于差分进化算法的大型公共场所人流引导路径规划方法的流程图。下面就流程图的内容分步描述整个算法的具体实施方式:
S1、执行初始化操作:
根据场景中可以行走的区域,随机生成大小为N的种群Q,在生成过程中,如果发现当前点不在可以行走的区域内时,则重新生成这个点,直至它落到可以行走的区域内,确保规划路径上的每个关键点都落在场景中可以行走的区域内。对于体育西路地铁站的场景,种群大小N设为30。用前后基因互相关联的编码方式为每个种群个体构建目标向量,目标向量中包含了路线信息和终点信息。随后,利用前文所述的计算目标函数值的方法计算种群个体目标向量的目标函数值。目标向量按如下方式被编码成一个M×(D+1)维的实数向量:
Pi,k=[θi,k,1,θi,k,2,...,θi,k,j,...,θi,k,D,ei,k]
Figure BDA0002954352820000111
其中,Xi表示种群中的第i个个体,θi,k,j表示第i个个体对应的路径方案中第k条路径的第j个关键点与前一个点连线与水平方向的夹角,Pi,k表示第i个个体对应的路径方案中的第k条路径,ei,k表示第i个个体选择的终点,θi,k,1对应的关键点的前一个点是路径的起点,θi,k,D对应的关键点的后一个点是ei,k,除了第一个关键点以外,其余关键点的前一个点均为关键点,i=1,2,......,N,N为种群个体的数量、k=1,2,……,M,M为一个路径方案中的路径数量、j=1,2,......,D,D为一条路径的关键点个数,关键点不包括起点和终点。
在体育西路地铁站场景中,本实施例将M的值设为16,D的值设为10,表示一个路径方案包括16条路径,每条路径上有10个关键点。在最终的路径方案中,并非所有的关键点都会发挥作用,即一条路径对应的10个关键点的后几个点可能是冗余的。在图2的场景示意图中,左上的入口归为一组,右上的入口归为一组,中央上方偏右的两个入口归为一组,其余5个入口每个入口单独为一组,这样总共有8组入口,场景内部的电梯口按上、下位置分为2组,共有16条路径。
在体育西路地铁站场景中,网格设为边长为1平方米的正方形。在仿真过程中,有如下假设:
(a)模仿行人的智能体持续以近似相同的频率根据泊松分布从入口进入场景;
(b)所有智能体的目的地均为某个电梯口;
(c)一半智能体选择场景上部的电梯口作为目的地,另一半选择场景下部的电梯口作为目的地。
S2、执行个体交叉操作,更新种群个体:
对于种群中的每个目标向量,使其与当代其他个体用差分进化算法的随机交叉方式生成一个测试向量,以将其他个体的信息结合到自身,在修正该测试向量使其合法后,计算测试向量的目标函数值,进行适度度评估,执行NSGA-II的基于非支配排序和拥挤距离的多目标选择操作,更新种群个体。
S21、种群里的每个个体Xi按如下方式构建测试向量Ui
Figure BDA0002954352820000121
其中,ui,m、xi,m、xb,m和xd,m分别是Ui、Xi、Xb和Xd的第m个元素,F和CR是差分进化的参数,F表示缩放比例因子,CR是交叉概率,两者均是范围在(0,1)之间的参数,Xb和Xd是在种群中除了Xi外随机选择的两个不同的个体。本实施例将F的值设为0.1、CR的设为0.9。
S22、按前文所述的方法计算测试向量的目标函数值;
S22、按前文所述的方法进行适应度评估:
S23、按前文所述的方法用多目标优化技术更新种群个体的目标向量。
S3、执行局部搜索操作,更新种群个体:
种群个体交叉操作完成后,执行局部搜索操作。局部搜索的目的是通过对个体对应的路径方案中的某条路径进行调整,以得到更好的方案。本实施中,考虑到每条路径对应的基因段实际上是一组前后相关联的关键点,因此将一整条路径作为局部搜索对象,而不是将路径上的单独一个点作为搜索对象,以提高搜索效率。
S31、局部搜索操作的测试向量按如下步骤得到:
对于当前种群的每一个种群个体,将其目标向量路径方案的值直接赋值给测试向量;
选定对测试向量路径方案熵值贡献最大的路径,应用该路径的起点和终点,按前后基因互相关联的编码方式,重新随机生成一条路径;
用生成的这条路径代替测试向量路径方案中对熵值贡献最大的那条路径,即淘汰掉带来最多安全隐患的路径。路径方案中的一条路径对路径方案熵值的贡献可以通过计算该路径方案中不含此条路径的其他路径的熵值来得到,仅用其他路径计算得到的熵值越小,说明这条路径对整个路径方案熵值的贡献越大。
S32、按前文所述的方法计算测试向量的目标函数值;
S33、按前文所述的方法进行适应度评估:
S34、按前文所述的方法用多目标优化技术更新种群个体的目标向量。
S4、重复执行S2、S3步骤,直到达到预设迭代数,或者所得方案在人流引导上的安全性和效率都满足了预设需求。
S5、将最后一次迭代得到的路径方案作为最终路径方案,利用贝塞尔曲线对路径方案中的路径进行曲线平滑等后期处理,得到可以实际应用的路径。
特别地,在进行人群行为仿真时,碰撞避免还可以用社会力模型,社会力模型更接近人的真实行为,但RVO2模型有较为完善的库,可以方便地实现并行运算,效率更高;路径导航算法还可以采用Dijkstra算法、D*算法等,但A*算法比Dijkstra算法和D*算法的效率高。在执行多目标选择操作更新种群个体时,还可以采用NSGA-II的改进版本NSGA-III算法,也可以采用线性加权的方式对两个目标函数值线性加权,将多目标选择操作转化为单目标操作,以便根据实际应用场景有侧重性地考虑路径规划方案的安全性或效率;此外,在对所得路径进行平滑处理时,还可以应用B-样条曲线对路径进行平滑。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (8)

1.一种基于差分进化算法的大型公共场所人流引导路径规划方法,其特征在于,所述路径规划方法包括以下步骤:
S1、执行初始化操作:随机生成大小为N的种群Q,所述种群Q包含N个种群个体;用前后基因互相关联的编码方式编码路径方案,为每个种群个体构建目标向量;基于人群行为仿真的结果计算所述目标向量的目标函数值;
S2、执行个体交叉操作:为每个所述种群个体构建测试向量;基于人群行为仿真的结果计算所述测试向量的目标函数值;进行适应度评估,执行多目标选择操作更新种群个体的目标向量;
S3、执行局部搜索操作:更新每个种群个体的测试向量;基于人群行为仿真的结果计算测试向量的目标函数值;进行适应度评估,执行多目标选择操作更新种群个体的目标向量;
S4、重复执行步骤S2、S3,直到达到预设迭代次数或路径方案的安全性和效率都满足预设需求;
S5、将最后一次迭代得到的路径方案作为最终路径方案,并对路径规划方案进行后期处理;
其中,所述步骤S1中目标向量按如下方式构建:
Pi,k=[θi,k,1,θi,k,2,...,θi,k,j,...,θi,k,D,ei,k]
Figure FDA0003641406290000011
其中,Xi表示种群中的第i个个体,θi,k,j表示第i个个体对应的路径方案中的第k条路径的第j个关键点与前一个点的连线与水平方向的夹角,Pi,k表示第i个个体对应的路径方案中的第k条路径,ei,k表示第i个个体在第k条路径的终点,i=1,2,......,N,N为种群个体的数量、k=1,2,......,M,M为路径方案中的路径数量、j=1,2,......,D,D为每条路径的关键点个数;
每个种群个体被编码成一个M×(D+1)维的实数向量:每个个体Xi包含M条路径;每条路径用一个(D+1)维的向量表示,按前后基因互相关联的编码方式对路径方案中的每一条路径进行编码:路径在编码中简单地用关键点序列表示,前D个实数是角的弧度值,表示该数对应的关键点与前一个点连线与水平线之间的夹角,第D+1个实数ei,k映射为具体目的地;θi,k,1对应的关键点的前一个点是路径的入口组的位置,θi,k,D对应的关键点的后一个点是ei,k对应的目的地,除了第一个关键点以外,其他关键点的前一个点均为关键点;所述关键点在路径上依次按预设距离排列。
2.根据权利要求1所述的基于差分进化算法的大型公共场所人流引导路径规划方法,其特征在于,所述目标函数值按如下方式计算:
f1=∑H(C)
f2=L/T
其中,C是以预设边长对场景进行网格划分后得到的某个方形区域的中心点;H(C)是不同路径的行人经过点C所在方形区域的频率的熵;f1是所有网格中心点熵的和,表示整个路径方案的熵;L是在进行评估的仿真结束时仍然留在场景内的人数,T是进行评估的仿真过程中进入场景的总人数,f2反映给定路径方案下人群通过场景的效率;两个目标函数值被合称为适应度。
3.根据权利要求1所述的基于差分进化算法的大型公共场所人流引导路径规划方法,其特征在于,所述步骤S2中测试向量按如下方式构建:
Figure FDA0003641406290000031
其中,ui,m是测试向量第i个向量中的第m个元素,xi,m是目标向量第i个向量中的第m个元素,m=1,2,......,M×(D+1),F和CR是范围在(0,1)之间的参数,F是缩放比例因子,CR是差分进化中个体交叉的概率,Xb和Xd是在种群中随机选择的除Xi外的两个不同的个体,xb,m和xd,m则是这两个个体的目标向量中的第m个元素。
4.根据权利要求1所述的基于差分进化算法的大型公共场所人流引导路径规划方法,其特征在于,所述步骤S3中更新测试向量按如下方式操作:
对每个种群个体,根据目标向量的路径方案计算测试向量;选定对测试向量路径方案熵值贡献最大的路径,应用该路径的起点和终点,按前后基因互相关联的编码方式,重新随机生成一条路径;用生成的这条路径代替测试向量路径方案中对熵值贡献最大的那条路径。
5.根据权利要求1所述的基于差分进化算法的大型公共场所人流引导路径规划方法,其特征在于,所述步骤S5中所述后期处理为平滑处理。
6.根据权利要求1至5任一所述的基于差分进化算法的大型公共场所人流引导路径规划方法,其特征在于,所述步骤S2和步骤S3中适应度评估,按如下方式操作:
比较所述种群个体目标向量和测试向量的适应度,适应度包括两个目标函数值,如果测试向量的两个目标函数值都小,即测试向量支配目标向量,则用测试向量取代目标向量;如果目标向量的两个目标函数值都小,即目标向量支配测试向量,则舍弃测试向量;如果测试向量有目标函数值比目标向量的目标函数值小,目标向量也有目标函数值比测试向量的目标函数值小,即目标向量和测试向量互不支配,此时将测试向量记录下来。
7.根据权利要求1至5任一所述的基于差分进化算法的大型公共场所人流引导路径规划方法,其特征在于,所述步骤S2和步骤S3中多目标选择操作,按如下方式操作:
将记录下来的测试向量加入当前种群,执行基于非支配排序和拥挤距离的多目标选择操作,逐个选择优秀的个体,使得当前种群的大小保持为N。
8.根据权利要求7所述的基于差分进化算法的大型公共场所人流引导路径规划方法,其特征在于,所述基于非支配排序和拥挤距离操作,具体过程如下:
将所有种群个体根据支配关系分配为指定个前沿,前一个前沿支配后一个前沿,同一个前沿内的种群个体互相不支配;如果一个种群个体的两个目标函数值都优于另一个种群个体,则称前者支配后者;
计算前沿中每个种群个体的拥挤距离,并按照拥挤距离从大到小排序;
优先选择靠前的前沿中的种群个体加入新种群,在同一个前沿中则优先将拥挤距离大的种群个体加入新种群,直到新种群大小为N。
CN202110219913.2A 2021-02-26 2021-02-26 基于差分进化算法的大型公共场所人流引导路径规划方法 Active CN112884229B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110219913.2A CN112884229B (zh) 2021-02-26 2021-02-26 基于差分进化算法的大型公共场所人流引导路径规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110219913.2A CN112884229B (zh) 2021-02-26 2021-02-26 基于差分进化算法的大型公共场所人流引导路径规划方法

Publications (2)

Publication Number Publication Date
CN112884229A CN112884229A (zh) 2021-06-01
CN112884229B true CN112884229B (zh) 2022-12-20

Family

ID=76054826

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110219913.2A Active CN112884229B (zh) 2021-02-26 2021-02-26 基于差分进化算法的大型公共场所人流引导路径规划方法

Country Status (1)

Country Link
CN (1) CN112884229B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113673695B (zh) * 2021-07-07 2023-07-21 华南理工大学 基于新型特征自动构造的人群行为规则自动提取方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105430707A (zh) * 2015-11-03 2016-03-23 国网江西省电力科学研究院 一种基于遗传算法的无线传感网多目标优化路由方法
CN112231968A (zh) * 2020-09-09 2021-01-15 山东师范大学 基于深度强化学习算法的人群疏散仿真方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AT514588B1 (de) * 2013-08-29 2015-02-15 Tech Universität Wien Verfahren zur Steuerung eines Fahrzeuges
CN106647262B (zh) * 2016-11-28 2020-01-21 中国人民解放军国防科学技术大学 一种面向敏捷卫星多目标任务规划的差分进化方法
CN110598863A (zh) * 2019-09-05 2019-12-20 南宁师范大学 协同进化的多目标差分进化方法
CN111767789B (zh) * 2020-05-13 2024-03-26 北京交通大学 一种基于多载体智能引导的人群疏散方法及系统
CN112132312B (zh) * 2020-08-14 2022-08-23 蓝海(福建)信息科技有限公司 一种基于进化多目标多任务优化的路径规划方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105430707A (zh) * 2015-11-03 2016-03-23 国网江西省电力科学研究院 一种基于遗传算法的无线传感网多目标优化路由方法
CN112231968A (zh) * 2020-09-09 2021-01-15 山东师范大学 基于深度强化学习算法的人群疏散仿真方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于遗传算法的最短路径规划;林煦涵等;《电子世界》;20170423(第08期);179-180 *

Also Published As

Publication number Publication date
CN112884229A (zh) 2021-06-01

Similar Documents

Publication Publication Date Title
CN105070042B (zh) 一种交通预测的建模方法
CN108197739B (zh) 一种城市轨道交通乘客流量预测方法
Ciaffi et al. Feeder bus network design problem: A new metaheuristic procedure and real size applications
CN107464005B (zh) 一种用于车辆预约的拓展型路径规划方法
CN111562785B (zh) 一种群集机器人协同覆盖的路径规划方法及系统
CN108021686B (zh) 一种快速匹配电子地图中公交线路与道路网络的方法
CN101853294B (zh) 一种基于遗传算法的多模式多标准路径选择方法
CN111612252A (zh) 大规模应急设施的自动选址方法、装置及可读存储介质
CN111238521A (zh) 一种无人驾驶车辆的路径规划方法及系统
CN109270927A (zh) 道路数据的生成方法及装置
CN112013829A (zh) 基于多目标优化的多uav/ugv协同长时作业路径规划方法
CN106228265A (zh) 基于改进粒子群优化的总拖期运输计划调度算法
Dong et al. A novel chaotic particle swarm optimization algorithm for parking space guidance
CN112884229B (zh) 基于差分进化算法的大型公共场所人流引导路径规划方法
CN113516277A (zh) 一种基于路网动态定价的网联智能交通路径规划方法
CN115270506A (zh) 一种人群沿楼梯上行的通行时间预测方法及系统
CN115545359A (zh) 一种面向复杂建筑火灾的动态智能疏散方法及装置
JPH07129882A (ja) 道路交通信号制御適正化装置、およびそれを用いた道路交通流管制システム
CN114969234A (zh) 一种面向全量城市路网的向量表征方法
CN111340303A (zh) 基于新型混合蛙跳算法的旅行商路线规划方法
CN110598948A (zh) 针对旅游景点多点归源的全排列路径规划的提速方法
CN114217609A (zh) 一种满足多约束条件的电网线路走廊带智能划分方法
CN112070275A (zh) 一种轨道交通路网协同限流的方法和装置及设备
Lee et al. Comparative analysis of machine learning algorithms to urban traffic prediction
CN113705891B (zh) 基于mra-bas-bp算法的城市商办综合体建筑停车需求预测方法

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