2021
263
17012020
Roberta Hansen rhansen@fi.uba.ar
Universidad de Buenos Aires, Argentina
Graciela A González ggonzal@fi.uba.ar
Consejo Nacional de Investigaciones Científicas y Técnicas CONICET, Argentina
Abstract:
Based on existing feedback control methods such as OGY and Pyragas, alternative new schemes are proposed for stabilization of unstable periodic orbits of chaotic and hypercha- otic dynamical systems by suitable modulation of a control parameter. Their performances are improved with respect to: (i) robustness, (ii) rate of convergences, (iii) reduction of waiting time, (iv) reduction of noise sensitivity. These features are analytically investigated, the achievements are rigorously proved and supported by numerical simulations. The proposed methods result successful for stabilizing unstable periodic orbits in some classical discrete maps like 1-D logistic and standard 2-D Hénon, but also in the hyperchaotic generalized .-D Hénon-like maps.
Keywords: chaos control, proportional feedback control, delayed feedback control, parameter modulation, control of hyperchaos.
Article
Feedback control modulation for controlling chaotic maps
Roberta Hansen rhansen@fi.uba.ar
Graciela A González ggonzal@fi.uba.ar
Recepción: 17 Enero 2020
Revisado: 11 Diciembre 2020
Publicación: 01 Mayo 2021
Chaotic behavior is a very interesting nonlinear phenomena, but in many situations, it is desirable to be avoided, for example, when it restricts the operating range of electronic or mechanics devices. Moreover, this goal should be achieved with the only help of tiny perturbations properly chosen [2].
The consagrated idea of Ott, Grebogi, and Yorke [24] consists on turning the presence of chaos into an advantage. Indeed, the system may be stabilized on a particular unstable periodic orbit (UPO) embedded in a strange attractor by applying a small time-dependent feedback perturbation to some accessible parameter or variable system. The periodic orbit is preserved, but its stability is modified keeping the trajectory to stay close to the UPO. This control strategy is known as the OGY method.
A simple proportional feedback (SPF) control method basically consists on a per- turbation proportional to the difference between the current system value and an unstable fixed point (or one of the UPO’s component), being the OGY method [24] a particular case of it. It is well known that it can control chaos for some 1-D maps [2]. Its usefulness arises when controlling chaos in highly dissipative systems; as there are cases in which it loses validity because Poincaré sections changes in each iteration, a recursive proportional feedback (RPF) has been proposed [36]. Both SPF and RPF methods require the exact knowledge of the UPO to be stabilized and the linearized dynamics about it. This is not always at one’s disposal in real-world implementations. Moreover, the OGY method results very sensitive to nonlinearities and fluctuations of external noise, mainly for large orbit periods. Nevertheless, the OGY method holds good, and extended versions of it are still studied (for example, in [22], a two-controlling parameter extension is applied on a coupled 2-D chaotic maps). With the aim to overcome the limitations of the OGY method, Pyragas introduced a self-controlling delayed feedback [27], which, as the OGY, does not modify the original UPO, but it does not depend explicitly on it. Its discrete-time version results in a perturbation, which is proportional to the difference between the cur- rent system value and a previous one. This delayed feedback control (DFC) successfully controls chaotic behavior in a variety of experiments [6] (and its references). However, the stabilization capability of the DFC may be weaker than the OGY’s. The domain of system parameters for which stabilization can be achieved via DFC is limited, namely, the method fails for highly UPO’s [32]. The Pyragas method is renamed in [32] as “time-delay autosynchronization” (TDAS). Its extended versions (referred as ETDAS and EDFC) uses information from many previous states improving on stabilization objectives achievement. Recently published works aim to the efficient choice of the control parameters for the EDFC implementation. In [1], the parameter design without accurate information about the system is dealt with by employing a data-driven pole placement method. Moreover, in [4], a fixed point stabilization by means of EDFC methodology is warranted under a condition on the set of the eigenvalues of the Jacobian matrix of the map at the fixed point by rigorous mathematical approach although the extension of its results to a .- period UPO stabilization seems to remain an open problem. On the other hand, any UPO with an odd number of real Floquet multipliers greater than one (or with an odd number of real positive Floquet exponents) can never be stabilized by TDAS/ETDAS method. This restriction was pointed out in [34], and it is known as the odd number limitation (ONL). Later, ONL was also stated for stabilization of UPO’s through ETDAS [33]. Several alternative methods, as the oscillating approach [20], the adaptive version [26], the delayed feedback with periodic control gain [16], the delayed feedback control methods based on predictive state [21, 35], the learning control of time-delay chaotic systems [11], the act-and-wait concept [8], and the observed-based-delayed feedback control [12], have been proposed to overcome this drawback. More references embracing the discrete and continuous-time cases may be found in [17].
Assume a dynamical system given by a continuously differentiable function . :
(1)where r is a scalar parameter, and that for
., the system develops chaotic behavior having an infinite number of UPOs embedded in a strange attractor, A, within its basin of attraction, B The dynamics in A is ergodic, meaning that for almost all initial condition in , a system trajectory visits any small neighborhood of every point of the UPO’s.
A widely used technique for controlling chaotic behavior is modulation: by judi- ciously varying control parameters any system trajectory resulting from a randomly cho- sen initial condition,
, is driven to a given .-period UPO,
Suppose that r can be finely tuned in a small-range around r0, namely,
Then the objective is to stabilize system (1) at . by feedback control mod- ulation. Namely, the (control) parameter is affected by a control 
(2)under the requirement that the adjusted parameter remains within a range for which the system is chaotic in the absence of perturbations.
Published works on controlling chaos by modulation are mostly based on experimen- tal or numerical arguments, while analytical approaches usually concentrate on a fixed point stabilization problem [31]. DFC and SPF controls applied to 1-D maps are con- fronted in [6] by means of linear stability analysis and by implementation in an analog electric circuit. On the other hand, analytical approaches consider only additive control, and features related to control of chaos objectives are disregarded [16, 18, 27].
Both OGY and Pyragas methods have been extensively studied and analyzed in the literature on the subject like [3] and also revisited in [29, 30] among others, while several modifications on them have been introduced throughly more than twenty years ago. The present work relies on modifications of these control methods by modulation techniques improving some aspects on their performances related to robustness, rate of convergences, reduction in the waiting time, reduction in noise sensitivity, among others. Moreover, all of these features are rigorously proved. The approach to stabilize any UPO of the chaotic map (1) by modulation of the parameter . is investigated both analytically and numerically. The possibility of relaxing some requirements of the OGY and Pyragas methods is explored, and alternative types of SPF and DFC methods are built. Different from the original OGY method, the SPF method does not require exact knowledge on the linearization data of the UPO to stabilize. Based on these ideas, modifications to the original Pyragas method [27] are introduced even in the extended version presented in [32], and alternative types of DFC and EDFC methods are also proposed. In all the cases, the convergence properties of the control strategies are proven taking into account control bounds as well as control performance issues of the considered schemes.
The article is organized as follows. Throughout the work, the problem is stated for
the general case of Eq. (1). The achievements of the proposed methods when applied to it are rigorously analyzed. By simplicity they are developed for 1-D maps at first. The generalized versions of the SPF and DFC methods, which involve a kind of switching controller, are introduced in Sections 2 and 3, respectively, and illustrated with the cele- brated logistic map. In Section 4, the issue of their extension to the .-dimensional case is established, and some meeting points with previously developed proposals are considered. Numerical simulations using the 2-D Hénon map illustrate the goals. In Section 5, an ad- hoc reformulation is designed to control the hyperchaotic .-D generalized Hénon-like maps. Finally, Section 6 contains the concluding remarks and lines of future research
This method outlined in [24] is briefly reviewed here. Assume that at time ., the trajectory falls close to the component p.. The linearized dynamics of (1) about pi ,and r0 is
(3)To force the system towards p, it is settled 
so from (3) it follows
(4)Equation (4) holds only when
hence, the required parameter perturbation,
is small, and the maximum parameter perturbation,
is proportional to . (with factor α.). When the trajectory is outside the .-neighborhood of
, the perturbation is not applied, and the system evolves at its nominal chaotic parameter 
Note that by ergodicity the control is eventually activated. The lapse while the control is off is known as waiting time. Then, for given ε, δ > 0, the control algorithm is
(5)Note that . is preserved by a parameter perturbation of SPF type, i.e.,
.
The control in Eq. (5) is turned on only at the end of each oscillation becoming very sensitive to nonlinearities and to fluctuation of external noises, mainly for large values of the orbit’s period, m. Namely, for every ., different control gain and waiting time come out, yielding to different performances for each
.
In order to improve these features, a sort of “switching control” is first proposed, which works as the OGY but applying the perturbation
for all 1 ≤ i ≤ m, i.e., each time the trajectory is close to any pi (a similar modification appears in [14] as part of a numerical implementation). As an immediate consequence, a net reduction on waiting time is obtained. The control algorithm is
(6)provided that
As a second consequence, this strategy
displays a notably better performance in presence of external noise. This fact is illustrated when applied to the well-known logistic map defined by
(7)
where 
and . is a control parameter. For
the map presents one chaotic attractor with B = [0,1] (except for small windows of periodicity). For fixed 
in that range, the parameter modulation is stated through the following controlled system:
(8)For the map (7),
The condition
makes the dynamics remain globally bounded, so, for
Figure 1 illustrates the comparison of applying the controls of Eqs. (5) and (6) about a 4-UPO of the map (7) for .. = 3.8. Even in the presence of noise, the shortening of the waiting time and the reduction in the control bound may be appreciated. Note that, once activated, the control stays on.
Here SPF modulation is selected from a set of control laws built by replacing each fixed gain α. with a coefficient β. adequately chosen:
(9)provided that
.This dynamics also preserves the .-UPO.
Once the control (9) is activated, Eq. (2) becomes
where
(10)The linear stability criterion for p. to be an asymptotically stable (a.s.) point, states the condition for β.. By applying chain rule, if
there exists a range
of βi values such that
(11)Under the stated conditions, convergence is formally proven.
Let the controlled system (2)and (9)with β. ∈
,
, for which Eq. (11)is valid and
Given δ > 0, there exists
such that for all
and for almost every initial condition x0 ∈ B, it verifies
for all k, and
converges to the m-periodic orbit
.
Proof.
. For all i, there exists σ > 0 such that |fβ′ . (p.)| < σ < 1, and for every ., there exists 0
such that
. Fixing
for almost every 
, the ergodicity of the uncontrolled system guarantees the existence of
for some i, then
Then by recursion it yields to
and the thesis follows.
The existence of a range
states the robustness of the method of Eq. (9),
which includes Eq. (6) (note that for
. By means of (9) the objective is fulfilled with smaller values for the control gain, or else, for the same control effort, 
a greater 
is allowed improving the waiting time to active the control. These facts are illustrated for the logistic map (7) for which Eq. (10) becomes
with
For r0. = 3.8, β.must verify
for all . to ensure the desired bound on the control effort and a globally bounded dynamics. As an example, the performances of controls of Eqs. (6) and (9) are compared when applied to stabilize the unstable fixed point
In Fig. 2(a), it is appreciated both: the neat reduction of control effort by changing the control gain 
by
, and an increase in the convergence time once the control is turned on. Figure 2(b) shows, for the similar control effort, the reduction in the waiting time explained above.
Verification of
-nearness becomes superfluous once K0 is detected.
The choosing of β. implies (but it is not equivalent to!):
(12)
Hence, a generalization of (9), which also includes (5), results in
(13)where
for some
for 
Note that the existence of K0 is warranted by system ergodicity.
Let the controlled system (2)and (13)with
m, for which Eq. (12)is valid and
there exists
. such that for all
, and for almost every initial condition x0 ∈ B, it verifies
converges to the m-periodic orbit
.
Proof. Let
such that
equivalent to
Therefore, there exists
such that
for
. Proceeding in the same way for
. values arise. Let
Without loss of generality, it is assumed that the condition
is verified, for the first time, by i = 1. Then
for
and so following:
By the same arguments it is obtained that 
and 
for
. Hence, control bounds and convergence to the periodic orbit follow.
Condition (11) relies on the existence of the range 
for each i.
Looking for a single common 
value, as referred in [3, Sect. 5.3.2], assumes a nonempty intersection of these intervals. However, this may not be the case in many examples as it is easily seen for the same 4-period UPO of Fig. 1 for which the ranges for .. and .. are (−11.0716, −2.325) and (8.2897, 20.9639), respectively.
The problem of finding a control depending just on one gain control for stabilizing an m-UPO is left unsolvable in [25]. In the terms of this work, it requires to find a value . such that Eq. (12) is satisfied with
for all . involving a nonlinear inequality whose solution is not warranted in any case.
For stabilizing the logistic map (7) to its unstable fixed point, p = 1 − 1/r0, additive forc- ing in the form of one time-delay linear perturbation,
was early introduced in [27]. Here DFC modulation is stated as follows:
(14)Control (14) vanishes when system (2) state attains ., then the fixed point is preserved. For the logistic map, the controlled system (8) and (14) yields to the two-dimensional dynamical system
which has P = [p, p] as a fixed point. The Jacobian matrix at P is a (2 x 2)-companion matrix
where
and 
The necessary and sufficient conditions for P to be a.s. [13] yields to a range
Fixing a . in this range and choosing an adequate . to assure the control effort to be bounded by ., the convergence of the trajectories to p is obtained. The resulting control performance of applying (14) maybe comparable to (or even better than) the one obtained by applying (9) if adequatecoefficients are chosen as shown in Fig. 3.
The idea of “switching control” that depends on .-different control gains adequately chosen, introduced in Section 2 for improving SPF modulation, is here applied to improve DFC modulation. The proposal of an extension of (14) to a .-UPO, ., of (1) is to set a “switching” control, as follows:
(15)
System (2) and (15) yields to a (m + 1)-dimensional system (with . switches) with 
as its state variables:
(16)Note that
is m-UPO of the free system (γ. = 0 for all i), and it is preserved when (15) is applied. The Jacobian matrix of system (16) at Pi is a
companion matrix given by
where
.
Unfortunately, in most cases for m > 1, there is no γi that warranties all of Ji eigenvalues to be of modulus less than one, and so, achieving UPO stabilization by this way. Even it is not possible for the 2-UPO of the logistic map.
The conditions for the orbit to be a.s. involve the product of the matrices, 
It should be found γi making this product have all eigenvalues of modulus less than one.
This is the equivalent of condition in Eq. (12) for the control (13). Moreover, it reduces to the (only) Jacobian matrix of the system for the fixed point stabilization (Eq. (14)).
As each .. is a companion matrix, the characteristic polynomial, 
of the product .. is obtained by using tools from [10], and it results
(17)Therefore, the control is proposed as
(18)if
fo. all
, being
for some
Existence of .. is a consequence of ergodicity. Its actual computation requires certain information on UPO’s location, but this does not violate the main feature of DFC method, i.e., exact knowledge of the UPO is not necessary
Introducing
and
the controlled system (2) and (18) is equivalent to
which has
UPO. Note that
and
. Next, following the same steps that in the proof of Proposition 2, the success of the proposed scheme is proved.
Let the controlled system (2)and (18)with γj, 1 ≤ j ≤ m, such that all the roots of
are of modulus less than one, and
Given δ > 0, there exists
such that for all
and for almost every initial condition x0 ∈ B, it verifies
and that .
converges to the m-UPO
Proof. The roots of 
are of modulus less than one, so
with the matrix norm induced by the Euclidean one. Let
And 0 < σ < 1 such that
which is the same as
Then there exists
such that
for 
Proceeding analogously for
it is obtained
and
. Without loss of generality, it is assumed that condition
is verified, for the first time, by i = 1. Then
where
. lies in the segment joining P1 and
By the same arguments it is obtained
(19)Hence, the convergence of
to the m-UPO
follows, and so the convergence of
to the m-UPO
. From Eq. (19) it also results
and so, the control bounds are obtained.
The successfulness of the controlled system (2) and (18) hinges on finding a range ofγi’s for orbit stability, i.e., for the roots of
in Eq. (17) to be within the unit circle It is easy to check that
is a necessary condition for it, so, the ONL remains for this modified DFC. The use of Jury test (as carried out in [18]) should be useful, but, even out of ONL, it is neither possible to obtain the γ.’s explicitly nor ensure a range finding in the general case. However, a first attempt to look for a solution is setting γi /= 0 for a fixed i and
for all 
in order to simplify the equation In this way,
results a (m − 1)-multiplicity root of
, and the other two are the roots of a quadratic function
where
. From the conditions on A and B. so that 
[13] it is deduced that a range
exists if and only
This strategy is illustrated in the logistic map (Eq. (7)) for the stabilization of m-UPO
. The last condition becomes
(20)The range
is settled on
(21)For r0 = 3.62 and 4-UPO,
condition (20) is accomplished. Taking, for instance,
, a range for .. is obtained from (21). Figure 4 shows this UPO being stabilized by control (18) for initial
condition x0 = 0.5 and 
= 0.05. The performance of the control u. is shown in the same figure. This resource does not work for r0> 3.625 because the Floquet multiplier of the 4-UPO is smaller than -3, and Eq. (20) is not accomplished anymore. Even so, 4-tuples
can be found for all the roots of 
to be within the unit circle. Indeed, by means of detailed heuristic search on the coefficient values of 
a solution is obtained. Figure 5 shows an example where r0 = 3.67. It is worth noting that, in this case, the control strategy (18) overcomes the Pyragas method, which fails above . 3.62 (see [32, p. 50]). Namely, by control (18) the 4-UPO is stabilized for an enlargement of r0 values: 3.625 ≤ r0 ≤ 3.67. Moreover, the control (18) also works successfully even in presence of noise as it is exhibited in Fig. 6.
The extension of the Pyragas method to stabilize the system in a m-UPO, known as “ex- tended time-delay autosynchronization system” (ETDAS) or “extended delayed feedback control” (EDFC), is defined as
(22)with 0 ≤ R < 1, which reduces to TDAS for R = 0. It was introduced to stabilize higher instabilities. Indeed, going back to the previous example, it is reported in [32] that with R = 0.5, stabilization of the 4-UPO of the logistic map can be maintained up to the parameter value
.
Based on control (18), a kind of extension like (22) is proposed to stabilize a system in a m-UPO that consists in replacing it by
(23)System (2) and (23) yields to a 2m-dimensional system (with . switches) with
as its state variables. For exam- ple, for the logistic map, the m = 4 case is given by
(24)Note that the 4-UPO of the free system (γi = 0 for all i)
is preserved when (23) is applied. The Jacobian matrix of system (24) at
is, in this case, a 8 × 8 matrix given by
where
and
and
the characteristic polynomial is given by
Applying the heuristic algorithm to search on the coefficient values of
, 4-tuples,
are found, and the 4-UPO for 3.67 < r0 ≤ 3.8 can be stabilized. Thus, the range of parameter values for which the 4-cycle-logistic map is stabilized, is enlarged. Figure 7 illustrates the orbit stabilization and the control performance for r0 = 3.8.
This section points out the required adjustments for applying the developed strategies in the .-dimensional case, and even more, they are compared to interesting background in the literature.
For SPF modulation of Section 2.3, Eq. (9) remains the same save that control gain is
and the distance is given by the Euclidean norm in
. The Jacobian of the controlled system at each
is given by
where
and
matrices, respectively.
The insertion of a switching controller as part of a SPF method appears in previous bibliography. Indeed, the multipoint OGY formula (see [3, Sect. 3.2.6]) may be inter- preted as an extension of Eq. (6) to the .-dimensional case. In turn, some modifications to overcome its limitations have been developed [3]. A two-level control method is proposed in [5], where the distance between the trajectory and the UPO is minimized in each step. The authors claim that its implementation on UPO’s requires a number of switchings greater than . to obtain successful control. In [25], it is proposed to choose feedback gains . such that the whole of the eigenvalues of
. be of modulus less than one. Note that this criterion may fail to detect parameters βi. for which the Jacobean matrix of the composition has the whole of its eigenvalues within the unit circle.
The stability condition stated in Eq. (11) becomes
for the .-dimensional case so forcing the trajectory to be closer to the UPO. This is verified if, for all
have the maximum singular value of modulus less than one [15], which is accomplished
has a basis of eigenvectors and the whole of its eigenvalues within the unit circle. If rank
pole- placement method [23] is a useful tool to obtain
. satisfying this condition.
Then Proposition 1 is straightforward generalized to the n-dimensional case. On the other hand, the generalized version of Proposition 2 to the .-dimensional case requires the product
. to have the whole of its eigenvalues within the unit circle, which is a weaker condition, but it yields to nonlinear equations in βi. Defining
the results of Section 2 are accordingly valid in the .-dimensional case.
Similarly, in the DFC modulation defined in Section 3.2, a vector
becomes the control gain for the n-dimensional version of Eq. (15). The stability criterion involves the product
matrix given by
where
are both n x n matrices, In is the (n × n)-identity matrix and 0n is the (n × n)-null matrix.
In the stability condition of DFC studied in [19], for the multidimensional case, the system is affected by an additive feedback control that depends only on a constant gain, 
The insertion of m different control gains
in Section 3.2 improves the proposals in [18, 19, 32], where
for all i, and which hinges on the eventual scenario, where nonempty overlapping γ.-ranges comes out. In addition, the proof of stability outlined in the present work avoids the computation of the product of the matri- ces Ji, and the characteristic polynomial of Eq. (17) is obtained by a more straightforward calculation.
The EDFC algorithm of Eq. (22) may also be formulated for the n-dimensional case. Its stability analysis yields to a 2mn-dimensional system (with . switches). As in the SPF modulation, controllability tools become useful to obtain adequate control gains. Defining .
the arguments on the validity of these methods are also straightforward generalized from Section 3, although a full description of them becomes tiresome.
In most cases, the parameter . affects only one component of the system given by f , say fj. Then, when looking for βi ∈ R (respectively γi) in these situations, a conductive
simplification can be made by setting all of its components equal to zero, except the .th one. In turn, for DFC strategy, it means that less delayed states are introduced, yielding to a decrease on the controlled system dimension (namely, on the . matrix dimension). Moreover, the (m + 1)-multiplication effect on the controlled system is avoided in these particular situations. This is the case of the two-dimensional Hénon map, where the second component of the state variable is just the delayed first one [7]:
(25)It is well known that for a = a0 = 1.4 and b = 0.3, the map posses a one-piece chaotic attractor. This map may be stabilized on its 2-UPO within the attractor by tuning the parameter . around the nominal value a0. This 2-UPO√is given by
with
being 
and 
The implementation of DFC modulation on (25) (confront with Eq. (16)) yields to the system
(26)where
(here
The stability condition is obtained from the product-matrix . w here
by looking for .. and .. that make the whole eigenvalues of . be of modulus less than one. Figure 8 shows the 2-UPO stabilization by this strategy. The performance of the control u. . γ..x. . z.), . = 1, 2, is appreciated in the same figure.
This methodology results useful for controlling hyperchaos. The applicability of stabi- lization methods on hyperchaotic systems is pointed out in [9]. Here the successfulness of a DFC-type algorithm for the stabilization of hyperchaotic generalized Hénon maps at a fixed point is analyzed. The .-dimensional system is defined [9, 28]:
(27)where xi,krepresents the state i at time . for i = 1, . . . , n. Fixing the parameter value ., the map is hyperchaotic for certain parameter range of . since the number of positive Lyapunov exponents is n − 1 [28]. It has two fixed points, , with x*1,2 = (−b−1± (b + 1)2 + 4a)/2, being x*1 within the hyperchaotic attractor stabilize (27) on it, an ad-hoc reformulation of control (18) is proposed involving n – 1 parameters γi. The n = 4 case is taken for illustration. When control is activated, the system results
(28)Note that although the number of control parameters is augmented, the algorithm remains not invasive, and the system dimension is not affected. By working out algebraic features of the characteristic polynomial, adequate sets of parameters are found. Formulation and procedure in the general n-case come out straightly. Figure 9 shows stabilization of (27) at the fixed point x*1. for a0 = 1.76 and b = 0.1. Note that the bigger size of . and the number of iteration ., needed for the system to achieve the controlled regime, is according to the bigger size of the embedded dimension .. The performance of the control, uk = y1(xk - yk) + y2(yk -zk) + y3(zk - wk), is displayed in the same figure.
It is well known that the advantage of any DFC method over SPF methods is just that the full knowledge of the UPO to be stabilized is not needed, instead, limitations on con- vergence issues arise. In this work, both SPF and DFC methods have been revisited, and, based on them, new modifications yielding to a “switching” SPF- and DFC-type strategies have been proposed. Throughout the paper, the introduction of each modification has been fully argued, and the comparison to similar proposals previously published have been mentioned in detail. The convergence of the improved version has been rigorously proved and confronted to the original method without claiming their optimality over all other conceivable approaches. Additionally, control performance aspects, usually disregarded in published works dealing with DFC methods, have been also taken care of.
Each one of the developed control strategies consists of two stages depending on various control parameters. For a .-UPO, the activated control law depends on . control gains, so more data must be stocked for increasing ., but this helps to find a successful control strategy. Under stated hypothesis, UPO stabilization is assured by choosing each control gain within a certain range of values. Among them, it is plausible to look for those values, which make the control strategy verify required conditions on waiting time and control magnitude bound. The right choice of the parameter . is crucial, and a bound on it is stated for which convergence and desired control magnitude are guaranteed.
Even though the improved DFC does not overcome the ONL, it does work success- fully when applied to an UPO presenting high instability with negative Floquet multipli- ers (and for which DFC fails) if adequate control gains γ. are found. The insertion of . control gains γ. together with the parameter control . into the EDFC method is also interesting: for a given UPO, the modified EDFC should require a smaller . than the original EDFC and so yielding to a greater speed of convergence.
The implementation of the proposed control strategies to stabilize any .-dimensional discrete-time systems in any .-UPO is also stated. For generalized Hénon-like systems, an ad-hoc control modulation is designed and applied to stabilize the hyperchaotic motion on a fixed point within its attractor. The hyperchaotic behavior is harder to be controlled than the classical chaos because more than one of the Lyapunov exponents are positive, which means more instability directions for dynamics to spread across the whole phase- space. This fact is critical in the .-D generalized Hénon maps case where the number of such instabilities is just n 1 [28], yielding a longer waiting time to achieve the controlled regime as shown in Fig. 9.
The feasibility and applicability of the these novel control modulation methods are all validated by the good simulation results.
Future work aims to introduce this kind of modifications into control schemes based on oscillating or periodic gain control or on predictive states. As the ONL restriction is overcome by these methods, improvements on their convergence and control parameter features result an attracting task.
The authors thank the reviewer for the suggestions and observations that improved the readability of this work.
