Abstract

In this paper, we formulate and investigate new switched HIV/AIDS models with drug treatment involving Caputo-fractional derivatives. Initially, due to the fractional derivative order related to the memory and hereditary effects and supposing that the model coefficients are time-varying parameters, we develop a Caputo-fractional order HIV/AIDS models with switching parameters and study their dynamics utilizing Lyapunov–Razumikhin technique. Furthermore, the results show that the fractional derivative () and the switching parameters are related to the critical threshold value ( or ) which ensures disease eradication under the condition of or . Then, a treatment compartment is introduced into the above model from the asymptomatic infected individuals until the full blown AIDS individuals. Novel sufficient conditions on the threshold value are derived to verify that the disease is eventually cleared as the critical threshold parameter is below unity. Finally, some simulations are employed to support the main results and one future research direction is presented.

1. Introduction

HIV/AIDS has been one of the most deadly diseases around the world since the first patients were documented in 1981. It has seriously affected human health and even life. According to clinical trials, the virus will go through life once infected with HIV. In the absence of drug treatment, the average life span of patients from HIV infection to AIDS is within 5–10 years [1]. With effective drug therapies, treated individuals can prolong their survival and improve their life quality but do not cure for HIV or AIDS. Thus, it is necessary to propose effective methods for prevention and control of AIDS.

In recent years, many epidemic issues have been studied by modelling mathematical models to reveal transmission of the disease or predict tendencies of the disease. In particular, Anderson et al. [2] investigated the dynamics of the initial HIV models using ordinary differential equations. The models have been improved by adding into various factors on the disease [35]. Okosun et al. [6] developed HIV/AIDS with treatment and screening of unaware infectives and investigated optimal control of the treatment. Pitchaimani and Monica [7] incorporated three time delays into a HIV-1 infection model and gave the existence of Hopf bifurcation of the model. Silva and Torres [8] introduced the fractional order into an HIV/AIDS model and surveyed local and uniform stability of its disease-free equilibrium. Huo et al. [9] added anti-HIV preventive vaccines into a fractional order HIV model and showed that the model has rich phenomenon of backward bifurcation with different dosages of vaccines and fractional derivative orders. In their models, they show that fractional order models are related to memory, history, or nonlocal property, which seems to better display complex behavior of real-world phenomena. Nowadays, fractional order models have been wide applied in many fields such as applied mathematics, engineering, economics, biology, and medicine [1012]. Modelling the dynamical behavior of the epidemic diseases by fractional derivative models has more effective than interorder modelling [1315]. For instance, Ullah et al. [16] studied the dynamics of fractional order derivative tuberculosis model and proved the existence and uniqueness of equilibria. Kheiri and Jafari [17] proposed a multipatch fractional order derivative epidemic model and derived relationship between the value of objective functional and the fractional derivative order. Khan and Atangana [18] formulated novel corona virus models with fractional derivative and studied the relationship between the fractional order parameters and the infected compartments. Naik et al. [19] modelled a fractional order HIV-1 model and verified the existence of equilibria and their asymptotical stability.

Seasonal changes in a fluctuating environment often have an effect on the transmission of an infectious disease, which naturally causes that the biological and environment parameters are varying in time. For instance, the spread of dengue disease (transmitted by Aedes aegypti) is fast in summer, and the spread is low in the winter. It has been shown that the transmission speed of childhood infections is the highest at the start of the school year and then declines gradually [20, 21]. The research works on the infection disease models affected due to seasonality (such as changes in host immunity or host behavior) have been carried out [2224]. In most of the epidemic models, the coefficients are traditionally assumed to be constant or smoothly varying functions. However, in some cases (such as vacation for school children), the models parameters are assumed to change abruptly in time [2528]. This type of parameter is called switching parameter, and the model is called the switched system. Since this type of epidemic model involves multiple subsystems and subject investigated is relatively complex, research work is relatively few [2931]. Thus, it is extremely significant to incorporate switching into epidemic models and study their dynamics.

This paper mainly develops and investigates switched HIV/AIDS models with drug treatment involving Caputo-fractional derivatives. Specifically, the models’ parameters under the influence of seasonality are supposed to be varying in time and change their function forms, and the derivative of the model is assumed to be fractional order with its memory property, which leads to a novel HIV/AIDS model. New global stability conditions of the disease-free equilibrium are derived to examine the dynamics of the disease based on Lyapunov–Razumikhin technique. Moreover, treatment strategies are proposed and evaluated to show that the drug therapy has important influence on the stability of the model.

This paper is organized as follows. The classic HIV/AIDS model is extended as a fractional switched HIV/AIDS model, and threshold conditions are obtained to illustrate that the disease-free equilibrium is globally asymptotically stable in Section 2. In Section 3, treatment strategies are applied into the above model, and sufficient conditions on the disease eradication are derived. Section 4 gives numerical examples to verify the proposed results. Conclusions are presented in Section 5.

