Next Article in Journal
Computing and Learning Year-Round Daily Patterns of Hourly Wind Speed and Direction and Their Global Associations with Meteorological Factors
Next Article in Special Issue
On the Exact Solution of Wave Equations on Cantor Sets
Previous Article in Journal
Consistency of Learning Bayesian Network Structures with Continuous Variables: An Information Theoretic Approach
Previous Article in Special Issue
Power-Type Functions of Prediction Error of Sea Level Time Series
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Active Control of a Chaotic Fractional Order Economic System

1
Department of Computer Engineering, Faculty of Engineering, University of Tunceli, 62100 Tunceli, Turkey
2
MI Group, Department of Mathematics, Faculty of Sciences and Techniques of Errachidia, Moulay Ismail University, 52000 Meknes, Morocco
3
Department of Mathematics, Faculty of Science, Firat University, 23100 Elazig, Turkey
*
Author to whom correspondence should be addressed.
Entropy 2015, 17(8), 5771-5783; https://doi.org/10.3390/e17085771
Submission received: 3 June 2015 / Revised: 30 July 2015 / Accepted: 3 August 2015 / Published: 11 August 2015
(This article belongs to the Special Issue Wavelets, Fractals and Information Theory I)

Abstract

:
In this paper, a fractional order economic system is studied. An active control technique is applied to control chaos in this system. The stabilization of equilibria is obtained by both theoretical analysis and the simulation result. The numerical simulations, via the improved Adams–Bashforth algorithm, show the effectiveness of the proposed controller.
MSC Classifications:
34H10; 26A33; 65P20

1. Introduction

Fractional calculus has 300-year history. However, applications of fractional calculus in physics and engineering have just begun. Many systems are known to display fractional order dynamics, such as viscoelastic systems, dielectric polarization and electromagnetic waves [1,2,3,4,5,6]. In recent years, the emergence of effective methods in differentiation and integration of non-integer order equations makes fractional order systems more and more attractive for the systems control community [7,8,9,10].
More recently, there has been a new trend to investigate the control and the dynamic behavior of fractional order chaotic systems. It has been shown that nonlinear chaotic systems may keep their chaotic behavior when their models become fractional [11,12,13].
In this paper, the aim is to control a chaotic fractional-order economic system, using a nonlinear active control method.
This paper is organized as follows: Some preliminaries about fractional calculus, the stability criterion and the numerical algorithm are given in Section 2. The fractional order economic system and its dynamical behavior are described in Section 3. The active control method and the numerical simulations are presented in Section 4. Concluding remarks are drawn in Section 5.

2. Preliminary Tools

2.1. Fractional Calculus

Historical introductions on fractional-order differential equations (FDEs) can be found in [3,4,5,6,14]. Commonly-used definitions for fractional derivatives are due to Riemann–Liouville and Caputo [15]. In what follows, Caputo derivatives are considered, taking the advantage that this allows for traditional initial and boundary conditions to be included in the formulation of the considered problem.
Definition 1. A real function f ( x ) , x > 0 , is said to be in the space C μ , μ R if there exits a real number λ > μ , such that f ( x ) = x λ g ( x ) , where g ( x ) C [ 0 , ) , and it is said to be in the space C μ m if and only if f ( m ) C μ for m I N .
Definition 2. The Riemann–Liouville fractional integral operator of order α of a real function f ( x ) C μ , μ - 1 , is defined as:
J α f ( x ) = 1 Γ ( α ) 0 x ( x - t ) α - 1 f ( t ) d t , α > 0 , x > 0 ,
and J 0 f ( x ) = f ( x ) . The operators J α have some properties, for α , β 0 and ξ - 1 :
  • J α J β f ( x ) = J α + β f ( x ) ,
  • J α J β f ( x ) = J β J α f ( x ) ,
  • J α x ξ = Γ ( ξ + 1 ) Γ ( α + ξ + 1 ) x α + ξ .
Definition 3. The Caputo fractional derivative D α of a function f ( x ) of any real number α, such that m - 1 < α m , m I N , for x > 0 and f C - 1 m in terms of J α , is:
D α f ( x ) = 1 Γ ( m - α ) 0 x ( x - t ) m - α - 1 f ( m ) ( t ) d t
and has the following properties for m - 1 < α m , m I N , μ - 1 and f C μ m :
  • D α J α f ( x ) = f ( x ) ,
  • J α D α f ( x ) = f ( x ) - k = 0 m - 1 f ( k ) ( 0 + ) x k k ! , f o r x > 0 .

