Discrete Dynamics in Nature and Society

Volume 2018, Article ID 9314907, 12 pages

https://doi.org/10.1155/2018/9314907

## Optimal Control of Rumor Spreading Model with Consideration of Psychological Factors and Time Delay

Business School, University of Shanghai for Science and Technology, Shanghai 200093, China

Correspondence should be addressed to Liang’an Huo; ten.haey@nilouhouh

Received 22 January 2018; Revised 1 June 2018; Accepted 5 June 2018; Published 9 July 2018

Academic Editor: Allan C. Peterson

Copyright © 2018 Liang’an Huo and Chenyang Ma. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Rumors have rapidly increasing influence on the society as well as individual life in the information age. How to control the spread of such rumors effectively has become an urgent problem to be solved. In this paper, we consider an optimal control of rumor spreading model with psychological factors and time delay. Firstly, we introduce a realistic optimal control of rumor spreading model with consideration of Holling-type II functional response and time delay. Secondly, by introducing two control strategies of both promoting scientific knowledge and releasing official information, we formulate an optimal control problem to minimize both the number of ignorant individuals and spreaders and the control cost. Thirdly, we prove the existence and the necessary conditions of optimal control strategies theoretically based on Pontryagin’s maximum principle. Our results indicate that the proposed control strategies are effective in reducing the number of spreaders and ignorant individuals and minimizing control cost.

#### 1. Introduction

With the help of network technology, rumors can widely spread in a second and lead to public opinion risks; they may lead to some unpredictable and irreparable consequences. Therefore, the impact and consequences of rumors should not be underestimated; it is a very urgent and meaningful subject to research on how to deal with rumors. Nowadays, mathematical modeling of rumor is one of the most important research areas; it has been studied by many scholars. Among them, the DK model proposed by Daley and Kendall in 1965 [1] is the early classical rumor model. In their model, the crowd was divided into three categories, the people who do not know the rumors, the people who know the rumors and spread rumors, and the people who know the rumor but do not transmit it. Zanette applied complex network theory to study the rumor spreading [2, 3]. Moreno et al. examined the dynamics of rumor spreading on scale-free networks [4]. Zhao et al. provided a good description of the amplitude effect caused by opinion leaders in the rumor spreading process [5]. Several authors considered suspected mechanism [6], incubation mechanism [7], forgetting mechanism [8], and hesitating mechanism [9] to the process of rumor spreading.

Most previous rumor spreading models mainly consider that rumor diffusion process meets the bilinear incidence rate. Rumor propagation is closely related to personal psychological quality. Thus, the bilinear incidence rate in real rumor spreading is not very appropriate. In recent years, some researchers have studied nonlinear incidence rate in rumor propagation process. For example, Huo established a rumor spreading model with consideration of Holling-type II functional response [10]. The Holling-type II functional response was used to describe the excellent explanations of the scientific knowledge. However, one significant drawback to Huo’s work was that the rumor spreading model did not consider the time delay in rumor spreading process. As is well known, individuals sometimes may not timely respond to rumors. Time delay is a common and inevitable phenomenon in nature; the development trend of the dynamic system depends on not only the present state, but also the past state. Time delays are often represented as latent period or immune period in epidemic model [11–13] and finite response period in neural network [14–16]. Similarly, time delays also exist in the rumor spreading process. For example, a person who is infected with rumor will not spread the rumor timely; it takes time for thinking and then for spreading rumor; there are little studies including time delay in rumor spreading model. Li et al. established a rumor model with time delay considering forgetting effect [17]. Li et al. proposed a time delay rumor propagation model with a saturated control function in emergencies [18, 19]. Meanwhile, rumor spreading emphasizes its timeliness; there is no doubt that inclusion of time delay in rumor spreading model makes it more realistic by describing the effects of latent period. Thus, considering the time delay when the rumor-infected individuals contact the rumor-susceptible individuals is necessary.

