# The Impact of Media Coverage and Emergency Strategies on the Rumor Spreading

Authors: Liang’an Huo; Yingying Cheng
Journal: Discrete Dynamics in Nature and Society
 (2018)

---

## Abstract

On the basis of nonlinear dynamic model, we research the propagation of rumors after emergencies, describe the impact of media coverage and emergency strategies by government on the transmission dynamics of information, and then obtain the basic reproduction number of rumor spreading. In order to overcome the limit of traditional methods of the static decision problem, the dynamic optimal control model for the rumor spreading is proposed based on the theorem of the optimal control. An optimal objective based on the maximum social utility is established and the optimal solution is acquired by using the Pontryagin Maximum Principle. Finally, numerical simulations show that the dynamic optimal control method has obvious superiority in modeling compared with the method without control. By means of the dynamic optimal control of media coverage and emergency strategies for government, the final rumor scale and the peak of spreaders can be effectively reduced.


---

## Body

## 1. Introduction

When emergencies occur, people's thirst for information is different from normal. If the propagation of information is not transparent, detailed, and controlled by the government and the media at this time, it may lead to people’s psychological and emotional tension, unnecessary public panic, and economic loss [1–4]. On the contrary, if there is a clear channel of information dissemination and effective information dissemination mode among the government, the media, and the public, this will be beneficial to reduce the social panic and prevent the expansion and dissemination of adverse effects in the event. At the same time, the information interactive platform can better establish the credibility of the government and the media.In recent years, many scholars have discussed the rumor spreading from different perspectives through qualitative analysis and theoretical modeling. Rapoport [5] first uses the epidemic model to study the spread of information. Goffman and Newill [6] compare the similarities between the spread of disease and rumor and use the stability of the epidemic model to explain the final propagation state of the rumor. Daley and Kendall [7] put forward a mathematical model of rumor spreading in 1960s. After that, scholars have expanded the classic D-K model. Based on summarizing the similarities and differences between the spread of infectious diseases and the mechanism of rumor diffusion, Klaus [8] emphatically expounds the reasonableness of using the epidemic model to explain the spread of rumors.With the development of dynamic theories, rumor spreading based on dynamics has attracted more and more attention. Dickinson [9] compared the spreading of epidemic with the rumor spreading mechanism. And they considered the propagation is a time-varying process, which can be subsumed under a general class of transient processes whose time-dependent evolution can be characterized exactly. According to the characteristics of audience, Thompson [10] divide the crowd into two categories, active and indifferent, set up a dynamic model of rumor spreading, and analyze the strategy of controlling the spread of rumors. Based on this, a large number of scholars [11–14] have expanded the rumor propagation model from the perspective of activity.Media reports are known as “subcutaneous injection” and “magic bullet” theory, and it is believed that media coverage has a great influence on public’s understanding and emergency management [15]. In the outbreak of large-scale infectious diseases, the media coverage will affect the final spread of disease; moderate coverage can limit the spread of infectious diseases and reduce the loss of infectious disease transmission. However, overreporting can cause panic and also bring disastrous consequences out of control [16, 17]. Referring to the theory of disease, Zhao et al. [18] built a model verifying that the media with high ratings but low credibility will cause widespread panic and cause a lot of rumors. Jane M Heffernan and Shannon Collinson [19] have applied the “media function” that affects the transmission rate in the mathematical epidemiology model, which allowed the effects of mass media during disease outbreaks to be studied in mathematical modeling literature. Jean M. Tchuenche and Chris T. Bauch [20] established a susceptible-infected-hospitalized-recovered model with vital dynamics, in which the media coverage of disease incidence and disease prevalence can influence people to reduce their contact rates. And the media function is merged into the model by using the exponential decreasing function skillfully. Luo et al. [21] explored the impact of positive news on rumor spreading in scale-free networks. It is concluded that most of the rumors related to hot events or emergencies can be propagated rapidly on the hotbed of online social networks. Huo et al. [22, 23] have made a great contribution in rumors spreading, and their results indicate that the media can be regarded as a good indicator in controlling the occurrence and spreading of rumor.When an emergency occurs, it will cause or may cause serious social harm. In order to avoid losses, some measures must be taken as soon as possible in response to the sudden incidents, such as natural disasters, accident disasters, public health events, and social security incidents. However, the main effect of media coverage is to popularize science education, improve the identification of the ignorant, reduce the probability of the ignorant believing rumors, and thus reduce the spreading rate of rumor dissemination, rather than prevent the spread of rumors. On the contrary, when rumors break out, relevant emergency strategies can not only prevent individuals from believing in rumors, but also reduce the spreader's enthusiasm for dissemination. Huo et al. [24], based on utility theory, proposed a model to describe the interplay between rumor spreading and authorities’ actions. Through numerical simulation, it is found that positive actions in the government can improve rumor management in emergencies and generate positive social effects. To study the spread of rumor in emergencies, Chen et al. [25] investigated a novel rumor spreading model with latent, constant recruitment, and varying total population. Li et al. [26] explored the impact of punishment of governments and sensitivity of individuals on the rumor spreading and then drew a conclusion that improving the punishment of government and increasing the sensitivity of individuals can effectively control the spreading of rumor.The spread of rumor not only caused the people to panic, but also brought huge economic losses. Therefore, many scholars have made great contributions to the immune control of rumors [27, 28]. However, controlling rumors is only the first goal. Saving control costs is equally important. Thus, optimal control theory has been widely applied in rumor governance in recent years. Based on the controlled heterogeneous node-based SIRS model, Huo et al. [29] considered an optimal control of rumor spreading model with psychological factors and time delay. Kandhway et al. [30] based on SIS and SIR models considered fixed and time varying information spreading rate and formulated two optimal control models, and then they proved the effectiveness of the optimal control strategy against the constant control strategy, a heuristic control strategy and no control. Huo et al. [31] established an optimal control model based on the SEIR model with latent period to study rumor management strategy. Inspired by the classical SIR model, Xu et al. [32] proposed a new rumor spreading model, which regards the transfer probability of nodes from susceptive to infected as a variable instead of a constant. Jeong et al. [33] proposed a rumor model with three control strategies for preventing the spread of rumor: (1) announcing the truth before ignorant receives rumor, (2) punishing spreaders, and (3) deleting information of the rumor in media, and considering optimal control problems to minimize the number of spreaders while minimizing the cost of three control strategies for preventing the spread of rumors.From the previous researches, the scholars mainly focus on the state of the final spreading, but ignore the influence factors and the control problems during the process. The individual’s spreading behavior will be determined by the amount of information obtained by an individual preference; meanwhile, the media and government measures will affect the rumor spreading. The facts indicate that the popularization of scientific knowledge, emergency strategies, authoritative information published by the media, expansion of information coverage, and guiding public opinion are the key to rumor governance in emergency management.Considering the positive impact of media coverage and the emergency strategies of government, we extended the classical DK model. On the basis of the optimal control theory, the optimal control strategy of rumor spreading is discussed, and a feasible management and control scheme is proposed for the control of the rumor in emergency management and obtainment of an emergency management strategy with maximization of social utility.The organization of this paper is as follows: the rumor spreading model is formulated in the next section. Section2 describes an analysis of the equilibrium, giving existence conditions and a kind of threshold condition for rumors outbreaks. Section 3 analyzes the optimal control problem in detail. We present some numerical simulations to illustrate the established results in Section 4. Finally, Section 5 outlines the work and some topics that are worthy of study.


