# Bifurcations of a New Fractional-Order System with a One-Scroll Chaotic Attractor

Authors: Xiaojun Liu; Ling Hong; Lixin Yang; Dafeng Tang
Journal: Discrete Dynamics in Nature and Society
 (2019)

---

## Abstract

In this paper, a new fractional-order system which has a chaotic attractor of the one-scroll structure is presented. Firstly, the stability of the equilibrium points of the system is investigated. And based on the stability analysis, the generation conditions of the one-scroll structure for the attractor are determined. In a commensurate-order case, bifurcations with the variation of a system parameter are investigated as derivative orders decrease from 0.99. In an incommensurate-order case, bifurcations with the variation of a derivative order are analyzed as other orders decrease from 1.


---

## Body

## 1. Introduction

Fractional calculus is a field of applied mathematics that deals with derivatives and integrals of arbitrary orders. Although the seeds of fractional derivatives were planted over 300 years ago, the development of fractional calculus is very slow at an early stage for the absence of geometrical interpretation and applications. In the recent several decades, it has been applied to almost every field of science, engineering, economics, secures communication, and so on [1–5].It is well known that fractional calculus is very suitable for the description of properties of various real materials. Meanwhile, fractional calculus provides an excellent instrument for the description of memory and hereditary properties of various materials and processes. This is the main advantage of fractional derivatives in comparison with classical integer-order models, in which such effects are in fact neglected. The advantages of fractional derivatives become apparent in modeling mechanical and electrical properties of real materials [6].Many new fractional-order systems were presented in recent years. Meanwhile, rich and complex dynamics, such as periodic solutions, chaos windows, all kinds of bifurcations, boundary, and interior crises were observed in these systems. For example, Chua system with a derivative order 2.7 makes chaos motion [7]. Chaotic dynamics of a damped van der Pol equation with fractional order is investigated in [8]. Chen studied the nonlinear dynamics and chaos in a fractional-order financial system [9]. A periodically forced complex Duffing’s oscillator was proposed, and chaos for the system was studied in detail [10]. Bifurcation, chaos control, and synchronization were investigated for a fractional-order Lorenz system with the complex variables [11]. In [12], boundary and interior crises were determined in a fractional-order Duffing system by a global numerical computation method. A proposed standard for the publication of new chaotic systems was studied in [13]. The authors investigated a simple chaotic flow with a plane of equilibria in [14].It is well known that bifurcation theory concerns the changes in qualitative or topological structures of limiting motions such as equilibriums, periodic solutions, homoclinic orbits, heteroclinic orbits, and invariant tori for nonlinear evolution equations as some relevant parameters in the equations vary. Generally, the subject can be traced back to the very earlier work of Poincare around 1892 [15]. Nowadays, it is a fundamental tool to analyze nonlinear problems which enables us to understand how and when a system organizes new states and patterns near the original “trivial” one when a control parameter crosses a critical value. For fractional-order systems, a bifurcation implies a qualitative or topological change in dynamics with a variation of a system parameter or derivative order, and bifurcation analysis becomes harder due to the nonlocal property of the operator of fractional calculus. Many references have studied the bifurcations of fractional-order systems [16–20]. However, these investigations mainly focus on the bifurcation of a fractional-order system as a system parameter or a derivative order varies. To our knowledge, few works concerns bifurcations with the variation of both a system parameter and a derivative order or with the variation of both a derivative order and other orders.Compared with integer-order chaotic systems, fractional-order chaotic systems with more complex dynamic characteristics and more system parameters can provide higher security for secure communication [21, 22]. In [23], the authors investigated the synchronization of a three-dimensional integer-order system. The differential equations of the system with simple structure were similar to those of the Lorenz system. It is well known that the Lorenz system has a chaotic attractor with double-scroll structure. The system in [23] has a chaotic attractor with only one-scroll and very abundant dynamic behaviors. Motivated by the above, in this paper, a corresponding fractional-order system is proposed and studied. Firstly, the stability of equilibrium points of the system is investigated. In a commensurate-order case, bifurcations with the variation of a system parameter are investigated as derivative orders decrease from 0.99. In an incommensurate-order case, bifurcations with the variation of a derivative order are analyzed as the other orders decrease from 1. Period-doubling and saddle-node bifurcations can be observed from the bifurcation diagrams by numerical simulations.The remainder of the paper is organized as follows. In Section2, the definitions of the fractional calculus and related preliminaries are given. A new fractional-order system with one-scroll attractor is presented in Section 3. In Section 4, bifurcations in the two cases of commensurate-order and incommensurate-order are analyzed, respectively. Conclusions of the paper are drawn in Section 5.