Scientific knowledge and official information have played a great role in the process of the rumor control. Most rumors and suspicions are caused by the lack of scientific knowledge and official information. Promoting scientific knowledge can enhance people’s ability to identify rumors indirectly; releasing official information can help people know the truth directly. For example, in March 2011, after the Tohoku earthquake, a buying frenzy for iodized salt erupted in China, people believed that iodized salt can protect people from nuclear radiation, and it has brought inestimable negative effects; however, the situation eased after the government promoting scientific knowledge and releasing official information. Therefore, this means that it is realistic to represent the actual dynamics of rumor spreading model including the effects of scientific knowledge and official information. Zhang et al. considered an 8-state ICSAR rumor propagation model with official rumor refutation [20]. Huo et al. present a way of analysis on interaction between two processes and a frame about the dissemination of scientific knowledge in order to prevent rumor transition [21]. Havakhor et al. argued that the extent of knowledge diffusion in the social media networks depends on the type of mechanism and the relative distribution of these three roles. [22]. Zhu et al. considered government intervention in the process of rumor spreading and proposed a rumor propagation model with government regulation and control [23]. Wang et al. built a knowledge transmission model by considering the self-learning mechanism to describe the dynamics of the knowledge transmission process [24]. Huo and Huang established a dynamic transmission 3D model to explore the influence of scientific knowledge in rumor spreading [25].

Recently, optimal control theory has been widely used to analyze the dynamic equation; an optimal control problem is described by suggesting an objective functional for optimizing the control cost and effect. A number of studies introduce vaccination and treatment on the spread of infectious diseases by using the control theory [26–28]. Due to the high similarity between rumor spreading model and epidemic model, various optimal control models were also proposed in rumor model. Among them, Kandhway formulated an optimal control problem, from the perspective of single campaigner, to maximize the spread of information when the campaign budget is fixed [29]. Dhar J et al. proposed media awareness as a control strategy for reducing the rumor spreading [30]. Huo et al. introduced control signals, such as science education and official media coverage attempt to convert lurkers and spreaders into stiflers [31].

Above all, in this paper, we will investigate the influence of scientific knowledge and official information on rumor spreading with psychological factors and time delay. To do this, we proposed a more realistic optimal control problem of rumor spreading model with consideration of Holling-type II functional response and time delay, in which time delay means a time that is used to think about whether to spread the rumor or not by one’s cognition or surrounding influence; two control strategies represent separately promoting scientific knowledge and releasing official information. Secondly, we use optimal control approach to minimize both the number of spreaders and ignorant individuals and the cost of two control strategies. Finally, we proved that the spreading of rumor could be controlled effectively by using the two control strategies together.

The organization of this paper is as follows. In Section 2, to prevent the spread of rumor, we introduce two control variables in rumor spreading model with time delay, and we use the optimal control techniques to find the optimal solution of the dynamics system. In Section 3, we perform some numerical simulations. In Section 4, we obtain the conclusions.

#### 2. The Model and Optimal Control Problem

In this section, we first consider a rumor spreading model with time delay. According to classic rumor model, we supposed rumor spread in a virtual community, where the population is divided into three different types: : ignorant individuals (referred to as ‘ignorant’ in this paper) who have never heard the rumor and can be easily infected by the rumor; : spreaders who have heard the rumor and passed it to others; and : removal individuals (referred to as ‘removal’ in this paper) who have heard the rumor but will not spread the rumor to others any more. They satisfy the condition .

Compared with traditional rumor spreading, we introduce the nonlinear incidence rate; we assume that rumor spreads between ignorant and spreaders, and the propagation meets Holling II functional response , in which represents the transmission rate and measures the inhibitory effect of . In the process of rumor spreading, on the one hand, due to the lack of the identification ability, one may be willing to expand the rumor; on the other hand, with the increase of spreaders, the government further intensifies the work of releasing official information; the effect of the blocking action of official information on individual’s choice whether to believe and spread the rumor is gradually increased. Thus, the contradictory psychological activity of individuals will lead to a decelerating spreading rate, which is coincident with Holling II functional response.

