Back

Predator-prey and the Hopf bifurcation

What you are seeing: the Rosenzweig-MacArthur predator-prey model with Holling Type II response: x˙=rx(1x/K)axy/(b+x)\dot{x} = r x (1 - x / K) - a x y / (b + x), y˙=eaxy/(b+x)dy\dot{y} = e a x y / (b + x) - d y. Prey (xx) grow logistically and are consumed by predators; predators (yy) grow by eating prey and die at rate dd. Increasing the prey carrying capacity KK destabilizes the coexistence equilibrium through a supercritical Hopf bifurcation at KH0.7K_H \approx 0.7. Below KHK_H the populations damp to a fixed point; above it they oscillate on a stable limit cycle.

Watch the phase portrait (x,y)(x, y) on the left and the time series on the right. Sliding KK from 0.4 to 2.0 sweeps through the Hopf: small KK gives damped spirals; large KK gives clean limit cycles whose amplitude grows as KKH\sqrt{K - K_H}.

Figure 1. Rosenzweig-MacArthur model with Hopf bifurcation. Method: RK4 with analytic equilibrium markers.
K1.50
speed3

WHAT TO TRY

  • Vary each control and watch the rail readouts respond.
  • Compare the diagnostic plot against the live scene.