2. Model Development

Noting that memory and hereditary properties of fractional order derivatives can make complex behavioral of epidemic models, we will extend the integer order HIV/AIDS models by introducing a fractional order HIV/AIDS model with switching parameters. Assume that the total individuals consist of the susceptible individuals , infected individuals with asymptomatic, infected individuals with symptomatic, and individuals diagnosed with AIDS , that is, . On the contrary, the variety of seasons can made biological parameters change abruptly in time. It is assumed the model coefficients are as switching parameters. Assume these parameters are controlled by a piecewise continuous (from the left) switching signal , in which is , and is the set of all switching rules. Assume that is a switched transmission rate between susceptible individuals and individuals with asymptomatic, is a switched transmission rate between the susceptible individuals and individuals with symptomatic, is a switched transmission rate between the susceptible individuals and individuals with individuals diagnosed with AIDS, is a switched transmission rate from individuals with asymptomatic moving to individuals with symptomatic, is a switched transmission rate from individuals with asymptomatic moving to individuals diagnosed with AIDS, is a switched transmission rate from individuals with symptomatic moving to individuals diagnosed with AIDS, and are mortality from diseases of individuals with asymptomatic, individuals with symptomatic, and individuals diagnosed with AIDS, respectively. denotes a recruitment rate of susceptible individuals. represents natural mortality rate of four individuals classes. Thus, the modified models is presented bywhere . Here, the initial conditions for system (1) are . All the parameters are assumed to be positive. For any values of the parameters of system (1), there exists a disease-free equilibrium , in which . In the following, we first give a positive invariant set . In this regard, many works have been done in the literatures [32, 33]. And, then, we investigate the global asymptotical stability of for system (1). Before giving the main results, we introduce the following definitions and lemmas of Riemann–Liouville and Caputo fractional derivative.

Definition 1 (see [34]). A Gamma function of is defined by

Definition 2 (see [34]). The Riemann–Liouville fractional integral of order of a function is defined bywhere is the Gamma function.

Definition 3 (see [34]). The Caputo derivative with fractional order of a function is defined by

Definition 4 (see [34]). The constant is an equilibrium point of the system , , and , if and only if .

Lemma 1 (see [34]). Let be such that and exist almost everywhere and let . Then, exists almost everywhere with

Lemma 2 (see [34]). Let and be continuous real-valued functions and nonnegative in . If is a nonnegative constant and satisfies the integral inequalitythen

Now, we give the following lemma to show the bounded of the solution for system (1):

Lemma 3. Assume that is any initial value for system (1), then the region is a set of positive invariant.

Proof. From the equations of system (1), it follows thatBased on the Laplace transform to (8), it follows thatSo, , which implies that is a set of positive invariant for system (1).
Therefore, we will derive the threshold value of system (1) by calculating the spectral radius of a next generation integral operator and investigate global stability of in the feasible region on the basis of the Lyapunov–Razumikhin method.

Theorem 1. Assume that is a solution of system (1). Suppose that andfor some , , , and . If , then the disease in system (1) dies out theoretically. In other words, the disease-free equilibrium is globally asymptotically stable.

Proof. Construct the following set of Lyapunov functions:where and .
By Lemma 1 and taking the derivative of along system (1), we haveNoting that and , it follows thatin which and .
Letting , equation (13) can be written asOn the contrary, taking , we haveCombining equations (14) and (15), it follows thatAssume that and . Taking the fractional integral on both sides of (16), for , we haveBy Lemma 2, it leads toIn addition, it follows that, for , . Similarly, it can be shown that , . Generally, for , it can be shown thatSince , it can be deducted that . Note that , and , , and converge to zero exponentially, and hence, approaches . In other words, the disease in system (1) dies out theoretically.

Remark 1. For all , , and Theorem 1 is satisfied. This is because, for , , and it follows that and , which leads to , , and converge to zero, and approaches .

Next, biological and environmental parameters are periodic variation due to the seasonal changes, and we will assume that a switching rule for system (1) is periodic and satisfies the following conditions [35]: denotes one switching period, and with and . Take be the set of periodic switching rule and . Thus, the following results are presented.

Theorem 2. Assume that is a solution of system (1). Assume that the switching rule is periodic. Suppose that andin which, for , , , and . If , then the disease in system (1) dies out theoretically. In other words, the disease-free equilibrium is globally asymptotically stable.

Proof. Assume that the switching rule is periodic. By the proof of Theorem 1, for , it follows thatin which . According to the condition of , we can get . For some integer , . The sequence , . In general, for and , it follows that , where . By the facts that , we can deduct that , , and converge to zero, and hence, approaches , which means that the disease in system (1) dies out theoretically.

