Open AccessArticle

Fuzzy Synchronization of Chaotic Systems with Hidden Attractors

Jessica Zaqueros-Martinez

Gustavo Rodriguez-Gomez

^1,*

Esteban Tlelo-Cuautle

and

Felipe Orihuela-Espina

^1,3

Department of Computer Science, Instituto Nacional de Astrofísica, Óptica y Electrónica (INAOE), Tonantzintla, Puebla 72840, Mexico

Department of Electronics, Instituto Nacional de Astrofísica, Óptica y Electrónica (INAOE), Tonantzintla, Puebla 72840, Mexico

School of Computer Science, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK

Author to whom correspondence should be addressed.

Entropy 2023, 25(3), 495; https://doi.org/10.3390/e25030495

Submission received: 15 December 2022 / Revised: 9 February 2023 / Accepted: 13 February 2023 / Published: 13 March 2023

(This article belongs to the Section Complexity)

Download

Browse Figures

Review Reports Versions Notes

Abstract

Chaotic systems are hard to synchronize, and no general solution exists. The presence of hidden attractors makes finding a solution particularly elusive. Successful synchronization critically depends on the control strategy, which must be carefully chosen considering system features such as the presence of hidden attractors. We studied the feasibility of fuzzy control for synchronizing chaotic systems with hidden attractors and employed a special numerical integration method that takes advantage of the oscillatory characteristic of chaotic systems. We hypothesized that fuzzy synchronization and the chosen numerical integration method can successfully deal with this case of synchronization. We tested two synchronization schemes: complete synchronization, which leverages linearization, and projective synchronization, capitalizing on parallel distributed compensation (PDC). We applied the proposal to a set of known chaotic systems of integer order with hidden attractors. Our results indicated that fuzzy control strategies combined with the special numerical integration method are effective tools to synchronize chaotic systems with hidden attractors. In addition, for projective synchronization, we propose a new strategy to optimize error convergence. Furthermore, we tested and compared different Takagi–Sugeno (T–S) fuzzy models obtained by tensor product (TP) model transformation. We found an effect of the fuzzy model of the chaotic system on the synchronization performance.

Keywords:

synchronization; chaos; hidden attractor; fuzzy control; T–S fuzzy model

1. Introduction

Synchronization of chaotic systems has received substantial attention theoretically and experimentally in recent years. Such attention is justified by the potential applications in secure communication [1,2,3], the control of chemical reactions with the aim of determining the favorable conditions for practical implementation [4], the stability of the chaotic wind power system in finite time [5], the synchronization of chaotic finance systems with known and unknown systems parameters [6,7], authentication from brain signals [8], chaotic attitude synchronization and anti-synchronization for master–slave satellites under unknown moments of inertia and disturbance torques [9], regulation of glucose–insulin concentrations from a chaotic regime (an illness) to a disorder-free equilibrium [10], and the relation between a meteorological phenomenon and influenza pandemics [11].

The synchronization of dynamic systems is the problem of enforcing several systems to match and converge [12]. Chaotic systems are deterministic nonlinear dynamic systems that exhibit a prominent sensitivity to initial conditions, noise, and system parameters [13]. Such characteristics of chaotic systems make the problem of synchronizing these systems particularly challenging. In the case of chaotic systems, it is known that a chaotic solution is neither asymptotically stable, nor totally unstable, but it alternates between two or more attractors over time.

Moreover, the synchronization of chaotic systems with hidden attractors is even more challenging as they can lead to erroneous coupling behavior and instability. Besides, the system analysis can be complicated [14,15,16]; on the one hand, hidden attractors have been observed in systems without no equilibrium points, with no unstable equilibrium points, with one stable equilibrium point, or with a line equilibrium. On the other hand, the basins of attraction of the hidden attractors do not intersect with small neighborhoods of any equilibrium points (if any) and are located far from such points. To avoid or reduce the risk of a sudden jump to unwanted behavior, it is helpful to know the properties of hidden attractors.

In order to carry out the synchronizations, it is necessary to leverage numerical systems so as to approximate the solutions of the systems involved. The global error

e_{n} (t) = y_{n} - y (t_{n})

, where

y_{n}

is the approximation to the theoretical value

y (t_{n})

at time

t_{n}

, of the numerical method is upper bounded by the error made on the initial conditions, the truncation error, and the rounding-up error [17,18]. The inadequate choice of the integration step may result in the growing of the truncation error and/or the rounding error. Consequently, the system may be perturbed, yielding undesired behaviors or even making the system inoperable in the presence of hidden attractors.

The goal of this research is to afford a new methodology to carry out the synchronization of two identical chaotic systems with hidden attractors by means of fuzzy control and a special integration method that exploits the oscillating property of the chaotic systems for numerical integration.

Based on our previous expertise, we opted for complete and projective synchronization schemes. The complete synchronization employs the linear matrix inequality (LMI), whilst the projective synchronization employs parallel distributed compensation (PDC). The former proposal was applied to a set of chaotic systems with hidden attractors known in the literature.

The main contribution of this research is to provide a new methodology to synchronize chaotic systems with hidden attractors through complete and projective fuzzy controls, evidencing that those fuzzy control strategies combined with such special integration method boost robustness over other popular approaches and, hence, are more likely successful in the task of synchronizing chaotic systems with hidden attractors.

The convergence to zero of the error in the synchronization varies depending on the strategy used for the control. Projective synchronization depends on two parameters,

α

and

γ

. We further propose here a new strategy for selecting the control parameter

γ

, then we manipulated the scaling parameter

α

, as well as the initial conditions to modulate the speed of convergence according to the demands of the application. For this, we supported our findings with numerical simulations.

Another contribution is to evidence that fuzzy controls used in this work are robust in handling multistable systems with hidden attractors, and the fuzzy controls used are also robust to changes in the initial conditions. Another contribution is that the tested chaotic systems were chosen arbitrarily from the literature. Therefore, these systems are interchangeable with any other continuous integer-order chaotic system with nonlinear terms with a state variable in common. Few works in the literature emphasize the numerical methods used to simulate chaotic systems since they can generate problems such as superstability. Our proposal is to use a special fixed-step numerical method that exploits the oscillatory characteristic of chaotic systems and where the integration step is easy to choose. Furthermore, an analysis of the performance of the synchronizations using different Takagi–Sugeno (T–S) fuzzy models was conducted. The results showed an effect of the fuzzy model of the chaotic system on the synchronization performance.

2. Related Work

There are two configurations in synchronization: bidirectional or mutual coupling systems, where both systems must adapt to each other, and master–slave systems or unidirectional coupling in which one system must adapt to another [19,20]. In both configurations, the synchronization efforts can be categorized by the type of error being minimized [12], such as complete or identical synchronization, anti-synchronization, phase synchronization, and projective synchronization, among other synchronizations.

The literature reports several efforts regarding the synchronization of chaotic systems in the presence of hidden attractors. One of the most-used controls is feedback control with complete error [16,21,22,23], synchronizing with one or two state variables with the master–slave and bidirectional approach [20,24]. To use this control, it is necessary to adjust the coupling parameter to achieve synchronization. This can be done through bifurcation diagrams, but they are very computationally expensive. This control intentionally fails when synchronizing Chua’s system with hidden attractors to show that the existence of hidden attractors significantly affects synchronization [15]. Another control is sliding mode control (SMC) with complete error or anti-synchronization [25,26,27,28]. However, this control requires studying and understanding the theory, which has strong mathematical foundations. Most of the time, researchers focus on solving the problem of synchronizing systems with perturbations and unknown system parameters. Another used control is adaptive control with complete error or anti-synchronization, where they focus on the time of convergence of the error in the synchronization [29] or the problem where the system parameters are unknown [30,31,32,33,34]. Another control proposal is the fixed-time synchronization observer, which emphasizes the time of convergence of the error in the synchronization under complete error [35]. Finally, another control, back-stepping control, synchronizes two different systems with complete error by setting an extra variable [36].

Most of the controls reviewed are specific regarding the set of chaotic systems that they synchronize. For example, SMC works with perturbed systems and with unknown parameters. Other controls are interested in the analysis of the convergence time of the synchronization. In this article, we addressed the problem of synchronizing two identical chaotic systems with hidden attractors considering the complete and projective errors. To solve this problem, it is necessary to have a control, but the choice of this control is not straightforward in the presence of hidden attractors [16]. In addition, a general control was sought that guarantees the solution to a broad range of chaotic systems, avoiding ad hoc controls, which particularize the solution.

On the other hand, fuzzy control is an efficient and effective tool in nonlinear systems and has been successfully applied to chaotic synchronization before [37]. Fuzzy control theory is widely studied. Furthermore, the part that is used to perform the synchronization does not require a priori knowledge of fuzzy control, nor specific skills to use it [38]. Fuzzy control applied to the synchronization of chaotic systems is a general and unified control because it synchronizes a variety of chaotic systems, unlike controls such as linear feedback control, nonlinear feedback control, and impulsive control, which deal with one or two kinds of specified chaotic systems [37].

3. Global Error

A known result in the literature is how the global error

e_{n} (t)

of the numerical approximation of a solution of ordinary differential equations with initial values is bounded according to

| e_{n} (t) | < {A δ + B h^{p} + C h^{q}}

where

A, B, C

are constants that depend on the methods’ coefficients and the Lipschitz constant of ordinary differential equations and

δ

is the bound initial conditions [17,18].

The bounds of the global error are then given in terms of the error in the initial values

δ

, the discretization or truncation error, and the round-off error. The choice of the integration step ought to maintain a compromise between the later two. Besides, the error in the initial values in the floating point arithmetic ought to be of the same order as that of the integration method. Because of the sensitivity to the initial conditions, in the case of dynamic chaotic systems, the integration step h and the error in the representation are of utmost importance. If we choose an integration step that keeps the discretization error bounded, but which does not bound the growing of the round-off error, the later can alter the chaotic system; for example, the chaos may be lost (superstability) or make the system attracted to a different attractor, thus changing the system’s behavioral regime. Even the variable step numerical integration methods with error control exhibit difficulties when choosing the integration step to constrain the growth of the error

e (t)

, leading to the so-called computational chaos of systems that are not originally chaotic. For example, Skufca reported how the Runge–Kutta–Fehlberg method implemented with error control and a variable step (ode45 in Matlab) when applied to a system that models the coupling of two oscillators originates the computational chaos [39].

4. Special Numerical Methods

Whenever some characteristic or property of the dynamical systems is known, said property can be exploited by the numerical method. Special numerical methods are those that exploit some specific property present in the dynamical system. In such circumstances, the numerical approximation to the solution of the ordinary differential equations with some initial values can be enhanced capitalizing on such prior knowledge about the dynamical system. It follows that, since chaotic systems exhibit oscillations, it is convenient when integrating them to use a method that leverages such an oscillatory property. Therefore, in this research, we opted for using the numerical method developed by Gautschi [40]. This is a well-known method based on trigonometric polynomials that effectively exploits this oscillatory property.

Gautschi’s method has been described and analyzed in detail in [17,40], but briefly, it is an explicit method of trigonometric order

q = 1

and algebraic order

p = 2

, given by

y_{n + 2} - y_{n + 1} = h (β_{1} f_{n + 1} + β_{0} f_{n})

where

\begin{matrix} β_{0} (ν) & = & - \frac{1}{2} (1 + \frac{1}{12} ν^{2} + \frac{1}{120} ν^{4} + \dots) \\ β_{1} (ν) & = & \frac{3}{2} (1 - \frac{1}{4} ν^{2} + \frac{1}{120} ν^{4} + \dots) \end{matrix}

with

ν = ω h

, the integration step h constant, and the truncation error

\frac{5}{12} h^{3} (y^{(3)} + ω^{2} y^{(1)})

. If the frequency

ω

is unknown, it is safer to underestimate it than to overestimate it [17].

5. Fuzzy Control

Fuzzification translates crisp inputs into fuzzy values for the inference system. The fuzzy inference system applies the fuzzy rules and generates the system outputs in fuzzy form. The counterpart, defuzzification, proceeds to recover crisp values for the outputs.

5.1. Fuzzification

Here, the fuzzy inference system employs the T–S fuzzy dynamic model [41]. The T–S model is described by a set of rules: IF–THEN. These rules are capable of exactly describing nonlinear systems within some region of interest. Nonlinearities are expressed as inferred fuzzy outputs by indicating the fuzzy membership function in the premise and the associated coefficients in the consequence [42].

Consider the continuous-time chaotic system given by

\dot{x} = f (x (t))

(1)

where

x (t) = {[x_{1} (t), x_{2} (t), \dots, x_{n} (t)]}^{T} \in R^{n}

is the state vector and

f

is a nonlinear function of appropriate dimensions. The T–S fuzzy model is composed of a set of rules of the form:

\begin{matrix} Rule i : & IF z_{1} (t) is F_{i 1} and \dots and z_{p} (t) is F_{i p} \\ THEN \dot{x} = A_{i} x (t) + b_{i}, i = 1, 2, 3, \dots, I \end{matrix}

(2)

where

z_{1} (t) \dots z_{p} (t)

are the premise state variables,

F_{i, j}

(

j = 1, \dots, p

) are fuzzy sets,

A_{i}

is a constant system matrix with the appropriate size, and

b_{i} \in R^{n}

is a constant vector; it represents the bias generated by the exact fuzzy model. The number of rules I depends on the variables exhibiting the nonlinearities.

5.2. Defuzzification

Using the singleton fuzzifier, fuzzy inference by the product, and defuzzification by the weighted average, the final outputs are calculated according to

\dot{x} = \frac{\sum_{i = 1}^{I} ω_{i} (z (t)) {A_{i} x (t) + b_{i}}}{\sum_{i = 1}^{I} ω_{i} (z (t))}

(3)

where

z (t) = {[z_{1} (t) z_{2} (t) \dots z_{p} (t)]}^{T},

with

w_{i} (z (t)) = \prod_{j = 1}^{p} F_{i j} (z_{j} (t))

in which

F_{i j} (z_{j} (t))

is the degree of membership of

z_{j} (t)

F_{i j}

, with

\sum_{i = 1}^{I} ω_{i} (z (t)) > 0

, and

ω_{i} (z (t)) \geq 0

i = 1, 2, 3, \dots, I

By employing

μ_{i} (z (t)) = ω_{i} (z (t)) / \sum_{i = 1}^{I} ω_{i} (z (t))

instead of

ω_{i} (z (t))

, System (3) is rewritten as

\dot{x} = \sum_{i = 1}^{I} μ_{i} (z (t)) {A_{i} x (t) + b_{i}}

(4)

where

\sum_{i = 1}^{I} μ_{i} (z (t)) = 1

with

μ_{i} (z (t)) \geq 0

(

i = 1, 2, \dots, I

When modeling a chaotic system using a T–S fuzzy surrogate, the goal is to build a T–S model following (2), representing exactly the nonlinear system in (1). The vector function

f (x (t))

is expressed as the inferred fuzzy output

\sum_{i = 1}^{I} μ_{i} (z (t)) {A_{i} x (t) + b_{i}}

in (4). There exist several ways to implement fuzzy modeling [38,42,43]. Here, we chose the one described in [42,43]. Because each chaotic test system has nonlinear terms with a state variable in common, we only need two rules to model them in a fuzzy manner. Other chaotic systems may require a larger number of rules. For instance, in the transformed Rossler system, it is not possible to extract two linear terms from the premise variable [43], and therefore, more than two rules are required for the fuzzification of this system.

5.3. Fuzzy Modeling of Chaotic Systems

In this paper, we considered the fuzzy modeling of chaotic systems of integer order with hidden attractors, as listed in Table 1. Since the derivatives that appear in chaotic systems are of integer order, it is said that the chaotic systems are of integer order. The first column reports the name of the system and the author of the system. In the second column, the mathematical model is listed. In the third column, the parameters of the system are listed; in case these do not come from the seminal work, we cite the source from which we took them. In the last column, it is indicated if the chaotic system has equilibrium points. Below, we present the fuzzy modeling of each of the systems.

The chaotic system given by Chua [44,49] was one of the earliest to be investigated with hidden attractors [50]. This three-dimensional dynamical system is known to exhibit nonlinearities [49]. For the Chua system with hidden attractors, assume that

x_{1} \in [- d, d]

and

d > 0

. Then, the fuzzy model that exactly represents the chaotic system under

x_{1} \in [- d, d]

is given by

\begin{matrix} Rule 1 : IF x_{1} (t) is F_{1} THEN \dot{x} (t) = A_{1} x (t) + b_{1} \\ Rule 2 : IF x_{1} (t) is F_{2} THEN \dot{x} (t) = A_{2} x (t) + b_{2} \end{matrix}

(5)

where

x = {[x_{1}, x_{2}, x_{3}]}^{T}

b_{1} = b_{2} = {[0, 0, 0]}^{T}

d = 1.1

F_{1} (x_{1}) = \frac{1}{2} (1 - \frac{ϕ (x_{1})}{d})

F_{2} (x_{1}) = 1 - F_{1} (x_{1})

, and

A_{1} = (\begin{matrix} α (d - 1) & α & 0 \\ 1 & - 1 & 1 \\ 0 & - β & - γ \end{matrix}) A_{2} = (\begin{matrix} - α (d + 1) & α & 0 \\ 1 & - 1 & 1 \\ 0 & - β & - γ \end{matrix})

ϕ (x_{1}) = \{\begin{matrix} \frac{ψ (x_{1})}{x_{1}} & x_{1} \neq 0 \\ m_{0} & x_{1} = 0 \end{matrix} sat (x_{1}) = \{\begin{matrix} - 1 & if x_{1} \leq - 1 \\ x_{1} & if ∣ x_{1} ∣ < 1 \\ 1 & if x_{1} \geq 1 \end{matrix}

The final output of the Chua fuzzy system is inferred by

\dot{x} = \sum_{i = 1}^{2} F_{i} (x_{1} (t)) {A_{i} x (t) + b_{i}} .

The chaotic behavior of the Chua system around the hidden attractor is illustrated in Figure 1a.

Analogous to the Chua system, the Sprott A system assumes that

x_{2} \in [- d, d]

and

d > 0

, then the fuzzy model is

\begin{matrix} Rule 1 : IF x_{2} (t) is F_{1} THEN \dot{x} (t) = A_{1} x (t) + b_{1} \\ Rule 2 : IF x_{2} (t) is F_{2} THEN \dot{x} (t) = A_{2} x (t) + b_{2} \end{matrix}

where

x = {[x_{1}, x_{2}, x_{3}]}^{T}

b_{1} = b_{2} = {[0, 0, 2]}^{T}

d = 5.0

F_{1} (x_{2}) = \frac{1}{2} (1 + \frac{x_{2}}{d})

F_{2} (x_{2}) = \frac{1}{2} (1 - \frac{x_{2}}{d})

, and

A_{1} = (\begin{matrix} 0 & 1 & 0 \\ - 1 & 0 & d \\ 0 & - d & 0 \end{matrix}) A_{2} = (\begin{matrix} 0 & 1 & 0 \\ - 1 & 0 & - d \\ 0 & d & 0 \end{matrix})

The final output of the Sprott A fuzzy system is inferred by

\dot{x} = \sum_{i = 1}^{2} F_{i} (x_{2} (t)) {A_{i} x (t) + b_{i}}

and the phase space of the Sprott A system is shown in Figure 1b.

For the NE₆ system, assume that

x_{3} \in [- d, d]

with

d > 0

, then the fuzzy model is

\begin{matrix} Rule 1 : IF x_{3} (t) is F_{1} THEN \dot{x} (t) = A_{1} x (t) + b_{1} \\ Rule 2 : IF x_{3} (t) is F_{2} THEN \dot{x} (t) = A_{2} x (t) + b_{2} \end{matrix}

where

x = {[x_{1}, x_{2}, x_{3}]}^{T}

b_{1} = b_{2} = {[0, 0, - 0.75]}^{T}

d = 3.0

F_{1} (x_{3}) = \frac{1}{2} (1 + \frac{x_{3}}{d})

F_{2} (x_{3}) = \frac{1}{2} (1 - \frac{x_{3}}{d})

, and

A_{1} = (\begin{matrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ - d & - 1 - d & 0 \end{matrix}) A_{2} = (\begin{matrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ d & - 1 + d & 0 \end{matrix})

The final output of the NE₆ fuzzy system is inferred by

\dot{x} = \sum_{i = 1}^{2} F_{i} (x_{3} (t)) {A_{i} x (t) + b_{i}}

and the phase space of this system is shown in Figure 1c.

The LE₄ system can be fuzzy modeled as follows. Assume that

x_{2} \in [- d, d]

with

d > 0

, then the fuzzy model is

\begin{matrix} Rule 1 : IF x_{2} (t) is F_{1} THEN \dot{x} (t) = A_{1} x (t) + b_{1} \\ Rule 2 : IF x_{2} (t) is F_{2} THEN \dot{x} (t) = A_{2} x (t) + b_{2} \end{matrix}

where

x = {[x_{1}, x_{2}, x_{3}]}^{T}

b_{1} = b_{2} = {[0, 0, 0]}^{T}

d = 0.8

F_{1} (x_{2}) = \frac{1}{2} (1 + \frac{x_{2}}{d})

F_{2} (x_{2}) = \frac{1}{2} (1 - \frac{x_{2}}{d})

, and

A_{1} = (\begin{matrix} 0 & 1 & 0 \\ - 1 & 0 & d \\ - d (a + b) & 0 & - d \end{matrix}) A_{2} = (\begin{matrix} 0 & 1 & 0 \\ - 1 & 0 & - d \\ d (a + b) & 0 & d \end{matrix})

The final output of the LE₄ fuzzy system is inferred by

\dot{x} = \sum_{i = 1}^{2} F_{i} (x_{2} (t)) {A_{i} x (t) + b_{i}}

and its phase space is shown in Figure 1d.

6. Fuzzy Synchronization

The general master–slave feedback control utilized to synchronize chaotic systems is summarized in Figure 2. In this paper, we used fuzzy control due to the characteristics that we explained in Section 1. The two fuzzy controls used to synchronize chaotic systems are described below.

6.1. Complete Fuzzy Synchronization

Let us consider the fuzzy representation of the chaotic system (2) as the master system and the system (6) as the slave system:

\begin{matrix} Rule 1 : IF y_{1} (t) is F_{1} THEN \dot{y} (t) = A_{1} y (t) + b_{1} + B u (t) \\ Rule 2 : IF y_{1} (t) is F_{2} THEN \dot{y} (t) = A_{2} y (t) + b_{2} + B u (t) \end{matrix}

(6)

where

y = {[y_{1}, y_{2}, y_{3}]}^{T}

are state variables and

A_{1}

A_{2}

b_{1}

b_{2}

F_{1}

F_{2}

, and d have the same values as the master system (2).

B

is an input matrix, and

u (t)

is the control.

Defuzzification is given by

\begin{matrix} \dot{y} = \sum_{i = 1}^{2} μ_{i} (y_{1} (t)) {A_{i} y (t) + b_{i}} + B u (t), μ_{i} (y_{1} (t)) = F_{i} (y_{1} (t)) . \end{matrix}

The error is

e (t) = y (t) - x (t)

. The error can be estimated from

\begin{matrix} \dot{e} (t) = \sum_{i = 1}^{2} μ_{i} (y_{1} (t)) {A_{i} y (t) + b_{i}} - \sum_{i = 1}^{2} μ_{i} (x_{1} (t)) {A_{i} x (t) + b_{i}} + B u (t) . \end{matrix}

(7)

The goal of the synchronization is to design the fuzzy control:

\begin{matrix} u (t) = - \sum_{i = 1}^{2} μ_{i} (y_{1} (t)) {C_{i} y (t) + b_{i}} + \sum_{i = 1}^{2} μ_{i} (x_{1} (t)) {C_{i} x (t) + b_{i}} \end{matrix}

(8)

such that

∥e (t)∥ \to 0

when

t \to \infty

. Substituting the fuzzy control (8) in the error (7), one obtains

\begin{matrix} \dot{e} (t) = \sum_{i = 1}^{2} μ_{i} (y_{1} (t)) {(A_{i} - B C_{i}) y (t)} - \sum_{i = 1}^{2} μ_{i} (x_{1} (t)) {(A_{i} - B C_{i}) x (t)} \end{matrix}

The gain matrices

C_{i}

are calculated during the design process. With the idea of linearization given in [38,42], if there exists gain matrices such that

\begin{matrix} {(A_{1} - B C_{1}) - (A_{2} - B C_{2})}^{T} \times {(A_{1} - B C_{1}) - (A_{2} - B C_{2})} = 0 \end{matrix}

then the total error of the system becomes linear when

\dot{e} (t) = G e (t)

according to the fuzzy control (8), where

G = A_{1} - B C_{1} = A_{2} - B C_{2}

. Moreover, if

G < 0

, the the error is asymptotically stable. This convergence was proven in [42].

6.2. Projective Fuzzy Synchronization

In the projective synchronization, the drive system is scaled to a scalar factor

α

. Let the master system (1) and the slave system be

\dot{y} = g (y (t), u (x (t), y (t)))

. If there exist a constant

α \neq 0

such that

{lim}_{t \to \infty} ∥y (t) - α x (t)∥ = 0

, then the synchronization is referred to as projective, and

α

is called the scaling factor.

In particular, for the case study at hand, our master system is given by its T–S fuzzy modeling (9) and the slave system (11).

\begin{matrix} Rule 1 : IF x_{1} (t) is F_{1} & THEN \dot{x} (t) = A_{1} x (t) \\ Rule 2 : IF x_{1} (t) is F_{2} & THEN \dot{x} (t) = A_{2} x (t) \end{matrix}

(9)

The final output of the fuzzy system (9) is inferred from

\dot{x} = \sum_{i = 1}^{2} μ_{i} (x_{1} (t)) A_{i} x (t)

(10)

where

μ

satisfies the same conditions as (4).

\begin{matrix} Rule 1 : IF x_{1} (t) is F_{1} & THEN \dot{y} (t) = A_{1} y (t) + u (t) \\ Rule 2 : IF x_{1} (t) is F_{2} & THEN \dot{y} (t) = A_{2} y (t) + u (t) \end{matrix}

(11)

The final output of the slave fuzzy system (11) is given by

\dot{y} = \sum_{i = 1}^{2} μ_{i} (x_{1} (t)) {A_{i} y (t) + u (t)}

(12)

Let the error

e (t) = y (t) - α x (t)

. Substituting the final outputs of the fuzzy systems (10) and (12), we obtain the derivative of the error:

\dot{e} (t) = \sum_{i = 1}^{2} μ_{i} (x_{1} (t)) {A_{i} e (t) + u (t)}

(13)

By means of PDC, which creates a controller fuzzy system with rule premises identical to those of the plant fuzzy system [51], the next fuzzy control is obtained:

\begin{matrix} Rule 1 : IF x_{1} (t) is F_{1} THEN & u (t) = - γ A_{1} [y (t) - α x (t)] - γ [y (t) - α x (t)] \\ Rule 2 : IF x_{1} (t) is F_{2} THEN & u (t) = - γ A_{2} [y (t) - α x (t)] - γ [y (t) - α x (t)] \end{matrix}

where

γ

is a control parameter.

The final output of the fuzzy control is given by:

\begin{matrix} u (t) = \sum_{i = 1}^{2} μ_{i} (x_{1} (t)) {- γ A_{i} [y (t) - α x (t)] - γ [y (t) - α x (t)]} \end{matrix}

(14)

Substituting the fuzzy control (14) in the error (13), one obtains

\dot{e} (t) = \sum_{i = 1}^{2} μ_{i} (x_{1} (t)) [(1 - γ) A_{i} - γ I] e (t)

(15)

where

I

is the identity matrix.

A condition for the convergence of the error (15) is that, if there exists a positive definite symmetric constant matrix P and constant

c > 0

, such that

[(1 - γ) A_{i} - γ I)]^{T} P + P [(1 - γ) A_{i} - γ I)] \leq - c I

i = 1, 2

, then the error (15) is globally exponentially stable.

Another convergence condition for the error (15) concerns the characteristic values

λ_{j}

j = 1, \dots, n

of the symmetric matrices

(1 - γ) (A_{i}^{T} + A_{i}) - 2 γ I

i = 1, 2

. If every

λ_{j}

has a strictly negative real part, then the error (15) is globally exponentially stable, which implies that (9) and (11) can asymptotically achieve projective synchronization. For the proof of these conditions, the reader is referred to the original work [52].

7. Projective Synchronizations with an Emphasis on Error Convergence

For projective synchronization, it is possible to find the value of the control parameter

γ

analytically so that the convergence of the synchronization error to zero can be made faster or slower according to the application. Through matrices

(1 - γ) (A_{i}^{T} + A_{i}) - 2 γ I

i = 1, 2

, the characteristic values are found as a function of

γ

, and its range of variation is determined. Parameter

γ

is selected, ensuring the characteristic values have a strictly negative real part, and so, the second condition stated at the end of Section 6.2 is met. By means of the standard deviation of the characteristic values

σ (λ)

λ = (λ_{1}, \dots, λ_{n})

, found in the valid range of

γ

, the set of characteristic values that gives us the sought convergence (fast or slow). When

σ (λ) = 0

, the set of associated characteristic values is very similar and not small in magnitude (due to the spectrum of variation of the values), and the convergence error quickly goes to zero. Then, determining the

λ^{*} = arg min σ (λ) = 0

, in the range of variation of

γ

, the value of

λ^{*}

is found so that the error converges as fast as possible. Otherwise,

λ^{*} = arg max σ (λ)

is found.

The aforementioned procedure was carried out using known strategies to determine the maxima and minima of a scalar function with a vector argument. However, when the expressions as a function of

γ

are polynomials of a degree greater than 1, finding

λ^{*}

is laborious in practice. In these cases (such as the ones discussed in this article), it is recommended to obtain the plot of

γ

vs.

σ (λ (γ))

to determine the desired

λ^{*}

and the range of variability of

γ

8. Simulation Results

In this section, we present the results of the numerical simulations of the chaotic systems with hidden attractors in Table 1 using the proposed fuzzy control approaches. We used the complete fuzzy and projective fuzzy synchronizations described in Section 6.1 and Section 6.2 to synchronize two identical chaotic systems with hidden attractors and different initial conditions. In the case of projective synchronization, we applied the strategy described in Section 7.

We further studied the rate of convergence by both control strategies applied to Chua’s system across a new set of initial conditions. The different initial conditions that we generated followed the methodology proposed by Danca [53]. With projective synchronization, we fixed the control parameter

γ

to the value that affords the fastest convergence, and we varied the scaling parameter and the different initial conditions. In each simulation, we looked for the first iteration from which each synchronization error was less than

2 %

for all subsequent iterations. With complete synchronization, we also varied the initial conditions and looked for the iterations that satisfy each synchronization error to be less than

2 %

. Finally, we compared the iterations found in order to analyze the convergence of the errors through them.

In addition, we analyzed different T–S fuzzy models and we found an effect of the fuzzy model of the tested chaotic system on the synchronization performance.

All the simulations were carried out on a personal computer with machine precision

2.2204 \times 10^{- 16}

in Fortran with the GNU Fortran 10.3.0 compiler for Windows. No compiler options and optimizations were selected. The graphs were made with matplotlib 3.5.1.

In all numerical simulations, the differential equations were solved by using the Gautschi integration method [40] based on trigonometric polynomials. Unlike other fixed-step integration methods such as Euler or Runge–Kutta, the Gautschi method is a special integration method that takes the oscillations of chaotic systems into account. Gautschi’s method has been shown to exhibit accuracy closest to the variable-step method with control error ode45 and requires fewer evaluations of the chaotic system than Runge–Kutta 4 (RK4) [54,55].

In order to employ the Gautschi integration method, we ought to establish the system frequency and the integration step. If the system frequency is unknown, it is possible and advisable to underestimate it [40]. In the simulations carried out here, we used an underestimated frequency equal to 1 and an integration step of

1 / 128

. With these values, we obtained the expected attractors, which we confirmed with the phase spaces of the test systems reported in the source articles. The simulation time was 150 time units for the systems Chua, NE₆, and LE₄ and 250 time units the Sprott A system. Despite their popularity, time units may be ambiguous, and we, therefore, report our simulations in terms of iterations instead. The total number of iterations was calculated as (final time—initial time)/integration step. In our case, 150 time units correspond to 19,200 iterations and 250 time units correspond to 32,000 iterations.

8.1. Results of Fuzzy Complete Synchronization

We synchronized each test system with another identical system, but with different initial conditions. First, we used the fuzzy models of each chaotic system with hidden attractors given in Section 5.3. Then, we applied the fuzzy control presented in Section 6.1.

The gain matrices

C_{1}

and

C_{2}

of each synchronization were obtained by solving the LMIs using the Matlab R2021b LMI toolbox (The MathWorks, Inc. Natick, MA, USA). Choose the input matrix

B

as the identity matrix.

Upon every synchronization, the partial errors converged to zero as the simulation progressed, that is all synchronizations were successful. The proposed control is robust and capable of synchronizing chaotic systems with hidden attractors.

8.1.1. Chua

Two Chua systems with hidden attractors were synchronized under error-complete and initial conditions

x_{0} = {[- 3.7727, - 1.3511, 4.6657]}^{T}

for the master system and

y_{0} = {[0.01, 0, 0]}^{T}

for the slave system [15]. In order to emphasize the difficulty of selecting a control to synchronize systems with hidden attractors, these systems were considered with the given initial conditions and feedback control to obtain synchronization failure [15].

To carry out the synchronization the gain matrices for the Chua chaotic system of integer order with hidden attractors with the chosen initial conditions were calculated. The matrices (16) were obtained. With these matrices, the error was asymptotically stable.

C_{1} = (\begin{matrix} 0.3364 & 1.1820 & 0 \\ 1.1820 & - 0.1250 & - 1.3841 \\ 0 & - 1.3841 & 0.1237 \end{matrix}), C_{2} = (\begin{matrix} - 4.3145 & 1.1820 & 0 \\ 1.1820 & - 0.1250 & - 1.3841 \\ 0 & - 1.3841 & 0.1237 \end{matrix})

(16)

Substituting the gain matrices in the fuzzy control (8) and using the Gautschi integration method, we synchronized the two Chua chaotic systems of integer order with hidden attractors. The phase space is shown in Figure 3a. It can be appreciated that the slave system (dotted magenta line) successfully followed the master system (solid green line).

Figure 3b presents the evolution of the error in the three state variables in the full simulation. Leveraging the theoretical guarantees, since

∥e (t)∥ \to 0

when

t \to \infty

, it was concluded that the systems were successfully synchronized using the proposed fuzzy control plotted in Figure 3c. Further, all synchronization errors were lower than

2 %

from Iteration 12,301. This case asserts the feasibility of fuzzy control for synchronizing chaotic systems in the presence of hidden attractors. This fuzzy control overcomes the difficulties presented by chaotic systems with hidden attractors such as unexpected behaviors, sudden changes, or even system inoperability.

8.1.2. Sprott A

Analogous to the previous synchronizations over the Chua system, we proceeded to synchronize the Sprott A system. The initial conditions used were

x_{0} = {[0, 5, 0]}^{T}

for the master system [56] and

y_{0} = {[1.1, 6.1, - 0.1]}^{T}

for the slave system. The calculated gain matrices for this system are

\begin{matrix} C_{1} = C_{2} & = (\begin{matrix} 0.5 & 0 & 0 \\ 0 & 0.5 & 0 \\ 0 & 0 & 0.5 \end{matrix}) \end{matrix}

In Figure 3d, the phase space where the slave system (dotted line) mimics the master system (solid line) is shown. Figure 3e plots the three synchronization errors, all of which converged to 0. Furthermore, from Iteration 2792, all were less than

2 %

. Finally, Figure 3f depicts the behavior of the control in the three state variables. We thus concluded that the synchronization was successful.

8.1.3. NE₆

In the complete fuzzy synchronization of the NE₆ system, we used as the initial conditions

x_{0} = {[0, 3, - 0.1]}^{T}

for the master system [47] and

y_{0} = {[2.4, 1.4, - 0.6]}^{T}

for the slave system. The gain matrices obtained with this system are

C_{1} = (\begin{matrix} 0.5 & 0.5 & - 1.5 \\ 0.5 & 0.5 & - 1.5 \\ - 1.5 & - 1.5 & 0.5 \end{matrix}), C_{2} = (\begin{matrix} 0.5 & 0.5 & 1.5 \\ 0.5 & 0.5 & 1.5 \\ 1.5 & 1.5 & 0.5 \end{matrix})

Figure 3g shows the phase space, in which the slave system (dotted line) mimics the master system (solid line). Figure 3h plots the three synchronization errors, each converging to zero. From Iteration 1885, the three errors were less than

2 %

. In Figure 3i, the behavior of the fuzzy control used is presented. The controls

u_{1}

and

u_{2}

have the same behavior; hence, the graphs overlap. For all of the above, this synchronization was also successful.

8.1.4. LE₄

The initial conditions used to synchronize this system were

x_{0} = {[0.2, 0.7, 0]}^{T}

for the master system [48] and

y_{0} = {[- 0.1, 0.3, 0.1]}^{T}

for the slave system. The gain matrices obtained are

C_{1} = (\begin{matrix} 0.5 & 0 & - 1.84 \\ 0 & 0.5 & 0.4 \\ - 1.84 & 0.4 & - 0.3 \end{matrix}), C_{2} = (\begin{matrix} 0.5 & 0 & 1.84 \\ 0 & 0.5 & - 0.4 \\ 1.84 & - 0.4 & 1.3 \end{matrix})

Figure 3j presents the phase space. The correct synchronization of the slave system (dotted line) with the master system (solid line) can be appreciated. Figure 3k shows the convergence of the errors to 0. From Iteration 3383, all errors were less than

2 %

. In Figure 3l, the behavior of the fuzzy control used is graphed. In this sense, synchronization was successful.

8.2. Results of Fuzzy Projective Synchronization

In this section, we present the results of the projective fuzzy synchronization for the test systems in Table 1. We considered the fuzzy models of the systems reported in Section 5.3. We applied the control in Section 6.2.

Fuzzy projective control (presented in Section 6.2) has the constraint that none of the chaotic systems involved, neither the master nor the slave, can have constant terms. Therefore, the only systems that can be synchronized with this control of those in Table 1 are Chua and LE₄.

Here, we present the results of applying the strategy given in Section 7 to choose

γ

and that allows modulating the speed of convergence.

8.2.1. Chua System

With those parameters presented in Section 5.3, we calculated the symmetric matrices

(1 - γ) (A_{i}^{T} + A_{i}) - 2 γ I

i = 1, 2

, which are

\begin{matrix} (\begin{matrix} 1.69124 (1 - γ) - 2 γ & 9.4562 (1 - γ) & 0 \\ 9.4562 (1 - γ) & - 2 (1 - γ) - 2 γ & - 11.0732 (1 - γ) \\ 0 & - 11.0732 (1 - γ) & 0.0104 (1 - γ) - 2 γ \end{matrix}) \end{matrix}

\begin{matrix} (\begin{matrix} - 35.516 (1 - γ) - 2 γ & 9.4562 (1 - γ) & 0 \\ 9.4562 (1 - γ) & - 2 (1 - γ) - 2 γ & - 11.0732 (1 - γ) \\ 0 & - 11.0732 (1 - γ) & - 0.0104 (1 - γ) - 2 γ \end{matrix}) \end{matrix}

Subsequently, we looked for the characteristic values of each matrix to determine the range of values in which

γ

must lie to satisfy that the symmetric matrices have strictly negative characteristic values. The equation to find

γ

is nonlinear. Clearing

γ

analytically can be cumbersome. Hence, instead, we opted to plot the standard deviation of the characteristic values by varying the parameter

γ

in the interval

[0.88, 1]

in Figure 4a. Such an interval can be found with trial and error with educated guesses; the sought

γ

values are those complying with the required convergence conditions.

Following the selection of the interval of the

γ

values that ensure convergence, we now present the synchronization under the lower and upper values of the function

σ (λ (γ))

within the interval where the convergence requirements are met:

γ = 0.88

and

γ = 1.0

, respectively, in the case of Chua. As the process involved some trial and error, we further show the synchronization for nearby values:

γ = 0.89

and

γ = 0.99

. The same initial conditions used for the complete synchronization were used here:

x_{0} = {[- 3.7727, - 1.3511, 4.6657]}^{T}

for the master system and

y_{0} = {[0.01, 0, 0]}^{T}

for the slave system [15]. The scaling factor used was

α = 0.5

. The errors of each state variable of the synchronizations are plotted in Figure 5a–c. The error in the synchronizations with

γ = 0.99

and

γ = 1

converged to zero in fewer iterations than with

γ = 0.88

and

γ = 0.89

in each error. It can be appreciated how the value of

γ

affected the speed of convergence: faster with

γ = 0.99

and

γ = 1

and slower with

γ = 0.88

and

γ = 0.89

, with implications for the domain application.

8.2.2. LE₄

When synchronizing this system projectively, we used the same initial conditions as for complete synchronization:

x_{0} = {[0.2, 0.7, 0]}^{T}

for the master system [48] and

y_{0} = {[- 0.1, 0.3, 0.1]}^{T}

for the slave system. We chose the scaling factor

α = 0.5

The symmetric matrices of this system are

\begin{matrix} (\begin{matrix} - 2 γ & 0 & - 3.68 (1 - γ) \\ 0 & - 2 γ & 0.8 (1 - γ) \\ - 3.68 (1 - γ) & 0.8 (1 - γ) & - 1.6 (1 - γ) - 2 γ \end{matrix}) \end{matrix}

\begin{matrix} (\begin{matrix} - 2 γ & 0 & 3.68 (1 - γ) \\ 0 & - 2 γ & - 0.8 (1 - γ) \\ 3.68 (1 - γ) & - 0.8 (1 - γ) & 1.6 (1 - γ) - 2 γ \end{matrix}) \end{matrix}

Using again the strategy described in Section 7, when calculating the characteristic values of the matrices, we also found polynomials of the third degree, as happened with the Chua system. The plots

γ

vs.

σ (λ (γ))

for

γ \in [0.7, 1.0]

are shown in Figure 4b. Analogously, as we did with Chua, this interval was found by trial and error following educated guesses. The lower and upper values of the function where the convergence requirements were met occurred at

γ = 0.77

and

γ = 1.0

, respectively. Again, because of the trial and error process, we also synchronized with nearby values:

γ = 0.78

and

γ = 0.99

Figure 5d–f plot the error for each state variable, and again, we show how the value of

γ

affects the speed of convergence;: faster with

γ = 0.99

and

γ = 1.0

and slower with

γ = 0.70

and

γ = 0.71

, in agreement with the strategy in Section 7.

8.3. A New Set of Initial Conditions

In this section, we analyze numerically the error convergence of both the complete and projective synchronizations of the Chua system with different initial conditions, and further, we investigate whether the scaling factor has an impact on the projective synchronizations.

To design further numerical simulations, like [53], first, we found the equilibrium points of the chaotic system; this system has three: two unstable

[6.58831, 0.0028364, - 6.58547]

[- 6.58831, - 0.0028364, 6.58547]

and one stable

[0, 0, 0]

, plotted in Figure 1a. Then, we verified numerically the solutions that start from neighborhoods of the unstable equilibrium points, which are attracted to the point of stable equilibrium or tend to infinity. Finally, we visualize the hidden attractors by trial and error by choosing the starting points for the numerical simulation outside the neighborhoods of the unstable equilibrium points. The list of initial conditions used is listed in Appendix A. These simulations were also carried out using Gautschi integration with an underestimated frequency of 1 and an integration step of 1/128. Again, the simulated time spanned 150 time units, equivalent to 19,200 iterations per simulation.

Afterwards, we conducted both the complete and the projective synchronizations. For the case of complete synchronization, we also used the gain matrices given in (16). For the projective synchronization, we chose to simulate scaling factors

α

equal to

0.5

, 1, 2, 3, 4, and 8. The control parameter was fixed to

γ = 0.99

in all projective simulations, which gave the fastest convergence, as estimated with the strategy given in Section 7.

In principle, the complete synchronization and the projective synchronization with

α = 1

are equivalent. However, the design of the respective control differed. We emphasize here such differences in our simulations.

All simulations were successful. Of course, some converged faster than others. To objectively establish this, we looked for the first iteration from which each synchronization error was less than

2 %

for each of the 30 newly found initial conditions. Figure 6 shows the box plot. Comparing the complete and projective synchronizations, there was a difference of 10,000 iterations between them, approximately, representing ∼50% of all 19,200 iterations.

A paired-samples right-tailed t-test was conducted to compare the results between the complete and the projective

α = 1

synchronizations. There was a significant difference between the scores for the complete (mean (

μ

) ± standard deviation (

σ

μ \pm σ = 13, 573 \pm 2032.9

) and projective (

μ \pm σ = 1641.2 \pm 256.85

) conditions;

t (29) = 32.5046, p < 0.001

We conducted a repeated-measures ANOVA test (Matlab’s R2021b ranova) to analyze the effect of the scaling factor

α

on the six tested projective synchronizations

α

= 0.5, 1, 2, 3, 4, 8. We obtained

F (5, 140) = 22.05

p < 0.001

. The post hoc comparison was conducted using a multiple comparison test (Matlab’s R2021b multcompare), where a Dunn and Sidak’s adjusted alpha revealed significant differences between

α = 0.5

, and the remaining

α

p < 0.05

, as well as

α = 8

and the others

α

p < 0.05

. The anterior can be appreciated in Figure 6.

8.4. Alternative T–S Fuzzy Models

In the synchronizations presented so far, we used the T–S fuzzy models provided in Section 5.3. However, for a given dynamical system, the fuzzy model is not unique [42,43]. A proposal to systematically obtain other T–S fuzzy models is the tensor product (TP) model transformation [57,58,59]. In order to establish the effect of the fuzzy model on the performance of the synchronization, using TP model transformation [58], we obtained two further variants of the T–S fuzzy models from each test system and compared their performances. These fuzzy models are the sum normalization non-negativeness (SNNN) type and the close-to-normality (CNO) type. In this article, we call the T–S fuzzy models presented in Section 5.3 Base, the SNNN type SNNN, and the CNO type CNO.

All synchronizations with the different T–S fuzzy models were successful with the errors converging to zero. To objectively compare the simulations, we used the synchronization error and looked for the first iteration that satisfied that each error was less than 2% for all subsequent iterations.

The results of these simulations are reported in Table 2. The results showed that, with the alternative T–S fuzzy models in complete fuzzy synchronization, the performances of the simulations improved or were equal. With the Chua system, an improvement of 79.7% was obtained. For the Sprott A system, it was 4.6%. With the NE₆ and LE₄ systems, the performances of Base and the alternatives SNNN and CNO matched. In the case of projective fuzzy synchronization, this was only carried out with the Chua and LE₄ systems due to the constraints of the fuzzy projective control, as explained in Section 8.2. With the Chua system, an improvement of 23.2% was obtained, while with the LE₄ system, 8.5% more iterations were required with the alternative SNNN, and the results of the alternative CNO and Base coincided.

9. Discussion

In this research, multistable chaotic systems [60] with hidden attractors were considered. When the hidden attractors of the systems are unknown, unexpected behaviors could be more frequent. The challenge of synchronizing them can be blamed on the fact that, currently, there is no standard strategy for the identification of all possible attractors, thus escaping unwanted behaviors. Determining the basins of attraction is a daunting task and highly time-consuming in most cases. The fuzzy controls used in this work evidenced robustness in handling multistable systems with hidden attractors. Specifically, we studied systems without equilibrium points (Sprott A and NE₆), with a line of equilibrium points (LE₄) and with a stable equilibrium point (Chua). All the synchronization efforts were successful, even in the case of the Chua system, where the synchronization through feedback control must take special care in the selection of its coupling parameter to avoid leaving the hidden attractor [15].

The tested chaotic systems were chosen arbitrarily from the literature. These systems are interchangeable with any other continuous integer-order chaotic system with nonlinear terms with a state variable in common. Chaotic maps and continuous systems of fractional order were excluded in this study. Chaotic systems with delay or with unknown parameters and/or initial conditions were not considered either. These latter would require a different design for the fuzzy controls used.

Regarding the initial conditions, 30 new initial conditions for the Chua system were found and presented in Section 8.3 and Appendix A. The corresponding synchronizations were performed, evidencing that the fuzzy controls used are also robust to changes in the initial conditions.

In this work, the Gautschi special method of numerical integration was used. This method exploits the oscillatory characteristic of chaotic systems. This is a fixed-step method that can use integration steps up to

10^{- 1}

and resources almost comparable to Euler’s method, providing competitive results to a fixed-step method with error control [54,55]. Few works in the literature pay attention to this fundamental aspect in numerical simulations of chaotic systems, but even variable-step methods with error control can present problems [39].

In traditional controls, pole assignment is a standard strategy to maximize the convergence rate (of the control). However, in fuzzy projective synchronization, this strategy is unfeasible because the equations of the fuzzy control (14) and the dynamics of the error (15) do not allow moving the poles in a simple manner. Instead, the values of the gamma parameter that satisfy the convergence conditions are sought. Therefore, a new strategy was proposed to optimize the convergence of the error in the fuzzy projective synchronizations through the analysis of the control parameters statistics.

Every dynamical system has infinitely many fuzzy models that can approximate it. Considering this, we also showed that the synchronization performance was influenced by the choice of this model.

10. Conclusions and Future Work

We studied the feasibility of synchronizing chaotic dynamic systems with hidden attractors using fuzzy control. We exemplified the success of such an approach with a case study on the Chua system, a case where a previous attempt using a different control strategy was difficult. We did not intend to send the message that fuzzy control is somehow superior to other forms of controls, but instead, we just intended to highlight that, under some circumstances, this strategy might present specific advantages. For instance, to apply the fuzzy controls presented here, it suffices to fuzzify the system (analogously to Section 5.1 and Section 5.2) and then continue with Section 6 to calculate the gain matrices or the control parameter. Importantly, the control design in this manner does not depend on the initial conditions. This is in contrast to the effort by [16] on the coupling variant of the problem for this same system.

We showed the flexibility of the approach by designing controls over two different synchronization errors; complete and projective. We stated the conditions that provide theoretical guarantees for the convergence to the solution. These conditions were previously reported in the literature. Here, we went one step further by empirically showing the practical consequences.

For projective synchronization, we proposed a new strategy for error convergence optimization. This strategy capitalizes on the theoretical conditions of the control through the parameter

γ

. Consequently, the proposed strategy can be applied to any chaotic system that is synchronized with projective fuzzy control.

We focused on four case studies. However, when reviewing the literature, we did not find theoretical limitations for this approach to be applicable to other cases of chaotic systems with hidden attractors. Moreover, for our case studies, the results behaved as predicted by theory. Hence, we believe it may be possible to extend our fuzzy control strategy to other chaotic systems.

Author Contributions

Conceptualization, J.Z.-M., G.R.-G., E.T.-C. and F.O.-E.; formal analysis, G.R.-G. and F.O.-E.; investigation, J.Z.-M.; methodology, J.Z.-M., G.R.-G., E.T.-C. and F.O.-E.; software, J.Z.-M. and G.R.-G.; supervision, G.R.-G., E.T.-C. and F.O.-E.; validation, G.R.-G.; writing—original draft, J.Z.-M.; writing—review and editing, G.R.-G., E.T.-C. and F.O.-E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Mexican Research Council: Consejo Nacional de Ciencia y Tecnología (CONACYT), Grant Number 731641, and the APC was partially funded by the University of Birmingham (C127.10028. 65661).

Data Availability Statement

A free Fortran implementation of the synchronizations can be downloaded at https://github.com/Jessica-ZM/Fuzzy_Synchronizations#fuzzy_synchronizations.

Acknowledgments

We would like to thank the Reviewers (blinded to us) for providing comments and suggestions, which improved this paper, and we especially express our gratitude to one of the Reviewers for the suggestion to test (and compare) different T–S fuzzy models.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

We found a new set of boundary conditions following the methodology suggested in [53]. These are listed in Table A1. These were used for the results reported in Section 8.3. Figure A1 presents the phase spaces of the first and eighth master systems in Table A1.

Figure A1. Some hidden attractors of master systems with new initial conditions. (a) First initial condition and (b) eighth initial condition of the master system.

Table A1. Set of 30 newly found initial conditions.

	$x_{1}$	$x_{2}$	$x_{3}$	$y_{1}$	$y_{2}$	$y_{3}$
1	4.5000	1.0000	−5.0000	0.0100	0.0100	0.0100
2	4.4000	0.9000	−4.9000	0.0100	0.0100	0.0100
3	4.3000	0.8000	−4.8000	0.0100	0.0100	0.0100
4	4.2000	0.7000	−4.7000	0.0100	0.0100	0.0100
5	4.1000	0.6000	−4.6000	0.0100	0.0100	0.0100
6	4.0000	0.5000	−4.5000	0.0100	0.0100	0.0100
7	3.9000	0.4000	−4.4000	0.0100	0.0100	0.0100
8	3.8000	0.3000	−4.3000	0.0100	0.0100	0.0100
9	3.7727	1.3511	−4.6657	0.0100	0.0100	0.0100
10	3.5000	0.1500	−6.1000	0.0100	0.0100	0.0100
11	3.4000	0.2500	−6.5000	0.0100	0.0100	0.0100
12	3.2000	0.3500	−6.3000	0.0100	0.0100	0.0100
13	3.1800	0.2000	−6.3000	0.0100	0.0100	0.0100
14	3.1500	0.2500	−6.5000	0.0100	0.0100	0.0100
15	3.1300	0.0500	−6.6600	0.0100	0.0100	0.0100
16	3.1000	0.3500	−6.4000	0.0100	0.0100	0.0100
17	3.0700	0.1500	−6.0000	0.0100	0.0100	0.0100
18	−3.0700	−0.1500	6.0000	0.0100	0.0100	0.0100
19	−3.1300	−0.0500	6.6600	0.0100	0.0100	0.0100
20	−3.1500	−0.2500	6.5000	0.0100	0.0100	0.0100
21	−3.1800	−0.2000	6.3000	0.0100	0.0100	0.0100
22	−3.2000	−0.3500	6.3000	0.0100	0.0100	0.0100
23	−3.3000	−0.0700	7.0000	0.0100	0.0100	0.0100
24	−3.4000	−0.2500	6.5000	0.0100	0.0100	0.0100
25	−3.5000	−0.1500	6.1000	0.0100	0.0100	0.0100
26	−3.6000	−0.9000	5.2000	0.0100	0.0100	0.0100
27	−3.7727	−1.3511	4.6657	4.0000	0.5000	−4.5000
28	−3.9000	−0.4000	4.4000	3.1000	0.3500	−6.4000
29	−4.0000	−0.5000	4.5000	3.3000	0.0700	−7.0000
30	−4.5000	−1.0000	5.0000	4.2000	0.7000	−4.7000

References

Ouannas, A.; Bendoukha, S.; Volos, C.; Boumaza, N.; Karouma, A. Synchronization of Fractional Hyperchaotic Rabinovich Systems via Linear and Nonlinear Control with an Application to Secure Communications. Int. J. Control. Autom. Syst. 2019, 17, 2211–2219. [Google Scholar] [CrossRef]
Khalaf, A.J.M.; Abdolmohammadi, H.R.; Ahmadi, A.; Moysis, L.; Volos, C.; Hussain, I. Extreme multi-stability analysis of a novel 5D chaotic system with hidden attractors, line equilibrium, permutation entropy and its secure communication scheme. Eur. Phys. J. Spec. Top. 2020, 229, 1175–1188. [Google Scholar] [CrossRef]
Chaurra-Gutierrrez, F.A.; Rodriguez-Gomez, G.; Feregrino-Uribe, C.; Tlelo-Cuautle, E.; Guillen-Fernandez, O. Secure Communication Scheme Based on Projective Synchronization of Hyperchaotic Systems. In Cybersecurity: A New Approach Using Chaotic Systems; Abd El-Latif, A.A., Volos, C., Eds.; Springer International Publishing: Cham, Switzerland, 2022; pp. 109–156. [Google Scholar] [CrossRef]
Bodale, I.; Oancea, V.A. Chaos control for Willamowski–Rössler model of chemical reactions. Chaos Solitons Fractals 2015, 78, 1–9. [Google Scholar] [CrossRef]
Wang, C.; li Zhang, H.; Fan, W.H.; Ma, P. Finite-time function projective synchronization control method for chaotic wind power systems. Chaos Solitons Fractals 2020, 135, 109756. [Google Scholar] [CrossRef]
Zhao, X.; Li, Z.; Li, S. Synchronization of a chaotic finance system. Appl. Math. Comput. 2011, 217, 6031–6039. [Google Scholar] [CrossRef]
Idowu, B.A.; Vaidyanathan, S.; Sambas, A.; Olusola, O.I.; Onma, O.S. A New Chaotic Finance System: Its Analysis, Control, Synchronization and Circuit Design. In Nonlinear Dynamical Systems with Self-Excited and Hidden Attractors; Pham, V.T., Vaidyanathan, S., Volos, C., Kapitaniak, T., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 271–295. [Google Scholar] [CrossRef]
Zaqueros-Martinez, J.; Rodriguez-Gomez, G.; Tlelo-Cuautle, E.; Orihuela-Espina, F. Synchronization of Chaotic Electroencephalography (EEG) Signals. In Cybersecurity: A New Approach Using Chaotic Systems; Abd El-Latif, A.A., Volos, C., Eds.; Springer International Publishing: Cham, Switzerland, 2022; pp. 83–108. [Google Scholar] [CrossRef]
Alsaade, F.W.; Yao, Q.; Bekiros, S.; Al-zahrani, M.S.; Alzahrani, A.S.; Jahanshahi, H. Chaotic attitude synchronization and anti-synchronization of master-slave satellites using a robust fixed-time adaptive controller. Chaos Solitons Fractals 2022, 165, 112883. [Google Scholar] [CrossRef]
Munoz-Pacheco, J.M.; Posadas-Castillo, C.; Zambrano-Serrano, E. The Effect of a Non-Local Fractional Operator in an Asymmetrical Glucose-Insulin Regulatory System: Analysis, Synchronization and Electronic Implementation. Symmetry 2020, 12, 1395. [Google Scholar] [CrossRef]
Oluwole, O.S.A. Dynamic Regimes of El Niño Southern Oscillation and Influenza Pandemic Timing. Front. Public Health 2017, 5, 301. [Google Scholar] [CrossRef] [Green Version]
Bendoukha, S.; Abdelmalek, S.; Ouannas, A. Secure Communication Systems Based on the Synchronization of Chaotic Systems. In Mathematics Applied to Engineering, Modelling, and Social Issues; Chapter Secure Communication Systems Based on the Synchronization of Chaotic Systems; Springer International Publishing: Cham, Switzerland, 2019; pp. 281–311. [Google Scholar] [CrossRef]
Çavuşoğlu, Ü.; Panahi, S.; Akgül, A.; Jafari, S.; Kaçar, S. A new chaotic system with hidden attractor and its engineering applications: Analog circuit realization and image encryption. Analog Integr. Circuits Signal Process. 2019, 98, 85–99. [Google Scholar] [CrossRef]
Blagov, M.V.; Kuznetsov, N.V.; Leonov, G.A.; Yuldashev, M.V.; Yuldashev, R.V. Simulation of PLL with impulse signals in MATLAB: Limitations, hidden oscillations, and pull-in range. In Proceedings of the 2015 7th International Congress on Ultra Modern Telecommunications and Control Systems and Workshops (ICUMT), Brno, Czech Republic, 6–8 October 2015; pp. 85–90. [Google Scholar] [CrossRef]
Kuznetsov, N.V.; Leonov, G.A.; Mokaev, T.N.; Seledzhi, S.M. Hidden attractor in the Rabinovich system, Chua circuits and PLL. AIP Conf. Proc. 2016, 1738, 210008. [Google Scholar] [CrossRef]
Kiseleva, M.A.; Kudryashova, E.V.; Kuznetsov, N.V.; Kuznetsova, O.A.; Leonov, G.A.; Yuldashev, M.V.; Yuldashev, R.V. Hidden and self-excited attractors in Chua circuit: Synchronization and SPICE simulation. Int. J. Parallel Emergent Distrib. Syst. 2018, 33, 513–523. [Google Scholar] [CrossRef]
Lambert, J.D. Computational Methods in Ordinary Differential Equations; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 1973. [Google Scholar]
Herinci, P. Discrete Variable Methods in Ordinary Differential Equations; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1962. [Google Scholar]
Pecora, L.M.; Carroll, T.L. Synchronization of chaotic systems. Chaos Interdiscip. J. Nonlinear Sci. 2015, 25, 097611. [Google Scholar] [CrossRef] [Green Version]
Volos, C.K.; Pham, V.T.; Azar, A.T.; Stouboulos, I.N.; Kyprianidis, I.M. Synchronization Phenomena in Coupled Dynamical Systems with Hidden Attractors. In Nonlinear Dynamical Systems with Self-Excited and Hidden Attractors; Pham, V.T., Vaidyanathan, S., Volos, C., Kapitaniak, T., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 375–401. [Google Scholar] [CrossRef]
Tian, H.; Wang, Z.; Zhang, P.; Chen, M.; Wang, Y. Dynamic Analysis and Robust Control of a Chaotic System with Hidden Attractor. Complexity 2021, 2021, 8865522. [Google Scholar] [CrossRef]
Rajagopal, K.; Kingni, S.T.; Kom, G.H.; Pham, V.T.; Karthikeyan, A.; Jafari, S. Self-Excited and Hidden Attractors in a Simple Chaotic Jerk System and in Its Time-Delayed Form: Analysis, Electronic Implementation, and Synchronization. J. Korean Phys. Soc. 2020, 77, 145–152. [Google Scholar] [CrossRef]
Lai, Q.; Kamdem Kuate, P.D.; Pei, H.; Fotsin, H. Infinitely Many Coexisting Attractors in No-Equilibrium Chaotic System. Complexity 2020, 2020, 8175639. [Google Scholar] [CrossRef] [Green Version]
Pm, G.; Kapitaniak, T. Synchronization in Coupled Multistable Systems with Hidden Attractors. Math. Probl. Eng. 2017, 2017, 5214235. [Google Scholar] [CrossRef] [Green Version]
Wei, Z.; Yousefpour, A.; Jahanshahi, H.; Erkin Kocamaz, U.; Moroz, I. Hopf bifurcation and synchronization of a five-dimensional self-exciting homopolar disc dynamo using a new fuzzy disturbance-observer-based terminal sliding mode control. J. Frankl. Inst. 2021, 358, 814–833. [Google Scholar] [CrossRef]
Vaidyanathan, S.; He, S.; Sambas, A. A new multistable double-scroll 4-D hyperchaotic system with no equilibrium point, its bifurcation analysis, synchronization and circuit design. Arch. Control. Sci. 2021, 31, 99–128. [Google Scholar] [CrossRef]
Wang, S.; Yousefpour, A.; Yusuf, A.; Jahanshahi, H.; Alcaraz, R.; He, S.; Munoz-Pacheco, J.M. Synchronization of a Non-Equilibrium Four-Dimensional Chaotic System Using a Disturbance-Observer-Based Adaptive Terminal Sliding Mode Control Method. Entropy 2020, 22, 271. [Google Scholar] [CrossRef] [Green Version]
Vaidyanathan, S.; Moroz, I.M.; Sambas, A. A new 4-D hyperchaotic system with no equilibrium, its multistability, offset boosting and circuit simulation. Arch. Control. Sci. 2020, 30, 575–597. [Google Scholar] [CrossRef]
Ramadoss, J.; Kengne, R.; Tokoue Ngatcha, D.; Kamdoum Tamba, V.; Rajagopal, K.; Motchongom Tingue, M. A Three-Dimensional Autonomous System with a Parabolic Equilibrium: Dynamical Analysis, Adaptive Synchronization via Relay Coupling, and Applications to Steganography and Chaos Encryption. Complexity 2022, 2022, 8362836. [Google Scholar] [CrossRef]
Zhang, X.; Wang, C. Multiscroll Hyperchaotic System with Hidden Attractors and Its Circuit Implementation. Int. J. Bifurc. Chaos 2019, 29, 1950117. [Google Scholar] [CrossRef]
Tamba, V.K.; Pham, V.T.; Hoang, D.V.; Jafari, S.; Alsaadi, F.E.; Alsaadi, F.E. Dynamic system with no equilibrium and its chaos anti-synchronization. Automatika 2018, 59, 35–42. [Google Scholar] [CrossRef]
Wang, X.; Pham, V.T.; Jafari, S.; Volos, C.; Munoz-Pacheco, J.M.; Tlelo-Cuautle, E. A New Chaotic System With Stable Equilibrium: From Theoretical Model to Circuit Implementation. IEEE Access 2017, 5, 8851–8858. [Google Scholar] [CrossRef]
Elsonbaty, A.; El-Sayed, A.M.A. Analytical study of global bifurcations, stabilization and chaos synchronization of jerk system with multiple attractors. Nonlinear Dyn. 2017, 90, 2637–2655. [Google Scholar] [CrossRef]
Azar, A.T.; Volos, C.; Gerodimos, N.A.; Tombras, G.S.; Pham, V.T.; Radwan, A.G.; Vaidyanathan, S.; Ouannas, A.; Munoz-Pacheco, J.M. A Novel Chaotic System without Equilibrium: Dynamics, Synchronization, and Circuit Realization. Complexity 2017, 2017, 7871467. [Google Scholar] [CrossRef]
Tian, H.; Wang, Z.; Zhang, H.; Cao, Z.; Zhang, P. Dynamical analysis and fixed-time synchronization of a chaotic system with hidden attractor and a line equilibrium. Eur. Phys. J. Spec. Top. 2022, 231, 2455–2466. [Google Scholar] [CrossRef]
Munoz-Pacheco, J.M.; Volos, C.; Serrano, F.E.; Jafari, S.; Kengne, J.; Rajagopal, K. Stabilization and Synchronization of a Complex Hidden Attractor Chaotic System by Backstepping Technique. Entropy 2021, 23, 921. [Google Scholar] [CrossRef]
Zhang, H.; Liu, D.; Wang, Z. Controlling Chaos: Suppression, Synchronization and Chaotification; Springer: London, UK, 2009. [Google Scholar] [CrossRef]
Tanaka, K.; Ikeda, T.; Wang, H.O. A unified approach to controlling chaos via an lmi-based fuzzy control system design. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 1998, 45, 1021–1040. [Google Scholar] [CrossRef]
Skufca, J.D. Analysis Still Matters: A Surprising Instance of Failure of Runge–Kutta–Felberg ODE Solvers. SIAM Rev. 2004, 46, 729–737. [Google Scholar] [CrossRef] [Green Version]
Gautschi, W. Numerical integration of ordinary differential equations based on trigonometric polynomials. Numer. Math. 1961, 3, 381–397. [Google Scholar] [CrossRef]
Takagi, T.; Sugeno, M. Fuzzy identification of systems and its applications to modeling and control. IEEE Trans. Syst. Man Cybern. 1985, SMC-15, 116–132. [Google Scholar] [CrossRef]
Lian, K.Y.; Chiu, C.S.; Chiang, T.S.; Liu, P. LMI-based fuzzy chaotic synchronization and communications. IEEE Trans. Fuzzy Syst. 2001, 9, 539–553. [Google Scholar] [CrossRef]
Lian, K.Y.; Chiang, T.S.; Chiu, C.S.; Liu, P. Synthesis of fuzzy model-based designs to synchronization and secure communications for chaotic systems. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2001, 31, 66–83. [Google Scholar] [CrossRef] [PubMed]
Chua, L.O. A zoo of strange attractors from the canonical Chua’s circuits. In Proceedings of the [1992] 35th Midwest Symposium on Circuits and Systems, Washington, DC, USA, 9–12 August 1992; Volume 2, pp. 916–926. [Google Scholar] [CrossRef]
Leonov, G.A.; Kuznetsov, N.V.; Kuznetsova, O.A.; Seledzhi, S.M.; Vagaitsev, V.I. Hidden Oscillations in Dynamical Systems. Wseas Trans. Syst. Control 2011, 6, 54–67. [Google Scholar]
Sprott, J.C. Some simple chaotic flows. Phys. Rev. E 1994, 50, R647–R650. [Google Scholar] [CrossRef] [PubMed]
Jafari, S.; Sprott, J.; Hashemi Golpayegani, S.M.R. Elementary quadratic chaotic flows with no equilibria. Phys. Lett. A 2013, 377, 699–702. [Google Scholar] [CrossRef]
Jafari, S.; Sprott, J. Simple chaotic flows with a line equilibrium. Chaos Solitons Fractals 2013, 57, 79–84. [Google Scholar] [CrossRef]
Chua, L.O. The Genesis of Chua’s Circuit. Arch. Fur Elektron. Und Uebertragungstechnik 1992, 46, 250–257. [Google Scholar]
Kuznetsov, N.; Leonov, G.; Vagaitsev, V. Analytical-numerical method for attractor localization of generalized Chua’s system. IFAC Proc. Vol. 2010, 43, 29–33. [Google Scholar] [CrossRef]
Niemann, D.; Li, J.; Wang, H.O.; Tanaka, K. Parallel distributed compensation for Takagi-Sugeno fuzzy models: New stability conditions and dynamic feedback designs. IFAC Proc. Vol. 1999, 32, 5374–5379. [Google Scholar] [CrossRef]
Nian, Y.; Zheng, Y. Generalized projective synchronization of chaotic systems based on Takagi-Sugeno fuzzy model. In Proceedings of the 2010 Seventh International Conference on Fuzzy Systems and Knowledge Discovery, Yantai, China, 10–12 August 2010; Volume 3, pp. 1291–1295. [Google Scholar] [CrossRef]
Danca, M.F. Hidden chaotic attractors in fractional-order systems. Nonlinear Dyn. 2017, 89. [Google Scholar] [CrossRef] [Green Version]
Pano-Azucena, A.D.; Tlelo-Cuautle, E.; Rodriguez-Gomez, G.; de la Fraga, L.G. FPGA-based implementation of chaotic oscillators by applying the numerical method based on trigonometric polynomials. AIP Adv. 2018, 8, 75217. [Google Scholar] [CrossRef] [Green Version]
Zaqueros-Martinez, J.; Rodríguez-Gómez, G.; Tlelo-Cuautle, E.; Orihuela-Espina, F. Trigonometric polynomials methods to simulate oscillating chaotic systems. AIP Conf. Proc. 2022, 2425, 420035. [Google Scholar] [CrossRef]
Pham, V.T.; Volos, C.; Kapitaniak, T. Systems Without Equilibrium. In Systems with Hidden Attractors: From Theory to Realization in Circuits; Springer International Publishing: Cham, Switzerland, 2017; pp. 51–63. [Google Scholar] [CrossRef]
Baranyi, P. The Generalized TP Model Transformation for T–S Fuzzy Model Manipulation and Generalized Stability Verification. IEEE Trans. Fuzzy Syst. 2014, 22, 934–948. [Google Scholar] [CrossRef] [Green Version]
Baranyi, P. TP-Model Transformation-Based-Control Design Frameworks; Springer: Cham, Switzerland, 2016. [Google Scholar] [CrossRef]
Baranyi, P. How to Vary the Input Space of a T–S Fuzzy Model: A TP Model Transformation-Based Approach. IEEE Trans. Fuzzy Syst. 2022, 30, 345–356. [Google Scholar] [CrossRef]
Dudkowski, D.; Jafari, S.; Kapitaniak, T.; Kuznetsov, N.V.; Leonov, G.A.; Prasad, A. Hidden attractors in dynamical systems. Phys. Rep. 2016, 637, 1–50. [Google Scholar] [CrossRef]

Figure 1. Hidden attractors of the test systems. (a) Chua system; the three equilibrium points are shown with black color to highlight the presence of the hidden attractor. (b) Sprott A system. (c) NE₆ system. (d) LE₄ system; the line of equilibrium

x_{3}

has is plotted in black.

x_{3}

has is plotted in black.

Figure 2. Synchronization scheme. Depiction of a master–slave feedback control where the output variables

y_{i}

are required to follow the references (set points) generated by the master system

x_{i}

. The error

e_{i}

is the difference between the output and the reference. The control

u_{i}

forces the output to follow the reference.

Figure 2. Synchronization scheme. Depiction of a master–slave feedback control where the output variables

y_{i}

are required to follow the references (set points) generated by the master system

x_{i}

. The error

e_{i}

is the difference between the output and the reference. The control

u_{i}

forces the output to follow the reference.

Figure 3. Results of complete synchronizations. (a) Phase space, (b) error, and (c) control of the Chua system. Analogously for the Sprott A system (d–f), NE₆ system (g–i), and LE₄ system (j–l). For Chua’s system, the synchronization errors were less than

2 %

from Iteration 12,301 onwards. For all other systems, this happened around Iteration 3000. If the units of time are seconds(s), 3000 iterations would be equivalent to

23.43

s, representing

15 %

of the total simulated time (150 s).

2 %

from Iteration 12,301 onwards. For all other systems, this happened around Iteration 3000. If the units of time are seconds(s), 3000 iterations would be equivalent to

23.43

s, representing

15 %

of the total simulated time (150 s).

Figure 4. Standard deviation of the characteristics’ values. (a) Standard deviation of the Chua system. (b) Standard deviation of the LE₄ system. Each one of the functions in its respective interval presents an absolute minimum and an absolute maximum.

Figure 5. Errors of projective synchronizations. (a) Error in first state variable, (b) in the second variable, and (c) in the third variable of the Chua system. Analogously for the LE₄ system (d–f). For the Chua system, the synchronizations with

γ = 0.99

and

γ = 1.00

reached zero faster than with

γ = 0.88

and

γ = 0.89

. For the LE₄ system, the fastest convergences were with

γ = 0.99

and

γ = 1.0

and slowest with

γ = 0.70

and

γ = 0.71

. The range of variation for

γ

was

[0.88, 1]

for Chua and

[0.70, 1]

for LE₄. The maximum difference between the synchronization errors was approximately

0.25

in Chua and

0.04

in LE₄.

γ = 0.99

and

γ = 1.00

reached zero faster than with

γ = 0.88

and

γ = 0.89

. For the LE₄ system, the fastest convergences were with

γ = 0.99

and

γ = 1.0

and slowest with

γ = 0.70

and

γ = 0.71

. The range of variation for

γ

was

[0.88, 1]

for Chua and

[0.70, 1]

for LE₄. The maximum difference between the synchronization errors was approximately

0.25

in Chua and

0.04

in LE₄.

Figure 6. Iteration at which error synchronization is less than

2 %

for all 30 of the newly found initial conditions. While the complete and projective synchronizations with

α = 1

are theoretically analogous, note that the control design differs and, hence, the observed differences.

Figure 6. Iteration at which error synchronization is less than

2 %

for all 30 of the newly found initial conditions. While the complete and projective synchronizations with

α = 1

are theoretically analogous, note that the control design differs and, hence, the observed differences.

Table 1. Chaotic systems with hidden attractors. In the first column, the name of each system is presented. In the second column, the mathematical model is given. In the third column, the parameters of the system are cited if they do not come from the seminal work. The fourth column indicates whether the systems have equilibrium points.

System Name	Mathematical Model	Parameters	Eq. Point
Chua	${\dot{x}}_{1} = α (x_{2} - x_{1} - ψ (x_{1}))$	$α = 8.4562$ ,	Yes
[44]	${\dot{x}}_{2} = x_{1} - x_{2} + x_{3}$	$β = 12.0732$ ,
	${\dot{x}}_{3} = - (β x_{2} + γ x_{3})$	$γ = 0.0052$ ,
	$ψ (x_{1}) = m_{1} x_{1} +$	$m_{0} = - 0.1768$ ,
	$(m_{0} - m_{1})$	$m_{1} = - 1.1468$
	$sat (x_{1})$	[45]
Sprott A	${\dot{x}}_{1} = x_{2}$		No
[46]	${\dot{x}}_{2} = - x_{1} + x_{2} x_{3}$
	${\dot{x}}_{3} = 1 - x_{2}^{2}$
NE₆	${\dot{x}}_{1} = x_{2}$		No
[47]	${\dot{x}}_{2} = x_{3}$
	${\dot{x}}_{3} = - x_{2} - x_{1} x_{3} -$
	$x_{2} x_{3} - a$	$a = 0.75$
LE₄	${\dot{x}}_{1} = x_{2}$		Line
[48]	${\dot{x}}_{2} = - x_{1} + x_{2} x_{3}$	$a = 4$
	${\dot{x}}_{3} = - a x_{1} x_{2} - b x_{1} x_{2} -$	$b = 0.6$
	$x_{2} x_{3}$

Table 2. Performance comparison of fuzzy synchronizations using different Takagi–Sugeno (T–S) fuzzy models. The name of the chaotic system is in the first column. The T–S fuzzy models are in the second column: Base refers to the T–S fuzzy models given in Section 5.3. Sum normalization non-negativeness (SNNN) and close-to-normality (CNO) are the variants obtained with the transformation of the tensor product (TP) model. The first iterations in which the synchronization errors satisfied the given comparison criteria using complete and projective fuzzy control are reported in the third and fourth columns, respectively.

System	T–S Fuzzy Model	Complete	Projective
Chua	Base	12,301	1670
	SNNN	2496	1281
	CNO	2496	1281
Sprott A	Base	2792
	SNNN	2663
	CNO	2663
NE₆	Base	1885
	SNNN	1885
	CNO	1885
LE₄	Base	3383	1771
	SNNN	3383	1923
	CNO	3383	1771

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zaqueros-Martinez, J.; Rodriguez-Gomez, G.; Tlelo-Cuautle, E.; Orihuela-Espina, F. Fuzzy Synchronization of Chaotic Systems with Hidden Attractors. Entropy 2023, 25, 495. https://doi.org/10.3390/e25030495

AMA Style

Zaqueros-Martinez J, Rodriguez-Gomez G, Tlelo-Cuautle E, Orihuela-Espina F. Fuzzy Synchronization of Chaotic Systems with Hidden Attractors. Entropy. 2023; 25(3):495. https://doi.org/10.3390/e25030495

Chicago/Turabian Style

Zaqueros-Martinez, Jessica, Gustavo Rodriguez-Gomez, Esteban Tlelo-Cuautle, and Felipe Orihuela-Espina. 2023. "Fuzzy Synchronization of Chaotic Systems with Hidden Attractors" Entropy 25, no. 3: 495. https://doi.org/10.3390/e25030495

APA Style

Zaqueros-Martinez, J., Rodriguez-Gomez, G., Tlelo-Cuautle, E., & Orihuela-Espina, F. (2023). Fuzzy Synchronization of Chaotic Systems with Hidden Attractors. Entropy, 25(3), 495. https://doi.org/10.3390/e25030495

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Fuzzy Synchronization of Chaotic Systems with Hidden Attractors

Abstract

1. Introduction

2. Related Work

3. Global Error

4. Special Numerical Methods

5. Fuzzy Control

5.1. Fuzzification

5.2. Defuzzification

5.3. Fuzzy Modeling of Chaotic Systems

6. Fuzzy Synchronization

6.1. Complete Fuzzy Synchronization

6.2. Projective Fuzzy Synchronization

7. Projective Synchronizations with an Emphasis on Error Convergence

8. Simulation Results

8.1. Results of Fuzzy Complete Synchronization

8.1.1. Chua

8.1.2. Sprott A

8.1.3. NE6

8.1.4. LE4

8.2. Results of Fuzzy Projective Synchronization

8.2.1. Chua System

8.2.2. LE4

8.3. A New Set of Initial Conditions

8.4. Alternative T–S Fuzzy Models

9. Discussion

10. Conclusions and Future Work

Author Contributions

Funding

Data Availability Statement

Acknowledgments

Conflicts of Interest

Appendix A

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

8.1.3. NE₆

8.1.4. LE₄

8.2.2. LE₄