ORIGINAL RESEARCH article

Front. Phys., 18 May 2023

Sec. Social Physics

Volume 11 - 2023 | https://doi.org/10.3389/fphy.2023.1194804

Dynamic analysis and optimal control of a stochastic information spreading model considering super-spreader and implicit exposer with random parametric perturbations

  • 1. School of Electronic and Information Engineering, University of Science and Technology Liaoning, Anshan, China

  • 2. School of Business Administration, University of Science and Technology Liaoning, Anshan, China

  • 3. School of Science, University of Science and Technology Liaoning, Anshan, China

Abstract

Environmental factors in social systems affect information spreading at all times. This paper proposes a stochastic S2EIR model that considers the presence of super-spreaders and implicit exposers in information spreading, as well as the stochastic perturbation of model parameters. The existence of a global positive solution using the formula is then demonstrated. Sufficient conditions for information disappearance and smooth distribution of information are calculated by using the Borel–Cantelli lemma and the strong law of large numbers. Furthermore, the optimal control strategy for the stochastic model is proposed using the Hamiltonian function. The results of the theoretical analysis are supported by numerical simulations and compared to the parameter variations of the deterministic model. The results of this study indicate that white noise facilitates the spread of information. The intensity of perturbation is proportional to the fluctuation of information spreading. Controlling random parameters can effectively facilitate the spread of information. For positive information, the randomness and complexity of the social system should be utilized to increase the spread of information. In contrast, for negative information, randomness in the social system should be suppressed to the greatest extent possible to limit information dissemination.

1 Introduction

Information is essential to the development of human society and is the primary means by which humans can comprehend society and improve their cognition. In general, the effects of various types of information on society vary. Information that is beneficial to social development, such as knowledge, innovation, and affinity, must be disseminated [1]; [2]; [3]; [4], while information that is harmful to social development, such as rumors and computer viruses, must be contained [5]; [6]; [7]. Therefore, it is valuable to investigate the characteristics and mechanisms of information dissemination.

Information dissemination is similar to the spread of infectious diseases. Daley and Kendall [8] first applied the classical infectious disease model to the spread of rumors in order to study the spread of harmful information, and they proposed the classical rumor spread DK model. The SI model [9], the SIS model [10], the SIR model [11], and other infectious disease models were then used as a research foundation to propose the SEIR model [12], the SPNR model [13], the ILSR model [14], etc.

The aforementioned models typically regard environmental factors as deterministic factors and do not account for them. However, there are many uncertain environmental factors in the real social environment. Random factors, such as government regulation of the media, sudden mass events, and outbreaks of epidemics or geological disasters, can interfere with the dissemination of information [15]. Therefore, researchers are beginning to consider the impact of uncertain environmental factors on the spread of information.

In recent years, extensive studies on stochastic factors in the field of virus propagation have been proposed. Shang [16] developed a low-dimensional system of nonlinear ordinary differential equations to model the mixed susceptible–infected (–recovered) SI(R) epidemics on a random network with general degree distributions. Zhang and Li [17] developed the stochastic SEIR model with migration and human awareness in complex networks and compared the distance between stochastic and deterministic systems at points of local disease equilibrium. Liu and Jiang [18] developed the stochastic SEIR model with asymptomatic infected individuals and calculated the specific form of the probability density near the quasi-endemic equilibrium of the stochastic system. Hussain et al. [19] proposed the stochastic SIRV model with general nonlinear incidence and vaccination and argued that stochastic fluctuations could limit the spread of disease. Wang et al. [20] designed the stochastic SICA model with standard incidence rates, provided sufficient conditions for HIV extinction and persistence, and discovered that increasing the intensity of stochastic perturbations could control the virus’s spread. Bobryk [21] believed that the spread rate in the SIR model must be positive, so he incorporated telegraphic noise, trichotomous noise, and bounded noise into the stochastic model and then discussed the effect of the stochastic perturbation on the stability behavior of disease-free equilibrium points. Kiouach and Sabbar [22] used Gaussian white noise and Lévy jumps to represent, respectively, parametric perturbations in stochastic models, and they calculated the disease’s mean persistence, ergodicity, and extinction. In addition, the application of infectious disease models with stochastic parameter perturbations to the novel coronavirus has been a hot topic over the past few years [23]; [24]; [25].