## 2. The Mathematical Model with the Impact of Media Coverage and Emergency Strategies

The spread of rumors will be influenced by media coverage. The media coverage improves personal discrimination ability through daily popular science information, thereby reducing the spreading rate of rumors. On the other hand, in order to reduce the loss caused by emergencies, it is necessary to take timely and effective emergency strategies in dealing with rumor spreading. Such measures include holding a press conference and releasing relevant facts for emergencies, which can effectively prevent the spread of rumors.Based on the study of the classic D-K model, this paper adds two factors, media coverage and emergency strategies by government, to the spread of rumors, analyzes the spread mechanism of the rumor under the influence of these two factors, and lays a theoretical foundation for further discussion of the control strategy of the rumors.


### 2.1. Model Hypothesis

( 1 ) Based on the classical D-K model for the classification and propagation of the population, this paper further assumes that the ignorant propagates rumors with a certain probability when the ignorant meets the spreader, influenced by the media reports, and the propagation probability is gradually reduced, and finally a certain stability value is reached.( 2 ) When the rumor spreaders meet with the stifler, the spreader loses interest in the spreading; when the two spreaders meet, each of them loses the power to spread the rumor.( 3 ) Rumors spread in a certain range of communities, the total number of people is variable, and individuals entering the community do not know any information.


### 2.2. Formulation of the Rumor Model

