# Computing and Controlling Basins of Attraction in Multistability Scenarios

Authors: John Alexander Taborda; Fabiola  Angulo
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

The aim of this paper is to describe and prove a new method to compute and control the basins of attraction in multistability scenarios and guarantee monostability condition. In particular, the basins of attraction are computed only using a submap, and the coexistence of periodic solutions is controlled through fixed-point inducting control technique, which has been successfully used until now to stabilize unstable periodic orbits. In this paper, however, fixed-point inducting control is used to modify the domains of attraction when there is coexistence of attractors. In order to apply the technique, the periodic orbit whose basin of attraction will be controlled must be computed. Therefore, the fixed-point inducting control is used to stabilize one of the periodic orbits and enhance its basin of attraction. Then, using information provided by the unstable periodic orbits and basins of attractions, the minimum control effort to stabilize the target periodic orbit in all desired ranges is computed. The applicability of the proposed tools is illustrated through two different coupled logistic maps.


---

## Body

## 1. Introduction

Complex bifurcation scenarios have been observed in nonlinear dynamic systems from virtually all areas of science, including a broad range of natural sciences, mechanical and electrical engineering, and economics and other areas of the social sciences [1–3]. Theoretical and applied researches have explained various bifurcation scenarios [4–6], and analytical, numerical, and experimental works have contributed to unraveling the complexity inherent to chaotic motion [3].Coupled chaotic maps are a set of special discrete-time dynamical systems that can describe chemical, epidemiological, physiological, biological, or engineering systems [7]. Interesting nonlinear phenomena have been reported in them. For example, in [8, 9], coupled logistic maps were analyzed and new scenarios for transition to chaos were found via the creation and destruction of multilayered tori. In those papers, novel routes to chaos were described, and authors found that, depending on the coupling constant value, the system approaches different periodic attractors.Control methods of coexisting attractors in multistability scenarios have been widely studied in the last decades. In [10], periodic signals were replaced by chaotic ones in order to eliminate multiple domains of attraction. In [11, 12], the influence of noise on preference and dominance of attractors in multistable systems was studied. In [13, 14], multistability was controlled using small periodic modulation of a system parameter. In [15], the basins of attraction (BA) were controlled using harmonic and stochastic modulation. In [16], an impulsive force was used in order to perturb one attractor of the system and to change its response such that the system response evolved to another attractor. In [17], control of multistability scenarios based on the selection of a particular attractor by periodic external modulation was presented. A complete report of control of multistability can be found in [18].In this paper, we prove a different technique to control domains of attraction in multistability scenarios which is called Fixed-Point Inducting Controller (FPIC). The FPIC is a feed-forward controller that forces the system solution to evolve to an existing desired attractor which cannot be always reached for the uncontrolled case because of initial conditions. This technique was initially proposed in [19] and successfully applied in [20, 21]. In particular, in [21], experiments showed the good performance of the controlled system, where the FPIC was used as a second control loop. However, in all previous works published until now, FPIC has been only used to stabilize unstable periodic orbits.The main contributions of this paper lie in numerical analysis and control design areas. The development of a novel methodology to compute and analyze the basins of attraction reinforced the numerical analysis. The methodology consists in decomposing the map of the system inp submaps, where p corresponds to the order of the periodic orbit to be controlled. To have the BA, the long-time response terms of this submap are depicted according to key colours. Apart from help to the control design, this way to compute the basin of attraction can be seen as the basins of attraction of period-1 orbits which are computed every sixth or fifth iteration for linear and nonlinear coupling, respectively, and the control goal can be thought of as the control of period-1 orbit. On the other hand, the methodology to design the controller based on FPIC technique belongs to control theory. When the proposed controller is used, the system is forced to evolve to a known periodic orbit that exists in the uncontrolled map. Moreover, by using the information obtained from bifurcation diagrams and basins of attraction, it is possible to compute the minimum control effort required to stabilize the target orbit in the defined region. In particular, the coexistence of periodic solutions in coupled logistic maps [8, 9, 22] is controlled by widening the basin of attraction of a period-p orbit that coexists with another one, and the minimum control effort is computed aided by the unstable period-2 orbit for the linear coupling and by unstable period-1 orbits for the nonlinear coupling case.The paper is organized as follows. In Section2, the coupled logistic maps are presented. The coexistence of period-6 and period-5 orbits is identified. In Section 3, the methodology to compute the basins of attraction is explained, and the procedure is applied to the linear and nonlinear coupling maps. In Section 4, the methodology to compute the controller is presented and applied to the considered systems and the minimum control effort required to guarantee monostabilization of the periodic orbit is computed. In Section 5, a brief discussion of the results is presented, and finally, in Section 6, the conclusions are given.


## 2. Coupled Logistic Maps

Coupled chaotic maps provide a source of bifurcation scenarios with nonlocal phenomena and coexistence of attractors. To design the proposed controller, two maps with different coupling mechanisms are chosen. For linear coupling map, there are two period-6 orbits that coexist in a range of the parameter set. Similarly, for nonlinear coupling case, two period-5 orbits coexist.


### 2.1. Coupled Logistic Maps: Linear Coupling

