Next Article in Journal
An Entropy-Based Approach to Path Analysis of Structural Generalized Linear Models: A Basic Idea
Next Article in Special Issue
Generalized Combination Complex Synchronization for Fractional-Order Chaotic Complex Systems
Previous Article in Journal
Averaged Extended Tree Augmented Naive Classifier
Previous Article in Special Issue
A Colour Image Encryption Scheme Using Permutation-Substitution Based on Chaos
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytic Exact Upper Bound for the Lyapunov Dimension of the Shimizu–Morioka System

1
Mathematics and Mechanics Faculty, St. Petersburg State University, 198504 Peterhof, Saint-Petersburg, Russia
2
National Research University Higher School of Economics, Soyusa Pechatnikov ul. 16, 190008 Saint-Petersburg, Russia
3
Department of Mathematical Information Technology, University of Jyväskylä, 40014 Jyväskylä, Finland
*
Author to whom correspondence should be addressed.
Entropy 2015, 17(7), 5101-5116; https://0-doi-org.brum.beds.ac.uk/10.3390/e17075101
Submission received: 26 May 2015 / Revised: 15 July 2015 / Accepted: 17 July 2015 / Published: 22 July 2015
(This article belongs to the Special Issue Recent Advances in Chaos Theory and Complex Networks)

Abstract

:
In applied investigations, the invariance of the Lyapunov dimension under a diffeomorphism is often used. However, in the case of irregular linearization, this fact was not strictly considered in the classical works. In the present work, the invariance of the Lyapunov dimension under diffeomorphism is demonstrated in the general case. This fact is used to obtain the analytic exact upper bound of the Lyapunov dimension of an attractor of the Shimizu–Morioka system.

1. Introduction

Two well-known definitions of Lyapunov exponents are the upper bounds of the exponential growth rate of the norms of linearized system solutions (LCEs) [1] and the upper bounds of the exponential growth rate of the singular values of the fundamental matrix of a linearized system (LEs) [2]. These definitions usually give the same values (e.g., for a “typical” regular system that satisfies the conditions of the Oseledets theorem). For a particular system, the LCEs and LEs may be different. Lyapunov [1] considered linearizations with bounded coefficients and introduced the class of regular linearizations: for regular linearization, the sum of LCEs equals the lower bound of the exponential growth rate of the normal fundamental matrix determinant (otherwise, the linearization is called irregular). Any regular linearization has exact LCEs (i.e., upper and lower bounds in the definition coincide), which are equal to LEs. In general, it can be shown that the ordered LEs are majorized by the ordered LCEs. More detailed discussions and examples can be found in [35]. Various characteristics of chaotic behavior are based on Lyapunov exponents. The sum of positive Lyapunov exponents is used [6,7] as the characteristic of the Kolmogorov–Sinai entropy rate [8,9]. Another measure of chaotic behavior is the Kaplan–Yorke or Lyapunov dimension, which is defined via the singular values or LEs [1012].
It is well known that a homeomorphism can change dimensions based on the space metric [13]. Thus, the question arises whether the Lyapunov dimension is invariant under changes (see, e.g., [14], “Is the Dimension of Chaotic Attractors Invariant under Coordinate Changes?”; and [1517]). Since the singular values and LEs are defined via the linearization of a system, the smooth changes of variables are of interest.
In the applied investigations, the invariance of the Lyapunov dimension under diffeomorphism is often used. However, in the case of irregular linearization, this fact is not strictly considered (see, e.g., the surveys [18,19]). In the present work, it is demonstrated in the general case that while singular values may not be invariant under changes of coordinates, the Lyapunov dimension of an attractor is invariant under diffeomorphism. This fact is used to obtain the exact formula of the Lyapunov dimension of an attractor of the Shimizu–Morioka system.

2. Lyapunov Dimension

