Back

Lyapunov spectrum of the Henon map

The Henon map xn+1 = 1 − a xn2 + yn, yn+1 = b xn is the canonical 2D chaos benchmark. At (a, b) = (1.4, 0.3) its attractor is a strange set with a positive largest Lyapunov exponent and a negative second exponent that sum exactly to ln|b|. Drag the handle on the parameter panel to explore how the attractor morphs and how the spectrum responds. The trace identity λ1 + λ2 = ln|b| is an exact dynamical invariant that the Benettin QR algorithm preserves to machine precision.

Figure 1. Henon attractor (left) and the (a, b) parameter panel (right). Method: 2D Henon iteration plus a Benettin QR algorithm on a 2x2 tangent frame with modified Gram-Schmidt re-orthonormalization at each step.

WHAT TO TRY

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