Abstract

Because the latent period and the infectious period of tuberculosis (TB) are very long, it is not reasonable to consider the time as constant. So this paper formulates a mathematical model that divides the latent period and the infectious period into n-stages. For a general n-stage stage progression (SP) model with bilinear incidence, we analyze its dynamic behavior. First, we give the basic reproduction number . Moreover, if , the disease-free equilibrium is globally asymptotically stable and the disease always dies out. If , the unique endemic equilibrium is globally asymptotically stable and the disease persists at the endemic equilibrium.

1. Introduction

Tuberculosis (TB) is one of the oldest recorded diseases of mankind. It is an disease caused by the infection of bacterium Mycobacterium tuberculosis. The disease is most commonly transmitted from a person suffering from infectious (active) tuberculosis to other persons by infected droplets created when the person with active TB coughs or sneezes. Most infected people do not develop the progressive disease. When the first time after being infected with Mycobacterium tuberculosis, the individual generally will experience a latent phase. TB progresses through a long latent period and a long infectious period. For this case, the infection can vary greatly in time. The progression of a typical TB infection can take four weeks to several years before the development into active disease, and a few individuals directly become infectious without experiencing latency. Moreover, most infected people do not develop the active disease in his or her life. In the infectious period, individual differences lead to different course. The longest infectious period is several decades while the shortest maybe only a few months. Moreover, the treatment of TB has side effects (sometimes quite serious) and takes varying time depending on the other various factors such as nutritional status and/or access to decent medical care and living conditions [1]. The progression of a TB infection goes through several distinct stages. Similarly, HIV virus has the long incubation and infectious periods (infection age, from 8 to 10 years). During the incubation period, the infectivity of infected people is varying depending on the time since infection. When symptom onset appears, AIDS population transmission rate depends on disease age (i.e., the time elapsed since the onset) [2]. Different from common infectious diseases, the time scale of TB or HIV/AIDS transmission is so long that demographic of the host population could affect transmission process. The classic compartmental models that postulate all the hosts are similar and imprecise to describe the spread of an infection. For explore the issue, many authors formulate staged progression (SP) models [1, 39] and delayed epidemic models [10]. In [5], the authors analyze a mathematical model for infectious diseases that progress through distinct stages within infected hosts and prove the global dynamics of the equilibria. Hyman et al. [4] created two SP models to show the impact of variations in infectiousness. In [10], the author formulate a delayed SIR epidemic model by introducing a latent period into susceptible and infectious individuals in incidence rate. Mathematical modeling has become an important tool in analyzing the epidemiological characteristics of infectious diseases and can provide useful control measures. In [1118], several variants and generalizations of the Beverton CHolt model (standard time-invariant, time-varying parameterized, generalized model or modified generalized model) have been investigated at the levels of stability, cycle-oscillatory behavior, permanence, and control through the manipulation of the carrying capacity. De la Sen et al. studied the impact of vaccination for infectious diseases. This paper considers the latent period and the infectious period to formulate a n-stage SP model with bilinear incidence (based on the model in [3]).

To formulate an SP model, the host population is divided into the following epidemiological classes or subgroups: the susceptible compartment S; the latent compartment (infected but not infectious), whose members are in the th stage of the disease progression, where ; the infectious compartment , whose members are in the th stage of the disease progression, where ; the treated compartment . denotes the total population. Here we assume that the latent period is averagely divided into stages and the infectious period is averagely divided into stages. We also assume that hosts in the treated compartment are noninfectious due to inactivity. Using Figure 1, we formulate the following model: The incidence form is , where the force of infection is the bilinear incidence; is the incubation period, and is the average incubation period from the th stage to th stage, for ; is the infectious period, and is the average infectious period from the th stage to th stage, for is the recruitment rate of susceptible; is per-capita treatment rates, respectively; is the per-capita natural death rate (and hence is an average lifespan of the healthy individuals in the population); is the per capita disease-induced death rate. , , , , , , are positive constants. Because in the model (1.1), the first equations do not contain variable , so dynamical behaviors of the model (1.1) are equivalent to the model (1.3).

Let , then using (1.3) we have This implies that . Therefore, the model can be studied in the feasible region

2. The Basic Reproduction Number

The disease-free equilibrium is obtained by setting the right side of each of the differential equations equal to zero in system (1.3).