Time delay is introduced in the system as follows: at time only a percentage of ignorant, who leave the ignorant class time unit ago, enter the spreader class at the rate of , a Holling-type II saturated incidence function. It represents the proportion of individuals who acquire infection by contacting spreaders and become spreaders at time , where is the removal rate during thinking time ; it represents a thinking time that is used to think about whether to spread the rumor or not by one’s cognition or surrounding influence. measures the individual psychological inhibition effect from the behavioral change of ignorant when their number increases. This incidence rate is more realistic than the usual bilinear incidence rate because it includes the individual psychological factors in the process of rumor spreading.

Finally, we assume there are two ways for spreader to become removal; spreader contact the removal ; because the removal knows the rumor, spreader loses motivation to spread rumor with a certain probability .

Based on the classical D-K model, Cintron [32] proposes a deterministic dynamic model. He thinks that people are enthusiastic about passing on the word as long as it is news; once they meet with others who already know the rumor, it is no longer exciting to spread it. In this paper, we suppose that two spreaders meet, they spread information to each other. It is possible that one feels that the other person already knows the information, so they continue to spread the rumors with little motivation; only the initiating spreader loses motivation to spread rumor again and becomes a stifler; we denote by the parameter the spreader’s contact rate.

The model is formulated as the following system of delayed differential equations:where is the number of newly coming people; is the immune rate of a newly coming people; and is natural removal rate.

Based on the above analysis, we know that scientific knowledge and official information have significant effect on controlling rumor spreading. In this paper, we introduce two control strategies to decrease the number of ignorant and spreaders. Promoting scientific knowledge can help the ignorant identify information when rumor occurs; as a result, they will not believe the rumors. Therefore, we assume that ignorant become the removal directly under the action of . However, scientific knowledge plays a weak role in influencing the spreaders, because scientific knowledge is not aimed at one specific rumor, and spreader chooses to believe the rumor with preconceived notions, and it is difficult to change their thought. In order to study the subject conveniently, we assume that has no effect on spreader. Similarly, when rumor occurred, the government released official information against the rumor, which can help the spreaders learn the rumor truth; as a result, they will not spread the rumor again. Therefore, we assume that spreaders become the removal directly under the action of . However, official information has no direct effect on ignorant compared with scientific knowledge. In order to study the subject conveniently, we assume that has no effect on the ignorant.

Hence, the optimal control model becomeswith initial conditions

In order to set an optimal control problem, firstly, we make the following notational conventions. Let be a given constant and define the control set:It indicates an admissible control set; in this optimal problem, we assume a restriction on the control variable such as , because vaccination of all the susceptible individuals at one time is impossible.

The purpose of optimal control aims to spend the least control cost to reduce negative effect caused by the rumor and maximize social utility. As a result, everybody knows the truth and they all become the removal, and they will not believe and spread the rumor easily. In the model, the physical meaning of the control variable in this problem is that if the low levels of the numbers of spreaders and ignorant individuals are reached, we get more removal individuals. In case of no control, the number of spreaders and ignorant individuals increases while the number of removal individuals decreases. The prefect time for control strategies can bring the number of infected individuals down to a small level, the number of ignorant individuals begins to be built again, and more individuals recover from the infection.

Now, we consider an optimal control problem to minimize the objective (cost) functional given bywhere , denote positive weight parameters that balance the size of the term , . , are positive weight parameters which are associated with the control , , and the square of the control variable reflects the severity of the side effects. and represent maximizing the intensity of promoting scientific knowledge and releasing official information separately.

##### 2.1. Existence of an Optimal Control

The existence of the optimal control pair without delay can be obtained by using the result in Fleming and Rishel [33]. On the basis of the literature, the principle with delay given in [34] provides necessary conditions for an optimal control problem. According to these papers, the optimal control problem (2) and (3) has the Lagrangian

Theorem 1. * An optimal control pair , exists so that*

