Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2017 / Article

Research Article | Open Access

Volume 2017 |Article ID 9563862 | 10 pages |

Delay Induced Hopf Bifurcation of an Epidemic Model with Graded Infection Rates for Internet Worms

Academic Editor: Marco Spadini
Received01 Sep 2017
Revised25 Nov 2017
Accepted03 Dec 2017
Published24 Dec 2017


A delayed SEIQRS worm propagation model with different infection rates for the exposed computers and the infectious computers is investigated in this paper. The results are given in terms of the local stability and Hopf bifurcation. Sufficient conditions for the local stability and the existence of Hopf bifurcation are obtained by using eigenvalue method and choosing the delay as the bifurcation parameter. In particular, the direction and the stability of the Hopf bifurcation are investigated by means of the normal form theory and center manifold theorem. Finally, a numerical example is also presented to support the obtained theoretical results.

1. Introduction

In the wake of developments in computer technology and communication technology, there is a rapid increase in computer viruses which has brought about huge financial losses [13]. Therefore, it is extremely urgent to analyze and protect computers against viruses. In order to understand the spread law of computer viruses over the Internet and in view of the high similarity between computer viruses and biological viruses, many computer virus propagation models have been developed and analyzed. For example, see [413] and the cited references therein.

All the models above assume that the infected computer has no infectivity. This is inconsistent with the fact that an infected computer which is in latency can also infect other computers through file copying or file downloading. Based on this, Yang et al. proposed some models [1417], by taking into account the fact that a computer immediately possesses infectivity once it is infected. However, these models make an assumption that the exposed computers and the infectious computers have the same infectivity. This is not consistent with the reality, because the infection rate of the exposed computers is less than that of the infectious ones. Thus, Wang et al. [18] proposed the following SEIQRS model with graded infection rates for Internet worms:where , , , , and denote the numbers of the susceptible, the exposed, the infectious, the quarantined, and the recovered computers at time in the Internet, respectively. is the recruitment of the susceptible computers; and are the contact rates of the exposed computers and the infectious computers, respectively; is the natural death rate of the computers; and are the death rates of the exposed computers and the infectious computers due to the attack of worms, respectively; is the quarantined rate of the infectious computers; is the average cured time; , , , , and are the state transition rates. Wang et al. [18] investigated local and global stability of system (1).

It should be pointed out that Wang et al. [18] neglect the fact that it needs a period to clean the worms in the exposed, the infectious, and the quarantined computers for the antivirus software. Time delays cause a stable equilibrium to become unstable and cause Hopf bifurcation phenomenon for a dynamical system. The occurrence of Hopf bifurcation means that the state of worm prevalence changes from an equilibrium point to a limit cycle, which is not welcomed in networks. Therefore, it is meaningful to investigate the effect of time delays on stability of dynamical systems.

As far as we know, there have been some researches on Hopf bifurcation of delayed computer virus models. For example, Feng et al. [19] studied the Hopf bifurcation of a delayed SIRS viral infection model in computer networks by taking the time delay due to the latent and temporary immune period as the bifurcation parameter. Dong et al. [9] proposed a delayed SEIR computer virus model with multistate antivirus and studied the Hopf bifurcation of the model by choosing the delay where the infectious nodes use antivirus software to clean the viruses as the bifurcation parameter. Zhang and Bi [20] investigated the Hopf bifurcation of a delayed computer virus propagation model with infectivity in latent period. For some other research works on the Hopf bifurcation of computer virus models one can refer to [2124]. Specially, Zhang et al. [24] studied the existence and properties of the Hopf bifurcation of a computer virus model with antidote in vulnerable system by regarding the time delay due to the period that the infected computers use to reinstall system as a bifurcation parameter. Motivated by the work above and considering the effect of the time delay on system (1), we consider the following delayed worm propagation model:where is the time delay due to the period that the antivirus software uses to clean the worms in the exposed, the infectious, and the quarantined computers.

The object of this paper is to study the existence and properties of the Hopf bifurcation of system (2). The remainder of this paper is organized as follows. The existence of a Hopf bifurcation is discussed by choosing the delay as the bifurcation parameter in Section 2. Properties of the Hopf bifurcation are studied by means of the normal form theory and the center manifold theorem. An example together with its numerical simulations is also presented in order to illustrate the effectiveness of our obtained theoretical results.

2. Stability of the Viral Equilibrium and Existence of Hopf Bifurcation

By a simple computation, we know that if and , then system (2) has a unique endemic equilibrium , where , , , , andThen, we obtain the characteristic equation of the linearized system at :whereFor , (4) becomeswith

According to Routh-Hurwitz criteria when the condition is satisfied, that is, when (8)–(12) hold, then is locally asymptotically stable when .Multiplying on both sides of (4), one can getFor , we assume that is the root of (13); thenwhereThen,withSince , (16) can be transformed intoIt equalswhere

Based on the discussion about the distribution of the root of (19) in [25, 26], we obtain the expression of , saySubstitute (21) into (16), then the expression of can be obtained, sayHence, we have

In order to obtain the main results in this paper, we suppose that ((23)) has at least one positive root . For , we have Differentiating (13) with respect to , we obtain where Further,

Therefore, if condition , , holds, then . Based on the Hopf bifurcation theorem in [27], we have the following.

Theorem 1. For system (2), if conditions , , and hold, then the endemic equilibrium is locally asymptotically stable when ; a Hopf bifurcation occurs at the endemic equilibrium when and a family of periodic solutions bifurcate from the endemic equilibrium near .

3. Direction and Stability of the Hopf Bifurcation

In this section, we describe the direction, stability, and period of the Hopf bifurcation of system (2) from the endemic equilibrium based on the normal form theory and center manifold theorem given by Hassard et al. [27]. Let , , , , and , and normalize the time delay by . And let , , so that the value is the Hopf bifurcation value. Define the space of continuous real valued functions as . Then, the delayed system (2) can be transformed into the functional differential equation in :where , , , and are given, respectively, bywith

By the Riesz representation theorem, there exists a function whose components are of bounded variation for such thatfor . In view of (29), we can choose with being the Dirac delta function.

For , defineThen, system (28) is then equivalent toFor ), defineand a bilinear inner productwhere . Then (from here onwards we refer by ) and are adjoint operators. Since are the eigenvalues of , they are also the eigenvalues of .

Suppose that and are the eigenvectors for and corresponding to and , respectively. Then, we haveFrom (37), we get such that and .

Proceeding in the same manner as Hassard et al. [27] and the similar computation process as that in [2830], we obtain with and can be obtained by the following two equations:whereThen, one can obtainThus, based on the properties of the Hopf bifurcation discussed in [27], we have the following.

Theorem 2. For system (2), if , then the Hopf bifurcation is supercritical (subcritical); if , then the bifurcated periodic solutions are stable (unstable); if  , then the period of the bifurcated periodic solutions increases (decreases).

4. Numerical Simulation

In this section, we perform some numerical simulations to support and explain our obtained results. We choose a set of parameters as follows: , , , , , , , , , , , , and . Then, system (2) becomes

Then, we get , , , , and the unique endemic equilibrium . Further we have