Consider a continuously-differentiable mapping F defined on an open set U ⊂ ℝn. Denote by TxF the Jacobian matrix of the mapping F at the point x. The continuous differentiability of F implies that:
F ( x + h ) F ( x ) = ( T x F ) h + o ( h )
By assumption, KU is a bounded invariant set: F K = K.
Denote by αj(A) the singular numbers of (n×n)-matrix A, such that α1(A) ≥ α2(A) ≥ … ≥ αn(A). Consider the relation:
ω d ( A ) = α 1 ( A ) α 2 ( A ) α j ( A ) α j + 1 s ( A ) ,
where d = j + s, s ∈ [0, 1], j ∈ ℕ, j ∈ [1, n].
We also write a well-known Horn inequality:
ω d ( A B ) ω d ( A ) ω d ( B ) .
Definition 1. We say that a local Lyapunov dimension of the mapping F at the point xK is a number:
dim L ( F , x ) = j + s ,
where j is the largest integer number from [0, n], such that:
α 1 ( T x F ) α 2 ( T x F ) α j ( T x F ) 1
and a number value s ∈ [0, 1) is as follows:
α 1 ( T x F ) α 2 ( T x F ) α j ( T x F ) α j + 1 s ( T x F ) = 1.
By definition, dimL(F, x) = 0 if α1(TxF) < 1, and dim L(F, x) = n if α1(TxF)α2(TxF) … αn(TxF) 1.
Definition 2. [20,21] The Lyapunov dimension of the mapping F on the set K is a number:
dim L ( F , K ) = sup K dim L ( F , x ) .
Definition 3. A local Lyapunov dimension of a one-parameter group of mappings Ft at the point xK is a number:
dim L x = lim t sup + dim L ( F t , x ) .
Definition 4. The Lyapunov dimension of mappings Ft on the set K is a number:
dim L K = sup K dim L x .
Consider a differential equation:
d x d t = f ( x ) , x n
with a continuously-differentiable vector-function f(x). Suppose that for any initial value x0, there exists a solution of Equations (3)x(t, x0), defined on t ∈ [0, +). Here, x(0, x0) = x0.
Denote by Ft(x0) = x(t, x0) a shift operator along the solutions of Equations (3) and suppose that the set K ⊂ ℝn is bounded and invariant: FtK = K, ∀t ∈ ℝ1. TxFt is the Jacobian matrix of the mapping Ft at the point x.
Consider a differential equation:
d y d t = g ( y ) , y n
with a continuously-differentiable vector function g(y). Suppose that for any initial value y0, there exists a solution of Equations (4)y(t, y0), defined on t ∈ [0, +). Here, y(0, y0) = y0.
Denote by F ˜ t ( y 0 ) = y ( t , y 0 ) a shift operator along the solutions of Equations (4) and suppose that the set K ˜ n is bounded and invariant: F ˜ t K ˜ = K ˜, ∀t ∈ ℝ1. T y F ˜ t is the Jacobian matrix of the mapping F ˜ t at the point y.
Consider a transformation of coordinates of system (3) Q : xy:
y = Q ( x ) .
Suppose that Q is a diffeomorphism between the set K and the set K ˜, i.e., xK, y K ˜.
Hence, for a certain point y, there exists an inverse mapping Q1 : yx that is continuously differentiable in the neighborhood of y.
x = Q 1 ( y ) .
Under the above assumptions, we state the following assertions.
Lemma 1. If for T x F t sup x K ω d 0 ( T x F t ) < 1, then for the Lyapunov dimension of mapping Ft of the set K, the following estimate:
dim L ( F t , K ) j + s ,
is satisfied. Here, j + s = d < d0, s = sup x K i = 1 j μ i | μ j + 1 |, 0 < s1, j is the largest integer: j ∈ [1; n]; µi = LEi(x), LEi(x) are Lyapunov exponents (LEs), i = 1 j μ i > 0, µj+1 < 0.
Proof. Consider an arbitrary number s 0 : 0 < sup x K i = 1 j μ i | μ j + 1 | < s 0 1.
Recall that:
μ i = lim sup t 1 t ln | α i ( T x F t ) | i = 1 k μ i lim sup t 1 t ln ω t ( T x F t ) .
Here, ω k ( T x F t ) = α 1 ( T x F t ) α 2 ( T x F t ) α k ( T x F t ).
Using the well-known relation ω d ( T x F t ) = ω j 1 s ( T x F t ) ω j + 1 s ( T x F t ), we obtain:
lim sup t 1 t ln ( ω j 1 s 0 ( T x F t ) ω j + 1 s 0 ( T x F t ) ) lim sup t 1 t ln ( ω j 1 s 0 ( T x F t ) ) + lim sup t 1 t ln ( ω j + 1 s 0 ( T x F t ) ) ( 1 s 0 ) i = 1 j μ i + s 0 i = 1 j + 1 μ i = μ j + 1 ( s 0 i = 1 j μ i μ j + 1 ) < 0.
It follows that the inequality:
lim sup t 1 t ln ω d 0 ( T x F t ) < 0
is satisfied for:
ω d 0 ( T x F t ) < 1 , x K .
Thus, the assertion of the lemma sup x K ω d 0 ( T x F t ) < 1 holds true.
From the definition of the Lyapunov dimension, we obtain that for all d < d0, Relation (5) is satisfied. □
Theorem 1. The Lyapunov dimension of the mapping Ft of the set K is invariant with respect to the diffeomorphism Q : K K ˜. Namely
dim L ( F t , K ) = dim L ( F ˜ t , K ˜ ) .
Proof. For the mappings Ft and F ˜ t ( t > 0 ), Relation (1) is satisfied.
Then, the continuous differentiability of F ˜ t and the properties of diffeomorphism Q imply the following chain of relations:
T y 0 F ˜ t ( y ˜ y 0 ) + o ( h y ) = F ˜ t ( y ˜ ) F ˜ t ( y 0 ) = Q [ F t ( Q 1 ( y ˜ ) ) ] Q [ F t ( Q 1 ( y 0 ) ) ] = T y 0 Q ( y ˜ y 0 ) + o ( h y ) ,
where = y0 + hy, x ˜ = x 0 + h x, Q1(y0) = x0, Q 1 ( y ˜ ) = x ˜,
T y 0 Q = Q x | F ( x 0 ) T x F t | Q 1 ( y 0 ) Q 1 y | y 0 .
The continuous differentiability of the inverse mapping Q1 (property of diffeomorphism) implies the following relation:
Q 1 y | y 0 = ( Q x ) 1 | Q 1 ( y 0 ) = x 0 .
Taking into account (9) and (8) in Relation (7), we obtain:
T y 0 F ˜ t = Q x | F ( x 0 ) T x 0 F t ( Q x ) 1 | x 0 .
Thus, for the nonsingular matrix Q x ( d e t ( Q x ) 0 ), Relation (10) is satisfied.
Denote ω ¯ d ( T u F t ) = sup u K ω d ( T u F t ), Q x = Q x, ( Q x ) 1 = ( Q x ) 1. Consider the following relations:
T y F ˜ T = T y F ˜ k τ = ( T F ˜ ( k 1 ) τ F ˜ t ) ( T F ˜ ( k 2 ) τ F ˜ t ) ( T y F ˜ t ) ,
where τ = T k, k ∈ ℕ.
With provision for (10), Relation (11) can be rewritten as:
T y F ˜ k τ = ( Q F k τ T F ( k 1 ) τ F t ( Q F ( k 1 ) τ ) 1 ) ( Q F ( k 1 ) τ T F ( k 2 ) τ F t ( Q F ( k 2 ) τ ) 1 ) ( Q F ( x ) T y F t Q x 1 ) .
Applying Horn Inequality (2) to the last relation, we obtain the estimate:
ω d ( T y F ˜ k τ ) ω d ( Q F k τ ) ω d ( T F ( k 1 ) τ F t ) ω d ( T F ( k 2 ) τ F t ) ω d ( T y F t ) ω d ( Q x 1 ) .
This implies that:
ω ¯ d ( T y F ˜ k τ ) ω ¯ d ( Q F k τ ) [ ω ¯ d ( T x F t ) ] k ω ¯ d ( Q x 1 ) .
In the last inequality, the following estimates ω ¯ d ( Q x 1 ) m 1 , ω ¯ d ( Q F k τ ) m 2 , ω ¯ d ( T x F t ) < 1 , where m1, m2 ∈ ℝ +, are satisfied.
Repeating all of the above operations for the inverse transformation x = Q1(y), we obtain relations similar to (10):
T x 0 F t = ( Q x ) 1 | F ( x 0 ) T y 0 F ˜ t Q x | x 0
and (12):
ω ¯ d ( T x F k τ ) ω ¯ d ( ( Q F k τ ) 1 ) [ ω ¯ d ( T y F ˜ t ) ] k ω ¯ d ( Q x ) ,
where ω ¯ d ( Q x ) l 1, ω ¯ d ( ( Q F k τ ) 1 ) l 2 ω ¯ d ( T y F ˜ t ) < 1, l1, l2 ∈ ℝ+. Suppose that M = max(m1, m2, l1, l2).
Then, Relations (12) and (13) take the form:
ω ¯ d ( T y F ˜ k τ ) M [ ω ¯ d ( T x F t ) ] k .
and:
ω ¯ d ( T x F k τ ) M [ ω ¯ d ( T y F ˜ t ) ] k ,
respectively. Taking into account the assertion of the lemma, we obtain:
sup y K ˜ ω d ( T y F ˜ t ) | t 0 , lim t sup y K ˜ ω d ( T y F ˜ t ) = 0 ,
and:
sup x K ω d ( T x F t ) | t 0 , lim t sup x K ω d ( T x F t ) = 0 ,
as k → ∞.
Denote dim L ( F t , K ) = D 1, dim L ( F ˜ t , K ˜ ) = D 2.
Suppose that D1 < D2 and that d is an arbitrary number: D1 < d < D2. Consider the transformation Q. Then, from (14) and the definition of the Lyapunov dimension, it follows that D2 ≤ d, but this contradicts the initial assumption. Consequently, the inequality must be valid:
D 2 D 1 .
Suppose now that for inverse transformation Q1, D2 < d < D1. Then, from (15), it follows that D1 ≤ d. This is in contrast with the above assumption. Consequently:
D 1 D 2 .
Relations (16) and (17) imply the assertion of the theorem:
dim L ( F t , K ) = dim L ( F ˜ t , K ˜ ) .
In addition, this implies that:
lim t sup y K ˜ ω d ( T y F ˜ t ) = 0 lim t sup x K ω d ( T x F t ) = 0.
Thus, if Q is a diffeomorphism, then the change of coordinates of dynamical system y = Q(x) preserves the Lyapunov dimension of mapping Ft on compact sets.

