Evacuation Model of Emotional Contagion Crowd Based on Cellular Automata

Crowd evacuation under emergency is an important task of world public security research and practice. In order to describe the microemotional contagion of evacuation individuals, a cellular automata-based evacuation model of emotional contagion crowd based on the classical SIS model of infectious diseases is proposed in this paper. Firstly, the state of evacuation individual is deﬁned as “emotional susceptible” and “emotional infective.” Then, a dynamic model considering emotional contagion is established with cellular automata. Based on the models of static ﬂoor ﬁeld and dynamic ﬂoor ﬁeld, the emotion updating rules and state updating rules are constructed. The inﬂuence of perception domain radius on pedestrian evacuation process is analyzed through experiments. The conclusion can provide evacuation guidance for evacuation individuals. The comparative experiment results show that the improved model can reﬂect the movement characteristics of evacuation individuals eﬀectively. The evacuation eﬃciency of the whole system is also eﬀectively improved due to the consideration of emotional contagion and evacuation strategy.


Introduction
With the development of social progress and frequent economic, political, and cultural exchanges, people's activities are increasingly rich. In recent years, various accidents caused by the crowd have happened frequently. erefore, how to grasp the characteristics of crowd evacuation behavior in emergency is the key to formulate effective guidance plan. In the past few decades, many researchers have mainly used experimental and simulation methods to study crowd evacuation behavior. e experimental methods mainly include the questionnaire method after the real event and controlled evacuation experiment with personnel participation. Afterwards, questionnaire method has information lag. Controllable experiments are often costly, and it is difficult to simulate dangerous scenarios, such as fire and terrorist attacks. erefore, building a crowd evacuation simulation model in line with the actual situation has gradually become an important method to analyze the characteristics of crowd evacuation behavior, which has been widely used [1][2][3]. Simulation model method is to design reasonable building facilities and environment and reshape the basic attributes of the crowd. en, take the crowd or individual as the research object to carry out the simulation of normal or emergency evacuation process. e evacuation simulation model can be used to evaluate the effectiveness of guidance measures, help the site management to make emergency plans, and provide theoretical support for emergency decision-making.
At present, the research of indoor crowd evacuation simulation model is mainly divided into macromodel and micromodel. e macrosimulation model usually adopts the classical graph theory method. e evacuation decisionmaking problem is transformed into an evacuation route selection problem. It usually adopts the minimum cost maximum flow method, the queuing theory method, and the shortest path method [4][5][6][7]. Micromodels usually take individuals as objects. It can identify many typical characteristics of individual movement or human evacuation behavior. Common models include cellular automata model, social force model, lattice gas model, and visual obstacle model [8][9][10][11][12][13]. Compared with other models, cellular automata (CA) model has powerful complex computing function and highly dynamic characteristics. It has a strong ability to simulate the spatiotemporal dynamic evolution of space complex system and requires relatively less computation. erefore, it has been widely used in crowd evacuation research. Cellular automata were proposed by Blue et al. [14] in 1997, and then it began to be applied in pedestrian evacuation. e basic CA model represents the environment site as a grid composed of cells. e evacuation individuals move from one cell to the adjacent cell according to certain rules. Liu et al. [15] established an improved CA evacuation model considering the influence of the density of people around the exit on the evacuation behavior. Li et al. [16] proposed an extended CA model, considering the effect of sudden hazard. Zhang et al. [17] established a multiagent decision-making framework based on CA model to study the evacuation problem of stadium with or without obstacles. Yang et al. [18] present a new static floor field method based on cellular automaton for simulations of evacuation processes. Zheng et al. [19] proposed an extended CA model to study the effect of three-dimensional smoke diffusion on crowd evacuation. Geng et al. [20] proposed an extended cellular automata model to simulate the individual evacuation information under adverse line of sight conditions with uncertainty.
In emergencies, individual emotion plays an important role in evacuation decision-making. e emotion of the evacuation individual is dynamic and can be perceived and spread. Franovetter [21] proposed the concept of "threshold theory" and proposed that every event has an emotional threshold. e OCEAN model [22] proposed by McCrae et al. is the most widely used personality model, and the model describes individual personality from five vector dimensions. Helbing [23] simulated the dynamic escape behavior of evacuating individuals and concluded that escape behavior was a behavioral manifestation of emotion and confirmed that emotion among many people was contagious. Wang et al. [24] investigated the correlations between evacuation crowd characteristics and the evacuation emotion in subway emergencies by questionnaire. Shi et al. [25] proposed a new deterministic model of impatience, which described the impatience during pedestrian evacuation qualitatively and quantitatively. Wang et al. [26] considered that emotional contagion has an impact on the behavior of both individuals and groups and formulated the control problem of emotional contagion as an optimal control problem of crowd evacuation. e human study of group emotion communication originated from infectious diseases. ere are two kinds of scientific epidemic models to explore how diseases spread among groups, such as SIR model and SIS model [27,28]. SIR model is to study the transmission relationship among "susceptible," "infective," and "removal." Fu et al. [29] modified the epidemiological SIR model and proposed a new (CA-SIRS) emotional contagion model. e SIS model considers that it is difficult for individuals to achieve complete immunity when considering the infection of some diseases. erefore, there are only "susceptible persons" and "infected persons" in the process of infection. e SIS model describes that an individual can change from a susceptible state to an infected state. After cure, he is still in a susceptible state and has the probability of being infected. However, the total number of infected people in SIS model should remain unchanged, and it is not suitable to describe the microemotional contagion of pedestrian evacuation individuals. erefore, this paper combines SIS model with CA model to build a crowd evacuation model, considering emotional contagion under emergency. en, the evacuation process is simulated.