## 2. Fractional Derivatives and Preliminaries

### 2.1. Definitions

Fractional calculus can be considered as a generalization of integration and differentiation. The operatorDtqa of fractional calculus can be defined by(1)Dtqa=dqdtqRq>01Rq=0∫atdτqRq<0,where q denotes the derivative order and R(q) corresponds to the real part of q. The numbers a and t represent the limits of the operator.In general, three definitions of fractional derivative are used frequently, namely, the Grunwald-Letnikov definition, the Riemann-Liouville, and the Caputo definitions [6, 24].The Grunwald-Letnikov definition (GL) derivative with fractional orderq can be described by(2)DtqaGL=limh→0⁡1hq∑j=0t-a/h-1qqjft-jh,where the symbol · represents the integer part.The Riemann-Liouville (RL) definition is(3)DtqaRLft=dndtn1Γn-q∫atfτt-τq-n+1dτ,n-1<q<n,where Γ(·) denotes the gamma function.The Caputo (C) fractional derivative is defined as follows:(4)DtqaCft=1Γn-q∫att-τn-q-1fnτdτ,n-1<q<n.For a fractional differential equation which is defined by Caputo derivatives, the initial condition takes on the same form as those for the integer-order ones, which can be measured easily in applications. For this reason, the Caputo derivative will be adopted in the rest of the paper.


### 2.2. Numerical Methods

Due to the nonlocal property of the operator of fractional calculus, it is not easy to obtain the numerical solutions for a fractional differential equation. Generally speaking, two approximation methods are frequently used, namely, an improved version of Adams-Bashforth-Moulton algorithm based on the predictor-correctors scheme and the frequency domain approximation [25–28]. For the accuracy [29], we will employ the improved predictor-corrector algorithm to solve a fractional differential equation in this paper.In order to get the approximate solution of a fractional-order chaotic system by the improved predictor-corrector algorithm, the following equation is considered:(5)dqxdtq=ft,x,0≤t≤Txk0=x0kk=0,1,⋯,q-1,which is equivalent to the Volterra integral equation (6)xt=∑k=0q-1x0ktkk!+1Γq∫0tt-τq-1fτ,xtdτ.Now, set h=T/N, tj=jh,(j=0,1,⋯,N). The corrector formula for (6) can thus be discretized as follows:(7)xhtn+1=∑k=0q-1x0ktn+1kk!+hqΓq+2ftn+1,xhptn+1+hqΓq+2∑j=0nαj,n+1ftj,xhtj,where predicted values xh(tn+1) are determined by the following formula:(8)xhptn+1=∑k=0q-1x0ktn+1kk!+1Γq∑j=0nβj,n+1ftj,xhtj,and(9)αj,n+1=nq+1-n-qn+1q,j=0n-j+2q+1+n-jq+1-2n-j+1q+1,1≤j≤n1,j=n+1.βj,n+1=hqqn-j+1q-n-jq,1≤j≤nThe error estimate of this approach is maxj=0,1,⋯,N⁡x(tj)-xh(tj)=O(hP), where p=min⁡(2,1+q).


### 2.3. The Stability of a Fractional-Order System

For fractional-order systems, the stability analysis of equilibrium points is complex and difficult due to the nonlocal property of fractional calculus. Here, the definitions of commensurate-order and incommensurate-order fractional-order systems will be given firstly.Definition 1.
For a fractional-order system, which can be described bydqx/dtq=f(x), where x=(x1,x2,⋯,xn)T is the state vector, q=(q1,q2,⋯,qn)T is the fractional derivative orders vector, and qi>0. The fractional-order system is commensurate-order when all the derivative orders satisfy q1=q2=⋯=qn; otherwise it is an incommensurate-order system [30].In order to investigate the stability of equilibrium points for fractional-order systems, the following lemma is used frequently.Lemma 2.
For a commensurate fractional-order system, the equilibrium points of the system are asymptotically stable if all the eigenvalues at the equilibriumE∗ satisfy the following condition:(10)arg⁡eigJ=arg⁡λj>π2q,j=1,2,⋯,n,where J is the Jacobian matrix of the system evaluated at the equilibria E∗ [31].