In the field of information spread, researchers have focused on the impact of stochastic factors on models of rumor spread. Cheng et al. [26] argued that individual activities play a significant role in rumor spread and that individual activities are affected by many uncertainties. As a result, they constructed the IaIdSaSdRaRd model and analyzed the dynamical behavior of the stochastic model based on this argument. Zhang and Zhu [27] concluded that the spread rate among individuals in rumor spread models is usually disturbed by random factors, and increasing the noise intensity can effectively control the spread of rumors. Tong et al. [28] constructed a stochastic rumor spread model with media coverage and age structure based on a deterministic model and confirmed the existence of a global positive solution to the model and the necessary conditions for rumor disappearance and smooth distribution. Zhou et al. [29] argued that information intervention influences rumor spread; consequently, they developed a stochastic SIRZ rumor spread model considering the information intervention mechanism and proposed an optimal control strategy for rumor control. Yue and Huo [30] developed a stochastic SICR rumor spreading model that accounted for media coverage and science education; the study’s results demonstrated that media coverage could inhibit rumor spread. Jain et al. [31] proposed the S1S2I model with expert interaction on a homogeneous network and discovered that noise perturbation caused persistent rumor spread. Huo and Dong [32] developed a stochastic ISR model with white noise media coverage, and the results of their study demonstrated that media coverage was inversely proportional to rumor spread. In addition, Mena et al. [33] analyzed the spread characteristics by developing an optimal control model based on the analysis of extinction and persistence of uncertain spread models.

The aforementioned scholars have conducted numerous studies on the influence of random factors on disease spread and rumor spread. In contrast, not all information is harmful to society in the real world, and positive information must be disseminated. Some beneficial information to society may be disseminated by individuals with authority, also known as super-spreaders. Then, it is more likely that the exposers who are in contact with the super-spreaders will comprehend and continue to spread the information. Conversely, the exposers who are not in contact with the super-spreaders are less likely to do so. The former is referred to as explicit spreaders, whereas the latter is referred to as implicit spreaders. This paper develops a stochastic S2EIR model that takes super-spreaders and exposers into account. Afterward, the existence of global positive solutions is demonstrated, suitable parameters are selected as control variables after calculating sufficient conditions for information disappearance and smooth information distribution, and the validity of the proposed theorem is validated through numerical simulation. Different from previous studies, this article considers that individuals with strong influence have played a significant role in promoting information spreading while also considering that some individuals are unwilling to spread information even if they are aware of it. Therefore, the common phenomena of “super-spreader” and “asymptomatic infection” in virus spreading are applied to the information spreading model. In addition, for each individual, the contact rate with information spreaders and the spread rate of being an information spreader may be affected by random perturbations in the environment. So, it is more practical to add random factors to the information spreading model. Meanwhile, the optimal control strategy containing random parameters of this paper is based on the optimal values calculated by the control variables.

The remainder of the paper is organized as follows. In Section 2, the stochastic S2EIR model taking super-spreaders and exposers into account is constructed. Section 3 proves the existence of a global positive solution. In Section 4, sufficient conditions for the disappearance of information are outlined. Section 5 provides sufficient conditions for stationary information distribution. The existence of optimal control and the optimal control strategy are discussed in Section 6. In Section 7, numerical simulation is used to analyze the impact of stochastic disturbance intensity on information spread and the effect of the control strategy. In the last section, conclusions are provided.

2 The model

This paper examines an open virtual community whose population size fluctuates over time t. The total population size can be represented by N(t), where all populations can be divided into five categories: 1) the receptive population S(t) who do not have access to information but are receptive to it; 2) the explicitly exposed population E1(t) that has contact with the super-spreaders; 3) the implicitly exposed population E2(t) that does not have contact with the super-spreaders; 4) the spread-information population I(t); and 5) the non-receptive population R(t) that is no longer interested in the information.

The model flow diagram is presented in Figure 1.

FIGURE 1

In the model of this paper, the implicit exposers are rather special. People who are not influenced by the information’s super-spreaders typically do not easily comprehend the information’s content since it is spread by super-spreaders. Then, these people demonstrate a lack of interest in the information. Therefore, despite being exposed to the information, these individuals choose not to spread it.

The parameters in the

S

2

EIR

