Control of the Lorenz Equations. Richard B. Choroszucha
|
|
- Jerome Cox
- 6 years ago
- Views:
Transcription
1 Control of the Lorenz Equations Richard B. Choroszucha 26.IV.2010
2 Contents List of Figures iii List of Tables iv 1 Lorenz Equations Spotlighting Chaos Analysis of Lorenz Equations Equilibrium Points Linearisation Lyapunov Argument Controllability and Observability Control Techniques Input State Linearisation Sliding Mode Controller Time-Delayed Feedback Control Conclusions Concluding Remarks Areas for Further Research References 18 i
3 CONTENTS Choroszucha ii A Derivation 20 B Code 26 B.1 Open Loop Simulations B.2 Input State Linearisation B.3 Sliding Mode Control B.4 Time-Delayed Feedback Control B.4.1 Time-Step Algorithm B.4.2 Delayed Differential Equation
4 List of Figures 1.1 Relative Tolerance , x 1 (0) = [ ] T Relative Tolerance , x 1 (0) = [ ] T Relative Tolerance , x 2 (0) = [ ] T Relative Tolerance , x 2 (0) = [ ] T States using Input State Solution Trajectories Control States using SMC Solution Trajectories Control States using TDFC Solution Trajectories Control iii
5 List of Tables 1.1 Declaration of Variables in Lorenz Equation Variable Declarations for TDFC iv
6 Chapter 1 Lorenz Equations The Lorenz equations are a simplified model of convective incompressible air flow between two horizontal plates with a temperature difference, subject to gravity. The motivation for these equations were to spotlight why weather is unpredictable, despite being a deterministic system. where The non-dimensionalised Lorenz Equations are: ẋ = σ(y x) (1.1) ẏ = (r z)x y (1.2) ż = xy bz, (1.3) x : Fluid Motion y : Horizontal Temperature z : Vertical Temperature σ : Prandtl Number [22] r : Rayleigh Number [22] b : Geometric Factor [23] τ : Non-dimensionalised time. Table 1.1: Declaration of Variables in Lorenz Equation The derivation of (1.1)-(1.3) are very involved and provide nothing in the way of insight due to approximations made for a simplified model that exhibits chaos. The derivation has been left for the readers review in Appendix A. For this paper, σ = 10, r = 28 and b = 8/3, perhaps the most common choice. 1
7 CHAPTER 1. LORENZ EQUATIONS Choroszucha Spotlighting Chaos Chaos is the strong dependence of the solution trajectories on the initial conditions, or the exponential divergence in solution trajectories for small differences in the initial conditions. For great reviews of chaos and attractors, please see [22, ] and [11]. To highlight chaos in the Lorenz Equations the author has provided the following four figures of the open loop (no control) Lorenz Equations to demonstrate how tiny differences in the initial conditions and tolerances on the integrator can propagate through the system causing drastically different solution trajectories. The open loop simulations were performed using MATLAB s ode45. Figures 1.1 and 1.2, pages 3 and 4, respectively, have the same initial conditions, but the Relative Tolerance on the error for the integrator has changed. To facilitate comparison, the author has placed each figure on a separate page at the top so one may flip back and forth quickly, much like a blink comparator. As can be seen by Figures 1.1 and 1.2, the solution trajectories are similar until about τ = 17.83, then after that they no longer share any resemblance. To call out the dependence on the initial conditions, Figure 1.3 and 1.4, pages 5 and 6, respectively, are simulations using the same tolerances, but with a small deviation from the original initial condition x 1 (0). As was the case for x 1 (0), x 2 (0) solutions remains close to the x 1 (0) solutions until τ 20, then are completely incomparable for both sets of tolerances. For code please see Appendix B.1. For the remainder of the paper, all calculations will be carried out with x(0) = [ ] T and RelTol= with all other settings being left as the default. This begs the questions, which has more of an effect, the error or initial conditions, and how can we control such a maniacal mathematical model?
8 CHAPTER 1. LORENZ EQUATIONS Choroszucha 3 Figure 1.1: Relative Tolerance , x 1 (0) = [ ] T
9 CHAPTER 1. LORENZ EQUATIONS Choroszucha 4 Figure 1.2: Relative Tolerance , x 1 (0) = [ ] T
10 CHAPTER 1. LORENZ EQUATIONS Choroszucha 5 Figure 1.3: Relative Tolerance , x 2 (0) = [ ] T
11 CHAPTER 1. LORENZ EQUATIONS Choroszucha 6 Figure 1.4: Relative Tolerance , x 2 (0) = [ ] T
12 Chapter 2 Analysis of Lorenz Equations 2.1 Equilibrium Points Setting the left-hand side of (1.1)-(1.3) to zero and solving, the equilibrium points are xe = x e y e z e = ± b(r 1) ± b(r 1) (r 1). 2.2 Linearisation Linearising the system through the Jacobian matrix, one obtains σ σ 0 (D x f)( x e ) = r z 1 x (2.1) y x b x e This results in the characteristic equation (λi (D x f)( x e )) = λ 3 + (b + σ + 1)λ 2 + (b + x 2 e bσ rσ + z e σ)λ +(bσ brσ + x 2 eσ + x e y e σ + bz e σ). (2.2) One can solve for the roots of the characteristic polynomial evaluated at the equilibria to obtain the local stability. 7
13 CHAPTER 2. ANALYSIS OF LORENZ EQUATIONS Choroszucha Lyapunov Argument Choosing the standard Lyapunov function candidate, making the substitutions x = x x e, y = y y e and z = z z e V = 1 2 (a(x x e) 2 + c(y y e ) 2 + d(z z e ) 2 ) (2.3) V = a(x x e )ẋ + c(y y e )ẏ + d(z z e )ż (2.4) = a(x x e )(σ((y y e ) (x x e ))) + c(y y e )((r (z z e ))(x x e ) (y y e )) +d(z z e )((x x e )(y y e ) b(z z e )). This leads to a natural definition of the coefficients to be a = 1, c = 1 and d = 1. One can σ then complete the square to arrive at ( V = (x x e ) r + 1 ) ( 2 ( ) ) 2 r (y y e) 1 (y y e ) 2 b(z z e ) 2. (2.5) 2 Implying that there are no unstable orbits in the Lorenz system. For the logic behind this statement, a Lyapunov function candidate was tested, such that V ( x) > 0 unless x = x e V (x e ) = 0, V (x) < 0 unless x = xe in which case V ( x e ) = 0. Using LaSalle s Invariance Theorem from [14], one can then find that this implies all trajectories are stable in the sense of Lyapunov. One may also show stronger stability, such as global stability of the origin, [22, 315]. 2.4 Controllability and Observability Augmenting the Lorenz Equations with controls on all variables ẋ = σ(y x) + α 1 u x (2.6) ẏ = (r z)x y + α 2 u y (2.7) ż = xy bz + α 3 u z (2.8) Y = c 1 x + c 2 y + c 3 z (2.9) one may check the controllability and observability of the system using [10]. Y is the output or measurements. Assuming this scheme, where α i, c i 0, the system is completely controllable and observable. However, in future works this may not be a good assumption to make because it may be impractical to control and measure some states.
14 Chapter 3 Control Techniques 3.1 Input State Linearisation Following the method outlined in [15], because a general case is quite complicated, and letting 0 g(x) = x, (3.1) 0 to apply Input State Linearisation we must satisfy that {g, ad g Y } is involutive and ad 2 fg 0 Y y σ Y x + x Y z = 0 (3.2) = 0 (3.3) These can be solved by inspection, (3.2) implies that the output cannot depend on y and (3.3) implies a solution for Y to be Y = 1 σ x2 2z + c 1, (3.4) where c 1 is a constant of integration and according to [15], this constant is determined by the control goal. This implies that z 1 = 1 σ x2 2z + c 1 (3.5) z 2 = 2x 2 + 2bz (3.6) z 3 = 4σx(y x) + 2b(xy bz). (3.7) Selecting the synthetic input v = β 1 z 1 + β 2 z 2 + β 3 z 3, (3.8) 9
15 CHAPTER 3. CONTROL TECHNIQUES Choroszucha 10 and substituting into the control formula u = 1 L g L 2 f Y (v L3 fy ). (3.9) This is where the treatment of this paper and [15] differs, the author has implemented a tracking scheme for x through modifying c 1, namely x = 5 sin(2πt) , is the equilibrium point in the x direction. For the synthetic control law, the author has chosen to place the linearised system s eigenvalues at -1, meaning the gains are β 1 = 1, β 2 = β 3 = 3. The control is turned on at τ = 20. The results of this simulation are close, but not exact, the frequency is higher than 1 Hz, and the amplitude is just shy of being 5. This can be seen in the following figures. Figure 3.1: States using Input State It should be noted in Figure 3.2 that it has not settled to a circle, it is closer to something like x 4 + y 4 = c 4.
16 CHAPTER 3. CONTROL TECHNIQUES Choroszucha 11 Figure 3.2: Solution Trajectories Figure 3.3 shows the control needed. Note, it never settles once activated. Figure 3.3: Control The code can be found in Appendix B.2. It should be noted, that if controls are present on all states then the nonlinearities can be canceled directly so that any trajectory is possible. 3.2 Sliding Mode Controller Sliding Mode Control (SMC) is a very promising control scheme. [23] has shown that SMC can be used to drive to an arbitrary point. Using a very elementary sliding surface for control of y, the author is able to do the same in a much less complicated manner by effectively changing r, but without the added x term multiplying the control as was discussed in [15]. S = y y e (3.10)
17 CHAPTER 3. CONTROL TECHNIQUES Choroszucha 12 The system goes rapidly to the desired point once the control is activated. The control is u 2 = ηsign(s) ẏ + 0.5sin(2πt) + y d. (3.11) Figure 3.4: States using SMC It should be noted in Figure 3.5 that the square is one equilibrium of the system and the + is the new equilibrium defined by y d = y e + 2. (3.12) As can be seen by the y vs τ plot of Figure 3.5, the state converges to 10.5 when y e 8.48.
18 CHAPTER 3. CONTROL TECHNIQUES Choroszucha 13 Figure 3.5: Solution Trajectories Figure 3.6 shows the control needed. Note, it settles to an approximately constant input. Figure 3.6: Control 3.3 Time-Delayed Feedback Control Time-Delayed Feedback Control (TDFC) is the dominant method in control of chaotic systems, however, it only guarantees the solution to settle onto an unstable periodic orbit (UPO, [16]) with a specified period, not a desired trajectory. TDFC is a proportional feedback controller with the feedback being delayed by a specified time. The general form for TDFC is where x(t) = f( x) + K x(t τ), (3.13)
19 CHAPTER 3. CONTROL TECHNIQUES Choroszucha 14 x(t) : States f( x) : Arbitrary Function of States K : Gain Matrix τ : Desired Period. Table 3.1: Variable Declarations for TDFC Some sources, e.g. [3] and [21], claim that the control will tend to zero for chaotic system using TDFC. There are two possible approaches one can take, either using this algorithm: 1. Declare initial conditions, separate time into a vector 2. Integrate system from time t i = t j to t f = t j+1 (a) When one period has passed and the Poincaré section condition is satisfied, implement control. (b) (Optional) Implement the optimal gain matrix [3] has proposed in the appendix. 3. Set initial conditions for next time step as being the final conditions from the previous time step. 4. Repeat until finished. or the method suggested in [19] using MATLAB s built in dde23 solver, the code for the first method can be found in Appendix B.4.1 and the dde23 has not been implemented properly, thus not shown here. Using a gain matrix K = 5I 3 3 (3.14) and a period of τ = 5, the following figures show that a periodic orbit is achieved and the control does go to zero.
20 CHAPTER 3. CONTROL TECHNIQUES Choroszucha 15 Figure 3.7: States using TDFC Figure 3.8: Solution Trajectories
21 CHAPTER 3. CONTROL TECHNIQUES Choroszucha 16 Figure 3.9: Control
22 Chapter 4 Conclusions 4.1 Concluding Remarks Theoretically, this analysis suggests that one may control the weather through controlling convection and temperatures. Realistically, controlling the full nonlinear partial differential equations would be incredibly difficult. It would take an amasing control scheme and possibly a very large amount of energy for input, or if time was not much of an issue, one could capitalise on the term that Edward Lorenz coined, the Butterfly Effect, whereby small inputs could drastically change the system, however, in acting on the system, the Butterfly Effect will ripple through, causing further issues. In a nutshell, this is a chicken and egg type problem. 4.2 Areas for Further Research Some tempting areas of research are: 1. A more robust tracking control for all states. 2. A more general control form and the corresponding Input-State Linearisation. 3. A mixture between TDFC and Input-State Linearisation controls. 4. Fuzzy control of chaotic systems. 5. Observers for chaotic systems. 17
23 References [1] John D. Anderson. Aerodynamics, Fundamentals of. Anderson Series. McGraw-Hill, Boston, 4th edition, [2] Serge Bielawski, Dominique Derozier, and Pierre Glorieux. Controlling unstable periodic orbits by a delayed continuous feedback. Phys. Rev. E, 49(2):R971 R974, Feb [3] James D. Biggs and Colin R. McInnes. Time-delayed feedback control in astrodynamics. Journal of Guidance, Control and Dynamics, 32(6): , November-December [4] S. Boccaletti, C. Grebogi, Y.-C. Lai, H. Mancini, and D. Maza. The control of chaos: Theory and applications. Physics Reports, 329:2000, [5] James Ward Brown and Ruel V. Churchill. Complex Variables and Applications. Mc- Graw Hill, Boston, 7th edition, [6] Jens Christian Claussen. Floquet stability analysis of ott-grebogi-yorke and difference control. New Journal of Physics, 10(6):063006, [7] B. Fiedler, V. Flunkert, M. Georgi, P. Hövel, and E. Schöll. Refuting the odd-number limitation of time-delayed feedback control. Phys. Rev. Lett., 98(11):114101, Mar [8] Xinping Guan and Cailian Chen. Adaptive fuzzy control for chaotic systems with h tracking performance. Fuzzy Sets and Systems, 139(1):81 93, [9] Richard Haberman. Applied Partial Differential Equations. Pearson - Prentice Hall, Upper Saddle River, 4th edition, [10] R. Hermann and A Krener. Nonlinear controllability and observability. IEEE Transactions on In Automatic Control, 22(5): , October [11] Robert C. Hilborn. Chaos and Nonlinear Dynamics. Oxford University Press, 2nd edition, [12] Chi-Chuan Hwang, Rong-Fong Fung, Jin-Yuan Hsieh, and Woei-Jong Li. A nonlinear feedback control of the lorenz equation. International Journal of Engineering Science, 37(14): ,
24 REFERENCES Choroszucha 19 [13] W. Just, B. Fiedler, M. Georgi, V. Flunkert, P. Hövel, and E. Schöll. Beyond the odd number limitation: A bifurcation analysis of time-delayed feedback control. Phys. Rev. E, 76(2):026210, Aug [14] Hassan Khalil. Nonlinear Systems. Prentice-Hall, Inc., Upper Saddle River, NJ, 3rd edition, [15] Chen Liqun and Liu Yanzhu. Control of the lorenz chaos by the exact linearization*. Applied Mathematics and Mechanics, 19(1):67 73, Jan [16] Hiroyuki Nakajima. On analytical properties of delayed feedback control of chaos. Physics Letters A, 232(3-4): , [17] Hiroyuki Nakajima and Yoshisuke Ueda. Limitation of generalized delayed feedback control. Physica D: Nonlinear Phenomena, 111(1-4): , [18] Edward Ott, Celso Grebogi, and James A. Yorke. Controlling chaos. Phys. Rev. Lett., 64(11): , Mar [19] Claire M. Postlethwaite and Mary Silber. Stabilizing unstable periodic orbits in the lorenz equations using time-delayed feedback control. Phys. Rev. E, 76(5):056214, Nov [20] K. Pyragas. Analytical properties and optimization of time-delayed feedback control. Phys. Rev. E, 66(2):026207, Aug [21] Kestutis Pyragas. Delayed feedback control of chaos. Philosophical Transactions of the Royal Society, 364: , July [22] Steven H. Strogatz. Nonlinear Dynamics and Chaos. Westview Press, Cambridge, [23] Shyi-Kae Yang, Chieh-Li Chen, and Her-Terng Yau. Control of chaos in lorenz system. Chaos, Solitons & Fractals, 13(4): , 2002.
25 Appendix A Derivation My derivation follows [11, ], but with added commentary using [9], [1] and [5] to expound upon the finer points in what [11] calls a development. 20
26 APPENDIX A. DERIVATION Choroszucha 21
27 APPENDIX A. DERIVATION Choroszucha 22
28 APPENDIX A. DERIVATION Choroszucha 23
29 APPENDIX A. DERIVATION Choroszucha 24
30 APPENDIX A. DERIVATION Choroszucha 25
31 Appendix B Code B.1 Open Loop Simulations 1 f u n c t i o n lorenzopen ( ) 2 3 sigma =10; 4 r =28; 5 b=8/3; 6 %t c =[0 5 0 ] ; 7 t c = 0 : : 5 0 ; 8 i c = [ ] ; 9 10 o d e o p t s=o d e s e t ( R e l T o l, 1E 6); 11 [ t x]= ode45 tc, i c, o d e o p t s ) ; 12 f i g u r e ( 1 ) 13 s u b p l o t ( 2, 2, 1 ) 14 p l o t 3 ( x ( :, 1 ), x ( :, 2 ), x ( :, 3 ), k ) 15 x l a b e l ( x ) 16 y l a b e l ( y ) 17 z l a b e l ( z ) 18 s u b p l o t ( 2, 2, 2 ) 19 p l o t ( t, x ( :, 1 ), k ) 20 x l a b e l ( \ t a u ) 21 y l a b e l ( x ) 22 s u b p l o t ( 2, 2, 3 ) 23 p l o t ( t, x ( :, 2 ), k ) 24 x l a b e l ( \ t a u ) 25 y l a b e l ( y ) 26 s u b p l o t ( 2, 2, 4 ) 27 p l o t ( t, x ( :, 3 ), k ) 28 x l a b e l ( \ t a u ) 29 y l a b e l ( z ) 30 s a v e a s ( gca, o p e n l o o p A l l - R e l T o l -6 - x p n g ) f i g u r e ( 2 ) 33 p l o t 3 ( x ( :, 1 ), x ( :, 2 ), x ( :, 3 ), k ) 34 x l a b e l ( x ) 35 y l a b e l ( y ) 36 z l a b e l ( z ) 37 s a v e a s ( gca, o p e n l o o p - R e l T o l -6 - x p n g ) f u n c t i o n dx=lrwhcw ( t, x ) dx=z e r o s ( 3, 1 ) ; 42 dx (1)= sigma ( x(2) x ( 1 ) ) ; 43 dx (2)= r x(1) x(2) x ( 1 ) x ( 3 ) ; 44 dx (3)= x ( 1 ) x(2) b x ( 3 ) ; 45 end 46 end B.2 Input State Linearisation 1 f u n c t i o n i s l i n e a r i s a t i o n ( ) 2 3 numpoints =10000; 26
32 APPENDIX B. CODE Choroszucha sigma =10; 6 r =28; 7 b=8/3; 8 t f =100; 9 10 beta=[ 1 3 3]; t t o t=l i n s p a c e ( 0, t f, numpoints ) ; 13 nt=numel ( t t o t ) ; 14 i c =[ ] ; 15 o d e o p t s=o d e s e t ( R e l T o l, 1E 9); 16 [ t x]= ode45 t t o t, i c, o d e o p t s ) ; f i g u r e ( 1 ) 19 s u b p l o t ( 2, 2, 1 ) 20 p l o t 3 ( x ( :, 1 ), x ( :, 2 ), x ( :, 3 ), k ) 21 x l a b e l ( x ) 22 y l a b e l ( y ) 23 z l a b e l ( z ) 24 s u b p l o t ( 2, 2, 2 ) 25 p l o t ( t, x ( :, 1 ), k ) 26 x l a b e l ( \ t a u ) 27 y l a b e l ( x ) 28 s u b p l o t ( 2, 2, 3 ) 29 p l o t ( t, x ( :, 2 ), k ) 30 x l a b e l ( \ t a u ) 31 y l a b e l ( y ) 32 s u b p l o t ( 2, 2, 4 ) 33 p l o t ( t, x ( :, 3 ), k ) 34 x l a b e l ( \ t a u ) 35 y l a b e l ( z ) 36 s a v e a s ( gca, i s f l - a l l. p n g ) f i g u r e ( 2 ) 39 p l o t 3 ( x ( :, 1 ), x ( :, 2 ), x ( :, 3 ), k ) 40 x l a b e l ( x ) 41 y l a b e l ( y ) 42 z l a b e l ( z ) 43 s a v e a s ( gca, i s f l. p n g ) f i g u r e ( 3 ) 46 p l o t ( t ( 1 : ( end 1)), d i f f ( x ( :, 4 ) ). / d i f f ( t ), k ) 47 x l a b e l ( \ t a u ) 48 y l a b e l ( u ) 49 s a v e a s ( gca, i s f l - c o n t r o l. p n g ) function dx=sonofthemorning ( t, x ) dx=z e r o s ( 4, 1 ) ; 54 u=z e r o s ( 3, 1 ) ; 55 F=[ sigma ( x(2) x ( 1 ) ) r x(1) x(2) x ( 1 ) x ( 3 ) x ( 1 ) x(2) b x ( 3 ) ] ; 56 i f ( t >20) f =2 p i ; 59 ampl =500; 60 x1d=(ampl s i n ( f t ) ) ; 61 x1ddot=ampl f c o s ( f t ) ; 62 x1dddot= ampl f ˆ2 s i n ( f t ) ; 63 x1ddddot= ampl f ˆ3 c o s ( f t ) ; 64 xd=[x1d x1ddot x1dddot ] ; d=(2 sigma b ) / ( b sigma ) ; 67 mu=d x1d ˆ2; z (1,1)=1/ sigma x (1)ˆ2 2 x (3)+mu; 70 z (2,1)= 2 x (1)ˆ2+2 b x ( 3 ) ; 71 z (3,1)= 4 sigma x ( 1 ) ( x(2) x (1))+2 b ( x ( 1 ) x(2) b x ( 3 ) ) ; 72 v=beta ( z ) ;%x1ddddot+ d xd ) ; 73 c=(b 2 sigma ) ; 74 D3Y=[4 sigma ( x(2) 2 x (1))+2 b x ( 2 ) 4 sigma x (2)+2 b x ( 1 ) 2 b ˆ 2 ] ; 75 numer=d3y F ; 76 denom=x (1)ˆ2 ( 4 sigma+2 b ) ; 77 u(2)=1/denom ( v numer ) ; 78 e l s e 79 u=z e r o s ( 3, 1 ) ; 80 end 81 dx=z e r o s ( 3, 1 ) ; 82 dx ( 1 : 3 ) =F+u ; 83 dx (4)=u ( 2 ) ; 84 end end B.3 Sliding Mode Control
33 APPENDIX B. CODE Choroszucha 28 1 function elemsmc ( ) 2 3 numpoints =10000; 4 5 e t a =0.25; 6 sigma =10; 7 r =28; 8 b=8/3; 9 e q p o i n t=s q r t ( b ( r 1 ) ) ; 10 eqvect =[ eqpoint eqpoint r 1 ] ; t f =100; 13 t t o t=l i n s p a c e ( 0, t f, numpoints ) ; 14 nt=numel ( t t o t ) ; 15 i c =[ ] ; 16 o d e o p t s=o d e s e t ( R e l T o l, 1E 9); 17 [ t x]= ode45 t t o t, i c, o d e o p t s ) ; f i g u r e ( 1 ) 20 s u b p l o t ( 2, 2, 1 ) 21 p l o t 3 ( x ( :, 1 ), x ( :, 2 ), x ( :, 3 ), k ) 22 x l a b e l ( x ) 23 y l a b e l ( y ) 24 z l a b e l ( z ) 25 s u b p l o t ( 2, 2, 2 ) 26 p l o t ( t, x ( :, 1 ), k ) 27 x l a b e l ( \ t a u ) 28 y l a b e l ( x ) 29 s u b p l o t ( 2, 2, 3 ) 30 p l o t ( t, x ( :, 2 ), k ) 31 x l a b e l ( \ t a u ) 32 y l a b e l ( y ) 33 s u b p l o t ( 2, 2, 4 ) 34 p l o t ( t, x ( :, 3 ), k ) 35 x l a b e l ( \ t a u ) 36 y l a b e l ( z ) 37 s a v e a s ( gca, smc - a l l. p n g ) f i g u r e ( 2 ) 40 p l o t 3 ( x ( :, 1 ), x ( :, 2 ), x ( :, 3 ), k ) 41 hold on 42 p l o t 3 ( x ( end, 1 ), x ( end, 2 ), x ( end, 3 ), k +, M a r k e r S i z e, 1 2 ) 43 p l o t 3 ( e q v e c t ( 1 ), e q v e c t ( 2 ), e q v e c t ( 3 ), k s, M a r k e r S i z e, 1 2 ) 44 x l a b e l ( x ) 45 y l a b e l ( y ) 46 z l a b e l ( z ) 47 s a v e a s ( gca, s m c. p n g ) f i g u r e ( 3 ) 50 p l o t ( t ( 1 : ( end 1)), d i f f ( x ( :, 4 ) ). / d i f f ( t ), k ) 51 x l a b e l ( \ t a u ) 52 y l a b e l ( u ) 53 s a v e a s ( gca, smc - c o n t r o l. p n g ) function dx=sonofthemorning ( t, x ) dx=z e r o s ( 4, 1 ) ; 58 u=z e r o s ( 3, 1 ) ; 59 F=[ sigma ( x(2) x ( 1 ) ) ( r x ( 3 ) ) x(1) x ( 2 ) x ( 1 ) x(2) b x ( 3 ) ] ; i f ( t >20) 62 yd=eqpoint +2; 63 S=x(2) yd ; 64 u(2)= ( e t a ) s i g n ( S) F( 2 ) s i n (2 p i t )+yd ; 65 end dx=z e r o s ( 3, 1 ) ; 68 dx ( 1 : 3 ) =F+u ; 69 dx (4)=u ( 2 ) ; 70 end 71 end B.4 Time-Delayed Feedback Control B.4.1 Time-Step Algorithm 1 f u n c t i o n stepalgorithm ( ) 2 3 c l o s e a l l ; 4 pause ( 0. 1 ) ; 5 6 f o l d e r n a m e=d a t e s t r ( now, 3 0 ) ; 7 note={ L o o k i n g a t d o u b l e o r b i t. } ; 8 numpoints =10000;
34 APPENDIX B. CODE Choroszucha type= s t e p ; 11 sigma =10; 12 r =28; 13 b=8/3; 14 %t c =[0 5 0 ] ; 15 t f =100; tauarr =10; 18 K=1.5 eye ( 3, 3 ) ;%+0.1 rand ( 3, 3 ) ; f o r o1=1 21 tau=t a u a r r ( o1 ) ; 22 t t o t=l i n s p a c e ( 0, t f, numpoints ) ; 23 nt=numel ( t t o t ) ; 24 x t o t=z e r o s ( nt, 6 ) ; 25 i c =[ ] ; 26 x t o t ( 1, : ) = i c ; 27 t c =[ t t o t ( 1 ) t t o t ( 2 ) ] ; 28 o d e o p t s=o d e s e t ( R e l T o l, 1E 9); 29 f o r a1 =1:( nt 1) 30 i f (mod( a1,100)==0) 31 f p r i n t f ( % d / % d \ n, a1, nt ) ; 32 end 33 [ t x]= ode45 tc, i c, o d e o p t s ) ; 34 x t o t ( a1 +1,:)= x ( end, : ) ; 35 i c=x ( end, : ) ; 36 i f ( a1 =nt 1) 37 t c =[ t t o t ( a1+1) t t o t ( a1 + 2 ) ] ; 38 end 39 end s a v e C h a o t i c P l o t s ( type, t t o t, xtot, f a l s e, foldername, note, tau,k) ; 42 end function dx=sonofthemorning ( t, x ) dx=z e r o s ( 6, 1 ) ; 47 f o r i 1 =1:3 48 e ( i 1,1)= calcu ( t, x, i 1 ) ; 49 end 50 u= K e ; dx=z e r o s ( 3, 1 ) ; 53 dx (1)= sigma ( x(2) x (1))+ u ( 1 ) ; 54 dx (2)= r x(1) x(2) x ( 1 ) x (3)+u ( 2 ) ; 55 dx (3)= x ( 1 ) x(2) b x (3)+u ( 3 ) ; 56 dx ( 4 : 6 ) =u ; 57 end f u n c t i o n myu=calcu ( t, x, a2 ) 60 i f ( t>=tau ) 61 myu=x ( a2) i n t e r p 1 ( t t o t, x t o t ( :, a2 ), t tau ) ; 62 e l s e 63 myu=x ( a2 ) ; 64 end 65 end 66 end B.4.2 Delayed Differential Equation 1 f u n c t i o n TDFCdde ( ) 2 3 sigma =10; 4 r =28; 5 b=8/3; 6 %t c =[0 5 0 ] ; 7 t f =100; 8 e q p o i n t=s q r t ( b ( r 1 ) ) ; 9 10 tau =10; 11 K=5 eye ( 3, 3 ) ; t f=tau 1 0 ; 14 t t o t=l i n s p a c e ( 0, t f, ) ; 15 nt=numel ( t t o t ) ; 16 x t o t=z e r o s ( nt, 6 ) ; 17 i c =[ ] ; 18 x t o t ( 1, : ) = i c ; o d e o p t s=o d e s e t ( R e l T o l, 1E 9); 21 ddeopts=d d e s e t ( R e l T o l, 1E 3); p e r T i c =[ t t o t ( 1 ) t t o t (1)+ tau ] ; 24 [ pert h e r s t o r y ]= ode45 [ t t o t ( 1 ) tau+t t o t ( 1 ) ], i c ( 1 : 3 ), o d e o p t s ) ; 25 s o l=dde23 [ min ( t t o t ) max( t t o t ) ], ddeopts ) ; 26 x t o t=d e v a l ( s o l, t t o t ) ; 27
35 APPENDIX B. CODE Choroszucha f i g u r e ( 1 ) 29 p l o t ( x t o t ( :, 1 ), x t o t ( :, 2 ), k ) 30 x l a b e l ( x ) 31 y l a b e l ( y ) 32 s a v e a s ( gca, s p r i n t f ( b - xy - % s - % d. p n g, uwhat, rho ) ) f i g u r e ( 3 ) 35 p l o t ( t t o t ( 1 : ( end 1)), d i f f ( x t o t ( :, 5 ) ). / d i f f ( t t o t ), k ) 36 x l a b e l ( t [ s ] ) 37 y l a b e l ( u _ x ) 38 s a v e a s ( gca, s p r i n t f ( b - ux - % s - % d. p n g, uwhat, rho ) ) f i g u r e ( 4 ) 41 p l o t ( t t o t ( 1 : ( end 1)), d i f f ( x t o t ( :, 6 ) ). / d i f f ( t t o t ), k ) 42 x l a b e l ( t [ s ] ) 43 y l a b e l ( u _ y ) 44 s a v e a s ( gca, s p r i n t f ( b - uy - % s - % d. p n g, uwhat, rho ) ) 45 pause ( 1 ) ; 46 c l o s e a l l ; function dx=sonofthemorning ( t, x ) dx=z e r o s ( 3, 1 ) ; 51 dx (1)= sigma ( x(2) x ( 1 ) ) ; 52 dx (2)= r x(1) x(2) x ( 1 ) x ( 3 ) ; 53 dx (3)= x ( 1 ) x(2) b x ( 3 ) ; 54 end function dx=daughterofthemorning ( t, x, xtau ) dx=z e r o s ( 6, 1 ) ; 59 e=x (1:3) xtau ( 1 : 3 ) ; 60 u= K e ; r=x (1)ˆ2+ x ( 2 ) ˆ 2 ; 63 dx (1)= sigma ( x(2) x ( 1 ) ) ; 64 dx (2)= r x(1) x(2) x ( 1 ) x (3)+u ; 65 dx (3)= x ( 1 ) x(2) b x ( 3 ) ; 66 dx (4)=u ; 67 end f u n c t i o n myu=calcu ( t, x, a2 ) 70 i f ( t>=tau ) 71 myu=x ( a2) i n t e r p 1 ( t t o t, x t o t ( :, a2 ), t tau ) ; 72 e l s e 73 myu=0; 74 end 75 end function r e t=herherstory ( t ) 78 r e t ( 1 : 3 ) = i n t e r p 1 ( pert, h e r s t o r y, mod( t, tau ) ) ; 79 r e t ( 4 ) = 0 ; 80 end 81 end
Chaos Control for the Lorenz System
Advanced Studies in Theoretical Physics Vol. 12, 2018, no. 4, 181-188 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2018.8413 Chaos Control for the Lorenz System Pedro Pablo Cárdenas Alzate
More informationarxiv: v2 [nlin.cd] 8 Sep 2012
An analytical limitation for time-delayed feedback control in autonomous systems Edward W. Hooton 1 1, 2, and Andreas Amann 1 School of Mathematical Sciences, University College Cork, Ireland 2 Tyndall
More informationA Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term
ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More informationEdward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology
The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More informationResearch Article Adaptive Control of Chaos in Chua s Circuit
Mathematical Problems in Engineering Volume 2011, Article ID 620946, 14 pages doi:10.1155/2011/620946 Research Article Adaptive Control of Chaos in Chua s Circuit Weiping Guo and Diantong Liu Institute
More informationLecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:
Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationSimplest Chaotic Flows with Involutional Symmetries
International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries
More informationDynamical analysis and circuit simulation of a new three-dimensional chaotic system
Dynamical analysis and circuit simulation of a new three-dimensional chaotic system Wang Ai-Yuan( 王爱元 ) a)b) and Ling Zhi-Hao( 凌志浩 ) a) a) Department of Automation, East China University of Science and
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationControlling a Novel Chaotic Attractor using Linear Feedback
ISSN 746-7659, England, UK Journal of Information and Computing Science Vol 5, No,, pp 7-4 Controlling a Novel Chaotic Attractor using Linear Feedback Lin Pan,, Daoyun Xu 3, and Wuneng Zhou College of
More informationChaos Control of the Chaotic Symmetric Gyroscope System
48 Chaos Control of the Chaotic Symmetric Gyroscope System * Barış CEVHER, Yılmaz UYAROĞLU and 3 Selçuk EMIROĞLU,,3 Faculty of Engineering, Department of Electrical and Electronics Engineering Sakarya
More informationGLOBAL CHAOS SYNCHRONIZATION OF UNCERTAIN SPROTT J AND K SYSTEMS BY ADAPTIVE CONTROL
GLOBAL CHAOS SYNCHRONIZATION OF UNCERTAIN SPROTT J AND K SYSTEMS BY ADAPTIVE CONTROL Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University Avadi,
More informationADAPTIVE CONTROL AND SYNCHRONIZATION OF HYPERCHAOTIC NEWTON-LEIPNIK SYSTEM
ADAPTIVE CONTROL AND SYNCHRONIZATION OF HYPERCHAOTIC NEWTON-LEIPNIK SYSTEM Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University Avadi, Chennai-600
More informationHYBRID CHAOS SYNCHRONIZATION OF HYPERCHAOTIC LIU AND HYPERCHAOTIC CHEN SYSTEMS BY ACTIVE NONLINEAR CONTROL
HYBRID CHAOS SYNCHRONIZATION OF HYPERCHAOTIC LIU AND HYPERCHAOTIC CHEN SYSTEMS BY ACTIVE NONLINEAR CONTROL Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical
More informationTime-delay feedback control in a delayed dynamical chaos system and its applications
Time-delay feedback control in a delayed dynamical chaos system and its applications Ye Zhi-Yong( ), Yang Guang( ), and Deng Cun-Bing( ) School of Mathematics and Physics, Chongqing University of Technology,
More informationCONTROLLING CHAOTIC DYNAMICS USING BACKSTEPPING DESIGN WITH APPLICATION TO THE LORENZ SYSTEM AND CHUA S CIRCUIT
Letters International Journal of Bifurcation and Chaos, Vol. 9, No. 7 (1999) 1425 1434 c World Scientific Publishing Company CONTROLLING CHAOTIC DYNAMICS USING BACKSTEPPING DESIGN WITH APPLICATION TO THE
More informationADAPTIVE CONTROL AND SYNCHRONIZATION OF A GENERALIZED LOTKA-VOLTERRA SYSTEM
ADAPTIVE CONTROL AND SYNCHRONIZATION OF A GENERALIZED LOTKA-VOLTERRA SYSTEM Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University Avadi, Chennai-600
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationGLOBAL CHAOS SYNCHRONIZATION OF UNCERTAIN LORENZ-STENFLO AND QI 4-D CHAOTIC SYSTEMS BY ADAPTIVE CONTROL
GLOBAL CHAOS SYNCHRONIZATION OF UNCERTAIN LORENZ-STENFLO AND QI 4-D CHAOTIC SYSTEMS BY ADAPTIVE CONTROL Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical
More information6.2 Brief review of fundamental concepts about chaotic systems
6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3-variable model that is a prototypical example of chaos theory. These equations were derived as a simplification
More informationThe projects listed on the following pages are suitable for MSc/MSci or PhD students. An MSc/MSci project normally requires a review of the
The projects listed on the following pages are suitable for MSc/MSci or PhD students. An MSc/MSci project normally requires a review of the literature and finding recent related results in the existing
More informationDynamical behaviour of a controlled vibro-impact system
Vol 17 No 7, July 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(07)/2446-05 Chinese Physics B and IOP Publishing Ltd Dynamical behaviour of a controlled vibro-impact system Wang Liang( ), Xu Wei( ), and
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationMULTISTABILITY IN A BUTTERFLY FLOW
International Journal of Bifurcation and Chaos, Vol. 23, No. 12 (2013) 1350199 (10 pages) c World Scientific Publishing Company DOI: 10.1142/S021812741350199X MULTISTABILITY IN A BUTTERFLY FLOW CHUNBIAO
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationA Novel Hyperchaotic System and Its Control
1371371371371378 Journal of Uncertain Systems Vol.3, No., pp.137-144, 009 Online at: www.jus.org.uk A Novel Hyperchaotic System and Its Control Jiang Xu, Gouliang Cai, Song Zheng School of Mathematics
More informationFinite-time hybrid synchronization of time-delay hyperchaotic Lorenz system
ISSN 1746-7659 England UK Journal of Information and Computing Science Vol. 10 No. 4 2015 pp. 265-270 Finite-time hybrid synchronization of time-delay hyperchaotic Lorenz system Haijuan Chen 1 * Rui Chen
More informationADAPTIVE DESIGN OF CONTROLLER AND SYNCHRONIZER FOR LU-XIAO CHAOTIC SYSTEM
ADAPTIVE DESIGN OF CONTROLLER AND SYNCHRONIZER FOR LU-XIAO CHAOTIC SYSTEM WITH UNKNOWN PARAMETERS Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 533 Homeworks Spring 07 Updated: Saturday, April 08, 07 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. Some homework assignments refer to the textbooks: Slotine
More informationGLOBAL CHAOS SYNCHRONIZATION OF HYPERCHAOTIC QI AND HYPERCHAOTIC JHA SYSTEMS BY ACTIVE NONLINEAR CONTROL
GLOBAL CHAOS SYNCHRONIZATION OF HYPERCHAOTIC QI AND HYPERCHAOTIC JHA SYSTEMS BY ACTIVE NONLINEAR CONTROL Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical
More informationGenerating a Complex Form of Chaotic Pan System and its Behavior
Appl. Math. Inf. Sci. 9, No. 5, 2553-2557 (2015) 2553 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/090540 Generating a Complex Form of Chaotic Pan
More informationTHE CONTROL OF CHAOS: THEORY AND APPLICATIONS
S. Boccaletti et al. / Physics Reports 329 (2000) 103}197 103 THE CONTROL OF CHAOS: THEORY AND APPLICATIONS S. BOCCALETTI, C. GREBOGI, Y.-C. LAI, H. MANCINI, D. MAZA Department of Physics and Applied Mathematics,
More informationPH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract
PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Mr. Benjamen P. Reed (110108461) IMPACS, Aberystwyth University January 31, 2014 Abstract A set of three coupled ordinary
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationComplete Synchronization, Anti-synchronization and Hybrid Synchronization Between Two Different 4D Nonlinear Dynamical Systems
Mathematics Letters 2016; 2(5): 36-41 http://www.sciencepublishinggroup.com/j/ml doi: 10.11648/j.ml.20160205.12 Complete Synchronization, Anti-synchronization and Hybrid Synchronization Between Two Different
More informationIntroduction Knot Theory Nonlinear Dynamics Topology in Chaos Open Questions Summary. Topology in Chaos
Introduction Knot Theory Nonlinear Dynamics Open Questions Summary A tangled tale about knot, link, template, and strange attractor Centre for Chaos & Complex Networks City University of Hong Kong Email:
More informationDiscussion of the Lorenz Equations
Discussion of the Lorenz Equations Leibniz Universität Hannover Proseminar Theoretische Physik SS/2015 July 22, 2015 (LUH) Lorenz Equations July 22, 2015 1 / 42 Outline 1 2 3 4 5 6 7 8 (LUH) Lorenz Equations
More informationDynamical Behavior And Synchronization Of Chaotic Chemical Reactors Model
Iranian Journal of Mathematical Chemistry, Vol. 6, No. 1, March 2015, pp. 81 92 IJMC Dynamical Behavior And Synchronization Of Chaotic Chemical Reactors Model HOSSEIN KHEIRI 1 AND BASHIR NADERI 2 1 Faculty
More informationMultistability in the Lorenz System: A Broken Butterfly
International Journal of Bifurcation and Chaos, Vol. 24, No. 10 (2014) 1450131 (7 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414501314 Multistability in the Lorenz System: A Broken
More informationCONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS
International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems
More informationGeneralized projective synchronization of a class of chaotic (hyperchaotic) systems with uncertain parameters
Vol 16 No 5, May 2007 c 2007 Chin. Phys. Soc. 1009-1963/2007/16(05)/1246-06 Chinese Physics and IOP Publishing Ltd Generalized projective synchronization of a class of chaotic (hyperchaotic) systems with
More informationADAPTIVE CHAOS SYNCHRONIZATION OF UNCERTAIN HYPERCHAOTIC LORENZ AND HYPERCHAOTIC LÜ SYSTEMS
ADAPTIVE CHAOS SYNCHRONIZATION OF UNCERTAIN HYPERCHAOTIC LORENZ AND HYPERCHAOTIC LÜ SYSTEMS Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University
More informationMATH 415, WEEK 12 & 13: Higher-Dimensional Systems, Lorenz Equations, Chaotic Behavior
MATH 415, WEEK 1 & 13: Higher-Dimensional Systems, Lorenz Equations, Chaotic Behavior 1 Higher-Dimensional Systems Consider the following system of differential equations: dx = x y dt dy dt = xy y dz dt
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationThe Existence of Chaos in the Lorenz System
The Existence of Chaos in the Lorenz System Sheldon E. Newhouse Mathematics Department Michigan State University E. Lansing, MI 48864 joint with M. Berz, K. Makino, A. Wittig Physics, MSU Y. Zou, Math,
More informationChaos in multiplicative systems
Chaos in multiplicative systems Dorota Aniszewska 1 and Marek Rybaczuk 2 1 Institute of Materials Science and Applied Mechanics Wroclaw University of Technology 50-370 Wroclaw, Smoluchowskiego 25 (e-mail:
More informationON STABILIZING N-DIMENSIONAL CHAOTIC SYSTEMS
International Journal of Bifurcation and Chaos, Vol. 3, No. 2 (23) 473 48 c World Scientific Publishing Company ON STABILIZING N-DIMENSIONAL CHAOTIC SYSTEMS LAURENT LAVAL and NACER K. M SIRDI Laboratoire
More informationAdaptive feedback synchronization of a unified chaotic system
Physics Letters A 39 (4) 37 333 www.elsevier.com/locate/pla Adaptive feedback synchronization of a unified chaotic system Junan Lu a, Xiaoqun Wu a, Xiuping Han a, Jinhu Lü b, a School of Mathematics and
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationResearch Article Hopf Bifurcation Analysis and Anticontrol of Hopf Circles of the Rössler-Like System
Abstract and Applied Analysis Volume, Article ID 3487, 6 pages doi:.55//3487 Research Article Hopf Bifurcation Analysis and Anticontrol of Hopf Circles of the Rössler-Like System Ranchao Wu and Xiang Li
More informationBidirectional Partial Generalized Synchronization in Chaotic and Hyperchaotic Systems via a New Scheme
Commun. Theor. Phys. (Beijing, China) 45 (2006) pp. 1049 1056 c International Academic Publishers Vol. 45, No. 6, June 15, 2006 Bidirectional Partial Generalized Synchronization in Chaotic and Hyperchaotic
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationIntroduction to Dynamical Systems Basic Concepts of Dynamics
Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic
More informationSPATIOTEMPORAL CHAOS IN COUPLED MAP LATTICE. Itishree Priyadarshini. Prof. Biplab Ganguli
SPATIOTEMPORAL CHAOS IN COUPLED MAP LATTICE By Itishree Priyadarshini Under the Guidance of Prof. Biplab Ganguli Department of Physics National Institute of Technology, Rourkela CERTIFICATE This is to
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps
More informationHYBRID CHAOS SYNCHRONIZATION OF UNCERTAIN LORENZ-STENFLO AND QI 4-D CHAOTIC SYSTEMS BY ADAPTIVE CONTROL
HYBRID CHAOS SYNCHRONIZATION OF UNCERTAIN LORENZ-STENFLO AND QI 4-D CHAOTIC SYSTEMS BY ADAPTIVE CONTROL Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical
More informationADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS
Letters International Journal of Bifurcation and Chaos, Vol. 12, No. 7 (2002) 1579 1597 c World Scientific Publishing Company ADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS A. S. HEGAZI,H.N.AGIZA
More informationUSING DYNAMIC NEURAL NETWORKS TO GENERATE CHAOS: AN INVERSE OPTIMAL CONTROL APPROACH
International Journal of Bifurcation and Chaos, Vol. 11, No. 3 (2001) 857 863 c World Scientific Publishing Company USING DYNAMIC NEURAL NETWORKS TO GENERATE CHAOS: AN INVERSE OPTIMAL CONTROL APPROACH
More informationSynchronizing Chaotic Systems Based on Tridiagonal Structure
Proceedings of the 7th World Congress The International Federation of Automatic Control Seoul, Korea, July 6-, 008 Synchronizing Chaotic Systems Based on Tridiagonal Structure Bin Liu, Min Jiang Zengke
More informationA State-Space Model for a Nonlinear Time-Delayed Feedback Loop
A for a Nonlinear Time-Delayed Feedback Loop Karl Schmitt Advisors: Jim Yorke, Rajarshi Roy, Tom Murphy AMSC 663 October 14, 2008 1 / 21 Goal To implement an alternative, discrete time model for coupled
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationA New Hyperchaotic Attractor with Complex Patterns
A New Hyperchaotic Attractor with Complex Patterns Safieddine Bouali University of Tunis, Management Institute, Department of Quantitative Methods & Economics, 41, rue de la Liberté, 2000, Le Bardo, Tunisia
More informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
More informationOUTPUT REGULATION OF THE SIMPLIFIED LORENZ CHAOTIC SYSTEM
OUTPUT REGULATION OF THE SIMPLIFIED LORENZ CHAOTIC SYSTEM Sundarapandian Vaidyanathan Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University Avadi, Chennai-600 06, Tamil Nadu, INDIA
More informationTHE ACTIVE CONTROLLER DESIGN FOR ACHIEVING GENERALIZED PROJECTIVE SYNCHRONIZATION OF HYPERCHAOTIC LÜ AND HYPERCHAOTIC CAI SYSTEMS
THE ACTIVE CONTROLLER DESIGN FOR ACHIEVING GENERALIZED PROJECTIVE SYNCHRONIZATION OF HYPERCHAOTIC LÜ AND HYPERCHAOTIC CAI SYSTEMS Sarasu Pakiriswamy 1 and Sundarapandian Vaidyanathan 1 1 Department of
More informationA new four-dimensional chaotic system
Chin. Phys. B Vol. 19 No. 12 2010) 120510 A new four-imensional chaotic system Chen Yong ) a)b) an Yang Yun-Qing ) a) a) Shanghai Key Laboratory of Trustworthy Computing East China Normal University Shanghai
More informationExperimental Characterization of Nonlinear Dynamics from Chua s Circuit
Experimental Characterization of Nonlinear Dynamics from Chua s Circuit John Parker*, 1 Majid Sodagar, 1 Patrick Chang, 1 and Edward Coyle 1 School of Physics, Georgia Institute of Technology, Atlanta,
More informationNonlinear System Analysis
Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control
More informationCalculating Fractal Dimension of Attracting Sets of the Lorenz System
Dynamics at the Horsetooth Volume 6, 2014. Calculating Fractal Dimension of Attracting Sets of the Lorenz System Jamie Department of Mathematics Colorado State University Report submitted to Prof. P. Shipman
More informationSimple conservative, autonomous, second-order chaotic complex variable systems.
Simple conservative, autonomous, second-order chaotic complex variable systems. Delmar Marshall 1 (Physics Department, Amrita Vishwa Vidyapeetham, Clappana P.O., Kollam, Kerala 690-525, India) and J. C.
More informationADAPTIVE STABILIZATION AND SYNCHRONIZATION OF HYPERCHAOTIC QI SYSTEM
ADAPTIVE STABILIZATION AND SYNCHRONIZATION OF HYPERCHAOTIC QI SYSTEM Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR Dr. SR Technical University Avadi, Chennai-600 062,
More informationADAPTIVE CONTROLLER DESIGN FOR THE ANTI-SYNCHRONIZATION OF HYPERCHAOTIC YANG AND HYPERCHAOTIC PANG SYSTEMS
ADAPTIVE CONTROLLER DESIGN FOR THE ANTI-SYNCHRONIZATION OF HYPERCHAOTIC YANG AND HYPERCHAOTIC PANG SYSTEMS Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical
More informationGLOBAL CHAOS SYNCHRONIZATION OF PAN AND LÜ CHAOTIC SYSTEMS VIA ADAPTIVE CONTROL
GLOBAL CHAOS SYNCHRONIZATION OF PAN AND LÜ CHAOTIC SYSTEMS VIA ADAPTIVE CONTROL Sundarapandian Vaidyanathan 1 and Karthikeyan Rajagopal 2 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical
More informationEen vlinder in de wiskunde: over chaos en structuur
Een vlinder in de wiskunde: over chaos en structuur Bernard J. Geurts Enschede, November 10, 2016 Tuin der Lusten (Garden of Earthly Delights) In all chaos there is a cosmos, in all disorder a secret
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationHandout 2: Invariant Sets and Stability
Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state
More informationDelayed Feedback Controller for Stabilizing Subsynchronous Oscillations in Power Systems
Delayed Feedback Controller for Stabilizing Subsynchronous Oscillations in Power Systems Yaar Küçükefe 1 and Adnan Kaypmaz 2 1 National Power Enerji, Tekirda, Turkey yasar.kucukefe@ieee.org 2 stanbul Technical
More informationA GALLERY OF LORENZ-LIKE AND CHEN-LIKE ATTRACTORS
International Journal of Bifurcation and Chaos, Vol. 23, No. 4 (2013) 1330011 (20 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127413300115 A GALLERY OF LORENZ-LIKE AND CHEN-LIKE ATTRACTORS
More informationAPPM 2460 CHAOTIC DYNAMICS
APPM 2460 CHAOTIC DYNAMICS 1. Introduction Today we ll continue our exploration of dynamical systems, focusing in particular upon systems who exhibit a type of strange behavior known as chaos. We will
More informationSYNCHRONIZATION CRITERION OF CHAOTIC PERMANENT MAGNET SYNCHRONOUS MOTOR VIA OUTPUT FEEDBACK AND ITS SIMULATION
SYNCHRONIZAION CRIERION OF CHAOIC PERMANEN MAGNE SYNCHRONOUS MOOR VIA OUPU FEEDBACK AND IS SIMULAION KALIN SU *, CHUNLAI LI College of Physics and Electronics, Hunan Institute of Science and echnology,
More informationConstructing a chaotic system with any number of equilibria
Nonlinear Dyn (2013) 71:429 436 DOI 10.1007/s11071-012-0669-7 ORIGINAL PAPER Constructing a chaotic system with any number of equilibria Xiong Wang Guanrong Chen Received: 9 June 2012 / Accepted: 29 October
More informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More information3. Controlling the time delay hyper chaotic Lorenz system via back stepping control
ISSN 1746-7659, England, UK Journal of Information and Computing Science Vol 10, No 2, 2015, pp 148-153 Chaos control of hyper chaotic delay Lorenz system via back stepping method Hanping Chen 1 Xuerong
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationK. Pyragas* Semiconductor Physics Institute, LT-2600 Vilnius, Lithuania Received 19 March 1998
PHYSICAL REVIEW E VOLUME 58, NUMBER 3 SEPTEMBER 998 Synchronization of coupled time-delay systems: Analytical estimations K. Pyragas* Semiconductor Physics Institute, LT-26 Vilnius, Lithuania Received
More informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationStrange dynamics of bilinear oscillator close to grazing
Strange dynamics of bilinear oscillator close to grazing Ekaterina Pavlovskaia, James Ing, Soumitro Banerjee and Marian Wiercigroch Centre for Applied Dynamics Research, School of Engineering, King s College,
More informationImpulsive Stabilization for Control and Synchronization of Chaotic Systems: Theory and Application to Secure Communication
976 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: FUNDAMENTAL THEORY AND APPLICATIONS, VOL. 44, NO. 10, OCTOBER 1997 Impulsive Stabilization for Control and Synchronization of Chaotic Systems: Theory and
More informationConverse Lyapunov theorem and Input-to-State Stability
Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts
More informationChaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel Nonlinear dynamics course - VUB
Chaos Lendert Gelens KU Leuven - Vrije Universiteit Brussel www.gelenslab.org Nonlinear dynamics course - VUB Examples of chaotic systems: the double pendulum? θ 1 θ θ 2 Examples of chaotic systems: the
More informationMore Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n.
More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n. If there are points which, after many iterations of map then fixed point called an attractor. fixed point, If λ
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More information