Advanced
Chaos in PID Controlled Nonlinear Systems
Chaos in PID Controlled Nonlinear Systems
Journal of Electrical Engineering and Technology. 2015. Jul, 10(4): 1843-1850
Copyright © 2015, The Korean Institute of Electrical Engineers
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : November 20, 2014
  • Accepted : February 23, 2015
  • Published : July 01, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Günyaz Ablay
Corresponding Author: Dept. of Electrical and Electronic Engineering, Abdullah Gül University, Turkey. (gunyaz.ablay@agu.edu.tr)

Abstract
Controlling nonlinear systems with linear feedback control methods can lead to chaotic behaviors. Order increase in system dynamics due to integral control and control parameter variations in PID controlled nonlinear systems are studied for possible chaos regions in the closed-loop system dynamics. The Lur’e form of the feedback systems are analyzed with Routh’s stability criterion and describing function analysis for chaos prediction. Several novel chaotic systems are generated from second-order nonlinear systems including the simplest continuous-time chaotic system. Analytical and numerical results are provided to verify the existence of the chaotic dynamics.
Keywords
1. Introduction
Chaos control and chaotification (anti-control of chaos) studies are providing significant effect on many novel applications, including high-performance circuits and devices, chemical reactions, biological systems, liquid mixing, secure communications, and critical decision-making in political, economic and military events [1] . Today, it seems that there are three ways to find threedimensional autonomous chaotic systems: chaotification, searching via computer programming and encountering chaos in simulations-experiments by chance [2] . Chaotification is a process that creates chaos from a non-chaotic system, or improves a chaotic system to present a different type of chaos for various potential applications in biological, medical, electronic, mechanical, optical, and chemical systems. The chaotification process usually involves designing a control mechanism. Various chaos control methods, e.g. OGY method [3] , delayed feedback control [4 , 5] , washout filter [6 , 7] , periodic perturbations [8 - 10] , stochastic perturbations [11] , and nonlinear control methods [12] can be used to alter dynamics of chaotic systems, but in this study chaotification means generating a chaotic structure from a non-chaotic system.
In recent years, some studies have been devoted to chaotification in discrete-time and continuous-time dynamical systems. Most of these studies have been done in the endeavor of achieving chaotification in discrete maps, and have showed that chaotification in discrete-time systems is very promising [13 - 19] . On the other hand, for continuous-time dynamical systems, some feedback control based studies recently have been available in the literature with some successful works, including time-delayed feedback and impulsive control approaches [20 - 35] . The chaotification process is different in discrete-time and continuous-time systems since their dynamic properties are completely different, e.g. a discrete chaotic system can have all positive Lyapunov exponents but a continuous chaotic system needs only a positive Lyapunov exponent to achieve stretch and fold mechanism. It is a common way to use trial-error approach through parameter tuning and numerical simulation in order to get chaotic systems [17 , 25 , 36] . In this work, the common classical control tools, i.e. PID control, Routh’s stability criterion and describing function analysis, will be used for chaotification in a systematic way.
The main goal of the study is to investigate and show that chaotic dynamics can occur in PID controlled continuous-time nonlinear systems. PID controller is the most common control approach used in industry [37] . It is also an effective method for controlling chaos [38 - 40] and chaos synchronization [41 , 42] . On the other hand, variations in the PID control parameters can result in generation of chaos from a non-chaotic nonlinear system. Since the integral term of the PID controller increases order of a system, it is possible to produce chaotic systems from second-order autonomous nonlinear systems. In this way, some new chaotic systems, including the simplest continuous-time chaotic system, are presented in this paper. These new chaotic systems pose high potentials for various applications in certain scientific, engineering and industrial fields, e.g. secure communications and circuit designs.
The paper is organized as follows. Section 2 presents novel chaotic dynamics under PID controlled nonlinear systems, and conclusion of the work is given in Section 3.
2. Chaos in Feedback Systems
Feedback control systems are the most common and successful approaches for solving practical control problems. On the other hand, controller parameters can be adjusted in such a way that unexpected behaviors including bifurcations and chaos can occur. Consider a second-order nonlinear system described by
PPT Slide
Lager Image
where a and b are the system parameters, the function
PPT Slide
Lager Image
represents the nonlinear part of the system and u is a feedback controller. Let the feedback controller be [37]
PPT Slide
Lager Image
where k 1 , k 2 and k 3 are the controller gains. By substituting the controller (2) into the system (1) and introducing new state variables as x = σ ,
PPT Slide
Lager Image
, and
PPT Slide
Lager Image
, we get
PPT Slide
Lager Image
where the new system parameters are kd = k 1 + a , kp = k 2 + b and ki = k 3 , and the function g (⋅) is the nonlinear term of the feedback system. The system (3) can be represented by the combination of the linear term G (·) and the nonlinear term g (⋅) , known as Lur’e form depicted in Fig. 1 . Note that the nonlinear parts of the chaotic systems to be generated in this work will depend on the system output
PPT Slide
Lager Image
, while many other nonlinear terms can also be considered (e.g., system (10)).
PPT Slide
Lager Image
A feedback system in the Lur’e form.
The system (3) can exhibit chaotic behaviors depending on its parameters and initial conditions. Chaotic systems include an infinite number of unstable periodic orbits embedded in an ergodic chaotic attractor [3] . This means that a chaotic motion is a type of limit cycle which is perturbed by the stretching and folding mechanism. By considering the process of periodic-doubling route to chaos, occurrence of a stable limit cycle is an important sign to chaotic behavior. Existence of such a limit cycle in nonlinear systems can be predicted with the describing function analysis that is a powerful tool for analyzing nonlinear systems [43] . The condition for the existence of a limit cycle is given by
PPT Slide
Lager Image
where N (⋅) is the describing function of the nonlinear term. Eq. (4) states that each intersection point of the curve G ( ) and the trajectory of −1/ N ( A ) corresponds to a limit cycle, which is illustrated in Fig. 2 . A specific method for analyzing chaotic systems in Lur’e form is based on describing function analysis proposed by Genesio and Tesi [44] . Genesio-Tesi procedure [45] for determining chaotic behavior of a system in Lur’e form necessitates: (a) Existence of an interacting stable limit cycle and unstable fixed point, and (b) suitable filtering effect along the system. This procedure gives approximate necessary conditions for the existence of a homoclinic orbit , whose perturbation can result in chaos [45 , 46] .
PPT Slide
Lager Image
Nyquist plot for stable limit cycle prediction.
A fast and effective way for finding existence of an interacting stable limit cycle and unstable fixed point can be obtained by using Routh’s stability condition and Nyquist plot of the linear subsystem G (·). In the presence of the nonlinear term g (⋅) , if a stable limit cycle or chaotic behavior exists for the system (3), then
PPT Slide
Lager Image
  • (ii) describing functionN(⋅) of the nonlinear termg(⋅) must intersect with the Nyquist plot of linear termG(·) as depicted inFig. 2.
