Dynamics and probability density function of a stochastic COVID-19 epidemic model with nonlinear incidence

In this paper, a stochastic SEIQ model with nonlinear incidence rate is proposed. The existence and uniqueness of global positive solution of the stochastic model are proved. Then, by constructing some Lyapunov functions, we derive a sufficient condition for the ergodic stationary distribution when R0s is greater than one. By solving a four-dimensional Fokker-Planck equation, we get the exact expression of log-normal probability density function of stationary distribution for the stochastic model. Sufficient condition for the extinction of the exposed and infected population is also provided. Finally, numerical simulations are presented to verify the above theoretical results.