model can be interpreted as follows:

  • 1) The number of individuals in a social system fluctuates over time. Therefore, in this paper, the author defines B as the number of individuals moving into the social system and defines μ as the emigration rate that is out of the social system due to force majeure factors.

  • 2) When information begins to spread through the social system, the receptive population will contact the spreaders with a probability of α. The proportion of super-spreaders is defined as m. Since super-spreaders will speed up the information spread, the spread rate of explicit exposers is defined as α(1 + m), and the spread rate of implicit exposers is represented by α.

  • 3) After the exposers make sense of the information, a portion of the explicit exposers who are interested in the information becomes the information spreaders with a probability of β. Additionally, the explicit exposers who are not interested in the information become implicit exposers with a probability of ϵ.

  • 4) Due to the different degrees of information acceptance among the exposers, the two types of exposers become the non-receptive population with probabilities of γ1 and γ2. Information is generally time-sensitive, and some spreaders lose interest in the information after a period of time. As a result, they become the non-receptive population with a probability of λ.

In addition, the uncertain factors in social systems are commonly referred to as environmental noise. It is not scientific to study the spread of information while ignoring random environmental noise fluctuations. Incorporating environmental noise into deterministic models is more representative of how information spreads in real society. The random factors added to the spread models mainly include three classical approaches: 1) Introducing Gaussian white noise into deterministic parameter perturbation models [34]. 2) Random perturbation encompassing the positive endemic equilibrium of deterministic models [35]. 3) Alternating between regimes based on the probability of Markov chains [36]. Since random perturbations in the environment may affect the contact rate of each individual with the information spreader and the spread rate of becoming an information spreader, this paper uses Gaussian white noise to generate random perturbations of α and β, and the parameters of random perturbation are expressed as follows:

Here, Wi(i = 1, 2) are independent standard Brownian motions and represent the intensities of Wi(i = 1, 2), respectively. In this paper, W1 and W2 represent the relationship without mutual influence between α and β, respectively. In other words, in real scenarios, the conversion rate of information disseminators may be influenced by individual understanding and information acceptability, but it is not affected by the natural contact rate, and there is a state of mutual independence between the two.

The stochastic perturbation parameters are introduced into the deterministic model to construct a stochastic S2EIR model driven by Gaussian white noise, and the stochastic model can be represented as follows:

3 Existence of the global and positive solution

In the rest of this paper, let be a complete probability space with a filtration satisfying the usual conditions. In addition, while contains all Pnull sets, it is increasing and right continuous [37]. It also can be denoted as follows:

The existence of the global solution is the basis of analyzing the dynamic behavior of the stochastic system (Eq. 2). At the same time, according to the actual situation, a positive value is required for the dynamic model of information transmission. The stochastic system (Eq. 2) can be proved as global and positive by Theorem 3.1.

Theorem 3.1

The existence of a unique positive solutionof the stochastic system (Eq. 2) is satisfied by any given initial value. The probability of the solution is 1 and remains in.

ProofThe existence of a unique local positive solution of the stochastic system (Eq. 2) on t ∈ [0, τe), which is based on the coefficients of the deterministic system, is locally Lipschitz continuous of any given initial value . τe is the explosion time [38]. It is required to have τe = a.s. to show this solution globally. The stopping time τ+ can be defined byLet set inf ∅ = (∅ denotes the empty set). It is easy to get τ+τe. So, if τ+ = a.s. is proved, then τe = and a.s. for all t ≥ 0. We assume that τ+ < ; then, T > 0 is existence such that P(τ+ < T) > 0. We define function V: by V(X) = InSE1E2I. Using the formula, we calculate the differential of V along the solution trajectories of the stochastic system (2). For ω ∈ (τ+ < T) and for all t ∈ [0, τe), we getPositivity of X(t) implies thatwhereSo, we haveNote that some components of X(τ+) equal 0. Thereby,Letting tτ+ in the system (Eq. 8), we haveAccording to Eqs. 89, it can be obtained that Eq. 10 is less than or equal to −. Meanwhile, for any given initial value, and S(u), E1(u), E2(u), I(u) in Eq. 10 belong to a positive invariant set and are bounded. Therefore, S(u), E1(u), E2(u), I(u) are greater than 0 and greater than −, and then, Eq. 10 is greater than −. This result is contradictory. In addition, the result obtained by Eq. 10 rejects the original hypothesis τ+ < . Thus, τ+ = .

4 Disappearance of the information

Theorems 4.1 and 4.2 give the condition for the disappearance of the information. The condition is expressed by intensities of noises and parameters of the deterministic system. In the stochastic S2EIR model built in this paper, information disappearance needs to meet the following conditions: 1) all exposers affected by the super-spreaders disappeared, and 2) all spreaders of information disappeared. If any of the two aforestated conditions is satisfied, information can disappear in the social system.