According to the basic assumptions of the model, an open and mixed population in a certain area can be divided into 3 categories: the ignorant (nodes who have never heard rumors), the spreaders (nodes who know and spread rumors), and the stiflers (nodes who understand the truth of rumors and stop spreading rumors). Similar to the traditional models, at any time each and every node in propagation network is assumed to be in one of three possible categories. LetIk(t) (Sk(t),Rk(t)) denote the relative densities of ignorant (spreaders, stiflers) with degree k at time t, respectively, with the normalization conditions Ik(t)+Sk(t)+Rk(t)=1. And in order to incorporate the dimensions of massive news coverage, we introduce the information variable M(t) to describe the cumulative density of media coverage in other region at time t. This is because the contact between different types of individuals leads to the spread of rumors and then changes in the type of individuals, as shown in Figure 1.Figure 1
Schematic diagram of rumor spreading under the influence of media coverage and emergency strategies.Based on the idea of differential equation modeling, the D-K model is further extended to establish a nonlinear dynamic model considering the impact of media coverage and emergency strategies by government. According to the above assumptions, one has the following system of3n+1 ordinary differential equations (ODEs):(1)dIktdt=Λ-β1-β2kMtα+kMtkIkt∑k′Pk′∣kSk′t-ηIktdSktdt=β1-β2kMtα+kMtkIkt∑k′Pk′∣kSk′t-ηSkt-δSkt-εkSkt∑k′Pk′∣kSk′t+Rk′tdRktdt=εkSkt∑k′Pk′∣kSk′t+Rk′t+δSkt-ηRktdMtdt=m0+ω∑k′PkSk′t-λMtwhereΛ: immigration rate, the individual enters a particular community at a constant rate;α: the half-saturation constant, which reflects the impact of media coverage during rumor spreading;η: removal rate, each type of individual moves out of the group for some reason;β1: the probability that the ignorant accepts and spreads rumors when the ignorant connects with the spreaders. β1IS indicates that some of the ignorants choose to trust the spreaders’ information and enter the group;β2: the maximum impact of media on the probability of transmission when ignorants are exposed to spreaders. β2kM/α+kM expresses that the ignorant is in contact with the spreaders; the probability of the media’s interference with the spread of rumor, with the increase of the number of spreaders, the intensities of the media, and the gradual increase in the media's interference in believing and disseminating untrue information to the individual, eventually tends to a maximum β2, where α represents saturation coefficient and measures the degree of media coverage affecting rumor spreading. Considering that media coverage can slow down the spread and diffusion of rumor, it cannot completely prevent the spreading, so we assume that β1≥β2;ε: stifling rate, the probability that a spreader transfers into a stifler when the spreader contacts with a spreader or a stifler;δ: emergency efficiency, by means of isolation and other compulsory administrative measures, the government loses the spreading ability of rumor spreaders;ω: rate of media coverage being implemented, the influence is bidirectional. While the media is affecting rumor spreading, the development of rumors in turn affects the frequency of media coverage. And the parameter ω is a proportional constant with which media coverage is being implemented based on the current rumor spread;λ: depletion rate of media coverage due to ineffective implementation;m0: density level of media coverage in this region based on the spread of rumor in another region. Even if there is no rumor case in the area, we must give warnings and nip the impossible situation in the bud.According to the previous studies,Pk′∣k=k′P(k′)/〈k〉 indicates the degree correlations so that a node of degree k is connected to a node of degree k′, and 〈k〉=∑kkP(k) is the average number of contacts each node can make. The media coverage function kM(t)/α+kM(t) is similar to the Michaelis –Menten function P1(M)=CM(t)/1+DM(t), (C,D≥0) (see [34, 35] and the references cited therein), which is a linear function of current prevalence of rumor and measures the impact of media coverage on a node in group k. It is related to the number of nodes he/she contacts and current density level of media coverage.


### 2.3. Model Analysis

This paper analyzes the mathematical model based on the stability related theory of differential equations, and due to the fact that the model depicts the dynamic propagation process of rumor in the crowd, the parameters involved are assumed to be nonnegative.For the sake of simplicity, we setΘ(t)=∑k′Pk′∣kSk′(t)=1/〈k〉∑k′k′P(k′)Sk′(t) representing the probability that a randomly selected neighbor of a given node is a spreader while Φ(t)=∑k′Pk′∣kRk′(t)=1/〈k〉∑k′k′P(k′)Rk′(t) denoting the probability that a randomly selected neighbor of a given node is a stifler. Then, system (1) can be rewritten as(2)dSktdt=β1-β2kMtα+kMtk1-Skt-RktΘt-ηSkt-δSkt-εkSktΘt+ΦtdRktdt=εkSktΘt+Φt+δSkt-ηRktdMtdt=m0+ω∑k′PkSk′t-λMtWith no loss of generality, we assume thatβ1=β2=β, then the system (2) will be transformed into the following system (3).(3)αβα+kMtk1-Skt-RktΘt-ηSkt-δSkt-εkSktΘt+Φt=0εkSktΘt+Φt+δSkt-ηRkt=0m0+ω∑k′PkSk′t-λMt=0Obviously, system (3) has a unique rumor-free equilibrium E0=0,0,…,0,m0/λ1×(2n+1).Furthermore, the existence of rumor spreading equilibrium pointE∗ will be discussed. From system (3), a simple calculation shows that(4)E∗=S1∗t,S2∗t,…,Sn∗t,R1∗t,R2∗t,…,Rn∗t,m0+ω∑kSk∗tPkλ1×2n+1(5)Sk∗t=αβηkΘHRk∗t=αβηkΘ·δ+kεΘ+ΦHM∗=m0+ω∑kSk∗tPkλwhere H=αη+kηM∗+αβkΘ·η+δ+kε(Θ+Φ).Following the recipe of van den Driessche and Watmough [34], System (4) can be written as dx/dt=F(x)-υ(x) and x=(Sk,Rk)T, where x1,x2,…,xn represent S1,S2,…,Sn, and xn+1,xn+2,…,x2n denote R1,R2,…,Rn, respectively. Applying the next generation matrix method in [34],where the rate of appearance of new spreaders is(6)Fx=αβkα+kMt1-Skt-RktΘt02n×1and the transfer rate of individuals out of compartments is(7)Vx=εkSktΘt+Φt+δSkt+ηSkt-εkSktΘt+Φt-δSkt+ηRkt2n×1the Jacobian matrices ofF(x) and υ(x) at the rumor-free equilibrium E0=0,0,…,0,m0/λ1×(2n+1) are as follows:(8)F=DFE0=F110002n×2n(9)V=DVE0=δ+η0-δηE2n×2nwhere(10)F11=αβk1α+1Mt1P11α+1Mt2P2⋯1α+1MtnPn2α+2Mt1P12α+2Mt2P2⋯2α+2MtnPn⋮⋮⋱⋮nα+nMt1P1nα+nMt2P2⋯nα+nMtnPn(11)V-1=1δ+η0δδ+η1η∙E2n×2nAccording to the concept of next generation matrix and reproduction number given in [34], the reproduction number of system (2) equals R0=ρFV-1, where ρFV-1 denotes the spectral radius of the matrix FV-1. Then,(12)R0=ρFV-1=αβλη+δk∑k=1nk2Pkαλ+km0The basic reproductive numberR0 describes the average number of new spreaders at the beginning of the spreading process. And according to [34], if R0<1, then there is no more rumor to spread in the networks; while if R0>1, then the rumor will persist in the networks.