## 2.1. Definitions

Fractional calculus can be considered as a generalization of integration and differentiation. The operatorDtqa of fractional calculus can be defined by(1)Dtqa=dqdtqRq>01Rq=0∫atdτqRq<0,where q denotes the derivative order and R(q) corresponds to the real part of q. The numbers a and t represent the limits of the operator.In general, three definitions of fractional derivative are used frequently, namely, the Grunwald-Letnikov definition, the Riemann-Liouville, and the Caputo definitions [6, 24].The Grunwald-Letnikov definition (GL) derivative with fractional orderq can be described by(2)DtqaGL=limh→0⁡1hq∑j=0t-a/h-1qqjft-jh,where the symbol · represents the integer part.The Riemann-Liouville (RL) definition is(3)DtqaRLft=dndtn1Γn-q∫atfτt-τq-n+1dτ,n-1<q<n,where Γ(·) denotes the gamma function.The Caputo (C) fractional derivative is defined as follows:(4)DtqaCft=1Γn-q∫att-τn-q-1fnτdτ,n-1<q<n.For a fractional differential equation which is defined by Caputo derivatives, the initial condition takes on the same form as those for the integer-order ones, which can be measured easily in applications. For this reason, the Caputo derivative will be adopted in the rest of the paper.


## 2.2. Numerical Methods

Due to the nonlocal property of the operator of fractional calculus, it is not easy to obtain the numerical solutions for a fractional differential equation. Generally speaking, two approximation methods are frequently used, namely, an improved version of Adams-Bashforth-Moulton algorithm based on the predictor-correctors scheme and the frequency domain approximation [25–28]. For the accuracy [29], we will employ the improved predictor-corrector algorithm to solve a fractional differential equation in this paper.In order to get the approximate solution of a fractional-order chaotic system by the improved predictor-corrector algorithm, the following equation is considered:(5)dqxdtq=ft,x,0≤t≤Txk0=x0kk=0,1,⋯,q-1,which is equivalent to the Volterra integral equation (6)xt=∑k=0q-1x0ktkk!+1Γq∫0tt-τq-1fτ,xtdτ.Now, set h=T/N, tj=jh,(j=0,1,⋯,N). The corrector formula for (6) can thus be discretized as follows:(7)xhtn+1=∑k=0q-1x0ktn+1kk!+hqΓq+2ftn+1,xhptn+1+hqΓq+2∑j=0nαj,n+1ftj,xhtj,where predicted values xh(tn+1) are determined by the following formula:(8)xhptn+1=∑k=0q-1x0ktn+1kk!+1Γq∑j=0nβj,n+1ftj,xhtj,and(9)αj,n+1=nq+1-n-qn+1q,j=0n-j+2q+1+n-jq+1-2n-j+1q+1,1≤j≤n1,j=n+1.βj,n+1=hqqn-j+1q-n-jq,1≤j≤nThe error estimate of this approach is maxj=0,1,⋯,N⁡x(tj)-xh(tj)=O(hP), where p=min⁡(2,1+q).


## 2.3. The Stability of a Fractional-Order System

For fractional-order systems, the stability analysis of equilibrium points is complex and difficult due to the nonlocal property of fractional calculus. Here, the definitions of commensurate-order and incommensurate-order fractional-order systems will be given firstly.Definition 1.
For a fractional-order system, which can be described bydqx/dtq=f(x), where x=(x1,x2,⋯,xn)T is the state vector, q=(q1,q2,⋯,qn)T is the fractional derivative orders vector, and qi>0. The fractional-order system is commensurate-order when all the derivative orders satisfy q1=q2=⋯=qn; otherwise it is an incommensurate-order system [30].In order to investigate the stability of equilibrium points for fractional-order systems, the following lemma is used frequently.Lemma 2.
For a commensurate fractional-order system, the equilibrium points of the system are asymptotically stable if all the eigenvalues at the equilibriumE∗ satisfy the following condition:(10)arg⁡eigJ=arg⁡λj>π2q,j=1,2,⋯,n,where J is the Jacobian matrix of the system evaluated at the equilibria E∗ [31].