This system is described by(1)xk+1=fxk+ϵyk-xk,yk+1=fyk+ϵxk-yk,where ϵ∈(-2,2) is the coupling constant and fz(k)=rz(1-z), r∈(1,4), is the parameter associated with the nonlinear part of (1). Since the system is symmetrical, there is an invariant line Δ and the restriction of the 2D map to Δ reduces it to a 1D map (the logistic map). The dynamics on this invariant set help us to study the dynamics and bifurcations of the 2D system. Moreover, symmetrical trajectories are generated by symmetrical initial conditions leading to symmetrical basins of attraction [22]. Novel routes to chaos through torus-breakdown mechanism of this system were reported in [9] and different dynamics were characterized in the parameter region given by (ϵ,r).Figure1 shows the dynamic behavior of (1) when r is varied and two symmetrical initial conditions are considered. Note that a pair of coexisting period-6 orbits are identified in these diagrams for r∈[3.045,3.065], which makes the main difference between Figures 1(a) and 1(b). The dynamical behavior in the rest of the interval is very similar and other differences cannot be easily seen. Table 1 shows the specific values of the two coexisting period-6 orbits. It can be observed that the solutions are mutually symmetrical in the sense that (x,y)=(y,x).Table 1
Two coexisting period-6 orbits for the linear coupling maps whenϵ=0.2725 and r=3.065. For the sake of simplicity, the periodic orbit values are presented with 8 decimal digits.
First period-6 orbit
Second period-6 orbit
(
x
∗
,
y
∗
)
1
=
(
0.79744636,0.20236605
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.20236605,0.79744636
)
(
x
∗
,
y
∗
)
2
=
(
0.33291675,0.65689339
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.65689339,0.33291675
)
(
x
∗
,
y
∗
)
3
=
(
0.76896860,0.60251974
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.60251974,0.76896860
)
(
x
∗
,
y
∗
)
4
=
(
0.49915797,0.77939325
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.77939325,0.49915797
)
(
x
∗
,
y
∗
)
5
=
(
0.84261193,0.45063017
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.45063017,0.84261193
)
(
x
∗
,
y
∗
)
6
=
(
0.29965625,0.86559446
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.86559446,0.29965625
)Figure 1
Bifurcation diagrams of (1) with ϵ=0.2725 and two symmetrical initial conditions.
(a)
(
x
0
,
y
0
)
=
(
0.1,0.5
)
(b)
(
x
0
,
y
0
)
=
(
0.5,0.1
)


### 2.2. Coupled Logistic Maps: Nonlinear Coupling

This system is defined by(2)xk+1yk+1=fxk+ϵfyk-fxkfyk+ϵfxk-fyk,where fx(k) and fy(k) were defined before, ϵ∈(-2,2), and r∈(1,4). Figure 2 shows an interesting coexistence scenario for this system. Two mutually symmetrical period-5 cycles coexist in it when ϵ=1.5 and r∈[2.775,2.815].Figure 2
Bifurcation diagrams of (2) with ϵ=1.5 and two symmetrical initial conditions.
(a)
(b)


## 2.1. Coupled Logistic Maps: Linear Coupling

This system is described by(1)xk+1=fxk+ϵyk-xk,yk+1=fyk+ϵxk-yk,where ϵ∈(-2,2) is the coupling constant and fz(k)=rz(1-z), r∈(1,4), is the parameter associated with the nonlinear part of (1). Since the system is symmetrical, there is an invariant line Δ and the restriction of the 2D map to Δ reduces it to a 1D map (the logistic map). The dynamics on this invariant set help us to study the dynamics and bifurcations of the 2D system. Moreover, symmetrical trajectories are generated by symmetrical initial conditions leading to symmetrical basins of attraction [22]. Novel routes to chaos through torus-breakdown mechanism of this system were reported in [9] and different dynamics were characterized in the parameter region given by (ϵ,r).Figure1 shows the dynamic behavior of (1) when r is varied and two symmetrical initial conditions are considered. Note that a pair of coexisting period-6 orbits are identified in these diagrams for r∈[3.045,3.065], which makes the main difference between Figures 1(a) and 1(b). The dynamical behavior in the rest of the interval is very similar and other differences cannot be easily seen. Table 1 shows the specific values of the two coexisting period-6 orbits. It can be observed that the solutions are mutually symmetrical in the sense that (x,y)=(y,x).Table 1
Two coexisting period-6 orbits for the linear coupling maps whenϵ=0.2725 and r=3.065. For the sake of simplicity, the periodic orbit values are presented with 8 decimal digits.
First period-6 orbit
Second period-6 orbit
(
x
∗
,
y
∗
)
1
=
(
0.79744636,0.20236605
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.20236605,0.79744636
)
(
x
∗
,
y
∗
)
2
=
(
0.33291675,0.65689339
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.65689339,0.33291675
)
(
x
∗
,
y
∗
)
3
=
(
0.76896860,0.60251974
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.60251974,0.76896860
)
(
x
∗
,
y
∗
)
4
=
(
0.49915797,0.77939325
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.77939325,0.49915797
)
(
x
∗
,
y
∗
)
5
=
(
0.84261193,0.45063017
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.45063017,0.84261193
)
(
x
∗
,
y
∗
)
6
=
(
0.29965625,0.86559446
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.86559446,0.29965625
)Figure 1
Bifurcation diagrams of (1) with ϵ=0.2725 and two symmetrical initial conditions.
(a)
(
x
0
,
y
0
)
=
(
0.1,0.5
)
(b)
(
x
0
,
y
0
)
=
(
0.5,0.1
)


## 2.2. Coupled Logistic Maps: Nonlinear Coupling

This system is defined by(2)xk+1yk+1=fxk+ϵfyk-fxkfyk+ϵfxk-fyk,where fx(k) and fy(k) were defined before, ϵ∈(-2,2), and r∈(1,4). Figure 2 shows an interesting coexistence scenario for this system. Two mutually symmetrical period-5 cycles coexist in it when ϵ=1.5 and r∈[2.775,2.815].Figure 2
Bifurcation diagrams of (2) with ϵ=1.5 and two symmetrical initial conditions.
(a)
(b)