2.2. Stability Criterion

In order to investigate the dynamics and to control the chaotic behavior of a fractional order dynamic system:
D t α X ( t ) = F ( X ( t ) ) ,
where X ( t ) = ( x 1 , x 2 , . . . , x n ) T R n , t > t 0 , t [ 0 , T ] , α ( 0 , 1 ) and F : T × R n R n is continuous in X. We will need the following indispensable stability theorem ( see Figure 1).
Theorem 1 (See [16]) For a given commensurate fractional order system (3), the equilibria can be obtained by calculating F ( X ) = 0 . These equilibrium points are locally-asymptotically stable if all of the eigenvalues λ of the Jacobian matrix J = F X at the equilibrium points satisfy:
| a r g ( λ ) | > π 2 α .
Figure 1. Stability region of the fractional order system (3).
Figure 1. Stability region of the fractional order system (3).
Entropy 17 05771 g001

2.3. The Adams–Bashforth–Moulton Algorithm

We recall here the improved version of Adams–Bashforth–Moulton algorithm [17] for the fractional-order systems. Consider the fractional order initial value problem:
D t α x = f ( x ( t ) ) 0 t T , x ( k ) ( 0 ) = x 0 ( k ) , k = 0 , 1 , . . . , m - 1 .
It is equivalent to the Volterra integral equation:
x ( t ) = k = 0 [ α ] - 1 x 0 ( k ) t k k ! + 1 Γ ( α ) 0 t ( t - s ) α - 1 f ( s , x ( s ) ) d s .
Diethelm et al. have given a predictor-corrector scheme (see [17]), based on the Adams–Bashforth–Moulton algorithm, to integrate Equation (6). By applying this scheme to the fractional order system (5), and setting:
h = T N , t n = n h , n = 0 , 1 , . . . , N ,
Equation (6) can be discretized as follows:
x h ( t n + 1 ) = k = 0 [ α ] - 1 x 0 ( k ) t k k ! + h α Γ ( α + 2 ) f ( t n + 1 , x h p ( t n + 1 ) ) + h α Γ ( α + 2 ) j = 0 n a j , n + 1 f ( t j , x h ( t j ) ) ,
where:
a j , n + 1 = n α + 1 - ( n - α ) ( n + 1 ) α , j = 0 , ( n - j + 2 ) α + 1 + ( n - j ) α + 1 - 2 ( n - j + 1 ) α + 1 , 1 j n , 1 , j = n + 1 ,
and the predictor is given by:
x h p ( t n + 1 ) = k = 0 [ α ] - 1 x 0 ( k ) t k k ! + 1 Γ ( α ) j = 0 n b j , n + 1 f ( t j , x h ( t j ) ) ,
where b j , n + 1 = h α α ( ( n + 1 ) - j ) α - ( n - j ) α . The error estimate of the above scheme is:
max j = 0 , 1 , . . , N | x ( t j ) - x h ( t j ) | = O ( h p ) ,
in which p = m i n ( 2 , 1 + α ) .

3. A Fractional Order Economic System

We consider a 3D system of fractional order autonomous differential equations; this system can be interpreted as an idealized macroeconomic model with foreign capital investments [18]. It can be described by:
D α x = m y + p x ( d - y 2 ) , D α y = - x + c z , D α z = s x - r y .
where α ( 0 , 1 ] and the state variables, x, y and z, are the savings of households, the gross domestic product (GDP) and the foreign capital inflow, respectively. Furthermore, the fractional derivation is considered with respect to time. Positive parameters represent corresponding ratios: m is the marginal propensity to saving, p is the ratio of capitalized profit, d is the value of the potential GDP, c is the output/capital ratio, s is the capital inflow/savings ratio and r is the debt refund/output ratio.

3.1. Dynamical Behavior