Remark 2. Even though threshold value of some subsystems in system (1) is greater than one, is globally asymptotically stable as long as threshold value of system (1) is less than one.

3. Treatment Strategies

In this section, we incorporate a treatment compartment into the fractional order HIV/AIDS model with switching parameters. Assume that infected individuals with asymptomatic receive drug treatment and the proportion of effective treatment is , infected individuals with symptomatic receive drug treatment and the proportion of effective treatment is , and individuals diagnosed with AIDS receive drug treatment and the proportion of effective treatment is . Therefore, a new fractional order HIV/AIDS model with switching parameters and treatment compartment can be written as follows:

According to Definition 4, the disease-free equilibrium of system (22) is , in which . For any initial value of system (22), the region is a positive invariant set. In the following, we analysis the global asymptotical stability of for system (22).

Theorem 3. Assume that is a solution of system (22). Suppose that andfor some , and . If , then the disease in system (22) dies out theoretically. In other words, the disease-free equilibrium is globally asymptotically stable.

Proof. Construct the following set of Lyapunov functions for system (22):where and .
According to system (22), we havein which and .
Letting , we can obtainTaking and combining equations (23) and (26), we haveIt follows thatLetting and and taking the fractional integral to both sides of (28), for , it follows thatFurthermore, we deriveThen, it follows that, for , . Similarly, it can be shown that , . Generally, for , it can be shown thatSince , it can be obtained that . Note that , and hence, , , , and converge to zero exponentially, and approaches . Therefore, the disease in system (22) dies out theoretically.
In addition, we study the dynamics of system (22) when the switching rule is periodic. The following results are given.

Theorem 4. Assume that is a solution of system (22). Assume that the switching rule is periodic. Suppose that andin which, for , and . If , then the disease in system (22) dies out theoretically. In other words, the disease-free equilibrium is globally asymptotically stable.

Proof. Assume that the switching rule is periodic. By the proof of Theorem 3, it follows that, for ,in which . Noting that , we can get . For some integer , . The sequence , . In general, for and , it follows that , where . Note that , and hence, , , , and converge to zero exponentially, and approaches . Therefore, the disease in system (22) dies out theoretically.

4. Numerical Simulations

In this section, numerical examples are presented to substantiate our theoretical results introduced. Suppose that , , and one period is with and . Most of the parameter values are taken from those in [2, 7, 8, 17]. Assume that constant parameter values are and . Take the switching parameters , , , , , , , , , , , , , , and with initial values , , , and . Taking the fractional order parameter , it follows that the threshold value in the subsystem is 1.6087; the threshold value in the subsystem is . It can be calculated by Theorem 2. Figure 1 shows that the disease can be cleared under the condition of , which sustains the result given in by Theorem 2. Moreover, when taking a different fractional order , Figure 2 describes the individuals behaviors. It is observed from Figure 2 that the dynamics of the asymptomatic infected individuals and the symptomatic infected individuals increase as the fractional order increases, but the susceptible individuals decrease. Therefore, the fractional order and switching parameters have a great impact on the disease eradication.

On the contrary, take the treatment parameter , and . The other parameters are the same as Figure 1. It can be obtained that by Theorem 4. Figure 3 shows that the disease dies out when , which consists of the result given in by Theorem 4. Moreover, Figure 4 describes the individuals’ behaviors with a different fractional order . Compared with the higher values of the fractional order, it can be found that the treatment is more effective when taking small value.

5. Conclusions

Research on HIV/AIDS epidemic models involving fractional order derivative , switching parameters, and treatment compartment are becoming one of the key areas in mathematical theory of epidemiology. As we all known, the research works on the switched HIV/AIDS epidemic models are very few [3638]. This paper has studied the global asymptotical stability of new fractional order HIV/AIDS models with switching parameters and treatment compartment. By constructing Lyapunov–Razumikhin technique, new threshold values are obtained to ensure eradication of the disease. More specifically, when or , the disease-free equilibrium of system (1) is globally asymptotically stable; in other words, the disease in system (1) dies out theoretically. Moreover, it has been shown that fractional order and switching parameters are related to threshold values, which have a significant effect on extinction of disease. On the contrary, when incorporating treatment compartments into the fractional order HIV/AIDS model with switching parameters, a new HIV/AIDS model is developed. The obtained results show that the disease can be cleared theoretically when the threshold value is less than one. In addition, numerical results show when the fractional order increases, the dynamics of the asymptomatic infected individuals and the symptomatic infected individuals increase while the susceptible individuals decrease. One future work is to investigate the endemic case for fractional order switched HIV/AIDS models.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 11847081), Key Scientific Research Projects of Henan Province (Grant nos. 17A110039 and 19A110007), Fundamental Research Funds for the Henan Provincial Colleges and Universities in Henan University of Technology (Grant no. 2017QNJH18), and High-Level Personal Foundation of Henan University of Technology (no. 2017BS009).