Basic Theory
2.1. Cellular Automata Model. Cellular automata are a kind of dynamic system with discrete time and space. It uses discrete space layout and discrete time interval to divide the cell into finite states. e evolution of the individual state of the cell is only related to its current state and the state of a local neighborhood. Each cell scattered in the regular grid takes a finite discrete state, follows the same action rules, and updates synchronously according to the determined local rules. e basic elements of cellular automata are as follows: (1) Space: the spatial lattice of cells in space can be onedimensional, two-dimensional, or multidimensional. (2) State set: there are only two different states, which can be coded as 0 and 1, respectively, or represented by different colors. (3) Neighbor: there is a cell with radius of 1 around a cell, which can affect the state of the cell at the next moment. e most common neighbors are Moore type and von Neumann type. (4) According to the state of the cell and its eight neighboring cells, the dynamic function or state transfer equation of the cell state at the next moment is determined.

Typical Models of Infectious
Diseases. Emotional contagion can be studied with reference to the disease transmission model. ere are two kinds of traditional infectious disease models: SIR model and SIS model. SIR model refers to three types of individual status: "susceptible," "infective," and "removal." S, I, and R, respectively, represent the proportion of the three groups in the total group. SIR model does not calculate the new birth rate and death rate in the process of infection, so S + I + R � 1. e model reflects the macroprocess of human infection: when infectious diseases occur, the total number of three groups of people is constant. Once "susceptible" (S) contacted with "infective" (I) will be infected with probability; "infective" (I) will change to "removal" (R) with a certain probability; "removal" is equivalent to vaccine injector and will no longer be infected with disease immunity.
e other classical SIS model considers that when some diseases are infected, it is difficult for individuals to achieve complete immunity. erefore, there are only "susceptible" and "infective" in the process of infection. SIS emotional contagion model describes that an individual changes from an easily infected state to an infected state. After being cured, he is still an easily infected individual and will enter the infectious system again. e differential system dynamics equation of the model is shown in equation (1). e state transition diagram is shown in Figure 1.
where β is the daily contact infection rate; c is the daily cure rate. Since the SIS model adopts the macroscopic model of infection dynamics expressed by differential equation, it is not suitable for describing the microemotional contagion of evacuation individuals. In addition, the SIS model assumes that the total number of people does not change during the process of infection. However, when pedestrians evacuate, once they leave the exit, it is a safe area. at is, the number of people in the dangerous area at every moment is dynamic. Based on the above, this paper combines the classical SIS model with cellular automata to build a dynamic pedestrian evacuation model, considering emotional contagion in emergencies.

Evacuation Model of Emotional Contagion
Crowd Based on Cellular Automata 3.1. Basic Cell Definition of Cellular Automata. In this paper, cellular automata model is used to simulate the evacuation situation in independent space. e space is L in length and W in width. It is divided into cell lattices of 0.4 × 0.4 for one person to stand. Each cell is represented by a set of number pairs, such as cell(x, y). e distance between adjacent cells is cl � 0.4m. e number of exits is m. e width is el i , where i ∈ [1, 2, . . . , m]. When the initial state t � 0, pedestrians are randomly distributed in the space, each occupying a cell. After the evacuation, Moore neighborhood rule with radius of 1 is adopted in the neighborhood of the cell, as shown in Figure 2. e time step is Δt � 0.4s. Each time, the pedestrians determine the next move strategy according to the transfer probability of the cell in each direction until they reach the exit and escape. Burstedde et al. [30] introduced the concepts of static floor field and dynamic floor field in the CA model. Many scholars added friction, inertia, force, obstacle, and other factors to its expansion. e probability of pedestrian moving to its adjacent cell lattice is determined by static floor field, dynamic floor field, and their interaction. is paper analyzes the transition probability of static floor field and dynamic floor field, respectively. In this paper, the model is established on the basis of static field and dynamic field.

