Abstract
Outbreak and large-scale of the infectious diseases have caused enormous economic losses to all countries in the world. Constructing a network model which could reflect the transmission dynamics of the epidemics and investigating their transmission laws have a significant meaning in the precaution and control of the epidemics. In this article, a fractional-order SIS epidemic network model is proposed. First, an expression of the basic reproduction number is deduced. Second, applying the Lyapunov function, the stability of the equilibrium points about the infectious model is analyzed in detail. Finally, an example is present to verify the theoretical analysis. Furthermore, on account of the fractional-order coefficient, its influence on the transmission dynamics is also exhibited.
1. Introduction
Epidemics, especially the appearance and outbreaks of novelty infectious diseases, have turned into a major global public health problem. The destruction of natural environment and the rapid change of social structures not only provide opportunities for the spread of existing epidemics but also lead to the mutation of viruses that appeared in the past. Therefore, establishing the disease network models, revealing their transmission laws, and even predicting their transmission trends are of great importance to seek the optimal strategy for the precaution and control of the epidemics.
In the last decade, network models [1, 2], which regard individuals in a population as nodes in the networks, and the contact relationships between individuals as the edges between nodes have been studied in the fields of Internet, economy, biology, and so on. The most representative work to study the transmission dynamics of epidemics in different network structures is to establish the mathematical differential equation of infectious disease. The most influential results are the SIS (susceptible-infected-susceptible) and SIR (susceptible-infected-recovered) models through the average field theory established in [3]. Meanwhile, by analyzing the stability of the SIR network model with nonmonotonic incidence, the correlation between transmission threshold and network equilibrium point is revealed [4, 5]. People’s knowledge of disease information greatly affects the threshold of disease transmission, and media is a very effective means of information transmission [6–8]. Therefore, Zhan et al. proposed an adaptive network driven by disease information and analyzed the dynamic behavior of disease transmission when the network topology and information evolve simultaneously [9]. Li et al. [10] also studied the influence of information acquisition on the dynamics of disease transmission under a random network. The more information obtained about the disease, the smaller the number of basic regeneration and the smaller the number of patients. Some scholars even analyzed the influence of network degree distribution and node weight distribution on transmission threshold [11, 12]. These studies reveal that not only the topology of the social networks but also the information diffusion circulating among the crowds have great influence on the epidemic transmission.
It is widely known that the fractional-order differential model only needs to use a few parameters to get a better effect [13]. Since the fractional-order systems have memory function, fractional-order differential equations have been taken advantage to analyze the epidemic dynamics, which are in accordance with the transmission characteristics of infectious diseases. Based on the fractional Lyapunov stability theory, some scholars have analyzed the equilibrium point stability and transmission dynamics of the fractional-order neural networks and the fractional HIV/AIDS epidemic systems [14–16]. Guo [17] studied a fractional-order SIR model in detail, and some theoretical results about the existence and stability of the unique positive solution were obtained. Furthermore, Rostamy et al. discussed the existence of multiple equilibrium points in the SIR model, indicating that selecting appropriate fractional-order parameters can expand the stability region of the equilibrium point [18]. Moreover, when the delay parameter equaled to the critical value, Hopfield bifurcation would appear. Rihan et al. [19] analyzed the dynamic behavior of the time-delay fractional-order infectious models with nonlinear incidence rate, based on the analytical methods for integer-order time-delay systems [20, 21]. From the above research, it can be found that establishing fractional infectious network models is very necessary to analyze the transmission laws of infectious diseases.
Since the basic reproduction number is a critical parameter for a disease transmission model, stability analyses about the equilibrium points were presented in [22]. According to the basic reproduction number, the major objective of this manuscript is to analyze the stability for a fractional SIS epidemic network model at the equilibrium points, including the disease-free equilibrium and the endemic equilibrium. When is satisfied, the system is asymptotically stable at the disease-free equilibrium, which reveals that the epidemic will eventually die out. Otherwise, when , there is only one endemic equilibrium point, and a certain amount of individuals will become infectious.
2. Preliminaries
First, some definitions of the Caputo fractional differential operator are introduced:
Definition 1 (see [23]). Define the Caputo fractional derivative as follows:where .
Remark 1. Supposing the fractional-order parameter , thenUsually, the Caputo fractional derivative is abbreviated to for simplicity.
Consider the following Caputo derivative equation:where is the initial value.
Definition 2. ([24]) For system (3), if and only if , then is an equilibrium point.
For a fractional-order epidemic model, maybe there are two kinds of equilibrium points, a disease-free equilibrium point and an endemic equilibrium point, which mean that the system could achieve and maintain the states by itself or under control.
Lemma 1 (see [25]). If is a continuous differentiable function, then for ,
Lemma 2 (see [26]). If is a continuous differentiable function, thenwhere , andfor any time instant .
3. Descriptions of the Epidemic Network Model
It is well known that the spread of an infectious disease can be thought of as the spread of a virus on a social network, which depends on both person-to-person contact and the transmission power of the virus. Taking into account the effects of birth rate, death rate, and transmission rate on the transmission of infectious diseases and combining with the fractional-order characteristics, in this article, a fractional-order SIS epidemic model on a degree uncorrelated network is proposed as follows:where is the Caputo derivative, for degree, represents the relative density of the susceptible (the infected) at time , and denotes the birth rate, denotes the natural death rate, denotes the infection rate, and denotes the recovery rate, respectively. denotes the ratio that a randomly selected edge starting from a node with degree points direct to an infected node with degree . Since in a degree uncorrelated network, denotes the ratio that a randomly selected edge starting from a node with degree points direct to a node with degree , which is in proportion to ; then, and could be rewritten as , where denotes the average degree of the whole network.
Let the right side of equation of (7) be zero as follows:and define be the disease-free equilibrium with , then we can get that
Similarly, (if it exists) denotes the endemic equilibrium of system (7), which expresses that the disease would exist in crowd. Therefore, the endemic equilibrium has the formwhere .
Since the basic reproduction number is a critical parameter to analyze the existence of the endemic equilibrium point, subsequently, an expression of the basic reproduction number for system (7) will be derived.
Define a function , then
To be sure that if the function has positive solutions for , the following condition must be satisfied:
Hence, for system (7), let ; if is satisfied, there is only disease-free equilibrium. If is satisfied, there is only one endemic equilibrium.
4. Stability Analysis
Subsequently, we will discuss the stability of system (7) at the disease-free equilibrium and endemic equilibrium, respectively, which is the basis in the research field of epidemic dynamics.
The basic reproduction number is an threshold parameter to characterize the initial stage of an infectious disease. It represents the average number of people infected by an infectious person over the course of illness in a population that is all susceptible. If , the infectious cannot grow. Conversely, if , the disease can invade the population. Thus, in this study, we will analyze the stability of equilibrium based on the basic reproduction number.
4.1. Stability Analysis at the Disease-Free Equilibrium
Theorem 1. For system (7), when , the disease-free equilibrium is asymptotically stable.
Proof. For system (7), construct a Lyapunov functionApplying Lemma 2, thenSubmitting equation (7) into the inequality (14), we can obtain thatSince and , then the inequality (15) can be transformed intoAs and , the inequality (16) could be written aswith .
It is quite clear thatthenThus, choose ; if is satisfied, then , and for system (7), the disease-free equilibrium is globally asymptotically stable, which means the epidemic will be extinct at last.
4.2. Stability Analysis of the Endemic Equilibrium
For discussing the stability of the endemic equilibrium, define a Volterra-type Lyapunov function given bythen
Since equation always holds, then
Hence, the endemic equilibrium is stable.
Remark 2. It is observed that when is satisfied, the disease will disappear ultimately. When is satisfied, the disease may exist and it will not spread to all individuals.
Remark 3. In this case, for system (7), no matter or , since at the disease-free equilibrium the eigenvalues of Jacobian matrix are and , we can also derive that the disease-free equilibrium is asymptotic stable.
5. An Example
In this part, an example is given to verify the theorem. For system (7), and the parameters are valued as given in Table 1.
From Figures 1 and 2 with , we can derive that system (7) is stable, and the smaller is, the slower the disease decays. The disease-free equilibrium is , which is consistent with the simulations.
From Figure 2, we can conclude that when belongs to some range, the disease would propagate in some individuals and will not die out. For example, choose ; at this moment, , and the infectious will be convergent to 0.018 approximately, and the disease will always exist, which implied that there exists an endemic equilibrium point at this moment. Choose , then , and the disease-free equilibrium is asymptotically stable. It means that the epidemic will eventually die out, which is verified from Figure 3. Unfortunately, the convergence time is a little long.
Although we have not theoretically analyzed the effects of the fractional parameter on disease transmission, the impacts are shown in Figure 4 via simulations. We discover that the smaller the parameter is, the slower the disease converges. From Figure 4, we can find that for , means that the system becomes a integer-order system, and it converges faster than the fractional-order parameter . Compared with the theoretical results in [4, 5], stability of the equilibrium points are all analyzed and some conditions are obtained. But only in this study, the influence of the fractional-order parameter on the stability is discussed and simulated.
Remark 4. From the above analysis, not only the influence of fractional-order on the equilibrium is considered but also the influence of network topology on the stability of the equilibrium is analyzed. Compared with the existing results about the stability analysis of the infectious models, it has considered more parameters which affect the disease transmission characters, but in regard to the fractional-order parameter, it lacks the qualitative analysis, which is also the main difficulty and challenge in the field of the epidemic transmission.
6. Conclusions
Epidemic dynamics of a fractional SIS network model has been investigated in this study. Based on the basic reproduction number and the Lyapunov function, stability of the equilibrium points have been analyzed. When is satisfied, the disease-free equilibrium is asymptotically stable and the epidemic will be extinct ultimately. When is satisfied, there is only one stable endemic equilibrium. At this moment, the epidemic would exist in some amount of individuals. An example is present to verify the theoretical achievements finally. Besides, the impact of fractional parameter to transmission dynamics is also exhibited by the simulation pictures.
Noted that the network topology has a tremendous impact on the transmission dynamics of the disease, the fractional-order parameter also affects the convergence rate of the disease. Therefore, according to the distribution of the population and the existing transmission of the disease, some control strategies are proposed. For example, to control the epidemic of COVID-19, measures such as restricting access, sealing off the cities, and isolating the infective and the close contacts have been taken with a good preventive effect. Therefore, in next step, we would construct a special fractional-order model about the COVID-19 based on the real and concrete transmission data and put forward an optimal control method to make sure the epidemic die out in a finite time. The comparison of the simulation and the real data would be made, which has a positive guiding role for the prevention and control of the COVID-19.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (61775198) and the Science and Technology Project in Henan Province (202102310203).