## 3. Methodology to Compute the Basins of Attraction

In this section, a methodology to compute the basins of attraction (BA) is developed. By using this methodology, it is possible to find the minimum control effort required to stabilize a desired periodic orbit.


### 3.1. Computation of the Basins of Attractions in Linear Coupling Maps

The BA corresponds to the set of initial conditions whose long-time response (LTR) approaches the attractor. The first and second period-6 orbits were defined in Table1. All numerical results associated with the second period-6 orbit are differentiated from the first one with an upper bar. In order to compute the BA and after applying FPIC control, we associate with each point k of the first period-6 orbit a section Bk that is given by |xLTR-x∗|<δ, such as what is illustrated in Figure 3. A similar procedure is developed to compute the BA of the second period-6 orbit, whose sections are defined as B¯k.Figure 3
Colour key to compute the BA.To compute the LTR terms, we start from any initial condition inside the considered region (0<x1<1 and 0<x2<1) and compute its successive iterations until the solution approaches the attractor. The general map to describe the evolution of the system is P, and it is given by (3)P≔x0⟼x1⟼x2⟼x3⟼x4⟼x5⟼x6⟼x7⟼x8⟼x9⟼⋯xm,where m is high enough. The map P is decomposed in six submaps Pi. Each Pi is formed by the set of points given by x(i-1+6j), where i=1⋯6 and j∈Z++{0}. If the LTR terms of Pi lie inside some Bk (or B¯k when the BA for the second period-6 orbit is computed), it is coloured according to the colour key displayed in Figure 3. If the LTR terms of Pi lie outside these regions, then it is gray coloured. This procedure is repeated for a lot of initial conditions.The BA of the first period-6 orbit (second period-6 orbit) can be computed and depicted only usingPi (P¯i) submaps for a fixed i, and its diagram may have up to seven colors. The gray coloured points correspond to the set of initial conditions whose LTR terms diverge from the first period-6 orbit (second period-6 orbit), and the other points (depicted with six different colors according to Figure 3) correspond to the initial conditions whose LTR terms lie inside some Bk (B¯k) region. The BA computed for the first period-6 orbit (second period-6 orbit) in this way is noted as BA(Pi) (BA(P¯i)). Figure 4(a) shows the diagram associated with BA(P1). All BA(Pi) have the same information but are coloured in different way. Figure 4(b) shows BA(P¯1). These BA let us compute the minimum control effort required to stabilize the periodic orbits.Figure 4
(a)BA(P1). (b) BA(P¯1). r=3.065 and ϵ=0.2725.
(a)
(b)


### 3.2. Computation of the Basins of Attraction for Nonlinear Coupling Maps

The procedure explained before is applied to the nonlinear coupling maps.BA(Pi) and BA(P¯i) of the two period-5 orbits are also symmetrical. Figures 5(a) and 5(b) show BA(P1) and BA(P¯1), respectively. The remaining diagrams BA(Pi) and BA(P¯i) for i=2⋯5 have similar behavior. As we said before, all BA(Pi) carry the same information but are coloured in different way.Figure 5
(a)BA(P1). (b) BA(P¯1). r=2.8 and ϵ=1.5.
(a)
(b)


## 3.1. Computation of the Basins of Attractions in Linear Coupling Maps

The BA corresponds to the set of initial conditions whose long-time response (LTR) approaches the attractor. The first and second period-6 orbits were defined in Table1. All numerical results associated with the second period-6 orbit are differentiated from the first one with an upper bar. In order to compute the BA and after applying FPIC control, we associate with each point k of the first period-6 orbit a section Bk that is given by |xLTR-x∗|<δ, such as what is illustrated in Figure 3. A similar procedure is developed to compute the BA of the second period-6 orbit, whose sections are defined as B¯k.Figure 3
Colour key to compute the BA.To compute the LTR terms, we start from any initial condition inside the considered region (0<x1<1 and 0<x2<1) and compute its successive iterations until the solution approaches the attractor. The general map to describe the evolution of the system is P, and it is given by (3)P≔x0⟼x1⟼x2⟼x3⟼x4⟼x5⟼x6⟼x7⟼x8⟼x9⟼⋯xm,where m is high enough. The map P is decomposed in six submaps Pi. Each Pi is formed by the set of points given by x(i-1+6j), where i=1⋯6 and j∈Z++{0}. If the LTR terms of Pi lie inside some Bk (or B¯k when the BA for the second period-6 orbit is computed), it is coloured according to the colour key displayed in Figure 3. If the LTR terms of Pi lie outside these regions, then it is gray coloured. This procedure is repeated for a lot of initial conditions.The BA of the first period-6 orbit (second period-6 orbit) can be computed and depicted only usingPi (P¯i) submaps for a fixed i, and its diagram may have up to seven colors. The gray coloured points correspond to the set of initial conditions whose LTR terms diverge from the first period-6 orbit (second period-6 orbit), and the other points (depicted with six different colors according to Figure 3) correspond to the initial conditions whose LTR terms lie inside some Bk (B¯k) region. The BA computed for the first period-6 orbit (second period-6 orbit) in this way is noted as BA(Pi) (BA(P¯i)). Figure 4(a) shows the diagram associated with BA(P1). All BA(Pi) have the same information but are coloured in different way. Figure 4(b) shows BA(P¯1). These BA let us compute the minimum control effort required to stabilize the periodic orbits.Figure 4
(a)BA(P1). (b) BA(P¯1). r=3.065 and ϵ=0.2725.
(a)
(b)