Static Floor Field.
e static floor field (SFF) in cellular automata is the distance from each cell to the effective exit. Let S ji (t) denote the shortest distance from the neighborhood cell i of evacuation individual j to the exit at time t. Let the location coordinate of evacuation individual j be C(x, y), and the location coordinate of neighborhood cell i be C(x i , y i ). ere are m effective exits in the evacuation space. e midpoint coordinate of each exit is C(x e k , y e k ). According to the static floor field of cellular automata, the transition probability formula is shown in the following equation: [1,m] x e k − x i 2 + y e k − y i [1,m] x e k − x i 2 + y e k − y i where the second term of the formula represents the minimum value of the neighborhood cell distance among all exits. e first term represents the maximum value of the minimum value of the neighborhood cell distance among all exits. e formula is a two-level progressive relation. e first layer is to calculate the shortest distance between eight cells in the neighborhood and all effective exits. e second layer is to compare the deviation between the shortest distance of eight cells in the neighborhood and the maximum of the shortest distance. If it is equal to the maximum value of the shortest distance, then S ji (t) � 0. If the deviation from the maximum value of the shortest distance is greater, then the cell has the shortest distance from the exit. at is, the greater the S ji (t) value is, the greater the probability of the cell being selected.

Dynamic Floor Field.
In the emergency situation, the state of the evacuation individual changes after the emotion reaches a certain threshold. us, it causes the change of behavior strategy. e behavior execution process of evacuation individual j at time t is shown in Figure 3. erefore, the dynamic field of this model mainly includes emotion updating rules and state updating rules.
where PE j (t − 1) is the emotional intensity of evacuation individual j at time t − 1; Δe j (t) is the emotion increment of evacuation individual j at time t; ξ is the perceptual attenuation coefficient; and the value range is [0, 1]. It indicates that there is a certain attenuation in the superposition of emotion increment of evacuation individuals. τ is the random emotional disturbance of individual j at time t. \Equation (3) shows the renewal mechanism of emotion and the evolution mechanism of infection. e emotional intensity PE j (t) is based on the self-emotional intensity at time t − 1 and superimposed with the emotional influence of others in the surrounding perception area.

