Stability and Hopf Bifurcation for a Delayed SLBRS Computer Virus Model

By incorporating the time delay due to the period that computers use antivirus software to clean the virus into the SLBRS model a delayed SLBRS computer virus model is proposed in this paper. The dynamical behaviors which include local stability and Hopf bifurcation are investigated by regarding the delay as bifurcating parameter. Specially, direction and stability of the Hopf bifurcation are derived by applying the normal form method and center manifold theory. Finally, an illustrative example is also presented to testify our analytical results.


Introduction
For the purpose of studying the prevalence of computer viruses and posing effective strategies of controlling their spread across the Internet and considering the similarity between the propagation of computer virus across the Internet and that of biological viruses across a population, many scholars have proposed computer virus models by properly modifying some epidemiological models [1][2][3][4][5][6]. In [2], Han and Tan proposed a computer virus model using an SIRS model and investigated the local stability and Hopf bifurcation for the endemic state of the model. In [4], Mishra and Pandey proposed an SEIRS model for the transmission of worms in computer network through vertical transmission and they obtained the threshold value determining whether the worms die out. All the models above assume that an infected computer in which the virus resides is in latency and cannot infect other computers.
Recently, Yang et al. [7] proposed the following SLBRS computer virus model by taking into account the fact that a computer immediately possesses infection ability once it is infected. Consider where ( ) denotes the percentage of computers that have no immunity at time ; ( ) denotes the percentage of infected computers that are latent at time ; ( ) denotes the percentage of infected computers that are breaking out at time ; and ( ) denotes the percentage of uninfected computers that have temporary immunity at time . is the constant rate for the external computers connecting to the Internet and that for the internal computers disconnecting from the Internet. is the rate of contact for the virus-free computer. is the rate coefficient from class to class , is the rate coefficient from class to class , and is the rate coefficient from class to 2 The Scientific World Journal class . Yang et al. [7] studied global stability of the model by constructing a Lyapunov function for the model. It is well known that there may be a time lag when the computer uses antivirus software to clean the virus and computer virus models with time delay have been studied by many scholars [1,2,6,[8][9][10][11][12][13][14][15]. In [8], Feng et al. investigated the Hopf bifurcation of delayed SIRS computer virus propagation model by regarding the latent period of computer virus as a bifurcation parameter. In [10], Liu investigated the Hopf bifurcation of a delayed SEIQRS model for the transmission of malicious objects in computer network by regarding the time delay due to the temporary immunity period as a bifurcation parameter. To the best of our knowledge, few papers deal with the research of Hopf bifurcation of system (1). Motivated by the work above and considering that there may be a time lag when the computer uses antivirus software to clean the virus, we incorporate time delay into system (1) and investigate the Hopf bifurcation of the following delayed SLBRS system: where is the time delay due to the period that a computer uses antivirus software to clean the virus.
The main purpose of this paper is to investigate the effect of delay on system (2). This paper is organized as follows. In Section 2, sufficient conditions for local stability and existence of local Hopf bifurcation are obtained by analyzing the distribution of the roots of the associated characteristic equation. In Section 3, direction and stability of the Hopf bifurcation are studied by applying the normal form theory and center manifold theorem. In Section 4, some numerical simulations are presented to verify the theoretical results.
If condition (H 2 ) holds, we know that (12) has at least a positive root 0 such that (6) has a pair of purely imaginary roots ± 0 . For 0 , the corresponding critical value of the delay is where 0 = − 0 0 , 2 = 0 2 − 1 1 + 2 0 , Taking the derivative with respect to on both sides of (6), we obtain Thus, According to the Hopf bifurcation theorem in [16], we have the following results.   (24)
Following the algorithms given in [16] and using similar computation process in [17], we can get the following important coefficients: with 20 ( ) = 20 (0) where 1 and 2 can be determined by the following equations, respectively, with (1) The Scientific World Journal 5 Then, we can get the following coefficients: In conclusion, we have the following results.

Numerical Simulation
In this section, we present a numerical example to justify the theoretical results above. We use the same coefficients as those used in [7] and consider the following system: from which we obtain the basic reproduction number 0 = 9.944 > 1 and the unique positive equilibrium * (0.1006, 0.1160, 0.2762, 0.5022) of system (36). Further, we obtain 0 = 2.3091, 0 = 58.5472. According to Theorem 1 in Section 2, we know that, when ∈ [0, 58.5472), the positive equilibrium * (0.1006, 0.1160, 0.2762, 0.5022) of system (36) is locally asymptotically stable. As can be seen from Figures 1 and 2, when = 55.2386 < 0 , the positive equilibrium * (0.1006, 0.1160, 0.2762, 0.5022) of system (36) is locally asymptotically stable. However, when the delay passes through the critical value of 0 = 58.5472, the positive equilibrium * will lose stability and a Hopf bifurcation occurs. This property can be illustrated by Figures 3 and 4. In addition, by some complex computations and from (35), we obtain 2 = −0.8143 < 0, 2 = −0.9001 < 0, and 2 = 0.0466 > 0. Therefore, we can conclude that the Hopf bifurcation of system (36) is subcritical, the bifurcating periodic solutions are stable, and the period of the periodic solutions decreases according to Theorem 2 in Section 3.

Conclusions
A delayed SLBRS computer virus model has been studied in this paper based on the model proposed in [7]. Compared with literature [7], we mainly consider the effect of the time delay due to the period that the computers use antivirus software to clean the virus on the model. The main results are presented in terms of local stability and local Hopf bifurcation. Sufficient conditions for local stability and existence of local Hopf bifurcation are obtained by analyzing the distribution of the roots of the associated characteristic equation. It has been shown that when the delay is suitably small ( < 0 ), the computer virus model is locally asymptotically stable. In such case, the propagation of the computer virus can be easily controlled. However, once the delay passes though the critical value 0 , a Hopf bifurcation occurs and a family of periodic solutions bifurcates from the positive equilibrium of the model. This is not welcome in networks because the computer virus will be out of control in this case. Further, direction and stability of the Hopf bifurcation have also been determined by using the normal form method and center manifold theory. Numerical simulations are also included to verify the analytical predictions.