Abstract

In this paper, a new fractional-order discrete noninvertible map of cubic type is presented. Firstly, the stability of the equilibrium points for the map is examined. Secondly, the dynamics of the map with two different initial conditions is studied by numerical simulation when a parameter or a derivative order is varied. A series of attractors are displayed in various forms of periodic and chaotic ones. Furthermore, bifurcations with the simultaneous variation of both a parameter and the order are also analyzed in the three-dimensional space. Interior crises are found in the map as a parameter or an order varies. Thirdly, based on the stability theory of fractional-order discrete maps, a stabilization controller is proposed to control the chaos of the map and the asymptotic convergence of the state variables is determined. Finally, the synchronization between the proposed map and a fractional-order discrete Loren map is investigated. Numerical simulations are used to verify the effectiveness of the designed synchronization controllers.

1. Introduction

In the recent several decades, chaos is an attractive phenomenon in nonlinear dynamical systems, which has been extensively analyzed and studied deeply. It is well known that chaos was first detected in continuous nonlinear systems. Its characteristics and the existence in discrete dynamical maps have also been the interesting topics. Many discrete maps with chaotic attractors have been proposed, such as the Logistic map, Hénon map, and Lozi map [15]. With the rapid development of fractional calculus, many works including chaos, control, and synchronization for fractional-order continuous systems have been reported [618]. It should be pointed out that the fractional discrete maps were not paid enough attention and properly explored until recently [19].

In 1974, Diaz and Olser first put forward the fractional difference [20]. Up to now, fractional-order discrete maps have obtained more and more attention. In [21], a discrete fractional Hénon map was introduced, and its chaotic behavior was discussed. Dynamics, stabilization, and synchronization for several fractional-order maps, such as the Ikeda map, Loren map, and Lozi map, were studied in [2228]. The discrete fractional calculus can avoid the tedious information or calculation error of the numerical discretization for the continuous ones due to the nonlocal property [29]. Therefore, more and more discrete maps with fractional operators need to be presented, and more abundant and complex dynamics behaviors need to be explored. Besides, it is well known that fractional-order discrete maps are sensitive not only to the small disturbance of parameters and initial conditions but also to the variation of fractional orders [30], which are the unique advantages of fractional-order systems. For this reason, a fractional-order discrete map is more suitable for data encryption and secure communications. Furthermore, fractional-order discrete maps have simple forms and rich dynamics, which are good for system analysis and numerical computation. Based on these, investigation of a new fractional-order discrete map including dynamics, stabilization, and synchronization is necessary and important for the development of fractional calculus.

In [31, 32], a two-dimensional noninvertible map with cubic order nonlinearity, which was taken as a chaotic cryptosystem, was proposed and studied. The evolution of attractors and their basins have been analyzed deeply and explained thoroughly. A noncyclic chaotic attractor for the map was displayed in [33]. Based on these, we extend the map to the fractional case and study its dynamics. The stability of the equilibrium points for the map is examined. By the bifurcation graphs and phase diagrams, the dynamics of the fractional-order discrete map with two different initial conditions is displayed as a parameter or a derivative order varies. Furthermore, bifurcations with the simultaneous variation of both a parameter and the order are also analyzed in the three-dimensional space. Interior crises occur in the map with the variation of a parameter or the order.

The main motivation of our work is to know whether bifurcations and chaos, which the integer-order discrete map possesses, also exist in the fractional-order counterpart. In fact, these dynamics behaviors do exist in the fractional-order map, and multifaceted complex dynamics is observed by means of the numerical simulations. For a chaotic system, control and synchronization are very important for its application in practical problems. In our work, we are also interested in studying the control and synchronization for the fractional-order map. Based on the stability theory of fractional-order discrete maps, a stabilization controller is proposed to control the chaos of the map. The synchronization between the proposed map and a fractional-order discrete Loren map is studied and realized.

2. Discrete Fractional Calculus