## 2.1. Model Hypothesis

( 1 ) Based on the classical D-K model for the classification and propagation of the population, this paper further assumes that the ignorant propagates rumors with a certain probability when the ignorant meets the spreader, influenced by the media reports, and the propagation probability is gradually reduced, and finally a certain stability value is reached.( 2 ) When the rumor spreaders meet with the stifler, the spreader loses interest in the spreading; when the two spreaders meet, each of them loses the power to spread the rumor.( 3 ) Rumors spread in a certain range of communities, the total number of people is variable, and individuals entering the community do not know any information.


## 2.2. Formulation of the Rumor Model

According to the basic assumptions of the model, an open and mixed population in a certain area can be divided into 3 categories: the ignorant (nodes who have never heard rumors), the spreaders (nodes who know and spread rumors), and the stiflers (nodes who understand the truth of rumors and stop spreading rumors). Similar to the traditional models, at any time each and every node in propagation network is assumed to be in one of three possible categories. LetIk(t) (Sk(t),Rk(t)) denote the relative densities of ignorant (spreaders, stiflers) with degree k at time t, respectively, with the normalization conditions Ik(t)+Sk(t)+Rk(t)=1. And in order to incorporate the dimensions of massive news coverage, we introduce the information variable M(t) to describe the cumulative density of media coverage in other region at time t. This is because the contact between different types of individuals leads to the spread of rumors and then changes in the type of individuals, as shown in Figure 1.Figure 1
Schematic diagram of rumor spreading under the influence of media coverage and emergency strategies.Based on the idea of differential equation modeling, the D-K model is further extended to establish a nonlinear dynamic model considering the impact of media coverage and emergency strategies by government. According to the above assumptions, one has the following system of3n+1 ordinary differential equations (ODEs):(1)dIktdt=Λ-β1-β2kMtα+kMtkIkt∑k′Pk′∣kSk′t-ηIktdSktdt=β1-β2kMtα+kMtkIkt∑k′Pk′∣kSk′t-ηSkt-δSkt-εkSkt∑k′Pk′∣kSk′t+Rk′tdRktdt=εkSkt∑k′Pk′∣kSk′t+Rk′t+δSkt-ηRktdMtdt=m0+ω∑k′PkSk′t-λMtwhereΛ: immigration rate, the individual enters a particular community at a constant rate;α: the half-saturation constant, which reflects the impact of media coverage during rumor spreading;η: removal rate, each type of individual moves out of the group for some reason;β1: the probability that the ignorant accepts and spreads rumors when the ignorant connects with the spreaders. β1IS indicates that some of the ignorants choose to trust the spreaders’ information and enter the group;β2: the maximum impact of media on the probability of transmission when ignorants are exposed to spreaders. β2kM/α+kM expresses that the ignorant is in contact with the spreaders; the probability of the media’s interference with the spread of rumor, with the increase of the number of spreaders, the intensities of the media, and the gradual increase in the media's interference in believing and disseminating untrue information to the individual, eventually tends to a maximum β2, where α represents saturation coefficient and measures the degree of media coverage affecting rumor spreading. Considering that media coverage can slow down the spread and diffusion of rumor, it cannot completely prevent the spreading, so we assume that β1≥β2;ε: stifling rate, the probability that a spreader transfers into a stifler when the spreader contacts with a spreader or a stifler;δ: emergency efficiency, by means of isolation and other compulsory administrative measures, the government loses the spreading ability of rumor spreaders;ω: rate of media coverage being implemented, the influence is bidirectional. While the media is affecting rumor spreading, the development of rumors in turn affects the frequency of media coverage. And the parameter ω is a proportional constant with which media coverage is being implemented based on the current rumor spread;λ: depletion rate of media coverage due to ineffective implementation;m0: density level of media coverage in this region based on the spread of rumor in another region. Even if there is no rumor case in the area, we must give warnings and nip the impossible situation in the bud.According to the previous studies,Pk′∣k=k′P(k′)/〈k〉 indicates the degree correlations so that a node of degree k is connected to a node of degree k′, and 〈k〉=∑kkP(k) is the average number of contacts each node can make. The media coverage function kM(t)/α+kM(t) is similar to the Michaelis –Menten function P1(M)=CM(t)/1+DM(t), (C,D≥0) (see [34, 35] and the references cited therein), which is a linear function of current prevalence of rumor and measures the impact of media coverage on a node in group k. It is related to the number of nodes he/she contacts and current density level of media coverage.


## 2.3. Model Analysis