## 3.2. Computation of the Basins of Attraction for Nonlinear Coupling Maps

The procedure explained before is applied to the nonlinear coupling maps.BA(Pi) and BA(P¯i) of the two period-5 orbits are also symmetrical. Figures 5(a) and 5(b) show BA(P1) and BA(P¯1), respectively. The remaining diagrams BA(Pi) and BA(P¯i) for i=2⋯5 have similar behavior. As we said before, all BA(Pi) carry the same information but are coloured in different way.Figure 5
(a)BA(P1). (b) BA(P¯1). r=2.8 and ϵ=1.5.
(a)
(b)


## 4. Methodology to Control the Basins of Attraction

In this section, domains of attraction in coexistence scenarios are controlled using FPIC technique. This technique has been proven to stabilize unstable and chaotic systems [19–21]. In those papers, analytical or numerical values of the steady-state control inputs were needed to guarantee the stabilization of the equilibrium point because FPIC was used as a second control loop. On the contrary, in this work, FPIC technique is used to control the BA in coupled logistic maps to guarantee stabilization of period-6 (for linear coupling) and period-5 (for nonlinear coupling) orbits, using the values of the periodic orbit to be controlled.


### 4.1. Fixed-Point Inducting Control

Theorem 1.
Consider a discrete-time system(4)zk+1=gzk.It is assumed that period-p solutions coexist in the system and two or more of them are stable. A stable period-p solution is noted in a compact form as z∗, and it is composed of p points such that it can be expressed as (5)z∗=z1∗,z2∗,…,zp∗=x1∗,x2∗,…,xn∗1,x1∗,x2∗,…,xn∗2,…,x1∗,x2∗,…,xn∗p.The objective of the controller is to expand the domain of attraction of a given period-p orbit, by forcing the system to evolve to this specific solution. Equation (6) shows the control law based on FPIC, with E as the control parameter. This controller stabilizes the target periodic orbit and expands its attraction domain as E increases. Hence,(6)
z
k
+
1
=
1
E
+
1
g
z
k
+
E
E
+
1
z
j
∗
i
f
k
mod
⁡
p
=
0
g
z
k
i
f
k
mod
⁡
p
≠
0
,where j is defined depending on the analysis of the controlled system.Proof.
The periodic orbits of the controlled system are the same as the uncontrolled one. This is because afterp iterations g(zj∗)=zj∗ (by definition of periodic orbit) and thus considering only the upper part of (6), it can be seen that the periodic orbits remain unchanged. Hence, when E is high enough, the controller transfers the state instantly to some Bk when kmod⁡p=0, irrespective of the state value. In this way, the new state (that is near to the periodic solution) is the initial condition for the next iteration obtaining convergence to the periodic solution because the periodic orbit is stable.Remarks. (i) From a practical point of view, the idea is not to fix E to a high value because the control effort increases as E does. To control systems whose initial conditions are away from the periodic orbit can lead to oversizing of the control effort. For example, for the systems analyzed in this work, the solution quickly diverges when the initial conditions are outside the considered region and the controller cannot be used.(ii) The minimum value ofE that should be put in (6) can be found based on numerical analysis, where the unstable period orbits that coexist in the dynamical system play an important role.


### 4.2. Control of Period-6 Orbits in Linear Coupling Maps

Let us consider the map defined by (1) with r=3.065 and ϵ=0.2725. We will stabilize the second period-6 orbit using FPIC control. In order to achieve the goal, we use the second point (x¯∗,y¯∗)2 as the target control; that is, x¯2∗≔(0.6569,0.3329). Taking into account (6), the linear coupling map is perturbed every sixth iteration according to the control law defined by(7)
x
k
+
1
y
k
+
1
=
1
E
+
1
f
x
k
+
ε
y
k
-
x
k
f
y
k
+
ε
x
k
-
y
k
+
E
E
+
1
x
¯
∗
y
¯
∗
if
k
mod
⁡
6
=
0
f
x
k
+
ε
y
k
-
x
k
f
y
k
+
ε
x
k
-
y
k
if
k
mod
⁡
6
≠
0
.Figure6 shows BA(P¯1) diagram for different E values. Colors were assigned as defined above. The elements of BA(P¯1) that tends to B¯2 increase as parameter E increases. Comparing Figures 4(b) and 6, it can be seen that the green coloured region enhances forming clear narrow contours. Now, the BA is only formed by two sets of points. The first one includes all points whose LTR terms end in B¯2 which are green coloured. The second one includes all points whose LTR terms are not inside any B¯i section for i=1,2,3,4,5,6, which are gray coloured. This implies that no long-time term of the submap P¯1 lies in B¯i section, for i=1,3,4,5,6. The green contours expand as E increases, widening the BA and diminishing the gray coloured parts.Figure 6
(a)BA(P¯1) for E=0.85. (b) BA(P¯1) for E=1.5. (c) BA(P¯1) for E=1.7. r=3.065 and ϵ=0.2725 in all cases.
(a)
(b)
(c)Figure7 shows BA(P1) when x2∗≔(0.3329,0.6569) is used as the reference in the FPIC controller. The same values of E have been used.Figure 7
(a)BA(P1) for E=0.85. (b) BA(P1) for E=1.5. (c) BA(P1) for E=1.7. r=3.065 and ϵ=0.2725 for all cases.
(a)
(b)
(c)


#### 4.2.1. Selection of Minimum Control Effort

