一种基于多粒子群协同演化的人车混合疏散仿真优化方法

一种基于多粒子群协同演化的人车混合疏散仿真优化方法
【专利摘要】本发明提供了一种基于多粒子群协同演化的人车混合疏散仿真优化方法，主要通过粒子群优化算法解决大型建筑物和路网集成环境中的人车混合疏散问题，利用粒子的运动过程模拟疏散个体，包括人员、车辆的混合疏散过程，并且采用基于信息素的多粒子群通信机制模拟人员和车辆之间的交互和信息共享，克服了粒子群算法只考虑了粒子之间的交互而忽略了环境对粒子运动过程的影响等不足，能够更好地模拟人车混合疏散过程中不同交通对象的竞争和协作，形成满足多个目标需求的疏散优化方案，提供合理高效的决策依据。
【专利说明】
-种基于多粒子群协同演化的人车混合疏散仿真优化方法
技术领域
[0001] 本发明属于智能计算和运筹学的交叉应用领域，具体地说是一种基于多粒子群协 同演化的人车混合疏散仿真优化方法，主要通过粒子群优化算法解决大型建筑物和路网集 成环境中的人车混合疏散问题，利用粒子的运动过程模拟疏散个体，包括人员、车辆的混合 疏散过程，并且采用基于信息素的多粒子群通信机制模拟人员和车辆之间的交互和信息共 享。
【背景技术】
[0002] 大型建筑物和路网集成环境下的人车混合疏散问题就是将紧急情况下聚集在大 型建筑物中的大量人员从事发地建筑撤离并进入周边路网，与路网中的移动对象共同在路 网疏散直到完全离开危险区域或到达安全区域。对于单独的建筑物人员疏散或者路网上的 车辆疏散问题，国内外研究者已经提出了大量的方法，但是紧急事件的影响范围一般是由 较小的范围逐渐扩大至较大范围内的移动对象，包括建筑物及其周边路网上的人员及车 辆。目前国内外对紧急情况下建筑物内的人员疏散研究较多，仿真模型和疏散策略研究也 相对成熟，但基于建筑物和交通路网的多模式混合交通疏散的研究相对不足。
[0003] 粒子群优化算法（Particle Swarm Optimization，PS0)是由Kennedy于 1995年提 出的一种基于群智能的优化算法。它源于对鸟类捕食行为的模拟，假设一群鸟要捜索一定 范围内的一块食物，且在整个捜索区域中只有一块食物。所有的鸟都不知食物在何处，但它 们知道自己当前位置距离食物还有多远。那么找到食物的最简单有效的方法就是捜寻并移 动到目前离食物最近的鸟的周围区域，并重复上述过程，最终所有的鸟都可W通过运种整 体协作的方式到达食物所在的位置。
[0004] 粒子群算法的运行机理不是依赖个体的自然进化规律，而是对生物群体的社会行 为进行模拟，它源于对生物群体行为和人类社会行为的研究。在生物群体中，个体与个体之 间、个体与群体之间都存在相互影响、相互作用，表现为一种群体中的信息共享机制。粒子 群算法就是受到运样一种社会行为的启发而提出的，个体间通过信息共享可W相互借鉴经 验，W促进整个群体的演化。由于粒子群算法思想简单、控制参数少、实现容易，已经成为研 究多目标优化问题的重要方法，被广泛应用于模式识别、信号处理、决策支持、系统控制、模 糊系统控制、神经网络训练等领域。
[0005] 虽然粒子群优化算法具有简单快速的优点，但是有两个主要的缺点限制其性能： 一方面，PSO对于疏散问题的研究主要集中于仿真，即粒子模拟疏散个体的行为，但是应用 粒子群算法对疏散过程进行优化的研究不多，其难点在于PSO模拟的一次疏散过程为一个 疏散方案，而本发明的目的是根据疏散目标得到较优的疏散方案，因此需要利用粒子群模 拟多次W逼近最优方案，该过程的演化需要用到疏散网络上的信息素；另一方面，传统粒子 群只考虑了粒子间的学习，而忽略了粒子所处的网络对粒子运动的影响。

【发明内容】