First, Theorem 4.1 gives the condition for the disappearance of information caused by the disappearance of the exposers affected by the super-spreaders.

Theorem 4.1

: For any given initial value,holds a.s. Furthermore,. Then,E1(t) tend to 0 exponentially a.s., where.

Proofuse the formula to calculate the differentiation of E1(t) in the stochastic system (2), and d ln E1(t) can be written asThus, ln E1(t) can be denoted asHere, Φ1(t) is a continuous local martingale. The quadratic variation of Φ1(t) can be denoted asBy exponential martingale inequality [38], it can be known thatwhere 0 < c < 1, k is a random integer. Using the Borel–Cantelli lemma, it is easy to know that the random integer k0(ω) exists such that for k > k0, for almost all ω ∈ Ω, . Therefore, for all t ∈ [0, k], we haveThen, it can be obtained thatnoting thatSubstituting Eq. 18 into Eq. 17, ln E1(t) can be written asHence, for k − 1 ≤ tk, can be obtained asBy the strong law of large numbers to the Brownian motion, let k and then t. It can be known that .Therefore,Finally, let c → 0. Then, can be obtained asNext, Theorem 4.2 gives the condition for the disappearance of information caused by the disappearance of the spreaders.

Theorem 4.2

For any given initial value, ,holds a.s. Furthermore,. Then,I(t) tend to 0 exponentially a.s., where.

Proofuse the formula to calculate the differentiation of I(t) in the stochastic system (Eq. 2), and d ln I(t) can be written asThus, ln I(t) can be denoted asWe denotewhere Φ2(t) is a continuous local martingale. The quadratic variation of Φ2(t) can be denoted asSimilar to Theorem 4.1, for all t ∈ [0, k], one can obtainThen, it can be obtained thatIt can be noted thatSubstituting Eq. 29 into Eq. 28, ln I(t) can be written asHence, for k − 1 ≤ tk, can be obtained asBy the strong law of large numbers to the Brownian motion, let k and then t. It can be known thatFinally, let c → 0. Then, can be obtained asRemark 4.1: andare decreasing inand. The information will disappear eventually ifandare large enough, whereand.

5 A sufficient condition for the stationary distribution

Theorem 5.1 gives the unique stationary distribution of the existence of the stochastic system (Eq. 2). This also means stability in a stochastic sense.

Theorem 5.1

If, the stochastic system (Eq. 2) with initial conditionand the following conditions are satisfied:where

Then, the stationary distribution π exists, and the solution of the stochastic system (Eq. 2) is ergodic.

By the information-existence equilibrium point, where,,, and . Then, it can be obtained that

Proofdefine a function V:whereThe differential L operator to Θ1 can be calculated asAccording to , it is easy to get thatand then, LΘ1 can be expressed aswhere and .By simple calculation, one can getDue to , it is easy to obtain thatSimilarly, the differential L operator to Θ2 can be calculated asAccording to , it is easy to get thatand then, LΘ2 can be expressed aswhere , , and .By simple calculation, one can getDue to , it is easy to obtain thatNext, the differential L operator to Θ3 can be calculated asAccording to , it is easy to get thatand LΘ3 can be obtained aswhere and I* >0.By simple calculation, one can getDue to , it is easy to obtain thatFinally, the differential L operator to Θ4 can be calculated asSubstituting Eqs 43, 48, 53, and 54 into Eq. 37, we getBy Eq. 34, the ellipsoidlies entirely in . According to [37], it is easy to know that the stochastic system (Eq. 2) has a stable stationary distribution. □

Remark 5.1: ByTheorem 5.1, there existso that the solution of the stochastic system (Eq. 2) fluctuates aroundE*. Moreover, the difference between the deterministic system and stochastic system (Eq. 2) decreases with the decrease in the values ofσ1, σ2.

6 The stochastic optimal control model

Based on the stochastic model described earlier, this paper proposes a control objective to promote the mass spread of information, considering that information, such as science, technology, and knowledge, plays a positive role in social development. Therefore, the two proportional constants α and β in the model are changed into control variables α(t) and β(t), respectively, which control the contact rate of the receptive and spread-information population and the proportion of the explicit exposers converted into the spread-information population.

Hence, the objective function can be proposed asand the objective function satisfies the state system as

