Abstract
Due to the heterogeneity of contact structure, it is more reasonable to model on networks for epidemics. Because of the stochastic nature of events and the discrete number of individuals, the spread of epidemics is more appropriately viewed as a Markov chain. Therefore, we establish stochastic SIRS models with vaccination on networks to study the mean and variance of the number of susceptible and infected individuals for large-scale populations. Using van Kampen's system-size expansion, we derive a high-dimensional deterministic system which describes the mean behaviour and a Fokker-Planck equation which characterizes the variance around deterministic trajectories. Utilizing the qualitative analysis technique and Lyapunov function, we demonstrate that the disease-free equilibrium of the deterministic system is globally asymptotically stable if the basic reproduction number R (0) < 1; and the endemic equilibrium is globally asymptotically stable if R (0) > 1. Through the analysis of the Fokker-Planck equation, we obtain the asymptotic expression for the variance of the number of susceptible and infected individuals around the endemic equilibrium, which can be approximated by the elements of principal diagonal of the solution of the corresponding Lyapunov equation. Here, the solution of Lyapunov equation is expressed by vectorization operator of matrices and Kronecker product. Finally, numerical simulations illustrate that vaccination can reduce infections and increase fluctuations of the number of infected individuals and show that individuals with greater degree are more easily infected.