Next Article in Journal
Entropy of Shortest Distance (ESD) as Pore Detector and Pore-Shape Classifier
Next Article in Special Issue
A Method for Choosing an Initial Time Eigenstate in Classical and Quantum Systems
Previous Article in Journal
Quantum Contextuality with Stabilizer States
Previous Article in Special Issue
A Unification between Dynamical System Theory and Thermodynamics Involving an Energy, Mass, and Entropy State Space Formalism
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Increase in Switching Systems

1
Centro de Investigación Operativa, Universidad Miguel Hernández. Avda. de la Universidad s/n, Elche 03202, Spain
2
Fachbereich Mathematik, Johan Wolfgang Goethe Universität. Frankfurt am Main 60054, Germany
*
Author to whom correspondence should be addressed.
Entropy 2013, 15(6), 2363-2383; https://doi.org/10.3390/e15062363
Submission received: 10 May 2013 / Revised: 3 June 2013 / Accepted: 3 June 2013 / Published: 7 June 2013
(This article belongs to the Special Issue Dynamical Systems)

Abstract

:
The relation between the complexity of a time-switched dynamics and the complexity of its control sequence depends critically on the concept of a non-autonomous pullback attractor. For instance, the switched dynamics associated with scalar dissipative affine maps has a pullback attractor consisting of singleton component sets. This entails that the complexity of the control sequence and switched dynamics, as quantified by the topological entropy, coincide. In this paper we extend the previous framework to pullback attractors with nontrivial components sets in order to gain further insights in that relation. This calls, in particular, for distinguishing two distinct contributions to the complexity of the switched dynamics. One proceeds from trajectory segments connecting different component sets of the attractor; the other contribution proceeds from trajectory segments within the component sets. We call them “macroscopic” and “microscopic” complexity, respectively, because only the first one can be measured by our analytical tools. As a result of this picture, we obtain sufficient conditions for a switching system to be more complex than its unswitched subsystems, i.e., a complexity analogue of Parrondo’s paradox.

1. Introduction

The time-switched dynamics of two one-dimensional, dissipative affine maps,
f ± ( x ) = θ ± x ± 1 0 < θ - , θ + < 1 , θ - θ +
was studied by the authors in [1], where they showed that the topological entropy of the resulting non-autonomous dynamical system coincides with the topological entropy of the autonomous dynamical system generating the switching (or control) sequence. In other words, the switched dynamics cannot be more complex than the switching sequence generator.
The situation envisaged in [1] is special in several regards.
(i)
The state space is R .
(ii)
Both the forward and backward dynamics are amenable to detailed analysis.
(iii)
The non-autonomous pullback attractor consists of singleton component sets.
In this follow-up paper we continue studying this question in an extended setting, namely, we consider this time switching between difference equations (called constituent maps) that have pullback attractors with nontrivial component sets. Our scope is to separate the particular results from the general ones, and so better understand the intricacies of switching and non-autonomous dynamics [2,3,4]. Specifically, in this paper:
(i’)
The state space is R d or a proper subset of it.
(ii’)
The constituent maps are not specified, except for the fact that they are supposed to have attractors.
(iii’)
The component sets of the pullback attractor of the systems under switching are supposed to be uniformly bounded.
Both from an instrumental and a conceptual point of view, the main challenge is introduced by the generalization (iii’). To begin with, working with attractors with nontrivial component sets instead of point-valued components sets requires using the Hausdorff distance. Contrarily to what happens in the latter case, the Hausdorff distance between component sets of a set-valued attractor is not continuous, in general, but only upper-continuous when a control sequence converges to another one [2]. This technical shortcoming is behind Assumption 1 in Section 3 (“Results using measurability”), and Assumption 3 in Section 5 (“Results using continuity”). Also the generality of (ii’) will be limited in different ways (injectivity of a certain map Φ in Section 3, or Assumption 2 in Section 4) in order to derive sharper results.
On the conceptual side, new features related to the extended geometry of set-valued attractors manifest themselves. Indeed, the observable complexity is going to be the result of two contributions. One, the “macroscopic” complexity, comes from the trajectory segments connecting different component sets of the pullback attractor (think that each component set has been lumped to a point). The other one, the “microscopic” complexity, comes from trajectory segments within the component sets. No geometrical meaning should be attached to this denomination since the (in general, uncountably many) component sets might be packed in a state space region of a size comparable to the attractor itself. The rationale for distinguishing two sorts of complexity is rather that our analytical tools will only be able to resolve the trajectories up to the precision set by the component sets. An offshoot of this picture is that the entropy of the coarse-grained dynamics (or “macroscopic entropy”) is a lower bound of the entropy of the switched dynamics. This fact will again allow us to relate the entropy of trajectories and control sequences, once we have shown that the latter coincide with the macroscopic entropy.
Thus, the main purpose of the present paper is the study of the macroscopic complexity in switched dynamics, as measured by the macroscopic topological entropy. In Section 3 and 5 we prove that, under certain provisos, the macroscopic topological entropy coincides with the topological entropy of the switching sequence generator. In this case, the complexity of the switched dynamics (a mixture of macroscopic and microscopic complexities) is certainly greater than the complexity of the switching sequence. This result can be brought in connection with Parrondo’s paradox, i.e., the emergence of new properties via switching [5,6,7,8]. Indeed, it suffices that the entropy of the switching sequence generator is higher than the entropies of the constituent maps for the switched dynamics to be more complex than the constituent dynamics.
This paper is organized as follows. Section 2 sets the mathematical framework (switching systems, Hausdorff distance and pullback attractors) of the paper. Additional materials on the Hausdorff distance have been collected in the Appendix for the reader’s convenience. A few results for further reference, but also interesting on their own, are proved as well in this section and in the Appendix. The main theoretical results of the paper are derived in Section 3 (Theorem 3) and Section 5 (Theorem 4) depending on whether the switched dynamics complies with Assumption 1 (Section 3), or with Assumption 3 (Section 5). In between, Section 4 scrutinizes a property needed in Theorem 3 and reformulates it as Assumption 2 to be used in the proof of Theorem 4. Some numerical simulations illustrating our theoretical results make up Section 6. The highlights of all these sections are summarized in the Conclusion, followed by the references and the Appendix.

2. Preliminaries

This section introduces the background on switching and non-autonomous dynamical systems needed to make this paper self-contained. The interested readers are referred to the books [2,3], the review [4], and the papers [9,10,11,12,13].

2.1. Switching Systems

