This paper investigates the modified function projective synchronization between fractional-order chaotic systems, which are partially linear financial systems with uncertain parameters. Based on the stability theory of fractional-order systems and the Lyapunov matrix equation, a controller is obtained for the synchronization between fractional-order financial chaotic systems. Using the controller, the error systems converged to zero as time tends to infinity, and the uncertain parameters were also estimated so that the phenomenon of parameter distortion was effectively avoided. Numerical simulations demonstrate the validity and feasibility of the proposed method.
National Natural Science Foundation of China71232004Fundamental Research Funds for the Central Universities of ChinaCDJSK10011. Introduction
In recent years, study on chaos is one of the most interesting research topics in real and physical systems [1–4]. In 1985, chaotic behavior was discovered in financial systems [5]. In fact, chaotic phenomena often appear in financial systems [6], such as the supernormal shock of the international oil prices, stock market crash, and financial crisis. Owing to its randomness and unpredictability, these chaotic phenomena usually lead to financial crisis and instability in the financial system. Since financial risks come from uncertainties, it is important to introduce unknown parameters in financial systems [7, 8]. Therefore, chaotic behaviors and uncertainties in financial systems must be taken into account [9]. It is necessary to investigate chaos control strategies for financial systems to address financial crisis and other related problems. Hence, it is important to achieve a synchronized and healthy development of financial markets.
Since Pecora and Carroll [10] proposed synchronization between two chaotic systems in 1990, chaotic synchronization has been extensively and intensively studied in various fields [11, 12]. Different types of chaotic synchronization have been reported, including generalized synchronization [13], lag synchronization [14], function projective synchronization [15], modified projective synchronization [16], and modified function projective synchronization (MFPS) [17, 18]. The modified function projective synchronization (MFPS) means that the drive and response systems could be synchronized up to a scaling function matrix, but not a constant matrix. Recent studies on complicated financial systems have shown satisfactory results using the nonlinear method [19]. Chaos synchronization in fractional-order financial systems has also been studied in recent years [20–22].
The synchronization of financial systems indicates that the financial systems in two different areas are to maintain synchronized development by the appropriate control conditions. That is, the drive system and the response system can be interpreted as a virtual economic target and a controlled objective, respectively. Hence, the goal of synchronization is to control a specific objective to a virtual economic target using the proposed controller. The dynamical behaviors of a financial system are more complex. Therefore, financial systems are not always completely synchronized, and some complex synchronization methods [20–22], such as function projective synchronization and MFPS, should be considered. Meanwhile, because of the complexity of the financial environment, this study highlights its theoretical value and ensures the comprehensiveness of the theory.
Also, in practice, the parameter distortion phenomenon frequently appears in fractional-order chaotic systems. In other words, parameters may be uncertain or drift with time in the chaotic synchronization between the response and drive systems because of various kinds of interferences. Many projective synchronization methods of chaotic systems with unknown parameters have been proposed in recent years [23]. Hence, achieving synchronization and identifying parameters are very important in financial chaotic systems.
In this study, we investigated fractional-order financial systems proposed by Chen [24]. A new controller is proposed to achieve synchronization. The controller makes the drive and response systems remain asymptotically stable. Based on the stability theory of fractional-order systems and the Lyapunov matrix equation, the MFPS was realized for partially linear and fractional-order (PLFO) financial chaotic systems, and the unknown parameters were also estimated so that the phenomenon of parameter distortion was avoided. Numerical simulation results showed that the proposed method effectively eliminated chaos and stabilized two financial systems.
The rest of the paper is organized as follows. In Section 2, the definition of MFPS, the PLFO chaotic system, and a useful lemma are presented. In Section 3, the scheme of MFPS is presented for PLFO chaotic systems with uncertain parameters. Numerical simulations are presented in Section 4. Finally, a conclusion is drawn in Section 5.
2. Preliminaries
The definitions of PLFO chaotic system and coupled PLFO chaotic system were given in [25]. We will give the definition of PLFO chaotic system with uncertain parameters and coupled PLFO chaotic system with uncertain parameters as follows.
Definition 1.
The fractional-order chaotic systems with uncertain parameters are defined as follows: (1)dqudtq=Mz,θ·u,dqudtq=fu,z.It is called an uncertain PLFO chaotic system, where the state vector is u=(u1,u2,…,un)T∈Rn, z∈R is a variable, and f:R×Rn→R is a differentiable function. Its order is subject to 0<q≤1. θ∈Rn represents the vector of uncertain parameters. The coefficient matrix M(z,θ)∈Rn×n is dependent on the variable z and the uncertain parameter θ.
System (1) can be written as(2)dqu1dtq=m11z,θu1+m12z,θu2+⋯+m1nz,θun,dqu2dtq=m21z,θu1+m22z,θu2+⋯+m2nz,θun,⋮dqundtq=mn1z,θu1+mn2z,θu2+⋯+mnnz,θun,dqzdtq=fu,z,where mij(z,θ) is a coefficient of uj in the ith differential equation, i,j=1,2,…,n. Hence, M(z,θ) can be described by(3)Mz,θ=m11z,θm12z,θ⋯m1nz,θm21z,θm22z,θ⋯m2nz,θ⋮⋮⋮⋮mn1z,θmn2z,θ⋯mnnz,θ.Obviously, system (1) is an uncertain partially linear and integer-order chaotic system with q=1.
Definition 2.
Consider two uncertain PLFO chaotic systems (1), which can be described as (4)dqurdtq=Mz,θ·ur,dqzdtq=fur,z,dqusdtq=M~z,θ~·us+ψ.The system is coupled by the variable z. The subscripts r and s in system (4) represent the drive and response systems, respectively. ψ∈Rn is a controller. θ~ is the estimated value of unknown parameters. The coefficient matrix M~(z,θ~)∈Rn×n is dependent on z and θ~. System (4) is the coupled PLFO chaotic system with uncertain parameters [26].
The synchronization error vector of the MFPS is defined as (5)et=us-ht·ur,where h(t)=diagh1t,h2t,…,hnt, hi(t) are continuously differentiable functions, and h(t) is a scaling function matrix. e(t)=(e1(t),e2(t),…,en(t))T∈Rn is an error state vector.
Definition 3.
System (4) is said to be MFPS with respect to the scaling function matrix h(t) if there exists a vector controller ψ such that (6)limt→∞et=limt→∞us-ht·ur=0,where · is the Euclidean norm. It implies that the MFPS between the drive and response system (4) can be achieved or the error dynamic systems are globally asymptotically stable.
Remark 4.
If h1(t)=h2(t)=⋯=hn(t), the MFPS converts to function projective synchronization. If h1(t),h2(t),…,hn(t) represent identical or different constants, the MFPS is simplified to modified projective synchronization (such as antisynchronization and complete synchronization). If h1(t)=h2(t)=⋯=hn(t)=0, the synchronization problem is transformed into a control problem.
Lemma 5 (see [27]).
Presume that the fractional-order autonomous system is(7)dqedtq=Ae·e,e0=e0,where qi∈(0,1](i=1,2,…,n) are the orders of the fractional derivative and q=(q1,q2,…,qn)T, A(e)∈Rn×n is a coefficient matrix which is a polynomial matrix and depends on the state vector e, and e=(e1,e2,…,en)T∈Rn is an n-dimension state vector. System (7) is asymptotically stable if and only if (8)argλ≥πα2,α=maxq1,q2,…,qn,where λ is an arbitrary eigenvalue of A(e). In this case, each state vector decays toward 0, such as t-α. Furthermore, the system is stable if and only if arg(λ)≥πα/2 and the critical eigenvalues that satisfy arg(λ)=πα/2 have geometric multiplicity one.
Then, the MFPS between the drive and response system (4) is transformed into the analysis of the asymptotical stability of zero solution of the error system (7).
Theorem 6 (see [28]).
Given the autonomous system (7), there exist a real symmetric positive definite matrix P and a positive definite matrix Q such that the MFPS between the drive and response system (4) can be achieved if PA(e)+(A(e))HP=-Q, where H is the conjugate transpose of a matrix and PA(e)+(A(e))HP=-Q is called a continuous Lyapunov matrix equation.
Proof.
Assume that λ is one of the eigenvalues of the polynomial matrix A(e) and the corresponding nonzero eigenvector is β; that is,(9)Aeβ=λβ.For (9), the Hermitian transpose is (10)AeβT¯=λ¯βH.Multiplying the left side of (9) by βHP, we can obtain (11)βHPAeβ=λβHPβ.Then, multiplying the right side of (10) by Pβ, we can also obtain(12)βHAeHPβ=λ¯βHPβ.From (11) and (12), we obtain (13)βHPAe+AeHPβ=λ+λ¯βHPβ.Since (PA(e)+(A(e))HP)H=PA(e)+(A(e))HP, then PA(e)+(A(e))HP is a Hermitian matrix. At the same time, since PA(e)+(A(e))HP=-Q, P is a real symmetric positive definite matrix, and Q is a positive definite matrix, then βHPβ>0 and βH(-Q)β<0, so that(14)λ+λ¯=βH-QββHPβ<0.From (14), we can obtain arg(λ)>π/2≥πα/2(α≤1). According to Lemma 5, the equilibrium point of system (7) is asymptotically stable. That is,(15)limt→∞et=limt→∞us-ht·ur=0,which indicates that the MFPS between the drive and response system (4) can be achieved. The proof is completed.
Remark 7.
Since A(e) is a polynomial matrix, it is difficult to solve for the eigenvalues. However, based on Theorem 6, the equilibrium point in system (7) is asymptotically stable if and only if PA(e)+(A(e))HP=-Q is a continuous Lyapunov matrix equation.
3. Designing the Lyapunov Matrix Equation
Next, we introduce a scheme of MFPS for coupled PLFO chaotic systems with uncertain parameters (7).
Theorem 8.
For a given scaling function matrix h(t), we can select a controller (16) and a synchronization error system (17) as follows:(16)ψ=-htM~z,θ~ur-Mz,θe+Ke+dqhturdtq,(17)dqedtq=M~z,θ~-Mz,θ+Ke,where K=diag{k1,k2,…,kn}, ki<0, 0<q≤1. The error system (17) can be equivalently written as(18)dqe1dtq=A11eθ11+A12eθ12+⋯+A1l1eθ1l1+k1e1,dqe2dtq=A2l1+1eθ2l1+1+A2l1+2eθ2l1+2+⋯+A2l2eθ2l2+k2e2,⋮dqendtq=Anln-1+1eθnln-1+1+Anln-1+2eθnln-1+2+⋯+Anlneθnln+knen,where eθij=θ~ij-θij (i=1,2,…,n, j=1,2,…,l1,l1+1,…,ln-1,ln-1+1,…,ln, where l1 stands for the number of unknown parameters in the first error equation, and ln-ln-1 stands for the number of unknown parameters in the nth error equation) and Aij is a coefficient of eθij in the ith error equation. Therefore, we can obtain the systems with uncertain parameters as (19)dqeθ11dtq=-A11e1+k1′eθ11,dqeθ12dtq=-A12e1+k2′eθ12,⋮dqeθ1l1dtq=-A1l1e1+kl1′eθ1l1,dqeθ2l1+1dtq=-A2l1+1e2+kl1+1′eθ2l1+1,⋮dqeθ2l2dtq=-A2l2e2+kl2′eθ2l2,⋮dqeθnlndtq=-Anlnen+kln′eθnln,where kj′<0(j=1,2,…,ln).We can achieve the MFPS of system (7) through the error system of (18) and the estimate of the system with unknown parameters (19). That is, the error system and the estimate of unknown parameters are asymptotically stable.
Proof.
Assume that e(t)=us-h(t)·ur; we have(20)dqetdtq=dqusdtq-dqhturdtq=M~z,θ~us+ψ-dqhturdtq.We select ψ=-h(t)M~(z,θ~)ur-M(z,θ)e+Ke+dq(h(t)ur)/dtq; then it has(21)dqetdtq=M~z,θ~us+ψ-dqhturdtq=M~z,θ~us-htM~z,θ~ur-Mz,θe+Ke+dqhturdtq-dqhturdtq=M~z,θ~-Mz,θ+Ke.Hence, we obtain a controller (16) and a synchronization error system (17).
From (18) and (19), the error system and the estimate of unknown parameters can be simply written as (22)dqedtqdqeθdtq=Aeeθ,where e=(e1,e2,…,en)T, eθ=(eθ11,eθ12,…,eθ1l1,eθ2(l1+1),…,eθ2l2,…,eθnln)T, and (23)A=k10⋯0A11⋯A1l10⋯0⋯0⋯00k2⋯00⋯0A2l1+1⋯A2l2⋯0⋯0⋮⋮⋱⋮⋮⋮⋮⋮⋮⋮⋯⋮⋮⋮00⋯kn0⋯00⋯0⋮Anln-1+1⋯Anln-A110⋯0k1′⋯00⋯0⋯0⋮0⋮⋮⋮⋮⋮⋱⋮⋮⋮⋮⋯⋮⋯⋮-A1l10⋯00⋯kl1′0⋯0⋯0⋯00-A2l1+1⋯00⋯0kl1+1′⋯0⋯0⋯0⋮⋮⋮⋮⋮⋮⋮⋮⋱⋮⋮⋮⋮⋮0-A2l2⋯00⋯00⋯kl2′⋯0⋯0⋮⋮⋯⋮⋮⋯⋮⋮⋯⋮⋱⋮⋯⋮00⋮-Anln-1+10⋮00⋯0⋯kln-1+1′⋯0⋮⋮⋯⋮⋮⋯⋮⋮⋮⋮⋯⋮⋱⋮00⋯-Anln0⋯00⋯0⋯0⋯kln′.Hence, we can obtain PA+AHP=-Q, where P is a real symmetric positive definite matrix and Q is a positive definite matrix. We usually choose P=E, where E is an identity matrix. Then, (24)Q=diag-2k1,-2k2,…,-2kn,-2k1′,…,-2kln′.According to Theorem 6, PA+AHP=-Q is a continuous Lyapunov matrix equation. Meanwhile, we can realize the MFPS of system (7) through controller (16). Therefore, limt→∞e(t)=0, which indicates that the MFPS within system (17) can be achieved. The proof is completed.
Remark 9.
In the proof of Theorem 8, it is easy to see that the scaling function matrix h(t) has no direct effect on dqe/dtq and dqeθ/dtq. The results verify the validity and feasibility of the proposed control strategy.
Remark 10.
Based on Theorem 8, we can choose the scaling function matrix to realize a variety of synchronizations by the method. For example, if h1(t)=h2(t)=⋯=hn(t) are continuously differentiable functions, the MFPS converts to function projective synchronization, and if h(t)=diagh1t,h2t,…,hnt, hi(t)(i=1,2,…,n) being real constants, the MFPS converts to modified projective synchronization (such as complete synchronization and antisynchronization).
Remark 11.
This method has universal applicability if and only if the chaotic system is a PLFO chaotic system (1).
4. Numerical Simulations
Huang and Li [29] proposed a dynamic model of finance comprising three first-order differential equations. Then, following the dynamic model, Chen [24] proposed a fractional-order financial system to describe the running of a financial system as follows:(25)dq1x1dtq1=nx3+mx1x2-ax1,dq2x2dtq2=1-bx2-x12,dq3x3dtq3=px1-cx3.It has three nonlinear equations of fractional-order chaotic systems. The state variables x1, x2, x3 represent the interest rate, the investment demand, and the price index, respectively. Changes in x1 are mostly influenced by two factors: contradictions from the investment market, which indicate a surplus between investment and savings, and structural adjustment from prices of goods. The changing rate of x2 is proportional to the rate of investment but inversely proportional to the cost of investment and interest rates. Changes in x3 are controlled by a contradiction between supply and demand in commercial markets and are affected by inflation rates. Particularly, the negative interest rate x1, negative investment demand x2, and negative price index x3 mean that the banks’ deposit rates are below the level of inflation, the investment market presents overinvestment, and the rate of deflation is increasing, respectively.
The constants a, b, c are representatives of the saving amount, the per-investment cost, and the elasticity of demand of the commercial markets. The fractional-order financial system (25) will exhibit chaotic behavior when qi=q≥0.85,i=1,2,3, and n=1, m=1, p=-1, a=3, b=0.1, and c=1. When q=0.9, the three-dimensional phase diagram of system (25) is illustrated in Figure 1 and the largest Lyapunov exponent is L=0.14. The results highlight the chaotic attractor of fractional-order financial system (25).
Three-dimensional phase diagram of fractional-order financial system (25).
It is obvious that system (25) is a PLFO chaotic system. Therefore, the drive and response systems can be constructed as follows:(26)dq1x1dtq1=nx3+mx1x2-ax1,dq3x3dtq3=px1-cx3,dq2x2dtq2=1-bx2-x12,dq1x4dtq1=n~x5+m~x4x2-a~x4+ψ1,dq3x5dtq3=p~x4-c~x5+ψ2,where ψi(i=1,2) is the controller. In the coupled system (26), the drive system evolves independently, while the response system is governed by the drive system through the variable x2. That is, the drive system is (27)dq1x1dtq1=nx3+mx1x2-ax1,dq3x3dtq3=px1-cx3,dq2x2dtq2=1-bx2-x12,and the response system is (28)dq1x4dtq1=n~x5+m~x4x2-a~x4+ψ1,dq3x5dtq3=p~x4-c~x5+ψ2.System (26) becomes a coupled PLFO chaotic system with uncertain parameters. Based on Theorem 8 and qi(i=1,2,3)=q, we have ur=x1x3, M(x2,θ)=mx2-anp-c, M~(x2,θ~)=m~x2-a~n~p~-c~, h(t)=h1(t)00h2(t), K=k100k2, e=e1e2=x1-h1(t)x4x3-h2(t)x5, em=m~-m, ea=a~-a, en=n~-n, ep=p~-p, and ec=c~-c. According to (16), we can describe the controller as(29)ψ=ψ1ψ2=-htM~x2,θ~ur-Mx2,θe+Ke+dqhturdtq.From (17) to (19), we can obtain the error systems as(30)dqe1dtq=x2e1em+-e1ea+e2en+k1e1,dqe2dtq=e1ep+-e2ec+k2e2.It shows that A11=x2e1, A12=-e1, A13=e2, A24=e1, and A25=-e2. Therefore, the unknown parameters systems are given by(31)dqemdtq=-x2e12+k1′em,dqeadtq=e12+k2′ea,dqendtq=-e1e2+k3′en,dqepdtq=-e1e2+k4′ep,dqecdtq=e22+k5′ec,where K′=diagk1′,k2′,k3′,k4′,k5′(ki′<0).
The initial conditions are x1(0)=-0.5, x2(0)=-0.2, x3(0)=8, x4(0)=3, x5(0)=-1, and q=0.9 and the estimated parameters have initial conditions n~(0)=-1, m~(0)=3, a~(0)=-0.5, p~(0)=3, c~(0)=3.5, and ki=-10(i=1,2) and ki′=-20(i=1,2,3,4,5). Let the scaling function factors h1(t)=0.1cos(2t) and h2(t)=t3; then the simulation results are shown in Figures 2 and 3. Obviously, in Figure 2, the errors converge to zero as time goes to infinity, which implies that the MFPS of system (26) is achieved with the scaling function factors. Figure 3 shows that the estimated values of unknown parameters n~, m~, a~, p~, and c~ converge to n=1, m=1, a=3, p=-1, and c=1, respectively, as t→∞. Therefore, the uncertain parameters n, m, a, p, and c are identified. Additionally, with the scaling function factors h1(t)=0.1cos(2t) and h2(t)=t3, phase portrait of the drive and response systems is illustrated in Figure 4.
MFPS errors of system (30) with h1(t)=0.1cos(2t) and h2(t)=t3.
Unknown parameters of system (31) with h1(t)=0.1cos(2t) and h2(t)=t3.
Phase portrait of the drive and response systems in x1-x3(x4-x5) plane with h1(t)=0.1cos(2t) and h2(t)=t3.
5. Conclusion
In this paper, the MFPS was investigated for the coupled uncertain fractional-order financial system. Based on the stability criterion of the fractional-order system and the continuous Lyapunov matrix equation, a synchronization controller was proposed. Then, using an uncertain financial system, we verified the validity of the proposed controller. The unknown parameters were also estimated so that the phenomenon of parameter distortion was effectively avoided. The results showed that the controller was feasible and effective. In theory, the proposed approach can realize a synchronized and healthy development of financial markets. In the case of possessing homogeneous investment demand, we can keep profits with price index synchronous development so that interference factors can be avoided.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the State Key Program of the National Natural Science of China (no. 71232004) and Fundamental Research Funds for the Central Universities of China (no. CDJSK1001).
ChianA. L.ZorottoF. A.RempelE. L.RogersC.Attractor merging crisis in chaotic business cycles2005243869875ChianA. C.-L.RempelE. L.RogersC.Complex economic dynamics: Chaotic saddle, crisis and intermittency2006295119412182-s2.0-3364653097810.1016/j.chaos.2005.08.218Zbl1142.91652De CesareL.SportelliM.A dynamic IS-LM model with delayed taxation revenues200525123324410.1016/j.chaos.2004.11.044MR2123640Zbl1110.910202-s2.0-13544265300FantiL.ManfrediP.Chaotic business cycles and fiscal policy: an IS-LM model with distributed tax collection lags200732273674410.1016/j.chaos.2005.11.024MR2280116Zbl1133.914822-s2.0-33749578125GrandmontJ.-M.On endogenous competitive business cycles1985535995104510.2307/1911010MR809902Zbl0593.90009DayR. H.1994BostonMIT PressMaC.WangX.Hopf bifurcation and topological horseshoe of a novel finance chaotic system20121727217302-s2.0-8005236167810.1016/j.cnsns.2011.05.029Zbl1241.91008KnightF. H.Risk, uncertainty and profit20054682690WangZ.HuangX.ShenH.Control of an uncertain fractional order economic system via adaptive sliding mode201283838810.1016/j.neucom.2011.11.0182-s2.0-84856604260PecoraL. M.CarrollT. L.Synchronization in chaotic systems199064882182410.1103/PhysRevLett.64.821MR10382632-s2.0-0343689904ChenY.LiX.Function projective synchronization between two identical chaotic systems20071858838882-s2.0-3454774787210.1142/S0129183107010607ChenY.AnH.LiZ.The function cascade synchronization approach with uncertain parameters or not for hyperchaotic systems20081971961102-s2.0-3914913616810.1016/j.amc.2007.07.036Zbl1136.93410YangS. S.DuanC. K.Generalized synchronization in chaotic systems1998910170317072-s2.0-003218241710.1016/S0960-0779(97)00149-5Zbl0946.34040WangL.YuanZ.ChenX.ZhouZ.Lag synchronization of chaotic systems with parameter mismatches201116298799210.1016/j.cnsns.2010.04.029MR2725820Zbl1221.372262-s2.0-77956056801DuH.ZengQ.WangC.Function projective synchronization of different chaotic systems with uncertain parameters200837233540254102-s2.0-4734911293510.1016/j.physleta.2008.06.036MaS.-J.ShenQ.HouJ.Modified projective synchronization of stochastic fractional order chaotic systems with uncertain parameters2013731-2931002-s2.0-8487978426010.1007/s11071-013-0770-6Zbl1281.34095LiZ. G.XuW.Adaptive modified projective synchronization and parameter estimation for chaotic systems with uncertain parameters20102713036MR2798117DuH.ShiP.A new robust adaptive control method for modified function projective synchronization with unknown bounded parametric uncertainties and external disturbances201685135536310.1007/s11071-016-2690-8MR3510618Zbl1349.932102-s2.0-84961154715PanI.DasS.DasS.Multi-objective active control policy design for commensurate and incommensurate fractional order chaotic financial systems201539250051410.1016/j.apm.2014.06.005MR32825912-s2.0-84922608499XieC.XuY.TongD.Chaos synchronization of financial chaotic system with external perturbation2015Art. ID 731376, 710.1155/2015/731376MR33374042-s2.0-84928473675731376DingK.XuX.Mixed synchronization of chaotic financial systems by using linear feedback control2016Art. ID 3269248, 710.1155/2016/3269248MR3579213Zbl067608352-s2.0-850060625043269248WangZ.HuangX.Synchronization of a chaotic fractional order economical system with active controlProceedings of the 2011 International Conference on Advanced in Control Engineering and Information Science, CEIS 2011August 2011chn5165202-s2.0-8405519005210.1016/j.proeng.2011.08.098WangC.ZhangH.-L.FanW.-H.Generalized dislocated lag function projective synchronization of fractional order chaotic systems with fully uncertain parameters20179814212-s2.0-8501491760310.1016/j.chaos.2017.02.010ChenW. C.Nonlinear dynamics and chaos in a fractional-order financial system20083651305131410.1016/j.chaos.2006.07.0512-s2.0-38649117916HanM.XuM.-L.A hybrid prediction model of multivariate chaotic time series based on error correction2013621210.7498/aps.62.1205101205102-s2.0-84880117204XuD.LiZ.Controlled projective synchronization in nonpartially-linear chaotic systems2002126139514022-s2.0-003660012410.1142/S0218127402005170MatignonD.Stability result for fractional differential equations with applications to control processing1997963968ZhouP.ZhuW.Function projective synchronization for fractional-order chaotic systems201112281181610.1016/j.nonrwa.2010.08.008MR2736173Zbl1209.340652-s2.0-78249280987HuangD. S.LiH. Q.1993ChengduHouse of Sichuan University