The minimum control effortEcr is defined as the minimum value of E that forces the system to evolve to the controlled orbit. The selection of Ecr can be computed by studying the behavior of unstable periodic orbits (UPOs) around the periodic orbits.The metamorphosis of controlledBA(P¯1) or BA(P1) is influenced by two unstable period-2 solutions of the system that are schematized in Figure 8. The fixed points of these period-2 orbits are named V1, V2, V3, and V4.Figure 8
Sketch of unstable period-2 solutionsVi in coexistence zone.For the controlled system, the pointsVi are hard to cover by the expansion of BA(P1). These points are always on the frontiers between the green and gray regions as it is shown in Figure 9. To find Ecr, we consider Δi region around each point Vi. Ecr is the minimum value such that all points in Δi0 regions converge to Vi points. Figure 9 shows examples of Δi regions when E=1.5. The neighborhood near Vi requires more control effort to converge to B2.Figure 9
V
i points and Δi sections in controlled linear map with E=1.5. (a) V1=(0.81326798,0.33518226). (b) V2=(0.33518226,0.81326798). (c) V3=(-0.07696298,0.57288468).
(a)
(b)
(c)Δ
i can be expressed as (Δx,Δy)≔(υx+δsin⁡θ,υy+δcos⁡θ), where δ is fixed to a small value and θ is varied in the range [0,2π]. Figure 10 illustrates the behavior of BA(P1) around Δi when E is varied. The convergence of all points to Δi is reached when Ecr≈1.78.Figure 10
Selection of critical valueEcr. Diagram E versus θ to select minimum value of E that guarantees total convergence of P1 to x∗=0.3329 and y∗=0.6568.
(a)
V
1
(b)
V
2
(c)
V
3


#### 4.2.2. Selection of the Reference

Similar results can be obtained if we choose other equilibrium points to perturb the mapPi. Table 3 summarizes the control effort requirements. Note that the critical control effort Ecr varies depending on the used reference. Coexistence phenomenon can be controlled with less effort if (x∗,y∗)=(0.8426,0.4506) or (x¯∗,y¯∗)=(0.4506,0.8426) are used to control the map.


### 4.3. Control of Period-5 Orbits in Nonlinear Coupling Maps

Now, we control two coexisting period-5 cycles of (2) when r=2.8 and ϵ=1.5. The reference is introduced every fifth iteration as is presented in(8)
x
k
+
1
y
k
+
1
=
1
E
+
1
f
x
k
+
ε
f
y
k
-
f
x
k
f
y
k
+
ε
f
x
k
-
f
y
k
+
E
E
+
1
x
¯
∗
y
¯
∗
if
k
mod
⁡
5
=
0
f
x
k
+
ε
f
y
k
-
f
x
k
f
y
k
+
ε
f
x
k
-
f
y
k
if
k
mod
⁡
5
≠
0
.According to Table2, we use (x¯∗,y¯∗)1≔(0.3526,0.7918) to control the system. Figure 11 presents BA(P¯1) of the controlled system when three different values of E are considered.Table 2
Two coexisting period-5 orbits for nonlinear coupling maps whenϵ=1.5 and r=2.8. For the sake of simplicity, the periodic orbit values are presented with 8 decimal digits.
First period-5 orbit
Second period-5 orbit
(
x
∗
,
y
∗
)
1
=
(
0.79175845,0.35258819
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.35258819,0.79175845
)
(
x
∗
,
y
∗
)
2
=
(
0.72790518,0.37290574
)
(
x
-
∗
,
y
-
∗
)
2
=
(
0.37290574,0.72790518
)
(
x
∗
,
y
∗
)
3
=
(
0.70487469,0.50446287
)
(
x
-
∗
,
y
-
∗
)
3
=
(
0.50446287,0.70487469
)
(
x
∗
,
y
∗
)
4
=
(
0.75867944,0.52373859
)
(
x
-
∗
,
y
-
∗
)
4
=
(
0.52373859,0.75867944
)
(
x
∗
,
y
∗
)
5
=
(
0.79131428,0.41974570
)
(
x
-
∗
,
y
-
∗
)
5
=
(
0.41974570,0.79131428
)Table 3
Summary of FPIC requirements in coexistence control of two period-6 cycles.
Maps
x
∗
x
-
∗
E
cr
P
1
,
P
-
1
(
0.7974,0.2023
)
(0.2023, 0.7974)
3.371
P
2
,
P
-
2
(
0.3329,0.6568
)
(0.6568, 0.3329)
1.774
P
3
,
P
-
3
(
0.7689,0.6025
)
(0.6025, 0.7689)
2.507
P
4
,
P
-
4
(
0.4991,0.7793
)
(0.7793, 0.4991)
1.492
P
5
,
P
-
5
(
0.8426,0.4506
)
(0.4506, 0.8426)
1.258
P
6
,
P
-
6
(
0.2996,0.8655
)
(0.8655, 0.2996)
4.112Figure 11
B
A
(
P
¯
1
) for E=0.25. (b) BA(P¯1) for E=0.45. (c) BA(P¯1) for E=0.65. r=2.8 and ϵ=1.5 for all cases.
(a)
(b)
(c)Selection of Minimum Control Effort. As in the previous case, the minimum control effort Ecr can be computed by studying the behavior of UPOs around the periodic orbits. Figure 12 shows an illustrative sketch of bifurcation scenario of the uncontrolled system. The unstable period-1 orbits are noted as N1 and N2. The metamorphosis of the controlled maps BA(P¯i) is influenced by unstable period-1 solutions of the system. We use points N1 and N2 to compute Ecr in each case. Table 4 presents the minimum value of parameter E to control coexistence phenomenon. All initial conditions belonging to the considered region converge to x¯∗ when E>Ecr. The lowest control effort that avoids coexistence phenomenon is reached when the map is controlled with the first point of the periodic orbit. If we induce the control using other points of the periodic orbit, then we need a greater control effort. Figure 13 shows the symmetrical behavior between BA(P1) and BA(P¯1) (compare with Figure 11). Critical control effort required in both cases has the same value.Table 4
Summary of FPIC requirements in coexistence control of two period-5 cycles.
Maps
x
∗
x
-
∗
E
cr
P
1
,
P
-
1
(
0.7918,0.3526
)
(0.3526,   0.7918)
0.6099
P
2
,
P
-
2
(
0.7279,0.3729
)
(0.3729,   0.7279)
0.6873
P
3
,
P
-
3
(
0.7049,0.5045
)
(0.5045,   0.7048)
1.2312
P
4
,
P
-
4
(
0.7588,0.5237
)
(0.5237,   0.7588 )
1.1892
P
5
,
P
-
5
(
0.7913,0.4197
)
(0.4197,   0.7913)
0.7504Figure 12
Sketch of unstable period-1 solutionsNi in coexistence zone.Figure 13
B
A
(
P
1
) for E=0.25. (b) BA(P1) for E=0.45. (c) BA(P1) for E=0.65. r=2.8 and ϵ=1.5 for all cases.
(a)
(b)
(c)


