Some Bifurcations of Codimensions 1 and 2 in a Discrete Predator–Prey Model with Non-Linear Harvesting

This paper delves into the dynamics of a discrete-time predator–prey system. Initially, it presents the existence and stability conditions of the fixed points. Subsequently, by employing the center manifold theorem and bifurcation theory, the conditions for the occurrence of four types of codimension 1 bifurcations (transcritical bifurcation, fold bifurcation, flip bifurcation, and Neimark–Sacker bifurcation) are examined. Then, through several variable substitutions and the introduction of new parameters, the conditions for the existence of codimension 2 bifurcations (fold–flip bifurcation, 1:2 and 1:4 strong resonances) are derived. Finally, some numerical analyses of two-parameter planes are provided. The two-parameter plane plots showcase interesting dynamical behaviors of the discrete system as the integral step size and other parameters vary. These results unveil much richer dynamics of the discrete-time model in comparison to the continuous model.