Emotional Perception Increment.
e emotional fluctuation of evacuation individuals is related to their own speed and the emotion of others around them. erefore, this model defines the emotional increment function of evacuation individual j at time t, as shown in the following equation: where Δf j (v, t) is the internal emotional change of evacuation individual j caused by speed at time t; Δφ j (t) is the emotional change of others perceived by individual j at time t.
(1) Definition of function Δf j (v, t): in the process of evacuation, if the evacuation speed is slower than the maximum expected speed, the negative emotions (irritability, tension, and anxiety) increase faster. If the walking speed of the evacuation individual is close to the maximum expected speed, the growth of their bad mood is relatively slow [24]. erefore, the emotional increment function caused by speed is shown in the following equation: where v j (t − 1) represents the speed of the previous movement of the evacuation individual j at time t; v max j is the maximum expected speed of individual j during evacuation; α is the velocity sensitive parameter of evacuation individual j at time t; and the value range is [0, 1].
is function can clearly express the influence of the difference between the previous moving speed and the maximum expected speed on the emotional fluctuation. If the moving speed is close to the maximum expected speed, it will cause less fluctuation of individual emotion. On the contrary, if the moving speed is far less than the maximum expected speed, the Δf j (v, t) value will be close to 1. In other words, the degree of emotion will increase.
(2) e definition of function Δφ j (t): Zhang et al. [31] proposed that there were circular perception domains and perception domain radius for evacuation individuals. e traditional perception domain discretization thinks that the emotional state of each individual in the perception area can be fully perceived by the individual, but the actual situation is not so. Geng et al. [20] considered that the influence of information communication in line of sight is different from that in nonline of sight. At the same time, this model considers that the scope of emotional perception is not limited to contact perception. For individuals, the emotions of the surrounding people will be perceived, especially the emotions in the visual area will be more intense through visual perception. erefore, this model divides the perception area into visual perception domain RV and nonvisual perception domain RI. It is considered that the emotion of others in the visual domain can be fully perceived, while the emotion in the nonvisual domain can be perceived through pushing, sound, action, etc., and the perception degree will decline to a certain extent. Li et al. used angle π/4 to represent the visual range of pedestrians in eight neighborhood directions [32]. Usually the horizontal angle of people's binocular vision is about 3π/4. π/2 is selected as the angle of visual area in this model, considering that people's field of vision will shrink under emotional fluctuations. en, the parameters of the visual perception domain RV and the nonvisual perception domain RI are defined as follows: is the angle radian between the forward direction dir j (t − 1) of the evacuation individual cell j and the horizontal axis of the coordinate. Its value can be (0, π/4, π/2, 3π/4, π, 5π/4, 6π/4, 7π/4) ③ dir j (t − 1) is the forward direction of evacuation individual cell j at the last moment. e corresponding relationship with θ j (t) is shown in Figure 4 ④ e visual area of evacuation individual is a sector with arc π/2 ⑤ R is the emotional perception domain radius of evacuation individuals e definition of visual perception domain RV: the sector with the cell coordinate cell j (x, y) of evacuation individual j at time t as the center, [θ j (t) − π/4, θ j (t) + π/4] as the angle, and R as the radius is the visual perception domain RV of the evacuation individual at time t. As shown in Figure 5, the green color area is the visual perception domain after discretization. e emotions of the evacuated individuals in this area can be perceived.
Definition of nonvisual perception domain RI: the circle with the cell coordinate cell j (x, y) of evacuation individual j at time t as the center, R as the radius, and removing the part of RV is the nonvisual perception domain RI of the evacuation individual at time t. As shown in Figure 5, the gray area is the nonvisual perception domain. e emotion of the evacuation individual is partially perceived. erefore, the visual perception domain and nonvisual perception domain, corresponding to different forward directions of evacuation individual cell j at time t, are shown in Figure 6.
According to the perception domain, the evacuation space is divided into three situations: exit perception area A 1 , adjacent edge area A 2 , and other area A 3 , as shown in Figure 7.
As shown in Figure 7, exit perception area A 1 is composed of two squares with exit endpoint as vertex and R as side length and a rectangle with length R and width el. When evacuation individuals are in this area, they are not affected by other people's emotions, because the exit is the most attractive. erefore, it can be considered that when in the exit perception area A 1 , the emotion of the evacuation individual no longer fluctuates, such as the location of n 1 . e perception area of evacuation individuals in the adjacent edge area A 2 is an incomplete circle, which needs to be simulated according to the incomplete rules, such as the location of n 2 and n 3 . When the evacuation individual is in the other area A 3 , the visual perception domain and nonvisual perception domain are complete, and the simulation calculation can be carried out according to the complete perception domain, such as the location of n 4 and n 5 .
Nilsson et al. used empirical experiments to analyze other people's influence and behavior choice during indoor evacuation. ey proved that the influence of close distance on evacuation individual's emotion and behavior is greater than that of long distance [33]. Fu et al. [29] also believed that the evacuation individuals were more sensitive to the emotional perception of nearby areas during the evacuation process. In this model, the influence of other people's emotion in visual perception domain and nonvisual perception domain on the evacuation individual is different. Combined with the distance factor, this model considers that Discrete Dynamics in Nature and Society the perceived degree of others' emotion in the visual perception domain is greater than that in the visual perception domain and is greater than that in the nonvisual perception domain. Taking Figure 8 as an example, evacuation individual j is more likely to be affected by n 6 , then by n 7 , and finally by n 8 .
In reference [29], the formula of distance as influence weight in emotion gain function is evacuation individual j to target cell k. is formula reflects that evacuation individual j is more sensitive to the emotion of close evacuation individual, and the emotion of others far away from evacuation individual j has less influence. In order to reflect the difference of perception degree in three different situations, this model redefines the distance weight as the following equation: e average change function of other people's emotion perceived by evacuation individual cell j at time t is shown in the following equation: where g 1 is the number of evacuees in the visual perception domain RV of evacuees cell j; g 2 is the number of evacuation individuals in the nonvisual perception domain RI of evacuation individual cell j; η is the degree to which the emotional difference of other people in the nonvisual domain is received by the evacuation individual j. Its value is not greater than the minimum distance weight in the visual perception domain, that is, 0 ≤ η ≤ min k∈RV dd jk . e first term of equation (7) represents the average increment of emotional change of evacuation individual j in the visual domain at time t. It also shows that the amount of emotional change is related to distance, and the emotional difference in close distance has a great influence on evacuation individual j. e latter term indicates that the emotional change of evacuation individual j in the nonvisual domain is the overall weakening due to the absence of specific situation.
In conclusion, the emotional increment function of evacuation individual j at time t is updated as shown in the following equation: e incremental function of emotional intensity indicates that the change of speed will cause the emotional fluctuation of evacuation individuals. When the travel speed is close to the maximum expected speed, the fluctuation of emotion is small. e greater the difference between the travel speed and the maximum expected speed, the greater the emotional fluctuation. e emotions of others around also affect the emotional fluctuation of evacuation individuals [29]. e emotional influence of others in the visual perception domain is greater than that of people in the nonvisual perception domain.
at is, the emotional influence of others in the visual perception domain is related to distance [24], and the emotion of evacuation individuals in close distance is easy to be perceived. It can be seen that the constructed incremental function of emotional intensity is more in line with the actual changes of individual emotions.

