Rabinovich–Fabrikant equations

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

The Rabinovich–Fabrikant equations are a set of three coupled ordinary differential equations exhibiting chaotic behavior for certain values of the parameters. They are named after Mikhail Rabinovich and Anatoly Fabrikant, who described them in 1979.

System description

The equations are:[1]

\dot{x} = y (z - 1 + x^2) + \gamma x \,
\dot{y} = x (3z + 1 - x^2) + \gamma y \,
\dot{z} = -2z (\alpha + xy), \,

where α, γ are constants that control the evolution of the system. For some values of α and γ, the system is chaotic, but for others it tends to a stable periodic orbit.

Danca and Chen[2] note that the Rabinovich–Fabrikant system is difficult to analyse (due to the presence of quadratic and cubic terms) and that different attractors can be obtained for the same parameters by using different step sizes in the integration.

Equilibrium points

File:Rabinovich-Fabrikant equilibrium point existence regions.svg
Graph of the regions for which equilibrium points \tilde{\mathbf{x}}_{1,2,3,4} exist.

The Rabinovich–Fabricant system has five hyperbolic equilibrium points, one at the origin and four dependent on the system parameters α and γ:[2]

\tilde{\mathbf{x}}_0 = (0,0,0)
\tilde{\mathbf{x}}_{1,2} = \left( \pm q_-, - \frac{\alpha}{q_-}, 1- \left(1-\frac{\gamma}{\alpha}\right)q_-^2 \right)
\tilde{\mathbf{x}}_{3,4} = \left( \pm q_+, - \frac{\alpha}{q_+}, 1- \left(1-\frac{\gamma}{\alpha}\right)q_+^2 \right)

where

q_{\pm} = \sqrt{ \frac{ 1 \pm \sqrt{ 1- \gamma \alpha \left( 1- \frac{3 \gamma}{4\alpha} \right) } }{2 \left(1- \frac{3\gamma}{4\alpha}\right) }}

These equilibrium points only exist for certain values of α and γ > 0.

γ = 0.87, α = 1.1

An example of chaotic behavior is obtained for γ = 0.87 and α = 1.1 with initial conditions of (−1, 0, 0.5).[3] The correlation dimension was found to be 2.19 ± 0.01.[4] The Lyapunov exponents, λ are approximately 0.1981, 0, −0.6581 and the Kaplan–Yorke dimension, DKY ≈ 2.3010[3]

γ = 0.1

Danca and Romera[5] showed that for γ = 0.1, the system is chaotic for α = 0.98, but progresses on a stable limit cycle for α = 0.14.

File:Rabinovich-Fabrikant LimitCicle.PNG
3D parametric plot of the solution of the Rabinovich-Fabrikant equations for α=0.14 and γ=0.1 (limit cycle is shown by the red curve)

See also

References

  1. Lua error in package.lua at line 80: module 'strict' not found.
  2. 2.0 2.1 Lua error in package.lua at line 80: module 'strict' not found.
  3. 3.0 3.1 Lua error in package.lua at line 80: module 'strict' not found.
  4. Lua error in package.lua at line 80: module 'strict' not found.
  5. Lua error in package.lua at line 80: module 'strict' not found.

External links