## 4.1. Fixed-Point Inducting Control

Theorem 1.
Consider a discrete-time system(4)zk+1=gzk.It is assumed that period-p solutions coexist in the system and two or more of them are stable. A stable period-p solution is noted in a compact form as z∗, and it is composed of p points such that it can be expressed as (5)z∗=z1∗,z2∗,…,zp∗=x1∗,x2∗,…,xn∗1,x1∗,x2∗,…,xn∗2,…,x1∗,x2∗,…,xn∗p.The objective of the controller is to expand the domain of attraction of a given period-p orbit, by forcing the system to evolve to this specific solution. Equation (6) shows the control law based on FPIC, with E as the control parameter. This controller stabilizes the target periodic orbit and expands its attraction domain as E increases. Hence,(6)
z
k
+
1
=
1
E
+
1
g
z
k
+
E
E
+
1
z
j
∗
i
f
k
mod
⁡
p
=
0
g
z
k
i
f
k
mod
⁡
p
≠
0
,where j is defined depending on the analysis of the controlled system.Proof.
The periodic orbits of the controlled system are the same as the uncontrolled one. This is because afterp iterations g(zj∗)=zj∗ (by definition of periodic orbit) and thus considering only the upper part of (6), it can be seen that the periodic orbits remain unchanged. Hence, when E is high enough, the controller transfers the state instantly to some Bk when kmod⁡p=0, irrespective of the state value. In this way, the new state (that is near to the periodic solution) is the initial condition for the next iteration obtaining convergence to the periodic solution because the periodic orbit is stable.Remarks. (i) From a practical point of view, the idea is not to fix E to a high value because the control effort increases as E does. To control systems whose initial conditions are away from the periodic orbit can lead to oversizing of the control effort. For example, for the systems analyzed in this work, the solution quickly diverges when the initial conditions are outside the considered region and the controller cannot be used.(ii) The minimum value ofE that should be put in (6) can be found based on numerical analysis, where the unstable period orbits that coexist in the dynamical system play an important role.


## 4.2. Control of Period-6 Orbits in Linear Coupling Maps

Let us consider the map defined by (1) with r=3.065 and ϵ=0.2725. We will stabilize the second period-6 orbit using FPIC control. In order to achieve the goal, we use the second point (x¯∗,y¯∗)2 as the target control; that is, x¯2∗≔(0.6569,0.3329). Taking into account (6), the linear coupling map is perturbed every sixth iteration according to the control law defined by(7)
x
k
+
1
y
k
+
1
=
1
E
+
1
f
x
k
+
ε
y
k
-
x
k
f
y
k
+
ε
x
k
-
y
k
+
E
E
+
1
x
¯
∗
y
¯
∗
if
k
mod
⁡
6
=
0
f
x
k
+
ε
y
k
-
x
k
f
y
k
+
ε
x
k
-
y
k
if
k
mod
⁡
6
≠
0
.Figure6 shows BA(P¯1) diagram for different E values. Colors were assigned as defined above. The elements of BA(P¯1) that tends to B¯2 increase as parameter E increases. Comparing Figures 4(b) and 6, it can be seen that the green coloured region enhances forming clear narrow contours. Now, the BA is only formed by two sets of points. The first one includes all points whose LTR terms end in B¯2 which are green coloured. The second one includes all points whose LTR terms are not inside any B¯i section for i=1,2,3,4,5,6, which are gray coloured. This implies that no long-time term of the submap P¯1 lies in B¯i section, for i=1,3,4,5,6. The green contours expand as E increases, widening the BA and diminishing the gray coloured parts.Figure 6
(a)BA(P¯1) for E=0.85. (b) BA(P¯1) for E=1.5. (c) BA(P¯1) for E=1.7. r=3.065 and ϵ=0.2725 in all cases.
(a)
(b)
(c)Figure7 shows BA(P1) when x2∗≔(0.3329,0.6569) is used as the reference in the FPIC controller. The same values of E have been used.Figure 7
(a)BA(P1) for E=0.85. (b) BA(P1) for E=1.5. (c) BA(P1) for E=1.7. r=3.065 and ϵ=0.2725 for all cases.
(a)
(b)
(c)


### 4.2.1. Selection of Minimum Control Effort