## 3. A New Fractional-Order System

In this section, a new fractional-order system which consist of three differential equations is proposed and can be denoted as follows:(11)Dq1x=ax-yDq2y=x-z-byDq3z=ax+4zy-2where x, y, and z are state variables of the system, a and b the system parameters, and q1, q2, and q3 derivative orders.When the derivative orders are selected asq1=q2=q3=q=0.99, the system parameters a=0.3 and b=0.02, and initial conditions (x0,y0,z0)=(-0.5,-1,1), system (11) is chaotic. In Figures 1(a) and 1(b), a chaotic attractor with a one-scroll on three-dimensional space and projected onto x-y plane are depicted. In this case, the corresponding Lyapunov exponents are λ1=0.117, λ2=0, and λ3=-10.49. When the derivative order q=0.985 and q=0.965, the corresponding attractors are displayed in Figures 1(c)–1(f).Figure 1
The chaotic attractor of system (11). (a) 3D plot of the attractor when q=0.99; (b) the attractor projected onto x-y plane when q=0.99; (c) 3D plot of the attractor when q=0.985; (d) the attractor projected onto x-y plane when q=0.985; (e) 3D plot of the attractor when q=0.965; (f) the attractor projected onto x-y plane when q=0.965.
(a) (b) (c) (d) (e) (f)The equilibrium points of system (11) can be calculated by solving the equations Dq1x=0, Dq2y=0, and Dq3z=0. The system contains two equilibriums, i.e.,(12)E10,0,0E22a-141-ab,a2a-141-ab,1-ab2a-141-ab.The equilibrium E2 exists when the system parameters satisfy the condition 1-ab≠0 or 1-ab≠a.The Jacobin matrix for system (11) evaluated at the equilibrium point (x∗,y∗,z∗) is given by(13)Jx∗,y∗,z∗=a-101-b-1a4z∗4y∗-2.Based on the matrix, the characteristic equation at the equilibriums E1 is (14)λ3+8+b-aλ2+8b-a+1-abλ+81-ab-a=0.By the Routh-Hurwitz test, the equilibrium point E1 has three roots with the negative real parts if and only if 8+b-a>0 and 0<81-ab-a<(8+b-a)[8(b-a)+1-ab]. In our case, when the values of the system parameters are taken as a=0.3, b=0.02, the corresponding eigenvalues for the equilibrium point E1 are λ1=-8.0047 and λ2,3=-0.1576±1.0089i. According to Lemma 2, the equilibrium point E1 is locally stable.The characteristic equation at equilibriumE2 is(15)λ3+b-a-4y-2λ2-4y-2b-a+ab-4z+1λ+4y-2ab-1-a1+4z=0.By the Routh-Hurwitz test, the equilibrium point E2 has three roots with the negative real parts if and only if (b-a)-4(y-2)>0 and(16)0<4y-2ab-1-a1+4z<b-a-4y-24y-2b-a+ab-4z+1.The corresponding eigenvalues for the equilibrium point E2 are λ1=-7.4038 and λ2,3=0.1029±0.9653i when the system parameters take the same values as before. Then based on Lemma 2, we can get that the fixed point E2 is unstable.In a nonlinear dynamical system, a saddle point is an equilibrium point on which the equivalent linearized model has at least one eigenvalue in the stable region and one eigenvalue in the unstable region. In the same system, a saddle point is called saddle point of index 1 if one of the eigenvalues is unstable and other eigenvalues are stable. Also, a saddle point of index 2 is a saddle point with one stable eigenvalue and two unstable eigenvalues. In chaotic systems, it is proved that scrolls are generated only around the saddle points of index 2. Moreover, saddle points of index 1 are responsible only for connecting scrolls [32, 33].From the above analysis we can see that the equilibrium pointE2 is a saddle point of index 2. In chaotic systems, it is proved that the scrolls of a chaotic attractor are generated only around the saddle points of index 2. Moreover, the saddle points of index 1 are responsible for connecting the scrolls. In the fractional-order system (11), the equilibrium E1 is not a saddle point of index 1. The necessary condition for the existence of double-scroll attractor in system (11) cannot be satisfied. Therefore, the chaotic attractor of system (11) is one scroll.Compared with an integer-order system, the derivative order is an important parameter for a fractional-order system. For system (11), the system parameters and the initial conditions are fixed. The phase diagrams with different values of the derivative order q are employed to demonstrate the behavior of system (11), as shown in Figure 2. From which it can be seen that system (11) converges to a fixed point as q=0.90, and it is period-1 for q=0.92, period-2 for q=0.95, and period-4 for q=0.955.Figure 2
The attractors for different values of the derivative orderq. (a) q=0.90; (b) q=0.92; (c) q=0.95; (d) q=0.955.
(a) (b) (c) (d)


