Control and synchronization of the financial systems with fractional-order are discussed in this paper. Based on the stability theory of fractional-order differential equations, Routh-Hurwitz stability condition, and by using linear control, simpler controllers are designed to achieve control and synchronization of the fractional-order financial systems. The proposed controllers are linear and easy to implement, which have improved the existing results. Theoretical analysis and numerical simulations are shown to demonstrate the validity and feasibility of the proposed method.
1. Introduction
Chaos, as a very interesting nonlinear phenomenon, has been intensively investigated in many fields over the last four decades. Since the pioneering work [1] that chaotic dynamics could exist in an economical model, research on the dynamical behavior of economical model has become one of the most interesting and important topics which have received increasing attention. Many continuous chaotic models have been proposed to study complex economic dynamics, such as the forced Vander-Pol model [2], the IS-ML model [3], Behrens-Feichtinger model [4], and Cournot-Puu model [5]. Just as all the other chaotic systems in engineering, financial chaotic system has complex dynamical behaviors and possess some special features, such as excessive sensitivity to initial conditions, the complex patterns of phase portraits, positive Lyapunov exponents, and bounded and fractal properties of the motion in the phase space. These features are inherent properties of the system itself, rather than caused by external disturbances, which denote some economic behaviors in the fields of finance, stocks, and social economics. In addition, as a nonlinear system, there exist many attractors with different topology, such as fixed points, limit cycles, quasiperiodic attractors, and chaotic attractors in chaotic financial systems, which make the contradiction in economic operation process complex and changeful. Chaotic behavior in financial systems is undesirable due to threatening the safety of investment. Therefore, to study the chaotic behaviors in nonlinear economical systems plays a very important role in decision-making by policy makers. Once the instability of a periodic solution, bifurcation, multiperiodic bifurcation, or similar phenomena appear in the economic chaotic systems, decision makers should take some measures to intervene. Some researches suggest that controlling chaos should improve the performance of a chaotic economy and stabilizing periodic solutions of a chaotic market model may increase economic efficiency [6–10]. In recent years, some authors have paid much attention to chaos control in economic models and obtained some results. For example, Ahmed et al. and Agiza used OGY method to control chaos in economic systems [11, 12], respectively. Pyragas [7] and Holyst and Urbanowicz [9] proposed the delayed feedback control, which has been widely used for controlling chaos in economic models. Sun et al. adopted impulsive control to control a financial model [13]. Du et al. applied phase space compression to control chaos in economic systems [14].
On the other hand, it has been found that many systems in interdisciplinary fields can be described by fractional differential equations, for example, dielectric polarization, electrode-electrolyte polarization, electromagnetic waves, viscoelastic systems, quantitative finance, and diffusion wave [15–17]. Recently, fractional-order systems have become an active research area, particularly in control and synchronization of chaotic systems. Interest has been growing in fractional calculus not only from physicists and engineers but also from researchers in life science and economics [18–20]. In fact, financial variables possessing long memories make fractional models more appropriate for dynamic behaviors in finance. Research on fractional-order financial models has a wider range of applications. In 2008, Chen studied nonlinear dynamics and chaos in a fractional-order financial system [21]. As he pointed out, one of the major differences between fractional-order and integer-order models is that fractional-order models possess memory; that is, the fractional-order model depends on the history of the system. The magnitude of the financial variables such as foreign exchange rates, gross domestic product, interest rates, production, and stock market prices can have very long memory; the reason for describing financial systems using a fractional nonlinear model is that it simultaneously possesses memory and chaos. From then on, researchers set out to investigate the fractional financial models, for instance, in [22] a sliding mode controller was designed for a fractional-order chaotic financial system and in [23] control of a fractional-order chaotic financial system by nonlinear feedback control was discussed. As we know, linear feedback control is especially attractive and has been successfully applied to practical implementation, which was adopted in [24–26] to realize control and synchronization of integer-order chaotic systems. However, there exists substantial difference between fractional-order differential systems and integer-order differential ones. Most of the properties, conclusions, and methods to deal with integer-order systems cannot be simply extended to the case of fractional-order ones. Therefore, results about fractional-order chaotic systems are much less than those of integer-order systems. To the best of our knowledge, there are few results about control and synchronization of fractional-order systems via linear feedback control [27, 28].
In this paper, based on the stability theory of fractional-order differential equations, Routh-Hurwitz stability condition, we investigate control and synchronization of fractional-order financial systems proposed by Chen via linear control. Sufficient conditions are established and easy to verify. Compared with sliding mode control and nonlinear feedback control used for discussing fractional-order financial systems in [22, 23], the linear control is economic and easy to implement, through which control and synchronization of fractional-order financial systems will be obtained only by choosing suitable feedback gains. The main job of this paper lies in two aspects. One is to remove chaotic phenomenon from fractional financial system by controlling, which makes prediction impossible in the financial world. The other is to realize harmonious and sustainable development between drive financial systems and response ones by investigating synchronization. The obtained results have a certain value to the theoretical guidance and application.
The remainder of this paper is organized as follows. In Section 2, preliminary results are presented and fractional-order financial system is described. In Section 3, some sufficient criteria for control of the fractional-order financial system are given. In Section 4, we discussed synchronization of the fractional-order financial system via linear error feedback. In Section 5, numerical simulations are given to illustrate the effectiveness of the main results. Finally, conclusions are drawn in Section 6.
2. Preliminaries and System Description
Fractional calculus is a generalization of integration and differentiation to a noninteger-order integrodifferential operator Dtα defined by Dtα={dαdtα,R(α)>0,1,R(α)=0,∫at(dτ)-α,R(α)<0.
Fractional calculus is being used in various fields gradually, such as biophysics, nonlinear dynamics, informatics, and control engineering [29]. There are some definitions for fractional derivatives. The commonly used definitions are Grunwald-Letnikov (GL), Riemann-Liouville (RL), and Caputo (C) definitions.
The Grunwald-Letnikov (GL) derivative with fractional-order α is given by αGDtαf(t)=limh→0fh(α)(t)=limh→0h-α∑i=0[(t-α)/h](-1)i(iα)f(t-ih),
where [·] means the integer part.
The Riemann-Liouvill (RL) fractional derivatives are defined by αRDtαf(t)=dndtn1Γ(n-α)∫atf(τ)(t-τ)(α-n+1)dτ,n-1<α<n,
where Γ(·) is the gamma function, Γ(τ)=∫0∞tτ-1e-tdt.
The Caputo (C) fractional derivative is defined as follows: αCDtαf(t)=1Γ(n-α)∫atf(n)(τ)(t-τ)(n-α-1)dτ,n-1<α<n.
It should be noted that the advantage of Caputo approach is that the initial conditions for fractional differential equations with Caputo derivatives take on the same form as those for integer-order differential, which have well-understood physical meanings. Therefore, in the rest of this paper, the notation D*α is chosen as the Caputo fractional derivative operator αCD.
The fractional-order chaotic financial system can be described by D*α1x=z+(y-a)x,D*α2y=1-by-x2,D*α3z=-x-cz,
where α1, α2, and α3 are the fractional-order, α1,α2,α3∈(0,1). System (2.5) describes the time variation of three state variables: the interest rate x, the investment demand y, and the price index z; a is the saving amount, b is the cost per investment, c is the elasticity of demand of commercial markets, parameters a, b, and c are positive real constants. When a=3, b=0.1, c=1, fractional-order (α1,α2,α3) is taken as (0.86,0.92,0.95), the largest Lyapunov exponent of system (2.5) is greater than 0 [23]. Figure 1 displays chaotic attractors of system (2.5).
Chaotic attractors of system (2.5) with order (0.86,0.92,0.95).
To obtain our results, the following lemma is presented.
Consider the following fractional-order system: D*α1x=f(x,y,z),D*α2y=g(x,y,z),D*α3z=h(x,y,z).
Lemma 2.1 (see [30]).
System (2.6) is asymptotically stable at the equilibrium points if |arg(λi(A))|>απ/2, α=max(α1,α2,α3), i=(1,2,3), for all eigenvalues λi of the Jacobian matrix J:
J=(∂f∂x∂f∂y∂f∂z∂g∂x∂g∂y∂g∂z∂h∂x∂h∂y∂h∂z).
3. Chaos Control3.1. Analysis of Equilibrium Points
Let z+(y-a)x=0,1-by-x2=0,-x-cz=0.
It is easy to obtain that system (2.5) has three equilibrium points if c-b-abc>0,P1=(0,1b,0),P2=((c-b-abc)c,(1+ac)c,(-1c)⋅(c-b-abc)c),P3=(-(c-b-abc)c,(1+ac)c,(1c)⋅(c-b-abc)c).
The Jacobian matrix of system (2.5) at equilibrium point (x*,y*,z*) is J=[y*-ax*1-2x*-b0-10-c].
The characteristic equation of the Jacobian matrix J is given by λ3+(a+b+c-y*)λ2+(ab+ac+bc-by*-cy*+2(x*)2+1)λ+abc-bcy+b+2c(x*)2=0.
When a=3, b=0.1, c=1, P1=(0,10,0). Substituting the parameters a, b, and c into (3.4), one obtains λ3-5.9λ2-6.6λ-0.6=0.
The roots of (3.5) are λ1=6.8730, λ2=-0.8730, and λ3=-0.1; obviously, λ1 is a positive real number, based on the lemma, arg(λ1)=0, arg(λ2)=arg(λ3)=π. So, the equilibrium point P1 is unstable.
When a=3, b=0.1, c=1, P2=(0.6,4,-0.6), P3=(-0.6,4,0.6). Substituting the parameters a, b, and c into (3.4), the characteristic equation of the Jacobian matrix J at equilibrium points P2,P3 is
λ3+0.1λ2+1.2λ+1.2=0.
The roots of (3.6) are λ1=-0.7256, λ2=0.3128+1.2474i, and λ3=0.3128-1.2474i. According to the lemma, λ1 is a negative real number, arg(λ1)=π, arg|(λ2,3)|=1.3251. α=2|arg(λ2,3)|/π=0.8436. Therefore, we can conclude that if the maximum fractional-order α among α1, α2, α3 is less than 0.8436, the equilibrium points P2, P3 are stable. On the contrary, the equilibrium points P2, P3 are unstable.
According to the above analysis, when the maximum fractional-order α among α1,α2, and α3 is less than 0.8436, there exist two stable equilibrium points; when α1=α2=α3=0.8436, system (2.5) will admit a limit cycle; when fractional-order α1, α2, and α3 are all greater than 0.8436, there are no stable equilibria, with all the equilibrium points being unstable, which implies that there may exist chaos for system (2.5). By calculating the values of Lyapunov exponents of system (2.5), it could be found that system (2.5) exhibits chaotic behaviors if maximum fractional-order α among α1, α2, andα3 is greater than or equal to 0.86 [23].
3.2. Chaos Control
In this subsection, linear state feedback controller is designed to control fractional-order chaotic financial system to its equilibrium.
The controlled fractional-order chaotic financial system is given by D*α1x=z+(y-a)x+k1(x-x*),D*α2y=1-by-x2+k2(y-y*),D*α3z=-x-cz+k3(z-z*),
where k1, k2, and k3 denote feedback gains and (x*,y*,z*) is the desired equilibrium point. Obviously, system (3.7) has one equilibrium point (x*,y*,z*).
The Jacobian matrix of system (3.7) at equilibrium point (x*,y*,z*) is J=[y*-a+k1x*1-2x*k2-b0-10k3-c].
The characteristic equation of the Jacobian matrix (3.8) is λ3+(a+b+c-k1-k2-k3-y*)λ2+(2(x*)2ab+ac+bc+k1k2+k2k3+k1k3-(b+c)k1+(y-a-c)k2-bcy*+(y-a-b)k3+1-by*+2(x*)2)λ-bck1+ck1k2+(cy*-ac-1)k2-(by*+ab+2(x*)2)k3+bk1k3+(a-y)k2k3-k1k2k3+abc+b+c-y*+2c(x*)2=0.
Our goal is to find suitable feedback gains such that all the state trajectories of system (3.7) are controlled to its equilibrium point, that is to say, roots of (3.9) should satisfy the conditions in the lemma.
Theorem 3.1.
When a=3, b=0.1, c=1, P1=(0,10,0); system (3.7) stabilizes to equilibrium point P1, if state feedback gains k1, k2, and k3 satisfy the following conditions:
-10≤k1-k3<-6,k1+k3<-6,k2<0.1.
Proof.
Substituting the parameters a, b, and c into (3.9), one obtains
(10λ+1-10k2)(λ2-(k1+k3+6)λ+k1k3-k1+7k3-6)=0.
It is very easy to obtain the roots of (3.11):
λ1=k2-0.1,λ2,3=12(k1+k3+6±(k1-k3+10)(k1-k3+6)).
Note that λ1 is a negative real number, λ2,3 are a pair of conjugate imaginary roots, and the real parts of imaginary root are negative, that is, arg(λ1)=π, arg(λ2,3)>π/2. Therefore, the trajectory of the controlled fractional-order system (3.7) is asymptotically stable at equilibrium point P1.
In Theorem 3.1, we designed three simple linear feedback controllers to ensure controlled system (3.7) stabilized to P1. In practice, two linear feedback controllers or single linear feedback controller in controlled system (3.7) will do the same thing. Then we have the following corollaries.
Corollary 3.2.
Controlled system (3.7) will approach asymptotically to P1 with one of the following conditions about feedback gains:
-10≤k1<-6,k2<0.1,k3=0,k1+k3<-6,-10≤k1-k3<-6,k2=0.
Corollary 3.3.
Controlled system (3.7) will approach asymptotically to P1 if feedback gains satisfy
-10≤k1<-6,k2=0,k3=0.
Theorem 3.4.
When a=3, b=0.1, c=1, P2=(0.6,4,-0.6); if feedback gains k1, k2, and k3 satisfy one of the following conditions, system (3.7) will approach and stabilizes to equilibrium point P2 asymptotically:
k1=0,k2<-0.9,k3=0,k1<121-62161220,k2=0,k3=0.
Proof.
(1) Substitute the parameters a=3, b=0.1, c=1, and k1=k3=0 into (3.9), one obtains
λ3+(110-k2)λ2+65λ+65=0.
According to the Routh-Hurwitz criterion, real parts of these eigenvalues λ1,2,3 of (3.16) are all negative if
110-k2>0,(110-k2)65>65,
that is, k2<-0.9. That implies that the trajectory of the controlled fractional-order system (3.7) is asymptotically stable at equilibrium point P2.
(2) Substituting the parameters a=3, b=0.1, c=1, and k2=k3=0 into (3.9), one obtains
λ3+(110-k1)λ2+(65-1110k1)λ+65-110k1=0.
By applying the Routh-Hurwitz criterion, if the following conditions for feedback gains are met,
110-k1>0,65-1110k1>0,(110-k1)(65-1110k1)>65-110k1,
then real parts of these eigenvalues λ1,2,3 are all negative. It follows that k1<(1/220)(121-62161). Thus, the trajectory of the controlled fractional-order system (3.7) is asymptotically stable at equilibrium point P2.
Remark 3.5.
Actually, we adopt single linear feedback controller to stabilize P2 of controlled system in Theorem 3.4, namely, stabilizing P2 by adding single linear feedback controller on the first state or the second state, but we cannot do it via adding single linear feedback controller on the third state; the reasons are described as follows.
Substituting the parameters a=3, b=0.1, c=1, and k1=k2=0 into (3.9), one obtains
λ3+(110-k3)λ2+(65+910k3)λ-1110k3+65=0.
By using the Routh-Hurwitz criterion, real parts of these eigenvalues λ1,2,3 are all negative if
110-k3>0,-1110k3+65>0,(110-k3)(65+910k3)>-1110k3+65.
By using simple calculation, such k3 is obviously absent. So the trajectory of the controlled fractional-order system (3.7) is not asymptotically stable at equilibrium point P2.
Theorem 3.6.
When a=3, b=0.1, c=1, P2=(-0.6,4,0.6); if feedback gains k1, k2, and k3 satisfy one of the following conditions, system (3.7) will approach and stabilize to equilibrium point P2 asymptotically:
k1=0,k2<-0.9,k3=0,k1<121-62161220,k2=0,k3=0.
Proof.
The proof is the same as that of Theorem 3.4 and so we omit it here.
4. Chaos Synchronization
In this section, we will investigate synchronization of fractional-order financial system (2.5). Three and two simple linear feedback controllers are designed to achieve synchronization, which simplify the existing synchronization schemes and reduce the synchronization cost. Drive system and response system are described as follows, respectively: D*α1xm=zm+(ym-a)xm,D*α2ym=1-bym-xm2,D*α3zm=-xm-czm,D*α1xs=zs+(ys-a)xs+u1,D*α2ys=1-bys-xs2+u2,D*α3zs=-xs-czs+u3,
where u1, u2, and u3 denote the external control inputs, to be designed later. It follows from systems (4.1) and (4.2) that the following error dynamical system is D*q1e1=(ym-a)e1+xse2+e3-u1,D*q2e2=(xm+xs)e1-be2-u2,D*q3e3=-e1-ce3-u3,
where e1=xm-xs, e2=ym-ys, and e3=zm-zs. Our aim is to find suitable control laws ui(i=1,2,3) for stabilizing the error dynamics system (4.3). To this end, the following theorem is proposed.
Theorem 4.1.
For any given initial conditions, synchronization between systems (4.1) and (4.2) will occur if control schemes are defined as follows:
u1=k1e1,u2=k2e2,u3=k3e3,
where k1, k2, and k3 are feedback gains and satisfy the following conditions:
k1>L-a,k2<-b-9L24(L-a-k1),k3<-c.
Proof.
Combining (4.3) with (4.4), the error system (4.3) is given by
D*q1e1=(ym-a-k1)e1+xse2+e3,D*q2e2=(xm+xs)e1-(b+k2)e2,D*q3e3=-e1-(c+k3)e3.
Error dynamical system (4.6) can be rewritten as the following matrix form:
[D*q1e1D*q2e2D*q3e3]=A[e1e2e3],
where
A=[ym-a-k1xs1xm+xs-(b+k2)0-10-(c+k3)].
Suppose that λ is one of the eigenvalues of matrix A and the corresponding noneigenvector is ε=(ε1,ε2,ε3)T, that is,
Aε=λε.
Taking conjugate transposal on both sides of (4.9), one obtains
(Aε)H=λ¯εH.
Equation (4.9) multiplied left by 1/2εH plus (4.10) multiplied right by 1/2ε, we have
εH(12A+12AH)ε=12(λ+λ̅)εHε.
Because a chaotic system has bounded trajectories, there exists a positive constant L, such that |x|<L, |y|<L. Thus, εH12(A+AH)ε=[ε1H,ε2H,ε3H][ym-a-k1xm+2xs20xm+2xs2-(b+k2)000-(c+k3)][ε1ε2ε3]≤[|ε1H|,|ε2H|,|ε3H|][L-a-k132L032L-(b+k2)000-(c+k3)][|ε1||ε2||ε3|].
From (4.11), we have
12(λ+λ̅)=εH((1/2)A+(1/2)AH)εεHε≤1εHε[|ε1H|,|ε2H|,|ε3H|]P[|ε1|,|ε2|,|ε3|]T,
where
P=[L-a-k132L032L-(b+k2)000-(c+k3)].
It is obvious that real parts of all eigenvalues λ are negative and matrix P should be negative definite, namely, the following inequalities hold:
L-a-k1<0,-(b+k2)(L-a-k1)-94L2<0,(b+k2)(c+k3)(L-a-k1)+94L2(c+k3)<0.
Simplifing the above inequalities, one has
k1>L-a,k2<-b-9L24(L-a-k1),k3<-c.
Therefore, based on stability theorem of fractional-order systems, error system (4.6) is asymptotically stable at the origin, which implies that synchronization between systems (4.1) and (4.2) will be achieved.
Based on the above analysis, it is easy to obtain that two linear feedback controllers could also achieve synchronization between systems (4.1) and (4.2). Then, we have the following corollary.
Corollary 4.2.
For any given initial condition, if control schemes are described as u1=k1e1, u3=k3e3 and feedback gains satisfy
L-a<k1<L-a+94bL2,k3<-c,
then the response system (4.2) can synchronize the drive system (4.1).
Proof.
The proof is similar to that of Theorem 4.1. After some computations, we have
12(λ+λ̅)=εH((1/2)A+1/2AH)εεHε≤1εHε[|ε1H|,|ε2H|,|ε3H|]P[|ε1|,|ε2|,|ε3|]T,
where
P=[L-a-k132L032L-b000-(c+k3)].
Matrix P must be negative definite; if the following inequalities hold:
L-a-k1<0,-b(L-a-k1)-94L2<0,b(c+k3)(L-a-k1)+94L2(c+k3)<0,
then one obtains
L-a<k1<L-a+94bL2,k3<-c.
Therefore, real parts of all eigenvalues λ are negative; according to the stability theorem of fractional-order systems, error system (4.6) is asymptotically stable. This means that the slave system (4.2) can asymptotically synchronize the master system (4.1).
Remark 4.3.
References [22, 23] presented control of the model by using sliding mode control and nonlinear control, respectively, but with no consideration of synchronization. In contrast, control cost of linear control is low.
Remark 4.4.
Response time of control and synchronization could be adjusted with suitable state feedback gains and error feedback gains in allowed limits.
5. Numerical Simulations
In this section, to verify theoretical results obtained in the previous section, the corresponding numerical simulations will be performed and an improved predictor-corrector algorithm is applied (see the appendix). In all simulations, fractional-order (α1, α2, and α3) is chosen as (0.86, 0.92, and 0.95) to ensure the existence of chaos in system (2.5).
Chaos Control
The parameters of system (2.5) are selected as a=3, b=0.1, and c=1. The initial values of controlled system (3.7) are chosen as (x(0),y(0),z(0))=(3,1,4), respectively. Based on Theorem 3.1, we chose state feedback gain (k1,k2,k3) as (-9,-5,-2), under these conditions, roots of (3.11) are λ1=-5.1, λ1=-2.5+0.866i, and λ1=-2.5-0.866i. Controlled system (3.7) is asymptotically stable at P1. Figure 2 shows that the controlled system (3.7) can be stabilized to P1. When state feedback gain (k1,k2,k3)=(0,-3,0) and (-2,0,0), roots (λ1,λ2,λ3) of (3.16) and (3.18) are (-2.8256,-0.1372-0.6371i,-0.1372+0.6371i) and (-2.8256,-0.1372+0.6371i,-0.1372-0.6371i), respectively. According to Theorem 3.4, controlled system (3.7) will stabilize at P2. Figures 3 and 4 display the simulation results, respectively. When the parameters of the controlled system (3.7) are selected as above, based on Theorem 3.6, P3 is stabilized. Figures 5 and 6 display the stabilization of the equilibrium point P3 for state feedback gain (0, −8, 0) and (−2, 0, 0), respectively.
The time response of the states for stabilizing system (3.7) to P1(k1=-9,k2=-5,k3=-2).
The time response of the states for stabilizing system (3.7) to P2(k1=0,k2=-3,k3=0).
The time response of the states for stabilizing system (3.7) to P2(k1=-2,k2=0,k3=0).
The time response of the states for stabilizing system (3.7) to P3(k1=0,k2=-8,k3=0).
The time response of the states for stabilizing system (3.7) to P3(k1=-2,k2=0,k3=0).
Chaos Synchronization
When the parameters of system (2.5) are chosen as a=3, b=0.1, and c=1, select the initial values of the drive and the response systems as (xm(0),ym(0),zm(0))=(4,-3,2) and (xs(0),ys(0),zs(0))=(-3,2,4), respectively. For error feedback gain (k1,k2,k3) = (12,9,-2), simulation result of the synchronization between systems (4.1) and (4.2) is shown in Figure 7. The synchronization error states between systems (4.1) and (4.2) and displayed in Figure 8. When error feedback gains k1=8, k2=0, and k3=-2, it can be seen that the derive system (4.1) and the response system (4.2) achieve the synchronization in Figure 9. Figure 10 displays the error state time response between systems (4.1) and (4.2).
Synchronization states of systems (4.1) and (4.2) (k1=12,k2=9,k3=-2).
Synchronization error states of systems (4.1) and (4.2) (k1=12,k2=9,k3=-2).
Synchronization states of systems (4.1) and (4.2) (k1=8,k2=0,k3=-2).
Synchronization error states of systems (4.1) and (4.2) (k1=8,k2=0,k3=-2).
6. Conclusion
In this paper, based on the stability theory of fractional-order systems and Routh-Hurwitz stability condition, some sufficient conditions for control and synchronization of the fractional-order chaotic financial system by linear feedback control have been derived. Finally, numerical simulations are provided to verify the effectiveness of the results obtained. The results are economical, reliable and efficient. It is noted that the method applied in the paper can also be extended to other fractional-order chaotic systems.
Appendix
An improved predictor-corrector algorithm [31] for fractional-order differential equations is presented in brief. In comparison with the classical one-step Adams-Bashforth-Moulton algorithm, the numerical approximation of the improved algorithm is more accurate and the computational cost is lower.
The following differential equation: dαxdtα=f(t,x),0≤t≤T,x(k)=x0(k),k=0,1,2,…,[α]-1,
is equivalent to the Volterra integral equation [32] x=∑k=0[α]-1x0(k)tkk!+1Γ(α)∫0tf(τ,x)(t-τ)1-αdτ.
Set h=T/N, tn=nh(n=0,1,2,…,N). Then the above equation can be discretized as follows: xh(tn+1)=∑k=0[α]-1x0(k)tkk!+hαΓ(α+2)f(tn+1,xhρ(tn+1))+hαΓ(α+2)∑j=0naj,n+1f(tj,xh(tj)),
where aj,n+1={nα+1-(n-α)(n+1)α+1,j=0,(n-j+2)α+1+(n-j)α+1-2(n-j+1)a+1,1≤j≤n,1,j=n+1,xhρ(tn+1)=∑k=0[α]-1x0(k)tn+1kk!+1Γ(α)∑j=0nbj,n+1f(tj,xh(tj)),
in which, bj,n+1=hα/α((n-j+1)α-(n-j)α).
The error estimate e is Max|x(tj)-xh(tj)|=O(hρ)(j=0,1,…,N), where ρ=Min(2,1+α).
Acknowledgments
The authors would thank the referee and the editor for their valuable comments and suggestions. This work is supported by the National Natural Science Foundation of China (no. 60974090), the Fundamental Research Funds for the Central Universities (no. CDJXS11172237), the Specialized Research Fund for the Doctoral Program of Higher Education of China (no. 20093401120001; no. 102063720090013), the Natural Science Foundation of Anhui Province (no. 11040606M12), and the Natural Science Foundation of Anhui Education Bureau (no. KJ2010A035).
StrotzR. H.McanultyJ. C.NainesJ. B.Goodwin's nonlinear theory of the business cycle: an electro-analog solution195321339041110.2307/1905446ZBL0050.36804ChianA. C.Nonlinear dynamics and chaos in macroeconomics20003360110.1142/S0219024900000723FantiL.ManfrediP.Chaotic business cycles and fiscal policy: an IS-LM model with distributed tax collection lags200732273674410.1016/j.chaos.2005.11.0242280116ZBL1133.91482FeichtingerG.1992Berlin, GermanySpringerPuuT.Chaos in duopoly pricing199116573581KopelM.Improving the performance of an economic system: controlling chaos1997732692892-s2.0-003123188710.1007/s001910050044PyragasK.Continuous control of chaos by self-controlling feedback199217064214282-s2.0-3425031555110.1016/0375-9601(92)90745-8HolystJ. A.HagelT.HaagG.WeidlichW.How to control a chaotic economy?19966131422-s2.0-003052700410.1007/BF01202371HolystJ. A.UrbanowiczK.Chaos control in economical model by time-delayed feedback method20002873-45875982-s2.0-003450976610.1016/S0378-4371(00)00395-2ChenL.ChenG.Controlling chaos in an economic model200737413493582-s2.0-3375116937710.1016/j.physa.2006.07.022AhmedE.El-MisieryA.AgizaH. N.On controlling chaos in an inflation-unemployment dynamical system19991091567157010.1016/S0960-0779(98)00192-11713508ZBL0958.91042AgizaH. N.On the analysis of stability, bifurcation, chaos and chaos control of Kopel map199910111909191610.1016/S0960-0779(98)00210-01706776ZBL0955.37022SunJ. T.QiaoF.WuQ.Impulsive control of a financial model200533542822882-s2.0-1254425602710.1016/j.physleta.2004.12.030ZBL1123.91325DuJ. G.HuangT. W.ShengZ. H.ZhangH. B.A new method to control chaos in an economic system201021762370238010.1016/j.amc.2010.07.0362733679ZBL1200.91195KoellerR. C.Polynomial operators, stieltjes convolution, and fractional calculus in hereditary mechanics1986583-425126484488210.1007/BF01176603ZBL0578.73040KoellerR. C.Applications of fractional calculus to the theory of viscoelasticity198451229429810.1115/1.3167616747787ZBL0544.73052HeavisideO.1971New York, NY, USAChelseaCarteaÁ.del-Castillo-NegreteD.Fractional diffusion models of option prices in markets with jumps200737427497632-s2.0-3375109307810.1016/j.physa.2006.08.071LaskinN.Fractional market dynamics20002873-448249210.1016/S0378-4371(00)00387-31802416ScalasE.GorenfloR.MainardiF.Fractional calculus and continuous-time finance20002841–437638410.1016/S0378-4371(00)00255-71773804ChenW. C.Nonlinear dynamics and chaos in a fractional-order financial system2008365130513142-s2.0-3864911791610.1016/j.chaos.2006.07.051DadrasS.MomeniH. R.Control of a fractional-order economical system via sliding mode201038912243424422-s2.0-7794991326610.1016/j.physa.2010.02.025Abd-ElouahabM. S.EddineHamriN.WangJ. W.Chaos control of a fractional-order financial system201020101827064610.1155/2010/270646ZBL1195.91185ChenY.WuX.GuiZ.Global synchronization criteria for a class of third-order non-autonomous chaotic systems via linear state error feedback control201034124161417010.1016/j.apm.2010.04.0132659663ZBL1201.93045ZhangJ. X.TangW. S.Control and synchronization for a class of new chaotic systems via linear feedback20095846756862-s2.0-7045023736110.1007/s11071-009-9509-9ZBL1183.70075SunM.JiaQ.TianL.A new four-dimensional energy resources system and its linear feedback control20093911011082-s2.0-6064908907510.1016/j.chaos.2007.01.125ZBL1197.93122XinB. G.ChenT.LiuY. Q.Synchronization of chaotic fractional-order windmt systems via linear state error feedback control20101085968510.1155/2010/859685OdibatZ. M.CorsonN.Aziz-AlaouiM. A.BertelleC.Synchronization of chaotic fractional-order systems via linear control2010201819710.1142/S02181274100254292656696ZBL1183.34095HilferR.2000SingaporeWorld ScientificMatignonD.Stability results of fractional differential equations with applications to control processingProceedings of the Institute for Mathematics and Computer Science, IEEE International Conference on Systems1996Lille, France963968DengW.Numerical algorithm for the time fractional Fokker-Planck equation200722721510152210.1016/j.jcp.2007.09.0152442403DiethelmK.FordN. J.Analysis of fractional differential equations2002265222924810.1006/jmaa.2000.71941876137ZBL1014.34003