This paper analyzes the mathematical model based on the stability related theory of differential equations, and due to the fact that the model depicts the dynamic propagation process of rumor in the crowd, the parameters involved are assumed to be nonnegative.For the sake of simplicity, we setΘ(t)=∑k′Pk′∣kSk′(t)=1/〈k〉∑k′k′P(k′)Sk′(t) representing the probability that a randomly selected neighbor of a given node is a spreader while Φ(t)=∑k′Pk′∣kRk′(t)=1/〈k〉∑k′k′P(k′)Rk′(t) denoting the probability that a randomly selected neighbor of a given node is a stifler. Then, system (1) can be rewritten as(2)dSktdt=β1-β2kMtα+kMtk1-Skt-RktΘt-ηSkt-δSkt-εkSktΘt+ΦtdRktdt=εkSktΘt+Φt+δSkt-ηRktdMtdt=m0+ω∑k′PkSk′t-λMtWith no loss of generality, we assume thatβ1=β2=β, then the system (2) will be transformed into the following system (3).(3)αβα+kMtk1-Skt-RktΘt-ηSkt-δSkt-εkSktΘt+Φt=0εkSktΘt+Φt+δSkt-ηRkt=0m0+ω∑k′PkSk′t-λMt=0Obviously, system (3) has a unique rumor-free equilibrium E0=0,0,…,0,m0/λ1×(2n+1).Furthermore, the existence of rumor spreading equilibrium pointE∗ will be discussed. From system (3), a simple calculation shows that(4)E∗=S1∗t,S2∗t,…,Sn∗t,R1∗t,R2∗t,…,Rn∗t,m0+ω∑kSk∗tPkλ1×2n+1(5)Sk∗t=αβηkΘHRk∗t=αβηkΘ·δ+kεΘ+ΦHM∗=m0+ω∑kSk∗tPkλwhere H=αη+kηM∗+αβkΘ·η+δ+kε(Θ+Φ).Following the recipe of van den Driessche and Watmough [34], System (4) can be written as dx/dt=F(x)-υ(x) and x=(Sk,Rk)T, where x1,x2,…,xn represent S1,S2,…,Sn, and xn+1,xn+2,…,x2n denote R1,R2,…,Rn, respectively. Applying the next generation matrix method in [34],where the rate of appearance of new spreaders is(6)Fx=αβkα+kMt1-Skt-RktΘt02n×1and the transfer rate of individuals out of compartments is(7)Vx=εkSktΘt+Φt+δSkt+ηSkt-εkSktΘt+Φt-δSkt+ηRkt2n×1the Jacobian matrices ofF(x) and υ(x) at the rumor-free equilibrium E0=0,0,…,0,m0/λ1×(2n+1) are as follows:(8)F=DFE0=F110002n×2n(9)V=DVE0=δ+η0-δηE2n×2nwhere(10)F11=αβk1α+1Mt1P11α+1Mt2P2⋯1α+1MtnPn2α+2Mt1P12α+2Mt2P2⋯2α+2MtnPn⋮⋮⋱⋮nα+nMt1P1nα+nMt2P2⋯nα+nMtnPn(11)V-1=1δ+η0δδ+η1η∙E2n×2nAccording to the concept of next generation matrix and reproduction number given in [34], the reproduction number of system (2) equals R0=ρFV-1, where ρFV-1 denotes the spectral radius of the matrix FV-1. Then,(12)R0=ρFV-1=αβλη+δk∑k=1nk2Pkαλ+km0The basic reproductive numberR0 describes the average number of new spreaders at the beginning of the spreading process. And according to [34], if R0<1, then there is no more rumor to spread in the networks; while if R0>1, then the rumor will persist in the networks.


## 3. Optimal Control Strategies of the Model