The initial conditions for system (Eq. 59) are satisfied:wherewhile U is the admissible control set. 0 and tf are the time interval. The control strength and importance of control measures are expressed as c1 and c2, respectively, which are the positive weight coefficients.

Theorem 6.1

There exists an optimal control pair (α*, β*) ∈ U, so the function is established as

ProofLet and

The following five conditions must be satisfied, and then, the optimal control pair is in existence.

  • i) The set of control variables and state variables is nonempty.

  • ii) The control set U 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 U.

  • v) There exist constants d1, d2 > 0 and ρ > 1 such that the integrand of the objective functional satisfied

It is clear that conditions (i)–(iii) are established. Then, condition (iv) can be easily established such thatNext, for any t ≥ 0, there is a positive constant M, which is satisfied when |X(t)| ≤ M; therefore,Let and ρ = 2. Then, condition (v) is established. Hence, the optimal control can be realized.

Theorem 6.2

There exist adjoint variablesδ1, δ2, δ3, δ4for the optimal control pair (α*, β*) that satisfy

The boundary conditions are as follows:

In addition, the optimal control pair (α*, β*) of the state system (Eq. 59) can be given by

ProofIn order to obtain the expression of the optimal control system and optimal control pair, we define a Hamiltonian function, which can be written asAccording to the Pontryagin maximum principle, the adjoint system can be written asand the boundary conditions of the adjoint system areThe optimal control formula can be written asThen, the optimal control pair (α*, β*) can be calculated based on Eq. 73 as

Remark 6.1: So far, the optimal control system that can be obtained includes the state system (Eq. 59) with the initial conditionsS(0) = S0, E1(0) = E1,0, E2(0) = E2,0, I(0) = I0and the adjoint system (Eq. 67) with boundary conditions with the optimization conditions. The optimal control system can be written asand

7 Numerical simulations

In this section, numerical simulation is conducted using the Runge–Kutta algorithm to validate the theorem of the stochastic system (2). Since the values of the parameters are not explicitly given in most recent studies, this section determines reasonable values for the model’s parameters by combining the value range of the basic reproduction number R0 and the basic conditions proposed by the theorem.

To observe the effect of external environmental factors on information spread and the effect of stochastic perturbations on changes in population characteristics in the deterministic model, the parameter values must satisfy the fundamental condition that information can spread in the social system, i.e., the basic regeneration number R0 > 1. Therefore, the parameter value can be B = 1, α = 0.5, m = 0.3, μ = 0.3, γ1 = 0.5, γ2 = 0.7, β = 0.6, ɛ = 0.2, λ = 0.2.

First, the perturbation strength can be σ = 0.001. Figure 2 shows the probability histogram of population S(t), E1(t), E2(t), I(t). As shown in Figure 2, the probability of all populations obeying the social system is stable. Figure 3 compares the density change trend between deterministic and stochastic systems about population S(t), E1(t), E2(t), I(t) over time. As shown in Figure 3, information spreads significantly better in the system with the addition of stochastic perturbation than in the deterministic system, indicating that stochastic environmental perturbations play a positive role in information spread. However, the information is in an unstable state as it spreads through the social system, and the population density of each group constantly fluctuates over time.

FIGURE 2

FIGURE 3

Second, the perturbation strength σ is reduced to 0.0001, while all other parameters remain unchanged. Figure 4 depicts the population S(t), E1(t), E2(t), I(t) probability histogram. Compared to Figure 2, the probabilities obeyed by all populations in the social system are stable, and their distributions have converged in Figure 4. Figure 5 compares the density change trend between deterministic and stochastic systems about population S(t), E1(t), E2(t), I(t) over time. Even though the perturbation intensity decreases in Figure 5, the stochastic system has greater information spread than the deterministic system. The information spread throughout the social system is not sufficiently stable as it is subject to fluctuations. In other words, information spread in a stochastic system produces fluctuations regardless of changes in perturbation strength. Due to the addition of random environmental factors, the densities of all populations are greater than those spread in a deterministic system.

FIGURE 4

FIGURE 5

Next, to observe the effects of different perturbation strengths on information spread, the trend plots of information spread over time for the stochastic system with perturbation strengths of 0.001 and 0.0001, respectively, are combined to be analyzed. As shown in Figure 6, as the intensity of the perturbation decreases, the fluctuation of information spread tends to become smoother, indicating that information spreads more easily in systems with random environmental factors and that controlling the random factors in the system can control the fluctuation of information spread effectively.