The minimum control effortEcr is defined as the minimum value of E that forces the system to evolve to the controlled orbit. The selection of Ecr can be computed by studying the behavior of unstable periodic orbits (UPOs) around the periodic orbits.The metamorphosis of controlledBA(P¯1) or BA(P1) is influenced by two unstable period-2 solutions of the system that are schematized in Figure 8. The fixed points of these period-2 orbits are named V1, V2, V3, and V4.Figure 8
Sketch of unstable period-2 solutionsVi in coexistence zone.For the controlled system, the pointsVi are hard to cover by the expansion of BA(P1). These points are always on the frontiers between the green and gray regions as it is shown in Figure 9. To find Ecr, we consider Δi region around each point Vi. Ecr is the minimum value such that all points in Δi0 regions converge to Vi points. Figure 9 shows examples of Δi regions when E=1.5. The neighborhood near Vi requires more control effort to converge to B2.Figure 9
V
i points and Δi sections in controlled linear map with E=1.5. (a) V1=(0.81326798,0.33518226). (b) V2=(0.33518226,0.81326798). (c) V3=(-0.07696298,0.57288468).
(a)
(b)
(c)Δ
i can be expressed as (Δx,Δy)≔(υx+δsin⁡θ,υy+δcos⁡θ), where δ is fixed to a small value and θ is varied in the range [0,2π]. Figure 10 illustrates the behavior of BA(P1) around Δi when E is varied. The convergence of all points to Δi is reached when Ecr≈1.78.Figure 10
Selection of critical valueEcr. Diagram E versus θ to select minimum value of E that guarantees total convergence of P1 to x∗=0.3329 and y∗=0.6568.
(a)
V
1
(b)
V
2
(c)
V
3


### 4.2.2. Selection of the Reference

Similar results can be obtained if we choose other equilibrium points to perturb the mapPi. Table 3 summarizes the control effort requirements. Note that the critical control effort Ecr varies depending on the used reference. Coexistence phenomenon can be controlled with less effort if (x∗,y∗)=(0.8426,0.4506) or (x¯∗,y¯∗)=(0.4506,0.8426) are used to control the map.


## 4.2.1. Selection of Minimum Control Effort

The minimum control effortEcr is defined as the minimum value of E that forces the system to evolve to the controlled orbit. The selection of Ecr can be computed by studying the behavior of unstable periodic orbits (UPOs) around the periodic orbits.The metamorphosis of controlledBA(P¯1) or BA(P1) is influenced by two unstable period-2 solutions of the system that are schematized in Figure 8. The fixed points of these period-2 orbits are named V1, V2, V3, and V4.Figure 8
Sketch of unstable period-2 solutionsVi in coexistence zone.For the controlled system, the pointsVi are hard to cover by the expansion of BA(P1). These points are always on the frontiers between the green and gray regions as it is shown in Figure 9. To find Ecr, we consider Δi region around each point Vi. Ecr is the minimum value such that all points in Δi0 regions converge to Vi points. Figure 9 shows examples of Δi regions when E=1.5. The neighborhood near Vi requires more control effort to converge to B2.Figure 9
V
i points and Δi sections in controlled linear map with E=1.5. (a) V1=(0.81326798,0.33518226). (b) V2=(0.33518226,0.81326798). (c) V3=(-0.07696298,0.57288468).
(a)
(b)
(c)Δ
i can be expressed as (Δx,Δy)≔(υx+δsin⁡θ,υy+δcos⁡θ), where δ is fixed to a small value and θ is varied in the range [0,2π]. Figure 10 illustrates the behavior of BA(P1) around Δi when E is varied. The convergence of all points to Δi is reached when Ecr≈1.78.Figure 10
Selection of critical valueEcr. Diagram E versus θ to select minimum value of E that guarantees total convergence of P1 to x∗=0.3329 and y∗=0.6568.
(a)
V
1
(b)
V
2
(c)
V
3


## 4.2.2. Selection of the Reference

Similar results can be obtained if we choose other equilibrium points to perturb the mapPi. Table 3 summarizes the control effort requirements. Note that the critical control effort Ecr varies depending on the used reference. Coexistence phenomenon can be controlled with less effort if (x∗,y∗)=(0.8426,0.4506) or (x¯∗,y¯∗)=(0.4506,0.8426) are used to control the map.


## 4.3. Control of Period-5 Orbits in Nonlinear Coupling Maps