In this section, we will recall the definition and related theories for the discrete fractional calculus. In the following, the symbol represents the order Caputo type delta fractional calculus of a function with [34], which is expressed as follows:where is the order, , and . In formula (1), the fractional sum of is defined as [35, 36]where . represents the falling function which is defined according to the Gamma function as

Generally speaking, the following method is employed to compute the numerical solutions for a fractional-order discrete map.

For an equation with the fractional calculus operator,the equivalent discrete integral one iswhere .

The following theorem is used to analyze the stabilization and synchronization for fractional discrete maps. For the proof of the theorem, please refer to the literature [37].

Theorem 1. The zero equilibrium of a linear fractional discrete system:where , is asymptotically stable iffor all the eigenvalues of .

3. A Fractional-Order Discrete Map

3.1. Description of the Map

Firstly, the two-dimensional discrete map with cubic nonlinearity in [3133] is described as follows:where are the state variables and parameters. The first-order difference of (8) is formulated as

By employing the Caputo-like delta difference given in (1) with the starting point , the corresponding fractional map is

Based on equations (4) and (5), we can obtainwhere is the discrete kernel function and . Based on this, the numerical solution for the fractional discrete map (10) can be obtained, which is as follows:

In the rest of the paper, the low limit is fixed as 0.

3.2. Stability of Equilibrium Points

Now, we turn to study the stability of equilibrium points for map (10). By the simple computation, we can get the three equilibrium points:when . The map has only one equilibrium point when . The Jacobian matrix of map (10) evaluated at an equilibrium point is

The corresponding eigenvalues for the equilibrium point are . In this paper, we only consider the case of map (10) with positive parameters. Therefore, the zero equilibrium point is unstable due to on the basis of Theorem 1.

For a fractional-order discrete map, the stability of a zero equilibrium point can be determined easily based on Theorem 1. Therefore, we will use a very simple method proposed in [38] for handling with the nonzero equilibrium points. For more special details about the method, please refer to Remark 2.5 in the literature [38].

In order to analyze the stability of the nonzero equilibrium points , let . Through the following variable transforms, we can get two newly maps with a zero equilibrium point:which are corresponding to , respectively. The Jacobian matrix of maps (16) and (17) evaluated at the zero equilibrium point is

The corresponding eigenvalues of for the zero equilibrium point are

When the parameters of map (10) are choosen as and the order is taken as , the corresponding eigenvalues are . By simple computation, we can obtainwhich means the equilibrium point is unstable in this case according to Theorem 1.

4. Dynamics Analysis

In this section, the dynamics of the fractional-order discrete map (10) with the variation of a parameter or the fractional order and the bifurcations with the simultaneous variation of both a parameter and the order will be analyzed in detail.

4.1. Dynamics as the Parameter Varies

When the order and parameter , the dynamics of map (10) with the variation of the parameter is analyzed. The bifurcation diagrams and the corresponding largest Lyapunov exponents (LLE) spectrums with two different initial conditions and are displayed in Figure 1. From which we can see that the dynamics of the map is abundant and shows a symmetry with different initial conditions in this case.

The evolution of the trajectories for different with is depicted in Figure 2. In the phase plane, there is a fixed point which means the map is period-1 for (Figures 2(a) and 2(b)). The map has a limit cycle attractor for , see Figure 2(c), which means a Hopf bifurcation occurs as the parameter increases from 1 to 1.5. The shape of the limit cycle changes as increases further (Figures 2(d) and 2(e)). The map keeps chaotic when varies from 1.75 to 2.2. In Figure 2(g), three small chaotic attractors appear in the phase plane when and convert to one large attractor when (Figure 2(h)). From Figures 2(i) and 2(j), we can see that the chaotic attractor becomes a large one suddenly, which implies that an interior crisis occurs as increases from 2 to 2.2.

The dynamics of the map with the initial condition , which is similar to Figure 2, is displayed in Figure 3. From the global dynamics perspective, the two chaotic attractors, which are depicted in Figures 2(i) and 3(i), collide with each other and convert to a large one (Figures 2(j) and 3(j)).

4.2. Dynamics as the Parameter Varies