## 4. Bifurcations of the New System

Bifurcations play a vital role in dynamics research for fractional-order systems. Therefore, in this section, bifurcation analysis is conducted to study the rich dynamical behavior of the fractional-order system (11) in the two cases of commensurate-order and incommensurate-order, respectively.


### 4.1. The Commensurate-Order Case

Firstly, with the system parameters fixed and the derivative order varying on the closed intervalq∈[0.9,1], bifurcation diagram for system (11) is depicted in Figure 3. Clearly, the evolution of the period-doubling scenario and saddle-node bifurcation can be observed from this figure. When the order satisfies q≤0.958, the route of leading to chaos for system (11) is period-doubling bifurcation. When q=0.978, the chaotic solutions disappear suddenly, and two new period-1 solutions appear, which means that the saddle-node bifurcation occurs as the order q varies.Figure 3
The bifurcation diagram of system (11) with the variation of the derivative order q.Secondly, bifurcations with the variation of the system parametera are studied for q=0.99 and b=0.02. In Figure 4, it can be seen that a series of period-doubling bifurcations occurs as the system parameter a decreases from 0.3 to 0.1. Meanwhile, the period-3 and chaos windows can also be observed from the bifurcation diagram.Figure 4
The bifurcation diagram of system (11) with the variation of the system parameter a.In order to get a better understanding of the dynamics of system (11), bifurcations with the variation of the parameter a for the different values of derivative order q are given. With the parameter a varying, the corresponding bifurcation diagrams with several specified values of the derivative order like q=0.98, q=0.97, q=0.96, q=0.95, q=0.94, and q=0.93 are plotted in Figure 5 for system (11). Compared these figures, it can be seen that the bifurcation structure of system (11) changes qualitatively with the variation of the parameter a and the order q. The area of the chaotic motion decreases with the augment of the periodic motion as the derivative order decreases. From Figures 5(a)–5(c), it is clear that the route leading to chaos is period-doubling bifurcation. Meanwhile, a typical period-doubling bifurcation, of which numerical solutions change from period-1 to period-2, can be seen clearly from Figures 5(d) and 5(e). The system becomes totally periodic as the order q decreases to 0.93; see Figure 5(f). In order to further discuss the bifurcations with the variation of both the system parameter a and derivative order q, a bifurcation diagram in the three-dimensional space is depicted in Figure 6. It can be seen that dynamics of system (11) becomes simple as the derivative order decreases from 0.98 to 0.93.Figure 5
The bifurcation diagrams of system (11) when the parameter a varied with the different values of the order q: (a) q=0.98; (b) q=0.97; (c) q=0.96; (d) q=0.95; (e) q=0.94; (f) q=0.93.
(a) (b) (c) (d) (e) (f)Figure 6
The bifurcation diagram of system (11) in three-dimensional space with the variation of both the parameter a and the order q.Thirdly, the bifurcation of system (11) with the variation of the parameter b is studied for q=0.99 and a=0.3. In Figure 7, there is a long time of the chaotic window for b≤-0.008. Meanwhile, the route leading to chaos is the period-doubling bifurcation for system (11).Figure 7
The bifurcation diagram of system (11) with the variation of the parameter b.Bifurcations with the variation of the parameterb are studied for different values of the derivative order q. When the derivative order q decreases from 0.98 to 0.93, the corresponding bifurcation diagrams are plotted for the fractional-order system (11) when the parameter b∈[-0.03,0.02]; see Figure 8, from which, it is clear that structure of dynamics of system (11) evolves as the order q varies. The interval of the chaotic motion increases and that of the periodic motion decreases. Meanwhile, system (11) is periodic completely when the order q=0.94, and period-4 and period-2 motions can be obtained from Figure 8(e). When the order q=0.93, only period-2 and period-1 motions exist; see Figure 8(f). In order to further discuss the bifurcations with the variation of both the system parameter b and derivative order q, a bifurcation diagram in the three-dimensional space is plotted in Figure 9. It can be seen that dynamics of system (11) becomes complex as the derivative order increases from 0.93 to 0.98.Figure 8
The bifurcation diagrams of system (11) when the parameter b varied with the different values of the order q: (a) q=0.98; (b) q=0.97; (c) q=0.96; (d) q=0.95; (e) q=0.94; (f) q=0.93.
(a) (b) (c) (d) (e) (f)Figure 9
The bifurcation diagram of system (11) in three-dimensional space with the variation of both the parameter b and the order q.


