Loading [MathJax]/jax/output/HTML-CSS/fonts/TeX/fontdata.js
A journal of IEEE and CAA , publishes high-quality papers in English on original theoretical/experimental research and development in all areas of automation
Volume 11 Issue 5
May  2024

IEEE/CAA Journal of Automatica Sinica

  • JCR Impact Factor: 15.3, Top 1 (SCI Q1)
    CiteScore: 23.5, Top 2% (Q1)
    Google Scholar h5-index: 77, TOP 5
Turn off MathJax
Article Contents
S. He, K. Sun, and H. Wang, “Dynamics of the fractional-order Lorenz system based on Adomian decomposition method and its DSP implementation,” IEEE/CAA J. Autom. Sinica, vol. 11, no. 5, pp. 1298–1300, May 2024. doi: 10.1109/JAS.2016.7510133
Citation: S. He, K. Sun, and H. Wang, “Dynamics of the fractional-order Lorenz system based on Adomian decomposition method and its DSP implementation,” IEEE/CAA J. Autom. Sinica, vol. 11, no. 5, pp. 1298–1300, May 2024. doi: 10.1109/JAS.2016.7510133

Dynamics of the Fractional-Order Lorenz System Based on Adomian Decomposition Method and Its DSP Implementation

doi: 10.1109/JAS.2016.7510133
More Information
  • Dear Editor,

    Dynamics and digital circuit implementation of the fractional-order Lorenz system are investigated by employing Adomian decomposition method (ADM). Dynamics of the fractional-order Lorenz system with derivative order and parameter varying is analyzed by means of Lyapunov exponents (LEs), bifurcation diagram, chaos diagram and phase diagram. Results show that the fractional-order Lorenz system has rich dynamical behaviors and it is a potential model for application. It is also found that the minimum order is affected by numerical algorithm and time step size. Finally, the fractional-order system is implemented on digital signal processor (DSP). The phase diagrams generated by the DSP are consistent with that generated by simulation.

    Introduction: The fractional-order Lorenz system with a new set of parameters is firstly analyzed by I. Grigorenko and E. Grigorenko [1], and they reported that the system can generate chaos when the total order is 2.91 by a numerical method they derived. Unfortunately, an error was found among the derived numerical method, thus the result in this letter is not reliable [2]. More recently, Jia et al. [3] analyzed dynamics of this system with order q = 0.7, 0.8 and 0.9 and implemented it in analog circuit by employing frequency domain method (FDM) [4]. However, whether this method accurately reflects chaotic characteristics in fractional-order chaotic system was questioned [5]. Another method for solving fractional-order chaotic systems is the Adams-Bashforth-Moulton algorithm (ABM) [6]. It can be used to analyze dynamics with continuous derivative order [7], and some researches of the fractional-order chaos are based on this algorithm [8]. But the calculation speed of this algorithm is very slow, and it consumes too many computer resources [9]. Meanwhile, ADM [10] is employed to obtain numerical solution of the fractional-order chaotic system for its high precision and fast speed of conver- gence [11]. In addition, based on ADM, LEs of the fractional-order system is calculated [12]. Furthermore, circuit design is essential for application of fractional-order chaotic systems. Although analog circuit implementation is widely reported by researchers, digital circuit realization of the fractional-order chaotic system has better flexibility and repeatability [13]. So, we focus on the dynamics of the fractional-order Lorenz system and its DSP implementation by employing ADM in this letter.

    Numerical solution for the fractional order Lorenz system: The fractional-order Lorenz system is presented by [1]

    {Dqt0x=a(yx)Dqt0y=cxxz+dyDqt0z=xybz (1)

    where a, b, c and d are the system parameters, and q is the derivative order. According to [1] and [3], we investigate dynamics and digital circuit realization of this system by fixing a=40, b=3, c=10, and varying d and q. Because ADM converges very fast [11], we choose the first 6 terms of ADM polynomial for the approximate solution in this letter, and the numerical solution of the fractional-order Lorenz system is denoted by

    {xn+1=6j=0κj1hjq/Γ(jq+1)yn+1=6j=0κj2hjq/Γ(jq+1)zn+1=6j=0κj3hjq/Γ(jq+1) (2)

    where h is the integration step-size, Γ() is the Gamma function, and kji() are defined as

    k01=xn,k02=yn,k03=zn (3)
    {κ11=a(κ02κ01)κ12=cκ01+dκ02κ01κ03κ13=bκ03+κ01κ02 (4)
    {κ21=a(κ12κ11)κ22=cκ11+dκ12κ01κ13κ11κ03κ23=κ11κ02+κ01κ12bκ13 (5)
    {κ31=a(κ02κ01)κ32=cκ21+dκ22κ01κ23κ11κ13Γ(2q+1)Γ2(q+1)κ21κ03κ33=κ01κ22+κ11κ12Γ(2q+1)Γ2(q+1)+κ21κ02bκ23 (6)
    {κ41=a(κ32κ31)κ42=cκ31+dκ32κ01κ33κ31κ03(κ21κ13+κ11κ23)Γ(3q+1)Γ(q+1)Γ(2q+1)κ43=κ01κ32+κ31κ02+bκ33+(κ21κ12+κ11κ22)Γ(3q+1)Γ(q+1)Γ(2q+1) (7)
    {κ51=a(κ42κ41)κ52=cκ41+dκ42κ01κ43(κ31κ13+κ11κ33)Γ(4q+1)Γ(q+1)Γ(3q+1)κ21κ23Γ(4q+1)Γ2(2q+1)κ41κ03κ53=κ01κ42+(κ31κ12+κ11κ32)Γ(4q+1)Γ(q+1)Γ(3q+1)+κ21κ22Γ(4q+1)Γ2(2q+1)+κ41κ02bκ43 (8)
    {κ61=a(κ52κ51)κ62=cκ51+dκ52κ01κ53(κ11κ43+κ41κ13)Γ(5q+1)Γ(q+1)Γ(4q+1)(κ21κ33+κ31κ23)Γ(5q+1)Γ(2q+1)Γ(3q+1)κ51κ03κ63=κ01κ52+κ51κ02bκ53+(κ11κ42+κ41κ12)Γ(5q+1)Γ(q+1)Γ(4q+1)+(κ21κ32+κ31κ22)Γ(5q+1)Γ(2q+1)Γ(3q+1). (9)

    According to (2), the chaotic sequences of the fractional-order Lorenz system are obtained with appropriate initial values. Meanwhile, (2) provides a necessary iterative for DSP implementation of the fractional-order Lorenz system.

    Dynamics with parameter varying: In this section, dynamics of the fractional-order Lorenz system with system parameter d and derivate order q varying are investigated. Parameter fixed dynamical analysis method and chaos diagram are used. Here we set N=20000 and h=0.01. Three cases are investigated.

    1) Fix d=25, and vary derivative order q from 0.75 to 1 with step size of 0.0005. The bifurcation diagram and LEs are shown in Fig. 1. It shows that the system generates chaos for 0.813q<1 excepting some periodic windows. Thus the minimum total order for the fractional-order Lorenz system to generate chaos is 2.439 and the corresponding phase diagram is shown in Fig. 2. In addition, the maximum Lyapunov exponent illustrates a decreasing trend as order q increasing.

    Figure  1.  Dynamics of the fractional-order Lorenz system with d=25 and q varying. (a) LEs; (b) Bifurcation diagram
    Figure  2.  Phase diagram of the fractional-order Lorenz system with d = 25 and q = 0.813

    2) Fix q=0.96 and vary d from 0 to 38 with step size of 0.1. When d decreases from 38, the system presents periodical states until it enters into chaos at d=32.1 by the period-doubling bifurcation as shown in Fig. 3(a). Chaos covers most of the range d[9.8,32.1] with several small periodic widows, like d[14.5,16.3][21.1,21.5]. Finally, the system becomes convergent at d=9.8 by a tangent bifurcation. To observe dynamics better, phase diagrams are presented in Fig. 4. When d =15, 21.5, and 37, the system is periodic, and the system is chaotic when d=20. It shows that the system presents different states with different parameter d.

    Figure  3.  Dynamics of the fractional-order Lorenz system with q = 0.96 and d varying. (a) LEs; (b) Bifurcation diagram
    Figure  4.  Phase diagrams of the fractional-order Lorenz system with q = 0.96 and d varying. (a) d=15; (b) d=20; (c) d=21.5; (d) d=37

    3) Vary q from 0.75 to 1 with step size of 0.0025 and vary d from 0 to 38 with step size of 0.38 simultaneously. The maximum Lyapunov exponent based chaos diagram in q-d parameter plane is illustrated in Fig. 5. In this figure, we only draw the case when the maximum Lyapunov exponent is larger than zero. According to Fig. 5, chaos exists in the range of d[10,32]. A high complexity region is observed within d[25,30] and q[0.8,0.97], which is favorable for practical application. So, the fractional-order Lorenz system is a good model for real application. The chaos diagram provides a parameter selection basis for the fractional-order Lorenz system in practical application. Compared with bifurcation analysis results based on FDM as shown in [3], results based on ADM are more detailed and accuracy. It also shows that we can analyze dynamics of the system with q varying continuously, but it is difficulty for FDM to do so.

    Figure  5.  Maximum Lyapunov exponent based chaos diagram.

    Discussion about the minimum order: Obviously, the minimum order for chaos is different for different system parameter. But it is also different when the numerical solution algorithm or time step size h is different [14]. Thus these two aspects are discussed as follows.

    1) Compared with other approaches, chaotic system has a much lower order if it is solved by ADM algorithm. The equilibrium point of this system is (0, 0, 0) and (±b(c+d),±b(c+d),c+d). When d = 25, the eigenvalues at (0, 0, 0) are λ1=45.6608,λ2=30.6608,λ3=3.0000, and the eigenvalues at (±105,±105,35) are λ1=25.2415, λ2=3.6207+17.8795i and λ3=3.620717.8795i. According to the stability theory, the lowest order q to generate chaos is q = 0.8726. It is not difficult to find out that ABM satisfies this result. However, FDM and ADM do not. According to [3], when q = 0.7, the system has rich dynamics and chaos still exists by applying FDM. According to Fig. 1, the minimum order of the system is q = 0.813 by applying ADM. Actually, the stability theory [15] is proposed to analyze fractional-order linear systems. For fractional-order nonlinear systems, Li et al. [16] proved that the stability theory does not always work when the specified matrix J(X) is time-varying. We think that it is more complex to analyze stability of fractional-order nonlinear system. Besides, although FDM and ADM do not satisfy the stability theory, they are widely used and accepted by researchers. In addition, it shows in [17] that different results of a fractional-order system may be achieved when simulations are performed based on different numerical methods. Since FDM and ADM can obtain chaos at a much lower order, they extend the parameter space of fractional-order chaotic systems.

    2) The effect of time step h should be further investigated. As for ADM, when h = 0.01, the lowest order to generate chaos is q = 0.813. We also find that the lowest order decrease with the decrease of the time step size h. As shown in Fig. 6, when h = 0.001, the lowest order is q = 0.505, and the lowest order is q = 0.402 for h = 0.0001. The system generates chaos with lower order when time step h is smaller, but more memory and computing resources are needed. It is not good for the real application. We think h = 0.01 is a suitable choice for general cases. However, the reason why the lowest order decreases with the decrease of time h needs further study.

    Figure  6.  Bifurcation diagrams of the fractional-order Lorenz system under different h. (a) h = 0.001; (b) h = 0.0001

    According to the discussion above, when a minimum order for chaos generation of a fractional-order chaotic system is presented, the certain set of parameters, numerical algorithm and time step size should also be specified.

    Digital circuit implementation: The DSP implementation method for the fractional-order chaotic systems is described in [13]. In this section, the fractional-order Lorenz system is realized according to (2)−(9). Here, the initial value is x0=[1 2 3]. Setting q = 0.8130, d = 25, the phase diagram is shown in Fig. 7(a). The corresponding MATLAB simulation result is illustrated in Fig. 2. Setting q = 0.96, d = 15, the phase diagram is shown in Fig. 7(b), and its MATLAB simulation counterpart is presented in Fig. 4(a). Setting q = 0.96 and varying d (d = 20 and d = 37), the phase diagrams are shown in Figs. 7(c) and 7(d). It can be seen that they consist with phase diagrams as shown in Figs. 4(b) and 4(d). It shows that the fractional-order Lorenz system is implemented on the DSP platform successfully. It lays a hardware foundation for the applications of the fractional-order Lorenz chaotic system.

    Figure  7.  Phase diagrams of the fractional-order Lorenz system recorded by the oscilloscope (a) q = 0.8130 and d = 25; (b) q = 0.96 and d = 15; (c) q = 0.96 and d = 20; (d) q = 0.96 and d = 37

    Conclusion: In this letter, based on ADM algorithm, we investigated the dynamics of the fractional-order Lorenz system. It shows that the fractional-order Lorenz system has rich dynamical characteristics. The system is more complex with smaller derivate order q as the maximum Lyapunov exponent decreases with the increase of q. The lowest order for chaos generation is different according to different numerical algorithms. The fractional-order Lorenz system has a much lower order for chaos if it is solved by ADM algorithm. Meanwhile, the lowest order for chaos is smaller when the time step size h is smaller. Finally, the system is implemented by employing DSP, and phase diagrams generated by the DSP device are consistent with the simulation results. Our further work will focus on real applications of the fractional-order Lorenz system.

    Acknowledgments: This work was supported by the National Natural Science Foundation of China (62061008, 62071496, 61901530).

  • [1]
    I. Grigorenko and E. Grigorenko, “Chaotic dynamics of the fractional Lorenz system,” Phys. Rev. Lett., vol. 91, no. 3, p. 34101, 2003. doi: 10.1103/PhysRevLett.91.034101
    [2]
    I. Grigorenko and E. Grigorenko, “Erratum: Chaotic dynamics of the fractional Lorenz system [Phys. Rev. Lett. 91, 034101(2003)],” Phys. Rev. Lett, vol. 96, no. 19, p. 199902, 2006. doi: 10.1103/PhysRevLett.96.199902
    [3]
    H. Jia, Z. Chen, and W. Xue, “Analysis and circuit implementation for the fractional-order Lorenz system,” Acta Phys. Sin., vol. 62, no. 14, p. 140503, 2013. doi: 10.7498/aps.62.140503
    [4]
    A. Charef, H. Sun, Y. Tsao, et al., “Fractal system as represented by singularity function,” IEEE Trans. Auto. Cont., vol. 37, no. 9, pp. 1465–1470, 1992. doi: 10.1109/9.159595
    [5]
    M. Wang, G. Sun, and Y. Wei, “Limitations of frequency domain approximation in the calculation of fractional order chaotic systems,” J. Herbin,Initit. Tech., vol. 43, no. 5, pp. 8–12, 2011.
    [6]
    F. Sadiki, K. Rasimi, A. Ibraimi, et al., “Approaches developed to solve fractional differential equations,” J. EJENS-Eur. J. Eng. Nat. Sci., vol. 8, p. 83, 2023.
    [7]
    K. Sun, X. Wang, and J. C. Sprott, “Bifurcations and chaos in fractional-order simplified Lorenz system,” Inter. J. Bifur. Chaos, vol. 20, no. 4, pp. 1209–1219, 2010. doi: 10.1142/S0218127410026411
    [8]
    R. Li and W. Chen, “Lyapunov-based fractional-order controller design to synchronize a class of fractional-order chaotic systems,” Nonlinear Dyn., vol. 76, no. 1, pp. 785–795, 2014. doi: 10.1007/s11071-013-1169-0
    [9]
    S. He, K. Sun, and H. Wang, “Solving of fractional-order chaotic system based on Adomian decomposition algorithm and its complexity analyses,” Acta Phys. Sin., vol. 63, no. 3, p. 30502, 2014. doi: 10.7498/aps.63.030502
    [10]
    G. Adomian, “A review of the decomposition method and some recent results for nonlinear equations,” Math. Comp. Model., vol. 13, no. 7, pp. 17–43, 1990. doi: 10.1016/0895-7177(90)90125-7
    [11]
    D. Cafagna and G. Grassi, “Bifurcation and chaos in the fractionalorder Chen system via a time-domain approach,” Int. J. Bifur. Chaos, vol. 18, no. 7, pp. 1845–1863, 2008. doi: 10.1142/S0218127408021415
    [12]
    R. Caponetto and S. Fazzino, “An application of Adomian decomposition for analysis of fractional-order chaotic systems,” Inter. J. Bifur. Chaos, vol. 23, no. 3, p. 1350050, 2013. doi: 10.1142/S0218127413500508
    [13]
    H. Cao, R. Chu, and Y. Cui, “Complex dynamical characteristics of the fractional-order cellular neural network and its DSP implementation,” Fractal and Fractional, vol. 7, no. 8, p. 633, 2023. doi: 10.3390/fractalfract7080633
    [14]
    D. Peng, K. Sun, S. He, et al., “What is the lowest order of the fractional-order chaotic systems to behave chaotically?” Chaos,Soliton Fract., vol. 119, pp. 163–170, 2019. doi: 10.1016/j.chaos.2018.12.022
    [15]
    D. Matignon, “Stability results for fractional differential equations with applications to control processing,” Comp. Eng. Sys. Appl., pp. 963–968, 1997.
    [16]
    L. Li, H. Peng, Q. Luo, et al., “Problem and analysis of stability decidable theory for a class of fractional order nonlinear system,” Acta Phys. Sin., vol. 62, no. 2, p. 20502, 2013. doi: 10.7498/aps.62.020502
    [17]
    M. S. Tavazoei and M. Haeri, “A proof for non existence of periodic solutions in time invariant fractional order systems,” Automatica, vol. 45, no. 8, pp. 1886–1890, 2009. doi: 10.1016/j.automatica.2009.04.001
  • Related Articles

    [1]Xiang Xu, Guang-Ren Duan.

    High-Order Fully Actuated System Models for Strict-Feedback Systems W

    ith Increasing Dimensions

    [J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(12): 2451-2462. doi: 10.1109/JAS.2024.124599
    [2]Nianyin Zeng, Xinyu Li, Peishu Wu, Han Li, Xin Luo. A Novel Tensor Decomposition-Based Efficient Detector for Low-Altitude Aerial Objects With Knowledge Distillation Scheme[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(2): 487-501. doi: 10.1109/JAS.2023.124029
    [3]Wei Chen, Guo-Ping Liu. Privacy-Preserving Consensus-Based Distributed Economic Dispatch of Smart Grids via State Decomposition[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(5): 1250-1261. doi: 10.1109/JAS.2023.124122
    [4]Haijing Wang, Jinzhu Peng, Fangfang Zhang, Yaonan Wang. High-Order Control Barrier Function-Based Safety Control of Constrained Robotic Systems: An Augmented Dynamics Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(12): 2487-2496. doi: 10.1109/JAS.2024.124524
    [5]Dingxin He, HaoPing Wang, Yang Tian, Yida Guo. A Fractional-Order Ultra-Local Model-Based Adaptive Neural Network Sliding Mode Control of n-DOF Upper-Limb Exoskeleton With Input Deadzone[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(3): 760-781. doi: 10.1109/JAS.2023.123882
    [6]Zhe Gao. Stable Model Order Reduction Method for Fractional-Order Systems Based on Unsymmetric Lanczos Algorithm[J]. IEEE/CAA Journal of Automatica Sinica, 2019, 6(2): 485-492. doi: 10.1109/JAS.2019.1911399
    [7]Amit S. Chopade, Swapnil W. Khubalkar, A.S. Junghare, M.V. Aware, Shantanu Das. Design and Implementation of Digital Fractional Order PID Controller Using Optimal Pole-Zero Approximation Method for Magnetic Levitation System[J]. IEEE/CAA Journal of Automatica Sinica, 2018, 5(5): 977-989. doi: 10.1109/JAS.2016.7510181
    [8]Wenjuan Gu, Yongguang Yu, Wei Hu. Artificial Bee Colony Algorithm-based Parameter Estimation of Fractional-order Chaotic System with Time Delay[J]. IEEE/CAA Journal of Automatica Sinica, 2017, 4(1): 107-113.
    [9]Shiming Chen, Xiaoling Chen, Zhengkai Pei, Xingxing Zhang, Huajing Fang. Distributed Filtering Algorithm Based on Tunable Weights Under Untrustworthy Dynamics[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(2): 225-232.
    [10]Kai Chen, Junguo Lu, Chuang Li. The Ellipsoidal Invariant Set of Fractional Order Systems Subject to Actuator Saturation: The Convex Combination Form[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(3): 311-319.
    [11]Kecai Cao, YangQuan Chen, Daniel Stuart. A Fractional Micro-Macro Model for Crowds of Pedestrians Based on Fractional Mean Field Games[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(3): 261-270.
    [12]Mojtaba Naderi Soorki, Mohammad Saleh Tavazoei. Constrained Swarm Stabilization of Fractional Order Linear Time Invariant Swarm Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(3): 320-331.
    [13]YangQuan Chen, Dingyü Xue, Antonio Visioli. Guest Editorial for Special Issue on Fractional Order Systems and Controls[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(3): 255-256.
    [14]Bingsan Chen, Chunyu Li, Benjamin Wilson, Yijian Huang. Fractional Modeling and Analysis of Coupled MR Damping System[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(3): 288-294.
    [15]Jiacai Huang, YangQuan Chen, Haibin Li, Xinxin Shi. Fractional Order Modeling of Human Operator Behavior with Second Order Controlled Plant and Experiment Research[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(3): 271-280.
    [16]Mingxiang Dai, Ying He, Xinmin Yang. Continuous-time System Identification with Nuclear Norm Minimization and GPMF-based Subspace Method[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(2): 184-191.
    [17]Xiaojuan Chen, Jun Zhang, Tiedong Ma. Parameter Estimation and Topology Identification of Uncertain General Fractional-order Complex Dynamical Networks with Time Delay[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(3): 295-303.
    [18]Cuihong Wang, Huanhuan Li, YangQuan Chen. H Output Feedback Control of Linear Time-invariant Fractional-order Systems over Finite Frequency Range[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(3): 304-310.
    [19]Zeeshan Alam, Liguo Yuan, Qigui Yang. Chaos and Combination Synchronization of a New Fractional-order System with Two Stable Node-foci[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(2): 157-164.
    [20]Fudong Ge, YangQuan Chen, Chunhai Kou. Cyber-physical Systems as General Distributed Parameter Systems: Three Types of Fractional Order Models and Emerging Research Opportunities[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(4): 353-357.

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(7)

    Article Metrics

    Article views (243) PDF downloads(28) Cited by()

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return