The fractional order is fixed as 0.98 and the parameter , and the dynamics of map (10) when is varied in the interval of [0.2, 1] with two initial conditions and are studied. Firstly, the bifurcation diagrams and the corresponding LLE spectrums with and are plotted in Figure 4. From which we can see that the dynamics of the map is different from that of the map as the parameter varies and lost the symmetry.

The phase diagrams for the map with are shown in Figure 5. The map is period-2 when (Figures 5(a) and 5(b)), and the periodic behavior persists for a long time until . When , the system is chaotic and the chaotic attractor is depicted in Figure 5(c). As increases from 0.87 to 0.88, the chaotic attractor becomes a large one, which means an interior crisis occurs.

The dynamics of the map with the initial condition is displayed in Figure 6. The map keeps chaotic as the parameter changes in the interval of [0.2, 1], and the chaotic attractor has different formations. A chaotic attractor consists of four small parts in the phase space when the parameter increases from 0.2 to 0.28 (Figures 6(a)6(c)) and converts to a whole one when (Figure 6(d)). The system is period-5 when , and the phase diagram and the discrete time evolution of the state variable are plotted in Figures 6(e) and 6(f). When the parameter increases from 0.87 to 0.88, the chaotic attractor becomes a large one suddenly (Figures 6(g) and 6(h)), which implies that an interior crisis occurs.

The stable region for the map in the parameter plane with is plotted in Figure 7 in order to give a guidance of choosing values of the parameters. From this figure, we can see that map (10) is chaotic when and .

4.3. Dynamics as the Order Varies

The parameters are fixed as and , the dynamics of map (10) is studied when the order is varied in this section. The bifurcation diagrams and the corresponding LLE spectrums with and are plotted in Figure 8. From which it can be seen that the dynamics of the map also shows a symmetry with different initial conditions in this case.

For different values of the , the phase diagrams with are shown in Figure 9. The map has period-1 attractor for (Figure 9(a)), and a limit cycle for (Figure 9(b)), which means a Hopf bifurcation, occurs as the order increases. The shape of the limit cycle changes as increases (Figure 9(c)). The map exists as a multicycle attractor for , see Figure 9(d). As the order increases to 0.84, the attractor becomes a chaotic one which is consisted by several small parts (Figure 9(e)). The small parts combine into one attractor when (Figure 9(f)). From Figures 9(g) and 9(h), it is clear that the chaotic attractor has three small parts in the phase plane when , and these parts become a whole one when . The chaotic attractor in Figure 9(i) becomes a large one (Figure 9(j)) when the order varies from 0.95 to 0.96.

The phase diagrams with initial condition as the order varies from 0.7 to 0.96 are shown in Figure 10, which are symmetric with those of the map with. From the global dynamics perspective, the two chaotic attractors, which are depicted in Figures 9(i) and 10(i), collide with each other and convert to a large one (Figures 9(j) and 10(j)).

4.4. Bifurcation with the Simultaneous Variation of Both Parameter and Order

In this section, the bifurcations of map (10) as a parameter and the order variation with the initial conditions and are studied. Firstly, the value of parameter is fixed as 0.95. The bifurcation diagram of the map is depicted in Figure 11, when the parameter and the order change simultaneously. Secondly, the value of parameter is set as 2.2. Figure 12 shows the bifurcation diagram of the map with the variation of the parameter and the order . From these figures, it can be seen that map (10) is periodic when the order is less than a certain threshold and appears chaotic behavior when the order is greater than the certain threshold.

In other words, the dynamics of map (10) becomes regular as the derivate order decreases from 1 to 0.5 and complex as the derivate order increases from 0.5 to 1.

5. Stabilization

The stabilization of map (10) will be studied in this section. From here, the case of the equilibrium to be at the origin is considered. For convenience, the controlled map (10) is rewritten as the following form:where and are the stabilization controllers.

Theorem 2. The fractional-order map (10) can be stabilized when the controllers are designed as the following form:

Proof. Through substituting (22) into (21), map (21) becomesMap (23) can be rewritten in the compact form:where . Based on Theorem 1, it is easy to see that the eigenvalues of satisfy the conditions , which implies that the chaos of map (10) can be controlled and the zero equilibrium of (24) is globally asymptotically stable.
In the numerical simulations, the values of the parameters are fixed as , and the fractional order . The controllers are used to stabilize map (10) when the iterations are chosen as . The stabilization results are displayed in Figure 13. It is clear that the state variables toward to zeros, which means the chaos of map (10) is stabilized and the results confirm the theoretical control results presented in Theorem 2.

6. Synchronization

Now, we consider the synchronization of map (10). Firstly, a fractional Lorenz map is taken as the drive system:where and state subscript denotes the drive system. Map (25) is chaotic when the parameters and the derivative order . For more details about the dynamics of map (25), please refer to the literature [23]. Map (10) with synchronization controllers and , which can be described by the following equations:where the subscript denotes the response system. The error state variables are defined as . If all the error states variables tend to 0 as the time , then maps (25) and (26) are synchronized. The following theorem is given to ensure the synchronization between the two maps can be realized.

Theorem 3. The drive and response maps (25) and (26) are synchronized when the controllers are designed as follows:

Proof. The error dynamical system with fractional Caputo difference isBy substituting controllers (27) into (28), we can obtain the error dynamical system:System (29) is rewritten in the compact form:where . It can be seen that the eigenvalues of the matrix satisfy the following stability condition:Based on Theorem 1, we can get that the zero equilibrium of (29) is globally asymptotically stable, which implies the two maps (25) and (26) are synchronized.
In the numerical simulations, the parameters are fixed as and the order is . The initial conditions of the two systems (25) and (26) are . The results of numerical simulations are depicted in Figure 14. From which, we can see that the error variables converge to zero rapidly as increases (Figures 14(a) and 11(b)). Meanwhile, the evolution of the state variables with time for the two maps (25) and (26) are synchronized under the designed controllers (27) (Figures 14(c) and 14(d)).

7. Conclusions and Discussion

A fractional-order discrete noninvertible map with cubic nonlinearity is proposed in this paper. Firstly, the stability of the equilibrium points for the map is analyzed. Secondly, the dynamics of the map with two different initial conditions is studied by numerical simulation. Bifurcation diagrams and phase plots are obtained as a parameter or the fractional order varies. A series of attractors of the map in different forms, including equilibrium points, limit cycles, and chaotic attractors, are plotted. Furthermore, bifurcations with the simultaneous variation of both parameter and order are also analyzed in the three-dimensional space. From the global dynamics perspective, interior crises occur in the map as a parameter or the order varies. Thirdly, based on the stability theory of fractional-order discrete maps, the chaos of the map is controlled by the stabilization controllers. Finally, the synchronization between the proposed map and a fractional-order discrete Loren map is investigated. Numerical simulations are implemented to verify the effectiveness of the designed controllers. The results obtained in this paper reveal that chaos really exists in the fractional-order formation for the map proposed in [3133]. More abundant local and global dynamics are found in the fractional-order map.

It is worth mentioning that the mechanism of interior crises occur in map (10) cannot be displayed from a global perspective due to the absence of effective global dynamics computation methods for fractional-order discrete maps. Therefore, developing the effective computation methods of global analysis for this kind of systems is our following work.

Data Availability

The data for the bifurcation diagrams used to support the findings of this study are included within the supplementary information files (available here). For the large of the data, we will supply the bifurcation data calculated by the software Matlab.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 11702194 and 11702195) and Natural Science Preparatory Study Foundation of Xi’an University of Posts and Telecommunications (no. 106/205020030).

Supplementary Materials

(1) discretesanweiIN1q098b and discretesanweiIN2q098b: the data of Figures 1(a) and 1(b). (2) discretesanweiIN1q098c and discretesanweiIN2q098c: the data of Figures 4(a) and 4(b). (3) discretesanweiqIN1 and discretesanweiqIN2: the data of Figures 8(a) and 8(b). (Supplementary Materials)