3. Lyapunov Dimension of the Shimizu–Morioka System

Consider the Shimizu–Morioka system [22] of the form:
x ˙ = y , y ˙ = x λ y x z , z ˙ = α z + x 2 ,
where α, λ are the positive parameters.
Using the diffeomorphism:
x x , y y , z z x 2 2 ,
System (18) can be reduced to the following system:
x ˙ = y , y ˙ = x λ y x z + x 3 2 , z ˙ = α z + x y + ( 1 + α 2 ) x 2 ,
where α, λ are the positive parameters of System (18). We say that System (20) is a transformed Shimizu–Morioka system.
The following assertion is valid [2327].
Theorem 2. Suppose that for the integers j ∈ [1, n] and s ∈ [0, 1], there exist a continuously differentiable function υ(x) and a nonsingular matrix S, such that:
λ 1 ( x , S ) + λ 2 ( x , S ) + + λ j ( x , S ) + s λ j + 1 ( x , S ) + υ ˙ ( x ) < 0 , x n .
Then, dimL K ≤ j + s.
Here, υ ˙ = L f υ, υ(x) is a Lyapunov function.
Further, for System (20), the following will be proven.
Theorem 3. Suppose that K is a bounded invariant set of System (20): (0, 0, 0) ∈ K and that the following relations:
λ 4 10 + 3 α 13 α , λ < 1 α α , 4 λ 8 + 15 α 8 α 2 24 α 3 2 α ( α + 1 )
are satisfied. Then:
dim L K = 3 2 ( λ + α ) λ + 4 + λ 2 .
Proof. Consider a nonsingular matrix:
S = ( 1 k 0 0 λ α 1 0 0 0 1 ) ,
where k = 1 1 α ( α λ ).
Then:
S J S 1 = ( α λ 1 k 0 k ( 3 2 x 2 z ) 1 k α x k ( α x + y ) x α ) ,
J = ( 0 1 0 3 2 x 2 + 1 z λ x ( 2 + α ) x + y x α ) ,
where J is the Jacobian matrix of the right-hand side of System (20), a = 2 + 2α − λ.
The characteristic polynomial of the matrix:
1 2 ( S J S 1 + ( S J S 1 ) * )
takes the form:
( α + Λ ) [ Λ 2 + λ Λ k 2 16 ( 3 x 2 2 z + 4 k 2 ) 2 k 2 4 ( a x + y ) 2 1 1 k 2 ] .
Here, * denotes a transposition sign. This implies that the eigenvalues of Matrix (24) are the numbers:
λ 2 = α , λ 1 , 3 = 1 2 λ ± 1 2 w 1 ( x , y , z ) ,
where w 1 ( x , y , z ) = k 2 ( a x + y ) 2 + k 2 4 ( 3 x 2 2 z + 4 k 2 ) 2 + ( 2 α λ ) 2. Obviously, λ1 ≥ λ2 ≥ λ3.
We choose the Lyapunov function to be of the form:
υ ( x , y , z ) = 1 s 4 4 + λ 2 V ( x , y , z ) ,
where:
V ( x , y , z ) = μ 1 ( 2 y 2 2 x y x 4 + 2 x 2 z ) + μ 2 x 2 4 α z + μ 3 ( z 2 x 2 z + x 4 4 + x y ) + μ 4 ( z 2 + y 2 x 4 4 x 2 ) .
Here, µ1, µ2, µ3, µ4 are running parameters, s ∈ [0, 1).
Inequality (21) takes the form:
2 ( λ 1 + λ 2 + s λ 3 + υ ˙ ) = ( 1 + s ) λ 2 α + ( 1 s ) 4 + λ 2 + w ( x , y , z ) + ( 1 s ) V ˙ 2 4 + λ 2 < 0 ,
where w ( x , y , z ) = k 2 ( a x + y ) 2 + 6 x 2 + k 2 4 ( 3 x 2 2 z ) 2 4 z.
Using the inequality:
p + q p + q 2 p , p > 0 , p + q > 0 ,
we obtain the estimate:
2 ( λ 1 + λ 2 + s λ 3 + υ ˙ ) = ( 1 + s ) λ 2 α + ( 1 s ) 4 + λ 2 + 1 s 2 4 + λ 2 ( w ( x , y , z ) + V ˙ ) < 0.
We have:
w + V ˙ = B 1 x 2 + B 2 x y + B 3 y 2 + C 1 x 4 + C 2 x 2 z + C 3 z 2 ,
where:
B 1 = 4 4 α + k 2 a 2 2 μ 1 + μ 3 , B 2 = 2 a k 2 4 α + 2 ( 2 + λ ) μ 1 + 2 μ 2 λ μ 3 , B 3 = k 2 2 ( 2 λ + 1 ) μ 1 + μ 3 2 λ μ 4 , C 1 = 9 4 k 2 + ( α + 1 ) μ 1 1 2 ( α + 1 ) μ 3 , C 2 = 3 k 2 + 2 ( 1 α ) μ 1 + ( 2 α + 1 ) μ 3 + ( 2 + α ) μ 4 , C 3 = k 2 2 α μ 3 2 α μ 4 .
Condition:
w + V ˙ 0
is satisfied for all x, y, z, if and only if the following inequalities hold true:
B 1 0 , 4 B 1 B 3 B 2 2 0 ,
C 1 0 , 4 C 1 C 3 C 2 2 0.
We represent the second inequality from (27) in the form:
4 B 1 B 3 ( B ˜ 2 + 2 μ 2 ) 2 ,
where B ˜ 2 = 2 a k 2 4 α + 2 ( 2 + λ ) μ 1 λ μ 3.
Inequality (29) is solvable with respect to µ2 if B1B3 0. This implies that (27) is equivalent to the system of inequalities:
B 1 0 , B 3 0.
The second inequality in (30) is valid under certain µ4 > 0 if 2 μ 1 1 2 λ + 1 ( k 2 + μ 3 ). Then, from (30), we obtain:
2 μ 1 4 4 α + k 2 α 2 + μ 3 , 2 μ 1 1 2 λ + 1 ( k 2 + μ 3 ) .
Now, we represent the second inequality from (28) as:
4 C 1 ( C ˜ 3 2 α μ 4 ) [ C ˜ 2 + ( 2 + α ) μ 4 ] 2 0 ,
where:
C ˜ 2 = 3 k 2 + 2 ( 1 α ) µ 1 + ( 2 α + 1 ) µ 3 , C ˜ 3 = k 2 2 α µ 3 .
Inequality (32) is solvable with respect to µ4 if the discriminant of its left-hand side is nonnegative. This implies that:
4 α 2 ( 2 + α ) 2 C 1 2 + 2 α 2 + α C 1 C ˜ 2 + C 1 C ˜ 3 0.
Since C1 0, for the last inequality to be valid, it is sufficient that:
4 α 2 ( 2 + α ) 2 C 1 + 2 α 2 + α C ˜ 2 + C ˜ 3 0
is satisfied. The latter is equivalent to:
2 μ 1 μ 3 k 2 ( 1 α ) 2 α .
For the sake of convenience, we assume μ 1 ˜ = 2 μ 1 k 2, μ 3 ˜ = μ 3 k 2. Then, from the first inequality in (28) and Inequalities (31) and (33), we obtain a system of linear inequalities with respect to μ 1 ˜ and μ 3 ˜:
μ 1 ˜ 4 α 4 α k 2 + a 2 + μ 3 ˜ ,
μ 1 ˜ 1 2 λ + 1 + ( 1 + μ 3 ˜ ) ,
μ 1 ˜ 9 2 ( α + 1 ) + μ 3 ˜ ,
μ 1 ˜ ( 1 α ) 2 α + μ 3 ˜ .
Suppose μ 3 ^ = 1 2 λ + 1 4 α 4 α k 2 a 2 1 1 2 λ + 1. Then, from (36) and (37), we obtain:
μ 1 ˜ μ 3 ^ 9 2 ( α + 1 ) μ 1 ˜ μ 3 ^ ( 1 α ) 2 α .
This implies two cases:
1 ) 9 2 ( α + 1 ) ( 1 α ) 2 α , if 3 7 2 α 3 + 7 2 2 ) ( 1 α ) 2 α 9 2 ( α + 1 ) , if 0 < α 3 7 2 .
From (34) and (36), it follows that:
4 α 4 α k 2 + a 2 + 9 2 ( α + 1 ) 0.
Then:
4 8 + 15 α 8 α 2 24 α 3 2 α ( α + 1 ) λ 4 + 8 + 15 α 8 α 2 24 α 3 2 α ( α + 1 ) .
From (34) and (37), it follows that:
4 α 4 α k 2 + a 2 + ( 1 α ) 2 α 0.
Then:
4 10 + 3 α 13 α λ 4 + 10 + 3 α 13 α .
This means that (26) is valid if (38) and the inequality:
λ 4 + 10 + 3 α 13 α
are satisfied. Thus, this implies that Inequality (26) is satisfied if and only if System (34)(37) is solvable with respect to μ 1 ˜ and μ 3 ˜ = μ 3 ^, where μ 1 ˜ 0, μ 3 ^ 0, μ 1 ˜ μ 3 ^.
In turn, (25) is satisfied if (26) is valid, i.e., when λ < 1 α α and:
s > 4 + λ 2 ( λ + 2 α ) λ + 4 + λ 2 .
Thus, under Condition (22), we obtain the matrix S, the function υ(x, y, z) and a number s, satisfying (39), such that (25) is valid.
Consequently, dimL K ≤ 2 + s for all s, satisfying (39). This implies that:
dim L K 3 2 ( λ + α ) λ + 4 + λ 2 .
At the same time, considering the linearization of System (20) in a neighborhood of the points (0, 0, 0), we obtain:
dim L ( 0 , 0 , 0 ) = 3 2 ( λ + α ) λ + 4 + λ 2 .
This relation and (40) yield Relation (23). □
For α = 0.375, λ = 0.81, classical Shimizu–Morioka System (18) possesses a strange self-excited local attractor of the Lorenz type (Figure 1a); for α = 0.191450, λ = 0.81, it changes to the self-excited local attractor of the Burke and Shaw type (Figure 2a). Transformed System (20) demonstrates similar types of local attractors under the same values of parameters α and λ (Figure 1b and Figure 2b). Here, numerical estimations of the Lyapunov dimension by the code from [28] (2.0749, Figure 1b; 2.0059, Figure 2b) are consistent with the analytic exact upper bound for the Lyapunov dimension given by (40) (the corresponding values are 2.2014 and 2.3251, respectively).
Thus, using the diffeomorphism (19), which saves the Lyapunov dimension of an invariant bounded set K, it is possible to state that for classical Shimizu–Morioka System (18), the set K has the same Lyapunov dimension as System (20).
Remark 1. An oscillation can generally be easily numerically localized if the initial data from its open neighborhood in the phase space lead to a long-term behavior that approaches the oscillation. Therefore, from a computational perspective, it is natural to suggest the following classification of attractors [2932], which is based on the simplicity of finding their basins of attraction in the phase space: an attractor is called a self-excited attractor if its basin of attraction intersects with any open neighborhood of an equilibrium; otherwise, it is called a hidden attractor. For a self-excited attractor, its basin of attraction is connected to an unstable equilibrium, and therefore (standard computational procedure), self-excited attractors can be localized numerically by the standard computational procedure: by constructing a solution using initial data from an unstable manifold in a neighborhood of an unstable equilibrium, observing how it is attracted and visualizing the oscillation. In contrast, the basin of attraction for a hidden attractor is not connected to any equilibrium. For example, hidden attractors are attractors in systems with no equilibria or with only one stable equilibrium (a special case of the multi-stability: the coexistence of attractors in multi-stable systems). Well-known examples of the hidden oscillations are nested limit cycles in the 16th Hilbert problem (see, e.g., [32,33]) and counterexamples to the Aizerman and Kalman conjectures on the absolute stability of nonlinear control systems [32,3436].