*Proof. *To prove the existence of an optimal control pair, we must check the following properties:

(I) The set of controls and corresponding state variables is nonempty.

(II) The control set is convex and closed.

(III) The right-hand side of the state system is bounded by a linear function in the state and control variables.

(IV) The integrand of the objective functional is convex on .

(V) There exist constants and such that the integrand of the objective functional satisfiesAccording to heorem 9.2.1 in Lukes [35], the existence of solution of system (2) is given with bounded coefficients, which gives condition (I). By the above definition, the control set is convex and closed, which gives condition (II). By definition, since the state system is linear in , the right side of system (3) satisfies condition (III). Using the boundedness of the solution, the integrand in the objective functional (5) is convex on (condition (IV) is proved). Next, let , , and ; we can havewhich satisfies condition (V).

We conclude that there exists an optimal control.

This completes the proof.

##### 2.2. Characterization of the Optimal Control

Next, let us derive a necessary condition for the optimal control strategies by using the Pontryagin’s maximum principle [34].

Theorem 2. *There exist adjoint variables , , that satisfywith boundary conditions**Furthermore, the optimal control variables are given as follows:*

*Proof. *To find the optimal control function for the optimal control problem, define the corresponding Hamiltonian asBy differentiating the Hamiltonian above with respect to respective states, we obtain the adjoint system as follows:and by the optimality conditions we find thatwhich impliesUsing the property of control set, we obtainso the optimal control pair is characterized as (12) and (13).

The optimality system consists of the state system coupled with the adjoint system with the initial and transversality conditions together with the characterization of the optimal control. The optimal control pair and the state are found by solving the following optimality system, which consists of the state system (2), the adjoint system (10), boundary conditions (3) and (11), and the characterization of the optimal control pair (12) and (13). This completes the proof.

#### 3. Numerical Simulations

To illustrate the theoretical analysis, in this section, we solve numerically the optimality system (2) and we present the results found. For the convenience purpose, we suppose that the number of people in the community will not increase over a short time, so we assume that , . Meanwhile, the number of people will not fall over short time; that is .

We note that the optimality system is a two-point boundary value problem, with separated boundary conditions at times t=0 and t=, according to [36], this involves use of an appropriate algorithm.

Let there exist a step size and integers with and . For reasons of programming, we consider m knots to left of 0 and right of , and we obtain the following partition:

Then, we have . Next we define the state and adjoint variables , and the controls , in terms of nodal points and .

*Algorithm*

*tep 1*

for i=,do*end* for

for ,doend for

*tep 2**end* for

*tep 3*

for i=,write*end* for.

Next, let , , separately. Furthermore, we choose the following initial values: , , . The weight constant values in the objective functional are , , , . We investigate and compare numerical results in the following three possible strategies for the control system. We also discuss the effect of different delays in control variables in** Section 3.4**.

##### 3.1. Only Promoting Scientific Knowledge Strategy ()

We suppose that the government uses only promoting scientific knowledge strategy to control rumor. With this strategy, only the control variable is used to optimize the objective function while the control variable is set to zero.

Next, let ; by calculation, it can satisfy optimality conditions. The trend of the change of , , both with and without control is shown in** Figures 1, 2**, and** 3** separately.** Figure 4** gives the corresponding optimal control variable .

**Figure 1** shows that, in the presence of a control , the number of ignorant () decreases fast compared to ‘without control’; under the effect of scientific knowledge, some ignorant become the removal.** Figure 2** shows that the number of spreaders () increases slowly with control compared to ‘without control’ at the beginning and the final size is much smaller.** Figure 3** shows that the number of the removal () increases fast with control compared to ‘without control’ and the final size is much larger. The dotted line denotes that more ignorant and spreaders become the removal when the control is implemented to the ignorant.** Figure 4** represents the optimal control variable . The curve has the maximum value in the beginning because scientific knowledge plays the most important role in the case of the high ignorant level; then it drops off because of the increasing number of spreaders and scientific knowledge has little effect on them.