Mechanism of Emotional Renewal.
e iterative process of evacuation strategy of evacuation individual j at time t is to judge the state update by emotion update and implement evacuation strategies corresponding to different states. e emotion of evacuation individual j at t time is related to the value of emotion, action strategy, state, and direction at time t − 1. It is also affected by the emotional intensity of others within the range of emotion perception. e change of emotion will cause the dynamic change of evacuation individual's state and behavior strategy. e specific state update process is shown in Figure 9.

State Update Mechanism.
Emotion gathered to a certain extent, resulting in the change of evacuation individual state. According to the SIS model, the evacuation state of the individual is defined as "emotional susceptible" and "emotional infective." e state transition rules follow SIS model. e initial definitions of "emotional susceptible" and "emotional infective" are as follows: when the emotional intensity PE j (t) of individual j exceeds the threshold λ emotion at time t, it is transformed into "emotional infective" by probability β. When the emotion fluctuation of "emotional infective" is lower than the threshold λ emotion , it is transformed into "emotional susceptible" by probability c. In the process of evacuation, S is the proportion of "emotional susceptible" and I is the proportion of "emotional infective." Its dynamic differential dynamic equation and system conservation are the same as equation (1).
Different evacuation individuals will adopt different evacuation strategies because of their different states. Literature [34] analyzed the 2010 German "Love Parade" stampede accident video and observed that the crowd presented two different behavioral characteristics, namely, "positive characteristics" and "negative characteristics." People with "positive characteristics" will actively seek exit, accelerate forward, and have obvious herd behavior; people with "negative characteristics" are slow, try to avoid the crowd, and have no obvious herd behavior. Based on this, this model maps pedestrian emotional intensity to pedestrian behavior. During the evacuation, the "emotional susceptible" (S) adopts the "calm" strategy BS, the speed is vs, and tends to choose the direction with low density ρ; while the "emotional infective" (I) adopts the "impulsive" strategy BI, the speed is vi, and tends to choose the direction with high density ρ.
PE j (t) is defined as the emotional intensity of evacuation individual j at time t, and the value range is (0, 1]. at is, the emotion of evacuation individual will not disappear in the whole process of evacuation. According to the different emotional intensity, the evacuation behavior is different. λ emotion is defined as the threshold value of emotional intensity change. It is the critical value from psychology to behavior, and the value range is (0, 1]. e mapping relationship between emotional intensity and individual behavior type is shown in Table 1.

