Stationary distribution and density function analysis of SVIS epidemic model with saturated incidence and vaccination under stochastic environments

In this article, we study the dynamical properties of susceptible-vaccinated-infected-susceptible (SVIS) epidemic system with saturated incidence rate and vaccination strategies. By constructing the suitable Lyapunov function, we examine the existence and uniqueness of the stochastic system. With the help of Khas’minskii theory, we set up a critical value [Formula: see text] with respect to the basic reproduction number [Formula: see text] of the deterministic system. A unique ergodic stationary distribution is investigated under the condition of [Formula: see text] . In the epidemiological study, the ergodic stationary distribution represents that the disease will persist for long-term behavior. We focus for developing the general three-dimensional Fokker–Planck equation using appropriate solving theories. Around the quasi-endemic equilibrium, the probability density function of the stochastic system is analyzed which is the main theme of our study. Under [Formula: see text] , both the existence of ergodic stationary distribution and density function can elicit all the dynamical behavior of the disease persistence. The condition of disease extinction of the system is derived. For supporting theoretical study, we discuss the numerical results and the sensitivities of the biological parameters. Results and conclusions are highlighted.