Now, we control two coexisting period-5 cycles of (2) when r=2.8 and ϵ=1.5. The reference is introduced every fifth iteration as is presented in(8)
x
k
+
1
y
k
+
1
=
1
E
+
1
f
x
k
+
ε
f
y
k
-
f
x
k
f
y
k
+
ε
f
x
k
-
f
y
k
+
E
E
+
1
x
¯
∗
y
¯
∗
if
k
mod
⁡
5
=
0
f
x
k
+
ε
f
y
k
-
f
x
k
f
y
k
+
ε
f
x
k
-
f
y
k
if
k
mod
⁡
5
≠
0
.According to Table2, we use (x¯∗,y¯∗)1≔(0.3526,0.7918) to control the system. Figure 11 presents BA(P¯1) of the controlled system when three different values of E are considered.Table 2
Two coexisting period-5 orbits for nonlinear coupling maps whenϵ=1.5 and r=2.8. For the sake of simplicity, the periodic orbit values are presented with 8 decimal digits.
First period-5 orbit
Second period-5 orbit
(
x
∗
,
y
∗
)
1
=
(
0.79175845,0.35258819
)
(
x
-
∗
,
y
-
∗
)
1
=
(
0.35258819,0.79175845
)
(
x
∗
,
y
∗
)
2
=
(
0.72790518,0.37290574
)
(
x
-
∗
,
y
-
∗
)
2
=
(
0.37290574,0.72790518
)
(
x
∗
,
y
∗
)
3
=
(
0.70487469,0.50446287
)
(
x
-
∗
,
y
-
∗
)
3
=
(
0.50446287,0.70487469
)
(
x
∗
,
y
∗
)
4
=
(
0.75867944,0.52373859
)
(
x
-
∗
,
y
-
∗
)
4
=
(
0.52373859,0.75867944
)
(
x
∗
,
y
∗
)
5
=
(
0.79131428,0.41974570
)
(
x
-
∗
,
y
-
∗
)
5
=
(
0.41974570,0.79131428
)Table 3
Summary of FPIC requirements in coexistence control of two period-6 cycles.
Maps
x
∗
x
-
∗
E
cr
P
1
,
P
-
1
(
0.7974,0.2023
)
(0.2023, 0.7974)
3.371
P
2
,
P
-
2
(
0.3329,0.6568
)
(0.6568, 0.3329)
1.774
P
3
,
P
-
3
(
0.7689,0.6025
)
(0.6025, 0.7689)
2.507
P
4
,
P
-
4
(
0.4991,0.7793
)
(0.7793, 0.4991)
1.492
P
5
,
P
-
5
(
0.8426,0.4506
)
(0.4506, 0.8426)
1.258
P
6
,
P
-
6
(
0.2996,0.8655
)
(0.8655, 0.2996)
4.112Figure 11
B
A
(
P
¯
1
) for E=0.25. (b) BA(P¯1) for E=0.45. (c) BA(P¯1) for E=0.65. r=2.8 and ϵ=1.5 for all cases.
(a)
(b)
(c)Selection of Minimum Control Effort. As in the previous case, the minimum control effort Ecr can be computed by studying the behavior of UPOs around the periodic orbits. Figure 12 shows an illustrative sketch of bifurcation scenario of the uncontrolled system. The unstable period-1 orbits are noted as N1 and N2. The metamorphosis of the controlled maps BA(P¯i) is influenced by unstable period-1 solutions of the system. We use points N1 and N2 to compute Ecr in each case. Table 4 presents the minimum value of parameter E to control coexistence phenomenon. All initial conditions belonging to the considered region converge to x¯∗ when E>Ecr. The lowest control effort that avoids coexistence phenomenon is reached when the map is controlled with the first point of the periodic orbit. If we induce the control using other points of the periodic orbit, then we need a greater control effort. Figure 13 shows the symmetrical behavior between BA(P1) and BA(P¯1) (compare with Figure 11). Critical control effort required in both cases has the same value.Table 4
Summary of FPIC requirements in coexistence control of two period-5 cycles.
Maps
x
∗
x
-
∗
E
cr
P
1
,
P
-
1
(
0.7918,0.3526
)
(0.3526,   0.7918)
0.6099
P
2
,
P
-
2
(
0.7279,0.3729
)
(0.3729,   0.7279)
0.6873
P
3
,
P
-
3
(
0.7049,0.5045
)
(0.5045,   0.7048)
1.2312
P
4
,
P
-
4
(
0.7588,0.5237
)
(0.5237,   0.7588 )
1.1892
P
5
,
P
-
5
(
0.7913,0.4197
)
(0.4197,   0.7913)
0.7504Figure 12
Sketch of unstable period-1 solutionsNi in coexistence zone.Figure 13
B
A
(
P
1
) for E=0.25. (b) BA(P1) for E=0.45. (c) BA(P1) for E=0.65. r=2.8 and ϵ=1.5 for all cases.
(a)
(b)
(c)


## 5. Discussion

The results presented above give rise to a new approach for coexistence control in multistable systems. The proposed methodology has significant differences with other methods reported in the literature, some of which are commented on below.(i) The Applicability of the Method. Our process is sequential, systematic, and applicable to any multistable system, even though we use low-dimensional discrete maps as illustrative examples. One of the most common methods to control multistability, which is based on slow harmonic modulation, has restrictions on its application because the attractor to be destroyed should be a focus type [23].(ii) The Requirement of Dynamic Behavior Knowledge. Fixed-point inducting control needs to know an approximate value of the periodic orbit to be stabilized, while all nonfeedback methods do not require a priori information of the system to be controlled [18]. However, this is not a major problem due to the simplicity of the required information and the multiple ways to obtain it, by either theoretical, empirical, or experimental sources.(iii) The Control Scheme, the Choice of the Attractor, and the Number of Parameters to Be Tuned. The most representative methods for multistability control can be classified into three categories:nonfeedback,feedback, andstochastic schemes [18]. We propose afeed-forward control configuration which has not been widely exploited. Particularly, the control law based on FPIC criterion has advantages with respect to the selection of the target orbit and the parameters setting. For example, when the coexistence control is realized by parameter modulation, the general form given by p=p0+pcsin⁡(2πfct+ψ0) has four parameters to be tuned [24]; then, the selection of a particular attractor depends on the right choice of these parameters. In contrast, to successfully apply FPIC control only it is needed to tune one parameter and to know the target orbit.


## 6. Conclusions

The coexistence of periodic solutions has been controlled using FPIC technique and the required control effort to obtain stabilization of the periodic orbits was computed.It has been proven that FPIC controller can be used to control BA. FPIC requires virtually no information about the system; just approximated knowledge of a periodic orbit is needed. Thus, this is a good control technique if the system is only known approximately, which is a fundamental requirement for other control methods.The applicability of the proposed tools has been illustrated in linear and nonlinear coupling logistic maps. However, the developed tools are generic and can be applied to any system with complex scenarios.



---

*Source: 313154-2015-12-03.xml*