Let s = ( , s - 2 , s - 1 , s 0 , s 1 , s 2 , ) S : = { - 1 , 1 } Z be the collection of all switching controls, which is a compact metric space with the metric
dist S ( s , s ) = n Z 2 - n s n - s n
In calculations though, the following, equivalent distance might be more convenient.
dist S ( s , s ) = 0 if s n = s n for all n Z , 2 - N if s n = s n for n < N , while s N s N and/or s - N s - N
The metric space ( S , dist S ) is compact; see, e.g., [14] (p. 207).
Consider the left shift operator σ defined on S ,
σ : ( , s n , s n + 1 , ) ( , s n + 1 , s n + 2 , )
Then σ is continuous with respect to the product topology on S , i.e., generates an autonomous dynamical system (group under composition) on S .
We consider discrete time dynamical systems generated by switching between two continuous mappings f ± 1 : R d R d , where the autonomous dynamical systems generated by each mapping are dissipative and have global attractors A ± . In other words, the (time) switched dynamics (or switching system) we are going to study in the following is described by the non-autonomous difference equation
x n + 1 = f s n ( x n )
in R d for different sequences s = ( s n ) n Z S . Sometimes one says that f ± 1 are the constituent maps of the switching system (4). For brevity, we shorten f ± 1 to f ± , x ± 1 to x ± , etc.
As an important example, let f ± 1 be affine maps,
f ± ( x ) = M ± x + b ±
where M ± are d × d -matrices, and x , b ± are column vectors. Linear systems correspond to b ± = 0 . Then
x n + k = i = 0 k - 1 M s n + i x n + i = 1 k - 1 M s n + i b s n + i = 2 k - 1 M s n + i b s n + 1 + M s n + k - 1 b s n + k - 2 + b s n + k - 1
for k 1 . The order in the matrix products i = 0 k - 1 M s n + i , etc., is right to left as the index i increases, i.e.,
i = 0 k - 1 M s n + i = M s n + k - 1 M s n + 1 M s n
Likewise if the matrices M ± are invertible, then from
x n - 1 = f s n - 1 - 1 ( x n ) = M s n - 1 - 1 ( x n - b s n - 1 )
we get
x n - k = i = 1 k M s n - i - 1 x n - i = 1 k M s n - i - 1 b s n - 1 - i = 2 k M s n - i - 1 b s n - 2 - M s n - k - 1 b s n - k
The same rule as in Equation (7) applies here to the products of inverse matrices.
As usual, let x denote the norm of x R d and M a compatible norm of the d × d matrix M.
Proposition 1
Consider the affine switched dynamics Equation (4) and Equation (5). Furthermore, suppose that μ : = max { M - , M + } < 1 , and β : = max { b + , b - } . Then
lim sup n x n β 1 - μ
for any x 0 R .
Proof.
From Equation (6), we have
x n + k μ k x n + ( μ k - 1 + μ k - 2 + + μ + 1 ) β = μ k x n + 1 - μ k 1 - μ β
where μ < 1 . Let now k . □
Linear systems with M - , M + < 1 are globally dissipative. The origin is a fixed point and, according to Equation (10) with β = 0 , all other orbits converge to it for each s S .

2.2. Hausdorff Metric and Sequences of Compact Subsets

Let K be the space of nonempty compact subsets of R d , which is a complete metric space with the Hausdorff metric
dist H ( A , B ) : = max { ρ ( A , B ) , ρ ( B , A ) }
where ρ ( A , B ) is the Hausdorff semi-distance defined by
ρ ( A , B ) : = max a A dist ( a , B ) dist ( a , B ) : = min b B a - b
At variance with dist H ( A , B ) , ρ ( A , B ) is not a metric because it is not always symmetric, and ρ ( A , B ) = 0 only implies A B . Intuitively speaking, for the distance between A and B to be small, both sets have to almost overlap, i.e., the difference set A Δ B = A B B A has to be “small”. Further results on the Hausdorff metric can be found in the Appendix.
Let B R ( a ) and B ¯ R ( a ) be the open and closed balls of radius R 0 and center a R d , respectively. Define K R = K B ¯ R ( 0 ) , i.e., the family of nonempty compact subsets of R d that are contained in B ¯ R ( 0 ) . Then, K R , dist H is a compact metric space.
The next proposition will be needed in Section 5. It is a particular case of Lemma A2, stated and proved in the Appendix.
Proposition 2
Let f : B ¯ R ( 0 ) R d be continuous, and A , B K R . Then
dist H ( f ( A ) , f ( B ) ) ω ( dist H ( A , B ) ) ,
where ω ( · ) is a continuous function with ω ( 0 ) = 0 .
In order to apply Proposition 2 to our constituent maps f ± , whose definition domain is in principle all R d , it suffices to consider their restrictions to B ¯ R ( 0 ) .
Further, we define a metric on the space K R : = K R Z of bi-infinite sequences A : = ( A n ) n Z of nonempty compact subsets of R d in B ¯ R ( 0 ) by
dist K R ( A , A ) : = n Z 2 - n dist H ( A n , A n )
for A = ( A n ) n Z , and A = ( A n ) n Z .

2.3. Skew Product Flows and Pullback Attractors

