2012 The "cord" attractor


Observability is an important concept (see here for understanding and interpreting results of several problems in nonlinear dynamics, such as analysis, synchronization, and control (see [1] [2] for instance). The reason for this is that the method’s performance and the reliability of the results depend on the “measured” variable. The "cord" system is of special interest for at least two reasons : i) It produces an attractor which is not equivalent to any known attractor, and ii) the system has two variables that provide poor observability of the dynamics.

The "cord" system is [3]

      \dot{x} = -y -z + a (F-x) \\[0.2cm]       
      \dot{y} = (x-1) \, y -bx \, z + G \\[0.2cm]       
      \dot{z} = bx \, y + (x-1) \, z   \, , 

and corresponds in fact to a variant of the Lorenz 84 system (Here) discovered "by error" by Luis Aguirre. The "cord" system is obtained by replacing the nonlinear terms of the first equation of the Lorenz 84 system by the linear counterparts. The symbolic observability coefficients of the "cord" system are \eta_x^2 = 0.81, \eta_y^2 = \eta_z^2 = 0.34. The good observability provided by the x-variable allowed us to perform a topological characterization of the cord attractor as detailed in [3].

The "cord" system is investigated for the parameter values as (a, F, G)=(0.258, 8, 1), with initial conditions (x_0, y_0, z_0)=(1,1,1). Parameter b is chosen as the bifurcation parameter. A bifurcation diagram (Fig. 1) reveals that one route to chaos for this system is the common period-doubling cascade. After the accumulation point, a chaotic regime is observed. Then, a long period-7 window is followed by a chaotic regime but there is no period-doubling cascade before returning to chaos. It seems that a boundary crisis thus takes place at the end of the period-7 window.

JPG - 15.9 ko
Fig. 1. Bifurcation diagram of the "cord" system.

This attractor is well bounded during the first part of the bifurcation diagram as always observed after a period-doubling cascade. The attractor boundaries become not well defined in the right part of the bifurcation diagram - the tails of the probability density function of visits to the boundary region become evanescent - after a boundary crisis (b \approx 4.046), leading to an attractor with over twice the size of the original one (Fig. 1).

Before the boundary crisis (beyond which it becomes nearly impossible to characterize the attractor due to its complexity), a typical chaotic attractor is shown in Fig. 2 for b=4.033. The main rotation is clearly seen in the y-z plane projection.

JPG - 20.3 ko
Fig. 2. The "cord" attractor for b=4.033.

A first-return map to a Poincaré section is shown in Fig. 3. It is a smooth unimodal map as commonly induced by a period-doubling cascade. The decreasing branch is layered, a signature announcing the crisis observed for slightly greater b-value.

JPG - 7.6 ko
Fig. 3. First-return map to a Poincaré section of the cord attractor.

[1] C. Letellier, L. A. Aguirre & J. Maquet, How the choice of the observable may influence the analysis of non linear dynamical systems, Communications in Nonlinear Science and Numerical Simulation, 11 (5), 555-576, 2006.

[2] C. Letellier & L. A. Aguirre, Interplay between synchronization, observability, and dynamics, Physical Review E, 82, 016204, 2010.

[3] C. Letellier & L. A. Aguirre, Required criteria for recognizing new types of chaos : Application to the “cord” attractor, Physical Review E, 85, 036204, 2012.

ATOMOSYD © 2007-2017 |  Suivre la vie du site  |  SPIP  |  scoty  |  MàJ . 06/02/2017