### 4.2. The Incommensurate-Order Case

In order to learn more about the characteristics of the new fractional-order system (11), the dynamics of system (11) with the variation of the different derivative orders will be investigated in this subsection. In the following work, the system parameters are taken as a=0.3 and b=0.02.Firstly, the bifurcation diagram with the derivative orderq1∈[0.9,1] as the others two derivative orders q2, q3 are both fixed and q2=q3=q=1 is plotted. A series of period-doubling bifurcations can be seen in Figure 10. Based on this, bifurcation diagrams versus the order q1 when the others two derivative orders decrease from 0.98 to 0.93 are given in Figure 11. For Figures 11(a)–11(c), it can be seen that the two branches of the bifurcation gradually decouple, and the area of the chaos gradually decreases with that of the period increases. For Figures 11(d)–11(f), it is clear that the fractional-order system (11) is periodic completely. A bifurcation diagram with the variation of both the order q1 and derivative order q in the three-dimensional space is plotted in Figure 12. It is clear that dynamics of system (11) becomes simple as the derivative order decreases from 0.98 to 0.93.Figure 10
The bifurcation diagram of system (11) with the variation of the order q1.Figure 11
The bifurcation diagram of system (11) when the order q1 varied with the different values of the order q: (a) q=0.98; (b) q=0.97; (c) q=0.96; (d) q=0.95; (e) q=0.94; (f) q=0.93.
(a) (b) (c) (d) (e) (f)Figure 12
The bifurcation diagram of system (11) in three-dimensional space with the variation of both the orders q1 and q.Secondly, the bifurcation diagram versus the orderq2∈[0.9,1] when the others two derivative orders q1, q3 are both fixed and q1=q3=q=1 is plotted in Figure 13. Meanwhile, a series of bifurcation diagrams versus the order q2 when the others two derivative orders decrease from 0.98 to 0.85 are given in Figure 14. By comparing Figure 14(a) and Figure 14(b), it is clear that the dynamics of system (11) becomes simple as the order q decreases. For Figures 14(b) and 14(c), it can be observed that the area of chaos increases and that of period decreases. From Figures 14(c) and 14(d), it can be seen that the two branches of the bifurcation gradually couple, and the area of the chaos gradually decreases with that of the period increases. For the rest of the figures in Figure 11, it can be seen that system (11) is totally periodic when the order q=0.90. Meanwhile, the typical period-doubling bifurcation can be seen in Figure 14(g). When the order q=0.85, only period-1 motion for system (11) exists. A bifurcation diagram with the variation of both the order q2 and derivative order q in the three-dimensional space is plotted in Figure 15. It is clear that dynamics of system (11) becomes simple as the derivative order decreases from 0.98 to 0.85.Figure 13
The bifurcation diagram of system (11) with the variation of the derivative order q2.Figure 14
The bifurcation diagram of system (11) when the order q2 varied with the different values of the order q: (a) q=0.98; (b) q=0.97; (c) q=0.96; (d) q=0.95; (e) q=0.94; (f) q=0.93; (g) q=0.90; (h) q=0.85.
(a) (b) (c) (d) (e) (f) (g) (h)Figure 15
The bifurcation diagram of system (11) in three-dimensional space with the variation of both the orders q2 and q.The bifurcation diagrams with the variation of the orderq3 when the other two orders q1, q2 decrease from 0.98 to 0.93 will not be given in here for the similarity.


## 4.1. The Commensurate-Order Case