Transfer
Strategy. In this model, cellular automata evacuation model is used. In each time step, the movement direction of the individual who does not leave the evacuation space is determined according to the next time movement probability. Let P ji (t) denote the probability of individual j transferring to cell C i at time t. e evacuation individual will choose the neighborhood cell with the maximum transition probability as the target and move at the speed v j . According to the theory of cellular automata, the transition probability is related to three basic characteristics such as the Discrete Dynamics in Nature and Society distance from neighborhood cell C i to the exit, whether it is an obstacle cell. At the same time, this model considers that the crowd density of the forward direction is also a factor that determines the selection of target cells. erefore, the cell transition probability is constructed as shown in the following equations: where k s is the influence coefficient of distance to the exit; k ρ is the influence coefficient of density. cp ji (t) is the state parameter of neighborhood cell i of evacuation individual j at time t. If it has been occupied by evacuation individuals, the value is 0; otherwise, it is 1, as shown below: where co ji (t) indicates whether the neighborhood cell i of evacuation individual. at time t is an obstacle. If it is an obstacle, the value is 0; otherwise, it is 1, as shown below: co ji (t) � 0, C x i , y i is barrier, where ρ ji (t) is the density around the neighborhood cell i of the evacuation individual j at time t. At present, all other evacuation individuals in the visible area of evacuation individuals will affect their choice. is model gives the density influence function of different types of groups combined with the decision types. erefore, the transfer probability function of different types of crowd is different.
(1) "Emotional susceptible" S adopts "calm" strategy BS. e maximum expected speed is v max j � 1 m/s. at is, when each step is updated, it moves 1 step. e target cell of "calm type" evacuation is 8 cells around, as shown in Figure 10.
At the same time, the density range of "calm type" evacuation individuals is the close density around the neighborhood cell i. As shown in Figure 10, the red cell is   8 Discrete Dynamics in Nature and Society evacuation individual j, and the blue cell is neighborhood cell i. When analyzing the density of i, evacuation individual j considers the density of the gray cell. Because the "calm" evacuation individuals described in this model choose the direction with low crowd density, the transfer probability formula of this kind of crowd is as follows: (2) "Emotional infective" I adopts "impulsive" strategy BI. e velocity is vi. ey tend to choose the direction with high crowd density ρ.
e maximum expected speed is v max j � 2m/s. When each step is updated, it moves 2 steps. en, the target cell of "impulsive" evacuation is two-step cell in each direction. Suppose the group of target cells is set G 1 , as shown in the blue cells in Figure 11(a). At this time, the moving speed is vi � v max j . If all the preferred eight cells are occupied, the target cell with the largest transition probability of eight cells in neighborhood is selected. Suppose the neighborhood cell is set G 2 , as shown in the blue cells in Figure 11(b). At this time, the moving speed is vi � 1m/s. e radius of density range of "impulsive" evacuation individuals is two steps, as shown in Figure 12.
Because the "impulsive" evacuation individuals described in this model choose the direction with high crowd density, the calculation method of the transition probability of this kind of crowd is as follows: Calculate P ji (t) i∈G 1 according to equation (9). If then calculate P ji (t) i∈G 2 .

Simulation Process.
According to the emotional intensity of individual evacuation, the individual evacuation is divided into two groups: "emotional susceptible" and "emotional infective." e "calm" strategy and "impulsive" strategy are adopted in evacuation. e emotional intensity of individual evacuation is spread by others and fluctuates, and the evacuation strategy of each individual evacuation is not single. Moreover, due to the influence of dynamic changes, the emotion of evacuation individuals affects their mobile behavior [34]. Considering that only when emotion accumulates to the time step of cell movement can pedestrian movement be driven, the improved model adopts parallel update mechanism with a unified time step Δt. In each time step Δt, the emotion is updated firstly to judge the change of the emotional state. en, the behavioral decisionmaking is updated. e simulation process of this model is as follows: ① Initialize the number of people and the layout of evacuation space; ② Calculate the SFF of each cell according to equation (2); ⑤ Update individual status: if the emotional intensity of "emotional susceptible" S is higher than the threshold λ emotion , it will turn into "emotional infective" I with probability β; if the emotional intensity of "emotional infective" I is lower than the threshold λ emotion , it will turn into "emotional susceptible" S with probability c. Realize the dynamic change of evacuation crowd state, as shown in Figure 13; ⑥ Calculate the individual transfer strategy according to equations (9)- (14) and select the target cell; ⑦ Solve conflicts: when multiple evacuees choose the same target unit, there will be conflicts. In this case, an evacuation individual is randomly selected to move. For simplicity, it is assumed that the probability of each individual being selected is equal [31]; ⑧ When all the evacuees in the space leave, the simulation ends; otherwise, jump to ④; ⑨ Every time the evacuation individual moves, pedometer + 1, until the end of the simulation.

