Codimension-Two Bifurcations of a Simplified Discrete-Time SIR Model with Nonlinear Incidence and Recovery Rates

In this paper, a simplified discrete-time SIR model with nonlinear incidence and recovery rates is discussed. Here, using the integral step size and the intervention level as control parameters, we mainly discuss three types of codimension-two bifurcations (fold-flip bifurcation, 1:3 resonance, and 1:4 resonance) of the simplified discrete-time SIR model in detail by bifurcation theory and numerical continuation techniques. Parameter conditions for the occurrence of codimension-two bifurcations are obtained by constructing the corresponding approximate normal form with translation and transformation of several parameters and variables. To further confirm the accuracy of our theoretical analysis, numerical simulations such as phase portraits, bifurcation diagrams, and maximum Lyapunov exponents diagrams are provided. In particular, the coexistence of bistability states is observed by giving local attraction basins diagrams of different fixed points under different integral step sizes. It is possible to more clearly illustrate the model's complex dynamic behavior by combining theoretical analysis and numerical simulation.