FIGURE 6

Finally, to verify the effect of the control strategy proposed in this chapter, the author kept other parameters constant and observed the density trends of the population E1(t) and I(t) over time when the optimal control strategy was used by controlling the stochastic parameters α and β. As shown in Figure 7, when the perturbation strength σ = 0.001, the densities of populations E1(t) and I(t) with the optimal control strategy for stochastic parameters α and β are superior to those without a control strategy. As shown in Figure 8, optimal control has a positive effect on information dissemination when the perturbation intensity σ = 0.0001 is reduced further. Furthermore, change the values of other parameters and compare the densities of E1(t) and I(t) change trends with different control strategies under different intensities of perturbation. As shown in Figure 9, no matter how the parameter values change, the optimal control strategy proposed in this paper can promote information dissemination effectively. At the same time, adopting optimal control strategies can suppress the fluctuations of information during its spreading.

FIGURE 7

FIGURE 8

FIGURE 9

According to the results of theoretical analysis and numerical simulation, random factors can promote the dissemination of information. For positive information, it is necessary to fully utilize the random phenomena in the social system to spread beneficial information for social development. For negative information, it is necessary to minimize the interference of random factors on the information to suppress the dissemination of information that is harmful to social development. Technological innovation information plays an important role in the development of enterprises in practical applications [3]. It is usually necessary to leverage the randomness of the social system for technological innovation information, enhance the range of exposure to technological innovation information in society, and promote the random dissemination of technological innovation information among different individuals so as to maximize the diffusion effect. On the contrary, the spread of rumor information hinders the development of society. During the outbreak of COVID-19 all over the world, the rumor that Shuanghuanglian oral liquid could treat COVID-19 was widely spread in China [39], which led to the citizens rushing to purchase the drug and caused social panic. It can be seen from this incident that the continuous emergence of We-media and unscrupulous merchants is a random phenomenon, and effective control of these random phenomena can quickly suppress the spread of rumors.

8 Conclusion

In this paper, stochastic factors in social systems are added to the deterministic model. Then, the stochastic S2EIR model with parameter perturbations is proposed, with two parameters represented by Gaussian white noise: the contact rate of super-spreaders and the conversion rate of information spreaders. The existence of a global positive solution is demonstrated, and sufficient conditions for information disappearance and smooth distribution of information are calculated. In addition, the optimal control strategy for the stochastic model is proposed. The effect of the probability density distribution and white noise perturbation on information dissemination is validated through numerical simulation. The dissemination trends of the information with different perturbation strengths are compared.

From the research presented in this paper, the following conclusions can be drawn: 1) White noise perturbation can promote information dissemination, and random environmental factors play a positive role in information dissemination. 2) As the intensity of perturbation increases, the randomness of the model and the fluctuation of the information spread trend become more obvious. 3) Information spread can be effectively controlled by controlling the stochastic parameters, and the optimal control strategy of this paper is based on the optimal values calculated by the control variables, which is different from previous studies.

There are many uncertainties in the actual social system, and it is more scientific to incorporate these uncertainties into deterministic models in order to develop stochastic models of dissemination. In this paper, the author considers the information spread phenomena of super-spreaders and implicit exposers and uses random perturbation terms to analyze the method of controlling information spread in the actual social system. The results indicate that for positive information, the randomness and complexity of the social system should be utilized to increase the spread of information. In contrast, for negative information, randomness in the social system should be suppressed to the greatest extent possible to limit information dissemination.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

Conceptualization: SK and XH; methodology: SK and YH; software: SK and YH; validation: SK, XH, and HL; formal analysis: SK and YH; investigation: SK and YH; data curation: SK and YH; writing—original draft preparation: SK; and writing—review and editing: XH and HL. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Funding

