CN107665282B - 基于异构情绪感染模型的人群疏散仿真方法及系统 - Google Patents
基于异构情绪感染模型的人群疏散仿真方法及系统 Download PDFInfo
- Publication number
- CN107665282B CN107665282B CN201710940633.4A CN201710940633A CN107665282B CN 107665282 B CN107665282 B CN 107665282B CN 201710940633 A CN201710940633 A CN 201710940633A CN 107665282 B CN107665282 B CN 107665282B
- Authority
- CN
- China
- Prior art keywords
- individual
- panic
- emotional
- value
- individuals
- 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.)
- Expired - Fee Related
Links
- 230000002996 emotional effect Effects 0.000 title claims abstract description 109
- 238000000034 method Methods 0.000 title claims abstract description 51
- 208000015181 infectious disease Diseases 0.000 title claims abstract description 50
- 238000004088 simulation Methods 0.000 title claims abstract description 28
- 230000008451 emotion Effects 0.000 claims abstract description 76
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000009877 rendering Methods 0.000 claims abstract description 15
- 230000000694 effects Effects 0.000 claims abstract description 14
- 230000008447 perception Effects 0.000 claims description 24
- 238000005315 distribution function Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 15
- 230000036651 mood Effects 0.000 claims description 14
- 230000001186 cumulative effect Effects 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 8
- 230000000007 visual effect Effects 0.000 claims description 7
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 238000009825 accumulation Methods 0.000 claims description 5
- 230000006735 deficit Effects 0.000 claims description 4
- 238000000638 solvent extraction Methods 0.000 claims description 4
- 230000004888 barrier function Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 claims description 3
- 208000001431 Psychomotor Agitation Diseases 0.000 claims description 2
- 206010038743 Restlessness Diseases 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 2
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 208000027534 Emotional disease Diseases 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了基于异构情绪感染模型的人群疏散仿真方法及系统，首先，设计了一种基于Durupinar模型的异构情绪感染计算方法，在该方法中，提出了一种基于分组的情绪感染计算方法，考虑个体间距离和关系两个因素对分组的影响，并利用人群分组算法对人群进行分组，研究分组对情绪感染计算的影响；在此基础上，提出基于特征的恐慌程度计算方法，根据个性、性别、年龄的特征差异将个体特殊化，并利用上述特征进行情绪恐慌程度计算。其次，将上述异构情绪感染模型与人群运动耦合，结合相对速度障碍法，利用情绪驱动人群运动。最后，使用真实感渲染的方法模拟人群疏散效果。本发明能较真实直观的观察人群运动中的情绪感染现象。
Description
技术领域
本发明涉及一种人群疏散仿真方法，尤其涉及基于异构情绪感染模型的人群疏散仿真方法及系统。
背景技术
运用人群运动仿真研究人群在突发情况下人群的运动特征与规律，能够帮助设计者对建筑物或者公共场所的人流疏散通道进行综合分析与评价，辅助建立紧急情况下人员的疏散方案，以及合理的建筑布局，以减少群体事件中拥挤、踩踏等不良事件的发生，有利于紧急情况下的人群疏散工作。在人群疏散过程中，个体往往会聚集成小的分组，而分组内部也存在个体差异，但现有的情绪感染模型大多没有考虑分组和个体差异对情绪的影响，个体之间的感染是同构的，即个体之间的感染速度相同，无法逼真地模拟人群运动过程。
发明内容
本发明的目的就是为了解决上述问题，提出了一种基于异构情绪感染模型的人群疏散仿真方法。该方法通过对情绪感染现象的仿真，实现了情绪感染下的人群运动，能够更好地反映现实情况下个体情绪的异构性，更真实地模拟人群运动过程，为解决公共安全问题提供了新的解决思路。
为了实现上述目的，本发明采用如下技术方案：
基于异构情绪感染模型的人群疏散仿真方法，包括以下步骤：
步骤(1)：输入人群信息、情绪信息和场景信息，分别对人群信息、情绪信息、场景信息初始化；
步骤(2)：基于Durupinar模型计算异构情绪感染程度；
步骤(3)：将步骤(2)得到的结果与相对速度障碍法RVO相耦合；
步骤(4)：结合真实感渲染平台对人群运动仿真结果进行真实感渲染。
所述步骤(2)采用社会比较理论的相似性计算方法，将个体之间的特征进行比较，并计算两者之间的相似度。利用一组有序的特征模拟每个个体，并对所有特征值加权线性求和来度量个体i，j之间的相似性s(i,j)，个体之间的相似性越高，越容易聚集在一起。
考虑个体的两个特征：距离和关系。其中，个体之间的距离越近，越容易聚在一起，而对于关系亲密性高的个体，必然会比无关系或亲密性低的个体更容易聚在一起。
在基于分组的情绪感染计算的基础上，进行基于特征的情绪恐慌程度计算。
在现实生活中，性别会影响个体情绪感染的强度，女性相对于男性更容易受到情绪感染；情绪感染的个体差异也与个性因素有关，有些人在人群中易被主流情绪影响，而有些人不易被影响；个体年龄越大，生活阅历越丰富，人们能够更好地调节情绪反应，更多地关注积极情绪，因此较少受到负面刺激的影响。因此，在人群恐慌情绪传染过程中，考虑个性、性别、年龄三个特征的影响。
所述人群信息，包括：人群的位置；所述情绪信息，包括：人群的情绪状态；所述场景信息，包括：场景的障碍物位置。
场景信息的初始化是输入障碍物信息并画出场景的Roadmap图。所述Roadmap图表示场景语义信息，通过Roadmap图将无障碍区域表示成图结构Graph，图结构Graph表示的是绕过障碍物的拓扑图；Roadmap图的顶点是场景内的随机点，Roadmap图的边是连接顶点的边，将任意两个顶点连接，同时在连接时要保证两个顶点之间的连线不经过任何障碍物；
人群信息的初始化是综合考虑场景的障碍物信息，把个体随机初始化在场景范围内并避开所有的障碍物；
情绪信息的初始化则是依据个体初始化情绪值，把所有个体按情绪状态分为已感染情绪个体和未感染情绪个体。
步骤(2)中，将个体划分为两种类型：S表示易感染情绪状态，I表示已感染情绪状态。将S状态个体的情绪值e设定为0，I状态个体的情绪值e设定在(0，1]之间。定义个体i的感知半径R，以个体i为圆心，以半径R画圆，圆形区域为个体i的感知范围，在感知范围内的其他个体为个体i的邻居。在S状态个体的感知范围内，当处于S状态个体与I状态个体相接触时，S状态个体将接收来自I状态个体的恐慌情绪，一旦S状态个体的情绪值e超过个体情绪值设定阈值，则S状态个体将会变为I状态个体，按照个性、性别、年龄特征计算感染情绪个体的情绪值e，所述感染情绪个体的情绪值e用于对人群运动速度的计算，并根据情绪值e对个体的情绪恐慌程度进行划分。
对于每个个体定义{e,d,D,T}四个变量，其中e是个体的情绪值，其取值范围为e∈[0,1]；d是个体的恐慌量；D是个体的恐慌累积量，通过d的累加获得；T是个体的恐慌阈值。
所述步骤(2)，包括：基于分组的情绪感染计算个体与个体之间的情绪影响因子和基于特征的情绪恐慌程度计算个体的情绪恐慌程度两个步骤。
所述步骤(2)中基于分组的情绪感染计算个体与个体之间的情绪影响因子，步骤如下：
步骤(201)：从人群中选取前k个个体，作为k个簇的各自的中心；
步骤(202)：利用个体之间相似度计算公式计算剩余个体到k个簇心的相似度，将剩余个体划分到相似度最大的簇中；
步骤(203)：根据聚类结果，重新计算k个簇各自的虚拟中心Ov，
步骤(204)：计算簇内所有元素Pi到虚拟簇心Ov的相似度s(Ov，Pi)，其中，Pi∈clusterm,1≤m≤k，m表示群组序号。选取与虚拟簇心相似度最大的元素Pi，将最大的元素Pi作为新的簇心。
步骤(205)：重复步骤(202)、步骤(203)和步骤(204)，直到k个簇心不再发生改变；
步骤(206)：人群划分为k个群组，个体i的所在分组为g(i)，其中，g(i)∈G，G＝{g(1),g(2),g(3)...g(k)}，G表示分组的集合。对于任意个体i，j，计算个体i，j两者之间的情绪影响因子we(i,j)，组内与组间个体的情绪影响因子不同。
所述步骤(202)中：
个体i，j之间的相似度s(i,j)计算公式如下：
其中，wd表示距离相似性特征的权重，wr表示关系相似性特征的权重；1/dis(i,j)表示个体i，j之间的距离特征值，个体i，j之间的距离dis(i,j)利用欧氏距离计算，rel(i,j)表示个体i，j之间的关系特征值，rel(i,j)＝U[0,1]，rel(i,j)均匀分布随机生成，rel(i,j)数值越大代表个体之间的关系越亲密。
所述步骤(203)的计算方法是对于每个簇，根据公式(2)(3)(4)，来产生一个虚拟簇心Ov，并将虚拟簇心Ov更新为簇心。
其中，1≤j≤agentnum，agentnum表示所有个体的总数，num表示该簇的个体数量。agent[i].posdx表示个体i的横坐标，agent[i].posdy表示个体i的纵坐标，agent[i].rel[j]表示个体i与个体j之间的关系值。
所述步骤(206)：个体i，j两者之间的情绪影响因子we(i,j)的计算公式：
所述步骤(2)在基于特征的情绪恐慌程度计算个体的情绪恐慌程度之前的步骤：
利用OCEAN个性模型，根据个性值N的大小，将人群个体划分为三类：敏感型、稳定型和冷静型；
定义个性感知因子PN，个性感知因子PN在(0，1)之间均匀随机产生，通过定义per1，per2的值来调整人群中不同个性的比例；
不同个性的人产生的恐慌程度是不同的，因此，定义个性恐慌因子fp，表示个性对个体情绪恐慌程度的影响。
利用正态分布函数
生成fp，fp∈(0,1)，μp和
表示p的期望和方差；μp和σp的取值受个性感知因子PN影响，如表1所示：
表1个性恐慌因子划分
| P<sub>N</sub>的值 | 个性分类 | 个性特质 | μ<sub>p</sub>的值 | σ<sub>p</sub>的值 |
| [0,per<sub>1</sub>] | 敏感型 | 冲动、不安 | 0.9 | 0.2 |
| (per<sub>1</sub>,per<sub>2</sub>] | 稳定型 | 随和、较冷静 | 0.55 | 0.2 |
| (per<sub>2</sub>,1.0] | 冷静型 | 冷静，保守 | 0.25 | 0.2 |
定义性别感知因子gen，性别感知因子gen在(0，1)之间均匀随机产生，通过定义gen1的值来调整人群中的性别比例；
由于性别的差异，不同性别的人产生的恐慌有所不同，定义性别恐慌因子fg，表示性别对个体情绪恐慌程度的影响。
利用正态分布函数
生成fg，fg∈(0,1)，μg和
表示g的期望和方差；μg和σg的取值受gen影响，如表2所示：
表2性别恐慌因子划分
| gen的值 | 性别 | μ<sub>g</sub>的值 | σ<sub>g</sub>的值 |
| [0,gen<sub>1</sub>] | 男 | 0 | 0.5 |
| (gen<sub>1</sub>,1.0] | 女 | 1.4 | 0.5 |
定义年龄感知因子age，年龄感知因子age在(0，1)之间均匀随机产生，通过设定age1，age2的值来调整人群中的年龄比例；
不同年龄阶段的人产生的恐慌是有所不同的，定义年龄恐慌因子fa表示年龄对个体情绪恐慌程度的影响；
利用正态分布函数
生成fa，fa∈(0,1)，μa和
表示a的期望和方差，μa和σa的取值受age影响，如表3所示：
表3年龄恐慌因子划分
| age的值 | 年龄段 | μ<sub>a</sub>的值 | σ<sub>a</sub>的值 |
| [0,age<sub>1</sub>] | 青年 | 0.9 | 0.2 |
| (age<sub>1</sub>,age<sub>2</sub>] | 中年 | 0.55 | 0.2 |
| (age<sub>2</sub>,1.0] | 老年 | 0.25 | 0.2 |
。
步骤(2)中，基于特征的情绪恐慌程度计算个体的情绪恐慌程度，计算过程为：
步骤(211)：对所有个体的个性感知因子PN、性别感知因子gen和年龄感知因子age的值进行随机初始化；
步骤(212)：利用per1、per2的值，根据表1初始化个性恐慌因子fp；利用gen1的值，根据表2初始化性别恐慌因子fg；利用age1、age2的值，根据表3初始化年龄恐慌因子fa；
步骤(213)：将所有个体划分为两类：已感染情绪的个体和易感染情绪个体；
步骤(214)：对于易感染情绪个体i，令个体的情绪值e＝0，个体的恐慌量d＝0，个体的恐慌累积量D＝0，利用
和
并计算个体i的恐慌阈值Ti；
步骤(215)：对于已感染情绪的个体j，令个体的恐慌累积量D＝0，个体的情绪值e＝1.0；利用
和
并根据公式(6)计算Ti，利用
和
并计算个体i的恐慌量di；
步骤(216)：在每个时间步t，对易感染情绪个体i的情绪值ei进行更新；
步骤(217)：根据个体的情绪值ei，利用表4将个体情绪划分为情绪状态。
表4情绪恐慌因子划分
| 情绪状态 | e<sub>i</sub>的值 |
| 无恐慌 | 0 |
| 轻度恐慌 | (0，0.4] |
| 中度恐慌 | (0.4，0.7] |
| 重度恐慌 | (0.7，1.0] |
。
所述步骤(214)中计算个体i的恐慌阈值Ti的公式
其中，
和
分别表示期望和方差，
是一个对数正态分布函数，也就是以e为底的函数。Ti通过对数正态分布函数计算可得。对数分布函数是对数为正态分布的任意随机变量的概率分布。
首先确定Ti是在[x1,x2]之间的任意随机数，把随机数Ti代入概率密度函数：
获得对数正态分布函数的概率密度y，在[0,y]之间产生一个随机数y1，而Ti的取值必须满足y1>y；令
x1＝10，x2＝70。
所述步骤(215)计算个体i的恐慌量di的公式：
其中，
和
分别表示期望和方差，
是一个对数正态分布函数，di通过对数正态分布函数计算可得；令
x1＝2，x2＝7。
所述步骤(216)：
步骤(2161)：利用欧氏距离公式计算个体i，j之间的距离dis(i,j)，如果dis(i,j)<R，并且在彼此的感知范围内，则两个个体视彼此为邻居；
步骤(2162)：计算在感知范围内的第j个体对易感染情绪个体i的恐慌传播值dj，根据公式(8)更新个体i接收的恐慌累积量D(i,R)；
其中，Neighbor(i)表示个体i的邻居。
步骤(2163)：如果Di<Ti，则易感染情绪个体i未被感染情绪，ei＝0；如果Di≥Ti，Ti是个体i的恐慌阈值，说明易感染情绪个体i被感染情绪，根据公式(9)计算个体的情绪值ei：
ei＝wp×fp+wg×fg+wa×fa (9)
其中，wp表示个性特征的权重，wg表示性别特征的权重，wa表示年龄特征的权重，且wp+wg+wa＝1，通过调节各权重的大小，来调节特征的影响程度。
所述步骤(3)中情绪感染模型与人群运动耦合的步骤如下：
步骤(31)：初始化个体参数：vpref、
vmax和amax；其中，vpref表示期望速度；
表示个体的速度惩罚因子，
的值与个体的活跃度成反比；vmax表示个体的最大速度；amax表示个体的最大加速度；
步骤(32)：在每个时间步t内，计算个体下一时刻的新速度；
步骤(33)：重复步骤(32)，直到个体到达目标位置。
所述步骤(32)的步骤：
步骤(321)：利用vmax，amax，根据公式(10)来限制个体的速度，得到在下一时刻的候选速度集合AVi(Vi)。
AVi(Vi)＝{V′i|||V′i||<Vi max∩||V′i-Vi||<amaxΔt} (10)
其中，个体下一时刻的候选速度V′i受最大速度以及最大加速度的限制，||V′i||<Vi max表示候选速度值不能超过个体的最大速度值，||V′i-Vi||<amaxΔt是由速度公式vt+1＝vt+at得到，即个体的候选速度也受最大加速度的限制。AVi表示下一时刻的合理速度集合；Vi表示个体当前的速度；V′i表示个体下一时刻的候选速度；vt+1表示个体下一时刻的速度；vt表示个体的当前速度；a表示个体的加速度。
步骤(322)：利用情绪值ei根据公式(11)，更新个体i的期望速度Vi pref(ei)，同时保证Vi pref(ei)满足公式(12)的要求：
Vi pref(ei)＝(1+ei)Vi pref (11)
||Vi pref(ei)||≤Vi max (12)
步骤(323)：根据公式(13)和(14)，利用期望速度Vi pref，在候选速度集合AVi(Vi)中选择使速度惩罚值penaltyi(V′i)最小的候选速度，作为下一时刻个体最佳的速度Vt+1。
其中，penaltyi(V′i)表示个体的速度惩罚值；
是个体的速度惩罚因子。tc′i(V′i)是个体i与周围个体的期望碰撞时间，||Vi pref-V′i||表示期望速度Vi pref与候选速度V′i的差值。个体的下一时刻速度Vt+1需要满足与其他个体碰撞时间最大，与期望速度的误差值最小两个条件。而情绪值ei通过改变期望速度Vi pref来影响个体的运动。
所述步骤(4)中的仿真系统是基于XNA技术开发的跨平台仿真系统，三维实时真实感渲染平台以Microsoft.NET Framework 4.0、XNA 4.0为平台。
基于异构情绪感染模型的人群疏散仿真系统，包括：存储器、处理器以及存储在存储器上并在处理器上执行的计算机指令，所述计算机指令在处理器上运行时，完成以下步骤：
步骤(1)：输入人群信息、情绪信息和场景信息，分别对人群信息、情绪信息、场景信息初始化；
步骤(2)：基于Durupinar模型计算异构情绪感染程度；
步骤(3)：将步骤(2)得到的结果与相对速度障碍法RVO相耦合；
步骤(4)：结合真实感渲染平台对人群运动仿真结果进行真实感渲染。
一种计算机可读存储介质，其上存储有计算机指令，所述计算机指令在处理器上运行时，完成以下步骤：
步骤(1)：输入人群信息、情绪信息和场景信息，分别对人群信息、情绪信息、场景信息初始化；
步骤(2)：基于Durupinar模型计算异构情绪感染程度；
步骤(3)：将步骤(2)得到的结果与相对速度障碍法RVO相耦合；
步骤(4)：结合真实感渲染平台对人群运动仿真结果进行真实感渲染。
本发明的有益效果：本发明以社会比较理论为基础，提出了一种基于异构情绪感染模型的人群疏散仿真方法。为了实现该方法，本发明构建了基于Durupinar模型的异构情绪感染方法的人群运动框架：首先，在运动过程中，将K-Mediods聚类算法与社会比较理论结合，根据距离和关系对人群进行分组，利用分组结果研究组间和组内不同的情绪感染计算方法；其次，为了有效地体现情绪的个体差异性，根据个性、性别、年龄的差异将个体特殊化，计算个体之间不同的情绪恐慌程度。然后，将上述方法与相对速度障碍法耦合，利用情绪驱动人群运动。最后，使用真实感渲染的方法模拟人群疏散效果。该方法能够更加真实地模拟情绪感染下的人群运动，体现了个体情绪的异构性，提高了人群疏散效率，为解决公共安全问题提供了新的解决思路。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。
图1为本发明模型简化的流程图。
具体实施方式
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
如图1所示，基于异构情绪感染模型的人群疏散仿真方法，包括以下几个过程：
过程1：初始化。这里的初始化包括对场景信息，个体的位置信息和情绪的初始化。
首先对场景障碍物进行初始化，画出场景的roadmap图。本发明使用Roadmap表示场景语义信息，通过Roadmap将无障碍区域表示成图(Graph)结构，表示的是绕过障碍物的拓扑图。该拓扑图由顶点和边组成G＝(V,E)。其中V表示顶点的集合，E表示边的集合，每条边e(v1,v2)∈E，其中顶点v1和顶点v2是边e的两个端点。Roadmap的边是连接顶点的，将任意两个距离为r的顶点连接，同时在连接时要保证两个顶点之间的连线不经过任何障碍物。
个体位置的初始化是综合考虑场景信息和障碍物信息，把个体随机初始化在场景范围内并避开所有的障碍物。情绪的初始化则是对个体初始化情绪值，把所有个体的情绪分为已感染和易感两种状态。
过程1：模拟基于分组的情绪感染计算过程。
第一步：从人群中选取前k个个体，作为k个簇的各自的中心，其中k表示分组个数；；
第二步：利用公式(1)计算剩余个体到k个簇心的相似度，将这些个体划分到相似度最大的簇中；
其中，wd表示距离相似性特征的权重(非负)，wr表示关系相似性特征的权重(非负)；1/dis(i,j)表示个体i，j之间的距离特征值，rel(i,j)表示个体i，j之间的关系特征值。
第三步：根据聚类结果，重新计算k个簇各自的中心，计算方法是对于每个簇，根据公式(2)(3)(4)，来产生一个虚拟簇心Ov，并将虚拟簇心Ov更新为簇心。
其中，1≤j≤agentnum，agentnum表示所有个体的总数，num表示该簇的个体数量，agent[i].posdx表示个体i的横坐标，agent[i].posdy表示个体i的纵坐标，agent[i].rel[j]表示个体i与个体j之间的关系值。
第四步：利用公式(1)计算簇内所有元素到虚拟簇心的相似度s(Ov,Pi)，其中，Pi∈clusterm,1≤m≤k，m表示群组序号。选取与虚拟簇心相似度最大的元素Pi，将其作为新的簇心。
第五步：重复2、3、4步，直到k个簇心不再发生改变；
第六步：将人群划分为k个群组，个体i的所在分组为g(i)，其中g(i)∈G,G＝{g1,g2,g3...gk}。对于任意个体i，j，利用公式(5)计算两者之间的情绪影响因子we，即组内与组间个体的情绪影响因子不同。
过程2：模拟基于特征的情绪恐慌程度计算过程。
第一步：对所有个体的PN、gen、age的值进行随机初始化。
第二步：在一个场景中，不同内在状态的个体具有不同的情绪恐慌程度。本发明定义个体性别影响因子fg，年龄影响因子fa，个性影响因子fp，其取值设置在(0，1)内。结合实际场景，调整per1，per2的值控制不同个性的个体数量，调整age1，age2的值控制不同年龄段的个体数量，调整gen1的值控制不同性别的个体数量，并利用正态分布函数对其初始化。
第三步：对未感染情绪个体进行初始化，其中情绪值e＝0，恐慌量d＝0，恐慌累积量D＝0，阈值T可通过μT和
根据公式(6)得到。
其中，μT和
分别表示T的期望和方差；
是一个对数正态分布函数。
第四步：对已感染情绪的个体进行初始化，其中情绪值e＝1.0，恐慌累积量D＝0，阈值T可通过μT和
根据公式(6)得到，恐慌量d可通过公式(7)计算得到。
其中，μd和
分别表示d的期望和方差；
是一个对数正态分布函数。
第五步：对于易感个体，计算个体的邻居。邻居的定义是以个体i为圆心画圆，同时定义个体的可视域，所有在圆内并且满足可视域范围的都是个体i的邻居。在实验中，个体与其邻居的距离通过半径R进行控制。根据现实生活中个体的影响范围，我们设置R<＝40。在本次实验中，我们设置R＝30。
在运动过程中，利用欧氏距离公式计算个体j与i之间的距离dis，若dis<R，并且j在i的可视域范围内，则j即为i的邻居。个体接收感知范围内的邻居的恐慌量，并根据公式(8)更新恐慌累积量D(i,R)。
在公式(8)中，Neighbor(i)表示个体i的邻居，D(i,R)表示个体i接收邻居恐慌量的累加值。邻居包括已感染个体和未感染个体，其中未感染个体的恐慌量d为0，已感染个体的恐慌量d通过对数正态分布函数计算得到。如果D≥T，说明个体被感染情绪，情绪值得到更新。这里的T是个体的阈值，通过对数正态分布函数得到。
第六步：根据第三步计算得到的阈值T和第五步计算得到的恐慌累加量D，如果D>＝T，说明个体被感染情绪，根据公式(9)计算个体的情绪值。
ei＝wp×fp+wg×fg+wa×fa (9)
其中，wp，wg，wa分别表示个性、性别、年龄特征的权重，且wp+wg+wa＝1，通过调节各权重的大小，来调节特征的影响程度。
第七步：如果D<T，则ei＝0，个体继续进行第五步和第六步的操作，直到个体到达目标位置。
第八步：对于已感染情绪个体，在每个时间步t内其情绪值e不再变化。
过程4：人群运动计算，在这一过程中，将情绪感染加入人群运动计算，控制整个人群的运动速度。个体的运动方向为出口所在方向。运动通过路径图进行导航，每一个图的顶点可看作一个临时目标点，个体的运动方向即向着临时目标点运动。
第一步：初始化个体参数：vpref、
vmax和amax。
其中，vpref是个体的期望速度，vmax和amax分别是是个体的最大速度和个体的最大加速度，
表示个体的速度惩罚因子，
的值与个体的活跃度成反比。
第二步：利用最大速度vmax，最大加速度amax，根据运动学和动力学公式(10)来限制个体的速度，得到在下一时刻的合理速度集合AV。
AVi(Vi)＝{V′i|||V′i||<Vi max∩||V′i-Vi||<amaxΔt} (10)
其中，Vi表示个体i的当前速度，V′i表示个体i下一时刻的候选速度，AVi表示个体的所有候选速度集合，所有符合条件的速度都属于AVi。||V′i-Vi||表示候选速度相对于当前速度的变化值的大小，t表示时间步，amaxt表示单位时间内最大速度变化值。在公式中，候选速度的大小受最大速度Vmax和最大加速度amax的制约。因此，||V′i||应小于最大速度值Vmax，且其速度变化值||V′i-Vi||也应小于amaxt。
第三步：利用情绪值e根据公式(11)，更新个体i的期望速度vpref，同时保证vpref满足公式(12)的要求；个体的期望速度值vpref，它是一个矢量，它的大小等于期望速度在目标方向上的值。个体的期望速度和情绪有关，个体情绪值越高，期望速度值越大。但是期望速度值并不是无限大，它始终小于最大速度值。
Vi pref(ei)＝(1+ei)Vi pref (11)
||Vi pref(ei)||≤Vi max (12)
第四步：根据公式(13)和(14)，利用期望速度vpref，在合理速度集合AV中选择使速度惩罚值penalty最小的速度，作为下一时刻个体最佳的速度Vt+1。
其中，penaltyi是个体i的惩罚度量，tc′i(V′i)是个体i与周围个体的期望碰撞时间，
是个体的速度惩罚因子，
的值与个体的活跃度成反比。个体的penaltyi是由期望速度与候选速度的绝对偏差以及期望避障时间共同决定。个体的最佳运动速度就是使得penaltyi最小的候选速度。
第五步：如果个体没有到达目标，则继续更新情绪和速度，否则个体停止运动。
过程5：结合真实感渲染技术，将人群运动计算结果进行渲染得到人群疏散仿真动画。
具体地，在真实感渲染平台中导入人群运动路径，生成人群运动的动画效果。该真实感渲染平台是基于XNA/MonoGame平台实现的支持跨平台的实时人群运动仿真系统。采用XNA Game Studio2013，以Microsoft Visual Studio2013为平台进行人群运动仿真。已生成的路径文件导入到该仿真平台，可以更加直观的观察人群运动效果，分析人群运动的情况。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (9)
1.基于异构情绪感染模型的人群疏散仿真方法，其特征是，包括以下步骤：
步骤(1)：输入人群信息、情绪信息和场景信息，分别对人群信息、情绪信息、场景信息初始化；
步骤(2)：基于Durupinar模型计算异构情绪感染程度；
步骤(3)：将步骤(2)得到的结果与相对速度障碍法RVO相耦合；
步骤(4)：结合真实感渲染平台对人群运动仿真结果进行真实感渲染；
所述步骤(2)，包括：基于分组的情绪感染计算个体与个体之间的情绪影响因子和基于特征的情绪恐慌程度计算个体的情绪恐慌程度两个步骤；
所述基于分组的情绪感染计算个体与个体之间的情绪影响因子包括：
将人群划分为k个群组，个体i的所在分组为g(i)，其中g(i)∈G,G＝{g1,g2,g3...gk}；对于任意个体i，j，计算两者之间的情绪影响因子we：
所述基于特征的情绪恐慌程度计算个体的情绪恐慌程度包括：
对于易感个体，计算个体的邻居，邻居的定义是以个体i为圆心画圆，同时定义个体的可视域，所有在圆内并且满足可视域范围的都是个体i的邻居，个体与其邻居的距离通过半径R进行控制；
在运动过程中，利用欧氏距离公式计算个体j与i之间的距离dis，若dis<R，并且j在i的可视域范围内，则j即为i的邻居；个体接收感知范围内的邻居的恐慌量，更新恐慌累积量D(i,R)：
其中，Neighbor(i)表示个体i的邻居，D(i,R)表示个体i接收邻居恐慌量的累加值；邻居包括已感染个体和未感染个体，其中未感染个体的恐慌量d为0，已感染个体的恐慌量d通过对数正态分布函数计算得到；如果D≥T，说明个体被感染情绪，情绪值得到更新；T是个体的阈值，通过对数正态分布函数得到；
根据阈值T和恐慌累加量D，如果D>＝T，说明个体被感染情绪，计算个体的情绪值；
ei＝wp×fp+wg×fg+wa×fa
其中，wp，wg，wa分别表示个性、性别、年龄特征的权重，且wp+wg+wa＝1，通过调节各权重的大小，来调节特征的影响程度；其中fp表示个性对个体情绪恐慌程度的影响，fg为性别对个体情绪恐慌程度的影响，fa为年龄对个体情绪恐慌程度的影响；
如果D<T，则ei＝0，个体继续进行更新恐慌累积量和计算个体的情绪值的操作，直到个体到达目标位置；
对于已感染情绪个体，在每个时间步t内其情绪值e不再变化。
2.如权利要求1所述的基于异构情绪感染模型的人群疏散仿真方法，其特征是，所述人群信息，包括：人群的位置；所述情绪信息，包括：人群的情绪状态；所述场景信息，包括：场景的障碍物位置；
场景信息的初始化是输入障碍物信息并画出场景的Roadmap图；所述Roadmap图表示场景语义信息，通过Roadmap图将无障碍区域表示成图结构Graph，图结构Graph表示的是绕过障碍物的拓扑图；Roadmap图的顶点是场景内的随机点，Roadmap图的边是连接顶点的边，将任意两个顶点连接，同时在连接时要保证两个顶点之间的连线不经过任何障碍物；
人群信息的初始化是综合考虑场景的障碍物信息，把个体随机初始化在场景范围内并避开所有的障碍物；
情绪信息的初始化则是依据个体初始化情绪值，把所有个体按情绪状态分为已感染情绪个体和未感染情绪个体。
3.如权利要求1所述的基于异构情绪感染模型的人群疏散仿真方法，其特征是，步骤(2)中，将个体划分为两种类型：S表示易感染情绪状态，I表示已感染情绪状态；将S状态个体的情绪值e设定为0，I状态个体的情绪值e设定在(0，1]之间；定义个体i的感知半径R，以个体i为圆心，以半径R画圆，圆形区域为个体i的感知范围，在感知范围内的其他个体为个体i的邻居；在S状态个体的感知范围内，当处于S状态个体与I状态个体相接触时，S状态个体将接收来自I状态个体的恐慌情绪，一旦S状态个体的情绪值e超过个体情绪值设定阈值，则S状态个体将会变为I状态个体，按照个性、性别、年龄特征计算感染情绪个体的情绪值e，所述感染情绪个体的情绪值e用于对人群运动速度的计算，并根据情绪值e对个体的情绪恐慌程度进行划分；
对于每个个体定义{e,d,D,T}四个变量，其中e是个体的情绪值，其取值范围为e∈[0,1]；d是个体的恐慌量；D是个体的恐慌累积量，通过d的累加获得；T是个体的恐慌阈值。
4.如权利要求3所述的基于异构情绪感染模型的人群疏散仿真方法，其特征是，所述步骤(2)中基于分组的情绪感染计算个体与个体之间的情绪影响因子，步骤如下：
步骤(201)：从人群中选取前k个个体，作为k个簇的各自的中心；
步骤(202)：利用个体之间相似度计算公式计算剩余个体到k个簇心的相似度，将剩余个体划分到相似度最大的簇中；
步骤(203)：根据聚类结果，重新计算k个簇各自的虚拟中心Ov，
步骤(204)：计算簇内所有元素Pi到虚拟簇心Ov的相似度s(Ov，Pi)，其中，Pi∈clusterm,1≤m≤k，m表示群组序号；选取与虚拟簇心相似度最大的元素Pi，将最大的元素Pi作为新的簇心；
步骤(205)：重复步骤(202)、步骤(203)和步骤(204)，直到k个簇心不再发生改变；
步骤(206)：人群划分为k个群组，个体i的所在分组为g(i)，其中，g(i)∈G，G＝{g(1),g(2),g(3)...g(k)}，G表示分组的集合；对于任意个体i，j，计算个体i，j两者之间的情绪影响因子we(i,j)，组内与组间个体的情绪影响因子不同。
5.如权利要求1所述的基于异构情绪感染模型的人群疏散仿真方法，其特征是，所述步骤(2)在基于特征的情绪恐慌程度计算个体的情绪恐慌程度之前的步骤，还包括：
利用OCEAN个性模型，根据个性值N的大小，将人群个体划分为三类：敏感型、稳定型和冷静型；
定义个性感知因子PN，个性感知因子PN在(0，1)之间均匀随机产生，通过定义per1，per2的值来调整人群中不同个性的比例；
不同个性的人产生的恐慌程度是不同的，因此，定义个性恐慌因子fp，表示个性对个体情绪恐慌程度的影响；
利用正态分布函数
生成fp，fp∈(0,1)，μp和
表示p的期望和方差；μp和σp的取值受个性感知因子PN影响，如表1所示：
表1个性恐慌因子划分
定义性别感知因子gen，性别感知因子gen在(0，1)之间均匀随机产生，通过定义gen1的值来调整人群中的性别比例；
由于性别的差异，不同性别的人产生的恐慌有所不同，定义性别恐慌因子fg，表示性别对个体情绪恐慌程度的影响；
利用正态分布函数
生成fg，fg∈(0,1)，μg和
表示g的期望和方差；μg和σg的取值受gen影响，如表2所示：
表2性别恐慌因子划分
定义年龄感知因子age，年龄感知因子age在(0，1)之间均匀随机产生，通过设定age1，age2的值来调整人群中的年龄比例；
不同年龄阶段的人产生的恐慌是有所不同的，定义年龄恐慌因子fa表示年龄对个体情绪恐慌程度的影响；
利用正态分布函数
生成fa，fa∈(0,1)，μa和
表示a的期望和方差，μa和σa的取值受age影响，如表3所示：
表3年龄恐慌因子划分
。
6.如权利要求5所述的基于异构情绪感染模型的人群疏散仿真方法，其特征是，步骤(2)中，基于特征的情绪恐慌程度计算个体的情绪恐慌程度，计算过程为：
步骤(211)：对所有个体的个性感知因子PN、性别感知因子gen和年龄感知因子age的值进行随机初始化；
步骤(212)：利用per1、per2的值，根据表1初始化个性恐慌因子fp；利用gen1的值，根据表2初始化性别恐慌因子fg；利用age1、age2的值，根据表3初始化年龄恐慌因子fa；
步骤(213)：将所有个体划分为两类：已感染情绪的个体和易感染情绪个体；
步骤(214)：对于易感染情绪个体i，令个体的情绪值e＝0，个体的恐慌量d＝0，个体的恐慌累积量D＝0，利用
和
并计算个体i的恐慌阈值Ti；
步骤(215)：对于已感染情绪的个体j，令个体的恐慌累积量D＝0，个体的情绪值e＝1.0；利用
和
并根据公式(6)计算Ti，利用
和
并计算个体i的恐慌量di；
其中，
和
分别表示期望和方差，
是一个对数正态分布函数，也就是以e为底的函数；步骤(216)：在每个时间步t，对易感染情绪个体i的情绪值ei进行更新；
步骤(217)：根据个体的情绪值ei，利用表4将个体情绪划分为情绪状态；
表4情绪恐慌因子划分
。
7.如权利要求1所述的基于异构情绪感染模型的人群疏散仿真方法，其特征是，所述步骤(3)中情绪感染模型与人群运动耦合的步骤如下：
步骤(31)：初始化个体参数：vpref、
vmax和amax；其中，vpref表示期望速度；
表示个体的速度惩罚因子，
的值与个体的活跃度成反比；vmax表示个体的最大速度；amax表示个体的最大加速度；
步骤(32)：在每个时间步t内，计算个体下一时刻的新速度；
步骤(33)：重复步骤(32)，直到个体到达目标位置。
8.基于异构情绪感染模型的人群疏散仿真系统，其特征是，包括：存储器、处理器以及存储在存储器上并在处理器上执行的计算机指令，所述计算机指令在处理器上运行时，完成如权利要求1-7任一项所述的基于异构情绪感染模型的人群疏散仿真方法。
9.一种计算机可读存储介质，其特征是，其上存储有计算机指令，所述计算机指令在处理器上运行时，完成如权利要求1-7任一项所述的基于异构情绪感染模型的人群疏散仿真方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710940633.4A CN107665282B (zh) | 2017-10-11 | 2017-10-11 | 基于异构情绪感染模型的人群疏散仿真方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710940633.4A CN107665282B (zh) | 2017-10-11 | 2017-10-11 | 基于异构情绪感染模型的人群疏散仿真方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107665282A CN107665282A (zh) | 2018-02-06 |
| CN107665282B true CN107665282B (zh) | 2021-06-25 |
Family
ID=61098122
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710940633.4A Expired - Fee Related CN107665282B (zh) | 2017-10-11 | 2017-10-11 | 基于异构情绪感染模型的人群疏散仿真方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107665282B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109408889B (zh) * | 2018-09-21 | 2022-08-12 | 同济大学 | 基于信息熵的宏观人群恐慌度量方法及其应用 |
| CN109460591A (zh) * | 2018-10-24 | 2019-03-12 | 山东师范大学 | 基于优化积极情绪感染的人群疏散仿真方法及装置 |
| CN109740192B (zh) * | 2018-12-12 | 2023-04-07 | 山东师范大学 | 基于阿诺德情绪模型的人群疏散仿真方法及系统 |
| CN111144006B (zh) * | 2019-12-27 | 2023-04-28 | 山东师范大学 | 基于知识和情绪双重驱动的人群疏散仿真方法及系统 |
| CN111339668B (zh) * | 2020-02-28 | 2022-05-10 | 西南交通大学 | 基于情绪认知的人群疏散可视化方法 |
| CN111639809B (zh) * | 2020-05-29 | 2023-07-07 | 华中科技大学 | 基于领导者和恐慌情绪的多智能体疏散仿真方法和系统 |
| CN112347626B (zh) * | 2020-10-28 | 2022-10-11 | 山东师范大学 | 一种人群疏散中恐慌情绪的最优化干预仿真方法及系统 |
| CN113177347A (zh) * | 2021-05-10 | 2021-07-27 | 清华大学 | 一种人群运动模拟方法、装置、设备及存储介质 |
| CN113239566B (zh) * | 2021-05-27 | 2021-11-09 | 大连理工大学 | 基于智能体的融合情绪化的灾后三级层次仿真系统及方法 |
| CN117236147B (zh) * | 2023-07-20 | 2024-07-02 | 中国地质大学（武汉） | 一种高密度人群疏散仿真方法、装置及存储介质 |
| CN117994767B (zh) * | 2023-12-21 | 2024-09-10 | 中南大学 | 一种基于视觉信号的高原列车乘员疏散引导智能调控系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8655939B2 (en) * | 2007-01-05 | 2014-02-18 | Digital Doors, Inc. | Electromagnetic pulse (EMP) hardened information infrastructure with extractor, cloud dispersal, secure storage, content analysis and classification and method therefor |
| CN106485327A (zh) * | 2016-09-08 | 2017-03-08 | 同济大学 | 一种灾害条件下人群踩踏行为演化的知识发现方法 |
| CN106530375A (zh) * | 2016-09-28 | 2017-03-22 | 山东师范大学 | 一种个性化情绪感染的人群动画生成方法 |
| CN106909730A (zh) * | 2017-02-21 | 2017-06-30 | 山东师范大学 | 基于同伦映射算法的建筑物三维模型仿真方法及系统 |
-
2017
- 2017-10-11 CN CN201710940633.4A patent/CN107665282B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8655939B2 (en) * | 2007-01-05 | 2014-02-18 | Digital Doors, Inc. | Electromagnetic pulse (EMP) hardened information infrastructure with extractor, cloud dispersal, secure storage, content analysis and classification and method therefor |
| CN106485327A (zh) * | 2016-09-08 | 2017-03-08 | 同济大学 | 一种灾害条件下人群踩踏行为演化的知识发现方法 |
| CN106530375A (zh) * | 2016-09-28 | 2017-03-22 | 山东师范大学 | 一种个性化情绪感染的人群动画生成方法 |
| CN106909730A (zh) * | 2017-02-21 | 2017-06-30 | 山东师范大学 | 基于同伦映射算法的建筑物三维模型仿真方法及系统 |
Non-Patent Citations (3)
| Title |
|---|
| 基于个性化情绪感染的人群动画生成方法;曹梦晓 等;《计算机科学》;20170630;第44卷(第6期);第306-311，316页 * |
| 曹梦晓 等.基于个性化情绪感染的人群动画生成方法.《计算机科学》.2017,第44卷(第6期), * |
| 融入关系分组的人群运动仿真;柳广鹏等;《小型微型计算机系统》;20160831;第37卷(第8期);第1735-1740页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107665282A (zh) | 2018-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107665282B (zh) | 基于异构情绪感染模型的人群疏散仿真方法及系统 | |
| CN111639809B (zh) | 基于领导者和恐慌情绪的多智能体疏散仿真方法和系统 | |
| CN113762530B (zh) | 面向隐私保护的精度反馈联邦学习方法 | |
| CN113536996B (zh) | 一种基于大量真实人群运动视频的人群运动仿真方法 | |
| CN107993237A (zh) | 一种基于窄带约束的几何活动轮廓模型图像局部分割方法 | |
| Hong et al. | Fast volume preservation for a mass-spring system | |
| CN110262855B (zh) | 车联网中基于背景信息的成员推测攻击原型系统 | |
| CN107704667B (zh) | 模拟集群性的人群运动仿真方法、装置和系统 | |
| Xue et al. | Crowd queuing simulation with an improved emotional contagion model | |
| CN107423480A (zh) | 一种人群运动仿真方法及系统 | |
| Wong et al. | Guidance path scheduling using particle swarm optimization in crowd simulation | |
| CN109740192B (zh) | 基于阿诺德情绪模型的人群疏散仿真方法及系统 | |
| CN114818360B (zh) | 人群应急疏散场景下的疏散出口设置方法及装置 | |
| CN106530375A (zh) | 一种个性化情绪感染的人群动画生成方法 | |
| CN110956684A (zh) | 基于残差网络的人群运动疏散仿真方法及系统 | |
| CN112257258B (zh) | 基于Boids集群运动模型的人群运动仿真方法 | |
| CN104217100B (zh) | 一种基于agent社会圈子网络的语言竞争模型的建模仿真方法 | |
| CN111125911B (zh) | 一种基于分区双层网络的语言竞争微观仿真模型的建模仿真方法 | |
| Xiang et al. | Using SIR model to simulate emotion contagion in dynamic crowd aggregation process | |
| Lee et al. | Online indoor fire evacuation system | |
| CN114548497B (zh) | 一种实现场景自适应的人群运动路径规划方法及系统 | |
| Baig et al. | Realistic modeling of agents in crowd simulations | |
| Yu et al. | Understanding and Modeling Collision Avoidance Behavior for Realistic Crowd Simulation | |
| JP4814846B2 (ja) | 対象物衝突模擬装置、その方法、そのプログラム、及びその記録媒体 | |
| CN112560326B (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210625 Termination date: 20211011 |