When m = 0 . 02 , p = 0 . 4 , c = 50 , d = 1 , r = 0 . 1 and s = 10 . The system (10) has three real equilibria E 0 ( 0 , 0 , 0 ) , E 1 ( 0 . 024 , 2 . 4 , 4 . 8 × 10 - 4 ) and E 2 ( - 0 . 024 , - 2 . 4 , - 4 . 8 × 10 - 4 ) .
At the equilibrium point E 0 , the Jacobian matrix of System (10) is given by:
J | E 0 = p d m 0 - 1 0 c s - r 0 ,
The eigenvalues of above matrix are given by:
λ 1 = 1 . 6761839 , λ 2 = - 0 . 6380920 + 2 . 5984529 i a n d λ 3 = - 0 . 6380920 - 2 . 5984529 i .
Hence, the equilibrium point E 0 is unstable. At the equilibrium point E 1 and E 2 , the Jacobian matrix of System (10) is given by:
J | E 1 , 2 = p d - 5 . 76 p m - 0 . 115 p 0 - 1 0 c s - r 0 ,
The eigenvalues of the above matrix are given by:
λ 1 = - 2 . 8852598 , λ 2 = 0 . 4906299 + 2 . 7503585 i a n d λ 3 = 0 . 4906299 - 2 . 7503585 i .
Here, λ 1 is a negative real number and λ 2 and λ 3 are a pair of complex conjugate eigenvalues with positive real parts. Therefore, the equilibrium points E 1 and E 2 are unstable. According to Theorem (4), System (10) exhibits chaotic behavior for α α m i n = 0 . 8876170531 . We stress here that the case α = 1 was studied in [18].

3.2. Numerical Simulations

In order to confirm the chaotic behavior of System (10), numerical simulations were conducted for α = 0 . 9 and the selected initial conditions ( x 0 , y 0 , z 0 ) = ( 0 . 05 , 0 . 1 , 0 . 02 ) . The time histories of the state variables, x, y and z, are graphically presented in Figure 2, and the phase diagrams are shown is Figure 3, while the chaotic attractor is plotted in Figure 4.
Figure 2. The time histories of variables (a) x(t), (b) y(t) and (c) z(t) for α = 0.9.
Figure 2. The time histories of variables (a) x(t), (b) y(t) and (c) z(t) for α = 0.9.
Entropy 17 05771 g002
Figure 3. Phase portraits: (a) x − y, (b) x − z and (c) y − z for System (10) when α = 0.9.
Figure 3. Phase portraits: (a) x − y, (b) x − z and (c) y − z for System (10) when α = 0.9.
Entropy 17 05771 g003
Figure 4. Chaotic attractor xyz for System (10) when α = 0.9.
Figure 4. Chaotic attractor xyz for System (10) when α = 0.9.
Entropy 17 05771 g004

4. Active Control of the Fractional Order Chaotic System

In this section, we investigate the problem of chaos control of the fractional chaotic System (10). In order to control it towards equilibrium points E 0 , E 1 and E 2 , as in [19], we assume that the controlled fractional order autonomous system is given by:
D α x = m y + p x ( d - y 2 ) + U 1 , D α y = - x + c z + U 2 , D α z = s x - r y + U 3 .
where U j ( t ) ( j = 1 , 2 , 3 ) are external active control inputs, which will be suitably-determined later. We prove the following result:
Theorem 2. Starting from any initial condition, an equilibrium point E i of system (11) is asymptotically stable when the controller U j , j = 1 , 2 , 3 , is active, for α α m i n .
Proof. As a Lyapunov candidate function associated with System (11), we consider the quadratic function defined by:
V ( t , ( X ( t ) - E i ) ) = ( x - x i ) 2 + ( y - y i ) 2 + ( z - z i ) 2 ,
where X = ( x , y , z ) T and E i = ( x i , y i , z i ) T is an equilibrium point. Note that V is a positive-definite function on R 3 . From system (11), we have:
D α V ( t , ( X ( t ) - E i ) ) = - 2 V ( t , ( X ( t ) - E i ) ) < 0 .
According to the Lyapunov theory, the equilibrium point E i is asymptotically stable.
To stabilize the chaotic orbits in (10) to its equilibrium E 0 (respectively, E 1 or E 2 ), we need to add the following active controllers:
  • For E 0 :
    U 1 : = - 1 . 4 x + 0 . 98 y + 0 . 4 x y 2 , U 2 : = - y - 49 . 9 z , U 3 : = - 10 x - z .
  • For E 1 :
    U 1 = 0 . 904 x + 1 . 02608 y - 0 . 002304 + 0 . 4 x y 2 + 1 . 92 x y + 0 . 0096 y 2 , U 2 = - y - 49 . 9 z , U 3 = - 10 x - z .
  • For E 2 :
    U 1 = 0 . 904 x + 1 . 02608 y + 0 . 002304 + 0 . 4 x y 2 - 1 . 92 x y - 0 . 0096 y 2 , U 2 = - y - 49 . 9 z , U 3 = - 10 x - z .