Based on the analysis of the previous section, we find that it is not enough for rumor management to discuss the dynamic process of rumor, and the optimal control techniques are of great use in optimal strategies to control various kinds of rumors. Based on the optimal control theory and numerical simulation method, we further expand the system (1) to explore the impact of the two factors such as media coverage and emergency strategies on the final spread of rumor.In real life, when rumors break out, some control measures will be taken by officials including precontrol measures (such as popularization of science education, Microblog hot tweens, press conference, and other media coverage measures) and postcontrol measures (by means of isolating and compulsive punishments of spreaders to transform them into stiflers).On one hand, when rumors prevail, the media should increase the broadcast to suppress them; on the other hand, when the rumor is controlled, the amount of media coverage should be reduced accordingly. So, we assume media coverage is a dynamically adjusted process with the density of the spreaders in the current environment (the range of rumor spreading). Meanwhile, when there are mass spreaders in the environment, a series of emergency strategies by government must be taken to reduce the harm of rumors, such as forced isolation, enhanced punishment, and cut off sources of rumors spreading. But taking into account the control costs and other factors, the emergency strategies by government can be slackened a little when the rumor situation slows down. From this perspective, emergency measures are also a dynamic adjustment process.In the actual model, we select the Lebesgue square-integrable control functionsum(t), uδ(t) to represent the control variables of media coverage and emergency strategies by government, respectively, and the system (1) can be rewritten as a new control system:(13)dIktdt=Λ-β1-β2kMtα1-umt+kMtkIkt∑k′Pk′∣kSk′t-ηIktdSktdt=β1-β2kMtα1-umt+kMtkIkt∑k′Pk′∣kSk′t-ηSkt-δ1-uδtSkt-εkSkt∑k′Pk′∣kSk′t+Rk′tdRktdt=εkSkt∑k′Pk′∣kSk′t+Rk′t+δ1-uδtSkt-ηRktdMtdt=m0+ω∑k′PkSk′t-λMtk=1,2,…,nwith initial conditions(14)Ik0=Ik0≥0,Sk0=Sk0≥0,Rk0=Rk0≥0,I10,…,In0,S10,…,Sn0,R10,…,Rn0T∈Ω~where Ω~=I1,…,In,S1,…,Sn,R1,…,RnT∈R+3n∣Ik+Sk+Rk=1,k=1,2,…,n.The objective of optimal control is to minimize the negative effects of rumor dissemination, in the meantime, and maximize the positive social utility. This paper selects the following aspects to consider the issue of social utility of rumor spreading:(1)
Maximize the number of ignorants and stiflers.(2)
Try to minimize the number of rumor spreaders.(3)
Minimize the controlling costs of variable control in the process.Now, we choose an objective function to be(15)maxJum,uδ=∫0tf∑kAkIkt+Rkt-BkSkt-Ckuk2t+uk2tdtwhere0,tf denotes the time period during a control strategy that is imposed on system (13). The parameters Ak,Bk,Ck are the weight coefficients of the nonspreading group, spreading group and process control and can be regarded as the balance coefficients between different variables in time tf, and all social utility is converted into a function with the goal of cost. um(t)=ummax, uδ(t)=uδmax represent the intensity of media coverage and emergency strategies by government, respectively; it has the greatest impact on the system. Ckum2,Ckuδ2 represent the cost of media coverage and emergency strategies by government, respectively. Since then, we can take the Lagrangian as(16)LIkt,Skt,Rkt,um,uδ=AkIkt+Rkt-BkSkt-Ckum2+uδ2The optimal control variableum(t),uδ(t) is determined by the following equation:(17)Jum∗,uδ∗=max⁡Jum,uδ∣um,uδ∈UwhereU=(um,uδ)∣measurable:0≤um≤1,0≤uδ≤1,t∈[0,tf],k=1,2,…,n is the control set, um(t),  uδ(t) are the measurable variables, and [0,tf] expresses the time period during which a control strategy is imposed on system (13), which indicates an admissible control.According to the conclusion of the existence of the optimal control solution by Fleming and Rishel [35], Theorem 1 demonstrates the existence of the optimal control solution in detail.Theorem 1.
Under the limit of (17), system (13) has an optimal control um∗,uδ∗, that is, max⁡J(um,uδ)∣um,uδ∈U=J(um∗,uδ∗).Proof.
In order to prove the conclusion of Theorem1, according to the theory of Fleming and Rishel [35] (Theorem 4.1, PP 68-69), the following conditions must be satisfied:①
The control set and corresponding state variable are not empty.②
The control set is a closed set of convexity by definition.③
The upper bound of the solution of system (13) is a linear function. Each right-hand side of the state system is continuous and is bounded by a sum of the bounded control and the state. Furthermore, it can be written as a linear function of the control variate um(t),uδ(t) with coefficients depending on time and state.④
The integral of the system (13) is concave function in the control set, c2-c1(umk+uδk) is the upper bound of the integral function, and c2>0,c1>0,k>1.
According to Theorem 9.2.1 of Lukes [36], the points on the boundary of system (13) are within the control set. It is known from definition that control set U is a closed set of convexity. According to Theorem 1, the solution of the system (13) has the supremum and satisfies condition (3); the item Ak(Ik(t)+Rk(t))-BkSk(t)-Ck(um2+uδ2) in the objective function is a concave function on the control set U, and there are some constants c2>0,c1>0,k>1 making(18)AkIkt+Rkt-BkSkt-Ckum2+uδ2≤c2-c1umk+uδk
Then system (13) satisfies all four conditions and has optimal control. That is, there exist um∗,uδ∗ making maxJ(um,uδ)∣um,uδ∈U=J(um∗,uδ∗) set up. This completes the proof.Theorem1 proves the existence of optimal control under the condition of the restrictive system (13). We discuss the necessary conditions for optimal control by means of the Pontryagin Maximum Principle [37].Pontryagin Maximum Principle requires the use of auxiliary functions to explore the optimal value problem. The differential equation system (13) will be used as the constraint condition of the objective function. After discussing the existence of the optimal control, we use Pontryagin Maximum Principle to solve the auxiliary variables of the differential equation, the corresponding boundedness condition, and the properties of the control variables um(t),uδ(t).Based on the above analysis, the Hamiltonian for the control problem is constructed as(19)HIkt,Skt,Rkt,umt,uδt,λ1kt,λ2kt,λ3kt,λ4kt=LIkt,Skt,Rkt,umt,uδt+∑k=1nλ1ktIktdt+λ2ktSktdt+λ3ktRktdt+λ4ktMtdt=LIkt,Skt,Rkt,Mt,umt,uδt+∑k=1nλ1ktΛ-β1-β2kMα1-umt+kMkIkt∑k′Pk′∣kSk′t-ηIkt+λ2ktβ1-β2kMα1-umt+kMkIkt∑k′Pk′∣kSk′t-ηSkt-δ1-uδtSkt-εkSkt∑k′Pk′∣kSk′t+Rk′t+λ3ktεkSkt∑k′Pk′∣kSk′t+Rk′t+δ1-uδtSkt-ηSkt∑k′Pk′∣k+λ4ktm0+ω∑k′PkSk′t-λMtwhereλ1k(t),λ2k(t),λ3k(t),λ4k(t) are the adjoint functions to be determined suitable, to simplify, with λ1k,λ2k,λ3k,λ4k, respectively.Next, we will derive the necessary conditions for the optimal control strategy by means of the Pontryagin Maximum Principle [37].Theorem 2.
LetIk∗(t),Sk∗(t),Rk∗(t),M∗(t) be optimal state solutions associated with the optimal control variable um∗(t),uδ∗(t) for the controlled system (13) with the corresponding initial conditions, where k=1,2,…,n. Then, there exist adjoint functions λ1k(t),λ2k(t), λ3k(t), and λ4k(t), which satisfy(20)dλ1ktdt=-Ak+λ1ktαβkΘα1-umt+kM-η-λ2ktαβkΘα1-umt+kMdλ2ktdt=Bk+λ1ktαβk2PkIktα1-umt+kMk-λ2ktαβk2PkIktα1-umt+kMk-εkΘt+Φt-εkSktkPkk-η-δ1-uδt-λ3ktεkΘt+Φt+εkSktkPkk+δ1-uδtdλ3ktdt=-Ak+λ2ktεkSktkPkk-λ3ktεkSktkPkk+ηdλ4ktdt=λ1ktαβk2IktΘα1-umt+kM2-λ4kt·λwith transversality conditions,(21)λ1ktf=λ2ktf=λ3ktf=λ4ktf=0,k=1,2,…,nFurthermore, the optimal control is given as follows:(22)um∗t=min⁡max⁡0,αβλ1k-λ2k2Ck1-umα+kM2,ummax(23)uδ∗t=min⁡max⁡0,δSktλ2k-λ3k2Ck,uδmaxand k=1,2,…,n.Proof.
To determine the adjoint equations and the transversality conditions, we use the Hamiltonian (5.6). By using the necessary condition for optimal control problems and differentiating the Hamiltonian (5.6) withIk(t), Sk(t), and Rk(t) at(24)Δ=Ikt,Skt,Rkt,Mt,umt,uδt∣Ikt=Ik∗t,Skt=Sk∗t,Rkt=Rk∗t,Mt=M∗t,umt=um∗,uδt=uδ∗we obtain(25)dλ1ktdt=-∂∂IktHIkt,Skt,Rkt,Mt,umt,uδt,λ1kt,λ2kt,λ3kt,λ4ktΔdλ2ktdt=-∂∂SktHIkt,Skt,Rkt,Mt,umt,uδt,λ1kt,λ2kt,λ3kt,λ4ktΔdλ3ktdt=-∂∂RktHIkt,Skt,Rkt,Mt,umt,uδt,λ1kt,λ2kt,λ3kt,λ4ktΔdλ4ktdt=-∂MtHIkt,Skt,Rkt,Mt,umt,uδt,λ1kt,λ2kt,λ3kt,λ4ktΔwhere k=1,2,…,n.
Hamiltonian minimizing condition: at the interior points, we have(26)∂∂umtHIkt,Skt,Rkt,Mt,umt,uδt,λ1kt,λ2kt,λ3kt,λ4ktΔ=-2Ckum+αβλ1k-λ2k1-umkα+kM2=0∂∂uδtHIkt,Skt,Rkt,Mt,umt,uδt,λ1kt,λ2kt,λ3kt,λ4ktΔ=-2Ckuδ+δSktλ2k-λ3k=0which implies(27)um∗t=0,αβλ1k-λ2k2Ck1-umα+kM2<0αβλ1k-λ2k2Ck1-umα+kM2,0≤αβλ1k-λ2k2Ck1-umα+kM2≤1ummaxαβλ1k-λ2k2Ck1-umα+kM2>1(28)uδ∗t=0,δSktλ2k-λ3k2Ck<0δSktλ2k-λ3k2Ck,0≤δSktλ2k-λ3k2Ck≤1uδmax,δSktλ2k-λ3k2Ck≥1that is,(29)um∗t=min⁡max⁡0,αβλ1k-λ2k2Ck1-umα+kM2,ummax(30)uδ∗t=min⁡max⁡0,δSktλ2k-λ3k2Ck,uδmaxk=1,2,…,n
The transversality conditions are obvious, namely,(31)λ1ktf=λ2ktf=λ3ktf=0,k=1,2,…,n.This completes the proof.Therefore, according to the analysis above, we can easily obtain the following optimality system:(32)dIktdt=Λ-αβkα1-um∗+kMkIkt∑k′Pk′∣kSk′t-ηIktdSktdt=αβkα1-um∗+kMkIkt∑k′Pk′∣kSk′t-ηSkt-δ1-uδ∗Skt-εkSkt∑k′Pk′∣kSk′t+Rk′tdRktdt=εkSkt∑k′Pk′∣kSk′t+Rk′t+δ1-uδ∗Skt-ηSktdMtdt=m0+ω∑k′PkSk′t-λMtdλ1ktdt=Ak+λ1ktαβkΘα1-um+kM-η-λ2ktαβkΘα1-um+kMdλ2ktdt=Bk+λ1ktαβk2PkIktα1-um+kMk-λ2ktαβk2PkIktα1-um+kMk-εkΘt+Φt-εkSktkPkk-η-δ1-uδ-λ3ktεkΘt+Φt+εkSktkPkk+δ1-uδdλ3ktdt=-Ak+λ2ktεkSktkPkk-λ3ktεkSktkPkk+ηwhere um(t) and uδ(t) are given by (27)-(30), respectively. The initial conditions are Ik(0)=I0, Sk(0)=S0, Rk(0)=R0, λik(tf)=0,i=1,2,3.  k=1,2,…,n. Obviously, the differential equation system satisfies the Lipschitz continuity condition, and then we get the only optimal control variable um∗(t) and uδ∗(t).