Define φ: N 0 × S × R d R d by φ ( 0 , s , x 0 ) = x 0 and
φ ( n , s , x 0 ) : = f s n - 1 f s 1 f s 0 ( x 0 ) n 1
Then φ is a continuous cocycle mapping with respect to σ, i.e.,
φ ( n + k , s , x 0 ) = φ ( n , σ k s , φ ( k , s , x 0 ) )
for all n , k 0 , and ( σ , φ ) is a discrete time skew product flow (non-autonomous dynamical system) on S × R d .
An entire solution of a discrete time skew product flow ( σ , φ ) is a mapping χ : S R d such that
χ ( σ n s ) = φ ( n - k , σ k s , χ ( σ k s ) )
for all s S and n , k Z with k n . In particular,
χ ( σ n s ) = φ ( n , s , χ ( s ) ) for all n 0
A pullback attractor is a family of nonempty compact subsets, A = { A ( s ) , s S } K which is φ-invariant, i.e.,
φ ( n , s , A ( s ) ) = A ( σ n s ) n 0
and pullback attracts, i.e.,
dist H φ ( n , σ - n s , D ) , A ( s ) 0 f o r n
for every nonempty bounded subset D R d . The A ( s ) are called the component sets of the attractor A .
We will assume that the switched dynamics (4) has a pullback attractor A = { A ( s ) : s S } such that the A ( s ) are nonempty, uniformly bounded compact subsets of R d . This means that there is R > 0 such that A ( s ) B ¯ R ( 0 ) for every s S . If the constituent maps are both affine, then it is inferred from Proposition 1 that one may choose R = β 1 - μ . A sufficient condition ensuring the existence of such a pullback attractor is that the two unswitched systems have a common bounded, positively invariant absorbing set.
Proposition 3
[2](Theorem 3.34) The set-valued mapping s A ( s ) is upper semi-continuous in ( K R , dist H ) , i.e.,
ρ A ( s ) , A ( s * ) 0 a s dist S ( s , s * ) 0
Counterexamples show that, in general, we cannot replace the Hausdorff semi-distance here by the Hausdorff metric, but in special cases we can do that, e.g., when the pullback attractor consists of singleton sets as in [1] (see Example 1 below), since then obviously
ρ A ( s ) , A ( s * ) = dist H A ( s ) , A ( s * ) = A ( s ) - A ( s * )
See Case 2 of the numerical simulations (Section 4.2) for other example.
A general, sufficient condition for the continuity of the map s A ( s ) in ( K R , dist H ) is provided by the following result.
Proposition 4
Suppose that dist H φ ( n , σ - n s , D ) , A ( s ) 0 uniformly in s for some nonempty, bounded set D R d . Then the map s A ( s ) is continuous in ( K R , dist H ) , i.e.,
dist H A ( s ) , A ( s * ) 0 a s dist S ( s , s * ) 0
Proof.
By the property (18) of pullback attractors, for all ε > 0 there exists an N such that
dist H φ ( n , σ - n s , D ) , A ( s ) < ε 2 dist H φ ( n , σ - n s * , D , A ( s * ) ) < ε 2
for all n N , where D R d is bounded and fixed for the time being. By hypothesis, N depends on ε but not on either s or s * . The triangle inequality then yields
dist H ( A ( s ) , A ( s * ) ) dist H A ( s ) , φ ( N , σ - N s , D ) + dist H φ ( N , σ - N s , D ) , φ ( N , σ - N s * , D + dist H φ ( N , σ - N s * , D , A ( s * ) ) < ε + dist H φ ( N , σ - N s , D ) , φ ( N , σ - N s * , D
where
φ ( N , σ - N s , D ) = f s - 1 f s - N + 1 f s - N ( D ) φ ( N , σ - N s * , D ) = f s - 1 * f s - N + 1 * f s - N * ( D )
Take now dist S ( s , s * ) < 2 - ( N + 1 ) so as s k = s k * for 0 ≤ k N , hence φ ( N , σ - N s , D ) = φ ( N , σ - N s * , D ) , and
dist H φ ( N , σ - N s , D ) , φ ( N , σ - N s * , D = 0
We conclude from Equation (20) that if dist S ( s , s * ) < 2 - ( N + 1 ) [where N = N ( ε ) is such that Equation (19) holds], then
dist H ( A ( s ) , A ( s * ) < ε
for all ε > 0 . □
Remark 1
It follows from Proposition 3 that the pullback attractor is also a forward attractor, i.e., with
dist H φ ( n , s , D ) , A ( σ n s 0 as n

3. Results Using Measurability

Remember that K R = K R Z is endowed with the metric (13). Let Φ: S K R be the map defined by Φ ( s ) = ( A ( σ n s ) ) n Z . We show first that Φ is Borel measurable.
By definition the Borel sigma-algebras of the product spaces S and K R are generated by the corresponding cylinder sets or just cylinders. If p , q Z , p q , the cylinders of S = { - 1 , + 1 } Z have the form
[ a p , . . . , a q ] p q = { ( s n ) n Z S : s n = a n for p l q }
where a p , . . . , a q { - 1 , + 1 } , and the cylinders of K R have the form
[ K p , . . . , K q ] p q = { ( A n ) n Z K R : A n = K n for p n q }
where K p , . . . K q K R . If p = q , then we simplify the notation to [ a p ] p or [ a ] p , and analogously for the cylinders of K R .
If, as before, σ is the (left) shift on S , and Σ is the shift on K R , then σ n [ a p , . . . , a q ] p q = [ a p , . . . , a q ] p - n q - n and, similarly, Σ n [ K p , . . . , K q ] p q = [ K p , . . . , K q ] p - n q - n for all n Z .
Lemma 1
Φ : S K R is Borel measurable.
Proof.
Since the cylinder sets build a semi-algebra of the product sigma-algebras they generate, it suffices to prove that Φ - 1 [ K p , . . . , K q ] p q is Borel measurable for every p , q Z , p q , and every K p , . . . K q K R (see Theorem 1.1 in [15]). Note that Φ - 1 [ K p , . . . , K q ] p q = unless [see Equation (17)] K p = A ( s ) for some s S , and K p + n = A ( σ n s ) = φ ( n , s , K p ) for n = 1 , . . . , q - p . In any case,
Φ - 1 [ K p , . . . , K q ] p q = Φ - 1 ( [ K p ] p [ K p + 1 ] p + 1   . . .   [ K q ] q ) = Φ - 1 [ K p ] p Φ - 1 [ K p + 1 ] p + 1   . . . .   Φ - 1 [ K q ] q
Since the map s A ( s ) is upper semi-continuous (hence Borel measurable), the sets Φ - 1 [ K p + n ] p + n , 0 n q - p , are Borel measurable and so is their intersection as well. □
Suppose now that Φ is one-to-one and let us explore when Φ - 1 is also Borel measurable, resulting in a Borel bimeasurable mapping. For p , q Z , p q ,
Φ [ a p , . . . , a q ] p q = ( A ( σ n [ a p , . . . , a q ] p q ) ) n Z
where
A ( σ n [ a p , . . . , a q ] p q ) = { A ( s ) : s σ n [ a p , . . . , a q ] p q } = { A ( s ) : s [ a p , . . . , a q ] p - n q - n } = { A ( s ) : s [ a p ] p - n [ a p + 1 ] p - n + 1   . . .   [ a q ] q - n } = { A ( s ) : s [ a p ] p - n } . . . { A ( s ) : s [ a q ] q - n } = A ( [ a p ] p - n ) . . .   A ( [ a q ] q - n )
Note that
A ( [ a ] p ) = { A ( s ) : s p = a }
is an uncountable union of compact sets.
Assumption 1.
The sets A ( [ + 1 ] p ) and A ( [ - 1 ] p ) are Borel measurable for every p Z .
Continuity or closedness of the mapping s A ( s ) are obvious sufficient conditions for Assumption 1 to hold. (As a matter of fact, by the closed map lemma, if s A ( s ) is continuous then it is closed because S is compact and K R is a Hausdorff topological space.) For example, in the affine, one-dimensional case studied in [1] A ( s ) is a singleton, i.e., A ( s ) = { χ ( s ) } , for all s S . That all the sets A ( [ a ] p ) are compact (hence Borel measurable) follows in this case from the continuity of s χ ( s ) . In Section 5 we will study with more detail the consequences of assuming the mapping s A ( s ) continuous.
Lemma 2
If Φ : s ( A ( σ n s ) ) n Z is one-to-one and Assumption 1 holds, then Φ is Borel bimeasurable.
Proof.
We only need to prove that Φ transforms cylinders into Borel measurable sets. In view of Equation (21) this boils down to showing that A ( σ n [ a p , . . . , a q ] p q ) is Borel measurable for every p , q Z , p q .
From Equation (22) and Assumption 1, it follows that A ( σ n [ a p , . . . , a q ] p q ) is a finite intersection of Borel measurable sets, hence it is Borel measurable. □
Consider the diagram
S σ S Φ Φ K R Σ K R
then it is straightforward to check that this diagram commutes, i.e., Φ σ = Σ Φ . In order to derive that h ( σ ) = h ( Σ Φ ( S ) ) , it suffices that Φ : S Φ ( S ) is a bimeasurable bijection (see Corollary 8.6.1 (iv) in [15]). This being the case, Lemma 3.2 yields the following result.
Theorem 2
Let Φ : S Φ ( S ) be a bijection and suppose that Assumption 1 holds. Then h ( σ ) = h ( Σ Φ ( S ) ) .
Other conditions leading also to h ( σ ) = h ( Σ Φ ( S ) ) will be discussed in Section 5.
Notice that h ( Σ Φ ( S ) ) corresponds to what we called macroscopic entropy in the Introduction. Indeed, since (i) Φ ( s ) is a trajectory in A , a set whose “points” are the component sets A ( s ) of A ; and (ii) χ ( s ) A ( s ) implies χ ( σ n s ) A ( σ n s ) for all n Z , it holds that h ( Σ Φ ( S ) ) measures the complexity of the trajectories ( χ ( σ n s ) ) n Z up to the precision set by the distinct component sets A ( σ n s ) . To relate the macroscopic entropy, h m a c r o h ( Σ Φ ( S ) ) , to the entropy of the switched dynamics (referred to as microscopic entropy in the Introduction), h m i c r o , let E K R be the set of entire orbits, and consider the commutative diagram
E Σ E i i Φ ( S ) Σ Φ ( S )
where i is the inclusion ( χ ( σ n s ) ) n Z ( A ( σ n s ) ) n Z . Thus Σ Φ ( S ) is a factor of (or semi-conjugate to) Σ E on account of i being measure-preserving ( i - 1 B = B for all measurable B Φ ( S ) ) and onto. It follows that h ( Σ E ) h ( Σ Φ ( S ) ) h m a c r o , i.e., h m a c r o is a lower bound of the topological entropy of the switched dynamics h m i c r o h ( Σ E ) . To keep with the physical identification of entropies, set h ( σ ) h c o n t r o l for the topological entropy (or complexity) of the control sequence generator. From
h m a c r o h m i c r o
and Theorem 2,
h c o n t r o l = h m a c r o ,
we conclude the following result on the relation between the topological entropies of the control sequence generator, h ( σ ) , and the ensuing switched dynamics, h m i c r o .
Theorem 3
Under the hypotheses of Theorem 2, h c o n t r o l h m i c r o .
This theorem provides sufficient conditions for a complexity version of Parrondo’s paradox. Indeed, it suffices that h ( f + ) , h ( f - ) < h ( σ ) (along with the assumptions of Theorem 2), where h ( f ± ) is the topological entropy of f ± , for the complexity of the switched dynamics to exceed the complexity of the constituent maps.
In the special case of pullback attractors consisting of singletons, like in [1], the inclusion i in the diagram (24) becomes the identity, a trivial isomorphism, and h m a c r o = h m i c r o .

4. Conditions for the Injectivity of Φ

According to Theorem 2, the injectivity of Φ (along with Assumption 1) is instrumental for h ( σ ) and h ( Σ ) to coincide. We explore next sufficient conditions for the injectivity of Φ.
If Φ is not injective, then there are control sequences s s * such that Φ ( s ) = Φ ( s * ) , i.e., A ( σ n s ) = A ( σ n s * ) for all n Z . Then
f s n ( A ( σ n s ) ) = A ( σ n + 1 s ) = A ( σ n + 1 s * ) = f s n * ( A ( σ n s * ) ) = f s n * ( A ( σ n s ) )
and, similarly,
f s n * ( A ( σ n s * ) ) = f s n ( A ( σ n s * ) )
for all n Z . It follows that
f + ( A ( σ n s ) = f - ( A ( σ n s ) , f + ( A ( σ n s * ) = f - ( A ( σ n s * )
for those n Z such that s n s n * (a nonempty set by hypothesis). This being the case, there are several conditions that prevent Equation (27) from occurring, thus guaranteeing the injectivity of Φ.
Assumption 2.
There is no A ( s ) A such that f + ( A ( s ) ) = f - ( A ( s ) ) or, equivalently, f + ( A ( s ) ) f - ( A ( s ) ) for all A ( s ) A .
This condition can be reworded as follows: f + and f - distinguish all pullback attractors of the switching dynamics.
Example 1
Let d = 1 , and f ± as in Equation (1). In this case, the pullback attractor A = { A ( s ) , s S } consists of singleton sets, i.e. A ( s ) = { χ ( s ) } , with χ ( s ) [ - 1 1 - θ - , 1 1 - θ + ] = A (see [1]). Let us check whether Assumption 2 holds:
f + ( χ ( s ) ) = f - ( χ ( s ) ) θ - θ + and χ ( s ) = 2 θ - - θ +
but 2 θ - - θ + [ - 1 1 - θ - , 1 1 - θ + ] because θ - + θ + < 2 . Since Assumption 2 holds, Φ is injective.
Example 2
Let A ± be the attractors of f ± , i.e., A + = A ( ( + 1 ) n Z ) , and A - = A ( ( - 1 ) n Z ) . If A + = A - , then Assumption 2 does not hold since f + ( A + ) = A + , and f - ( A - ) = A - . Think, for instance, of two linear dynamics f ± ( x ) = M ± x having only 0 as a fixed point. Therefore, in this case the topological entropy of the switched dynamics may be different from the topological entropy of the control sequence.
Now suppose that at least one of the maps f ± is invertible and define g : = f + - 1 f - if f + is invertible, or g : = f - - 1 f + otherwise; if both maps are invertible, then either choice is the inverse of the other one. Then
f + ( A ( s ) ) = f - ( A ( s ) ) A ( s ) = g ( A ( s ) ) .
Therefore, Assumption 2 can be particularized in the following way.
Assumption 2 for Invertible Mappings.
In case that one of the maps f ± is invertible, there is no g-invariant A ( s ) A .
Example 3
Suppose that M ± are two d × d -matrices, with M + being invertible, and f ± ( x ) = M ± x , so that g ( x ) = G x with G = M + - 1 M - . Thus, when applying Assumption 2 for invertible mappings, we are asking about the existence of invariant sets of known linear maps. Linear maps can be decomposed as actions of scalings (expansion/contraction, projection, reflection) and/or rotations. Therefore, the existence of G-invariant pullback attractors A ( s ) boils down to well-known geometrical facts.
  • The only invariant set of an (in general anisotropic) expansion/contraction is { 0 } .
  • If a linear map is a projection onto a subspace V, then any set contained in V is invariant.
  • If a linear map is a reflection, then any set symmetric under reflections (like a star-like shaped object) is invariant.
  • If a linear map is a rotation by an angle α, then any set symmetric under that rotation is invariant, in particular any ball centered at the origin of radius r 0 and dimension n or 0 n d - 1 if contained in the orthogonal complement of the rotation axis.
Since the origin is invariant in any case, it might happen that A ( s ) = { 0 } for some s S (for example, for ( + 1 ) n Z or ( - 1 ) n Z , since { 0 } A + , A - , where A ± are the attractors of f ± ). In this case, the version of Assumption 2 for invertible mappings would not hold.
From the preceding discussion we conclude the following result.
Lemma 3
Suppose that the maps f ± are such that Assumption 2 holds. Then the map Φ: S K R defined by Φ ( s ) : = ( A ( σ n s ) ) n Z is one-to-one.

5. Results Using Continuity

The main technical difference between the setting of [1] and the present general setting is that the map s A ( s ) is not necessarily continuous (Proposition 3). It is, however, worth following the same path as in [1] to see how far we can go and under which conditions. The first condition is obvious.
Assumption 3.
The set-valued mapping s A ( s ) is continuous in ( K R , dist H ) , i.e.,
dist H A ( s ) , A ( s * ) 0 a s dist S ( s , s * ) 0
Proposition 4 provides a sufficient condition for Assumption 3 to hold.
We are going to show in this section that if Assumptions 2 (Section 4) and 3 are fulfilled, then Φ Φ ( S ) is a homeomorphism.
Lemma 4
Under Assumption 3, the map Φ: S K R defined by Φ ( s ) = ( A ( σ n s ) ) n Z is uniformly continuous.
Proof.
We have to show that for any ε > 0 , there exists N N such that
dist K R ( A ( σ n s ) ) n Z , ( A ( σ n s * ) ) n Z < ε
whenever dist S ( s , s * ) < 2 - ( N + 1 ) (i.e., whenever s n = s n * for n N ). The uniformity follows then from the compactness of S .
Consider first the second term in the decomposition
dist K R ( ( A ( σ n s ) ) n Z , ( A ( σ n s * ) ) n Z = n N dist H A ( σ n s ) , A ( σ n s * ) 2 n + n N + 1 dist H A ( σ n s ) , A ( σ n s * ) 2 n
Since A ( σ n s ) , A ( σ n s * ) B ¯ R ( 0 ) for any n Z , we have
dist H A ( σ n s ) , A ( σ n s * ) 2 R ,
hence
n N + 1 dist H A ( σ n s ) , A ( σ n s * ) 2 n 4 n = N + 1 1 2 n - 1 = 4 2 N - 1 < ε 2
if N N I I : = log ( 8 / ε ) log 2 + 2
We use next Assumption 1 to show that the first term in Equation (29) with N = N I I can also be made smaller than ε / 2 . Indeed, for each n with n N I I there exist by Equation (28) N n N such that
dist S ( σ n s , σ n s * ) < 2 - ( N n + 1 ) dist H A ( σ n s ) , A ( σ n s * ) < ε 6
Set N I : = max { N n : n N I I } . Then
dist S ( s , s * ) < 2 - ( N I + N I I + 1 ) dist S ( σ n s , σ n s * ) < 2 - ( N n + 1 )
for n N I I since
( σ n s ) k = ( σ n s * ) k for k N n s k + n = s k + n * for k N n
where n N I I and k N n N I . Therefore if dist S ( s , s * ) < 2 - ( N I + N I I + 1 ) , then
n N I I dist H A ( σ n s ) , A ( σ n s * ) 2 n < ε 6 n N I I 1 2 n < ε 2
All in all, we conclude from Equation (29), Equation (31) and Equation (30) with N = N I I that if dist S ( s , s * ) < 2 - ( N I + N I I + 1 ) , then dist K R ( A ( σ n s ) ) n Z , ( A ( σ n s * ) ) n Z < ε . □
Remark.
The above proof exploits the special structure of the model and provides insight into it. The result, in fact, follows from the continuity of the mapping Φ and the compactness of the metric space ( S , dist S ) .
If the correspondence s Φ ( s ) : = ( A ( σ n s ) ) n Z is invertible, then given a sequence ( A n ) n Z Φ ( S ) there exists a unique control sequence s = ( s n ) n Z such that A n + 1 = f s n ( A n ) for every n. Therefore, s n can be determined from the knowledge of A n and A n + 1 .
Lemma 5
Suppose that Assumption 2 holds so that the mapping Φ is one-to-one. Then, Φ - 1 : Φ ( S ) S is continuous.
Proof.
Given N > 0 , we have to prove that there exists δ > 0 such that
dist K R ( A n ) n Z , ( A n * ) n Z < δ dist S s , s * < 2 - N
where ( A n ) n Z , ( A n * ) n Z Φ ( S ) , and s = Φ - 1 ( ( A n ) n Z ) , s * = Φ - 1 ( ( A n * ) n Z ) .
From
dist K R ( A n ) n Z , ( A n * ) n Z = n N dist H ( A n , A n * ) 2 n + n > N dist H ( A n , A n * ) 2 n < δ
it follows
dist H ( A n , A n * ) < 2 N δ for n N
Consider now two neighboring sequences in Φ ( S ) ,
. . . A - N , A - N + 1 , . . . , A 0 , . . . , A N - 1 , A N , . . . , . . . A - N * , A - N + 1 * , . . . , A 0 * , . . . , A N - 1 * , A N * , . . . ,
of distance δ apart. From these sequences we determine the corresponding switching sequence under Φ - 1 :
. . . s - N , s - N + 1 , . . . , s 0 , . . . , s N - 1 , . . . , . . . s - N * , s - N + 1 * , . . . , s 0 * , . . . , s N - 1 * , . . . .
Note that A n = A ( σ n s ) , where s = ( s n ) n Z , and similarly A n * = A ( σ n s * ) , where s * = ( s n * ) n Z .
Now suppose that, contrarily to what we need to prove, there is an n, - N n N - 1 , such that s n * = - s n . For brevity assume s n = + 1 , i.e.
A n + 1 = f + ( A n ) , and A n + 1 * = f - ( A n * ) .
Then, by the triangle inequality,
dist H ( f + ( A n * ) , f - ( A n * ) ) dist H ( f + ( A n * ) , f + ( A n ) ) + dist H ( f + ( A n ) , f - ( A n * ) ) .
We claim that the right hand side of Equation (34) can be made arbitrarily small by choosing δ in Equation (32) sufficiently small.
First of all, by Proposition 2 and Equation (32)
dist H ( f + ( A n * ) , f + ( A n ) ) < ω ( dist H ( A n * , A n ) )
where ω ( · ) is continuous and ω ( 0 ) = 0 . As for the second term on the right hand side of Equation (34), use Equation (33) and Equation (32) to derive
dist H ( f + ( A n ) , f - ( A n * ) ) = dist H ( A n + 1 , A n + 1 * ) < 2 N δ
Therefore,
dist H ( f + ( A n * ) , f - ( A n * ) ) < ω ( dist H ( A n * , A n ) ) + 2 N δ
where
ω ( dist H ( A n * , A n ) ) 0 when δ 0
This proves our claim.
It follows from Equation (35) and Equation (36) that f + ( A n * ) = f - ( A n * ) , in contradiction to Assumption 2. Hence we conclude that s n * = s n for - N n N - 1 , i.e., that dist S s , s * < 2 - N , where s = Φ - 1 ( ( A n ) n Z ) , and s * = Φ - 1 ( ( A n * ) n Z ) .
Notice that we did not assume Φ to be continuous in Lemma 5. With Lemmas 4 and 5 we derive a different version of Theorem 2.
Theorem 4
Suppose that Assumptions 2 and 3 hold. Then the mapping Φ: S Φ ( S ) defined as s ( A ( σ n s ) ) n Z is a homeomorphism, hence h ( σ ) = h ( Σ Φ ( S ) ) , i.e., h c o n t r o l = h m a c r o h m i c r o .

6. Numerical Simulations

In this section we illustrate the previous theoretical results with a few numerical simulations. As in [1] we resort to two-state Markov processes with transition probability matrices
P = p 1 - p q 1 - q
to generate the control sequences. We recall that if the transition probability matrix is irreducible, then the process is ergodic (and even there is a unique invariant measure whose metric entropy coincides with the topological entropy of the corresponding topological Markov chain).

6.1. Case 1

Consider f ± : R 2 R 2 with
f + ( x , y ) = ( 0 . 25 x , 0 . 75 y )
and the Hénon map
f - ( x , y ) = ( 1 - 1 . 4 x 2 + 0 . 3 y , x )
A common bounded, positively invariant absorbing set for Equation (38) and Equation (39) is the quadrilateral with vertices [16] (p. 614).
P 1 = ( - 1 . 33 , 0 . 42 ) , P 2 = ( 1 . 32 , 0 . 133 ) , P 3 = ( 1 . 245 , - 0 . 14 ) , P 4 = ( - 1 . 06 , - 0 . 5 )
Therefore, the switching system defined by Equation (38) and Equation (39) has a pullback attractor, which is also a forward attractor (Remark 1). We take advantage of this fact in the estimation of the topological entropy by taking a sample of long forward trajectories and deleting the first 100 points.
The topological entropy of these two maps in bits per iteration (i.e., taking logarithms to base 2) are
h ( f + ) = 0 , h ( f - ) = 0 . 67072 ± 0 . 00004
( h ( f - ) taken from [17]). Table 1 summarizes the numerical results. h c o n t r o l h ( σ ) ( = log λ , where λ is the largest positive eigenvalue of the admissibility matrix [1,15] denotes the topological entropy of the switching sequence. As above, h m i c r o denotes the topological entropy of the switched dynamics. Conditioned on the assumptions of Theorem 3 or 4, the inequality h c o n t r o l h m i c r o holds. Figure 1 shows the extrapolation technique we used based on entropy rates of order 3 L 7 (see [1] for more details). A total of 40 orbits were generated for each choice of the probabilities ( p , q ) = ( 0 . 7 , 0 . 2 ) , ( 0 . 25 , 1 ) , and ( 0 , 1 ) , with random initial conditions within the trapping region (40). Note that h ( f + ) , h ( f - ) < h c o n t r o l for ( p , q ) = ( 0 . 7 , 0 . 2 ) , ( 0 . 25 , 1 ) .
Table 1. Numerical results of Case 1 ( h c o n t r o l h m i c r o ).
Table 1. Numerical results of Case 1 ( h c o n t r o l h m i c r o ).
( p , q ) h control (bit/iter) h micro ± s . d . (bit/iter)
( 0 . 7 , 0 . 2 ) 1 1 . 237 ± 0 . 002
( 0 . 25 , 1 ) 0 . 694 0 . 702 ± 0 . 013
( 0 , 1 ) 0 0 . 000 ± 0 . 000
Figure 1. This plot illustrates the extrapolation technique used to estimate the values of h m i c r o in Table 1.
Figure 1. This plot illustrates the extrapolation technique used to estimate the values of h m i c r o in Table 1.
Entropy 15 02363 g001
It turns out that, except for the periodic switching ( p , q ) = ( 0 , 1 ) , the complexity of the switched dynamics is higher than the complexity of the control switching sequence, and also higher than the entropies of the constituent maps, see Equation (41). This provides an instance of Parrondo’s paradox with regard to the complexity as measured by the topological entropy.

6.2. Case 2

Consider the constituent affine maps
f ± ( x ) = M ± x + b ±
where x = ( x 1 , x 2 ) R 2 ,
M + = 1 4 5 - 4 2 - 1 M - = 1 4 5 - 6 3 - 4
and
b + = 3 2 b - = - 3 - 2
A special feature of the matrices M ± is that they are simultaneous diagonalizable. Indeed, in the common eigenvector base
v 1 = 2 1 v 2 = 1 1
the constituent maps are transformed to (we keep the notation x = ( x 1 , x 2 ) for the coordinates in the new base)
f + ( x 1 , x 2 ) = 1 4 3 0 0 1 x 1 x 2 + 1 1 = 3 4 x 1 + 1 1 4 x 2 + 1 , f - ( x 1 , x 2 ) = 1 4 2 0 0 - 1 x 1 x 2 + - 1 - 1 = 1 2 x 1 - 1 - 1 4 x 2 - 1 ,
Therefore, the two-dimensional switched dynamics (42) decomposes in the base (43) into two uncoupled, one-dimensional switched dynamics whose constituent maps are of the type (1) studied in [1]. Here the pullback attractors are singleton sets (hence the macroscopic and microscopic complexities coincide), Assumptions 2 and 3 are fulfilled, and h m i c r o = h m a c r o = h c o n t r o l , the latter equality holding by Theorem 4. Furthermore, h ( f + ) = h ( f - ) = 0 , thus h ( f + ) , h ( f - ) < h c o n t r o l for ( p , q ) = ( 0 . 7 , 0 . 2 ) , ( 0 . 25 , 1 ) , as in Case 1. Table 2 shows numerical evidence of this statement. Figure 2 shows the extrapolation technique based on entropy rates of orders 3 L 7 .
Table 2. Numerical results of Case 2 ( h c o n t r o l = h m i c r o ).
Table 2. Numerical results of Case 2 ( h c o n t r o l = h m i c r o ).
( p , q ) h control (bit/iter) h micro ± s . d . (bit/iter)
( 0 . 7 , 0 . 2 ) 1 1 . 002 ± 0 . 003
( 0 . 25 , 1 ) 0 . 694 0 . 687 ± 0 . 007
( 0 , 1 ) 0 0 . 000 ± 0 . 000
Figure 2. This plot illustrates the extrapolation technique used to estimate the values of h m i c r o in Table 2.
Figure 2. This plot illustrates the extrapolation technique used to estimate the values of h m i c r o in Table 2.
Entropy 15 02363 g002

7. Conclusions

Consider the time-switched dynamics associated with two maps f ± : R d R d having global attractors A ± . In this paper we studied the relation between the topological entropy of the control sequence generator and the topological entropy of the switched dynamics. Compared with the previous paper [1], we found some new technical obstacles.
First of all, we had to circumvent the shortcoming that the map s A ( s ) is not necessarily continuous in ( K R , dist H ) ; here s is a control sequence and A ( s ) is the corresponding component set of the pullback attractor, A = { A ( s ) , s S } . To this end we assumed in Section 3 (“Results using measurability”), (i) the measurability of certain unions of component sets (Assumption 1); and (ii) the injectivity of Φ : s ( A ( σ n s ) ) n Z . It follows then, Theorem 2, that Φ is a topological conjugacy between the shift on the control sequences, ( { + 1 , - 1 } Z , σ ) , and the shift on the trajectories of the switched dynamics up to the precision set by the component sets, ( K R , Σ ) . Being the entropy of coarse-grained trajectories, h ( Σ ) is a lower bound of h m i c r o , the topological entropy of the switched dynamics. We call h ( Σ ) the macroscopic entropy and write h m a c r o . Likewise, h ( σ ) is the entropy of the control sequence generator, thus we write h c o n t r o l . In sum, h c o n t r o l = h m a c r o h m i c r o [Equation (25) and Equation (26)] under the hypotheses (i) and (ii).
In Section 4 we surveyed sufficient conditions for the injectivity of Φ postulated in Section 3. We formulated our conclusions as Assumption 2, both in a general version, and in a special version when one of the maps f ± is invertible.
In Section 5 (“Results using continuity”) we hypothesized from the outset that the mapping s A ( s ) is continuous (Assumption 3) and paralleled the approach of [1]. In Theorem 4 we proved that Φ is a homeomorphism (or topological conjugacy), Assumptions 2 and 3 granted. Once again we concluded that h c o n t r o l = h m a c r o h m i c r o .
The main results of this paper, Theorems 3 and 4, are in line with what we call the complexity version of Parrondo’s paradox, namely, that the complexity of a switched dynamics may be higher than the complexity of its constituent maps. Moreover, they provide sufficient conditions as well, to wit: h ( f + ) , h ( f - ) < h c o n t r o l .
Finally, both possibilities h c o n t r o l < h m i c r o (Case 1) and h c o n t r o l = h m i c r o (Case 2) have been numerically illustrated in Section 6.

Acknowledgements

J.M.A. and A.G. were financially supported by the Spanish Ministerio de Economía y Competitividad, project MTM2012-31698. P.E.K. was financially supported by the Spanish Ministerio de Ciencia e Innovación, project MTM2011-22411, and Conserjería de Innovación, Ciencia y Empresa (Junta de Andalucía) under the Ayuda 2009/FQM314 and the Proyecto de Excelencia P07-FQM-02468.

Conflict of Interest

The authors declare no conflict of interest.

Appendix

Recall that the Hausdorff semi-distance and Hausdorff distance between two nonempty subsets X and Y of a metric space ( M , dist ) are denoted by ρ ( X , Y ) and dist H ( X , Y ) , respectively, while
diam X = sup x , x X dist ( x , x )
is called the diameter of a set X.
Let
B ε ( x ) = z M : dist ( z , x ) < ε , B ¯ ε ( x ) = z M : dist ( z , x ) ε
be the open and closed balls in ( M , dist ) with center at x and radius ε > 0 , respectively. The set of all points within a distance ε of X, i.e.,
E ε ( X ) = z X B ¯ ε ( z )
is called the ε-inflation (or ε-expansion) of X. It follows trivially from this definition that
X E ε ( Y ) and Y E ε ( X ) dist H ( X , Y ) ε
The converse of Equation (45) is in general false, except when X and Y are compact.
Proposition A1.
Let X and Y be compact subsets of ( M , dist ) . Then
dist H ( X , Y ) = ε X E ε ( Y ) and Y E ε ( X )
Proof.
If dist H ( X , Y ) = ε , then ρ ( X , Y ) ε and ρ ( Y , X ) ε . This amounts to
max x X dist ( x , Y ) ε , where dist ( x , Y ) : = min y Y dist ( x , y )
and
max y Y dist ( y , X ) ε , where dist ( y , X ) : = min x X dist ( x , y )
respectively. Therefore, see Equation (46), for every x X , there exists y Y such that x B ¯ ε ( y ) , and, see Equation (47), for every y Y , there exists x X such that y B ¯ ε ( x ) . □
It follows from Proposition A1 that
dist H ( X , Y ) ε X E ε ( Y ) and Y E ε ( X )
when X and Y are compact.
Among the different properties of the Hausdorff distance, we mention here only two of them. Let K ( M ) be the set of all nonempty compact subsets of M.
  • If ( M , dist ) is complete, then ( K ( M ) , dist H ) is also a complete.
  • If ( M , dist ) is compact, then ( K ( M ) , dist H ) is also a compact.
In Section 2.2 we consider the metric space M = R d endowed with the Euclidean distance · , as well as the complete, compact metric space ( K R , dist H ) , where K R = K ( R d ) B ¯ R ( 0 ) .
Lemma A2.
Let ( M 1 , dist 1 ) be a compact metric space and ( M 2 , dist 2 ) a complete metric space. If f : M 1 M 2 is a continuous map, then there exists a real-valued function ω : [ 0 , ) [ 0 , ) with ω ( 0 ) = 0 such that, for any X , Y K ( M 1 ) ,
dist 2 , H ( f ( X ) , f ( Y ) ) ω ( dist 1 , H ( X , Y ) )
where dist 1 , H (resp. dist 2 , H ) is the Hausdorff distance in K ( M 1 ) (resp. K ( M 2 ) ).
The function ω here is called a modulus of continuity. Moduli of continuity are used to express in a convenient way both the continuity at a point and the uniform continuity of maps between metric spaces as, for instance, in Equation (49).
Proof of Lemma A2.
Let X and Y be compact subsets of ( M 1 , dist 1 , H ) such that dist 1 , H ( A , B ) δ . By Equation (48),
X i = 0 m B ¯ δ ( y i ) and Y j = 0 n B ¯ δ ( x j )
where x j X and y i Y . Thus,
f ( X ) i = 0 m f ( B ¯ δ ( y i ) ) and f ( Y ) j = 0 n f ( B ¯ δ ( x j ) )
In turn, f ( B ¯ δ ( y i ) ) and f ( B ¯ δ ( x j ) ) are compact in M 2 because f : M 1 M 2 is continuous. This being the case,
f ( B ¯ δ ( y i ) ) B ¯ r i ( f ( y i ) ) and f ( B ¯ δ ( x j ) ) B ¯ s j ( f ( x j ) )
where [see Equation (44)]
r i diam f ( B ¯ δ ( y i ) ) and s j diam f ( B ¯ δ ( x j ) )
for 1 i m and 1 j n . In sum,
f ( X ) i = 0 m B ¯ r i ( f ( y i ) ) and f ( Y ) j = 0 n B ¯ s j ( f ( x j ) )
Let ω f : [ 0 , + ] [ 0 , + ] be the modulus of continuity of the uniformly continuous map f of M 1 into M 2 , i.e.,
dist 2 ( f ( z ) , f ( z ) ) < ω f ( dist 1 ( z , z ) )
where
lim t 0 ω f ( t ) = ω f ( 0 ) = 0
Without loss of generality, ω f may be assumed to be continuous and even differentiable. Then
r i diam f ( B ¯ δ ( y i ) ) = max { dist 2 ( f ( z ) , f ( z ) ) : z , z B ¯ δ ( y i ) } max { ω f ( dist 1 ( z , z ) ) : z , z B ¯ δ ( y i ) }
for 1 i m , and, analogously,
s j max { ω f ( dist 1 ( z , z ) ) : z , z B ¯ δ ( x j ) }
for 1 j n . From Equation (51) we conclude that r i and s j can be made arbitrarily small by choosing δ sufficiently small.
All in all, given ε > 0 and X , Y M 1 compact such that dist H ( A , B ) < δ , we can choose δ sufficiently small so that r i , s j ε for 1 i m , and 1 j n . From Equation (50) we obtain
f ( X ) i = 0 m B ¯ ε ( f ( y i ) ) E ε ( f ( Y ) ) and f ( Y ) j = 0 n B ¯ ε ( f ( x j ) ) E ε ( f ( X ) )
hence dist 2 , H ( f ( X ) , f ( Y ) ) ε by Equation (45).
This, in fact, proves the uniform continuity (since K ( M 1 ) is compact) of the set-valued map K ( M 1 ) A f ( A ) K ( M 2 ) . Equation (49) is just an equivalent formulation of this fact. □

References

  1. Amigó, J.M.; Kloeden, P.E.; Giménez, A. Switching systems and entropy. J. Differ. Equ. Appl. 2013. [Google Scholar] [CrossRef]
  2. Kloeden, P.E.; Rasmussen, M. Nonautonomous Dynamical Systems; American Mathematical Society: Providence, RI, USA, 2011. [Google Scholar]
  3. Liberzon, D. Switching in Systems and Control; Birkhäuser: Boston, MA, USA, 2003. [Google Scholar]
  4. Shorten, R.; Wirth, F.; Mason, O.; Wulff, K.; King, C. Stability criteria for switched and hybrid systems. SIAM Rev. 2007, 49, 545–592. [Google Scholar] [CrossRef]
  5. Almeida, J.; Peralta-Salas, D.; Romera, M. Can two chaotic systems give rise to order? Phys. D: Nonlinear Phenom. 2005, 200, 124–132. [Google Scholar] [CrossRef]
  6. Boyarsky, A.; Góra, P.; Islam, M.S. Randomly chosen chaotic maps can give rise to nearly ordered behavior. Phys. D: Nonlinear Phenom. 2005, 210, 284–294. [Google Scholar] [CrossRef]
  7. Cánovas, J.S.; Linero, A.; Peralta-Salas, D. Dynamic Parrondo’s paradox. Phys. D: Nonlinear Phenom. 2006, 218, 177–184. [Google Scholar] [CrossRef]
  8. Harner, G.P.; Abbott, D. Parrondo’s paradox. Stat. Sci. 1999, 14, 206–213. [Google Scholar]
  9. Kloeden, P.E. Pullback attractors in nonautonomous difference equations. J. Differ. Equ. Appl. 2000, 6, 33–52. [Google Scholar] [CrossRef]
  10. Cheban, D.N.; Kloeden, P.E.; Schmalfuß, B. The relationship between pullback, forward and global attractors of nonautonomous dynamical systems. Nonlinear Dyn. Syst. Theory 2002, 2, 125–144. [Google Scholar]
  11. Kloeden, P.E. Pullback attractors of nonautonomous semidynamical systems. Stoch. Dyn. 2003, 3, 101–112. [Google Scholar] [CrossRef]
  12. Kloeden, P.E.; Siegmund, S. Bifurcations and continuous transitions of attractors in autonomous and nonautonomous systems. Int. J. Bifurcat. Chaos Appl. Sci. Eng. 2005, 15, 743–762. [Google Scholar] [CrossRef]
  13. Kloeden, P.E. Nonautonomous attractors of switching systems. Dyn. Syst. 2006, 21, 209–230. [Google Scholar] [CrossRef]
  14. Amigó, J.M. Permutation Complexity in Dynamical Systems; Springer: Berlin, Germany, 2010. [Google Scholar]
  15. Walters, P. An Introduction to Ergodic Theory; Springer: Berlin, Germany, 2000. [Google Scholar]
  16. Peitgen, H.O.; Jürgens, H.; Saupe, D. Chaos and Fractals: New Frontiers of Science, 2nd ed.; Springer: Berlin, Germany, 2004. [Google Scholar]
  17. Hirata, Y.; Mees, A.I. Estimating topological entropy via a symbolic data compression technique. Phys. Rev. E 2003, 67, 026205. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Amigó, J.M.; Kloeden, P.E.; Giménez, Á. Entropy Increase in Switching Systems. Entropy 2013, 15, 2363-2383. https://doi.org/10.3390/e15062363

AMA Style

Amigó JM, Kloeden PE, Giménez Á. Entropy Increase in Switching Systems. Entropy. 2013; 15(6):2363-2383. https://doi.org/10.3390/e15062363

Chicago/Turabian Style

Amigó, José M., Peter E. Kloeden, and Ángel Giménez. 2013. "Entropy Increase in Switching Systems" Entropy 15, no. 6: 2363-2383. https://doi.org/10.3390/e15062363

APA Style

Amigó, J. M., Kloeden, P. E., & Giménez, Á. (2013). Entropy Increase in Switching Systems. Entropy, 15(6), 2363-2383. https://doi.org/10.3390/e15062363

Article Metrics

Back to TopTop