Taking into account the above-described controllers, the equilibria E 0 , E 1 and E 2 are stabilized, and then, the chaos is controlled in system (10).
Now, we implement the improved Adams–Bashforth algorithm for numerical simulations (for t [ 0 , 100 ] and t [ 0 , 300 ] ). The unstable point E 0 has been stabilized, as shown in Figure 5, Figure 6 and Figure 7. We remark that the behaviors of x ( t ) , y ( t ) and z ( t ) start as chaotic; then, when the control is activated at t = 40 , the equilibrium point is rapidly stabilized. The equilibria E 1 and E 2 are stabilized in an analogous way (see Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13).
Figure 5. Time histories of System (11) for x signal at the equilibrium E0 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 5. Time histories of System (11) for x signal at the equilibrium E0 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g005
Figure 6. Time histories of system (11) for y signal at the equilibrium E0 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 6. Time histories of system (11) for y signal at the equilibrium E0 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g006
Figure 7. Time histories of System (11) for z signal at the equilibrium E0 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 7. Time histories of System (11) for z signal at the equilibrium E0 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g007
Figure 8. Time histories of System (11) for x signal at the equilibrium E1 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 8. Time histories of System (11) for x signal at the equilibrium E1 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g008
Figure 9. Time histories of System (11) for y signal at the equilibrium E1 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 9. Time histories of System (11) for y signal at the equilibrium E1 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g009
Figure 10. Time histories of System (11) for z signal at the equilibrium E1 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 10. Time histories of System (11) for z signal at the equilibrium E1 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g010
Figure 11. Time histories of System (11) for x signal at the equilibrium E2 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 11. Time histories of System (11) for x signal at the equilibrium E2 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g011
Figure 12. Time histories of System (11) for y signal at the equilibrium E2 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 12. Time histories of System (11) for y signal at the equilibrium E2 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g012
Figure 13. Time histories of System (11) for z signal at the equilibrium E2 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Figure 13. Time histories of System (11) for z signal at the equilibrium E2 with α = 0.9: (a) tmax = 100, (b) tmax = 300.
Entropy 17 05771 g013

5. Conclusion

In this paper, chaos control of a fractional-order chaotic economic system is studied. Furthermore, we have studied the local stability of the equilibria using the Matignon stability condition. Analytical conditions for nonlinear active control have been implemented. Simulation results have illustrated the effectiveness of the proposed control method.

Acknowledgments

The authors would like to thank the reviewers for their valuable comments and suggestions to improve the present work.

Author Contributions

