In this paper, we propose two control strategies for the Q–S synchronization of fractional-order discrete-time chaotic systems. Assuming that the dimension of the response system m is higher than that of the drive system n, the first control scheme achieves n-dimensional synchronization whereas the second deals with the m-dimensional case. The stability of the proposed schemes is established by means of the linearization method. Numerical results are presented to confirm the findings of the study.
1. Introduction
Since the concept of chaos synchronization was first proposed nearly three decades ago, dynamical systems exhibiting chaotic behavior have attracted considerable attention in a number of fields. In addition to the original continuous-time chaotic systems such as Lorenz and Chen systems, discrete-time chaotic systems or maps have been proposed and studied. Among the most widely used 2-component discrete systems of these are the Hénon map [1], the Lozi map [2], and the flow model [3]. Higher-dimensional systems have also been studied such as the generalized Hénon map [4] and the Stefanski map [5]. All of these maps are classified as integer–order signifying that the difference order is integer.
Although fractional calculus was studied by mathematicians as early as the nineteenth century, it wasn’t until recently that it gained popularity in applied science and engineering. For the longest time, the study and application of fractional calculus was limited to continuous time. However, most recently, researchers have diverted their attention to the discrete-time case and attempted to put together a complete theoretical framework for the subject. Perhaps one of the earliest works is that of [6]. Among the most interesting and relevant works on fractional discrete calculus in the last decade is that of [7] where the authors introduce a backward fractional difference operator. In [8], the author discusses the discrete-time difference counterparts of conventional Riemann and Caputo derivatives. More on the general notation of fractional discrete-time calculus can be found in [9, 10]. Furthermore, the numerical formulas corresponding to a fractional difference systems can be found in [11]. Recent studies have examined the stability conditions for fractional discrete-time systems including [12]. Most recently, some advances have been made in relation to chaotic fractional discrete-time systems and their applications. A number of fractional difference maps including the fractional logistic map [13], the fractional sine map [14], fractional Hénon map [15], and fractional cubic logistic map [16].
Chaos synchronization in its simplest form refers to the control of a response chaotic system to force its trajectory towards that of a drive [17]. In this case, the error defined as the difference between the states of the drive and those of the response towards zero as t→∞. Numerous other types of synchronization have been proposed in the literature, whereby a general function of the two state vectors is forced to zero instead. The synchronization of integer–order discrete-time systems subject to different definitions of the error has been studied in the literature such as hybrid synchronization [18], generalized synchronization [19], and inverse full state hybrid projective synchronization [20]. As for fractional-order discrete systems, the literature related to their synchronization is still very limited including synchronization of the logistic map [21], synchronization based on the stability condition [22], synchronization of linearly coupled fractional Hénon maps [23], and impulsive synchronization of fractional-order discrete-time hyperchaotic systems [24].
In this paper, we are concerned with the Q–S synchronization of fractional-order discrete-time chaotic systems with different dimensions and non–identical orders. In this scheme, two functions Q and S are used to condition the response and drive states, respectively, such that different values for Q and S lead to different types of synchronization. Q–S synchronization was first proposed by Yan in 2005 [25] for continuous–time dynamical systems. The author developed a backstepping approach with a strict feedback form to synchronize two identical systems. The author, then, extended the Q–S scheme to discrete-time systems [26]. In the years that followed, several algorithms were proposed for the Q–S synchronization of integer-order continuous-time systems including [27–29]. Several studies also looked at the fractional continuous-time case such as [30]. Finally, in [31, 32], the authors investigated Q–S synchronization between integer-order and fractional-order continuous-time systems with different dimensions.
In the next section of this paper, we will present the notation adopted in our study and formulate the Q–S synchronization problem for a fractional-order discrete-time chaotic drive-response pair. We denote the number of states for the drive and response by n and m, respectively. Sections 3 and 4 present the control laws for the m–dimensional and n–dimensional cases, respectively. Section 5 presents some numerical results related to two particular examples. Finally, Section 6 provides a general summary of the main findings of this study.
2. Problem Formulation
Consider the following drive and response maps described as(1)ΔaαCXt=AXt+υ-1+fXt+υ-1,∀t∈Na+1-α,ΔaβCYt=BYt+υ-1+gYt+υ-1+U,∀t∈Na+1-βwhere Xt=x1t,…,xntT and Yt=y1t,…,ymtT denote the drive and response state vectors, respectively, A∈Rn×n, B∈Rm×m, f:Rn→Rn, g:Rm→Rm are nonlinear functions and U=ui1≤i≤m is a vector controller to be determined by the control laws of the synchronization scheme. Note that Na denotes the set of natural numbers starting from a and 0<υ≤1.
The notation ΔaυCXt denotes the Caputo type delta difference of Xt defined over Na [8]. This operator can be defined as(2)ΔaυCXt=Δa-n-υΔnXt=1Γn-υ∑s=at-n-υt-s-1n-υ-1ΔsnXs,for υ∉N, t∈Na+n-υ, and n=υ+1. The operator ΔsnXt denotes the υ–th fractional sum of xi:Na→R defined in [7] as(3)Δa-υXt=1Γυ∑s=at-υt-s-1υ-1Xs,with υ>0, Na=a,a+1,a+2,…, and tυ is given in terms of the Gamma function Γ as(4)tυ=Γt+1Γt+1-υ.
For the purpose of this study, we will assume a=0 and thus the subscript may be ignored. With these notations in mind, we can go ahead and define the Q-S synchronization of our drive-response pair. The following definition is consistent with the literature reviewed earlier in this paper.
Definition 1.
The drive-response pair (1) is Q-S synchronized in the dth dimension if there exist a controller U=ui1≤i≤m and two functions Q:Rm→Rd and S:Rn→Rd such that(5)limt→+∞et≔QYt-SXt=0.
It is easy to observe from Definition 1 that, depending on our choice of the pair QYt,SXt, we may end up with one of many synchronization types, namely, (6)Yt,Xt⟶Complete synchronization,Yt,-Xt⟶Anti synchronization,Yt,θXt,θ∈R∗⟶Projective synchronization,Yt,MXt,M∈Rm×n⟶Matrix projective synchronization,MYt,Xt,M∈Rn×m⟶Inverse matrix projective synchronization,Yt,SXt⟶Generalized synchronization,QYt,Xt⟶Inverse generalized synchronization.
For reasons that will become clear later on, suppose that Q can be divided into two parts:(7)QYt=QYt+qYt,where Q is an invertible d×d matrix and q:Rd→Rd is a nonlinear function. In the following two sections, we will present two distinct synchronization schemes of dimensions d=m and d=n, respectively.
3. Q-S Synchronization in Dimension m
The Q–S synchronization error defined in (5) can be written in the following form:(8)et=QYt+qYt-SXt,where Q is an invertible m×m matrix, q:Rm→Rm, and S:Rn→Rm. The β order Caputo fractional difference of (8) is of the following form:(9)ΔβCet=ΔβCQYt+qYt-SXt=ΔβCQYt+ΔβCqYt-ΔβCSXt=QΔβCYt+ΔβCqYt-ΔβCSXt=QBYt+β-1+gYt+β-1+U+ΔβCqYt-ΔβCSXt=QBYt+β-1+QgYt+β-1+Q×U+ΔβCqYt-ΔβCSXt.This can be further simplified to(10)ΔβCet=B-Cet+β-1+R+Q×U,where C∈Rm×m is an appropriately chosen control matrix and(11)R=C-Bet+QBYt+QgYt+ΔβCqYt-ΔβCSXt.The following remark is important.
Remark 2.
Note that our aim in synchronization is to find a suitable controller U that forces the synchronization error to zero asymptotically. Since U is part of the response system, the function QYt hinders our process. To overcome this problem, we assumed that QYt can be decomposed into a linear part (matrix Q) and some other nonlinear function qYt. The reasoning behind this is that matrix Q is constant and thus may be factored out of the fractional difference operator.
Before stating the proposed control law and establishing its stability, it is important to state the following theorem, which is essential for our proof. Interested readers are referred to [33] for the proof of this result.
Theorem 3.
The zero equilibrium of the linear fractional-order discrete-time system:(12)ΔaυCet=Met+υ-1,where e(t)=e1(t),…,en(t)T, 0<υ≤1, M∈Rn×n and ∀t∈Na+1-υ, is asymptotically stable, if(13)λ∈z∈C:z<2cosargz-π2-υυandargz>υπ2,for all the eigenvalues λ of M.
With this stability result, we are now ready to present our synchronization scheme.
Theorem 4.
The drive-response pair (1) is globally Q–S synchronized in dimension m subject to(14)U=-Q-1×R,where Q-1 is the inverse of Q and the control matrix C is selected such that all the eigenvalues λ of B-C satisfy(15)-2β<λ<0.
Proof.
By substituting (14) into (10), we obtain a new formulation for the error system:(16)ΔβCet=B-Cet+β-1.Now, it is easy to see that subject to (15), all eigenvalues λ of matrix B-C satisfy (17)argλ=π>βπ2,and(18)λ<2cosargλ-π2-ββ.Therefore, by means of Theorem 3, we establish the global asymptotic stability of the zero solution. Consequently, the drive-response maps (1) are globally Q-S synchronized in dimension m.
4. Q-S Synchronization in Dimension n
Let us now move to establish Q-S synchronization in dimension n. In this section, we assume n<m. We define the following notations: (19)Q∈Rn×m,q:Rm→Rn,andS:Rn→Rn.We set the last m-n controller components to zero, i.e., (20)ui=0for n<i≤m.Then, the error system reduces to(21)ΔβCet=QBYt+β-1+QgYt+β-1+Q^×U^+ΔβCqYt-ΔβCSXt,where Q^ is an invertible n×n matrix and U^=u1,…,unT. This can be further simplified to(22)ΔCβet=A-Let+T+Q^×U^,where L∈Rn×n is an appropriate control matrix and(23)T=L-Aet+QBYt+QgYt+ΔβCqYt-ΔβCSXt.
To achieve Q–S synchronization between the drive-response maps (1) in dimension n, we propose choosing(24)U^=-Q^-1T.Substituting (24) into (22) yields the error dynamics(25)ΔβCet=A-Let+β-1.The following theorem can be proven in the exact same way as Theorem 4.
Theorem 5.
By selecting control matrix L such that all the eigenvalues λ of matrix A-L satisfy -2β<λ<0, the drive and response maps (1) are globally Q–S synchronized in n dimensions subject to (20) and (24).
5. Numerical Examples
In order to show the validity of the proposed control strategies, let us consider some numerical examples. The 2D fractional Hénon map given by(26)ΔαCx1t=x2t+α-1+1-a1x12t+α-1-x1t+α-1,ΔαCx2t=b1x1t+α-1-x2t+α-1,is chosen as the driving map. The fractional Hénon map (26) was proposed in [34]. The authors showed that this map has a chaotic attractor, for instance, when a1,b1=1.4,0.3 and α=0.984. In order to unify the notation, (26) can be rewritten as(27)ΔαCXt=AXt+fXt,where(28)A=-11b1-1andf=1-a1x12t0.The chaotic behavior of the drive map is depicted in Figure 1.
Phase space plot for the fractional Hénon map with a1,b1=1.4,0.3, α=0.984, and x10=x20=0.
As for the response map, we choose the 3D fractional-order generalized Hénon map [35], which is of the following form:(29)ΔβCy1t=-y1t+β-1-b2y3t+β-1+u1t+β-1,ΔβCy2t=b2y3t+β-1+y1t+β-1-y2t+β-1+u2t+β-1,ΔβCy3t=1+y2t+β-1-a2y32t+β-1-y3t+β-1+u3t+β-1.This system exhibits a chaotic behavior subject to certain conditions. We choose as an example the parameters a2,b2=0.99,0.2 and β=0.984, which has chaotic trajectories as shown in Figure 2. Map (29) can be rewritten as(30)ΔβCYt=BYt+gYt,where(31)B=-10-b21-1b201-1andg=001-a2y32t.
Phase portraits for the fractional generalized Hénon map with a2,b2=0.99,0.2, β=0.984 and y10,y20,y30=0.1,0.2,0.5.
We will take 3- and 2-dimensional Q–S synchronization cases separately following the control laws of Theorems 4 and 5, respectively.
Case 1.
Let us, first, consider the 3-dimensional Q–S synchronization case. Based on the approach stated in Section 3, the error system is given by(32)e1t,e2t,e3tT=Qy1t,y2t,y3t-Sx1t,x2t,where (33)Qy1t,y2t,y3t=y1t,y1t+y2t,y3t3+y2tT,and (34)Sx1t,x2t=x1t,x2t,x1t+x2t.We may rewrite Q in the following form:(35)Qy1t,y2t,y3t=Q×y1t,y2t,y3tT+qy1t,y2t,y3twhere(36)Q=100120003and(37)qy1t,y2t,y3t=0,0,y2ty3tT.It is easy to see that(38)Q-1=100-121200013.According to our approach described in Section 3, there exists a control matrix C such that all the eigenvalues B-C satisfy the condition of Theorem 4. We may simply choose the following:(39)C=000100010.We start by calculating R by means of (11) and then U by means of (14). The resulting error system is given by(40)ΔβCe1t=-e1t+β-1-b2e3t+β-1,ΔβCe2t=-e2t+β-1+b2e3t+β-1,ΔβCe3t=-e3t+β-1.The time evolution of the errors is depicted in Figure 3 for the initial values: (41)e10e20e30=0.10.31.6.The errors clearly converge to zero in sufficient time indicating that the drive map (26) and the response map (29) become Q–S synchronized in 3 dimensions.
The evolution of errors over time for Case 1.
Case 2.
As for the two-dimensional case, we follow the approach described in Section 4. The error system is given by(42)e1t,e2t,e3tT=Qy1t,y2t,y3t-Sx1t,x2t,where (43)Qy1t,y2t,y3t=y1t+2y2t+y32t,-y1t+3y3tT,and(44)Sx1t,x2t=x1t,x1tx2t.We can write(45)Qy1t,y2t,y3t=Q×y1t,y2t,y3tT+qy1t,y2t,y3t,where(46)Q=120-103,and(47)qy1t,y2t,y3t=y32t,0T.This yields(48)Q^=12-10andQ^-1=12-10.According to Theorem 5, the control matrix L can be selected as(49)L=01-1+b1-12,which clearly satisfies the condition of Theorem 5. Next, we calculate T as specified in (23) and U^ as in (24). The complete control vector U is formed by appending a zero to the end of U^. The resulting error system is given by(50)ΔβCe1t=-e1t+β-1,ΔβCe2t=e1t+β-1-12e2t+β-1.The convergence of the errors to zero is depicted in Figure 4 for the initial conditions: (51)e10e20=0.751.14.
The evolution of errors over time for Case 2.
6. Concluding Remarks
This paper has studied Q–S synchronization of fractional-order discrete-time chaotic systems of different dimensions. The Q–S scheme aims for the general definition of the error towards zero in finite time, thereby covering a range of different synchronization types. We assume an n-dimensional drive map and an m-dimensional response map with the condition m>n. Two different control laws have been derived and the asymptotic stability of their zero solution investigated through the linearization method. The two strategies correspond to the synchronization dimensions d=n and d=m.
Computer simulation results have been presented, whereby a 3-dimensional fractional-order generalized Hénon map response synchronized the 2-dimensional fractional-order Hénon drive. The two proposed schemes were utilized to achieve 2- and 3-dimensional synchronization. The errors have been shown to converge towards zero in sufficient time.
Data Availability
No external data has been used in this manuscript.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
HénonM.A two-dimensional mapping with a strange attractor1976501697710.1007/BF01608556MR0422932Zbl0576.58018LOZIR.Un atracteur étrange du type attracteur de H énon19783991010.1051/jphyscol:1978505ItohM.YangT.ChuaL. O.Conditions for impulsive synchronization of chaotic and hyperchaotic systems200111255156010.1142/S0218127401002262MR1830352Zbl1090.375202-s2.0-0035587560HitzlD. L.ZeleF.An exploration of the H\'enon quadratic map198514330532610.1016/0167-2789(85)90092-2MR793708Zbl0602.58025StefańskiK.Modelling Chaos and Hyperchaos with 3-D Maps199891839310.1016/S0960-0779(97)00051-92-s2.0-0031700210Zbl0934.37033MillerK. S.RossB.Fractional difference calculusProceedings of the Int. Symp. Univalent Func., Frac. Calc. Appl.1989Horwood, Chichester139152Ellis Horwood Ser. Math.AticiF. M.EloeP. W.Discrete fractional calculus with the nabla operator20093112MR2558828Zbl1189.39004AbdeljawadT.On Riemann and Caputo fractional differences20116231602161110.1016/j.camwa.2011.03.036MR2824747Zbl1228.260082-s2.0-79960975255HolmM.2011Lincoln, NebraskaUniversity of Nebraska-LincolnMR2873503GoodrichC.PetersonA. C.2015New York, NY, USASpringer International10.1007/978-3-319-25562-0Zbl1350.39001AnastassiouG. A.Principles of delta fractional calculus on time scales and inequalities2010523-455656610.1016/j.mcm.2010.03.055MR2658507Zbl1201.260012-s2.0-77954861707WyrwasM.MozyrskaD.GirejkoE.Stability of discrete fractional-order nonlinear systems with the nabla Caputo differenceProceedings of the 6th Workshop on Fractional Differentiation and Its Applications, FDA 2013February 2013France1671712-s2.0-84881031485WuG.BaleanuD.Discrete fractional logistic map and its chaos2014751-228328710.1007/s11071-013-1065-7MR3144852Zbl1281.34121WuG. C.BaleanuD.ZengS. D.Discrete chaos in fractional sine and standard maps20143785-648448710.1016/j.physleta.2013.12.010MR3149702Zbl1331.37048LiuY.Discrete Chaos in Fractional Henon Maps2014183170175MR3297306LiuC. X.HuangL. L.WuK. T.Chaos in discrete fractional cubic logistic map and bifurcation analysis201512105111PecoraL. M.CarrollT. L.Synchronization in chaotic systems199064882182410.1103/PhysRevLett.64.821MR1038263Zbl0938.370192-s2.0-0343689904OuannasA.AzarA. T.Abu-SarisR.A new type of hybrid synchronization between arbitrary hyperchaotic maps201786188718942-s2.0-8503238466410.1007/s13042-016-0566-3OuannasA.OdibatZ.Generalized synchronization of different dimensional chaotic dynamical systems in discrete time2015811-276577110.1007/s11071-015-2026-0MR3355066Zbl1347.34086OuannasA.GrassiG.Inverse full state hybrid projective synchronization for chaotic maps with different dimensions20162592-s2.0-8498838061410.1088/1674-1056/25/9/090503090503WuG. C.BaleanuD.Chaos synchronization of the discrete fractional logistic map2014102969910.1016/j.sigpro.2014.02.022WuG.-C.BaleanuD.XieH.-P.ChenF.-L.Chaos synchronization of fractional chaotic maps based on the stability condition201646037438310.1016/j.physa.2016.05.045MR3511762LiuY.Chaotic synchronization between linearly coupled discrete fractional Hénon maps20169033133172-s2.0-8495651507110.1007/s12648-015-0742-4MegherbiO.HamicheH.DjennouneS. d.BettayebM.A new contribution for the impulsive synchronization of fractional-order discrete-time chaotic systems20179031519153310.1007/s11071-017-3743-3MR3713942Zbl1380.37072YanZ.Q-S (lag or anticipated) synchronization backstepping scheme in a class of continuous–time hyperchaotic systems: a symbolic-numeric computation approach200515210.1063/1.1827431MR2133445YanZ.Q-S synchronization in 3D Hénon-like map and generalized Hénon map via a scalar controller2005342430931710.1016/j.physleta.2005.04.049Zbl1222.370932-s2.0-21244494302ManfengH.ZhenyuanX.A general scheme for Q-S synchronization of chaotic systems20086941091109910.1016/j.na.2007.06.038MR2426675Zbl1153.343382-s2.0-44649156300ZhaoJ.RenT.Q–S synchronization between chaotic systems with double scaling functions201062366567210.1007/s11071-010-9752-0Zbl1209.93080ZhaoJ.ZhangK.A general scheme for Q-S synchronization of chaotic systems with unknown parameters and scaling functions201021672050205710.1016/j.amc.2010.03.035MR2647073Zbl1192.651532-s2.0-77953137579OuannasA.AzarA. T.VaidyanathanS.On a simple approach for Q-S synchronisation of chaotic dynamical systems in continuous-time201781202710.1504/IJCSM.2017.083167MR3650805OuannasA.Abu-SarisR.A robust control method for QS synchronization between different dimensional integer-order and fractional-order chaotic systems201517OuannasA.AzarA. T.ZiarT.VaidyanathanS.Q-S chaos synchronization between fractional–order master and integer–order slave systemsProceedings of the 15th Int. Multi-Conf. Sys. Sig. Dev.2018CermákJ.GyõriI.NechvátalL. e.On explicit stability conditions for a linear fractional difference system201518365167210.1515/fca-2015-0040MR3351493HuT.Discrete Chaos in Fractional Henon Map2014522432248ShuklaM. K.SharmaB. B.Investigation of chaos in fractional order generalized hyperchaotic Henon map20177826527310.1016/j.aeue.2017.05.009