Experiment 1.
is experiment analyzes the influence of perception domain on the total evacuation time. In the improved model, the emotion perception domain of the evacuation individual is a circle with a radius R. Set a single exit room as the simulation scene. e size is 12m × 12m. e exit width is el 1 � 3cl. e size of each cell is cl × cl, cl � 0.4m. e initial number of people is N � 200, which is randomly distributed in the space and given the initial emotion value. In order to analyze the influence of perception domain radius R, set nn � R/cl (cl � 0.4m). en, nn � 1, 2, 3, 4, 5. Emotional intensity threshold λ emotion , proportion of emotional infected people I, emotional infection coefficient β, and emotional calm coefficient c and other parameters all affect the overall evacuation time of the system. erefore, three different scenarios with different parameters are set for simulation experiments. Scenario 1: β � 0.33, c � 0.67; scenario 2: β � 0.5, c � 0.5; scenario 3: β � 0.67, c � 0.33. Next, we analyze the impact of different perception domain radius nn on the total evacuation time under three scenarios, when the emotional threshold λ emotion � 0.2, 0.4, 0.6, 0.8. Each experiment runs 100 times, and the average of the total evacuation time is taken. e simulation results are shown in Figures 14-16. It can be seen from the figure that, when λ emotion � 0.2, there is no difference in the total evacuation time between scenario 1 and scenario 3, while the total evacuation time of the perceived radius nn � 4.5 is smaller at scenario 2; When λ emotion � 0.4, the total evacuation time of the perceived radius nn � 4.5 is less than nn � 1, 2, 3 in the three scenarios. When λ emotion � 0.6 and λ emotion � 0.8, the three scenarios show that the minimum total evacuation time appears when nn � 1, which generally shows that the larger the perception domain radius is, the larger the total evacuation time is. erefore, when λ emotion < 0.5, the larger the perceived radius nn is, the faster the total evacuation time is; when λ emotion > 0.5, the larger the perceived radius nn is, the longer the total evacuation time is. e values of β and c have no decisive influence on this conclusion. In order to verify this conclusion, we carried out simulation experiments on the same three scenarios with the initial numbers of people N � 300 and N � 400, and the results are consistent.

Experiment 2.
e perceived domain radius of individual evacuation has a direct impact on the emotion of individual evacuation. We carried out simulation experiments on the mean system emotions in scenario 1, scenario 2, and scenario 3 under different emotion thresholds. e experimental results are shown in Figure 17. In the figure, the horizontal axis is the perception domain radius nn � 1, 2 , 3, 4, 5. e vertical axis is the mean emotional value of the whole system, and the formula is as follows: , that is, the cumulative sum of the mean system emotion at each time divided by the total evacuation time.
e value of this parameter reflects the emotional level of the individual in the process of evacuation. e initial position and initial emotional value of the experiment were the same as that of experiment 1. e results show that, when λ emotion � 0.2 and λ emotion � 0.4, with the increase of perception domain radius, the mean system emotion increases; when λ emotion � 0.6 and λ emotion � 0.8, with the increase of perception domain radius, the mean system emotion decreases. e conclusion of this experiment is that λ emotion � 0.5 can be regarded as a cut-off point of emotional value. e initial people numbers of N � 300 and N � 400 are simulated, and the conclusion is consistent.

12
Discrete Dynamics in Nature and Society

Experiment 3.
e change of emotional value will lead to the change of evacuation individual state. is experiment analyzes the influence of different perceived domain radius of evacuation individuals on the state of evacuation individuals. According to the previous two experiments, λ emotion � 0.5 is the inflection point of the change of evacuation time and mean system emotional value. In order to observe the detailed changes of pedestrian status more clearly, λ emotion � 0.4, λ emotion � 0.5, and λ emotion � 0.6 are taken in this experiment respectively, to analyze the change of the proportion of state I in the system at each time. e experiment adopts the same scene as the previous experiment. e initial position and initial emotional value are the same. e initial number of people are N � 200 and N � 300 respectively. Take β � 0.5 and c � 0.5. e experimental results are shown in Figures 18 and 19.
It can be seen from the figure that no matter N � 200 or N � 300, the system will be transformed into group I eventually with the evacuation time when λ emotion � 0.4. When the perception domain radius nn � 1, the change speed is slower, and when nn � 2, 3, 4, 5, the change speed is faster. e proportion curve of state I of perception domain radius nn � 1 at λ emotion � 0.5 is different from other situations significantly. When N � 200, the proportion of I increases first and then decreases rapidly. In the last stage of evacuation, group I leaves first, and the last one is individual S. So, at the last moment, S � 1. When N � 300, the number of group I decreased slowly. is shows that for nn � 1, it is more beneficial to transform into group I when the number of people is small. But state group I and state group S always coexist in the system under the critical threshold no matter what kind of initial number. Under two initial numbers, when the perception domain radius nn � 2, 3, 4, 5, the evacuation individual status in the system is ultimately group I. From the curve change, when nn � 2, the change speed of I state group is slightly slower than that of nn � 3, 4, 5. When λ emotion � 0.6, whether the initial number is N � 200 or N � 300, with the change of evacuation time, the system will be all group S. It can be seen that the state change speed of perception radius nn � 1 is the slowest, followed by nn � 2. When N � 200, the change speed of group I of nn � 3 is slightly slower than that of nn � 4 and nn � 5, while when N � 300; the difference of curves of nn � 3, 4, 5 is small. erefore, when λ emotion ≤ 0.5, the larger the perceived domain radius of the evacuation individual is, the more favorable it is for the evacuation individual to change into state I. When λ emotion > 0.5, the larger the perceived domain radius of the evacuation individual, the easier the evacuation individual is to change to state S.

Experiment 4.
e improved model and CA model are compared in the same scenario to verify the effectiveness of the model. e basic scene of reference [24] is used as the simulation scene. e unit cell is cl × cl, cl � 0.4m. e single exit room size is 18cl × 14cl, and the exit width is el 1 � 2cl. When the initial state t � 0, 150 individuals were randomly distributed in the space, each occupying a cell. Moore's neighbor rule is used when evacuation individuals move. Other parameters of this model are set as λ emotion � 0.6, ξ � 0, β � 1, and c � 1. CA model does not consider emotional contagion, only according to the SFF rule evacuation model. After the evacuation starts, the evacuation individuals move to the exit according to different model moving rules. Each time step is 0.4s. Figures 20 and 21 record the density maps of the two models at different time steps, respectively. e cell density function is shown in equation (15). It follows Moore's neighbor rule, so num � 8 in the formula. cp ji (t) denotes whether there are people in the neighbor cell of cell j at time t. If cp ji (t) � 1, the neighbor cell has people; otherwise, the neighbor cell has no people. 14 Discrete Dynamics in Nature and Society As can be seen from the figure, the personnel movement trend of the improved model is consistent with the CA model. e pedestrian flow at the exit is arched, which is in line with the characteristics of individual movement flow. In the same time step, the evacuation individuals move faster in the improved model.
is is because the emotional contagion of evacuation individuals leads to the dynamic movement speed, which reduces the overall evacuation time.

Experiment 5.
In order to further verify that the improved model can reflect the movement characteristics of individual evacuation, the scene of reference [24] is used to carry out the experiment. N people to be evacuated are distributed in a single exit room with the size of 12m × 12m randomly. Each cell is cl × cl, cl � 0.4m. at is, there are 900 cells in a single exit room. Suppose the total number of people is N � 100, 200, 300, 400, 500, 600, 700, 800, and the exit width is 2, 3, 4, 5, 6, 7, and 8 respectively.  It can be seen from the figure that the trend presented by the two models is consistent. e wider the exit width is, the smaller the evacuation time is. e evacuation time of the improved model is shorter than that of the CA model under the same number of people and the same exit width. e reason is that the improved model takes into account the emotional contagion of individual evacuation, and different emotional groups adopt different evacuation strategies.

Conclusion
In this paper, the SIS model of infectious diseases is combined with cellular automata, and a crowd evacuation model considering emotional contagion is proposed. In the improved model, the dynamic movement rules of evacuation individuals are constructed. e action strategy, state transition rules, emotion update function, and linkage relationship among them are refined.
ree simulation experiments were carried out to study the influence of perception domain radius on the total evacuation time, the mean system emotion, and the proportion of state group I under different emotional threshold, infection coefficient, and calm coefficient. e results show that when the system emotion threshold λ emotion ≤ 0.5, the larger the perceived domain radius, the higher the mean system emotion value, the easier the evacuation individual is to change to state I and leave the site quickly. When λ emotion > 0.5, the larger the perceived domain radius, the lower the mean system emotion value, and the easier the evacuation individual is to change to state S, resulting in a longer evacuation time.
e conclusion can provide the following evacuation guidance for evacuation individuals. When λ emotion ≤ 0.5, the field of vision of the evacuation individual can be enlarged and the observation area is far away.
is is conducive to judging the critical situation, forming a convergent emotion, and then make positive action strategies. When λ emotion > 0.5, the evacuation individual pays attention to the emotion or behavior strategy of others nearby, which is more conducive to change personal emotion and state. If the observation scope is too large, it may lead to some misjudgments and adopt a calm strategy, resulting in prolonged evacuation time.
Two simulation experiments are compared with the CA model. e results show that the improved model can effectively simulate the process of individual evacuation, showing the characteristics of crowd conformity and exit arch. In addition, the comparative analysis also shows that emotional contagion has an impact on evacuation time. Because different emotional individuals choose different evacuation strategies, evacuation individuals move faster than the traditional CA model.
In summary, this model can not only realize the dynamic characteristics of individual state and emotion in the macroevacuation model but also analyze the microcharacteristics of individual state change and emotion fluctuation. is conclusion provides a basis for managers to grasp the characteristics of emotional contagion among people and put forward effective control strategies. Of course, the content of this chapter needs to be improved, and the simulation results need to be supported by empirical data. In the future, empirical experiments will be conducted to verify the significance of the conclusions in this chapter.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.  16 Discrete Dynamics in Nature and Society