Open AccessArticle

Hidden Homogeneous Extreme Multistability of a Fractional-Order Hyperchaotic Discrete-Time System: Chaos, Initial Offset Boosting, Amplitude Control, Control, and Synchronization

Amina-Aicha Khennaoui

¹,

Adel Ouannas

^2,3,

Stelios Bekiros

^4,5,6

Ayman A. Aly

⁷

Ahmed Alotaibi

⁷

Hadi Jahanshahi

^8,*

and

Hajid Alsubaie

⁷

Faculty of New Technologies of Information and Communication, University of Constantine 2, El Khroub 25016, Algeria

Department of Mathematics and Computer Science, University of Larbi Ben M’hidi, Oum El Bouaghi 04000, Algeria

Nonlinear Dynamics Research Center (NDRC), Ajman University, Ajman 346, United Arab Emirates

⁴

Faculty of Economics Management & Accountancy, University of Malta, MSD 2080 Msida, Malta

⁵

LSE Health, Department of Health Policy, London School of Economics and Political Science, London WC2A 2AE, UK

⁶

IPAG Business School, 184 Boulevard Saint-Germain, 75006 Paris, France

⁷

Department of Mechanical Engineering, College of Engineering, Taif University, Taif 21944, Saudi Arabia

⁸

Department of Mechanical Engineering, University of Manitoba, Winnipeg, MB R3T 5V6, Canada

Author to whom correspondence should be addressed.

Symmetry 2023, 15(1), 139; https://doi.org/10.3390/sym15010139

Submission received: 18 October 2022 / Revised: 11 December 2022 / Accepted: 26 December 2022 / Published: 3 January 2023

(This article belongs to the Special Issue Symmetry: Recent Developments in Engineering Science and Applications)

Download

Browse Figures

Versions Notes

Abstract

Fractional order maps are a hot research topic; many new mathematical models are suitable for developing new applications in different areas of science and engineering. In this paper, a new class of a 2D fractional hyperchaotic map is introduced using the Caputo-like difference operator. The hyperchaotic map has no equilibrium and lines of equilibrium points, depending on the values of the system parameters. All of the chaotic attractors generated by the proposed fractional map are hidden. The system dynamics are analyzed via bifurcation diagrams, Lyapunov exponents, and phase portraits for different values of the fractional order. The results show that the fractional map has rich dynamical behavior, including hidden homogeneous multistability and offset boosting. The paper also illustrates a novel theorem, which assures that two hyperchaotic fractional discrete systems achieve synchronized dynamics using very simple linear control laws. Finally, the chaotic dynamics of the proposed system are stabilized at the origin via a suitable controller.

Keywords:

fractional map; complexity; initial boosting attractors; amplitude control; chaos; synchronization

1. Introduction

While discrete fractional calculus was first introduced roughly fifty years ago, the non-integer-order derivative has been a mathematical idea since 1695 [1]. In particular, the discretization of continuous-time operators [2] in 1974 led to the first derivations of fractional differential equations. Due to the recent widespread interest in discrete fractional calculus, many difference operators, notably the Grünwald–Letnikov difference operator [3], have been proposed in the literature. Many efforts have been made to thoroughly analyze the dynamics of both classical systems and fractional systems since the discovery of chaos phenomena in dynamical systems [3,4]. Regarding the latter, various works have been conducted on the topic of studying chaotic behaviors in nonlinear maps defined by fractional order difference equations [5,6,7]. Wu et al. [8] proposed a fractional logistic map and studied the chaotic behavior based on the Caputo-left difference operator. One-dimensional fractional maps with chaotic attractors and quasi-periodic behaviors were suggested in [9]. Further evidence for chaotic attractors was provided by Khennaoui et al. [10] in three fractional maps, namely the Wang map, the Stefanski map, and the Rössler map. The dynamic behavior of a higher-dimensional fractional-order chaotic map was described by Peng et al. in [11] that same year, while Wang et al. examined the chaos and complexity of a fractional order higher-dimensional multi-cavity chaotic map in [12]. Memristive maps with integer and fractional orders have recently been proposed [13,14]. A memristor-based Rulkov neuron map was introduced in [15] by using a discrete memristor to replicate the magnetic effects, whereas [16] illustrated the dynamic features of a novel discrete memristive hyperchaotic map.

Compared with integer-order chaotic maps, fractional-order chaotic maps have more complex dynamic characteristics, i.e., where a parameter can increase the degree of freedom of the model. These additional degrees of freedom can also be employed to uncover the subtleties of ecological processes that occur in the real world. In order to better understand these complex dynamical characteristics, more discrete maps with fractional operators must be provided. Furthermore, one of the distinctive benefits of fractional-order maps is that they are sensitive to variations in fractional orders in addition to tiny perturbations in parameters and initial circumstances. Additionally, fractional-order discrete maps are advantageous for system analysis and numerical computation because they have straightforward forms and rich dynamics.

In recent years, coexisting attractors and hidden attractors have been reported in various chaotic fractional maps. For example, in [17], the presence of chaos in fractional discrete systems with rectangle-shaped and square-shaped fixed points was investigated; ref. [18] provided an illustration of the behavior of a discrete fractional-order Bonhoeffer-van der Pol oscillator. The coexistence of various attractors, such as chaotic and hyperchaotic attractors as well as quasi-periodic and periodic attractors, is highlighted by the map. Furthermore, Almatroud et al. [19] showed that a brand-new 2D fractional-order chaotic map included hidden extreme multistability. A fractional discrete system without fixed points is provided in [20]. To demonstrate the existence of chaotic hidden attractors in the system dynamics, the 0–1 test and the computation of the approximation entropy have been used. Later, using bifurcation diagrams, Lyapunov exponents, and a 0–1 test, the complex dynamic behaviors of a fractional map with hidden attractors were addressed; see [21]. The synchronization and control of the system dynamics are also suggested. A synchronized chaotic system aims to achieve the state of a chaotic master system. Hence,

t \to + \infty

gives us a synchronization error that converges toward zero. The design of synchronization schemes is influenced by a wide range of techniques. The chaos synchronization phenomenon was studied in many different ways, such as hybrid projection synchronization, inverse matrix projection synchronization, Q–S synchronization, complete synchronization, etc. [22,23,24].

In this study, a discrete-time new fractional system with the line of equilibrium and no equilibrium points was constructed with the following properties: homogeneous multistability, hidden dynamics, and amplitude control. This new fractional-order hyperchaotic map possesses not only hidden attractors but also coexisting attractors. The initial boosting phenomena were investigated in the proposed hyperchaotic map, as seen in Section 3, with an independent knob for partial amplitude control in Section 4. A suitable controller stabilizes the hyperchaotic system in Section 5. A novel theorem is presented in Section 6, which claims that two fractional hyperchaotic systems can be controlled linearly to show synchronized dynamics. Our conclusions with ideas for future work are presented in Section 7.

2. Description and Analysis of the Fractional Hyperchaotic Map

Recently, Bao et al. [25] proposed a new hyperchaotic map by introducing the nonlinear sinus function into a simple map, which is governed by Equation (1). As shown below, it has only one nonlinear term and two real constants controllers A and B:

\{\begin{matrix} x_{1} (n + 1) = x_{1} (n) + A x_{1} (n) sin x_{2} (n) + B, \\ x_{2} (n + 1) = x_{1} (n) + x_{2} (n) . \end{matrix}

(1)

Inspired by some relevant works on fractional-order discrete-time system(s) (FODTS) [26], we define the fractional-order case by introducing the

η

-Caputo difference operator given by:

\begin{matrix} ^{C} Δ_{a}^{η} g (s) = \frac{1}{η (m - η)} \sum_{τ = a}^{s - (m - η)} {(s - τ - 1)}^{(m - η - 1)} Δ_{τ}^{m} g (τ), \end{matrix}

(2)

where

s \in N_{a + m - η}

m = ⌈η⌉ + 1

and

η \notin N

. The corresponding fractional-order form is built as follows

\{\begin{matrix} ^{C} Δ_{a}^{η} x_{1} (s) = A x_{1} (s - 1 + η) sin x_{2} (s - 1 + η) + B, \\ ^{C} Δ_{a}^{η} x_{2} (s) = x_{1} (s - 1 + η), \end{matrix}

(3)

where

s \in N_{a + 1 - η}

and

η \in (0, 1]

are the fractional orders. The 2D FODTS (3) is symmetric with respect to the plane

S = {(x_{1}, x_{2}) \in R, x_{1} = 0, x_{2} = 0}

, and control parameters A and B, where it is always invariant under the transformation

(x_{1}, x_{2}, A, B) \to (- x_{1}, - x_{2}, - A, - B)

. Thus, the attractor of the proposed FODTS could display symmetric behavior.

The equilibrium points of the FODTS (3) are easily determined from:

\{\begin{matrix} A x_{1} sin x_{2} + B = 0, \\ x_{1} = 0 . \end{matrix}

(4)

Similar to the integer-order counterpart, two different cases are considered:

Case a: when

B \neq 0

Equation (4) has no solution, demonstrating that the fractional map (3) has no equilibrium point.

Case b: when

B = 0

, the fixed points of the FODTS (3) are given by

L = (0, C);

where C is a real number. In other words, the proposed system has line equilibrium points. The Jacobian matrix

(J a c)

for system (3), computed at L, is:

J a c = [\begin{matrix} A sin C & 0 \\ 1 & 0 \end{matrix}],

(5)

whose determinant is equal to zero

(det J a c = 0)

. According to the stability theorem of a fractional-order discrete-time system [19], the line of the equilibrium point is stable if

\{\begin{matrix} \frac{- T r (J a c)}{2} \geq \sqrt{D e t (J a c)}, \\ η > l o g_{2} \frac{\sqrt{T r {(J a c)}^{2} - 4 D e t (J a c)} - T r (J a c)}{2}, \\ D e t (J a c) > 0 \end{matrix}

(6)

The conditions

\frac{- T r (J a c)}{2} \geq \sqrt{D e t (J a c)}

and

D e t (J a c) > 0

cannot be satisfied at the same time; thus, system (3) has a high possibility of hidden attractors.

By definition, all chaotic attractors generated by the proposed FODTS (4) are hidden attractors.

3. Dynamical Analysis and Numerical Simulations

Nonlinear dynamical analysis tools, such as phase diagrams, bifurcation diagrams, and Lyapunov exponents are used to analyze the hidden dynamics. We shall first give the numerical formula of the proposed map. By applying the

η

fractional sum,

Δ_{τ}^{- η} g (τ) = \frac{1}{Γ (η)} \sum_{τ = a}^{s - γ} {(s - τ - 1)}^{(γ - 1)} g (τ)

, the numerical formula is designed as:

\{\begin{matrix} {x_{1}}_{n} = {x_{1}}_{0} + \frac{1}{Γ (η)} \sum_{j = 0}^{n - 1} \frac{Γ (n - 1 - j + η)}{Γ (n - j)} A {x_{1}}_{j} sin {x_{2}}_{j} + B, \\ {x_{2}}_{n} = {x_{2}}_{0} + \frac{1}{Γ (η)} \sum_{j = 0}^{n - 1} \frac{Γ (n - 1 - j + η)}{Γ (n - j)} {x_{1}}_{j} . \end{matrix}

(7)

As one can see from Equation (7), the proposed FODTS (3) has the memory effect. Hence, the states

{x_{1}}_{n}, {x_{2}}_{n}

depend on

x_{0}, x_{1}, x_{2}, \dots, x_{n - 1}

Using numerical simulations of the same bifurcation parameters for the integer-order map reported in [25], we present the dynamics of the suggested fractional-order map in order to find hidden chaotic attractors.

3.1. Hidden Attractors and Bifurcation Analysis

The control parameter B is assigned as

B = \pm 0.1

B = 0

, respectively, and the initial condition is selected as

{x_{1}}_{0} = \pm 1, {x_{2}}_{0} = \pm 2

, respectively. The properties of the coexisting symmetrical hidden attractors are illustrated by considering the bifurcation diagrams and Lyapunov exponents (LEs) for

A \in [2.3, 2.9]

and

A \in [- 2.9, - 2.3]

, respectively. The bifurcation diagram is obtained by plotting the local maxima of the state

x_{2}

in terms of the control parameter A, whereas the LEs are computed numerically using the Jacobian matrix algorithm. The LEs are computed numerically using the Jacobian matrix algorithm for fractional maps [27].

Case A: No Fixed Point

B = 0.1

and the initial values

{x_{1}}_{0} = 1

{x_{2}}_{0} = - 2

; the bifurcation diagram and the corresponding LEs are plotted in Figure 1a. These diagrams are obtained for

η = 0.98

and for

A \in [2.3, 2.9]

. From the LEs and bifurcation diagram in Figure 1a, it can be seen that the FODTS (3) gradually changes from the periodic state to a chaotic state at

A \in [2.665, 2.692] \cup [2.704, 2.85]

, to a hyperchaotic state for

A \in [0.8, 0.9286]

, through period-doubling. With the change of system parameter A, the LE diagram approves the dynamic behavior observed in the bifurcation diagram.

In contrast, when the control parameter B is fixed as

B = - 0.1

and the initial values are set to

{x_{1}}_{0} = - 1

{x_{2}}_{0} = 2

, while the control parameter A changes in the range

[- 2.9, - 2.3]

; the FODTS (3) also shows the period-doubling route to chaos. Namely, the FODTS (3) is chaotic when

A \in [2.665, 2.692] \cup [2.704, 2.85]

, where

m a x λ

is positive, and is in hyperchaos when

A \in [0.8, 0.9286]

; where

λ_{1}, λ_{2} > 0

. Basically, Figure 1 shows that under symmetrical initial values, the dynamic evolution of the FODTS caused by the symmetrical control parameters

A \in [2.3, 2.9]

and

A \in [- 2.9, - 2.3]

, respectively, had symmetrical similarity.

To observe the coexistence of hidden symmetric attractors in the FODTS (3), the phase portraits, as seen in Figure 1, are disposed of in Figure 2. Red trajectories were obtained for the constant parameter

B = 0.1

, fractional order

η = 0.98

, and the initial values

{x_{1}}_{0} = 1, {x_{2}}_{0} = - 2

; blue trajectories were obtained for

B = - 0.1

η = 0.98

{x_{1}}_{0} = - 1, {x_{2}}_{0} = 2

. For different values of A, the FODTS shows symmetrical hidden periodic attractors, symmetrical hidden chaotic attractors, and symmetrical hidden hyperchaotic attractors. Therefore, it is concluded that the FODTS has complex hidden dynamics. Moreover, the coexistence behavior indicates the sensitivity of the map (3) to the control parameters and initial values.

Case B: Line of the Equilibrium Point

As the system parameters and the initial state variables of the 2D fractional system are taken as

B = 0

{x_{1}}_{0} = 1, {x_{2}}_{0} = - 2

, the bifurcation diagrams along with the LEs diagram for the fractional order value

η = 0.95

are shown in Figure 3a. As illustrated in the previous section, when

B = 0

, the FODTS (3) has a line of equilibria. When the bifurcation A increases from 2.3 to 2.9, it can be seen from the bifurcation diagram that the FODTS starts from the period and goes into chaos via the cascade of period-doubling bifurcation and then turns into hyperchaos at interval

A \in [2.689, 2.734] \cup [2.77, 2.813]

. In contrast, when the parameter A is adjusted in the region

[- 2.9, - 2.3]

, the bifurcation of the state

x_{2}

and the LEs are numerically simulated as shown in Figure 3b. Similarly, it is easily viewed that the bifurcation diagram in Figure 3b had certain symmetrical similarities with the one reported in Figure 3a. Therefore, it proves that the coexistence behavior appears in the hidden attractors of the FODTS with the lines of the equilibrium points.

Without loss of generality, the hidden attractors at several critical values corresponding to Figure 3 are shown in Figure 4.

3.2. The Effect of Fractional Order $η$

In order to further understand the dynamics of the FODTS, the fractional order

η

is considered a bifurcation parameter, and the bifurcation diagram along with the LEs are derived as shown in Figure 5. Note that the blue diagram is obtained for the initial condition

(1, - 2)

and system parameters

A = 2.7, B = 0.1

, while the red diagram is obtained for

A = - 2.7

B = - 0.1

, with initial conditions

{x_{1}}_{0} = - 1

{x_{2}}_{0} = 2

. It is observed that the new FODTS produces more complex dynamics compared with the integer-order map (1). In particular, for the corresponding integer-order value

η = 1

, the FODTS is chaotic where the maximum LEs is bigger than 0. However, it becomes hyperchaotic as the fractional order value decreases, where both of the Lyapunov exponents

(L E_{1}, L E_{2})

become bigger than 0, which indicates that the dynamic characteristic of the FODTS is more complex. Moreover, when

η \in [0.8, 0.9286]

, the FODTS (3) is in hyperchaos, where

λ_{1}, λ_{2} > 0

with some unbounded region for the largest values of n. When

η \in [0.9286, 0.9739] \cup [0.9869, 1]

the FODTS is in chaos, where

m a x λ > 0

. The phase portraits for the fractional order values

η = 0.9

and

η = 1

, which are associated with the hidden hyperchaotic and hidden chaotic attractors, respectively, are illustrated at the top of Figure 5.

The bifurcation diagram and LEs are used to analyze the relationship between the control parameter A and the fractional order

η

. When the system parameter A is adjusted in the region

[2.2, 2.9]

, the fractional order

η

, respectively, is selected as

η = 0.9

η = 0.95

η = 1

. Three kinds of bifurcation diagrams are depicted in Figure 6, with different color regions. Correspondingly, the LEs are calculated and shown in Figure 6. The dynamic behavior of the FODTS depends on the value of the control parameter A and the value of

η

. Figure 6 shows that the states of the FODTS (3) go from a periodic to a chaotic state with the decrease of order

η

, as A increases and

η

decreases, it is observed that the chaotic band suddenly disappears and the system goes into infinity. The LE diagrams confirm that the system passes from periodic regions to chaotic and hyperchaotic regions as the system parameter A increases, and the fractional order value decreases.

3.3. Hidden Extreme Homogeneous Multistability

To reveal the extreme multistability of the FODTS (3), the maximum LEs and the bifurcation diagrams of the state variable

{x_{2}}_{n}

are calculated as shown in Figure 7, where the fractional order

η

is chosen as

0.98

. The parameters are selected as

A = 2.7

B = 0.1

, and the IV is fixed as

{x_{1}}_{0} = 1

, while

{x_{2}}_{0}

varies in the range

[0, 200]

. As shown in Figure 7, when

η = 0.98

, system (3) generates multiple chaotic attractors along the

x_{2}

axis, where the values of the LEs of the FODTS (3) are nearly the same, which indicates that all of these attractors have very close chaotic features. To exhibit the phenomenon of multistability, typical hyperchaotic attractors of the FODTS (3) with no equilibrium for a different IV

{x_{2}}_{0}

are simulated, as shown in Figure 8. As can be seen, different initial condition values

{x_{2}}_{0}

lead to different hidden hyperchaotic attractors of the same shape. This special phenomenon is known as homogeneous multistability, which is a new phenomenon in nonlinear systems [28]. The hidden homogeneous multistability of the FODTS (3) (with respect to small changes of the initial state

{x_{2}}_{0}

) is illustrated. Note that there are many other hidden hyperchaotic attractors with unified LEs for other IVs.

3.4. Initial Offset Boosting

In the following, we will discuss an unusual phenomenon of the new FODTS, namely the initial offset boosting by changing the IV

{x_{2}}_{0}

with the

2 π

period as

x_{2} 0 = - 2 + 2 K π

with the fixed order

η = 0.95

, where

K = - 1, 0, 1, 2

and

{x_{1}}_{0} = 1

. In Figure 9, we present the bifurcation diagrams of system (3) when

A \in [2.3, 2.9]

and

B = 0

. For four different values of IC

{x_{2}}_{0}

, the FODT systems experience four bifurcation diagrams with periods

2 π

. All of these diagrams have the same structure with the

2 π

regime of homogenous multistability. The same result can be obtained if we consider the second case

B \neq 0.1

Fix

A = 2.8

B = 0.1

, and let the fractional order

η = 0.95

; the hidden chaotic attractors are illustrated in Figure 10 for

{x_{2}}_{0} = - 2 + 2 K

. These hidden chaotic attractors have the same shape as regime

2 π

4. Amplitude Control Analysis

In this section, the amplitude control method is used to detect coexisting attractors of the FODTS (3). Amplitude control has potential application value in engineering applications because it can achieve stability and explain the unpredictable behaviors that occur in engineered systems [29,30,31].

It is very interesting that the proposed FODTS possesses the property of amplitude control, which means that the amplitudes of one or both states of the map are adjusted. By introducing a controlled parameter

c_{1}

into the variable

x_{2}

, respectively, the FODTS (3) can be changed to

\{\begin{matrix} ^{C} Δ_{a}^{η} x_{1} (s) = A x_{1} (s - 1 + η) sin (x_{2} (s - 1 + η)) + B, \\ ^{C} Δ_{a}^{η} x_{2} (s) = \frac{x_{1} (s - 1 + η)}{c_{1}}, \end{matrix}

(8)

which is identical to Equation (3) with transformation

x_{1} = u

and

x_{2} = c_{1} v

. Therefore, parameter

c_{1}

proportionally controls the amplitude of the FODTS states

x_{1}

and

x_{2}

according to

c_{1}

. In order to discuss the new fractional discrete-time system (8) in terms of controller

c_{1}

, we fix

A = 2.8

B = 0.1

η = 0.98

, and

{x_{1}}_{0} = 1; {x_{2}}_{0} = - 2

; while parameter

c_{1}

takes the values

{- 5, 2, 12, 13}

. The phase diagrams are depicted in Figure 11. These diagrams have the same shape structure as regime

c_{1}

. In particular, from Figure 11, it can be seen that the chaotic attractor of system (8) is larger as

c_{1}

increases. Nevertheless, there are other cases of coexisting attractors that depend on the control parameter

c_{1}

5. Complexity Analysis of the FODTS

In this section, the complexity of the FODTS (3) is analyzed using the approximate entropy (ApEn) [20]. An approximate entropy indicates how complex a system generated by a time series is. In general, a time series with the largest values of ApEn is considered as being more complex. ApEn can be calculated for both continuous and discrete systems and integer and fractional orders. Firstly, the calculating process of the

A p E n

algorithm is presented [32]. For the specific calculating method, we consider N points of the state

{(x_{1})}_{i = 1, . ., N}

; the specific steps are as follows:

Step 1. Construct a sequence of m vectors. For a given time series ${(x_{1})}_{i = 1, . ., N}$ , the m vector sequence ${X_{m} (1), \dots, X_{i + m - 1}}$ is constructed as

$X_{m} (1) = {x (i), x (i + 1), \dots, x (i + m - 1)}, 1 \leq i \leq B - m + 1 .$
Step 2. For each $1 \leq i \leq N - m + 1$ , define the following equation

$C_{r}^{m} (i) = \frac{n b r o f d [X (i), X (j)] \leq r}{N - m + 1},$

where $d [X (i), X (j)]$ is the distance between $X (i)$ and $X (j)$ given by

$d [X (i), X (j)] = m a x_{k = 0, \dots, m} |X (i + k - 1) - X (j + k - 1)|,$

and $r = 0.2 s t d (x)$ in which $s t d (x)$ presents the standard deviation of the data x.
Step 3. On the basis of $C_{r}^{m}$ , the average value is denoted to be

$ϕ^{m} (r) = \frac{1}{n - m - 1} \sum_{i = 1}^{n - m + 1} log C_{i}^{m} (r) .$

(9)
Step 4. The $A p E n$ is calculated as follows

$A p E n = ϕ^{m} (r) - ϕ^{m + 1} (r) .$

(10)

The

A p E n

analysis results of the proposed FODT system (3) with varying

η

and A are shown in Figure 12. It can be seen that the proposed system (3) can have a higher complexity with relatively larger parameter values a and fractional order

η

. The analysis results in Figure 12 indicate that the

A p E n

test cannot distinguish chaos and hyperchaos but it can be used as the parameter choice in the practical application.

6. Control Laws

This section deals with the stability analysis of hyperchaotic FODTS (3). It is our first objective to present the theorem below, which is used as the basis for analyzing the stabilization of fractional discrete systems. Next, we perform numerical simulations to demonstrate the effectiveness of the proposed control method.

Theorem 1

([33]). Let

x = 0

be an equilibrium point for the h-fractional discrete system

_{h}^{C} Δ_{a}^{η} x (s) = f (s + η h, x (s + η h)), s {(h N)}_{a + (1 - η) h} .

(11)

Assume that there exists a positive definite and decrescent scalar function

V (s, x (s))

, such that

_{h}^{C} Δ_{a}^{η} V (s, x (s)) \leq 0

Then the zero equilibrium point of the nonlinear h-fractional discrete system is asymptotically stable.

Lemma 1

([33]). For any discrete time

s \in {(h N)}_{a + (1 - η) h}

, we have

_{h}^{C} Δ_{a}^{η} x^{2} (s) \leq 2 x {(s + η h)}_{h}^{C} Δ_{a}^{η} x (s), 0 < η \leq 1 .

(12)

To obtain our results, the following theorem is presented.

Theorem 2.

The two-dimensional FODTS (3) is controlled under the following one-dimensional control law

C (s) = - (|A| + 1) x_{1} (s) - x_{2} (s) - B,

(13)

where

s \in {(h N)}_{a + (1 - η) h}

Proof.

The controlled FODTS can be described as

\{\begin{matrix} ^{C} Δ_{a}^{η} x_{1} (s) = A x_{1} (s - η + 1) sin (x_{2} (s - η + 1)) + B + C (s - η + 1), \\ ^{C} Δ_{a}^{η} x_{2} (s) = x_{1} (s - η + 1), \end{matrix}

(14)

Consequently, (14) takes the form:

\{\begin{matrix} ^{C} Δ_{a}^{η} x_{1} (s) = A x_{1} (s - η + 1) sin (x_{2} (s - η + 1)) - (|A| + 1) x_{1} (s - η + 1) - x_{2} (s - η + 1), \\ ^{C} Δ_{a}^{η} x_{2} (s) = x_{1} (s - η + 1), \end{matrix}

(15)

One can utilize the Lyapunov method by first letting the Lyapunov function,

V (s)

, in the form:

V = \frac{1}{2} x_{1}^{2} (s) + \frac{1}{2} x_{1}^{2} (s),

(16)

the adoption of the Caputo h-difference operator implies that

_{h}^{C} Δ_{a}^{η} V (s) = \frac{1}{2}_{h}^{C} Δ_{a}^{η} x_{1}^{2} (s) + \frac{1}{2}_{h}^{C} Δ_{a}^{η} x_{2}^{2} (s) .

(17)

By using Lemma 1, it follows that

\begin{matrix} _{h}^{C} Δ_{a}^{η} V & \leq & x_{1} {(s + η h)}_{h}^{C} Δ_{a}^{η} x_{1} (s) + x_{2} {(s + η h)}_{h}^{C} Δ_{a}^{η} x_{2} (s) \\ = & x_{1} (s - η + 1) [A x_{1} (s - η + 1) sin (x_{2} (s - η + 1)) - (|A| + 1) x_{1} (s - η + 1) - x_{2} (s - η + 1)] \\ x_{2} (s - η + 1) x_{1} (s - η + 1) \\ = & A x_{1}^{2} (s - η + 1) sin (x_{2} (s - η + 1)) - (|A| + 1) x_{1}^{2} (s - η + 1) \\ \leq & |A x_{1}^{2} (s - η + 1) sin (x_{2} (s - η + 1))| - (|A| + 1) x_{1}^{2} (s - η + 1) \\ \leq & |A| x_{1}^{2} (s - η + 1) |sin (x_{2} (s - η + 1))| - (|A| + 1) x^{2} (s - η + 1) \\ \leq & |A| x_{1}^{2} (s - η + 1) - (|A| + 1) x_{1}^{2} (s - η + 1) \\ = & - x_{1}^{2} (s - η + 1) \leq 0 \end{matrix}

Theorem 1 states that (14) has an asymptotically stable zero equilibrium. This implies that the chaotic trajectories of this proposed FODTS (3) are stabilized by the linear control law (13). □

As shown in Section 3, the FODTS (3) with

A = 2.8, B = 0.1

, and

η = 0.98

display chaotic behavior. Figure 13 depicts the evolution of states and the phase space plot of the controlled FODTS (3), respectively. As one can see, the FODTS (3) is completely controlled in a stable state. These plots confirm that the hidden chaotic dynamic of the FODTS is controlled effectively.

7. Synchronization

In this section, we will investigate the synchronization of the hyperchaotic FODTS (3). We present a new theorem that ensures that the dynamics of two chaotic FODT systems are synchronized through linear control laws of a very simple nature. The master discrete system is defined as

\{\begin{matrix} ^{C} Δ_{a}^{η} x_{1 m} (s) = A x_{1 m} (s - η + 1) sin (x_{2 m} (s - η + 1)) + B, \\ ^{C} Δ_{a}^{η} x_{2 m} (s) = x_{1 m} (s - η + 1), \end{matrix}

(18)

and the slave system as

\{\begin{matrix} ^{C} Δ_{a}^{η} x_{1 s} (s) = A x_{1 s} (s - η + 1) sin (x_{2 s} (s - η + 1)) + B + L_{1}, \\ ^{C} Δ_{a}^{η} x_{2 s} (s) = x_{1 s} (s - η + 1) + L_{2}, \end{matrix}

(19)

where

L_{1}

and

L_{2}

are synchronization controllers that have to be designed. Using the master (18) and slave discrete systems (19), the error dynamics can be derived as

\{\begin{matrix} e_{1} = x_{1 s} - x_{1 m}, \\ e_{2} = x_{2 s} - x_{2 m}, \end{matrix}

(20)

lim_{s \to + \infty} | e_{i} (s) | = 0

, then the synchronization scheme is achieved.

The following theorem, which we consider as being one of the most important outcomes of this research, is intended to highlight the important results obtained by the proposed synchronization scheme.

Theorem 3.

The master system (18) and the slave system (19) achieve synchronized dynamics, provided that the control law is selected as

\{\begin{matrix} L_{1} (s) = - (|A| + \frac{{(A l)}^{2}}{4}) e_{1}, \\ L_{2} (s) = - e_{1} - e_{2}, \end{matrix}

(21)

where

s \in {(h N)}_{a + (1 - η) h} .

Proof.

To establish an asymptotic convergence of the synchronization errors, given in (20), we start applying the Caputo-type fractional-order differences on (20), which yields:

\{\begin{matrix} _{h}^{C} Δ_{a}^{η} e_{1} = A [x_{1 s} (s - η + 1) sin (x_{2 s} (s - η + 1)) - x_{1 s} (s - η + 1) sin (x_{2 s} (s - η + 1))] + L_{1}, \\ _{h}^{C} Δ_{a}^{η} e_{2} = x_{1 s} (s - η + 1) - x_{1 m} (s - η + 1) + L_{2}, \end{matrix}

(22)

(21) is subsumed into (22) to yield the following new discrete system:

\{\begin{matrix} _{h}^{C} Δ_{a}^{η} e_{1} = A [x_{1 s} (s - η + 1) sin (x_{s} (s - η + 1)) - x_{1 s} (s - η + 1) sin (x_{2 s} (s - η + 1))] - (|A| + \frac{{(A l)}^{2}}{4}) e_{1}, \\ _{h}^{C} Δ_{a}^{η} e_{2} = - e_{2}, \end{matrix}

(23)

Let

V = \frac{1}{2} e_{1}^{2} (s) + \frac{1}{2} e_{2}^{2} (s) .

This implies

^{C} Δ_{a}^{η} V =

^{C} Δ_{a}^{η} e_{1}^{2} (s) +

^{C} Δ_{a}^{η} e_{2}^{2} (s),

and by using Lemma 1, we obtain

\begin{matrix} ^{C} Δ_{a}^{η} V & \leq & e_{1} {(s - η + 1)}^{C} Δ_{a}^{η} e_{1} (s - η + 1) + e_{2} {(s - η + 1)}^{C} Δ_{a}^{η} e_{2} (s - η + 1) \\ = & a e_{1} (t - η + 1) [x_{1 s} (s - η + 1) sin (x_{2 s} (s - η + 1)) - x_{1 s} (s - η + 1) sin (x_{2 s} (s - η + 1))] \\ - (|A| + \frac{{(A l)}^{2}}{4}) e_{1}^{2} - e_{2}^{2} (s - η + 1) \end{matrix}

one can write:

\begin{matrix} x_{1 s} (s - η + 1) sin (x_{2 s} (s - η + 1)) - x_{1 s} (s - η + 1) sin (x_{2 s} (s - η + 1)) \\ = x_{1 s} (s - η + 1) (sin (x_{2 s} (s - η + 1)) - sin (x_{2 m} (s - η + 1))) \\ - (x_{1 s} (s - η + 1) - x_{1 m} (s - η + 1)) sin (x_{2 m} (s - η + 1)) \end{matrix}

So, we have

\begin{matrix} ^{C} Δ_{a}^{η} V & \leq & |A| |e_{1} (t - η + 1)| |x_{1 s} (s - η + 1)| |(sin (x_{2 s} (s - η + 1)) - sin (x_{2 m} (s - η + 1)))| \\ + |A| e_{1}^{2} (s - η + 1) |sin (x_{2 m} (s - η + 1))| - (|A| + \frac{{(A)}^{2}}{4}) e_{1}^{2} - e_{2}^{2} (s - η + 1) \end{matrix}

since fractional chaotic maps have the ‘property of boundedness’, we can put

|x_{1 s} (s - η + 1)| \leq l

and

|(sin (x_{2 s} (s - η + 1)) - sin (x_{2 m} (s - η + 1)))| \leq |x_{2 s} (s - η + 1) - x_{2 m} (s - η + 1)|

so we have

\begin{matrix} ^{C} Δ_{a}^{η} V & \leq & |A| l |e_{1} (s - η + 1)| |(x_{2 s} (s - η + 1) - x_{2 m} (s - η + 1))| + |A| e_{1}^{2} (s - η + 1) \\ - (|A| + \frac{{(A l)}^{2}}{4}) e_{1}^{2} - e_{2}^{2} (s - η + 1) \\ = & |A| l |e_{1} (s - η + 1)| |e_{2} (s - η + 1)| - \frac{{(A)}^{2}}{4} e_{1}^{2} - e_{2}^{2} (s - η + 1) \\ = & - {(\frac{|A| l}{2} e_{1} - e_{2})}^{2} < 0 . \end{matrix}

According to the stability theorem (Theorem 1), error system (22) is stabilized to the origin. This means that the fractional slave discrete system (19) can asymptotically synchronize the fractional master discrete system (18). □

The synchronization errors between the hyperchaotic FODT master and slave systems are depicted in Figure 14 for the initial values

e_{1} (0) = 1, e_{2} (0) = - 2,

and

η = 0.98

. There is a convergence of errors into zero, which indicates that chaos synchronization has indeed been achieved.

8. Conclusions

A 2D fractional map with hidden multistability is presented using the Caputo-like difference operator. Bifurcation diagrams, Lyapunov exponents, and phase portraits have been used to study the dynamics of the system in detail. Multiple hidden attractors have been revealed in the new fractional map. In order to stabilize the chaotic trajectories of the fractional discrete-time system, a new theorem has been illustrated. A new theorem assures that two hyperchaotic fractional discrete-time systems can achieve synchronized dynamics via very simple linear control laws. Future work will involve leveraging the hardware implementation developed herein in order to develop cryptographic applications based on fractional maps with hidden chaotic attractors. It is being investigated how fractional maps with hidden attractors can be used to develop new types of fractional pseudo-number generators.

Author Contributions

Conceptualization, A.-A.K., A.O., S.B., A.A.A., A.A., H.J., and H.A.; methodology, A.-A.K., A.O., S.B., A.A.A., A.A., H.J., and H.A.; software, A.-A.K., A.O., S.B., A.A.A., A.A., H.J., and H.A.; validation, A.-A.K., A.O., S.B., A.A.A., A.A., H.J., and H.A.; writing—original draft preparation, A.-A.K., A.O., S.B., A.A.A., A.A., H.J., and H.A.; writing—review and editing, A.-A.K., A.O., S.B., A.A.A., A.A., H.J., and H.A.; supervision, A.-A.K., A.O., S.B., A.A.A., A.A., H.J., and H.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Taif University Researchers Supporting Project number (TURSP-2020/77), Taif University, Taif, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

Herrmann, R. Fractional Calculus—An Introduction for Physicists; World Scientific: Singapore, 2018. [Google Scholar]
Diaz, J.B.; Olser, T.J. Differences of fractional order. Math. Comput. 1974, 28, 185–202. [Google Scholar] [CrossRef] [Green Version]
Goodrich, C.; Peterson, A.C. Discrete Fractional Calculus; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
Zaslavsky, G.M.; Zaslavskij, G.M. Hamiltonian Chaos and Fractional Dynamics; Oxford University Press: Oxford, UK, 2005. [Google Scholar]
Wang, Y.; Liu, S.; Li, H. On fractional difference logistic maps: Dynamic analysis and synchronous control. Nonlinear Dyn. 2020, 102, 579–588. [Google Scholar] [CrossRef]
Kassim, S.; Hamiche, H.; Djennoune, S.; Bettayeb, M. A novel secure image transmission scheme based on synchronization of fractional-order discrete-time hyperchaotic systems. Nonlinear Dyn. 2017, 88, 2473–2489. [Google Scholar] [CrossRef]
Khennaoui, A.A.; Ouannas, A.; Bendoukha, S.; Grassi, G.; Lozi, R.P.; Pham, V.T. On fractional–order discrete–time systems: Chaos, stabilization and synchronization. Chaos Solitons Fractals 2019, 119, 150–162. [Google Scholar] [CrossRef]
Wu, G.C.; Baleanu, D. Discrete fractional logistic map and its chaos. Nonlinear Dyn. 2014, 75, 283–287. [Google Scholar] [CrossRef]
Alpar, O. Dynamics of a new generalized fractional one-dimensional map: Quasiperiodic to chaotic. Nonlinear Dyn. 2018, 94, 1377–1390. [Google Scholar] [CrossRef]
Khennaoui, A.A.; Ouannas, A.; Bendoukha, S.; Grassi, G.; Wang, X.; Pham, V.T.; Alsaadi, F.E. Chaos, control, and synchronization in some fractional-order difference equations. Adv. Differ. Equ. 2019, 1, 1–23. [Google Scholar] [CrossRef] [Green Version]
Peng, Y.; Sun, K.; Peng, D.; Ai, W. Dynamics of a higher dimensional fractional-order chaotic map. Phys. A Stat. Mech. Its Appl. 2019, 525, 96–107. [Google Scholar] [CrossRef]
Wang, L.; Sun, K.; Peng, Y.; He, S. Chaos and complexity in a fractional-order higher-dimensional multi-cavity chaotic map. Chaos Solitons Fractals 2020, 131, 109488. [Google Scholar] [CrossRef]
Rong, K.; Bao, H.; Li, H.; Hua, Z.; Bao, B. Memristive Hénon map with hidden Neimark-Sacker bifurcations. Nonlinear Dyn. 2022, 108, 4459–4470. [Google Scholar] [CrossRef]
Bao, H.; Li, H.; Hua, Z.; Xu, Q.; Bao, B.C. Sine Transform-Based Memristive Hyperchaotic Model with Hardware Implementation. IEEE Trans. Ind. Inform. 2022; in press. [Google Scholar] [CrossRef]
Li, K.; Bao, H.; Li, H.; Ma, J.; Hua, Z.; Bao, B. Memristive Rulkove neuron model with magnetic inclution effects. IEEE Trans. Ind. Inform. 2021, 18, 1726–1736. [Google Scholar] [CrossRef]
Bao, H.; Hua, Z.; Li, H.; Chen, M.; Bao, B. Discrete memristor hyperchaotic maps. IEEE Trans. Circuits Syst. Regul. Pap. 2021, 68, 4534–4544. [Google Scholar] [CrossRef]
Hadjabi, F.; Ouannas, A.; Shawagfeh, N.; Khennaoui, A.A.; Grassi, G. On Two-Dimensional Fractional Chaotic Maps with Symmetries. Symmetry 2020, 12, 756. [Google Scholar] [CrossRef]
Liu, T.; Mou, J.; Banerjee, S.; Cao, Y.; Han, X. A new fractional-order discrete BVP oscillator model with coexisting chaos and hyperchaos. Nonlinear Dyn. 2021, 106, 1011–1026. [Google Scholar] [CrossRef]
Amatroud, A.O.; Khennaoui, A.A.; Ouannas, A.; Pham, V.T. Infinite Line of Equilibrium in a Novel Fractional Map with Coexisting Attractors and Initial Offset Boosting. Available online: https://www.degruyter.com/document/doi/10.1515/ijnsns-2020-0180/pdf (accessed on 5 December 2022).
Pincus, S.M. Approximate entropy as a measure of system complexity. Proc. Natl. Acad. Sci. USA 1991, 88, 2297–2301. [Google Scholar] [CrossRef] [Green Version]
Ouannas, A.; Khennaoui, A.A.; Momani, S.; Grassi, G.; Pham, V.T. Chaos and control of a three-dimensional fractional order discrete-time system with no equilibrium and its synchronization. AIP Adv. 2020, 10, 045310. [Google Scholar] [CrossRef]
Ouannas, A.; Bendoukha, S. Generalized and inverse generalized synchronization of fractional–order discrete–time chaotic systems with non–identical dimensions. Adv. Differ. Equ. 2018, 2018, 303. [Google Scholar] [CrossRef]
Bendoukha, S.; Ouannas, A.; Wang, X.; Khennaoui, A.A.; Pham, V.T.; Grassi, G.; Huynh, V.V. The co-existence of different synchronization types in fractional-order discrete-time chaotic systems with non–identical dimensions and orders. Entropy 2018, 20, 710. [Google Scholar] [CrossRef] [Green Version]
Ouannas, A.; Khennaoui, A.A.; Zehrour, O.; Bendoukha, S.; Grassi, G.; Pham, V.T. Synchronisation of integer-order and fractional-order discrete-time chaotic systems. Pramana 2019, 92, 52. [Google Scholar] [CrossRef]
Bao, B.C.; Li, H.Z.; Zhu, L.; Zhang, X.; Chen, M. Initial-switched boosting bifurcations in 2D hyperchaotic map. Chaos Interdiscip. J. Nonlinear Sci. 2020, 30, 033107. [Google Scholar] [CrossRef]
Khennaoui, A.A.; Almatroud, A.O.; Ouannas, A.; Al-sawalha, M.M.; Grassi, G.; Pham, V.T.; Batiha, I.M. An Unprecedented 2-Dimensional Discrete-Time Fractional-Order System and Its Hidden Chaotic Attractors. Math. Probl. Eng. 2021, 2021, 6768215. [Google Scholar] [CrossRef]
Wu, G.C.; Baleanu, D. Jacobian matrix algorithm for Lyapunov exponents of the discrete fractional maps. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 95–100. [Google Scholar] [CrossRef]
Zhang, S.; Zeng, J.; Wang, X.; Zeng, Z. A novel no-equilibruim HR neuron model with hidden homogeneous extreme multistability. Chaos Solitons Frractals 2021, 145, 11761. [Google Scholar]
Li, C.; Sprott, J.C. Finding coexisting attractors using amplitude control. Nonlinear Dyn. 2014, 780, 2059–2064. [Google Scholar] [CrossRef]
Leutcho, G.D.; Wang, H.; Kengne, R.; Kengne, L.K.; Njitacke, Z.T.; Fozin, T.F. Symmetry breaking amplitude control and constant Lyapunov exponent based on single parameter snap flows. Eur. Phys. J. Spec. Top. 2021, 230, 1887–1903. [Google Scholar] [CrossRef]
Viayajumar, M.D.; Natiq, H.; Leutcho, G.D.; Rajagopal, K.; Jafri, S.; Hussain, I. Hidden and self-Excited Collective Dynamics of a New Multistable Hyper-Jerk System with Unique Equilibrium. Int. J. Bifurc. Chaos 2022, 32, 2250063. [Google Scholar] [CrossRef]
Cafagna, D.; Grassi, G. An effective method for detecting chaos in fractional-order systems. Int. J. Bifurc. Chaos 2010, 20, 669–678. [Google Scholar] [CrossRef]
Baleanu, D.; Wu, G.C.; Bai, Y.R.; Chen, F.L. Stability analysis of Caputoâ like discrete fractional systems. Commun. Nonlinear Sci. Numer. Simul. 2017, 48, 520–530. [Google Scholar] [CrossRef]

Figure 1. Bifurcation diagrams and LEs of the FODTS (3) with no equilibria for fractional order value

η = 0.98

: (a) for

A \in [2.3, 2.9]

;

B = 0.1

, and initial value

{x_{1}}_{0} = 1; {x_{2}}_{0} = - 2

; (b) for

A \in [- 2.9, - 2.3]

;

B = - 0.1

, and initial value

{x_{1}}_{0} = - 1; {x_{2}}_{0} = 2

Figure 1. Bifurcation diagrams and LEs of the FODTS (3) with no equilibria for fractional order value

η = 0.98

: (a) for

A \in [2.3, 2.9]

;

B = 0.1

, and initial value

{x_{1}}_{0} = 1; {x_{2}}_{0} = - 2

; (b) for

A \in [- 2.9, - 2.3]

;

B = - 0.1

, and initial value

{x_{1}}_{0} = - 1; {x_{2}}_{0} = 2

Figure 2. Hidden attractors of the FODTS (3) for the symmetrical initial conditions

({x_{1}}_{0}, {x_{2}}_{0}) = (\pm 1, \pm 2)

, fractional order

η = 0.98

, positive system parameters in red diagrams, and negative system parameters in blue diagrams: (a) hidden periodic attractors for

A = \pm 2.62

and

B = \pm 0.1

; (b) hidden chaotic attractors for

A = \pm 2.772

and

B = \pm 0.1

; (c) hidden hyperchaotic attractors for

A = \pm 2.85

and

B = \pm 0.1

Figure 2. Hidden attractors of the FODTS (3) for the symmetrical initial conditions

({x_{1}}_{0}, {x_{2}}_{0}) = (\pm 1, \pm 2)

, fractional order

η = 0.98

, positive system parameters in red diagrams, and negative system parameters in blue diagrams: (a) hidden periodic attractors for

A = \pm 2.62

and

B = \pm 0.1

; (b) hidden chaotic attractors for

A = \pm 2.772

and

B = \pm 0.1

; (c) hidden hyperchaotic attractors for

A = \pm 2.85

and

B = \pm 0.1

Figure 3. Bifurcation diagrams and Lyapunov exponents of the FODTS (3) with the lines of the equilibrium points in the

A - x_{2}

plane for the fractional order value

η = 0.98

: (a) for

A \in [2.3, 2.9]

;

B = 0

, and IV

x_{10} = 1; x_{20} = - 2

; (b) for

A \in [- 2.9, - 2.3]

;

B = 0

, and initial conditions

x_{10} = - 1; x_{20} = 2

Figure 3. Bifurcation diagrams and Lyapunov exponents of the FODTS (3) with the lines of the equilibrium points in the

A - x_{2}

plane for the fractional order value

η = 0.98

: (a) for

A \in [2.3, 2.9]

;

B = 0

, and IV

x_{10} = 1; x_{20} = - 2

; (b) for

A \in [- 2.9, - 2.3]

;

B = 0

, and initial conditions

x_{10} = - 1; x_{20} = 2

Figure 4. Hidden attractors of the FODTS (3) with the lines of the equilibrium points for the symmetrical IV

(x_{10}, x_{20}) = (\pm 1, \pm 2)

, fractional order

η = 0.98

, and positive system parameters in red, and negative system parameters in blue with

B = 0

: (a) hidden periodic attractors for

A = \pm 2.752

, (b) hidden hyperchaotic attractors for

A = \pm 2.797

, (c) hidden chaotic attractors for

A = \pm 2.867

Figure 4. Hidden attractors of the FODTS (3) with the lines of the equilibrium points for the symmetrical IV

(x_{10}, x_{20}) = (\pm 1, \pm 2)

, fractional order

η = 0.98

, and positive system parameters in red, and negative system parameters in blue with

B = 0

: (a) hidden periodic attractors for

A = \pm 2.752

, (b) hidden hyperchaotic attractors for

A = \pm 2.797

, (c) hidden chaotic attractors for

A = \pm 2.867

Figure 5. Hidden attractors, bifurcation diagrams, and Lyapunov exponents (LEs) versus

η

of the FODTS (3) for the symmetrical initial conditions

({x_{1}}_{0}, {x_{2}}_{0}) = (\pm 1, \pm 2)

, and positive control parameter

A = 2.7

B = 0.1

(red diagram), and negative control parameters

A = - 2.7

B = - 0.1

(blue diagram).

Figure 5. Hidden attractors, bifurcation diagrams, and Lyapunov exponents (LEs) versus

η

of the FODTS (3) for the symmetrical initial conditions

({x_{1}}_{0}, {x_{2}}_{0}) = (\pm 1, \pm 2)

, and positive control parameter

A = 2.7

B = 0.1

(red diagram), and negative control parameters

A = - 2.7

B = - 0.1

(blue diagram).

Figure 6. Different bifurcation diagrams and Lyapunov exponents of the FODTS (3) with the variation of control parameter A and fractional order

η

Figure 6. Different bifurcation diagrams and Lyapunov exponents of the FODTS (3) with the variation of control parameter A and fractional order

η

Figure 7. Bifurcation and LEs of the FODTS (3) versus IV

{x_{2}}_{0}

, for parameter values

A = 2.7

B = 0

and fractional order

η = 0.98

Figure 7. Bifurcation and LEs of the FODTS (3) versus IV

{x_{2}}_{0}

, for parameter values

A = 2.7

B = 0

and fractional order

η = 0.98

Figure 8. Five coexisting hidden hyperchaotic attractors for different initial conditions and control parameters

A = 2.7, B = 0.1

and order

η = 0.98

x_{20} = 1

Figure 8. Five coexisting hidden hyperchaotic attractors for different initial conditions and control parameters

A = 2.7, B = 0.1

and order

η = 0.98

x_{20} = 1

Figure 9. Offset boosting of the fractional map (3) for varying initial condition

y (0)

with

x (0) = 1

and fractional order

μ = 0.95

. Bifurcation diagrams as the control parameter a increases in

[2.3, 2.9]

and

b = 0.1

Figure 9. Offset boosting of the fractional map (3) for varying initial condition

y (0)

with

x (0) = 1

and fractional order

μ = 0.95

. Bifurcation diagrams as the control parameter a increases in

[2.3, 2.9]

and

b = 0.1

Figure 10. Four hidden chaotic attractors with

{x_{2}}_{0} = - 2 - 2 π

(grey),

{x_{2}}_{0} = - 2

(blue),

{x_{2}}_{0} = - 2 + 2 π

(red),

{x_{2}}_{0} = - 2 + 4 π

(magenta).

Figure 10. Four hidden chaotic attractors with

{x_{2}}_{0} = - 2 - 2 π

(grey),

{x_{2}}_{0} = - 2

(blue),

{x_{2}}_{0} = - 2 + 2 π

(red),

{x_{2}}_{0} = - 2 + 4 π

(magenta).

Figure 11. Rescaled phase diagrams under different total amplitude controllers

c_{1}

of FODTS with

{x_{1}}_{0} = 1, {x_{2}}_{0} = - 2

, and fractional order

η = 0.98

c_{1} = - 5

(green diagram),

c_{1} = 2

(red diagram),

c = 12

(magenta diagram),

c_{1} = 13

, (blue diagram).

Figure 11. Rescaled phase diagrams under different total amplitude controllers

c_{1}

of FODTS with

{x_{1}}_{0} = 1, {x_{2}}_{0} = - 2

, and fractional order

η = 0.98

c_{1} = - 5

(green diagram),

c_{1} = 2

(red diagram),

c = 12

(magenta diagram),

c_{1} = 13

, (blue diagram).

Figure 12. The approximate entropy ApEn of the FODTS (3) in the three-dimensional space with the variation of system parameter A and fractional order

η

, for

B = 0.1

Figure 12. The approximate entropy ApEn of the FODTS (3) in the three-dimensional space with the variation of system parameter A and fractional order

η

, for

B = 0.1

Figure 13. Evolution of states and the phase space plot of the controlled FODTS (3) with fractional order

η = 0.98

and system parameters

A = 2.8, B = 0.1

Figure 13. Evolution of states and the phase space plot of the controlled FODTS (3) with fractional order

η = 0.98

and system parameters

A = 2.8, B = 0.1

Figure 14. Evolution of synchronization error states with FODTS

η = 0.98

and system parameters

A = 2.8, B = 0.1

Figure 14. Evolution of synchronization error states with FODTS

η = 0.98

and system parameters

A = 2.8, B = 0.1

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

Khennaoui, A.-A.; Ouannas, A.; Bekiros, S.; Aly, A.A.; Alotaibi, A.; Jahanshahi, H.; Alsubaie, H. Hidden Homogeneous Extreme Multistability of a Fractional-Order Hyperchaotic Discrete-Time System: Chaos, Initial Offset Boosting, Amplitude Control, Control, and Synchronization. Symmetry 2023, 15, 139. https://doi.org/10.3390/sym15010139

AMA Style

Khennaoui A-A, Ouannas A, Bekiros S, Aly AA, Alotaibi A, Jahanshahi H, Alsubaie H. Hidden Homogeneous Extreme Multistability of a Fractional-Order Hyperchaotic Discrete-Time System: Chaos, Initial Offset Boosting, Amplitude Control, Control, and Synchronization. Symmetry. 2023; 15(1):139. https://doi.org/10.3390/sym15010139

Chicago/Turabian Style

Khennaoui, Amina-Aicha, Adel Ouannas, Stelios Bekiros, Ayman A. Aly, Ahmed Alotaibi, Hadi Jahanshahi, and Hajid Alsubaie. 2023. "Hidden Homogeneous Extreme Multistability of a Fractional-Order Hyperchaotic Discrete-Time System: Chaos, Initial Offset Boosting, Amplitude Control, Control, and Synchronization" Symmetry 15, no. 1: 139. https://doi.org/10.3390/sym15010139

APA Style

Khennaoui, A.-A., Ouannas, A., Bekiros, S., Aly, A. A., Alotaibi, A., Jahanshahi, H., & Alsubaie, H. (2023). Hidden Homogeneous Extreme Multistability of a Fractional-Order Hyperchaotic Discrete-Time System: Chaos, Initial Offset Boosting, Amplitude Control, Control, and Synchronization. Symmetry, 15(1), 139. https://doi.org/10.3390/sym15010139

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

Article Menu

Hidden Homogeneous Extreme Multistability of a Fractional-Order Hyperchaotic Discrete-Time System: Chaos, Initial Offset Boosting, Amplitude Control, Control, and Synchronization

Abstract

1. Introduction

2. Description and Analysis of the Fractional Hyperchaotic Map

3. Dynamical Analysis and Numerical Simulations

3.1. Hidden Attractors and Bifurcation Analysis

3.2. The Effect of Fractional Order $η$

3.3. Hidden Extreme Homogeneous Multistability

3.4. Initial Offset Boosting

4. Amplitude Control Analysis

5. Complexity Analysis of the FODTS

6. Control Laws

7. Synchronization

8. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

Article Menu

Hidden Homogeneous Extreme Multistability of a Fractional-Order Hyperchaotic Discrete-Time System: Chaos, Initial Offset Boosting, Amplitude Control, Control, and Synchronization

Abstract

1. Introduction

2. Description and Analysis of the Fractional Hyperchaotic Map

3. Dynamical Analysis and Numerical Simulations

3.1. Hidden Attractors and Bifurcation Analysis

3.2. The Effect of Fractional Order η

3.3. Hidden Extreme Homogeneous Multistability

3.4. Initial Offset Boosting

4. Amplitude Control Analysis

5. Complexity Analysis of the FODTS

6. Control Laws

7. Synchronization

8. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

3.2. The Effect of Fractional Order $η$