Firstly, with the system parameters fixed and the derivative order varying on the closed intervalq∈[0.9,1], bifurcation diagram for system (11) is depicted in Figure 3. Clearly, the evolution of the period-doubling scenario and saddle-node bifurcation can be observed from this figure. When the order satisfies q≤0.958, the route of leading to chaos for system (11) is period-doubling bifurcation. When q=0.978, the chaotic solutions disappear suddenly, and two new period-1 solutions appear, which means that the saddle-node bifurcation occurs as the order q varies.Figure 3
The bifurcation diagram of system (11) with the variation of the derivative order q.Secondly, bifurcations with the variation of the system parametera are studied for q=0.99 and b=0.02. In Figure 4, it can be seen that a series of period-doubling bifurcations occurs as the system parameter a decreases from 0.3 to 0.1. Meanwhile, the period-3 and chaos windows can also be observed from the bifurcation diagram.Figure 4
The bifurcation diagram of system (11) with the variation of the system parameter a.In order to get a better understanding of the dynamics of system (11), bifurcations with the variation of the parameter a for the different values of derivative order q are given. With the parameter a varying, the corresponding bifurcation diagrams with several specified values of the derivative order like q=0.98, q=0.97, q=0.96, q=0.95, q=0.94, and q=0.93 are plotted in Figure 5 for system (11). Compared these figures, it can be seen that the bifurcation structure of system (11) changes qualitatively with the variation of the parameter a and the order q. The area of the chaotic motion decreases with the augment of the periodic motion as the derivative order decreases. From Figures 5(a)–5(c), it is clear that the route leading to chaos is period-doubling bifurcation. Meanwhile, a typical period-doubling bifurcation, of which numerical solutions change from period-1 to period-2, can be seen clearly from Figures 5(d) and 5(e). The system becomes totally periodic as the order q decreases to 0.93; see Figure 5(f). In order to further discuss the bifurcations with the variation of both the system parameter a and derivative order q, a bifurcation diagram in the three-dimensional space is depicted in Figure 6. It can be seen that dynamics of system (11) becomes simple as the derivative order decreases from 0.98 to 0.93.Figure 5
The bifurcation diagrams of system (11) when the parameter a varied with the different values of the order q: (a) q=0.98; (b) q=0.97; (c) q=0.96; (d) q=0.95; (e) q=0.94; (f) q=0.93.
(a) (b) (c) (d) (e) (f)Figure 6
The bifurcation diagram of system (11) in three-dimensional space with the variation of both the parameter a and the order q.Thirdly, the bifurcation of system (11) with the variation of the parameter b is studied for q=0.99 and a=0.3. In Figure 7, there is a long time of the chaotic window for b≤-0.008. Meanwhile, the route leading to chaos is the period-doubling bifurcation for system (11).Figure 7
The bifurcation diagram of system (11) with the variation of the parameter b.Bifurcations with the variation of the parameterb are studied for different values of the derivative order q. When the derivative order q decreases from 0.98 to 0.93, the corresponding bifurcation diagrams are plotted for the fractional-order system (11) when the parameter b∈[-0.03,0.02]; see Figure 8, from which, it is clear that structure of dynamics of system (11) evolves as the order q varies. The interval of the chaotic motion increases and that of the periodic motion decreases. Meanwhile, system (11) is periodic completely when the order q=0.94, and period-4 and period-2 motions can be obtained from Figure 8(e). When the order q=0.93, only period-2 and period-1 motions exist; see Figure 8(f). In order to further discuss the bifurcations with the variation of both the system parameter b and derivative order q, a bifurcation diagram in the three-dimensional space is plotted in Figure 9. It can be seen that dynamics of system (11) becomes complex as the derivative order increases from 0.93 to 0.98.Figure 8
The bifurcation diagrams of system (11) when the parameter b varied with the different values of the order q: (a) q=0.98; (b) q=0.97; (c) q=0.96; (d) q=0.95; (e) q=0.94; (f) q=0.93.
(a) (b) (c) (d) (e) (f)Figure 9
The bifurcation diagram of system (11) in three-dimensional space with the variation of both the parameter b and the order q.


## 4.2. The Incommensurate-Order Case