[0006] 本发明的目的在于:为了克服已有技术的不足，即传统疏散模型大多注重疏散过 程的模拟，缺少对疏散过程的优化，从而提出一种基于信息素的协同演化粒子优化算法，一 方面使用多粒子群来模拟人员和车辆运两种类型的交通对象及其之间的交互，另一方面利 用疏散网络的信息素来得到最优疏散方案。多粒子群协同演化优化算法可W模拟不同群体 的运动过程，通过群体间的通信机制既体现环境对疏散过程的影响，又可W通过共享和交 互信息来达到种群间的相互作用和协同疏散，从而可W得到符合实际需求的多目标优化疏 散方案。
[0007] 为了达到上述目的，本发明采用的技术方案是:一种基于多粒子群协同演化的人 车混合疏散仿真优化方法，其特征在于，包括如下步骤： 步骤1.根据疏散场景初始化疏散网络拓扑图：定义AfetroWKW,巧为疏散网络，其中W =(1，2，…，W为疏散网络中的道路交叉点，定义M为车辆允许经过的节点，MHTfen知， Rank2, ??? , Rankn], N^ C i￥ ； 步骤2.设置模型参数:包括粒子群的粒子数量#、最大迭代次数7_^^?；惯性系数《、粒 子群学习因子Cl, C2, C3、信息素浓度初值化信息素浓度挥发系数芦，阳:1]; 步骤3.设置人员数量M、车辆数量#2,生成规模为M的粒子群1和规模为#2的粒子群2， 即粒子群1的粒子数量#= M，粒子群2的粒子数量#=化； 步骤4.初始化人员和车辆位置，W及信息素浓度;将代表人员和车辆的粒子随机放置 在疏散网络节点上，其中，代表车辆的粒子只能放置在属于W麻]节点上，将初始位置加入每 个粒子的路径，初始化所有粒子的速度为0,初始化疏散网络节点上的信息素浓度为步骤2 设置的化 步骤5.设置迭代计数器gen = 1; 步骤6.对于粒子群1和粒子群2，分别执行步骤7-步骤20; 步骤7.从当前粒子群中随机选择一个粒子作为第一个，设仁1; 步骤8.判断第A个粒子的状态是否为"成功撤离"，若是则执行步骤17,否则执行步骤 9; 步骤9.统计当前时刻疏散网络中每个节点的粒子数，记为numj(t)，计算t时刻所有 "未撤离"状态的粒子所覆盖的面积，记为areaU)，根据公式（1)计算淋1刻粒子群的密度c/ U):
(1) 步骤10.计算当前粒子群中所有粒子的适应值:根据每个粒子的当前位置，根据公式 (2)计算每个粒子的适应值：
(2) 式中，joartic沁为第A个粒子，边为粒子妓Ij出口m的距离； 步骤11.根据步骤10得到的所有粒子适应值，找出当前粒子群中适应值巧 最小的粒子所在的位置记为当前粒子群最优粒子位置根据第A个粒子的路径信息，找 出其路径上距离出口最近的位置记为粒子个体最优位置根据当前时刻疏散网络中每 (3) 个位置节点的信息素浓度，将最大信息素浓度所在位置记为巧。^ pAsr。?。。。； 击驢12.按公古（3)巧新笛Zr个紛子的巧底， KU)和KU + 1)分别为粒子在t时刻和t + 1时刻的速度，1U)为粒子在淋]刻的位置， W为步骤2设置的惯性系数，Cl, C2, C3为步骤2设置的学习因子，ii, n, n为(〇，1)之间的 随机数，Msr。。。。。分别为步骤11得到的粒子个体最优位置、粒子群最优粒子 位置、最大信息素浓度所在位置； 步骤13.更新第A个粒子的位置，并将该位置加入粒子的路径如，记录第A个粒子通 过上一位置节点的时间，W及该粒子从上一位置节点到最新节点的距离:根据步骤12得到 的速度按公式(4)更新第A个粒子在下一时刻的位置，将粒子移动到最新位置； 著稼+ 1)二錢左》夺巧安手0 (4) 1U)和Kt + 1)分别为粒子在t时刻和t + 1时刻的位置； 步骤14.根据步骤13得到的最新位置判断该位置是否为出口位置，若是，则代表第A个 粒子已到达出口，标记该粒子状态为"成功撤离"，否则执行步骤15; 步骤15.粒子计数加1，即A = A + 1，若^ <薦兑明还未遍历完该粒子群的所有粒子， 返回步骤9,否则表明完成了一个时间步的捜索，执行步骤16; 步骤16.判断当前粒子群中所有粒子的状态是否为"成功撤离"，若是则表明本轮迭代 完成，将本轮迭代的时间步数记为r，执行步骤17,否则返回步骤7; 步骤17.根据步骤13得到的当前粒子群中每个粒子的路径，根据步骤13得到的路径上 通过每个节点的时间和距离计算每个粒子通过整个路径的总时间和总距离，记为粒子疏散 时间和路径长度，按公式巧)计算信息素浓度：
(5) S(gen)为第gen次迭代位置节点i上的信息素浓度，，0为步骤2设置的信息素浓度挥发 系数，A诗为位置节点i在个体层次上的信息素浓度增量，这均根据公式(6)计算：
6) 其中，0为步骤2设置的信息素浓度初值，扣W (gen)为步骤17得到的第gen次迭代中第 A个粒子的路径长度，min扣^(1，…，gen)为第A个粒子在第1次到第gen次迭代中的最短 路径； 步骤18.计算疏散目标值:根据步骤17得到的粒子疏散时间，按照公式(7)计算所有粒 子的疏散时间总和巧，根据步骤9得到的各个时刻粒子密度，按照公式(8)计算总密度值巧， 巧和巧即为本发明所优化的两个疏散目标值：
(7) (8) t*为步骤17得到的第A个粒子的疏散时间，#为当前粒子群的粒子数量，c/( t)为步骤9 得到的时刻粒子群的密度，r为步骤16得到的本轮迭代的时间步数； 步骤19 .根据步骤18得到的疏散目标值凡和巧按照W下化reto支配关系确定当前粒 子群的化reto最优解集该化reto最优解集为迭代过程中满足凡和仍最小的若干个 疏散方案的集合； 在多目标优化问题中，Pareto支配是一个非常重要的概念，它是解的一种偏序关系，定 义如下： 化'6*〇支配:对于决策空间冲的任意两个决策向量a和6, a支配6当且仅当向量*3的 目标值都不大于向量6的目标值，且至少存在一个目标，向量a对应的目标值严格小于向 量6对应的目标值； 化reto最优解：对决策空间f中的任一向量a，口为化reto最优，当且仅当决策空间f 中存在化reto支配向量幻;的决策向量； 化reto最优解集：由决策空间f中的所有化reto最优解构成，并定义该集合对应的目 标向量集合为化reto最优前沿； 步骤20.根据步骤19得到的化reto最优解集t，对t中的每一个解，即疏散 方案，按照公式(9)更新其方案中所有粒子经过路径上节点位置的信息素浓度：
(9) 芦为步骤2设置的信息素浓度挥发系数，Ar/为位置节点i在全局层次上的信息素浓 度增量;Ar/按照公式（10)计算：
(10) 0为步骤2设置的信息素浓度初值，如为步骤13得到的第A个粒子的路径，诗为步 骤13得到的第A个粒子通过节点i的时间； 步骤21.粒子群间通信：每个粒子群中的粒子除了通过各自的信息素合作捜索出口 夕h还通过粒子群间的通信来共享信息素，W改进捜索过程;每个粒子群按照公式（11)~公 式（16)来更新自己的信息素：
(11) (12) (13) (14) (15) (16)
4,馬叫解]分别是粒子群I和粒子群2的通信因子，是衡量粒子群I和粒子群2相互影 响的参数;AU),泌U)分别为由步骤9得到的t时刻粒子群1和粒子群2的密度，r为步骤16 得到的本轮迭代的时间步数； 步骤22.迭代计数器洗础日1，即gen = gen +1; 步骤23.若肾^/?《7_#心是步骤2中设置的最大迭代次数，将所有粒子放置于步 骤4设置的各自的初始位置节点上，重复步骤6-22，否则，迭代结束执行步骤24; 步骤24.输出化reto最优解集。
[0008] 作为进一步优选的，所述步骤9中t时刻粒子群的密度(t)为当前时刻粒子群体 中未到达出口的粒子个数与未撤离粒子在疏散网络中所占面积之比。
[0009] 作为进一步优选的，所述步骤12中粒子速度更新方式为粒子通过对当前速度、粒 子个体最优位置、粒子群最优粒子位置、最大信息素浓度所在位置运几个因素的综合学习 得出。
[0010] 作为进一步优选的，所述步骤17中信息素浓度在个体层次上的增量由粒子自身迭 代过程中的个体最优路径决定。
[0011] 作为进一步优选的，所述步骤19中化reto最优解集为同时满足粒子群疏散时间总 和凡最短和粒子总密度值巧最小的疏散方案。
[0012] 作为进一步优选的，所述步骤21中行人粒子群1和车辆粒子群2间的通信规则为一 个粒子群通过引入另一个粒子群的信息素来改进捜索过程。
[0013] 作为进一步优选的，本方法的参数设置为：人员数量似1 = 2000,车辆数量化= 1000,最大迭代次数200,惯性系数W= 1，粒子群学习因子幻=3, C2= 3, C3= 2， 信息素浓度初值0= 1000,信息素浓度挥发系数= 0.3。
[0014] 本发明与现有技术相比具有W下优点： 1、本发明采用多粒子群算法模拟行人和车辆的疏散过程，能够更好地模拟人车混合疏 散过程中的人员群体和车辆行为，W及不同对象之间的交互，粒子的运动除了受到自身经 验和群体最优个体的影响之外，同时还向邻域中的最优个体学习，运种局部和全局并存的 学习机制既模拟了疏散过程中的从众行为、小群体现象，又加速了寻找安全出口的进程。
[0015] 2、本发明采用基于信息素的多粒子群协同演化算法优化人车混合疏散过程，各种 群在各自独立进化的同时相互间共享和交互信息，表现在各种群利用从外界获得的信息W 指导自身的捜索过程，同时还把自身寻优经验与其他种群分享，从而使得整个系统协同进 化，使得人员和车辆能够协作达到最优。
[0016] 3、本发明与基于单一粒子群的方法相比，不但能够模拟行人和车辆的疏散过程， 而且能够通过粒子群间的通信对疏散方案进行优化，且拥有较好的收敛速度和收敛效果。
[0017] 4、本发明与其他方法相比，能够得到更好的疏散方案，并且同时满足疏散时间和 疏散密度两个目标要求，从而使得疏散路径更加合理，疏散方案更加符合实际需求。
【附图说明】
[0018] 图1是本发明的实验流程图。
[0019] 图2是本发明中多粒子群协同演化算法流程图。
[0020] 图3是本发明的疏散网络图。
[0021 ]图4是本发明的密度示意图。
[0022] 图5是本发明中多粒子群协同演化算法收敛图。
[0023] 图6是本发明方法和粒子群优化算法的密度变化曲线对比图。
[0024] 图7是本发明方法和粒子群优化算法在疏散路径拥堵程度上的比较。
【具体实施方式】
[0025] 为了更好地理解本发明，下面结合实施例进一步阐明本发明的内容，但本发明的 内容不仅仅局限于下面的实施例。本领域技术人员可W对本发明作各种改动或修改，运些 等价形式同样在本申请所列权利要求书限定范围之内。
[0026] 图1和图2是本发明的实验流程图和多粒子群协同演化优化算法流程图，从中可W 看出具体实现过程如下： 步骤1.根据疏散场景初始化疏散网络拓扑图：定义AfetroWKW,巧为疏散网络，其中W =(1，2，…，W为疏散网络中的道路交叉点，定义M为车辆允许经过的节点，MHTfen知， Ranh，???，Rank。}，N', N。
[0027] 步骤2.设置模型参数:包括粒子群的粒子数量狀最大迭代次数肢性系数 W、粒子群学习因子幻，C2, C3、信息素浓度初值化信息素浓度挥发系数只，卢€巧4]。
[0028] 步骤3.设置人员数量化、车辆数量#2,生成规模为化的粒子群1和规模为始的粒子 群2,即粒子群1的粒子数量#=化，粒子群2的粒子数量#=始。
[0029] 步骤4.初始化人员和车辆位置，W及信息素浓度;将代表人员和车辆的粒子随机 放置在疏散网络节点上，其中，代表车辆的粒子只能放置在属于W脚节点上，将初始位置加 入每个粒子的路径，初始化所有粒子的速度为0，初始化疏散网络节点上的信息素浓度为步 骤2设置的化
[0030] 步骤5.设置迭代计数器gen = 1。
[0031] 步骤6.对于粒子群1和粒子群2,分别执行步骤7-步骤20。
[0032] 步骤7.从当前粒子群中随机选择一个粒子作为第一个，设A = 1。
[0033] 步骤8.判断第A个粒子的状态是否为"成功撤离"，若是则执行步骤17,否则执行 步骤9。
[0034] 步骤9.统计当前时刻疏散网络中每个节点的粒子数，记为numjU)，计算t时刻 所有"未撤离"状态的粒子所覆盖的面积，记为areaU)，根据公式（1)计算t时刻粒子群的 酱度t)。
[0035]
(I) 步骤10.计算当前粒子群中所有粒子的适应值:根据每个粒子的当前位置，根据公式 (2)计算每个粒子的适应值：
(2) 式中，jOarticJ饼为第A个粒子，垃为粒子A到出口m的距离。
[0036] 步骤11.根据步骤10得到的所有粒子适应值，找出当前粒子群中适应值/ 最小的粒子所在的位置记为当前粒子群最优粒子位置根据第A个粒子的 路径信息，找出其路径上距离出口最近的位置记为粒子个体最优位置根据当前时刻 疏散网络中每个位置节点的信息素浓度，将最大信息素浓度所在位置记为pAsra。。。。。
[0037] 步骤12.按公式(3)更新第A个粒子的速度：
(3) KU)和KU + 1)分别为粒子在t时刻和t + 1时刻的速度，ZU)为粒子在t时刻的位 置，W为步骤2设置的惯性系数，Cl, C2, C3为步骤2设置的学习因子，ii, n为(0，1)之 间的随机数，Aux Msr。。。。。分别为步骤11得到的粒子个体最优位置、粒子群最优 粒子位置、最大信息素浓度所在位置。
[0038] 步骤13.更新第A个粒子的位置，并将该位置加入粒子的路径如，记录第A个粒 子通过上一位置节点的时间，W及该粒子从上一位置节点到最新节点的距离:根据步骤12 得到的速度按公式(4)更新第A个粒子在下一时刻的位置，将粒子移动到最新位置； J佑手巧二黎皆+!鳴夺1) (4) Kt)和Kt + 1)分别为粒子在t时刻和t + 1时刻的位置。
[0039] 步骤14.根据步骤13得到的最新位置判断该位置是否为出口位置，若是，则代表 第A个粒子已到达出口，标记该粒子状态为"成功撤离"，否则执行步骤15。
[0040] 步骤15.粒子计数加1，即A = A + 1，若A <薦兑明还未遍历完该粒子群的所有 粒子，返回步骤9,否则表明完成了一个时间步的捜索，执行步骤16。
[0041] 步骤16.判断当前粒子群中所有粒子的状态是否为"成功撤离"，若是则表明本轮 迭代完成，将本轮迭代的时间步数记为r，执行步骤17,否则返回步骤7。
[0042] 步骤17.根据步骤13得到的当前粒子群中每个粒子的路径，根据步骤13得到的路 径上通过每个节点的时间和距离计算每个粒子通过整个路径的总时间和总距离，记为粒子 疏散时间和路径长度，按公式巧)计算信息素浓度：
(5) S(gen)为第gen次迭代位置节点i上的信息素浓度，声为步骤2设置的信息素浓度挥发 系数，这呀为位置节点i在个体层次上的信息素浓度增量，这种根据公式(6)计算： (6) 其中，0为步骤2设置的信息素浓度初值，扣A (gen)为步骤17得到的第gen次迭代中 第A个粒子的路径长度，min扣^(1，…，gen)为第A个粒子在第1次到第gen次迭代中的最 短路径。
[0043] 步骤18.计算疏散目标值:根据步骤17得到的粒子疏散时间，按照公式(7)计算所 有粒子的疏散时间总和巧，根据步骤9得到的各个时刻粒子密度，按照公式(8)计算总密度 值巧，巧和巧即为本发明所优化的两个疏散目标值：
(7) (8) 为步骤17得到的第A个粒子的疏散时间，#为当前粒子群的粒子数量，c/U)为步骤9得 到的时刻粒子群的密度，7?步骤16得到的本轮迭代的时间步数。
[0044] 步骤19.根据步骤18得到的疏散目标值巧和巧按照W下化reto支配关系确定当前 粒子群的化reto最优解集该化reto最优解集为迭代过程中满足巧和巧最小的若干 个疏散方案的集合。
[0045] 在单目标优化问题中，最优解通常是唯一确定的，而对多目标优化问题而言，由于 各个目标函数之间是冲突的，不可折衷的，因而多目标优化问题的最优解集通常是一个集 合，运是多目标优化问题相对于单目标优化问题的本质差别。
[0046] 在多目标优化问题中，Pareto支配(dominate)是一个非常重要的概念，它是解的 一种偏序关系，定义如下： 化reto支配:对于决策空间f中的任意两个决策向量a和6, a支配6当且仅当向量a 的目标值都不大于向量6的目标值，且至少存在一个目标，向量对应的目标值严格小于 向量6对应的目标值； 化reto最优解：对决策空间中的任一向量？，U'为化reto最优，当且仅当决策空间f 中存在化reto支配向量'巧的决策向量； 化reto最优解集：由决策空间f中的所有化reto最优解构成，并定义该集合对应的目 标向量集合为F*areto最优前沿(Pareto Optimal Rront)。
[0047] 步骤20.根据步骤19得到的化reto最优解集，对冲的每一个解(疏 散方案），按照公式(9)更新其方案中所有粒子经过路径上节点位置的信息素浓度：
(9) 芦为步骤2设置的信息素浓度挥发系数，金'r/为位置节点i在全局层次上的信息素浓 度增量;A r/按照公式（10)计算： (10) 0为步骤2设置的信息素浓度初值，如为步骤13得到的第A个粒子的路径，社为步骤 13得到的第A个粒子通过节点i的时间。
[004引步骤17和步骤20表明了疏散网络中信息素浓度的更新受到两个方面因素的影响， 一个是每个粒子根据自身迭代过程中的个体最优路径更新本路径的信息素，另一方面，根 据Pareto最优解集中的非支配解集更新该方案中的所有路径上的信息素。
[0049] 步骤21.粒子群间通信:每个粒子群中的粒子除了通过各自的信息素合作捜索出 口外，还通过粒子群间的通信来共享信息素，W改进捜索过程;每个粒子群按照公式（11)~ 公式（16)来更新自己的信息素：
(11) (12) (13) (14) V 斟扮别是粒子群1和粒子群2的通信因子，是衡量粒子群1和粒子群2相互影 响的参数;AU),泌U)分别为由步骤9得到的t时刻粒子群1和粒子群2的密度，r为步骤16 得到的本轮迭代的时间步数。
[0050] 步骤22.迭代计数器沪础日1，即洗n =洗n +1。
[0051] 步骤23.若洗n《7_倘义7_似心是步骤2中设置的最大迭代次数，将所有粒子放置 于步骤4设置的各自的初始位置节点上，重复步骤6-22，否则，迭代结束执行步骤24。
[0化2]步骤24.输出化reto最优解集。
[0053]为了验证本发明的有效性，我们将本发明与基本粒子群算法进行了仿真比较。在 仿真实验过程中，选取武汉体育中屯、体育馆及其周边区域作为实验场景。本发明的参数设 置为:人员数量M= 2000,车辆数量始=1000,最大迭代次数200,惯性系数W = 1，粒子群学习因子Ci= 3，C2= 3，C3= 2,信息素浓度初值0= 1000,信息素浓度挥发系数 只=0.3。
[0054] 图3是本发明所使用场景的疏散网络图，其中标记为圆点的节点为仅允许行人通 行，标记为方形的节点为允许行人或车辆通行。
[0055] 图4是本发明中关于密度定义的示意图，其中密度为场景中未疏散个体的数量与 未疏散个体所占面积之比，从图中可W看出，疏散个体越分散，密度越低，本发明W疏散时 间和密度作为优化目标，得到的优化方案既能够保证时间最短，又能降低疏散过程中拥堵 程度，从而得到安全、快速的疏散方案。
[0056] 图5是本发明中多粒子群协同演化算法的收敛图，从图中可W看出，随着迭代次数 的增加，疏散总时间和总密度两个目标值都在降低，说明算法能够保证收敛，且算法在150 代左右达到最优。
[0057] 图6是本发明方法和基本粒子群优化算法的密度变化比较，可W看出在疏散时间 上，本发明的方法比基本粒子群优化算法明显缩短，并且密度在整个疏散过程中保持较低 的水平，说明本发明的方法能够使疏散个体在疏散全过程中尽可能地分散，从而避免产生 拥堵。
[0058] 图7是本发明方法和基本粒子群优化算法在疏散路径上的比较，图中路段由细到 粗变化，代表路段的通过人数/车辆数由低到高，从图中标注区域可W看出本发明的方法明 显优于基本粒子群优化算法，说明本发明的方法得到的疏散路径更加合理。
[0059] 表1列出了本发明的多粒子群协同演化优化算法(CEPSO)和基本粒子群算法(PSO) 在不同迭代次数下的实验结果，可W看出，随着迭代次数的增加，本发明方法(CEPSO)和基 本粒子群算法(PSO)在总疏散时间、总密度和时间步上都趋于下降趋势，但本发明方法在运 几个评价指标上均明显优于本粒子群算法，说明本发明的方法具有更快的收敛速度。
[0060] 表1 CEPSO与PSO方法优化结果比较
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
【主权项】
1. 一种基于多粒子群协同演化的人车混合疏散仿真优化方法，其特征在于，包括如下 步骤： 步骤1.根据疏散场景初始化疏散网络拓扑图：定义Λ/eirorK#，幻为疏散网络，其中# ={1，2，…，/3}为疏散网络中的道路交叉点，定义屬为车辆允许经过的节点， Rankz, ··· , Rankn], C 步骤2.设置模型参数:包括粒子群的粒子数量#、最大迭代次数惯性系数《、粒 子群学习因子C1, Ω，C3、信息素浓度初值认信息素浓度挥发系数P[04]; 步骤3.设置人员数量M、车辆数量#2，生成规模为屬的粒子群1和规模为杨的粒子群2， 即粒子群1的粒子数量I=麗，粒子群2的粒子数量if =屬； 步骤4.初始化人员和车辆位置，以及信息素浓度;将代表人员和车辆的粒子随机放置 在疏散网络节点上，其中，代表车辆的粒子只能放置在属于节点上，将初始位置加入每 个粒子的路径，初始化所有粒子的速度为0,初始化疏散网络节点上的信息素浓度为步骤2 设置的仏 步骤5.设置迭代计数器= 1; 步骤6.对于粒子群1和粒子群2，分别执行步骤7-步骤20; 步骤7.从当前粒子群中随机选择一个粒子作为第一个，设A= 1; 步骤8.判断第1 个粒子的状态是否为"成功撤离"，若是则执行步骤17,否则执行步骤 9； 步骤9.统计当前时刻疏散网络中每个节点的粒子数，记为?)，计算?时刻所有"未 撤离"状态的粒子所覆盖的面积，记为area( ?)，根据公式（1)计算?时刻粒子群的密度忒?):(1) 步骤10.计算当前粒子群中所有粒子的适应值:根据每个粒子的当前位置，根据公式 (2)计算每个粒子的适应值：式中，pariide*为第女个粒子，为粒子女至Ij出口《的距离； 步骤11.根据步骤10得到的所有粒子适应值，找出当前粒子群中适应值/Xpariidei) 最小的粒子所在的位置记为当前粒子群最优粒子位置根据第1 个粒子的路径信息，找 出其路径上距离出口最近的位置记为粒子个体最优位置怂^，根据当前时刻疏散网络中每 个位置节点的信息素浓度，将最大信息素浓度所在位置记为; 击3娶11桉公忒审新笛Zr个鈴子的逋麽.(3) KU)和KU + 1)分别为粒子在?时刻和? + 1时刻的速度，ZU)为粒子在?时刻的位置， ?为步骤2设置的惯性系数，Cl，C2, C3为步骤2设置的学习因子，ii，r2, ?为(0，1)之间的 随机数，/^^、/^^、/^1#_?^分别为步骤11得到的粒子个体最优位置、粒子群最优粒子 位置、最大信息素浓度所在位置； 步骤13.更新第1 个粒子的位置，并将该位置加入粒子的路径paih，记录第1 个粒子通 过上一位置节点的时间，以及该粒子从上一位置节点到最新节点的距离:根据步骤12得到 的速度按公式(4)更新第1 个粒子在下一时刻的位置，将粒子移动到最新位置； '% + i 二赚+ ?% + 1> (4) 1U)和IU + 1)分别为粒子在?时刻和? + 1时刻的位置； 步骤14.根据步骤13得到的最新位置判断该位置是否为出口位置，若是，则代表第女个 粒子已到达出口，标记该粒子状态为"成功撤离"，否则执行步骤15; 步骤15.粒子计数加1，即I = 1 + 1，若I < Af说明还未遍历完该粒子群的所有粒 子，返回步骤9,否则表明完成了一个时间步的搜索，执行步骤16; 步骤16.判断当前粒子群中所有粒子的状态是否为"成功撤离"，若是则表明本轮迭代 完成，将本轮迭代的时间步数记为Γ，执行步骤17,否则返回步骤7; 步骤17.根据步骤13得到的当前粒子群中每个粒子的路径，根据步骤13得到的路径上 通过每个节点的时间和距离计算每个粒子通过整个路径的总时间和总距离，记为粒子疏散 时间和路径长度，按公式(5)计算信息素浓度：石(糾/3)为第ge/3次迭代位置节点i上的信息素浓度，P为步骤2设置的信息素浓度挥发 系数，Δ#为位置节点i在个体层次上的信息素浓度增量，根据公式(6)计算：其中，0为步骤2设置的信息素浓度初值，c/ist (ge/3)为步骤17得到的第ge/3次迭代中第女 个粒子的路径长度，min t/is*(l，···，ge/3)为第1 个粒子在第1次到第ge/3次迭代中的最短 路径； 步骤18.计算疏散目标值:根据步骤17得到的粒子疏散时间，按照公式(7)计算所有粒 子的疏散时间总和丹，根据步骤9得到的各个时刻粒子密度，按照公式(8)计算总密度值/?， 丹和/?即为本发明所优化的两个疏散目标值：以为步骤17得到的第1 个粒子的疏散时间，#为当前粒子群的粒子数量，t/U)为步骤9得 到的时刻粒子群的密度，7?步骤16得到的本轮迭代的时间步数； 步骤19.根据步骤18得到的疏散目标值F1和/?按照以下Pareto支配关系确定当前粒子 群的Pareto最优解集A2?_sei，该Pareto最优解集为迭代过程中满足/7I和/?最小的若干个疏 散方案的集合； 在多目标优化问题中，Pareto支配是一个非常重要的概念，它是解的一种偏序关系，定 义如下： Pare to支配:对于决策空间冲的任意两个决策向量a和/?，a：支配当且仅当向量a：的 目标值都不大于向量祕勺目标值，且至少存在一个目标，向量对应的目标值严格小于向量 M寸应的目标值； Pare to最优解:对决策空间冲的任一向量，'《为Pareto最优，当且仅当决策空间Γ中 存在Pareto支配向量_￡?的决策向量； Pareto最优解集：由决策空间冲的所有Pareto最优解构成，并定义该集合对应的目标 向量集合为Pareto最优前沿； 步骤20 .根据步骤19得到的Pareto最优解集对中的每一个解，即疏散 方案，按照公式(9)更新其方案中所有粒子经过路径上节点位置的信息素浓度：(9) P为步骤2设置的信息素浓度挥发系数，Ar/为位置节点i在全局层次上的信息素浓度 增量;Δ r/按照公式(10)计算：歧/步骤2设置的信息素浓度初值，押咖为步骤13得到的第1 个粒子的路径，$为步骤 13得到的第1个粒子通过节点i'的时间； 步骤21.粒子群间通信：每个粒子群中的粒子除了通过各自的信息素合作搜索出口 外，还通过粒子群间的通信来共享信息素，以改进搜索过程;每个粒子群按照公式（11)~公 式（16)来更新自己的信息素：馬(AefOg分别是粒子群1和粒子群2的通信因子，是衡量粒子群1和粒子群2相互影响 的参数;必⑴，必⑴分别为由步骤9得到的埘刻粒子群1和粒子群2的密度，7?步骤16得到 的本轮迭代的时间步数； 步骤22.迭代计数器gao加1，即gao = gao +1; 步骤23.若糾/3彡是步骤2中设置的最大迭代次数，将所有粒子放置于步 骤4设置的各自的初始位置节点上，重复步骤6-22，否则，迭代结束执行步骤24; 步骤24.输出Pareto最优解集。2. 根据权利要求1所述的基于多粒子群协同演化的人车混合疏散仿真优化方法，其特 征在于，所述步骤9中?时刻粒子群的密度c/U)为当前时刻粒子群体中未到达出口的粒子个 数与未撤离粒子在疏散网络中所占面积之比。3. 根据权利要求1所述的基于多粒子群协同演化的人车混合疏散仿真优化方法，其特 征在于，所述步骤12中粒子速度更新方式为粒子通过对当前速度、粒子个体最优位置、粒子 群最优粒子位置、最大信息素浓度所在位置这几个因素的综合学习得出。4. 根据权利要求1所述的基于多粒子群协同演化的人车混合疏散仿真优化方法，其特 征在于，所述步骤17中信息素浓度在个体层次上的增量由粒子自身迭代过程中的个体最优 路径决定。5. 根据权利要求1所述的基于多粒子群协同演化的人车混合疏散仿真优化方法，其特 征在于，所述步骤19中Pareto最优解集为同时满足粒子群疏散时间总和丹最短和粒子总密 度值/?最小的疏散方案。6. 根据权利要求1所述的基于多粒子群协同演化的人车混合疏散仿真优化方法，其特 征在于，所述步骤21中行人粒子群1和车辆粒子群2间的通信规则为一个粒子群通过引入另 一个粒子群的信息素来改进搜索过程。7. 根据权利要求1所述的基于多粒子群协同演化的人车混合疏散仿真优化方法，其特 征在于，本方法的参数设置为：人员数量#1 = 2000,车辆数量#2 = 1000,最大迭代次数Γ_ 侧Z = 200,惯性系数《= 1，粒子群学习因子〇 = 3，C2 = 3，C3 = 2,信息素浓度初值0 = 1000,信息素浓度挥发系数P = 0.3。
【文档编号】G06Q10/04GK106022510SQ201610306666
【公开日】2016年10月12日
【申请日】2016年5月11日
【发明人】宗欣露, 王春枝, 叶志伟, 刘伟, 徐慧, 陈宏伟, 蒋颖丽
【申请人】湖北工业大学