Abstract
With the consideration of mechanism of prevention and control for the spread of viral diseases, in this paper, we propose two novel virus dynamics models where state feedback control strategies are introduced. The first model incorporates the density of infected cells (or free virus) as control threshold value; we analytically show the existence and orbit stability of positive periodic solution. Theoretical results imply that the density of infected cells (or free virus) can be controlled within an adequate level. The other model determines the control strategies by monitoring the density of uninfected cells when it reaches a risk threshold value. We analytically prove the existence and orbit stability of semi-trivial periodic solution, which show that the viral disease dies out. Numerical simulations are carried out to illustrate the main results.
Similar content being viewed by others
1 Introduction
Viruses are the most abundant type of biological entities and are found in almost every ecosystem on Earth, which can infect all types of life forms, from animals and plants to bacteria and archaea. Many serious diseases, such as ebola, AIDS, avian influenza, and SARS, are caused by viruses. The control and hence eradication of diseases is one of the major concerns in the study of viral epidemiology. Since viral reproduction always involves host cells and uses the cellular machinery for the synthesis of their genome and other components, there have been two methodologies used to mimic these processes by mathematical models. Perhaps the earliest and simplest classical virus dynamics model was developed by Kermack et al. in [1] and Anderson et al. in [2, 3] as follows
where \(x\) denotes the density of uninfected cells and \(y\) denotes the density of infected cells, \(\lambda \) is the rate of production of uninfected cells, \(\delta \) is their per capita death rate, \(\beta \) is the rate of infection of uninfected cells, and \(\mu \) is the rate of disappearance of infected cells. Model (1) is very important in viral epidemiology, which has been studied by many authors. A more detailed description of model (1) and its dynamic behaviors may be found in [1–3] and the references therein.
Not long after that many similar models which include the dynamics of free virus or immune response have been used to describe the short-term dynamics of virus load during drug treatment and have helped to estimate virus turnover rates in vivo. We refer to some of them in [4–13] and the references therein. However, Bartholdy et al. [14] and Wodarz et al. [15] found that the turnover of free virus is much faster than that of infected cells, which allowed them to make a quasi-steady state assumption—that is, the amount of free virus is simply proportional to the number of infected cells. Therefore, the density of infected cells \(y\) in model (1) also can be considered as a measure of the free virus load.
In recent years, many researchers studied the evolution of virus dynamics models from different perspectives. These include the existence of the threshold value, which is the index for the persistence and extinction of a viral disease; the local or global stability of virus-free equilibrium and the virus equilibrium; the existence and stability of periodic solutions, to name just a few (see, e.g., [16–25]). Particularly, Nakata [26] studied the global dynamics of a viral infection model with a latent period and nonlinear function which denotes the incidence rate of the virus infection in vivo. Wang and Zhao [27] considered the investigation of the effects of periodic drug treatment on a standard within-host virus model. Ball et al. [28] examined the dynamics of strain replacement in a simple model that includes a convex trade-off between rapid virus reproduction and long-term host cell survival. Lang and Li [29] investigated the consequences of a more general CTL response and show that a sigmoidal response function gives rise to complex behaviors previously unobserved.
It is well known that exploring an effective and easily implementable control measure to keep the level of spread of diseases is significant both theoretically and practically. For example, of many strategies, the efficient ways for elimination or control of malignant neoplasm are still chemotherapy, radiotherapy, and immunotherapy in clinical treatment in recent years. The clinical data show that the course of these treatments for some viral diseases is relatively short but the density of various cells in vivo changed radically, which in nature is submitted to short temporary effects that are negligible compared to the process duration. These short-time perturbations are often assumed to be in the form of impulses in the modeling process.
Very recently, many works have been focused on the analysis of mathematical models described by ordinary differential equations with control effects [30–33]. Particularly, the state-dependent impulsive control strategies are applied widely to the prevention of spread of infectious disease due to its economic, high efficiency, and feasibility nature. This idea can be found in many other areas like agricultural production and fishery industry; where the control measures, such as catching, poisoning, releasing the natural enemy, and harvesting, are taken only when the number of population reaches an economic threshold value. Following this idea, the dynamic behaviors of population, epidemic and language models with state-dependent impulsive effects are considered, and the existence and stability of positive periodic solution by the Poincaré map, properties of the Lambert \(W\) function, and analog of Poincaré criterion are obtained [34–38].
Motivated by the above works, we propose, in this paper, two novel virus dynamics models with state feedback control strategies. The main purpose is to investigate the state feedback control strategies which govern whether the viral disease dies out or not, and further to examine how the state feedback control strategies affect the prevention and control of viral disease.
This paper is organized as follows. In the next section, we introduce two novel virus dynamics models where the state feedback control strategies are considered, some basic definitions, and an important lemma. The density of infected cells (or free virus) as control threshold value, and some sufficient conditions are presented in Sect. 3 for the existence and orbital stability of positive periodic solutions. In Sect. 4, the existence and stability of semi-trivial periodic solution are obtained, where the density of uninfected cells is the control threshold value. The numerical simulations are carried out in Sect. 5 for illustration. Some concluding remarks are presented in Sect. 6.
2 Model formulation and preliminaries
The basic reproductive rate for model (1) is \(R_0=\frac{\lambda \beta }{\mu \delta }\), which is defined as the average number of secondary infected cells generated by a single infected cell placed in an uninfected cell (or free virus). It is clear that (i) if \(R_0<1\), then model (1) admits only a global asymptotical stable virus-free equilibrium \((\frac{\lambda }{\delta },0)\); (ii) if \(R_0>1\), then model (1) admits an unstable trivial equilibrium \((\frac{\lambda }{\delta },0)\) and a unique globally asymptotically stable positive equilibrium \((x^*,y^*)\), where
Generally, the objective of treatment and control for viral disease is to reduce the density of free virus or infected cells in vivo. For all this, we propose, firstly, a novel virus dynamics model with state feedback control strategies, while regarding the density of free virus or infected cells as control threshold value. The control model is governed by the following state feedback impulsive differential equation:
The meaning of model (3) is as follows: when the density of infected cells (or free virus) \(y\) reaches the critical threshold value \(Y_c\) at time \(t_i(Y_c)\) at the \(i\)-th time, control strategies (such as chemotherapy, radiotherapy, immunotherapy, etc.) are taken and the densities of uninfected cells \(x\) and infected cells (or free virus) \(y\) immediately become \((1-p)x(t_i(Y_c))+\tau \) and \((1-q)y(t_i(Y_c))\), respectively, where \(p\), \(q\in (0,1)\), and \(\tau >0\).
Remark 1
It is clear that a priori time of control strategies depends on the density of infected cells (or free virus), which makes our control strategy a “state feedback control.”
On the other hand, from the dependencies between the densities of uninfected cells and infected cells (or free virus) in model (1), it is obvious that if the density of uninfected cells is high then infected cells (or free virus) are growing fast. Therefore, regarding the density of uninfected cells as control threshold value, we propose another novel model simulating immune boosting and viral suppressing as state feedback control strategies, and investigate its dynamical behaviors. We suppose that, therefore, when the density of uninfected cells \(x\) reaches the risk threshold value \(X_{r}\) at time \(t_i(X_r)\) at the \(i\)-th time, the control strategies are taken and the densities of uninfected cells \(x\) and infected cells (or free virus) \(y\) turn very suddenly to a great degree to \((1-p)x(t_i(X_r))+\tau \) and \((1-q)y(t_i(X_r))\), respectively, where \(p\), \(q\in (0,1)\), \(\tau >0\), and \((1-p)X_r+\tau <X_r\).
Under the assumption aforementioned, we come to a control model that is governed by the following ordinary differential equation with state feedback control strategies:
We assume, throughout this paper, that \(R_0=\frac{\lambda \beta }{\mu \delta }>1\). By the biological background, we only consider models (3) and (4) in the region \(\mathbb {R}_+^2=\{(x,y): x> 0, y> 0\}\), where the biology makes sense. Obviously, \(\mathbb {R}_+^2\) is divided into four domains with vertical isocline \(\frac{\mathrm {d}x}{\mathrm {d}t}=0\) and horizontal isocline \(\frac{\mathrm {d}y}{\mathrm {d}t}=0\), followed by
The global existence and uniqueness of solution for models (3) and (4) are guaranteed by the smoothness of the right-hand sides of models (3) and (4). For more details, we refer to [39].
Lemma 1
Solutions of models (3) and (4) with the initial value in the interior of \(\mathbb {R}_+^2\) at time \(t=t_0\ge 0\) are positive.
The proof of Lemma 1 is obvious, hence we omit it here.
Let \(\mathbb {S}\subset \mathbb {R}^2_+\) be a nonempty set and \(P_0\in \mathbb {R}^2_+\) be a point. The distance between \(P_0\) and \(\mathbb {S}\) is defined by \(\rho (P_0, \mathbb {S})=\inf _{P\in \mathbb {S}}|P-P_0|\). Let \(z(t)=(x(t),y(t))\) be a solution of model (3) starting from initial point \(z_0\in \mathbb {R}_+^2\) at \(t=t_0\). We define the positive orbit as follows:
Definition 1
(Orbital stability [40]) A trajectory \(O^+(z_0,t_0)\) is said to be orbitally stable if for any given \(\varepsilon >0\), there is a constant \(\delta \,{=}\,\delta (\varepsilon )>0\) such that for any other solution \(z^*(t)\) of model (3), \(\rho (z^*(t),O^+(z_0,t_0))\,{<}\,\varepsilon \) for all \(t>t_0\) when \(\rho (z^*(t_0),O^+(z_0,t_0))<\delta \).
Definition 2
(Orbitally asymptotical stability [40]) A trajectory \(O^+(z_0,t_0)\) is said to be orbitally asymptotically stable if it is orbitally stable, and there exists a constant \(\eta >0\) such that for any other solution \(z^*(t)\) of model (3), \(\lim _{t\rightarrow \infty }\rho (z^*(t),O^+(z_0,t_0))=0\) when \(\rho (z^*(t_0),O^+(z_0,t_0))<\eta \).
Next, for model (3), let
and
For any point \(A_0(x_0,Y_c)\in {\varGamma }_2\), consider \(O^+(A_0,t_0)\) starting from point \(A_0\) at time \(t_0\), then \(O^+(A_0,t_0)\) jumps to point \(A_0^+((1-p)x_0+\tau ,(1-q)Y_c)\) on section \({\varGamma }_1\) at \(t=t_0^+\) due to control effects \(\triangle x(t)=-px(t)+\tau \) and \(\triangle y(t)=-qY_c\). Then, trajectory \(O^+(A_0,t_0)\) intersects section \({\varGamma }_2\) at point \(A_1(x_1,Y_c)\) again. Therefore, we can define a Poincaré map on section \({\varGamma }_2\) as follows:
Furthermore, we can define \(x_2=F_1(x_1)=F_1^2(x_0)\), \(\ldots \), \(x_k=F_1(x_{k-1})=F_1(F_1^{k-1}(x_0))=F_1^k(x_0)\), \(\ldots \).
Definition 3
A trajectory \(O^+(A_0,t_0)\) of model (3) is said to be order-\(k\) periodic if there exists a positive integer \(k\ge 1\) such that \(k\) is the smallest integer for \(F_1^k(x_0)=x_0\).
Similarly, we can define two sections to model (4) by
and a Poincaré map on section \({\varGamma }_4\)
Next, we consider the following autonomous system with pulse effects:
where \(f\) and \(g\) are continuous differentiable functions defined on \(\mathbb {R}^2\) and \(\varphi \) is a sufficiently smooth function with \(\mathrm {grad}\varphi \ne 0\). Let \((\mu (t),\nu (t))\) be a positive \(T\)-periodic solution of system (8). The following result comes from Corollary 2 of Theorem 1 of [41].
Lemma 2
(Analog of Poincaré criterion) If the Floquet multiplier \(\mu \) satisfies \(|\mu |<1\), where
with
and \(f\), \(g\), \(\frac{\partial \xi }{\partial x}\), \(\frac{\partial \xi }{\partial y}\), \(\frac{\partial \eta }{\partial x}\), \(\frac{\partial \eta }{\partial y}\), \(\frac{\partial \varphi }{\partial x}\), and \(\frac{\partial \varphi }{\partial y}\) have been calculated at the point \((\mu (\tau _j),\nu (\tau _j))\), \(f_+=f(\mu (\tau _j^+),\nu (\tau _j^+))\), \(g_+=g(\mu (\tau _j^+),\nu (\tau _j^+))\), and \(\tau _j\) \((j\in N)\) is the time of the \(j\)-th jump, then \((\mu ,\nu )\) is orbitally asymptotically stable.
3 Analysis of the virus dynamics model (3)
Since virus equilibrium \((x^*,y^*)\) is globally asymptotically stable for model (1), any solution of model (3) without state feedback control strategy will eventually tend to \((x^*,y^*)\). If \(Y_c\le y^*\), we easily see that trajectory of model (3) with initial value \((x_0,y_0)\in {\varGamma }_1\) will intersect section \({\varGamma }_2\) infinitely many times, and when \(Y_c> y^*\), then trajectory of model (3) starting from point \((x_0,y_0)\in {\varGamma }_1\) may intersect section \({\varGamma }_2\) finitely many times. Therefore, in this section, we give some sufficient conditions for the existence and stability of positive periodic solutions in two cases of \(Y_c\le y^*\) and \(Y_c> y^*\), respectively.
3.1 Case \(Y_c\le y^*\)
The following result is on the existence of positive order-1 periodic solution.
Theorem 1
For any \(p\), \(q\in (0,1)\), and \(\tau >0\), model (3) admits a positive order-1 periodic solution.
Proof
For point \(A_0^+(x_0,(1-q)Y_c)\in {\varGamma }_1\) in domain \(I\!I\) with \(x_0\le (1-p)\frac{\mu }{\beta }\), in view of the properties (5) of the phase space of model (3), trajectory \(O^+(A_0^+,t_0)\) of model (3) will enter into domain \(I\!I\!I\) or in succession enter into domain \(I\!V\) and finally intersect \({\varGamma }_2\) at point \(A_1(x_1,Y_c)\). Therefore, we have \(x_1>\frac{\mu }{\beta }\). At point \(A_1\), trajectory \(O^+(A_0^+,t_0)\) jumps to point \(A_1^+((1-p)x_1+\tau ,(1-q)Y_c)\) on section \({\varGamma }_1\) due to state feedback control strategies. Furthermore, trajectory \(O^+(A_0^+,t_0)\) intersects section \({\varGamma }_2\) at point \(A_2(x_2,Y_c)\).
From the facts \(x_0\le (1-p)\frac{\mu }{\beta }\) and \(x_1>\frac{\mu }{\beta }\), then \(x_0<(1-p)x_1+\tau \). It follows that point \(A_0^+\) is left point \(A_1^+\). We claim that point \(A_1\) is right point \(A_2\). In fact, if point \(A_1\) is left point \(A_2\) or the two points coincide, then orbits \(\widetilde{A_0^+A_1}\) and \(\widetilde{A_1^+A_2}\) must intersect at a point \((x^0,y^0)\). This shows that there are two different solutions which start from this point. It contradicts with the uniqueness of solution for model (3). Therefore, by Poincaré map (6), it follows that
On the other hand, for the intersection \(B_0^+(\frac{\mu }{\beta },(1-q)Y_c)\) of line \(L_1:\beta x-\mu =0\) and section \({\varGamma }_1\), trajectory \(O^+(B_0^+,t_0)\) intersects section \({\varGamma }_2\) at point \(B_1(\hat{x}_1,Y_c)\) and then jumps to point \(B_1^+((1-p)\hat{x}_1+\tau ,(1-q)Y_c)\) on section \({\varGamma }_1\) and finally reaches point \(B_2(\hat{x}_2,Y_c)\) on section \({\varGamma }_2\) again. If there is a positive constant \(\tau ^*\) such that \((1-p)\hat{x}_1+\tau ^*=\frac{\mu }{\beta }\), then \(B_1^+\) coincides with point \(B_0^+\) for \(\tau =\tau ^*\); that is, point \(B_1\) coincides with point \(B_2\). Otherwise, point \(B_1^+\) is left point \(B_0^+\) due to \((1-p)\hat{x}_1+\tau <\frac{\mu }{\beta }\) for \(\tau \in (0,\tau ^*)\) and is right point \(B_0^+\) due to \((1-p)\hat{x}_1+\tau >\frac{\mu }{\beta }\) for \(\tau >\tau ^*\). Further, from the properties (5) of phase space of model (3), point \(B_2\) is right point \(B_1\) for any \(\tau \in (0,\tau ^*)\cup (\tau ^*,+\infty )\). Namely, \(\hat{x}_2>\hat{x}_1\) for this case.
Thus, from the above discussion, we get that
-
(i)
if \(\hat{x}_1=\hat{x}_2\), then model (3) has a positive order-1 periodic solution;
-
(ii)
if \(\hat{x}_2>\hat{x}_1\), then
$$\begin{aligned} \hat{x}_2-\hat{x}_1=F_1(\hat{x}_1)-\hat{x}_1>0. \end{aligned}$$(11)By (10) and (11), it follows that Poincaré map (6) has a fixed point. This amounts to saying that model (3) has a positive order-1 periodic solution. This completes the proof.
Let \((\phi (t),\psi (t))\) be a positive order-1 periodic solution of model (3) with period \(T\). On the orbital stability of solution \((\phi (t),\psi (t))\) of model (3), we have the following Theorem 2.
Theorem 2
If
then \((\phi (t),\psi (t))\) is orbitally asymptotically stable.
Proof
Suppose that \((\phi (t),\psi (t))\) intersects the sections \({\varGamma }_1\) and \({\varGamma }_2\) at points \(C^+((1-p)\phi (T)+\tau ,(1-q)Y_c)\) and \(C(\phi (T),Y_c)\), respectively. Comparing with system (8), we have
\(\xi (x,y)=-px+\tau \), \(\eta (x,y)=-qy\), \(\varphi (x,y)=y-Y_c\), \((\phi (T),\psi (T))=(\phi (T),Y_c)\), and \((\phi (T^+),\psi (T^+))=((1-p)\phi (T)+\tau ,(1-q)Y_c)\). Thus,
and
Furthermore, it follows from (13), (14), and (9) that
and
On the other hand, integrating both sides of the second equation of model (3) along the orbit \(\widetilde{C^+C}\) gives
By condition (12), we know that model (3) satisfies all conditions of Lemma 2. It then follows from Lemma 2 that the order-1 periodic solution \((\phi (t),\psi (t))\) of model (3) is orbitally asymptotically stable. This completes the proof.
Remark 2
Generally, the condition (12) of Theorem 2 is not easy to test since the expression of periodic solution \((\phi (t),\psi (t))\) is unknown. We note that, however, it is also weaker since the exponent term of condition (12) is less than 1.
The following corollary is a direct consequence of Theorem 2.
Corollary 1
Let \((\phi (t),\psi (t))\) be a positive order-1 periodic solution of model (3) with periodic \(T\). If
then \((\phi (t),\psi (t))\) is orbitally asymptotically stable.
3.2 Case \(Y_c>y^*\)
As already mentioned above, this case \(Y_c>y^*\) is complex since the virus equilibrium \((x^*,y^*)\) of model (1) is globally asymptotically stable. Additionally, since control measures are aimed at reducing the density of infected cells (or free virus) \(y\), we suppose that \((1-q)Y_c<y^*\). We now have the following result on the existence and orbital stability of positive order-k periodic solution of model (3) with \(Y_c>y^*\).
Theorem 3
For any \(p\), \(q\in (0,1)\), and \(\tau >0\), one of the following statements is valid.
-
(a)
If there is a positive constant \(\hat{x}=\hat{x}(Y_c)\in (0,x^*]\) such that trajectory \(O^+(E_0,t_0)\) of model (3) starting from point \(E_0(\hat{x},(1-q)Y_c)\) is tangent to line \(L: y=Y_c\) at point \((\frac{\mu }{\beta },Y_c)\) and \((1-p)\frac{\lambda }{\delta }+\tau <\hat{x}\), then model (3) has a positive order-1 or order-2 periodic solution, which is orbitally asymptotically stable.
-
(b)
Suppose that for any \(x\in (0,x^*]\), trajectory \(O^+(A_0,t_0)\) of model (3) starting from the initial point \(A_0(x_0,(1-q)Y_c)\) cuts line \(L: y=Y_c\) at point \(B(\hat{x}_0,Y_c)\), where \(\hat{x}_0\ge x^*\). Further, if \((1-p)\frac{\lambda }{\delta }+\tau <x^*\), then model (3) has a positive order-1 or order-2 periodic solution, which is orbitally asymptotically stable.
-
(c)
Suppose that for any \(x\in (0,x^*]\), trajectory \(O^+(F,t_0)\) of model (3) starting from point \(F(x,(1-q)Y_c)\) does not intersect with line \(L:y=Y_c\). Further, if \((1-p)\frac{\lambda }{\delta }+\tau <x^*\), then model (3) has no positive order-\(k\) \((k\ge 1)\) periodic solution and the virus equilibrium \((x^*,y^*)\) is globally asymptotically stable.
Proof
We first prove conclusion (a). Suppose that there is a positive constant \(\hat{x}=\hat{x}(Y_c)\in (0,x^*]\) such that trajectory \(O^+(E_0,t_0)\) of model (3) starting from point \(E_0(\hat{x},(1-q)Y_c)\) crosses section \({\varGamma }_1\) at point \(E_{0p}(\hat{x}_1,(1-q)Y_c)\) and is tangent to line \(y=Y_c\) at point \(E_1(\frac{\mu }{\beta },Y_c)\). Then trajectory starting from point \((x,(1-q)Y_c)\) with \(x\in (\hat{x},x^*)\) will tend to virus equilibrium \((x^*,y^*)\) and not intersect with section \({\varGamma }_2\). Moreover, since \((1-p)\frac{\lambda }{\delta }+\tau <\hat{x}\), then \((1-p)x+\tau <\hat{x}\) for any point \((x,Y_c)\in {\varGamma }_2\) and \(x<\frac{\lambda }{\delta }\). Therefore, for any \(E(x,Y_c)\in {\varGamma }_2\) and \(x<\frac{\lambda }{\delta }\), trajectory \(O^+(E,t_0)\) starting from point \(E(x,Y_c)\) will intersect with \({\varGamma }_2\) infinitely many times due to state feedback control strategies \(\triangle x=-px(t)+\tau \) and \(\triangle y(t)=-qY_c\).
Considering any two points \(D_i(x_i,Y_c)\) and \(D_j(x_j,Y_c)\) on section \({\varGamma }_2\), where \(x_i\), \(x_j\in (0,\frac{\delta }{\lambda })\), and \(x_i<x_j\), it is obviously that point \(D_i^+((1-p)x_i+\tau ,(1-q)Y_c)\) is left point \(D_j^+((1-p)x_j+\tau ,(1-q)Y_c)\) in view of control strategies \((1-p)x_i+\tau <(1-p)x_j+\tau <\hat{x}\). Further, trajectories \(O^+(D_i,t_0)\) and \(O^+(D_j,t_0))\) intersect section \({\varGamma }_2\) at points \(D_{i+1}(x_{i+1},Y_c)\) and \(D_{j+1}(x_{j+1},Y_c)\), respectively. We claim that
In fact, if (18) is false, that is \(x_{j+1}\ge x_{i+1}\), then point \(D_{j+1}(x_{j+1},Y_c)\) is right point \(D_{i+1}(x_{i+1},Y_c)\), or the two points coincide. So it follows that trajectories \(O^+(D_i,t_0)\) and \(O^+(D_j,t_0)\) intersect at a point \((\hat{x},\hat{y})\). This implies that there are two different solutions which start from the point \((\hat{x},\hat{y})\), which contradicts the uniqueness of solution. Inequality (18) is thus valid.
Suppose that trajectory \(O^+(C_0,t_0)\) of model (3) starting from point \(C_0(x_0,Y_c)\) (\(x_0<\frac{\lambda }{\delta }\)) on section \({\varGamma }_2\) jumps to point \(C_0^+((1-p)x_0+\tau ,(1-q)Y_c)\) on section \({\varGamma }_1\), and reaches section \({\varGamma }_2\) at point \(C_1(x_1,Y_c)\) again due to the fact \((1-p)\frac{\lambda }{\delta }+\tau <\tilde{x}\), where \(x_1\in (x^*,\frac{\lambda }{\delta })\) and then jumps to point \(C_1^+((1-p)x_1+\tau ,Y_c)\) at section \({\varGamma }_1\). At state \(C_1^+\), trajectory \(O^+(C_0,t_0)\) intersects section \({\varGamma }_2\) at point \(C_2(x_2,Y_c)\), where \(x_2\in (x^*,\frac{\lambda }{\delta })\). By the Poincaré map (6) of section \({\varGamma }_2\), it follows that \(x_1=F_1(x_0)\) and \(x_2=F_1^2(x_0)\). Repeating the above process, we have \(x_{n+1}=F_1^n(x_0)\) \((n=0,1,\ldots )\). Particularly, model (3) has a positive order-1 periodic solution when \(x_0=x_1\) and a positive order-2 periodic solution when \(x_0\ne x_1\) and \(x_0=x_2\).
Now, we consider the general situation for \(x_0\ne x_1\ne x_2\ne \cdots \ne x_{k}\) \((k>2)\). On the relation of \(x_0\), \(x_1\), and \(x_2\), there are the following two different cases.
-
(C1)
\(x_0<x_1\) In this case, we can get from (18) that \(x_1>x_2\). This results in the relation of \(x_0\), \(x_1\), and \(x_2\) to be one of the following two cases.
-
(i)
\(x_2<x_0<x_1\) In this case, \(x_3>x_1>x_2\) by (18). Repeating the above process, we have
$$\begin{aligned} x^*&< \dots <x_{2k}<\dots <x_0< x_1<\dots <x_{2k+1}\\&< \dots <\frac{\lambda }{\delta }. \end{aligned}$$ -
(ii)
\(x_0<x_2<x_1\) Similar to (i), we have
$$\begin{aligned} x_0&< x_2<\dots <x_{2k}<\cdots <x_{2k+1}<\dots \\&< x_1<\frac{\lambda }{\delta }. \end{aligned}$$ -
(C2)
\(x_0>x_1\). In this case, it follows from (18) that \(x_1<x_2\). This results in the relation of \(x_0\), \(x_1\), and \(x_2\) to be one of two cases.
-
(iii)
\(x_1<x_0<x_2\) In this case, \(x_2>x_1>x_3\) by (18). Repeating the above process, we have
$$\begin{aligned} x^*&< \dots <x_{2k+1}<\dots <x_0<\dots <x_{2k}\\&< \dots <\frac{\lambda }{\delta }. \end{aligned}$$ -
(iv)
\(x_1<x_2<x_0\) Similar to (iii), we have
$$\begin{aligned} x^*&< x_1<\dots <x_{2k+1}<\cdots <x_{2k}\\&< \dots <x_2<x_0. \end{aligned}$$
-
(i)
Further, in (i) of case (C1), we have \(\lim _{k\rightarrow \infty }x_{2k}=\theta _*\) and \(\lim _{k\rightarrow \infty }x_{2k+1}=\theta ^*\), where \(x^*<\theta _*<\theta ^*<\frac{d}{\lambda }\). Hence \(\theta _*=F_1(\theta ^*)\) and \(\theta ^*=F_1(\theta _*)\). So model (3) has an orbitally asymptotically stable positive order-2 periodic solution. Similarly, in (ii) of case (C1) and (iv) of case (C2), model (3) has an orbitally asymptotically stable positive order-1 periodic solution. In (iii) of case (C2), model (3) has an orbitally asymptotically stable positive order-2 periodic solution. This proves (a).
Next, we prove conclusion (b). If for any \(x_0\in (0,x^*]\), trajectory \(O^+(A_0,t_0)\) of model (3) starting from initial point \(A_0(x_0,(1-q)Y_c)\) cuts the line \(L: y=Y_c\) at point \(A_1(x_1,Y_c)\), where \(x^*<x_1<\frac{\lambda }{\delta }\), then for any \(E(x,Y_c)\in {\varGamma }_2\), trajectory \(O^+(E,t_0)\) will intersect with section \({\varGamma }_2\) infinitely many times due to state feedback control strategies and the fact \((1-p)\frac{\lambda }{\delta }+\tau <x^*\). Similar to conclusion (a), we can also obtain that model (3) has a positive order-1 or order-2 periodic solution, which is orbitally asymptotically stable. Conclusion (b) thus follows.
Finally, we turn to (c). If for any \(x\in (0,x^*]\), the trajectory \(O^+(G,t_0)\) of model (3) starting from point \(G(x,(1-q)Y_c)\) does not intersect with line \(L:y=Y_c\), then trajectory starting from point \((x,(1-q)Y_c)\) of section \({\varGamma }_1\) with \(x\in (0,x^*]\) will tend to virus equilibrium \((x^*,y^*)\) and not intersect with section \({\varGamma }_2\). Furthermore, any other trajectory intersects section \({\varGamma }_2\) at most finitely many times, and then tends to virus equilibrium \((x^*,y^*)\) due to \((1-p)\frac{\lambda }{\delta }+\tau <x^*\). In this case, model (3) has no positive order-\(k\) \((k\ge 1)\) periodic solution and virus equilibrium \((x^*,y^*)\) is globally asymptotically stable. This is (c). This completes the proof.
Remark 3
From conclusion (a) of Theorem 3, we note that \((1-p)\frac{\lambda }{\delta }<\hat{x}\) is a sufficient condition for system (3) to have a positive order-1 or order-2 periodic solution.
Remark 4
Conclusion (c) of Theorem 3 shows that the state feedback control strategies are invalid when the controlled strengths \(p\), \(q\), and \(\tau \) remain at a relatively low level and threshold value \(Y_c\) is greater than \(y^*\).
4 Analysis of the virus dynamics model (4)
Similar to Sect. 3, we also divided two cases to give some sufficient conditions for the existence and stability of periodic solutions for model (4).
4.1 Case of \(X_r\le x^*\)
The following result is on the extinct of \(y\).
Theorem 4
For any solution \((x(t),y(t))\) of model (4) with initial condition \((x_0,y_0)\in R_+^2\), if \(x_0\le X_r\) then \(\lim _{t\rightarrow \infty }y(t)=0\).
Proof
Suppose that \((x(t),y(t))\) be any solution of model (4) with initial value \((x_0,y_0)\in \mathbb {R}_+^2\) and \(x_0\le X_r\). Since any trajectory of model (4) starting from domain \(I\) will enter into domain \(I\!I\) , where \(I\) and \(I\!I\) are given by (5), solution \((x(t),y(t))\) must intersect with section \({\varGamma }_4\) at a point \(A(X_r,y_1)\), and then jump to section \(\gamma _3\) at point \(A^+((1-q)X_r+\tau ,(1-q)y_1)\) due to state feedback control strategies. Obviously,
So we consider only trajectories which start from points \((x_0,y_0)\in {\varGamma }_3\) and \(y_0\le \omega \).
Suppose that trajectory \(O^+(E_0,t_0)\) of model (4) starting from point \(E_0((1-p)X_r+\tau ,y_0)\) (where \(0<y_0\le \omega \)) first intersect section \({\varGamma }_4\) at point \(E_1(X_r,y_1)\) and then jumps to point \(E_1^+((1-p)X_r+\tau ,y_1^+)\) on section \({\varGamma }_3\) due to control strategies, and reaches section \({\varGamma }_4\) at point \(E_2(X_r,y_2)\) again, where \(y_1\), \(y_2\in (0,\frac{\lambda -\delta X_r}{\beta X_r})\). Repeating the above process, we can get two point sequences \(\{E_n^+((1-p)X_r+\tau ,y_n^+)\}\) and \(\{E_n(X_r,y_n)\}\), where \(y_n^+=(1-q)y_n\). The corresponding impulsive time sequences are marked by \(\{t_n\}\).
Integrating both sides of the first equation of model (4) from the orbit \(\widetilde{E_n^+E}_{n+1}\) we have
Further, integrating both sides of the second equation of model (4) from the orbit \(\widetilde{E_n^+E}_{n+1}\), we can get
This together with (19) shows that
From the above recursive formula, it can be easily shown that
This results in \(\lim _{n\rightarrow \infty }y_n=0\), where we used the facts \((1-p)X_r+\tau <X_r\) and \(\beta X_r-\mu <0\). Moreover, \(\lim _{t\rightarrow \infty }y(t)=0\). This completes the proof.
Let \(y\equiv 0\) for \(t\in [0,\infty )\) in model (4), we can get the following reduce model
It can be easy to calculate that model (20) has a \(T\) periodic solution which is given by
for all \(nT<t\le (n+1)T\), where
\(x(0)=(1-p)X_r+\tau \) and \(x(T)=X_r\). This means that model (4) has the following semi-trivial periodic solution for \(nT<t\le (n+1)T\) \((n=0,1,2,\ldots )\)
On the stability of this semi-trivial periodic solution, we have the following Theorem 5.
Theorem 5
For any \(p\), \(q\in (0,1)\), \(\tau >0\), and \(X_r\le x^*\), the semi-trivial periodic solution (21) is orbitally asymptotically stable.
Proof
Suppose that \((\widetilde{x}(t),\widetilde{y}(t))\) is a solution of of small-amplitude perturbation of periodic solution \((\phi (t),\psi (t))\) with initial value \((\widetilde{x}(0),\widetilde{y}(0))=((1-p)X_r+\tau ,\widetilde{y}_0)\), which first intersects section \({\varGamma }_4\) at point \((X_r,\widetilde{y}_1)\) and then jumps to point \(((1-p)X_r+\tau ,\widetilde{y}_1^+)\). Further, solution \((\widetilde{x}(t),\widetilde{y}(t))\) intersects section \({\varGamma }_4\) at point \((X_r,\widetilde{y}_2)\) again. Repeating the above process, we have two point sequences \(\{((1-p)X_r+\tau ,\widetilde{y}_n^+)\}\) and \(\{(X_r,\widetilde{y}_n)\}\), where \(\widetilde{y}_n^+=(1-q)\widetilde{y}_n\). Further, by Theorem 4, it is clear that \(\lim _{t\rightarrow \infty }\widetilde{y}(t)=0\). This shows that semi-trivial periodic solution (21) is orbitally asymptotically stable. This completes the proof.
4.2 Case of \(x^*<X_r<\frac{\lambda }{\delta }\)
From the properties (5) of the phase space of model (4), we know that there a point \(E_0((1-p)X_r+\tau ,\gamma )\) such that trajectory \(O^+(E_0,t_0)\) is tangent to the section \({\varGamma }_4\) at the intersection \(E_1(X_r,\frac{\lambda -\delta X_r}{\beta X_r})\) of section \({\varGamma }\) and isocline \(\frac{\mathrm {d}x}{\mathrm {d}t}=0\), where \(\gamma =\gamma (X_r,p,\tau )\).
If \((1-q)\frac{\lambda -\delta X_r}{\beta X_r}<\gamma \), then trajectory \(O^+(E,t_0)\) starting from point \(E(X_r,y)\) (where \(y\le \frac{\lambda -\delta X_r}{\beta X_r}\)) will intersect with \({\varGamma }_4\) infinitely many times due to control strategies \(\triangle x=-pX_r+\tau \) and \(\triangle y(t)=-qy(t)\).
For any two points \(E_i(X_r,y_i)\) and \(E_j(X_r,y_j)\) on section \({\varGamma }_4\), where \(y_i\), \(y_j\in (0,\frac{\lambda -\delta X_r}{\beta X_r})\), and \(y_i<y_j\), in view of control strategies, \(E_i^+((1-p)X_r+\tau ,(1-q)y_i)\) is on below \(E_j^+((1-p)X_r+\tau ,(1-q)y_j)\). Namely, \((1-q)y_i<(1-q)y_j\). It implies that
Similar to the discussion of Theorem 3, we give the following result on the existence and orbital stability of periodic solution.
Theorem 6
For any \(p\), \(q\in (0,1)\), \(\tau >0\), and \(X_r>x^*\), if
then model (4) admits an order-1 periodic solution which is orbitally asymptotically stable.
Remark 5
We want to make clear that the order-1 periodic solution in Theorem 6 may be a positive periodic solution, or it may also be a semi-trivial periodic solution by condition (22).
5 Numerical simulation and discussion
To illustrate the theoretical results and the feasibility of state control strategies, we perform numerical simulations for different control parameters \(p\), \(q\), \(\tau \), \(Y_c\), and \(X_r\) using the software MATLAB. We fixed parameters as those in papers [4, 8, 11] with \(\lambda =50\), \(\delta =0.1\), \(\beta =0.002\), and \(\mu =0.5\) in model (1). It is easy to calculate that model (1) has a unique globally asymptotically stable virus equilibrium \((x^*,y^*)=(250,50)\), which is illustrated here by the blue line in Fig. 1a.
Firstly, we choose the control parameters to be \(p=0.3\), \(q=0.5\), \(\tau =50\), and \(Y_c=45<y^*\), respectively. Using Theorem 1, we know that model (3) has a positive order-1 periodic solution \((\phi (t),\psi (t))\) which is shown by red line in Fig. 1a. At the same time, Fig. 1a also shows that the periodic solution starting from the initial point \((\phi (0),\psi (0))=(251.8965,22.5000)\). Further, by the condition of Corollary 1, it can be easily shown that
Therefore, it follows that \((\phi (t),\psi (t))\) is orbitally asymptotically stable by Corollary 1, which is shown in Fig. 1b.
Further, we fixed the values of \(q\), \(\tau \), and \(Y_c\) as mentioned before and choose \(p\) to be \(0.2\), \(0.4\), \(0.6\), and \(0.8\), respectively. Numerical simulations show that the density of \(y\) can be controlled within a certain limit and the period \(T\) of order-1 periodic solution for model (3) increases with the increase of immune strength \(p\), which is shown in Fig. 2a. The similar numerical results are presented in Fig. 2b.
Secondly, let \(p=0.3\), \(q=0.5\), \(\tau =50\), and \(Y_c=55>y^*\); it is easy to see that model (3) has a positive locally orbitally asymptotically stable order-1 periodic solution by parts (a) or (b) of Theorem 3, which is shown in Fig. 3. However, we only change model parameter \(Y_c=56\), and fixed \(x(0)=250\) and \(y(0)\in [5,35]\), it can be found that the attraction domains of positive order-1 periodic solution and virus equilibrium appear alternately from Fig. 4. It is different from the case of \(Y_c=55\). Therefore, the value of control threshold value \(Y_c\) plays an important role in the treatment of various diseases.
Further, we fixed \(Y_c=55\), \(q=0.5\), and \(\tau =50\) and change the control strength \(p\) to be \(0.35\), \(0.45\), \(0.55\), \(0.65\), and \(0.75\), respectively. Numerical simulations in Fig. 5 show that trajectories intersect with section \({\varGamma }_2\) finitely many times and then tend to virus equilibrium \((x^*,y^*)\) when \(p\) is small, but trajectories intersect with section \({\varGamma }_2\) infinitely many times and model has a orbitally asymptotically stable periodic solution with increase of parameter \(p\). The plots in Fig. 5 show also that control strength \(p\) has only a small effect on the periodic of treatment cycle. Similar results can be obtained if we choose \(q\) as a control parameter. However, if we fixed \(Y_c=55\), \(p=0.3\), \(q=0.5\), and change the control strength \(\tau \) to be \(120\), \(80\), \(40\), and \(0\), respectively. As shown in Fig. 6, the periodic of treatment cycle decrease as control strength \(\tau \) increases. This implies that when we reduce the concentration of infected cells (or free virus) we also should reasonably control the concentration of uninfected cells in the course of treating disease because the high concentrations of uninfected cells will promote the growth of infected cells (or free virus).
Thirdly, we discuss the existence and stability of semi-trivial periodic solution of model (4). Here, we choose \(p=0.3\), \(q=0.5\), \(\tau =5\), and \(X_r=250=x^*\). By Theorem 5, it follows that model (4) has a orbitally asymptotically stable semi-trivial periodic solution, which are shown in Fig. 7. Theoretical results and numerical simulations show clearly that infected cells (or free virus) \(y\) dies out if we control the density of uninfected cells \(x\) within a certain range. It also provides theoretical basis for finding a new measure to prevent the spread of virus disease. In fact, this is also consistent with the results in [42], where Ross proposed a mathematical model to study the spread between human beings and mosquitoes for malaria in earlier 1911. A concept of threshold density is introduced and it is concluded that “ \(\cdots \) in order to counteract malaria anywhere we need not banish Anopheles there entirely–we need only to reduce their numbers below a certain figure.”
Finally, fixing \(p=0.5\), \(\tau =10\), and \(X_r=350>x^*\), model (4) has an orbitally asymptotical stability of positive order-1 periodic solution when \(q=0.4385\) as shown in Fig. 8. At the same time, many numerical simulations show that the attraction domain of this periodic solution is very small; that is, the stability of this periodic solution is very sensitive to the initial value of solution. The plots in Fig. 9 show that the solution tends to virus equilibrium for \(q=0.43\), tends to positive order-1 periodic solution for \(q=0.4385\), and tends to the semi-trivial periodic solution for \(q=0.44\), respectively. It implies that control strength \(q\) is sensitive for the persistence and extinction of infected cells (or free virus) \(y\) for \(X_r>x^*\). This is the relevant contents of Theorem 6 and Remark 5.
6 Concluding remarks
The dynamic complexity of two novel virus epidemic models with state feedback control strategies are analyzed systemically in this paper. Firstly, choosing the density of infected cells or free virus as a control threshold value, some sufficient conditions on the existence and orbitally asymptotical stability of positive order-1 or order-2 periodic solution of model (3) are presented by the Poincaré map, the analog of Poincaré criterion and qualitative analysis method. This amounts to that we can control the density of infected cells or free virus at a low level over a long period of time by adjusting control strength. It is concluded that the state feedback control strategies are feasible and effective. Further, choosing the density of uninfected cells \(x\) as control threshold value, we obtain that the existence and orbitally asymptotical stability of periodic solution. In particular, the existence and orbitally asymptotical stability of semi-trivial periodic are obtained. This implies that infected cells or free virus \(y\) dies out if we control the density of cells \(x\) within a certain range. It provides theoretical basis for finding a new measure to prevent and control the spread of viral disease.
References
Kermack, W.O., McKendrick, A.G.: Contributions to the mathematical theory of epidemics. Proc. R. Soc. Lond. Ser. A 141, 94–122 (1933)
Anderson, R.M.: Epidemiology. In: Cox, F.E.G. (ed.) Modern parasitology. Blackwell Scientific Publications Ltd, Oxford (1993)
Anderson, R.M., May, R.M.: Population biology of infectious dieases: I. Nature 280, 361–367 (1979)
Bonhoeffer, S., Coffin, J.M., Nowak, M.A.: Human immunodeficiency virus drug therapy and virus load. J. Virol. 71, 3275–3278 (1997)
Bonhoeffer, S., May, R.M., Shaw, G.M., Nowak, M.A.: Virus dynamics and drug therapy. Proc. Natl. Acad. Sci. USA 94, 6971–6976 (1997)
Nowak, M.A., Bonhoeffer, S., Hill, A.M., Boehme, R., Thomas, H.C.: Viral dynamics in hepatitis B virus infection. Proc. Natl. Acad. Sci. USA 93, 4398–4402 (1996)
Nowak, M.A., Bangham, C.R.M.: Population dynamics of immune responses to persistent viruses. Science 272, 74–79 (1996)
Ji, Y., Min, L.Q., Zheng, Y., Su, Y.M.: A viral infection model with periodic immune response and nonlinear CTL response. Math. Comput. Simul. 80, 2309–2316 (2010)
Leenheer, P.D.E., Smit, H.L.: Virus dynamics: a global analysis. SIAM J. Appl. Math. 63, 1313–1327 (2003)
Spouge, J.L., Shrager, R.I., Dimitrov, D.S.: HIV-1 infection kinetics in tissue cultures. Math. Biosci. 138, 1–22 (1996)
Wang, K.F., Wang, W.D., Pang, H.Y., Liu, X.N.: Complex dynamic behavior in a viral model with delayed immune response. Phys. D 226, 197–208 (2007)
Coffin, J.M.: HIV population dynamics in vivo: implications for genetic variation, pathogenesis, and therapy. Science 267, 483–489 (1995)
Ho, D.D., Neumann, A.U., Perelson, A.S., Chen, W., Leonard, J.M., Markowitz, M.: Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection. Nature 373, 123–126 (1995)
Bartholdy, C., Christensen, J.P., Wodarz, D., Thomsen, A.R.: Persistent virus infection despite chronic cytotoxic T-lymphocyte activation in Gamma interferon-deficient mice infected with lymphocytic chroriomeningitis virus. J. Virol. 74, 10304–10311 (2000)
Wodarz, D., Christensen, J.P., Thomsen, A.R.: The importance of lytic and nonlytic immune responses in viral infections. Trends Immunol. 23, 194–200 (2002)
Chirove, F., Lungu, E.: Effects of replicative fitness on competing HIV strains. BioSystems 113, 28–36 (2013)
Feng, Z.L., Velasco-Hernandez, J., Tapia-Santos, B., A. Leite, M.C.: A model for coupling within-host and between-host dynamics in an infectious disease. Nonlinear Dyn. 68, 401–411 (2012)
Pang, J.H., Cui, J.A., Hui, J.: The importance of immune responses in a model of hepatitis B virus. Nonlinear Dyn. 67, 723–734 (2012)
Sharma, S., Samanta, G.P.: Dynamical behaviour of an HIV/AIDS epidemic model. Differ. Equ. Dyn. Syst. (2013). doi:10.1007/s12591-013-0173-7
Wang, S.L., Wang, S.L., Song, X.Y.: Hopf bifurcation analysis in a delayed oncolytic virus dynamics with continuous control. Nonlinear Dyn. 67, 629–640 (2012)
Wang, X., Tao, Y., Song, X.: A delayed HIV-1 infection model with Beddington-DeAngelis functional response. Nonlinear Dyn. 62, 67–72 (2010)
Yang, Y.P., Xiao, Y.N., Wang, N., Wu, J.H.: Optimal control of drug therapy: Melding pharmacokinetics with viral dynamics. BioSystems 107, 174–185 (2012)
Yuan, Z.H., Ma, Z.J., Tang, X.H.: Global stability of a delayed HIV infection model with nonlinear incidence rate. Nonlinear Dyn. 68, 207–214 (2012)
Yu, Y., Nieto, J.J., Torres, A., Wang, K.F.: A viral infection model with a nonlinear infection rate. Bound. Value Probl. 2009, Article ID 958016, 19 pages (2009)
Wang, T.L., Hu, Z.X., Liao, F.C.: Stability and Hopf bifurcation for a virus infection model with delayed humoral immunity response. J. Math. Anal. Appl. 411, 63–74 (2014)
Nakata, Y.: Global dynamics of a viral infection model with a latent periodic and Beddington-DeAngelis response. Nonlinear Anal. 74, 2929–2940 (2011)
Wang, Z., Zhao, X.Q.: A within-host virus model with periodic multidrug therapy. Bull. Math. Biol. 75, 543–563 (2013)
Ball, C.L., Gilchrist, M.A., Coombs, D.: Modeling within-host evolution of HIV: mutation, competition and strain replacement. Bull. Math. Biol. 69, 2361–2385 (2007)
Lang, J., Li, M.Y.: Stable and transient periodic oscillations in a mathematical model for CTL response to HTLV-I infection. J. Math. Biol. 65, 181–199 (2012)
Ghorai, A., Kar, T.K.: Biological control of a predator-prey system through provision of a super predator. Nonlinear Dyn. 74, 1029–1040 (2013)
Jana, S., Kar, T.K.: A mathematical study of a prey-predator model in relevance to pest control. Nonlinear Dyn. 74, 667–883 (2013)
Lv, Y.F., Yuan, R., Pei, Y.Z.: Dynamics in two nonsmooth predator-prey models with threshold harvesting. Nonlinear Dyn. 74, 107–132 (2013)
Zhang, Y.: Stability of discrete-time Markovian jump delay systems with delayed impulses and partly unknown transition transition probabilities. Nonlinear Dyn. 75, 101–111 (2014)
Nie, L.F., Teng, Z.D., Guo, B.Z.: A state dependent pulse control strategy for a SIRS epidemic system. Bull. Math. Biol. 75, 1697–1715 (2013)
Nie, L.F., Teng, Z.D., Torres, A.: Dynamic analysis of an SIR epidemic model with state dependent pulse vaccination. Nonlinear Anal. Real World Appl. 13, 1621–1629 (2010)
Nie, L.F., Teng, Z.D., Nieto, J.J., Jung, I.H.: Dynamic analysis of a two-Language competitive model with control strategies. Math. Probl. Eng. 2013, Article ID 654619, 13 pages (2013)
Jiang, G.R., Lu, Q.S.: The dynamics of a prey-predator model with impulsive state feedback control. Discrete Contin. Dyn. Syst. Ser. B 6, 1310–1320 (2006)
Tang, S.Y., Cheke, R.A.: State-dependent impulsive models of integrated pest management (IPM) strategies and their dynamic consequences. J. Math. Biol. 50, 257–292 (2005)
Lakshmikantham, V., Bainov, D.D., Simeonov, P.S.: Theory of Impulsive Differential Equations. World Scientific, Singapore (1989)
Hale, J., Kocak, H.: Dynamics and Bifurcations. Springer-Verlag, New York (1991)
Simeonov, P.S., Bainov, D.D.: Orbital stability of periodic solutions of autonomous systems with impulse effect. Int. J. Syst. Sci. 19, 2561–2585 (1988)
Ross, R.: The Prevention of Malaria, 2nd edn. Murray, London (1911)
Acknowledgments
The authors would like to thank antonymous referees for their constructive suggestions and comments that substantially improved the original manuscript. This research has been partially supported by the National Natural Science Foundation of China (Grant Nos. 11001235, 11271312, and 11261056), the China Postdoctoral Science Foundation (Grant Nos. 20110491750 and 2012T50836), and the Natural Science Foundation of Xinjiang (Grant No. 2011211B08).
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Nie, LF., Teng, ZD. & Jung, I.H. Complex dynamic behavior in a viral model with state feedback control strategies. Nonlinear Dyn 77, 1223–1236 (2014). https://doi.org/10.1007/s11071-014-1372-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11071-014-1372-7