Zakia Hammouch designed the research, Toufik Mekkaoui performed the numerical experiment, Haci Mehmet Baskonus and Hasan Bulut analysed the data. All authors have read and approved the final version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Belgacem, F.B.M.; Baskonus, H.M.; Bulut, H. Variational Iteration Method for Hyperchaotic Nonlinear Fractional Differential Equations Systems. In Advances in Mathematics and Statistical Sciences, Proceedings of the 3rd International Conference on Mathematical, Computational and Statistical Sciences, (MCSS’15), Dubai, United Arab Emirates, 22–24 February 2015; pp. 445–453.
  2. Miller, K.S.; Rosso, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  3. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  4. Nagy, A.M.; Sweilam, N.H. An efficient method for solving fractional Hodgkin–Huxley model. Phys. Lett. A 2014, 378, 1980–1984. [Google Scholar] [CrossRef]
  5. Sweilam, N.H.; Nagy, A.M.; El-Sayed, A.A. Second kind shifted Chebyshev polynomials for solving space fractional order diffusion equation. Chaos Solitons Fractals 2015, 73, 141–147. [Google Scholar] [CrossRef]
  6. Petras, I. Fractional-Order Nonlinear Systems: Modeling, Analysis and Simulation; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  7. Baleanu, D.; Guvenc, B.; Tenreiro-Machado, J.A. New Trends in Nanotechnology and Fractional Calculus Applications; Springer: New York, NY, USA, 2010. [Google Scholar]
  8. Bulut, H.; Belgacem, F.B.M.; Baskonus, H.M. Some New Analytical Solutions for the Nonlinear Time-Fractional KdV-Burgers-Kuramoto Equation. In Advances in Mathematics and Statistical Sciences, Proceedings of the 3rd International Conference on Mathematical, Computational and Statistical Sciences, (MCSS’15), Dubai, United Arab Emirates, 22–24 February 2015; pp. 118–129.
  9. Caponetto, R.; Dongola, G.; Fortuna, L. Fractional Order Systems: Modeling and Control Application; World Scientific: Singapore, Singapore, 2010. [Google Scholar]
  10. Caponetto, R.; Fazzino, S. An application of Adomian Decomposition Method for analysis of fractional-order chaotic systems. Int. J. Bifurc. Chaos 2013, 23, 1–7. [Google Scholar] [CrossRef]
  11. Tavazoei, M.; Haeri, M. A necessary condition for double scroll attractor existence in fractional-order systems. Phys. Lett. A 2007, 367, 102–113. [Google Scholar] [CrossRef]
  12. Wang, S.; Yu, Y. Application of Multistage Homotopy-perturbation Method for the Solutions of theChaotic Fractional Order Systems. Int. J. Nonlinear Sci. 2012, 13, 3–14. [Google Scholar]
  13. Zaslavsky, G. Hamiltonian Chaos and Fractional Dynamics; Oxford University Press: Oxford, UK, 2008. [Google Scholar]
  14. Podlubny, I. Fractional Differential Equations: Mathematics in Science and Engineering; Academic Press: Waltham, MA, USA, 1999. [Google Scholar]
  15. Caputo, M. Linear models of dissipation whose Q is almost frequency independent. J. R. Astral. Soc. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  16. Matignon, D. Stability results for fractional differential equations with applications to control processing. In Proceedings of IMACS Multiconference on Computational Engineering in Systems Applications, Lille, France, 9–12 July 1996; pp. 963–968.
  17. Diethelm, K.; Ford, N.; Freed, A.; Luchko, A. Algorithms for the fractional calculus: A selection of numerical method. Comput. Methods Appl. Mech. Eng. 2005, 94, 743–773. [Google Scholar] [CrossRef]
  18. Volos, C.K.; Ioannis, M.K.; Ioannis, N.S. Synchronization Phenomena in Coupled Nonlinear Systems Applied in Economic Cycles. WSEAS Trans. Syst. 2012, 11, 681–690. [Google Scholar]
  19. Hammouch, Z.; Mekkaoui, T. Control of a new chaotic fractional-order system using Mittag-Leffler stability. Nonlinear Stud. 2015. to appear. [Google Scholar]

Share and Cite

MDPI and ACS Style

Baskonus, H.M.; Mekkaoui, T.; Hammouch, Z.; Bulut, H. Active Control of a Chaotic Fractional Order Economic System. Entropy 2015, 17, 5771-5783. https://doi.org/10.3390/e17085771

AMA Style

Baskonus HM, Mekkaoui T, Hammouch Z, Bulut H. Active Control of a Chaotic Fractional Order Economic System. Entropy. 2015; 17(8):5771-5783. https://doi.org/10.3390/e17085771

Chicago/Turabian Style

Baskonus, Haci Mehmet, Toufik Mekkaoui, Zakia Hammouch, and Hasan Bulut. 2015. "Active Control of a Chaotic Fractional Order Economic System" Entropy 17, no. 8: 5771-5783. https://doi.org/10.3390/e17085771

APA Style

Baskonus, H. M., Mekkaoui, T., Hammouch, Z., & Bulut, H. (2015). Active Control of a Chaotic Fractional Order Economic System. Entropy, 17(8), 5771-5783. https://doi.org/10.3390/e17085771

Article Metrics

Back to TopTop