4. Conclusion

In the present work, the invariance of the Lyapunov dimension under diffeomorphism is strictly proven. For the attractor of the Shimizu–Morioka system, a smooth change of coordinates is suggested and a Lyapunov function constructed, which allow one to obtain the analytic exact upper bound for the Lyapunov dimension. Similar results for some other Lorenz-like systems can be found, e.g., in [28,37,38].

Acknowledgments

This work was supported by the Russian Scientific Foundation (14-21-00041, Sections 1–2) and the Saint-Petersburg State University (6.38.505.2014, Section 3).

Author Contributions

Gennady Leonov, Tatyana Alexeeva and Nikolay Kuznetsov did the analytical estimates. Tatyana Alexeeva and Nikolay Kuznetsov made numerical experiments and wrote the paper. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; nor in the decision to publish the results

References

  1. Lyapunov, A.M. The General Problem of the Stability of Motion; CRC Press: Boca Raton, FL, USA, 1992. [Google Scholar]
  2. Oseledec, V.I. Multiplicative Ergodic Theorem: Characteristic Lyapunov Exponents of Dynamical Systems. Trans. Moscow Math. Soc. 1968, 19, 179–210. [Google Scholar]
  3. Leonov, G.A.; Kuznetsov, N.V. Time-Varying Linearization and the Perron Effects. Int. J. Bifurcat. Chaos. 2007, 17, 1079–1107. [Google Scholar]
  4. Kuznetsov, N.V.; Mokaev, T.N.; Vasilyev, P.A. Numerical Justification of Leonov Conjecture on Lyapunov Dimension of Rossler Attractor. Commun. Nonlinear Sci. Numer. Simulat. 2014, 19, 1027–1034. [Google Scholar]
  5. Kuznetsov, N.V.; Alexeeva, T.; Leonov, G.A. Invariance of Lyapunov Characteristic Exponents, Lyapunov Exponents, and Lyapunov Dimension for Regular and Non-regular Linearizations 2014. arXiv:1410.2016v2.
  6. Millionschikov, V.M. A Formula for the Entropy of Smooth Dynamical Systems. Differencial’nye Uravenija 1976, 12, 2188–2192. [Google Scholar]
  7. Pesin, Y. Characteristic Lyapunov Exponents and Smooth Ergodic Theory. Russ. Math. Surv. 1977, 32, 55–114. [Google Scholar]
  8. Kolmogorov, A. On Entropy per Unit Time as a Metric Invariant of Automorphisms. Dokl. Akad. Nauk SSSR 1959, 124, 754–755. [Google Scholar]
  9. Sinai, Y. On the Notion of Entropy of Dynamical Systems. Dokl. Akad. Nauk SSSR 1959, 124, 768–771. [Google Scholar]
  10. Kaplan, J.L.; Yorke, J.A. Chaotic Behavior of Multidimensional Difference Equations. In Functional Differential Equations and Approximations of Fixed Points; Springer: Berlin, Germany, 1979; pp. 204–227. [Google Scholar]
  11. Douady, A.; Oesterle, J. Dimension de Hausdorff des Attracteurs. CR Acad. Sci. Paris 1980, 290, 1135–1138. [Google Scholar]
  12. Ledrappier, F. Some Relations between Dimension and Lyapounov Exponents. Commun. Math. Phys. 1981, 81, 229–238. [Google Scholar]
  13. Pontrjagin, L.; Schnirelmann, L. Sur une Propriété Métrique de la Dimension. Ann. Math. 1932, 156–162. [Google Scholar]
  14. Ott, E.; Withers, W.; Yorke, J. Is the Dimension of Chaotic Attractors Invariant under Coordinate Changes? J. Stat. Phys. 1984, 36, 687–697. [Google Scholar]
  15. Dettmann, C.; Frankel, N.; Cornish, N. Chaos and Fractals around Black Holes. Fractals 1995, 3, 161–181. [Google Scholar]
  16. Hoover, W.; Hoover, C. Why Instantaneous Values of the “Covarian” Lyapunov Exponents Depend upon the Chosen State-Space Scale. Fractals 2014, 20, 5–8. [Google Scholar]
  17. Sprott, J.; Hoover, W.; Hoover, C. Heat Conduction, and the Lack Thereof, in Time-Reversible Dynamical Systems: Generalized Nosé-Hoover Oscillators with a Temperature Gradient. Phys. Rev. E 2014, 89, 042914. [Google Scholar]
  18. Gelfert, K.; Motter, A. (Non)Invariance of Dynamical Quantities for Orbit Equivalent Flows. Commun. Math. Phys. 2010, 300, 411–433. [Google Scholar]
  19. Barreira, L.; Gelfert, K. Dimension Estimates in Smooth Dynamics: A Survey of Recent Results. Ergod. Theor. Dyn. Syst. 2011, 31, 641–671. [Google Scholar]
  20. Temam, R. Infinite-Dimensional Dynamical Systems in Mechanics and Physics, 2nd ed; Springer: New York, NY, USA, 1997. [Google Scholar]
  21. Hunt, B. Maximum Local Lyapunov Dimension Bounds the Box Dimension of Chaotic Attractors. Nonlinearity 1996, 9, 845–853. [Google Scholar]
  22. Shimizu, T.; Morioka, N. On the Bifurcation of a Symmetric Limit Cycle to an Asymmetric One in a Simple Model. Phys. Lett. A 1980, 76, 201–204. [Google Scholar]
  23. Leonov, G.A. On Estimations of the Hausdorff Dimension of Attractors. Vestnik Leningradskogo Gos. Univ. Ser. 1991, 1, 41–44. [Google Scholar]
  24. Leonov, G.A.; Boichenko, V.A. Lyapunov’s Direct Method in the Estimation of the Hausdorff Dimension of Attractors. Acta Appl. Math. 1992, 26, 1–60. [Google Scholar]
  25. Boichenko, V.A.; Leonov, G.A.; Reitmann, V. Dimension Theory for Ordinary Differential Equations; Vieweg Teubner Verlag: Wiesbaden, Germany, 2005. [Google Scholar]
  26. Leonov, G.A. Strange Attractors and Classical Stability Theory. Nonlinear Dyn. Syst. Theor. 2008, 8, 49–96. [Google Scholar]
  27. Leonov, G.A. Lyapunov Functions in the Attractors Dimension Theory. J. Appl. Math. Mech. 2012, 76, 129–141. [Google Scholar]
  28. Leonov, G.; Kuznetsov, N.; Mokaev, T. Homoclinic Orbits, and Self-excited and Hidden Attractors in a Lorenz-Like System Describing Convective Fluid Motion. Eur. Phys. J. Special Topics. 2015, 224, 1421–1458. [Google Scholar]
  29. Kuznetsov, N.V.; Leonov, G.A.; Vagaitsev, V.I. Analytical-Numerical Method for Attractor Localization of Generalized Chua’s System. IFAC Proc. Vol. 2010, 4, 29–33. [Google Scholar]
  30. Leonov, G.A.; Kuznetsov, N.V.; Vagaitsev, V.I. Localization of Hidden Chua’s Attractors. Phys. Lett. A 2011, 375, 2230–2233. [Google Scholar]
  31. Leonov, G.A.; Kuznetsov, N.V.; Vagaitsev, V.I. Hidden Attractor in Smooth Chua Systems. Physica D 2012, 241, 1482–1486. [Google Scholar]
  32. Leonov, G.A.; Kuznetsov, N.V. Hidden Attractors in Dynamical Systems. From Hidden Oscillations in Hilbert-Kolmogorov, Aizerman, and Kalman Problems to Hidden Chaotic Attractors in Chua Circuits. Int. J. Bifurcat. Chaos. 2013, 23, 1330002. [Google Scholar]
  33. Kuznetsov, N.V.; Kuznetsova, O.A.; Leonov, G.A. Visualization of Four Normal Size Limit Cycles in Two-Dimensional Polynomial Quadratic System. Differ. Equ. Dyn. Syst. 2013, 21, 29–34. [Google Scholar]
  34. Leonov, G.A.; Bragin, V.O.; Kuznetsov, N.V. Algorithm for Constructing Counterexamples to the Kalman Problem. Dokl. Math. 2010, 82, 540–542. [Google Scholar]
  35. Leonov, G.A.; Kuznetsov, N.V. Algorithms for Searching for Hidden Oscillations in the Aizerman and Kalman Problems. Dokl. Math. 2011, 84, 475–481. [Google Scholar]
  36. Bragin, V.O.; Vagaitsev, V.I.; Kuznetsov, N.V.; Leonov, G.A. Algorithms for Finding Hidden Oscillations in Nonlinear Systems. The Aizerman and Kalman Conjectures and Chua’s Circuits. J. Comput. Syst. Sci. Int. 2011, 50, 511–543. [Google Scholar]
  37. Leonov, G.A.; Kuznetsov, N.V. On Differences and Similarities in the Analysis of Lorenz, Chen, and Lu Systems. Appl. Math. Comput. 2015, 256, 334–343. [Google Scholar]
  38. Leonov, G.; Kuznetsov, N.; Korzhemanova, N.; Kusakin, D. Estimation of Lyapunov Dimension for the Chen and Lu Systems 2015. arxiv:1504.04726.