This work was supported by the National Natural Science Foundation of China (71472080), the Humanities and Social Sciences Research Projects of the Education Department of Liaoning Province of China (2020LNJC11), and the Natural Science Foundation of Liaoning Science and Technology Agency of China (2022-MS-356).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1.

    PangGWangXWangLHaoFLinYWanPet alEfficient deep reinforcement learning-enabled recommendation. IEEE Trans Netw Sci Eng (2023) 10:87186. 10.1109/TNSE.2022.3224028

  • 2.

    ZhuHWangYYanXJinZResearch on knowledge dissemination model in the multiplex network with enterprise social media and offline transmission routes. Physica A: Stat Mech its Appl (2022) 587:126468. 10.1016/j.physa.2021.126468

  • 3.

    HuYLiuHZhaoJTuLDynamic analysis of dissemination model of innovation ability of enterprise r&d personnel. Physica A: Stat Mech its Appl (2019) 531:121743. 10.1016/j.physa.2019.121743

  • 4.

    ShangYLie algebraic discussion for affinity based information diffusion in social networks. Open Phys (2017) 15:70511. 10.1515/phys-2017-0083

  • 5.

    SunHShengYCuiQAn uncertain sir rumor spreading model. Adv Difference Equations (2021) 2021:286. 10.1186/s13662-021-03386-w

  • 6.

    ChengYHuoLZhaoLDynamical behaviors and control measures of rumor-spreading model in consideration of the infected media and time delay. Inf Sci (2021) 564:23753. 10.1016/j.ins.2021.02.047

  • 7.

    DongNPLongHVGiangNLThe fuzzy fractional siqr model of computer virus propagation in wireless sensor network using caputo atangana–baleanu derivatives. Fuzzy Sets Syst (2022) 429:2859. 10.1016/j.fss.2021.04.012

  • 8.

    DaleyDJKendallDGEpidemics and rumours. Nature (1964) 204:1118. 10.1038/2041118a0

  • 9.

    ChenLSunJGlobal stability of an si epidemic model with feedback controls. Appl Maths Lett (2014) 28:535. 10.1016/j.aml.2013.09.009

  • 10.

    CaoBShanMZhangQWangWA stochastic sis epidemic model with vaccination. Physica A: Stat Mech its Appl (2017) 486:12743. 10.1016/j.physa.2017.05.083

  • 11.

    ZamanGHan KangYJungIHStability analysis and optimal vaccination of an sir epidemic model. Biosystems (2008) 93:2409. 10.1016/j.biosystems.2008.05.004

  • 12.

    EfimovDUshirobiraROn an interval prediction of Covid-19 development based on a seir epidemic model. Annu Rev Control (2021) 51:47787. 10.1016/j.arcontrol.2021.01.006

  • 13.

    JiangGLiSLiMDynamic rumor spreading of public opinion reversal on weibo based on a two-stage spnr model. Physica A: Stat Mech its Appl (2020) 558:125005. 10.1016/j.physa.2020.125005

  • 14.

    XiaYJiangHYuZGlobal dynamics of ilsr rumor spreading model with general nonlinear spreading rate in multi-lingual environment. Chaos, Solitons and Fractals (2022) 154:111698. 10.1016/j.chaos.2021.111698

  • 15.

    LiuQJiangDDynamics of a multigroup sirs epidemic model with random perturbations and varying total population size. Commun Pure Appl Anal (2020) 19:1089110. 10.3934/cpaa.2020050

  • 16.

    ShangYMixed si (r) epidemic dynamics in random graphs with general degree distributions. Appl Maths Comput (2013) 219:50428. 10.1016/j.amc.2012.11.026

  • 17.

    ZhangYLiYEvolutionary dynamics of stochastic seir models with migration and human awareness in complex networks. Complexity (2020) 2020:115. 10.1155/2020/3768083

  • 18.

    LiuQJiangDStationary distribution and probability density for a stochastic seir-type model of coronavirus (Covid-19) with asymptomatic carriers. Chaos, Solitons and Fractals (2023) 169:113256. 10.1016/j.chaos.2023.113256

  • 19.

    HussainGKhanAZahriMZamanGStochastic permanence of an epidemic model with a saturated incidence rate. Chaos, Solitons and Fractals (2020) 139:110005. 10.1016/j.chaos.2020.110005

  • 20.

    WangXWangCWangKExtinction and persistence of a stochastic sica epidemic model with standard incidence rate for hiv transmission. Adv Difference Equations (2021) 2021:26017. 10.1186/s13662-021-03392-y

  • 21.

    BobrykRStability analysis of a sir epidemic model with random parametric perturbations. Chaos, Solitons and Fractals (2021) 143:110552. 10.1016/j.chaos.2020.110552

  • 22.

    KiouachDSabbarYDynamic characterization of a stochastic sir infectious disease model with dual perturbation. Int J Biomathematics (2021) 14:2150016. 10.1142/S1793524521500169

  • 23.

    LiYWeiZDynamics and optimal control of a stochastic coronavirus (Covid-19) epidemic model with diffusion. Nonlinear Dyn (2021) 109:91120. 10.1007/s11071-021-06998-9

  • 24.

    DananeJAllaliKHammouchZNisarKSMathematical analysis and simulation of a stochastic Covid-19 lévy jump model with isolation strategy. Results Phys (2021) 23:103994. 10.1016/j.rinp.2021.103994

  • 25.

    Aslam NoorMRazaAArifMSRafiqMSooppy NisarKKhanIet alNon-standard computational analysis of the stochastic Covid-19 pandemic model: An application of computational biology. Alexandria Eng J (2022) 61:61930. 10.1016/j.aej.2021.06.039

  • 26.

    ChengYHuoLZhaoLRumor spreading in complex networks under stochastic node activity. Physica A: Stat Mech its Appl (2020) 559:125061. 10.1016/j.physa.2020.125061

  • 27.

    ZhangYZhuJDynamics of a rumor propagation model with stochastic perturbation on homogeneous social networks. J Comput Nonlinear Dyn (2022) 17. 10.1115/1.4053269

  • 28.

    TongXJiangHChenXLiJCaoZDeterministic and stochastic evolution of rumor propagation model with media coverage and class-age-dependent education. Math Methods Appl Sci n/a (2022) 46:712539. 10.1002/mma.8959

  • 29.

    ZhouYZhangJZhuCWangHModelling and analysis of rumour propagation based on stochastic optimal control. Alexandria Eng J (2022) 61:1286980. 10.1016/j.aej.2022.06.057

  • 30.

    YueXHuoLDynamical behavior of a stochastic sicr rumor model incorporating media coverage. Front Phys (2022) 10. 10.3389/fphy.2022.1010428

  • 31.

    JainADharJGuptaVStochastic model of rumor propagation dynamics on homogeneous social network with expert interaction and fluctuations in contact transmissions. Physica A: Stat Mech its Appl (2019) 519:22736. 10.1016/j.physa.2018.11.051

  • 32.

    HuoLDongYAnalyzing the dynamics of a stochastic rumor propagation model incorporating media coverage. Math Methods Appl Sci (2020) 43:690320. 10.1002/mma.6436

  • 33.

    MenaHPfurtschellerL-MRomero-LeitonJPRandom perturbations in a mathematical model of bacterial resistance: Analysis and optimal control. Math Biosciences Eng (2020) 17:447799. 10.3934/mbe.2020247

  • 34.

    JiangDShiNA note on nonautonomous logistic equation with random perturbation. J Math Anal Appl (2005) 303:16472. 10.1016/j.jmaa.2004.08.027

  • 35.

    BerettaEKolmanovskiiVShaikhetLStability of epidemic model with time delays influenced by stochastic perturbations. Mathematics Comput Simulation (1998) 45:26977. 10.1016/S0378-4754(97)00106-7

  • 36.

    DuNKonRSatoKTakeuchiYDynamical behavior of lotka–volterra competition systems: Non-autonomous bistable case and the effect of telegraph noise. J Comput Appl Maths (2004) 170:399422. 10.1016/j.cam.2004.02.001

  • 37.

    LahrouzAOmariLExtinction and stationary distribution of a stochastic sirs epidemic model with non-linear incidence. Stat Probab Lett (2013) 83:9608. 10.1016/j.spl.2012.12.021

  • 38.

    MaoXStationary distribution of stochastic population systems. Syst Control Lett (2011) 60:398405. 10.1016/j.sysconle.2011.02.013

  • 39.

    ZhouXQiuLHaoTScir rumor propagation model with the chord mechanism in social networks. Int J Mod Phys C (2022) 33:2250014. 10.1142/S0129183122500140

Summary

Keywords

information spreading, stochastic S2EIR model, super-spreader and implicit exposer, random parametric perturbations, optimal control

Citation

Kang S, Hou X, Hu Y and Liu H (2023) Dynamic analysis and optimal control of a stochastic information spreading model considering super-spreader and implicit exposer with random parametric perturbations. Front. Phys. 11:1194804. doi: 10.3389/fphy.2023.1194804

Received

27 March 2023

Accepted

03 May 2023

Published

18 May 2023

Volume

11 - 2023

Edited by

Alessandro Vezzani, National Research Council (CNR), Italy

Reviewed by

Cheng Hu, Xinjiang University, China

Yilun Shang, Northumbria University, United Kingdom

Updates

Copyright

*Correspondence: Xilin Hou,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics