Fuzzy Synchronization of Chaotic Systems with Hidden Attractors

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.

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\left(t\right)=y_n-y\left(t_n\right)$, where $y_n$ is the approximation to the theoretical value $y\left(t_n\right)$ 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, $\alpha$ and $\gamma$. We further propose here a new strategy for selecting the control parameter $\gamma$, then we manipulated the scaling parameter $\alpha$, 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\left(t\right)$ of the numerical approximation of a solution of ordinary differential equations with initial values is bounded according to

$|e_n\left(t\right)|<\{A\delta +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 $\delta$ 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 $\delta$, 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\left(t\right)$, 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({\beta }_1{f}_{n+1}+{\beta }_0{f}_n)$

where

$\begin{array}{ccc}\hfill {\beta }_0\left(\nu \right)& =& -\frac{1}{2}(1+\frac{1}{12}{\nu }^2+\frac{1}{120}{\nu }^4+\dots )\hfill \\ \hfill {\beta }_1\left(\nu \right)& =& \frac{3}{2}(1-\frac{1}{4}{\nu }^2+\frac{1}{120}{\nu }^4+\dots )\hfill \end{array}$

with $\nu =\omega h$, the integration step h constant, and the truncation error $\frac{5}{12}{h}^3(y^{\left(3\right)}+{\omega }^2{y}^{\left(1\right)})$. If the frequency $\omega$ 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

(1)$\dot{\mathbf{x}}=\mathbf{f}\left(\mathbf{x}\left(t\right)\right)$

where $\mathbf{x}\left(t\right)={[x_1\left(t\right),x_2\left(t\right),\dots ,x_n\left(t\right)]}^T\in {\mathbb{R}}^n$ is the state vector and $\mathbf{f}$ is a nonlinear function of appropriate dimensions. The T–S fuzzy model is composed of a set of rules of the form:

(2)$\begin{array}{cc}\mathrm{Rule}i:\hfill & \mathbf{IF}{z}_1\left(t\right)\mathrm{is}{F}_{i1}\mathrm{and}\dots \hspace{0.17em}\mathrm{and}{z}_p\left(t\right)\mathrm{is}{F}_{ip}\\ & \mathbf{THEN}\dot{\mathbf{x}}={\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{b}}_i,i=1,2,3,\dots ,I\end{array}$

where $z_1\left(t\right)\dots z_p\left(t\right)$ are the premise state variables, $F_{i,j}$ ($j=1,\dots ,p$) are fuzzy sets, ${\mathbf{A}}_i$ is a constant system matrix with the appropriate size, and ${\mathbf{b}}_i\in {\mathbb{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

(3)$\dot{\mathbf{x}}=\frac{{\sum }_{i=1}^I{\omega }_i\left(\mathbf{z}\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{b}}_i\}}{{\sum }_{i=1}^I{\omega }_i\left(\mathbf{z}\left(t\right)\right)}$

where $\mathbf{z}\left(t\right)={[z_1\left(t\right)z_2\left(t\right)\dots z_p\left(t\right)]}^T,$ with $w_i\left(\mathbf{z}\left(t\right)\right)={\prod }_{j=1}^p{F}_{ij}\left(z_j\left(t\right)\right)$ in which $F_{ij}\left(z_j\left(t\right)\right)$ is the degree of membership of $z_j\left(t\right)$ in $F_{ij}$, with ${\sum }_{i=1}^I{\omega }_i\left(\mathbf{z}\left(t\right)\right)>0$, and ${\omega }_i\left(\mathbf{z}\left(t\right)\right)\ge 0$, $i=1,2,3,\dots ,I$.

By employing ${\mu }_i\left(\mathbf{z}\left(t\right)\right)={\omega }_i\left(\mathbf{z}\left(t\right)\right)/{\sum }_{i=1}^I{\omega }_i\left(\mathbf{z}\left(t\right)\right)$ instead of ${\omega }_i\left(\mathbf{z}\left(t\right)\right)$, System (3) is rewritten as

(4)$\dot{\mathbf{x}}=\sum _{i=1}^I{\mu }_i\left(\mathbf{z}\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{b}}_i\}$

where ${\sum }_{i=1}^I{\mu }_i\left(\mathbf{z}\left(t\right)\right)=1$ with ${\mu }_i\left(\mathbf{z}\left(t\right)\right)\ge 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 $\mathbf{f}\left(\mathbf{x}\right(t\left)\right)$ is expressed as the inferred fuzzy output ${\sum }_{i=1}^I{\mu }_i\left(\mathbf{z}\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{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

(5)$\begin{array}{c}\hfill \mathrm{Rule}1:\mathbf{IF}{x}_1\left(t\right)\mathrm{is}{F}_1\mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_1\mathbf{x}\left(t\right)+{\mathbf{b}}_1\\ \hfill \mathrm{Rule}2:\mathbf{IF}{x}_1\left(t\right)\mathrm{is}{F}_2\mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_2\mathbf{x}\left(t\right)+{\mathbf{b}}_2\end{array}$

where $\mathbf{x}={[x_1,x_2,x_3]}^T$, ${\mathbf{b}}_1={\mathbf{b}}_2={[0,0,0]}^T$, $d=1.1$, $F_1\left(x_1\right)=\frac{1}{2}\left(1-\frac{\varphi \left(x_1\right)}{d}\right)$, $F_2\left(x_1\right)=1-F_1\left(x_1\right)$, and 

${\mathbf{A}}_1=\left(\begin{array}{ccc}\alpha (d-1)& \alpha & 0\\ 1& -1& 1\\ 0& -\beta & -\gamma \end{array}\right){\mathbf{A}}_2=\left(\begin{array}{ccc}-\alpha (d+1)& \alpha & 0\\ 1& -1& 1\\ 0& -\beta & -\gamma \end{array}\right)$

$\varphi \left(x_1\right)=\left\{\begin{array}{cc}\frac{\psi \left(x_1\right)}{x_1}& x_1\ne 0\\ m_0& x_1=0\end{array}\right.\mathrm{sat}\left(x_1\right)=\left\{\begin{array}{cc}-1\hfill & \mathrm{if}{x}_1\le -1\hfill \\ x_1\hfill & \mathrm{if}\mid x_1\mid <1\hfill \\ 1\hfill & \mathrm{if}{x}_1\ge 1\hfill \end{array}\right.$

The final output of the Chua fuzzy system is inferred by

$\dot{\mathbf{x}}=\sum _{i=1}^2{F}_i\left(x_1\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{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{array}{c}\hfill \mathrm{Rule}1:\mathbf{IF}{x}_2\left(t\right)\mathrm{is}{F}_1\mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_1\mathbf{x}\left(t\right)+{\mathbf{b}}_1\\ \hfill \mathrm{Rule}2:\mathbf{IF}{x}_2\left(t\right)\mathrm{is}{F}_2\mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_2\mathbf{x}\left(t\right)+{\mathbf{b}}_2\end{array}$

where $\mathbf{x}={[x_1,x_2,x_3]}^T$, ${\mathbf{b}}_1={\mathbf{b}}_2={[0,0,2]}^T$, $d=5.0$, $F_1\left(x_2\right)=\frac{1}{2}\left(1+\frac{x_2}{d}\right)$, $F_2\left(x_2\right)=\frac{1}{2}\left(1-\frac{x_2}{d}\right)$, and 

${\mathbf{A}}_1=\left(\begin{array}{ccc}0& 1& 0\\ -1& 0& d\\ 0& -d& 0\end{array}\right){\mathbf{A}}_2=\left(\begin{array}{ccc}0& 1& 0\\ -1& 0& -d\\ 0& d& 0\end{array}\right)$

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

$\dot{\mathbf{x}}=\sum _{i=1}^2{F}_i\left(x_2\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{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{array}{c}\hfill \mathrm{Rule}1:\mathbf{IF}{x}_3\left(t\right)\mathrm{is}{F}_1\mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_1\mathbf{x}\left(t\right)+{\mathbf{b}}_1\\ \hfill \mathrm{Rule}2:\mathbf{IF}{x}_3\left(t\right)\mathrm{is}{F}_2\mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_2\mathbf{x}\left(t\right)+{\mathbf{b}}_2\end{array}$

where $\mathbf{x}={[x_1,x_2,x_3]}^T$, ${\mathbf{b}}_1={\mathbf{b}}_2={[0,0,-0.75]}^T$, $d=3.0$, $F_1\left(x_3\right)=\frac{1}{2}\left(1+\frac{x_3}{d}\right)$, $F_2\left(x_3\right)=\frac{1}{2}\left(1-\frac{x_3}{d}\right)$, and 

${\mathbf{A}}_1=\left(\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ -d& -1-d& 0\end{array}\right){\mathbf{A}}_2=\left(\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ d& -1+d& 0\end{array}\right)$

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

$\dot{\mathbf{x}}=\sum _{i=1}^2{F}_i\left(x_3\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{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{array}{c}\hfill \mathrm{Rule}1:\mathbf{IF}{x}_2\left(t\right)\mathrm{is}{F}_1\mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_1\mathbf{x}\left(t\right)+{\mathbf{b}}_1\\ \hfill \mathrm{Rule}2:\mathbf{IF}{x}_2\left(t\right)\mathrm{is}{F}_2\mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_2\mathbf{x}\left(t\right)+{\mathbf{b}}_2\end{array}$

where $\mathbf{x}={[x_1,x_2,x_3]}^T$, ${\mathbf{b}}_1={\mathbf{b}}_2={[0,0,0]}^T$, $d=0.8$, $F_1\left(x_2\right)=\frac{1}{2}\left(1+\frac{x_2}{d}\right)$, $F_2\left(x_2\right)=\frac{1}{2}\left(1-\frac{x_2}{d}\right)$, and 

${\mathbf{A}}_1=\left(\begin{array}{ccc}0& 1& 0\\ -1& 0& d\\ -d(a+b)& 0& -d\end{array}\right){\mathbf{A}}_2=\left(\begin{array}{ccc}0& 1& 0\\ -1& 0& -d\\ d(a+b)& 0& d\end{array}\right)$

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

$\dot{\mathbf{x}}=\sum _{i=1}^2{F}_i\left(x_2\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{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:

(6)$\begin{array}{c}\hfill \mathrm{Rule}1:\mathbf{IF}{y}_1\left(t\right)\mathrm{is}{F}_1\mathbf{THEN}\dot{\mathbf{y}}\left(t\right)={\mathbf{A}}_1\mathbf{y}\left(t\right)+{\mathbf{b}}_1+\mathbf{B}\mathbf{u}\left(t\right)\\ \hfill \mathrm{Rule}2:\mathbf{IF}{y}_1\left(t\right)\mathrm{is}{F}_2\mathbf{THEN}\dot{\mathbf{y}}\left(t\right)={\mathbf{A}}_2\mathbf{y}\left(t\right)+{\mathbf{b}}_2+\mathbf{B}\mathbf{u}\left(t\right)\end{array}$

where $\mathbf{y}={[y_1,y_2,y_3]}^T$ are state variables and ${\mathbf{A}}_1$, ${\mathbf{A}}_2$, ${\mathbf{b}}_1$, ${\mathbf{b}}_2$, $F_1$, $F_2$, and d have the same values as the master system (2). $\mathbf{B}$ is an input matrix, and $\mathbf{u}\left(t\right)$ is the control.

Defuzzification is given by

$\begin{array}{c}\hfill \dot{\mathbf{y}}=\sum _{i=1}^2{\mu }_i\left(y_1\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{y}\left(t\right)+{\mathbf{b}}_i\}+\mathbf{B}\mathbf{u}\left(t\right),{\mu }_i\left(y_1\left(t\right)\right)=F_i\left(y_1\left(t\right)\right).\end{array}$

The error is $\mathbf{e}\left(t\right)=\mathbf{y}\left(t\right)-\mathbf{x}\left(t\right)$. The error can be estimated from

(7)$\begin{array}{c}\hfill \dot{\mathbf{e}}\left(t\right)=\sum _{i=1}^2{\mu }_i\left(y_1\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{y}\left(t\right)+{\mathbf{b}}_i\}-\sum _{i=1}^2{\mu }_i\left(x_1\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{x}\left(t\right)+{\mathbf{b}}_i\}+\mathbf{B}\mathbf{u}\left(t\right).\end{array}$

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

(8)$\begin{array}{c}\hfill \mathbf{u}\left(t\right)=-\sum _{i=1}^2{\mu }_i\left(y_1\left(t\right)\right)\{{\mathbf{C}}_i\mathbf{y}\left(t\right)+{\mathbf{b}}_i\}+\sum _{i=1}^2{\mu }_i\left(x_1\left(t\right)\right)\{{\mathbf{C}}_i\mathbf{x}\left(t\right)+{\mathbf{b}}_i\}\end{array}$

such that $∥\mathbf{e}\left(t\right)∥\to 0$ when $t\to \infty$. Substituting the fuzzy control (8) in the error (7), one obtains

$\begin{array}{c}\hfill \dot{\mathbf{e}}\left(t\right)=\sum _{i=1}^2{\mu }_i\left(y_1\left(t\right)\right)\left\{({\mathbf{A}}_i-\mathbf{B}{\mathbf{C}}_i)\mathbf{y}\left(t\right)\right\}-\sum _{i=1}^2{\mu }_i\left(x_1\left(t\right)\right)\left\{({\mathbf{A}}_i-\mathbf{B}{\mathbf{C}}_i)\mathbf{x}\left(t\right)\right\}\end{array}$

The gain matrices ${\mathbf{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{array}{c}\hfill {\{({\mathbf{A}}_1-\mathbf{B}{\mathbf{C}}_1)-({\mathbf{A}}_2-\mathbf{B}{\mathbf{C}}_2)\}}^T\times \{({\mathbf{A}}_1-\mathbf{B}{\mathbf{C}}_1)-({\mathbf{A}}_2-\mathbf{B}{\mathbf{C}}_2)\}=0\end{array}$

then the total error of the system becomes linear when $\dot{\mathbf{e}}\left(t\right)=\mathbf{G}\mathbf{e}\left(t\right)$ according to the fuzzy control (8), where $\mathbf{G}={\mathbf{A}}_1-\mathbf{B}{\mathbf{C}}_1={\mathbf{A}}_2-\mathbf{B}{\mathbf{C}}_2$. Moreover, if $\mathbf{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 $\alpha$. Let the master system (1) and the slave system be $\dot{\mathbf{y}}=\mathbf{g}(\mathbf{y}\left(t\right),\mathbf{u}(\mathbf{x}\left(t\right),\mathbf{y}\left(t\right)))$. If there exist a constant $\alpha \ne 0$ such that ${lim}_{t\to \infty }∥\mathbf{y}\left(t\right)-\alpha \mathbf{x}\left(t\right)∥=0$, then the synchronization is referred to as projective, and $\alpha$ 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).

(9)$\begin{array}{cc}\hfill \mathrm{Rule}1:\mathbf{IF}{x}_1\left(t\right)\mathrm{is}{F}_1& \mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_1\mathbf{x}\left(t\right)\hfill \\ \hfill \mathrm{Rule}2:\mathbf{IF}{x}_1\left(t\right)\mathrm{is}{F}_2& \mathbf{THEN}\dot{\mathbf{x}}\left(t\right)={\mathbf{A}}_2\mathbf{x}\left(t\right)\hfill \end{array}$

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

(10)$\dot{\mathbf{x}}=\sum _{i=1}^2{\mu }_i\left(x_1\left(t\right)\right){\mathbf{A}}_i\mathbf{x}\left(t\right)$

where $\mu$ satisfies the same conditions as (4).

(11)$\begin{array}{cc}\hfill \mathrm{Rule}1:\mathbf{IF}{x}_1\left(t\right)\mathrm{is}{F}_1& \mathbf{THEN}\dot{\mathbf{y}}\left(t\right)={\mathbf{A}}_1\mathbf{y}\left(t\right)+\mathbf{u}\left(t\right)\hfill \\ \hfill \mathrm{Rule}2:\mathbf{IF}{x}_1\left(t\right)\mathrm{is}{F}_2& \mathbf{THEN}\dot{\mathbf{y}}\left(t\right)={\mathbf{A}}_2\mathbf{y}\left(t\right)+\mathbf{u}\left(t\right)\hfill \end{array}$

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

(12)$\dot{\mathbf{y}}=\sum _{i=1}^2{\mu }_i\left(x_1\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{y}\left(t\right)+\mathbf{u}\left(t\right)\}$

Let the error $\mathbf{e}\left(t\right)=\mathbf{y}\left(t\right)-\alpha \mathbf{x}\left(t\right)$. Substituting the final outputs of the fuzzy systems (10) and (12), we obtain the derivative of the error:

(13)$\dot{\mathbf{e}}\left(t\right)=\sum _{i=1}^2{\mu }_i\left(x_1\left(t\right)\right)\{{\mathbf{A}}_i\mathbf{e}\left(t\right)+\mathbf{u}\left(t\right)\}$

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{array}{cc}\hfill \mathrm{Rule}1:\mathbf{IF}{x}_1\left(t\right)\mathrm{is}{F}_1\mathbf{THEN}& \mathbf{u}\left(t\right)=-\gamma {\mathbf{A}}_1[\mathbf{y}\left(t\right)-\alpha \mathbf{x}\left(t\right)]-\gamma [\mathbf{y}\left(t\right)-\alpha \mathbf{x}\left(t\right)]\hfill \\ \hfill \mathrm{Rule}2:\mathbf{IF}{x}_1\left(t\right)\mathrm{is}{F}_2\mathbf{THEN}& \mathbf{u}\left(t\right)=-\gamma A_2[\mathbf{y}\left(t\right)-\alpha \mathbf{x}\left(t\right)]-\gamma [\mathbf{y}\left(t\right)-\alpha \mathbf{x}\left(t\right)]\hfill \end{array}$

where $\gamma$ is a control parameter.

The final output of the fuzzy control is given by:

(14)$\begin{array}{c}\hfill \mathbf{u}\left(t\right)=\sum _{i=1}^2{\mu }_i\left(x_1\left(t\right)\right)\{-\gamma {\mathbf{A}}_i[\mathbf{y}\left(t\right)-\alpha \mathbf{x}\left(t\right)]-\gamma [\mathbf{y}\left(t\right)-\alpha \mathbf{x}\left(t\right)]\}\end{array}$

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

(15)$\dot{\mathbf{e}}\left(t\right)=\sum _{i=1}^2{\mu }_i\left(x_1\left(t\right)\right)[(1-\gamma ){\mathbf{A}}_i-\gamma \mathbf{I}]\mathbf{e}\left(t\right)$

where $\mathbf{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-\gamma )A_i-\gamma I){]}^{T}P+P[(1-\gamma )A_i-\gamma I)]\le -cI$, $i=1,2$, then the error (15) is globally exponentially stable.

Another convergence condition for the error (15) concerns the characteristic values ${\lambda }_j$, $j=1,\dots ,n$ of the symmetric matrices $(1-\gamma )(A_i^T+A_i)-2\gamma I$, $i=1,2$. If every ${\lambda }_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 $\gamma$ analytically so that the convergence of the synchronization error to zero can be made faster or slower according to the application. Through matrices $(1-\gamma )(A_i^T+A_i)-2\gamma I$, $i=1,2$, the characteristic values are found as a function of $\gamma$, and its range of variation is determined. Parameter $\gamma$ 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 $\sigma \left(\lambda \right)$, $\lambda =({\lambda }_1,\dots ,{\lambda }_n)$, found in the valid range of $\gamma$, the set of characteristic values that gives us the sought convergence (fast or slow). When $\sigma \left(\lambda \right)=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 ${\lambda }^{*}=argmin\sigma \left(\lambda \right)=0$, in the range of variation of $\gamma$, the value of ${\lambda }^{*}$ is found so that the error converges as fast as possible. Otherwise, ${\lambda }^{*}=argmax\sigma \left(\lambda \right)$ 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 $\gamma$ are polynomials of a degree greater than 1, finding ${\lambda }^{*}$ is laborious in practice. In these cases (such as the ones discussed in this article), it is recommended to obtain the plot of $\gamma$ vs. $\sigma \left(\lambda \right(\gamma \left)\right)$ to determine the desired ${\lambda }^{*}$ and the range of variability of $\gamma$.

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 $\gamma$ 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 $\mathbf{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 ${\mathbf{x}}_0={[-3.7727,-1.3511,4.6657]}^T$ for the master system and ${\mathbf{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.

(16)${\mathbf{C}}_1=\left(\begin{array}{ccc}0.3364& 1.1820& 0\\ 1.1820& -0.1250& -1.3841\\ 0& -1.3841& 0.1237\end{array}\right),{\mathbf{C}}_2=\left(\begin{array}{ccc}-4.3145& 1.1820& 0\\ 1.1820& -0.1250& -1.3841\\ 0& -1.3841& 0.1237\end{array}\right)$

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 $∥\mathbf{e}\left(t\right)∥\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 ${\mathbf{x}}_0={[0,5,0]}^T$ for the master system [56] and ${\mathbf{y}}_0={[1.1,6.1,-0.1]}^T$ for the slave system. The calculated gain matrices for this system are

$\begin{array}{cc}\hfill {\mathbf{C}}_1={\mathbf{C}}_2& =\left(\begin{array}{ccc}0.5& 0& 0\\ 0& 0.5& 0\\ 0& 0& 0.5\end{array}\right)\end{array}$

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 ${\mathbf{x}}_0={[0,3,-0.1]}^T$ for the master system [47] and ${\mathbf{y}}_0={[2.4,1.4,-0.6]}^T$ for the slave system. The gain matrices obtained with this system are

${\mathbf{C}}_1=\left(\begin{array}{ccc}0.5& 0.5& -1.5\\ 0.5& 0.5& -1.5\\ -1.5& -1.5& 0.5\end{array}\right),{\mathbf{C}}_2=\left(\begin{array}{ccc}0.5& 0.5& 1.5\\ 0.5& 0.5& 1.5\\ 1.5& 1.5& 0.5\end{array}\right)$

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 ${\mathbf{x}}_0={[0.2,0.7,0]}^T$ for the master system [48] and ${\mathbf{y}}_0={[-0.1,0.3,0.1]}^T$ for the slave system. The gain matrices obtained are

${\mathbf{C}}_1=\left(\begin{array}{ccc}0.5& 0& -1.84\\ 0& 0.5& 0.4\\ -1.84& 0.4& -0.3\end{array}\right),{\mathbf{C}}_2=\left(\begin{array}{ccc}0.5& 0& 1.84\\ 0& 0.5& -0.4\\ 1.84& -0.4& 1.3\end{array}\right)$

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 $\gamma$ 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-\gamma )(A_i^T+A_i)-2\gamma I$, $i=1,2$, which are

$\begin{array}{c}\left(\begin{array}{ccc}1.69124(1-\gamma )-2\gamma & 9.4562(1-\gamma )& 0\\ 9.4562(1-\gamma )& -2(1-\gamma )-2\gamma & -11.0732(1-\gamma )\\ 0& -11.0732(1-\gamma )& 0.0104(1-\gamma )-2\gamma \end{array}\right)\end{array}$

$\begin{array}{c}\left(\begin{array}{ccc}-35.516(1-\gamma )-2\gamma & 9.4562(1-\gamma )& 0\\ 9.4562(1-\gamma )& -2(1-\gamma )-2\gamma & -11.0732(1-\gamma )\\ 0& -11.0732(1-\gamma )& -0.0104(1-\gamma )-2\gamma \end{array}\right)\end{array}$

Subsequently, we looked for the characteristic values of each matrix to determine the range of values in which $\gamma$ must lie to satisfy that the symmetric matrices have strictly negative characteristic values. The equation to find $\gamma$ is nonlinear. Clearing $\gamma$ analytically can be cumbersome. Hence, instead, we opted to plot the standard deviation of the characteristic values by varying the parameter $\gamma$ in the interval $[0.88,1]$ in Figure 4a. Such an interval can be found with trial and error with educated guesses; the sought $\gamma$ values are those complying with the required convergence conditions.

Following the selection of the interval of the $\gamma$ values that ensure convergence, we now present the synchronization under the lower and upper values of the function $\sigma \left(\lambda \right(\gamma \left)\right)$ within the interval where the convergence requirements are met: $\gamma =0.88$ and $\gamma =1.0$, respectively, in the case of Chua. As the process involved some trial and error, we further show the synchronization for nearby values: $\gamma =0.89$ and $\gamma =0.99$. The same initial conditions used for the complete synchronization were used here: ${\mathbf{x}}_0={[-3.7727,-1.3511,4.6657]}^T$ for the master system and ${\mathbf{y}}_0={[0.01,0,0]}^T$ for the slave system [15]. The scaling factor used was $\alpha =0.5$. The errors of each state variable of the synchronizations are plotted in Figure 5a–c. The error in the synchronizations with $\gamma =0.99$ and $\gamma =1$ converged to zero in fewer iterations than with $\gamma =0.88$ and $\gamma =0.89$ in each error. It can be appreciated how the value of $\gamma$ affected the speed of convergence: faster with $\gamma =0.99$ and $\gamma =1$ and slower with $\gamma =0.88$ and $\gamma =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: ${\mathbf{x}}_0={[0.2,0.7,0]}^T$ for the master system [48] and ${\mathbf{y}}_0={[-0.1,0.3,0.1]}^T$ for the slave system. We chose the scaling factor $\alpha =0.5$.

The symmetric matrices of this system are

$\begin{array}{c}\left(\begin{array}{ccc}-2\gamma & 0& -3.68(1-\gamma )\\ 0& -2\gamma & 0.8(1-\gamma )\\ -3.68(1-\gamma )& 0.8(1-\gamma )& -1.6(1-\gamma )-2\gamma \end{array}\right)\end{array}$

$\begin{array}{c}\left(\begin{array}{ccc}-2\gamma & 0& 3.68(1-\gamma )\\ 0& -2\gamma & -0.8(1-\gamma )\\ 3.68(1-\gamma )& -0.8(1-\gamma )& 1.6(1-\gamma )-2\gamma \end{array}\right)\end{array}$

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 $\gamma$ vs. $\sigma \left(\lambda \right(\gamma \left)\right)$ for $\gamma \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 $\gamma =0.77$ and $\gamma =1.0$, respectively. Again, because of the trial and error process, we also synchronized with nearby values: $\gamma =0.78$ and $\gamma =0.99$.

Figure 5d–f plot the error for each state variable, and again, we show how the value of $\gamma$ affects the speed of convergence;: faster with $\gamma =0.99$ and $\gamma =1.0$ and slower with $\gamma =0.70$ and $\gamma =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 $\alpha$ equal to $0.5$, 1, 2, 3, 4, and 8. The control parameter was fixed to $\gamma =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 $\alpha =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 $\alpha =1$ synchronizations. There was a significant difference between the scores for the complete (mean ($\mu$) ± standard deviation ($\sigma$): $\mu \pm \sigma =13,573\pm 2032.9$) and projective ($\mu \pm \sigma =1641.2\pm 256.85$) conditions; $t\left(29\right)=32.5046,p<0.001$.

We conducted a repeated-measures ANOVA test (Matlab’s R2021b ranova) to analyze the effect of the scaling factor $\alpha$ on the six tested projective synchronizations $\alpha$ = 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 $\alpha =0.5$, and the remaining $\alpha$s, $p<0.05$, as well as $\alpha =8$ and the others $\alpha$s, $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 $\gamma$. 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.

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.

Footnotes

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.

Figures and Tables

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) NE6 system. (d) LE4 system; the line of equilibrium [Forumla omitted. See PDF.] has is plotted in black.

Figure 2. Synchronization scheme. Depiction of a master–slave feedback control where the output variables [Forumla omitted. See PDF.] are required to follow the references (set points) generated by the master system [Forumla omitted. See PDF.]. The error [Forumla omitted. See PDF.] is the difference between the output and the reference. The control [Forumla omitted. See PDF.] 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), NE6 system (g–i), and LE4 system (j–l). For Chua’s system, the synchronization errors were less than [Forumla omitted. See PDF.] 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 [Forumla omitted. See PDF.] s, representing [Forumla omitted. See PDF.] of the total simulated time (150 s).

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), NE6 system (g–i), and LE4 system (j–l). For Chua’s system, the synchronization errors were less than [Forumla omitted. See PDF.] 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 [Forumla omitted. See PDF.] s, representing [Forumla omitted. See PDF.] 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 LE4 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 LE4 system (d–f). For the Chua system, the synchronizations with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] reached zero faster than with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. For the LE4 system, the fastest convergences were with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] and slowest with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. The range of variation for [Forumla omitted. See PDF.] was [Forumla omitted. See PDF.] for Chua and [Forumla omitted. See PDF.] for LE4. The maximum difference between the synchronization errors was approximately [Forumla omitted. See PDF.] in Chua and [Forumla omitted. See PDF.] in LE4.

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 LE4 system (d–f). For the Chua system, the synchronizations with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] reached zero faster than with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. For the LE4 system, the fastest convergences were with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] and slowest with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. The range of variation for [Forumla omitted. See PDF.] was [Forumla omitted. See PDF.] for Chua and [Forumla omitted. See PDF.] for LE4. The maximum difference between the synchronization errors was approximately [Forumla omitted. See PDF.] in Chua and [Forumla omitted. See PDF.] in LE4.

Figure 6. Iteration at which error synchronization is less than [Forumla omitted. See PDF.] for all 30 of the newly found initial conditions. While the complete and projective synchronizations with [Forumla omitted. See PDF.] 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=\alpha (x_2-x_1-\psi \left(x_1\right))$	$\alpha =8.4562$,	Yes
[44]	${\dot{x}}_2=x_1-x_2+x_3$	$\beta =12.0732$,
	${\dot{x}}_3=-(\beta x_2+\gamma x_3)$	$\gamma =0.0052$,
	$\psi \left(x_1\right)=m_1{x}_1+$	$m_0=-0.1768$,
	$(m_0-m_1)$	$m_1=-1.1468$
	$\mathrm{sat}\left(x_1\right)$	[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$
NE6	${\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$
LE4	${\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
NE6	Base	1885
	SNNN	1885
	CNO	1885
LE4	Base	3383	1771
	SNNN	3383	1923
	CNO	3383	1771

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.

	${\mathit{x}}_1$	${\mathit{x}}_2$	${\mathit{x}}_3$	${\mathit{y}}_1$	${\mathit{y}}_2$	${\mathit{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

1. 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, pp. 2211-2219. [DOI: https://dx.doi.org/10.1007/s12555-018-0216-5]

2. 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, pp. 1175-1188. [DOI: https://dx.doi.org/10.1140/epjst/e2020-900225-9]

3. 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. Cybersecurity: A New Approach Using Chaotic Systems; Abd El-Latif, A.A.; Volos, C. Springer International Publishing: Cham, Switzerland, 2022; pp. 109-156. [DOI: https://dx.doi.org/10.1007/978-3-030-92166-8_6]

4. Bodale, I.; Oancea, V.A. Chaos control for Willamowski–Rössler model of chemical reactions. Chaos Solitons Fractals; 2015; 78, pp. 1-9. [DOI: https://dx.doi.org/10.1016/j.chaos.2015.06.019]

5. 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. [DOI: https://dx.doi.org/10.1016/j.chaos.2020.109756]

6. Zhao, X.; Li, Z.; Li, S. Synchronization of a chaotic finance system. Appl. Math. Comput.; 2011; 217, pp. 6031-6039. [DOI: https://dx.doi.org/10.1016/j.amc.2010.07.017]

7. 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. Nonlinear Dynamical Systems with Self-Excited and Hidden Attractors; Pham, V.T.; Vaidyanathan, S.; Volos, C.; Kapitaniak, T. Springer International Publishing: Cham, Switzerland, 2018; pp. 271-295. [DOI: https://dx.doi.org/10.1007/978-3-319-71243-7_12]

8. Zaqueros-Martinez, J.; Rodriguez-Gomez, G.; Tlelo-Cuautle, E.; Orihuela-Espina, F. Synchronization of Chaotic Electroencephalography (EEG) Signals. Cybersecurity: A New Approach Using Chaotic Systems; Abd El-Latif, A.A.; Volos, C. Springer International Publishing: Cham, Switzerland, 2022; pp. 83-108. [DOI: https://dx.doi.org/10.1007/978-3-030-92166-8_5]

9. 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. [DOI: https://dx.doi.org/10.1016/j.chaos.2022.112883]

10. 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. [DOI: https://dx.doi.org/10.3390/sym12091395]

11. Oluwole, O.S.A. Dynamic Regimes of El Niño Southern Oscillation and Influenza Pandemic Timing. Front. Public Health; 2017; 5, 301. [DOI: https://dx.doi.org/10.3389/fpubh.2017.00301]

12. Bendoukha, S.; Abdelmalek, S.; Ouannas, A. Secure Communication Systems Based on the Synchronization of Chaotic Systems. 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. [DOI: https://dx.doi.org/10.1007/978-3-030-12232-4_9]

13. Ç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, pp. 85-99. [DOI: https://dx.doi.org/10.1007/s10470-018-1252-z]

14. 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. 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. [DOI: https://dx.doi.org/10.1109/ICUMT.2015.7382410]

15. 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. [DOI: https://dx.doi.org/10.1063/1.4951991]

16. 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, pp. 513-523. [DOI: https://dx.doi.org/10.1080/17445760.2017.1334776]

17. Lambert, J.D. Computational Methods in Ordinary Differential Equations; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 1973.

18. Herinci, P. Discrete Variable Methods in Ordinary Differential Equations; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1962.

19. Pecora, L.M.; Carroll, T.L. Synchronization of chaotic systems. Chaos Interdiscip. J. Nonlinear Sci.; 2015; 25, 097611. [DOI: https://dx.doi.org/10.1063/1.4917383]

20. Volos, C.K.; Pham, V.T.; Azar, A.T.; Stouboulos, I.N.; Kyprianidis, I.M. Synchronization Phenomena in Coupled Dynamical Systems with Hidden Attractors. Nonlinear Dynamical Systems with Self-Excited and Hidden Attractors; Pham, V.T.; Vaidyanathan, S.; Volos, C.; Kapitaniak, T. Springer International Publishing: Cham, Switzerland, 2018; pp. 375-401. [DOI: https://dx.doi.org/10.1007/978-3-319-71243-7_17]

21. 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. [DOI: https://dx.doi.org/10.1155/2021/8865522]

22. 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, pp. 145-152. [DOI: https://dx.doi.org/10.3938/jkps.77.145]

23. Lai, Q.; Kamdem Kuate, P.D.; Pei, H.; Fotsin, H. Infinitely Many Coexisting Attractors in No-Equilibrium Chaotic System. Complexity; 2020; 2020, 8175639. [DOI: https://dx.doi.org/10.1155/2020/8175639]

24. Pm, G.; Kapitaniak, T. Synchronization in Coupled Multistable Systems with Hidden Attractors. Math. Probl. Eng.; 2017; 2017, 5214235. [DOI: https://dx.doi.org/10.1155/2017/5214235]

25. 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, pp. 814-833. [DOI: https://dx.doi.org/10.1016/j.jfranklin.2020.11.011]

26. 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, pp. 99-128. [DOI: https://dx.doi.org/10.24425/acs.2021.136882]

27. 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. [DOI: https://dx.doi.org/10.3390/e22030271]

28. 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, pp. 575-597. [DOI: https://dx.doi.org/10.24425/acs.2020.134678]

29. 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. [DOI: https://dx.doi.org/10.1155/2022/8362836]

30. Zhang, X.; Wang, C. Multiscroll Hyperchaotic System with Hidden Attractors and Its Circuit Implementation. Int. J. Bifurc. Chaos; 2019; 29, 1950117. [DOI: https://dx.doi.org/10.1142/S0218127419501177]

31. 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, pp. 35-42. [DOI: https://dx.doi.org/10.1080/00051144.2018.1491934]

32. 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, pp. 8851-8858. [DOI: https://dx.doi.org/10.1109/ACCESS.2017.2693301]

33. 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, pp. 2637-2655. [DOI: https://dx.doi.org/10.1007/s11071-017-3828-z]

34. 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. [DOI: https://dx.doi.org/10.1155/2017/7871467]

35. 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, pp. 2455-2466. [DOI: https://dx.doi.org/10.1140/epjs/s11734-022-00553-2]

36. 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. [DOI: https://dx.doi.org/10.3390/e23070921]

37. Zhang, H.; Liu, D.; Wang, Z. Controlling Chaos: Suppression, Synchronization and Chaotification; Springer: London, UK, 2009; [DOI: https://dx.doi.org/10.1007/978-1-84882-523-9]

38. 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, pp. 1021-1040. [DOI: https://dx.doi.org/10.1109/81.728857]

39. Skufca, J.D. Analysis Still Matters: A Surprising Instance of Failure of Runge–Kutta–Felberg ODE Solvers. SIAM Rev.; 2004; 46, pp. 729-737. [DOI: https://dx.doi.org/10.1137/S003614450342911X]

40. Gautschi, W. Numerical integration of ordinary differential equations based on trigonometric polynomials. Numer. Math.; 1961; 3, pp. 381-397. [DOI: https://dx.doi.org/10.1007/BF01386037]

41. Takagi, T.; Sugeno, M. Fuzzy identification of systems and its applications to modeling and control. IEEE Trans. Syst. Man Cybern.; 1985; SMC-15, pp. 116-132. [DOI: https://dx.doi.org/10.1109/TSMC.1985.6313399]

42. Lian, K.Y.; Chiu, C.S.; Chiang, T.S.; Liu, P. LMI-based fuzzy chaotic synchronization and communications. IEEE Trans. Fuzzy Syst.; 2001; 9, pp. 539-553. [DOI: https://dx.doi.org/10.1109/91.940967]

43. 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, pp. 66-83. [DOI: https://dx.doi.org/10.1109/3477.907565] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18244768]

44. Chua, L.O. A zoo of strange attractors from the canonical Chua’s circuits. Proceedings of the [1992] 35th Midwest Symposium on Circuits and Systems; Washington, DC, USA, 9–12 August 1992; Volume 2, pp. 916-926. [DOI: https://dx.doi.org/10.1109/MWSCAS.1992.271147]

45. 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, pp. 54-67.

46. Sprott, J.C. Some simple chaotic flows. Phys. Rev. E; 1994; 50, pp. R647-R650. [DOI: https://dx.doi.org/10.1103/PhysRevE.50.R647] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9962166]

47. Jafari, S.; Sprott, J.; Hashemi Golpayegani, S.M.R. Elementary quadratic chaotic flows with no equilibria. Phys. Lett. A; 2013; 377, pp. 699-702. [DOI: https://dx.doi.org/10.1016/j.physleta.2013.01.009]

48. Jafari, S.; Sprott, J. Simple chaotic flows with a line equilibrium. Chaos Solitons Fractals; 2013; 57, pp. 79-84. [DOI: https://dx.doi.org/10.1016/j.chaos.2013.08.018]

49. Chua, L.O. The Genesis of Chua’s Circuit. Arch. Fur Elektron. Und Uebertragungstechnik; 1992; 46, pp. 250-257.

50. Kuznetsov, N.; Leonov, G.; Vagaitsev, V. Analytical-numerical method for attractor localization of generalized Chua’s system. IFAC Proc. Vol.; 2010; 43, pp. 29-33. [DOI: https://dx.doi.org/10.3182/20100826-3-TR-4016.00009]

51. 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, pp. 5374-5379. [DOI: https://dx.doi.org/10.1016/S1474-6670(17)56915-3]

52. Nian, Y.; Zheng, Y. Generalized projective synchronization of chaotic systems based on Takagi-Sugeno fuzzy model. Proceedings of the 2010 Seventh International Conference on Fuzzy Systems and Knowledge Discovery; Yantai, China, 10–12 August 2010; Volume 3, pp. 1291-1295. [DOI: https://dx.doi.org/10.1109/FSKD.2010.5569112]

53. Danca, M.F. Hidden chaotic attractors in fractional-order systems. Nonlinear Dyn.; 2017; 89, [DOI: https://dx.doi.org/10.1007/s11071-017-3472-7]

54. 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. [DOI: https://dx.doi.org/10.1063/1.5038583]

55. 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. [DOI: https://dx.doi.org/10.1063/5.0081384]

56. Pham, V.T.; Volos, C.; Kapitaniak, T. Systems Without Equilibrium. Systems with Hidden Attractors: From Theory to Realization in Circuits; Springer International Publishing: Cham, Switzerland, 2017; pp. 51-63. [DOI: https://dx.doi.org/10.1007/978-3-319-53721-4_4]

57. Baranyi, P. The Generalized TP Model Transformation for T–S Fuzzy Model Manipulation and Generalized Stability Verification. IEEE Trans. Fuzzy Syst.; 2014; 22, pp. 934-948. [DOI: https://dx.doi.org/10.1109/TFUZZ.2013.2278982]

58. Baranyi, P. TP-Model Transformation-Based-Control Design Frameworks; Springer: Cham, Switzerland, 2016; [DOI: https://dx.doi.org/10.1007/978-3-319-19605-3]

59. 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, pp. 345-356. [DOI: https://dx.doi.org/10.1109/TFUZZ.2020.3038488]

60. Dudkowski, D.; Jafari, S.; Kapitaniak, T.; Kuznetsov, N.V.; Leonov, G.A.; Prasad, A. Hidden attractors in dynamical systems. Phys. Rep.; 2016; 637, pp. 1-50. [DOI: https://dx.doi.org/10.1016/j.physrep.2016.05.002]