In order to learn more about the characteristics of the new fractional-order system (11), the dynamics of system (11) with the variation of the different derivative orders will be investigated in this subsection. In the following work, the system parameters are taken as a=0.3 and b=0.02.Firstly, the bifurcation diagram with the derivative orderq1∈[0.9,1] as the others two derivative orders q2, q3 are both fixed and q2=q3=q=1 is plotted. A series of period-doubling bifurcations can be seen in Figure 10. Based on this, bifurcation diagrams versus the order q1 when the others two derivative orders decrease from 0.98 to 0.93 are given in Figure 11. For Figures 11(a)–11(c), it can be seen that the two branches of the bifurcation gradually decouple, and the area of the chaos gradually decreases with that of the period increases. For Figures 11(d)–11(f), it is clear that the fractional-order system (11) is periodic completely. A bifurcation diagram with the variation of both the order q1 and derivative order q in the three-dimensional space is plotted in Figure 12. It is clear that dynamics of system (11) becomes simple as the derivative order decreases from 0.98 to 0.93.Figure 10
The bifurcation diagram of system (11) with the variation of the order q1.Figure 11
The bifurcation diagram of system (11) when the order q1 varied with the different values of the order q: (a) q=0.98; (b) q=0.97; (c) q=0.96; (d) q=0.95; (e) q=0.94; (f) q=0.93.
(a) (b) (c) (d) (e) (f)Figure 12
The bifurcation diagram of system (11) in three-dimensional space with the variation of both the orders q1 and q.Secondly, the bifurcation diagram versus the orderq2∈[0.9,1] when the others two derivative orders q1, q3 are both fixed and q1=q3=q=1 is plotted in Figure 13. Meanwhile, a series of bifurcation diagrams versus the order q2 when the others two derivative orders decrease from 0.98 to 0.85 are given in Figure 14. By comparing Figure 14(a) and Figure 14(b), it is clear that the dynamics of system (11) becomes simple as the order q decreases. For Figures 14(b) and 14(c), it can be observed that the area of chaos increases and that of period decreases. From Figures 14(c) and 14(d), it can be seen that the two branches of the bifurcation gradually couple, and the area of the chaos gradually decreases with that of the period increases. For the rest of the figures in Figure 11, it can be seen that system (11) is totally periodic when the order q=0.90. Meanwhile, the typical period-doubling bifurcation can be seen in Figure 14(g). When the order q=0.85, only period-1 motion for system (11) exists. A bifurcation diagram with the variation of both the order q2 and derivative order q in the three-dimensional space is plotted in Figure 15. It is clear that dynamics of system (11) becomes simple as the derivative order decreases from 0.98 to 0.85.Figure 13
The bifurcation diagram of system (11) with the variation of the derivative order q2.Figure 14
The bifurcation diagram of system (11) when the order q2 varied with the different values of the order q: (a) q=0.98; (b) q=0.97; (c) q=0.96; (d) q=0.95; (e) q=0.94; (f) q=0.93; (g) q=0.90; (h) q=0.85.
(a) (b) (c) (d) (e) (f) (g) (h)Figure 15
The bifurcation diagram of system (11) in three-dimensional space with the variation of both the orders q2 and q.The bifurcation diagrams with the variation of the orderq3 when the other two orders q1, q2 decrease from 0.98 to 0.93 will not be given in here for the similarity.


## 5. Conclusions

In this paper, a new fractional-order system is presented. Firstly, the stability of the equilibriums is analyzed. Based on the stability analysis, the reason of the generation for the one scroll of the attractor is analyzed. Phase diagrams for the different values of the derivative order are obtained to show the rich dynamics of the new fractional-order system.Bifurcations of the new fractional-order system in commensurate-order and incommensurate-order cases are studied in detail. In the commensurate-order case, when the derivative order decreases from 0.99, bifurcations with the variation of a system parameter are investigated. In the incommensurate-order case, bifurcations with the variation of a derivative order when the other orders decreases from 1 are analyzed. Period-doubling and saddle-node bifurcations can be observed.What is more, it can be concluded that the dynamics of the new fractional-order system becomes periodic or simple when the derivative order approaches to 0 and chaotic or complex when the order approaches to 1 from a global perspective. These results obtained in this paper can be referenced for the bifurcation control of fractional-order systems. The generalization of the conclusion for other fractional-order systems will be our future work.



---

*Source: 8341514-2019-01-01.xml*