Figure 1. Self-excited local attractors in: (a) System (18) for α = 0.375, λ = 0.81; (b) System (20) for α = 0.375, λ = 0.81.
Figure 1. Self-excited local attractors in: (a) System (18) for α = 0.375, λ = 0.81; (b) System (20) for α = 0.375, λ = 0.81.
Entropy 17 05101f1
Figure 2. Self-excited local attractors in: (a) System (18) for α = 0.191450, λ = 0.81, “Burke and Shaw-like”; (b) System (20) for α = 0.191450, λ = 0.81, “Burke and Shaw-like”.
Figure 2. Self-excited local attractors in: (a) System (18) for α = 0.191450, λ = 0.81, “Burke and Shaw-like”; (b) System (20) for α = 0.191450, λ = 0.81, “Burke and Shaw-like”.
Entropy 17 05101f2

Share and Cite

MDPI and ACS Style

Leonov, G.A.; Alexeeva, T.A.; Kuznetsov, N.V. Analytic Exact Upper Bound for the Lyapunov Dimension of the Shimizu–Morioka System. Entropy 2015, 17, 5101-5116. https://0-doi-org.brum.beds.ac.uk/10.3390/e17075101

AMA Style

Leonov GA, Alexeeva TA, Kuznetsov NV. Analytic Exact Upper Bound for the Lyapunov Dimension of the Shimizu–Morioka System. Entropy. 2015; 17(7):5101-5116. https://0-doi-org.brum.beds.ac.uk/10.3390/e17075101

Chicago/Turabian Style

Leonov, Gennady A., Tatyana A. Alexeeva, and Nikolay V. Kuznetsov. 2015. "Analytic Exact Upper Bound for the Lyapunov Dimension of the Shimizu–Morioka System" Entropy 17, no. 7: 5101-5116. https://0-doi-org.brum.beds.ac.uk/10.3390/e17075101

Article Metrics

Back to TopTop