Besides, periodic solutions exist if ki = kdkp . The suitable filtering effect along the system has a qualitative and heuristic nature, but this can easily be solved with numerical simulations. That is to say, after achieving conditions (i) and (ii), a rough bifurcation diagram of the system (3) can be obtained for a selected range of parameters in order to search for the existence of chaos. In the following subsections, system (3) will be evaluated for various nonlinearities for the existence of chaos with the approach outlined above.
- 2.1 Quadratic nonlinearities
For a quadratic nonlinearity,
PPT Slide
Lager Image
, the system (3) can be rewritten as
PPT Slide
Lager Image
The system (6) has a fixed point at xe = (0, 0, 0), and linearization at the fixed point results in the characteristic equation s 3 + kds 2 + kps + ki = 0. When the system parameters are selected as kd = 2.6 , kp = 0.1 and ki = 1.81 , it can easily be shown that unstable fixed point (condition (i)) and stable limit cycle (condition (ii)) conditions are satisfied for the possibility of chaos. A bifurcation diagram exhibiting a period-doubling route to chaos for the system (6) is shown in Fig. 3 for the parameter ki versus maximum value of x . The numerical simulations are done using MATLAB programs. It is clear that there is an interleaving of chaos and order for 1.8 ≤ ki ≤ 2.1 . Fig. 4(a) shows the attractor of this new chaotic system.
PPT Slide
Lager Image
A bifurcation diagram for the quadratic system (6) for the parameter ki versus maximum value of x.
PPT Slide
Lager Image
Chaotic attractors for quadratic nonlinearities: (a) system (6), (b) system (9) and (c) system (10).
For further assessment of this new chaotic system, eigenvalues of the system are calculated as
PPT Slide
Lager Image
Hence, the fixed point xe is a saddle focus and satisfies the inequality | s 1 | > | σ | > 0 , which implies that there is a homoclinic chaos according to Shil’nikov theorem [47 - 49] . In addition, the chaotic attractor has a dissipative feature since the divergence of flows is described by
PPT Slide
Lager Image
Eq. (8) shows that all system orbits will be confined to a specific limit set of zero volume and asymptotic motion will converge to an attractor [50] . The dimension of the attractor can be estimated through Kaplan-Yorke dimension as DKY = 2.0076 (see Table 2 ).
Novel chaotic dynamics in feedback systems.
PPT Slide
Lager Image
Novel chaotic dynamics in feedback systems.
It is interesting to see that when kp = 0 in (6), a one of the simplest chaotic system (five-term) can be obtained. Namely,
PPT Slide
Lager Image
where kd = 2.38 and ki = 1.6 . The system (9) has a fixed point xe = (0, 0, 0) , and linearization at the fixed point yields three eigenvalues, s 1 = −2.61 and s 2,3 = 0.11 ± 0.77 i . Again, the fixed point xe is a saddle focus and satisfies the inequality | s 1 | > | σ | > 0 . Consequently, the system produces a homoclinic chaos. The chaotic attractor has a dissipative feature because
PPT Slide
Lager Image
System (9) is one of the simplest chaotic systems, since it has been proved that four-term autonomous systems cannot exhibit chaos [51] . While there are some other simple five-term chaotic systems (including one quadratic term [52] and two quadratic terms [53] ), this new chaotic system has a much larger chaos region (i.e., the existence of chaos in certain ranges of parameter variations, e.g. 1.55 ≤ ki ≤ 1.65 ) compared to the chaotic dynamic introduced in [52] .
Fig. 4(b) shows the attractor of the chaotic system (9). A bifurcation diagram of the system (9) for the parameter kp versus maximum value of x exhibiting a period-doubling route to chaos is presented in Fig. 5 . Complex solutions exist for kd = 2.38 and 1.1 ≤ ki ≤ 1.65 . The dimension of the attractor is DKY = 2.0179 , but the basin of attractor is small and initial conditions must be chosen carefully.
PPT Slide
Lager Image
Bifurcation diagram of the simplest chaotic system (9) for ki versus maximum value of the x.
The quadratic nonlinearity can also be considered as a product of two different state variables. When the nonlinearity is described by
PPT Slide
Lager Image
, another simple dissipative chaotic system can be introduced by
PPT Slide
Lager Image
where the system parameters are kd = 2.6 , kp = 0.53 and ki = 3.8 . The bifurcation diagram of the system (10) for ki versus xmax is presented in Fig. 6 . It is seen that chaotic solutions exist for 3.5 < ki < 4.5 . The fixed point of the system is at the origin and carries the features of the saddle focus since the eigenvalues are s 1 = −2.87 and s 2,3 σ ± = 0.13 ± 1.14 i . A homoclinic chaos exists since | s 1 | > | σ | > 0 . Fig. 4(c) shows the attractor of this novel chaotic system.
PPT Slide
Lager Image
Bifurcation diagram of the system (10) for ki versus maximum value of x.
- 2.2 Cubic/quintic nonlinearities
The nonlinear function of the system (3) can be composed of a simple cubic nonlinearity,
PPT Slide
Lager Image
where the system parameters are kd = 0.9 , kp = −9 and ki = 5 . System (11) has a fixed point xe = (0,0,0) , which is a saddle point. The unstable fixed point and stable limit cycle conditions (conditions (i) and (ii)) are satisfied for the possibility of chaos. For 3.6 ≤ ki ≤ 5.1 , a period-doubling route to chaos scheme can be observed in the numerical simulations and the bifurcation diagram given in Fig. 8 . The existence of the chaotic attractor ( Fig. 7a ) can be given by
PPT Slide
Lager Image
and observed in numerical simulations. Since divergence of flows is negative, the system (11) is dissipative and the dimension of the attractor is calculated as DKY = 2.0526.
PPT Slide
Lager Image
Chaotic attractors for various nonlinearities, (a) cubic (11), (b) exponential (14), (c) hyperbolic (15).
PPT Slide
Lager Image
Bifurcation diagram of the cubic system (11) for ki versus maximum value of x.
A more complicated form of the cubic nonlinearity, known as odd quartic function, can also yield a chaotic behavior,
PPT Slide
Lager Image
where the system parameters are kd = 1 , kp = −5 and ki = 4 . For the possibility of chaos, conditions (i) and (ii) are satisfied, and a bifurcation diagram of the system exhibiting a period-doubling route to chaos is displayed in Fig. 9 . It is seen form the bifurcation diagram that there is a wide chaos region (2.5 ≤ ki ≤ 5.1), but the continuity of this region is disrupted with periodic windows. Moreover, there is only a single fixed point xe = (0,0,0) , and it is a saddle focus since its eigenvalues are calculated as s 1 = −3.06 and s 2,3 σ ± = 1.03 ± 0.49 i . This system also produces a homoclinic chaos since | s 1 | > | σ | > 0 . The divergence of flows is
PPT Slide
Lager Image
and thus, the system is dissipative. The dimension of the attractor is DKY = 2.1228.
PPT Slide
Lager Image
Bifurcation diagram of the system (12) for ki versus maximum value of x.
Another dissipative chaotic system example can be obtained with a quintic nonlinearity,
PPT Slide
Lager Image
where the system parameters are kd = 1 , kp = −9 and ki = 10 . The bifurcation diagram of the system (13) is shown in Fig. 10 . Similar to the other chaotic systems, this system also has a saddle fixed point xe = (0,0,0) , and exhibits a homoclinic chaos. The dimension of the attractor is calculated as DKY = 2.0754 and Lyapunov exponents of the system are given in Table 2 . Compared to the other chaotic systems introduced in this section, the bifurcation diagram in Fig. 10 shows that the system (13) has the largest chaos region ( 5.5 ≤ ki ≤ 16).
PPT Slide
Lager Image
Bifurcation diagram of the quintic system (13) for ki versus maximum value of x.
- 2.3 Exponential and hyperbolic nonlinearities
The physical systems can also have exponential, hyperbolic and trigonometric nonlinearities. By considering a saturated exponential function, the reference system (3) becomes
PPT Slide
Lager Image
where the system parameters are kd = 1 , kp = −1 and ki = 7 . A bifurcation diagram of the system exhibiting a period-doubling route to chaos is depicted in Fig. 11 . There is an interleaving chaos and order for 0 ≤ ki ≤ 8.2 , interrupted with some periodic windows. The fixed point of the system (14) is at the origin and behaves as a saddle focus, and the system exhibits a homoclinic chaos. The system has a dissipative characteristics because
PPT Slide
Lager Image
, and the dimension of the attractor is found as DKY = 2.1304 . The attractor is shown in Fig.7b .
PPT Slide
Lager Image
Bifurcation diagram of the exponential system (14) for ki versus maximum value of x.
Now, let the system (3) contain a hyperbolic function,
PPT Slide
Lager Image
where the system parameters are kd = 1 , kp = −5 and ki = 5 . The bifurcation diagram of the system is illustrated in Fig. 12 . Similar to the other chaotic systems, the system (15) has a saddle fixed point xe = (0,0,0) and exhibits a homoclinic chaos. The dimension of the attractor of this dissipative chaotic system is calculated as DKY = 2.1667. The attractor is displayed in Fig. 7c .
PPT Slide
Lager Image
Bifurcation diagram of the hyperbolic system (15) for ki versus maximum value of x.
- 2.4 Characterization of the novel chaotic systems
The describing function and stable limit cycle prediction of some common nonlinear functions are provided in Table 1 . It is seen from Table 1 that it is not always possible to get stable limit cycle prediction for all nonlinearities. In addition, while some nonlinearities can yield a stable limit cycle, they may not produce a chaotic behavior under system (3). It is also interesting to see that the common hard nonlinearities, e.g. sign(·) function, cannot yield a stable limit cycle and chaotic behavior by considering the system (3).
Analysis of common nonlinearities and chaos prediction for the system (3).
PPT Slide
Lager Image
Analysis of common nonlinearities and chaos prediction for the system (3).
The characterization of the attractors of chaotic systems is an important aspect. The geometry of chaotic attractors can be characterized with the dimension of an attractor. One significant quantitative characterization method is the Kaplan-Yorke dimension (KYD) [54] . It represents an upper bound for the information dimension of the system. The KYD is calculated by DKY = d + ( λ 1 + ... + λd /| λ d+1 | where λ 1 ≤ λ d are the Lyapunov exponents (LEs) and d is the largest integer for which λ 1 + ... + λd ≥ 0. For the calculation of LEs, the Wolf algorithm [55] is used. It is seen from Table 2 that LEs of the introduced chaotic systems consist of positive, zero and negative exponents. Since the chaotic systems are dissipative, the sum of the LEs is negative. LEs can also be used for unpredictability measure of the chaotic dynamics, i.e. Kolmogolov-Sinai entropy (KSE) [56] calculations. KSE is defined by the sum of the positive LEs, and gives information about the level of the randomness of a chaotic system. The greater KSE means the higher randomness. Table 2 summarizes the LEs, KSE and KYD of the novel continuous-time chaotic systems introduced in this work.
3. Conclusion
This paper presented novel chaotic systems generated from PID controlled second-order nonlinear systems. The study set out to determine the possibility of chaotic behaviors in simple nonlinear systems controlled with PID systems. This research has shown that the integral control term of the PID controller increases system order and makes chaotification possible in second-order systems with simple nonlinearities. The chaos prediction approach is based on the Routh’s stability criterion and describing function analysis. Several new chaotic systems, including the simplest chaotic dynamic, are introduced, and their existences are shown with bifurcation diagrams and Lyapunov exponents. These novel chaotic systems have high potentials in many applications including secure communications, circuit designs and random number generations. This research also shows that bifurcations and chaos can occur in industrial systems regulated by the PID controllers.
Acknowledgements
This work was supported by Research Fund of the Abdullah Gül University. Project No: 4.
BIO
Günyaz Ablay He received his PhD degree in nuclear engineering from the Ohio State University, and his MS and BS degrees in electrical engineering from Firat University. His research interests include control theory and its applications, chaos theory and nuclear energy technologies.
References
Chen G. 2001 Wiley Encyclopedia of Electrical and Electronics Engineering John Wiley & Sons, Inc. “Chaos, Bifurcations, and Their Control,”
Lü J. , Chen G. , Cheng D. 2004 “A new chaotic system and beyond: The generalized Lorenz-like system,” Int. J. Bifurcation Chaos 14 (05) 1507 - 1537    DOI : 10.1142/S021812740401014X
Ott E. , Grebogi C. , Yorke J. A. 1990 “Controlling chaos,” Physical Review Letters 64 (11) 1196 - 1199    DOI : 10.1103/PhysRevLett.64.1196
Pyragas K. 2006 “Delayed feedback control of chaos,” Philos Trans Roy Soc A 364 2309 - 2334    DOI : 10.1098/rsta.2006.1827
Pyragas K. 1992 “Continuous control of chaos by self-controlling feedback,” Physics Letters A 170 (6) 421 - 428    DOI : 10.1016/0375-9601(92)90745-8
Tesi A. , Abed E. H. , Genesio R. , Wang H. O. 1996 “Harmonic balance analysis of period-doubling bifurcations with implications for control of nonlinear dynamics,” Automatica 32 (9) 1255 - 1271    DOI : 10.1016/0005-1098(96)00065-9
Zhou S. , Lin X. , Li H. 2011 “Chaotic synchronization of a fractional-order system based on washout filter control,” Comm. in Nonlinear Science and Numerical Simulation 16 (3) 1533 - 1540    DOI : 10.1016/j.cnsns.2010.06.022
Braiman Y. , Goldhirsch I. 1991 “Taming chaotic dynamics with weak periodic perturbations,” Phys. Rev. Lett. 66 (20) 2545 - 2548    DOI : 10.1103/PhysRevLett.66.2545
Lei Y. , Xu W. , Xu Y. , Fang T. 2004 “Chaos control by harmonic excitation with proper random phase,” Chaos, Solitons & Fractals 21 (5) 1175 - 1181    DOI : 10.1016/j.chaos.2003.12.086
Li R. , Xu W. , Li S. 2006 “Chaos controlling of extended nonlinear Liénard system based on the Melnikov theory,” Applied Mathematics and Computation 178 (2) 405 - 414    DOI : 10.1016/j.amc.2005.11.057
Fahy S. , Hamann D. R. 1992 “Transition from chaotic to nonchaotic behavior in randomly driven systems,” Phys. Rev. Lett. 69 (5) 761 - 764    DOI : 10.1103/PhysRevLett.69.761
Ablay G. 2009 “Sliding mode control of uncertain unified chaotic systems,” Nonlinear Analysis: Hybrid Systems 3 (4) 531 - 535    DOI : 10.1016/j.nahs.2009.04.002
Wang X. F. , Chen G. 1999 “On feedback anticontrol of discrete chaos,” International Journal of Bifurcation and Chaos 09 (07) 1435 - 1441    DOI : 10.1142/S0218127499000985
Wiggins S. 2003 Introduction to applied nonlinear dynamical systems and chaos Springer New York
Chen G. , Shi Y. 2006 “Introduction to anti-control of discrete chaos: theory and applications,” Philos Trans A Math Phys Eng Sci 364 (1846) 2433 - 2447    DOI : 10.1098/rsta.2006.1833
Shi Y. , Chen G. 2005 “Chaotification of discrete dynamical systems governed by continuous maps,” International Journal of Bifurcation and Chaos 15 (02) 547 - 555    DOI : 10.1142/S0218127405012351
Hu G. 2009 “Generating hyperchaotic attractors with three positive lyapunov exponents via state feedback control,” International Journal of Bifurcation and Chaos 19 (02) 651 - 660    DOI : 10.1142/S0218127409023275
Li Z. , Park J.-B. , Joo Y.-H. 2001 “Chaotifying continuous-time TS fuzzy systems via discretization,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications 48 (10) 1237 - 1243    DOI : 10.1109/81.956020
Codreanu S. 2002 “Desynchronization and chaotification of nonlinear dynamical systems,” Chaos, Solitons & Fractals 13 (4) 839 - 843    DOI : 10.1016/S0960-0779(01)00060-1
Zhang J. , Xu D. , Zhou J. , Li Y. 2012 “Chaotification of vibration isolation floating raft system via non-linear time-delay feedback control,” Chaos, Solitons & Fractals 45 (9-10) 1255 - 1265    DOI : 10.1016/j.chaos.2012.05.012
Kuznetsov A. S. , Kuznetsov S. P. 2013 “Parametric generation of robust chaos with time-delayed feed-back and modulated pump source,” Comm. in Non-linear Science and Numerical Simulation 18 (3) 728 - 734    DOI : 10.1016/j.cnsns.2012.08.006
Yu S. , Chen G. 2012 “Chaotifying continuous-time nonlinear autonomous systems,” Int. Journal of Bifurcation and Chaos 22 (09) 12 -
Şahin S. , Guzelis C. 2010 “Chaotification of Real Systems by Dynamic State Feedback,” IEEE Antennas and Propagation Magazine 52 (6) 222 - 233
Xie Q. , Han Z. , Zhang W. 2012 “Chaotification via system immersion,” Journal of Computational and Applied Mathematics 236 (7) 1775 - 1782    DOI : 10.1016/j.cam.2011.10.008
X. F. Wang , G. Chen , X. Yu 2003 Chaos Control Springer Berlin Heidelberg “Generating Chaos in Continuous-Time Systems via Feedback Control,” 179 - 204
Chen H.-K. , Lee C.-I. 2004 “Anti-control of chaos in rigid body motion,” Chaos, Solitons & Fractals 21 (4) 957 - 965    DOI : 10.1016/j.chaos.2003.12.034
Ge Z.-M. , Chang C.-M. , Chen Y.-S. 2006 “Anti-control of chaos of single time scale brushless dc motors and chaos synchronization of different order systems,” Chaos, Solitons & Fractals 27 (5) 1298 - 1315    DOI : 10.1016/j.chaos.2005.04.095
Akhmet M. U. , Fen M. O. 2013 “Replication of chaos,” Comm. in Nonlinear Science and Numerical Simulation 18 (10) 2626 - 2666    DOI : 10.1016/j.cnsns.2013.01.021
Raiesdana S. , Goplayegani S. M. H. 2013 “Study on chaos anti-control for hippocampal models of epilepsy,” Neurocomputing 111 54 - 69    DOI : 10.1016/j.neucom.2012.11.033
Saha P. , Banerjee S. , Chowdhury A. R. 2003 “On the study of control and anti-control in magneto-convection,” Physics Letters A 306 (4) 211 - 226    DOI : 10.1016/S0375-9601(02)01514-1
Moreno-Valenzuela J. 2013 “Adaptive anti control of chaos for robot manipulators with experimental evaluations,” Comm. in Nonlinear Science and Numerical Simulation 18 (1) 1 - 11    DOI : 10.1016/j.cnsns.2012.06.003
Li Y. , Tang W. K. S. , Chen G. 2005 “Hyperchaos evolved from the generalized Lorenz equation,” International Journal of Circuit Theory and Applications 33 (4) 235 - 251    DOI : 10.1002/cta.318
Zhou J. , Xu D. , Zhang J. , Liu C. 2012 “Spectrum optimization-based chaotification using time-delay feedback control,” Chaos, Solitons & Fractals 45 (6) 815 - 824    DOI : 10.1016/j.chaos.2012.02.015
Fossi D. O. T. , Woafo P. 2013 “Generation of complex phenomena in a simple electromechanical system using the feedback control,” Communications in Nonlinear Science and Numerical Simulation 18 (1) 209 - 218    DOI : 10.1016/j.cnsns.2012.06.021
Starkov K. , Chen G. 2004 “Chaotification of poly-nomial continuous-time systems and rational normal forms,” Chaos, Solitons & Fractals 22 (4) 849 - 856    DOI : 10.1016/j.chaos.2004.02.052
Yu S. , Chen G. 2012 “Anti-control of continuous-time dynamical systems,” Communications in Nonlinear Science and Numerical Simulation 17 (6) 2617 - 2627    DOI : 10.1016/j.cnsns.2011.10.001
Åström K. J. , Hägglund T. 2006 Advanced PID Control ISA
Das S. , Acharya A. , Pan I. 2014 “Simulation studies on the design of optimum PID controllers to suppress chaotic oscillations in a family of Lorenz-like multi-wing attractors,” Mathematics and Computers in Simulation 100 72 - 87    DOI : 10.1016/j.matcom.2014.03.002
Alfi A. 2012 “Chaos suppression on a class of uncertain nonlinear chaotic systems using an optimal H∞ adaptive PID controller,” Chaos, Solitons & Fractals 45 (3) 351 - 357    DOI : 10.1016/j.chaos.2012.01.001
Davendra D. , Zelinka I. , Senkerik R. 2010 “Chaos driven evolutionary algorithms for the task of PID control,” Computers & Mathematics with Applications 60 (4) 1088 - 1104    DOI : 10.1016/j.camwa.2010.03.066
Coelho L. dos S. , Grebogi R. B. 2010 “Chaotic synchronization using PID control combined with population based incremental learning algorithm,” Expert Systems with Applications 37 (7) 5347 - 5352    DOI : 10.1016/j.eswa.2010.01.022
Hsu C.-F. , Tsai J.-Z. , Chiu C.-J. 2012 “Chaos synchronization of nonlinear gyros using self-learning PID control approach,” Applied Soft Computing 12 (1) 430 - 439    DOI : 10.1016/j.asoc.2011.08.027
Atherton D. 2011 An Introduction to Nonlinearity in Control Systems Bookboon
Genesio R. , Tesi A. 1991 “Chaos prediction in nonlinear feedback systems,” Control Theory and Applications, IEE Proceedings D 138 (4) 313 - 320    DOI : 10.1049/ip-d.1991.0042
Genesio R. , Tesi A. 1992 “Harmonic balance methods for the analysis of chaotic dynamics in nonlinear systems,” Automatica 28 (3) 531 - 548    DOI : 10.1016/0005-1098(92)90177-H
Guckenheimer J. , Holmes P. 1983 Nonlinear oscillations, dynamical systems, and bifurcation of vector fields Springer-Verlag New York
Turaev D. V. , Shilnikov L. P. 1998 “An example of a wild strange attractor,” Sbornik Mathematics 189 (2) 291 - 314    DOI : 10.1070/SM1998v189n02ABEH000300
Shilnikov L.P. 1965 “A case of the existence of a denumerable set of periodic motions,” Sov. Math. Dokl. 6 163 - 166
Silva C. P. 1993 “Shil’nikov’s theorem-a tutorial,” IEEE Trans Circuits and Systems I: Fundamental Theory and Applications 40 (10) 675 - 682    DOI : 10.1109/81.246142
Hunt B. R. 2004 The Theory of Chaotic Attractors Springer
Fu Z. , Heidel J. 1999 “Non-chaotic behaviour in three-dimensional quadratic systems,” Nonlinearity 12 (3) 739 -    DOI : 10.1088/0951-7715/12/3/501
Sprott J. C. 2000 “Simple chaotic systems and circuits,” American Journal of Physics 68 (8) 758 - 763    DOI : 10.1119/1.19538
Munmuangsaen B. , Srisuchinwong B. 2009 “A new five-term simple chaotic attractor,” Physics Letters A 373 (44) 4038 - 4043    DOI : 10.1016/j.physleta.2009.08.068
Frederickson P. , Kaplan J. L. , Yorke E. D. , Yorke J. A. 1983 “The liapunov dimension of strange attractors,” Journal of Differential Equations 49 (2) 185 - 207    DOI : 10.1016/0022-0396(83)90011-6
Wolf A. , Swift J. B. , Swinney H. L. , Vastano J. A. 1985 “Determining Lyapunov exponents from a time series,” Physica D: Nonlinear Phenomena 16 (3) 285 - 317    DOI : 10.1016/0167-2789(85)90011-9
Frigg R. 2004 “In what sense is the Kolmogorov-Sinai Entropy a measure for chaotic behaviour? Bridging the gap between dynamical systems theory and communication theory,” British J. for the Philosophy of Science 55 (3) 411 - 434    DOI : 10.1093/bjps/55.3.411