Back

Duffing oscillator

What you are seeing: a ball in a double-well potential being pushed back and forth at a fixed frequency. Without the push the ball settles into one of the two wells; with a small push it rattles around in one well; with a strong enough push it occasionally jumps over the barrier between the wells. The Duffing equation x¨+δx˙x+x3=γcos(ωt)\ddot x + \delta\dot x - x + x^3 = \gamma\cos(\omega t) is the standard textbook example of a chaotic driven oscillator.

The left panel shows the phase portrait (x,x˙)(x, \dot x). Dots: xx vs x˙\dot x sampled once per drive period T=2π/ωT = 2\pi/\omega ("strobe", or Poincare section). A clean periodic orbit lays down a single dot; period-2 lays two; chaos sprays a fractal cloud. The right panel is the bifurcation diagram: for each value of γ\gamma on the horizontal axis, we plot the strobed xx values vertically. You can see the period-doubling cascade as γ\gamma crosses about 0.40.4.

Figure 1. Left: phase portrait of the Duffing oscillator with stroboscopic Poincare points (one per drive period). Right: bifurcation diagram across γ\gamma, showing the period-doubling cascade to chaos. Method: classical fourth-order Runge-Kutta from shared/js/engine/ode-rk.js, fixed step dt=T/200dt = T/200.
delta0.30
gamma0.500
omega1.20
speed0.6

WHAT TO TRY

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