Abstract
In this paper, we study the dynamic behavior of a stochastic tungiasis model for public health education. First, the existence and uniqueness of global positive solution of stochastic models are proved. Secondly, by constructing Lyapunov function and using It formula, sufficient conditions for disease extinction and persistence in the stochastic model are proved. Thirdly, under the condition of disease persistence, the existence and uniqueness of an ergodic stationary distribution of the model is obtained. Finally, the importance of public health education in preventing the spread of tungiasis is illustrated through the combination of theoretical results and numerical simulation.
1. Introduction
Tungiasis is a parasitic skin disease caused by female sand fleas invading the human epidermis. It is an endemic zoonosis in the tropics. The disease is mainly manifested in the feet. People walking barefoot are infected with fleas on their skin, and the body of the fleas can penetrate the skin and become parasitic [1–3]. A large number of reported infections indicate that the disease spreads in many countries in economically disadvantaged areas, such as Latin America, the Caribbean, and sub-Saharan Africa [4]. In these areas, many patients cannot walk properly due to the effects of the disease, leading to further poverty. However, it remains a neglected health problem in poor communities in countries affected by the disease [5]. This neglect has been verified due to the lack of epidemiological data in Kenya [6]. About 2.6 million Kenyans are known to have been affected by the disease, and 10 million are at risk of tungiasis [7]. Tungiasis is a common disease among people living in the mountains of Kenya, western Kenya, the Great Rift Valley, and coastal areas [2]. In 2010, the prevalence of Tungiasis among children aged 5 to 12 in the endemic area of Southern Muranga reached [8]. Dirty living environment is the main factor causing tungiasis [9]. Therefore, good personal hygiene and wearing shoes are the keys to control the spread of the disease. Thus, the introduction of public health education is essential in the preventive treatment of tungiasis. The public health department can attract people’s attention by reporting the transmission route, mode of transmission, and the change of the number of infected people to achieve the purpose of public health education. At the same time, public health programs such as free insecticide delivery, free shoe distribution, and immunization may affect the spread of the disease [10].
Mathematical models have proved to be a powerful tool for studying the spread of infectious diseases and are widely used in biology, medicine, and disease control. To control the spread of infectious diseases, scholars build mathematical models to study the dynamic behavior of infectious diseases. For example, Kahuru et al. [11] conducted stability analysis on the transmission dynamics of tungiasis in endemic areas. In [12], Lv et al. studied the dynamics and optimal control of the transmission of tungiasis. In particular, Nyang’ inja et al. [13] discussed a mathematical model of the impact of public health education on tungiasiswhere , , and , respectively, represent the number of susceptible people, public health educators, and infected people at time t. All parameters here are positive. is the constant input rate of the susceptible population. is the effective contact rate for disease transmission. is the recovery rate of the infected person. is the rate at which public health education is disseminated to susceptible populations in accordance with public health intervention strategies. is the natural mortality rate of various populations. Since the impact of public health education is not permanent and the control measures are taken gradually disappear, the educated population will be infected at a lower rate , where . represents the effect of public health education on reducing infection. is the fatality rate of the disease. The asymptotic behavior of the equilibrium point is discussed in [13].
However, the spread process of infectious diseases is bound to be affected by random disturbances, which is an essential factor to be considered in the establishment of mathematical models. Compared with deterministic models, the discussion of stochastic models is more consistent with the reality of life [14–16]. So far, the stochastic systems play an important role and have been widely used in many fields such as physics, chemistry, mechanics, economics and finance systems, aerospace engineering systems, and control theory [17–21]. For example, Wang et al. [17] developed a stochastic integral sliding mode control strategy for singularly perturbed Markov jump descriptor systems subject to nonlinear perturbation. The authors provided a solution to the stabilization problem for SPMJDSs by developing a novel stochastic integral SMC strategy. Cheng et al. [18] focused on static output feedback control for fuzzy Markovian switching singularly perturbed systems (FMSSPSs) with deception attacks and asynchronous quantized measurement output. The authors developed a more general scheme, where the partial information of FMSSPS is available in both quantizer and controller. Cheng et al. [19] analyzed the problem of hidden Markov model based control for periodic systems subject to singular perturbations and Lur’e cone-bounded nonlinearity. The highlight of their research lied that the hidden Markov model detector is forwarded to observe the fading channel mode, whose detection probabilities are generalized to be partially recognized. Inspired by [22–25], we assume that the random perturbation under the influence of white noise is proportional to the variables , and in model (1) to obtain a stochastic tungiasis model:where are the intensity of the white noise, is independent standard Brownian motion defined on a complete probability space with the filtration satisfying the usual conditions. When the noise intensity , model (2) becomes the deterministic model (1). Moreover, we assign .
This article is organized as follows. In Section 2, the existence and uniqueness of the global positive solution for the model (2) are proved. In Section 3, the sufficient conditions for disease extinction and strong persistence of the model (2) are obtained. In Section 4, the existence and uniqueness of an ergodic stationary distribution of the model (2) are obtained under the condition of disease persistence. In Section 5, the influence of random noise on disease transmission is analyzed by numerical simulation. Finally, some conclusions are given in Section 6.
2. Existence and Uniqueness of the Positive Solution
In this section, we will prove that there is a unique global positive solution to model (2).
Theorem 1. For any initial value , there exists a unique positive solution to model (2) for , and the solution will remain in with probability one.
Proof. Because the coefficients of model (2) satisfy the local Lipschitz conditions, there is a unique local positive solution on for any initial value , where is the explosion time [26]. And then, we show that a.s.. Let be sufficiently large so that lies within the interval . For , we define the stopping time asThen, we set . Clearly, is increasing when . Let , then we have a.s.. If a.s. is true, then a.s. Next, we declare that a.s. If this statement is false, for , such that . Thus, there is a positive integer such that for .
Define a -function : bywhere positive constants, and , are the undetermined coefficients. The nonnegativity of the function can be derived from for . By using It formula, we havewhereTaking , thenSo,Integrating the above equation, with both sides from 0 to and taking expectation, which can be obtained from Theorem 5.8(ii) in [27],Denote , then . For , one or more of is equal to or ; thus,Consequently,where is the indicator function of . Letting , thenwhich contradicts the hypothesis, so we get that . Therefore, the solution of the model (2) will not explode in a finite time with probability one.
3. Extinction and Persistence of the Disease
In this section, we discuss the problem of disease survival. First, we give the following lemma.
Lemma 1. Let be the solution of the model (2), then , And a.s..
Proof. The proof is similar to Lemma 2.2 in [28], so it is omitted here.
Then, we analyze the persistence and extinction of disease in model (2).
Theorem 2. If , then the disease of model (2) will be extinct at an exponential rate with probability one.
Proof. From the stochastic model (2),Calculating the integral from 0 to on both sides of (13) and dividing by ,So,From Lemma 1,From the model (2), by using It formula, we haveIntegrating (17) with both sides from 0 to , dividing by t, we can obtainCombining with (16) and Lemma 1, we havewhich implies thatThus, the disease will die out exponentially with probability one.
Theorem 3. If , then the disease will be strong persistent in the mean, andwhere is a positive constant, and it is defined in the proof.
Proof. Assignwhere are the undetermined coefficients. We use It formulawhereUsing inequality , and taking , we haveThus,Integrating (26), with both sides over the interval [0, t] and dividing by t lead toBy taking the limit inferior of both sides in (27) and combining with Lemma 1, one haswhich shows that the disease is strong persistent.
4. Existence of Ergodic Stationary Distribution
In this section, we investigate the conditions for the existence of a unique ergodic stationary distribution.
Lemma 2 (see [16]). Assume that there exists a bounded open domain with regular boundary , and(1)There is a positive number satisfied that (2)There exists a nonnegative -function V such that is negative for any Then, the Markov process has a ergodic stationary distribution , and it is unique.
Theorem 4. If , for any initial value , the solution of model (2) admits a unique stationary distribution .
Proof. First, the diffusion matrix of the model (2) is as followsLet be any bounded open domain in , that satisfiesfor all . So, condition (1) in Lemma 2 is satisfied.
And then, we are going to verify condition (2). Consider a nonnegative -function withwhere and is a sufficiently small constant and is a unique minimum value point of the -function .
Denotewhere , . are positive numbers, and we will show that later is a sufficiently small constant satisfyingwhere . is a sufficient large constant such that, , will be given as follows.
Using It formula, one haswhere . Similarly, we can getTherefore,where . And,From equations (36)–(39), we havewhereSinceFrom the above analysis, we can conclude that , . So, this tells us that condition (2) in Lemma 2 also satisfies. Hence, model (2) has a unique stationary distribution .
5. Numerical Simulations
The numerical analysis supporting the theoretical results in model (2) can be obtained by Milstein method [29]. The discretized equations are listed as follows:where the time increment and are independent distributed Gaussian random variables.
Next, some numerical simulations are given to verify our theorem.
Example 1. Consider model (2) with parameter values (some of the parameters are from [30, 31]): , and the initial point is . And, we can compute that , the conditions of Theorem 2 are satisfied. Theorem 2 tells us that diseases become extinct, see Figure 1. From Figure 1(c), we can see that compared with the deterministic model (1), random white noise will accelerate the rate of disease extinction.
(a)
(b)
(c)
Example 2. Consider model (2) with parameter values (some of the parameters are from [30, 31].): , and the initial point is . Simple calculation leads to . Therefore, the conditions of Theorem 3 hold. According to Theorem 3, the disease in model (2) has a strong persistence in the mean, see Figure 2. Comparing Figure 2(c) with Figure 2(d), we can see that the disease of model (2) goes to strong persistence in the mean. Meanwhile, under this condition, model (2) has a unique ergodic stationary distribution, see Figure 3.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
6. Conclusion
In this paper, we investigate the dynamics of a stochastic tungiasis model in public health education. Firstly, the existence and uniqueness of the global positive solution of model (2) are proved by constructing appropriate Lyapunov functions. Then, the sufficient conditions for disease extinction and persistence are analyzed. At the same time, the existence and uniqueness of the ergodic stationary distribution of model (2) are obtained under the condition of disease persistence. Through the combination of theoretical research and numerical simulation, we found that the smaller white noise can make the model (2) has a unique ergodic stationary distribution when . In comparison, the larger white noise can cause the extinction of disease in the model (2) when . Compared with the deterministic model, in the epidemic model, the combination of random white noise interference and public health education has a greater impact on the persistence and extinction of the disease, enriching the dynamic behavior of the epidemic model.
Data Availability
Data sharing is not applicable to this paper as no datasets were generated or analyzed during the current study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research work was funded by the National Natural Science Foundation of China (no.11871314), the Natural Sciences Foundation of Shanxi Province (no. 201901D111314) and the Natural Science Foundation of Datong City (no. 2020147).