## 4. Numerical Simulations

To complement the mathematical analysis carried out in the previous section, a series of numerical simulations were performed to support and extend our theoretical results on BA scale-free networks (p(k)=2m2k-3,m=3). In the first part, we analyze the sensitivity of the basic reproductive number R0 in terms of the model parameters. During the second part, we simulated and analyzed the model involving optimal control.From Figure2, it is easy to see that the basic reproduction number R0 decreases with increasing the model parameters δ and m0, expressed as inhibition of spreading, while it increases with the increase of the parameters α and λ, expressed as promotion of spreading. Comparing the first two pictures, we can see that the influence of m0 on R0 is greater than that of δ, which shows that prior media coverage is more effective than the postcontrol from the government. From the last two pictures, it is known that the influence of α on R0 is less than that of λ, which means that we should improve the effectiveness of campaigning (i.e., lower the depletion of media coverage λ) to reduce R0.Figure 2
The relationship between the basic reproduction numberR0 and the model parameters δ, m0, α, and λ on BA scale-free network.
(a) (b) (c) (d)First, the dynamic propagation process of the rumor without the optimal control is discussed. Figure3 depicts the total densities of rumor spreaders changing over time with R0<1, while Figure 4 with R0>1. From these two figures, it can be found that when the corresponding parameter is adjusted to make R0>1, the rumor propagator tends to a stable value. The rumor continues to spread and causes bad social influence. When the corresponding parameter is adjusted to make R0<1, the rumor spreaders will be reduced to 0, and the rumor will be well controlled.Figure 3
The impact of different values ofδ,m0,α,λ on the total density of spreaders when R0<1.
(a) (b) (c) (d)Figure 4
The impact of different values ofδ,m0,α,λ on the total density of spreaders when R0>1.
(a) (b) (c) (d)The influence of different valuesδ on the total density of spreaders is shown in Figures 3 and 4. We can see that the bigger the value of δ is (the more effective of the emergency strategies control means), the lower the total density of infected nodes is. The same operation is done for the parameters m0, α, and λ. From Figures 3(b) and 4(b), one obtains that the larger the value of m0 is, the more easily the rumor could be governed, which means that the rumor control and government center had better timely report when the rumor emerges in another neighboring region owing to rapid spatial spread, although there is no spreader found in the present region. From Figures 3(c) and 4(c), we can obtain that the larger the value of α is, the more difficult the rumor could be governed. And Figures 3(d) and 4(d) illustrate how the densities of spreaders change over time with different depletion rate of media coverage λ, and it can be concluded that the bigger the λ, the larger the final densities of spreaders.In this part, based on the theory of optimal control, combined with Pontryagin Maximum Principle, we respectively simulated the dynamic changes of four groups of people in different situations, such as Figures5 and 6. The two schemes are as follows: (1) without control: um=0, uδ=0, (2) with control: um≠0, uδ≠0. The parameters are set as shown in Table 1.Table 1
Relevant parameter values.
Parameters Λ α β 1 β 2 η ε δ ω m 0 λ A k B k C k  Values 0.0001 0.5 0.133 0.133 0.001 0.000005 0.5 5 0.05 2 100 100 200Figure 5
Evolution of densities of spreaders and stiflers with control or without control.
(a) (b)Figure 6
Time evolution of control variables ((a) forum, (b) for uδ).
(a) (b)Figure5 depicts the comparison of time evolution of spreaders and stiflers which applies both media coverage and emergency strategies by government to without control. As the diagram shows, in the presence of the optimal control variables um and uδ, the peak value of spreaders is greatly reduced and the density is reduced to 0 faster than without control. Compared with the scheme without control, the time to reach the stability of the stiflers is shorter, and the final size is smaller. These phenomena show that the media coverage can make the ignorant know the rumor and reduce the infection rate of the ignorant. And the emergency strategies by government can effectively curb the spread of passion, so that the rumors will disappear as soon as possible.Figure6 shows the time evolution of the control variables um and uδ. The maximum value at the beginning of the um shows that the initial input of the media coverage is relatively large. As time evolves, the role of the media begins to play a scale effect, and the event is gradually clear and its cost decreases correspondingly. In the early stage of propagation, with the increase of spreader, the control variable uδ gradually increased to the peak and then dropped to zero. This is because the government helps spreaders understand the truth by means of forcing segregation, holding press conferences, clarifying facts, and using other means, thus preventing spreaders from spreading rumors and transforming them into stiflers. As the density of the spreaders drops, the cost of the control variable is also corresponding reduction.


## 5. Conclusion

Based on the classical SIR model, an improved SIRM rumor spreading model is proposed. On the one hand, this paper introduces the media coverage into the rumor spreading model, and with the propaganda function of the media, the rumor can recognize the reality and thus manage the rumor from the source. On the other hand, the emergency strategies by government are also used to block the diffusion of spreaders. Then the dynamic equations are obtained and the threshold of rumor spreading is calculated. At the same time, for the purpose of maximizing the social utility and minimizing the cost of rumor governance, an optimal control problem has been formulated. The existence of an optimal control has been shown, and the corresponding optimality system has been derived. The proposed optimal control scheme can achieve a low level of spreaders at a low cost. Finally, the theoretical analysis is verified by numerical simulation.



---

*Source: 4137129-2018-12-09.xml*