If , , it is easy to deduce the disease-free equilibrium as follows:

Next, we derive the basic reproductive number of (1.3) by the method of next-generation matrix formulated in [19].

Let . Then the last equations of model (1.3) can be written as where By calculating the Jacobian matrices of and at the disease-free equilibrium , we have

Let be the identity matrix. Solving the matrix equation , we obtain that is a lower triangular matrix where where equals .

is the next-generation matrix for model (1.3). It then follows that the spectral radius of matrix is According to Theorem 2 in [19], the basic reproduction number of model (1.3) is Therefore, gives the number of secondary infectious cases produced by an infectious individual during his or her effective infectious period when introduced in a population of susceptible. If , then becomes unstable and the disease becomes endemic. Moreover, a unique endemic equilibrium exists in the interior of . Next, we prove the uniqueness of the endemic equilibrium when .

3. Equilibria

An equilibrium of (1.3) satisfies where is given in Section 1. The disease-free equilibrium exists for all positive parameter values. Next we consider the existence of endemic equilibria , , , , , . First, we introduce the definition and properties of M-matrices. Most of the texts on matrix theory can find them [20].

Definition 3.1. is a M-matrix if(1)Off-diagonal entries of are nonpositive, and(2) is positively stable, namely, all eigenvalues of have positive real parts.

Proposition 3.2. Properties of M-matrices(1), , , the spectral radius of .(2) is nonsingular and .(3)There exists such that for .According to the above, we know that

Proposition 3.3. The following holds for the matrix defined above.(1) is a M-matrix.(2) exists and is a nonnegative matrix.(3)There exists such that for .By Proposition 3.3, we know that Then, we obtain the result.

Theorem 3.4. If , then is the only equilibrium in . If , then a unique endemic equilibrium exists in the interior of .

Proof . The last equations of (3.1) can be written in the form or
Multiplying the row vector to (3.5), we have
Using , we obtain
Because , then It is clear that the equation has a unique positive solution in the interval when . Substitute into (3.1), and we obtain are unique positive solution. So the endemic equilibrium is unique when .

4. Stability of the Equilibria

In this section, we employ the direct Lyapunov method with a Lyapunov function of the form where is a properly selected constant, is the population of the th compartment, and is the equilibrium level, to study properties of this model. This function is referenced in many papers [2123], including the models with multiple parallel infectious stages [4, 24] and models with nonlinear incidence rates of different forms [2528].

Now we are ready to proceed to the global properties of the model.

Theorem 4.1. If , the disease-free equilibrium is globally asymptotically stable in .

Proof . Consider the function where and (; ;) are properly selected constants. Their values are as follows:
Because , the each part of should be less than 1, that is, < , < . Therefore, we can obtain that are greater than 0. This function is defined and continuous for all .
Its derivative along the solutions to the system (1.3) is
If , then . Note that, only if . The maximum invariant set in is the singleton . The global stability of when follows from LaSalle's invariance principle [29]. The global attractivity of and the Lyapunov function imply that is also locally stable, since otherwise will have a homoclinic orbit that is entirely contained in , contradicting that the largest compact invariant set in is . This establishes the global stability of when .

Theorem 4.2. If , the endemic equilibrium is globally asymptotically stable in .

Proof . Let us consider the function where This function is defined and continuous for all . By compute the derivative of along the solutions to the system (1.3), it follows that
Recalling that , and that , we obtain Applying the inequality we have Therefore, for all , provided that are positive, where the equality holds only on the straight line . It is easy to see that for both these systems, is the only equilibrium state on this line. Therefore, by Lyapunov-LaSalle asymptotic stability theorem [30, 31], the positive equilibrium state is globally asymptotically stable in .

5. Conclusion

According to the different length of the latent period and the infectious period of TB, in this paper, we proposed a general -stage SP model with bilinear incidence to study the transmission dynamics of TB. What to do to make the results more accurate and tally with the actual situation. We prove that the global dynamics are completely determined by the basic reproduction number . If , then the disease-free equilibrium is globally asymptotically stable and the disease always dies out. If , the unique endemic equilibrium is globally asymptotically stable in the interior of the feasible region, and the disease persists at the endemic equilibrium.

Acknowledgments

The work is supported by the National Sciences Foundation of China (10471040) and the National Sciences Foundation of Shanxi Province (2009011005-1).