PRAMANA
c Indian Academy of Sciences
— journal of physics
Vol. 80, No. 5 May 2013 pp. 825–835
Comparison of backstepping and modified active control in projective synchronization of chaos in an extended Bonhöffer–van der Pol oscillator K S OJO1,2,∗ , A N NJAH1 and S T OGUNJO2 1 Department
of Physics, University of Agriculture, P.M.B. 2240, Abeokuta, Nigeria of Physics, Federal University of Technology, P.M.B. 704, Akure, Nigeria ∗ Corresponding author. E-mail:
[email protected] 2 Department
MS received 4 September 2012; revised 20 November 2012; accepted 6 December 2012 Abstract. In this article, projective synchronization of double–scroll attractor of an extended Bonöffer–van der Pol oscillator (BVPO) is considered via the backstepping and active control techniques. In each synchronization scheme, a single control function is designed to achieve projective synchronization between two Bonhöffer–van der Pol oscillator evolving from different initial conditions. To obtain a single control function via the active control, the coefficient of the error dynamics is chosen such that the number of control functions is reduced from three to one, thereby, reducing control function complexity in design. The results show that the transient error dynamics convergence and synchronization time are achieved faster via the backstepping than that of the active control technique. However, the control function obtained via the active control is simpler with a more stable synchronization time and hence, it is more suitable for practical implementation. Numerical simulations are presented to confirm the effectiveness of the analytical results. Keywords. Modified active control; backstepping control; projective synchronization of chaos in an extended Bonöffer–van der Pol oscillator. PACS Nos 05.45.Pq; 05.45.Gg; 05.45.Ac
1. Introduction Chaos synchronization is one the critical issues in nonlinear dynamical science because of its various applications in physics, secure communication, chemical reactors, control theory, biological networks, artificial neural networks, etc. Particularly, in neuroscience research, synchronization plays a very important role in the analysis of migraine and epilepsy [1]. In recent decades, the synchronization of chaotic systems has been extensively investigated both theoretically and practically [2–7]. The study of chaos synchronization has led to the discovery of various types of synchronization such as DOI: 10.1007/s12043-013-0526-3; ePublication: 12 April 2013
825
K S Ojo, A N Njah and S T Ogunjo complete synchronization [2], phase synchronization [3], lag synchronization [4], generalized synchronization [5], anticipated synchronization [6], projective synchronization [7], etc. Projective synchronization was reported by Mainnieri and Rehacek [7] in partially linear coupled three-dimensional systems in which the response of the dynamical state variables synchronized up to a scaling constant with the drive dynamical variables. Projective synchronization was later extended to general class of chaotic systems and not only to partially linear systems and this is called generalized projective synchronization [8,9]. Projective synchronization is one of the most interesting synchronization scheme because of the proportional relation that exists between the dynamical state variables of the drive and response systems, a feature that can be used to extend binary digit to M-nary digit for fast communication [10]. Following the seminal work of Pecora and Carroll [2], various synchronization techniques have been developed in search of improved and effective methods of achieving stable synchronous state between identical and nonidentical chaotic systems. These methods include feedback control, adaptive control, sliding mode, backstepping control, active control [11–19] to list a few. Notable among these methods are the active control and backstepping techniques which have excellent performance in the synchronization of identical and nonidentical chaotic systems [20,21]. The ubiquitous application of active control techniques has encouraged researchers to introduce active control based on different stability criteria. For instance, Lei et al [22] introduced active control based on Lyapunov stability theory and Routh Hurwitz criteria which has the advantage of possible implementation and has been used to synchronize a few chaotic systems [22,23]. However, the active control technique always gives controllers which are as many as the dimensions of the system been synchronized. In this article, active control technique is used to design a controller for projective synchronization of three-dimensional autonomous chaotic system in such a way that the number of controllers is reduced from three to one thereby, reducing significantly the controller complexity and cost. Hence, this makes it feasible for practical implementation. Meanwhile, the backstepping technique has been recognized as a powerful design technique for stabilization, tracking and synchronization of chaotic systems. It has been reported in [24] that backstepping can guarantee global stability, tracking and transient performance of a broad class of strict-feedback nonlinear systems. Recently, it has been used for controlling and tracking hyperchaotic systems [25]. According to ref. [26], some of the advantages of the method include: applicability to a variety of chaotic systems irrespective of whether they contain external excitation or not; needs only one controller to realize synchronization of chaotic systems and finally there is no derivative in the controller. Zhang et al [25] states that the controller is singularity-free from nonlinear term of quadratic forms, gives flexibility to construct a control law which can be extended to higher-dimensional hyperchaotic systems, while ref. [27] adds that it requires less control effort compared to differential geometric methods. The Bonhöffer–van der Pol (or Fizhugh–Nagumo) oscillator model was derived from van der Pol oscillator to give more accurate and reliable description of nonlinear dynamical systems which can show a stable threshold phenomenon as well as stable oscillations. The Bonhöffer–van der Pol oscillator is closely related to Fizhugh–Hulley (FH) model of the squid giant axon, to the cat’s carotid sinus nerve [28]. Indeed, many studies on Bonhöffer–van der Pol oscillator have shown its various applications in medicine (see, for 826
Pramana – J. Phys., Vol. 80, No. 5, May 2013
Backstepping and modified active control in projective synchronization example, refs [28–31] and references therein). Recently, there has been a resurgent interest in the theoretical and experimental generation of multiscroll chaotic attractors [32]. This is partly due to many practical applications foreseen in such fields as digital and secure communication, synchronous prediction, random number generation, information systems, etc. [31]. An efficient secure communication model should exhibit appreciable synchronous performance even in chaotic state. Thus, exploring projective synchronization of multiscroll chaotic system will be a stimulating subject of research. However, it has received inadequate attention. A few reports on synchronization of multiscroll chaotic systems are limited to complete synchronization [33–35]. The goal of this paper is to design a single control function via active control and backstepping techniques and compare their performance in projective synchronization of two identical double-scroll chaotic attractors generated from an extended Bonhöffer–van der Pol oscillator. To the best of our knowledge, this problem has not been considered till now. The rest of this paper is organized as follows: In §2, a brief description of the system of the extended Bonhöffer–van der Pol oscillator is given. Sections 3 and 4 are devoted to projective synchronization of extended Bonhöffer–van der Pol oscillators via active control and backstepping techniques with numerical simulations. Section 5 deals with the comparison of active control and backstepping techniques while, §6 concludes the paper. 2. Description of the Bonhöffer–van der Pol oscillator model The normalized equations for the extended Bonhöffer–van der Pol oscillator are x˙ = −z + Ax + tanh Bx, y˙ = z − δy z˙ = x − y
(1)
where x and y are the state variables corresponding to the voltage across the capacitors and z corresponds to the value of fixed resistor. A and B are the control parameters of the oscillator. Extensive experimental and numerical study of the bifurcation and chaotic phenomenon of the oscillator described in (1) was recently carried out by Nshiuchi et al [36]. Besides, various dynamical behaviours observed in double-scroll chaotic attractors of Bonhöffer–van der Pol oscillator were reported for (A, B, δ) = (1.0, 1.0, 1.2). This attractor is shown in figure 1. 3. Projective synchronization of chaos in an extended Bonhöffer–van der Pol oscillators via active control technique 3.1 Design of a single active controller Let the system of the Bonhöffer–van der Pol oscillator be written as x˙1 = −z 1 + Ax1 + tanh Bx1 , y˙1 = z 1 − δy1 , z˙ 1 = x1 − y1 . Pramana – J. Phys., Vol. 80, No. 5, May 2013
(2) 827
K S Ojo, A N Njah and S T Ogunjo 10 8 6 4
y
2 0 -2 -4 -6 -8 -10 -10
-8
-6
-4
-2
0
2
4
6
8
10
x
Figure 1. Two-dimensional view of the double-scroll chaotic attractor of the extended Bonhöffer–van der Pol oscillator with A = 1.0, B = 1.0 and δ = 1.2.
Then the response system is x˙2 = −z 2 + Ax1 + tanh Bx2 + u 1 (t), y˙2 = z 2 − δy2 + u 2 (t), z˙ 2 = x2 − y2 + u 3 (t),
(3)
where u i (t), i = 1, 2, 3 are the control functions to be determined. Subtracting (2) from (3) we obtain the error dynamics as e˙1 = −e3 + Ae1 + tanh Bx2 − α tanh Bx1 + u 1 (t), e˙2 = e3 − δe2 + u 2 (t), e˙3 = e1 − e2 + u 3 (t),
(4)
where e1 = x2 − αx1 , e2 = y2 − αy1 and e3 = z 2 − αz 1 . The projective synchronization between the drive (2) and the response (3) reduces to asymptotic stability of the error system (4) at equilibrium. To achieve this, control functions are re-defined to eliminate terms in (4) which cannot be expressed as a linear term in e1 , e2 and e3 as follows: u 1 (t) = α tanh Bx1 − tanh Bx2 + v1 (t) u 2 (t) = v2 (t) u 3 (t) = v2 (t). Substituting (5) into (4) yields ⎛ ⎞ ⎛ ⎞⎛ ⎞ ⎛ ⎞ e˙1 A 0 −1 e1 v1 (t) ⎝ e˙2 ⎠ = ⎝ 0 −δ 1 ⎠ ⎝ e2 ⎠ + ⎝ v2 (t) ⎠ . 1 −1 0 e˙3 e3 v3 (t)
(5)
(6)
When eq. (6) is stabilized by the feedback v1 (t), v2 (t) and v3 (t), the errors will converge to zero as t → ∞ which implies that projective synchronization of systems (2) and (3) is 828
Pramana – J. Phys., Vol. 80, No. 5, May 2013
Backstepping and modified active control in projective synchronization achieved. To achieve this goal using the active control technique, a constant matrix C is chosen to control the error dynamics (4) such that the feedback v1 (t), v2 (t) and v3 (t) are ⎛ ⎞ ⎛ ⎞ v1 (t) e1 ⎝ v2 (t) ⎠ = C ⎝ e2 ⎠ , (7) e3 v3 (t) where C is a 3 × 3 matrix. There are various choices of the feedback C that would control the error dynamics in (6). We optimize the way this choice is made so that the problem of controller complexity is significantly reduced and the choice is ⎛ ⎞ −(λ + A) 0 1 0 0 0⎠. C=⎝ (8) 0 0 0 With this choice the control functions are reduced from three to one which is u 1 (t) = α tanh Bx1 − tanh Bx2 − e1 (λ + A) + e3 ,
(9)
where λ is a positive constant control parameter. 3.2 Numerical simulation results
20
20
15
15
10
10
5
5
0
0
y
y
To verify the effectiveness of the designed controller we used the fourth-order Runge– Kutta algorithm with initial conditions (x1 , y1 , z 1 ) = (1.0, −1.0, 1.0) and (x2 , y2 , z 2 ) =
-5
-5
-10
-10
-15
-15
-20 -20 -15 -10 -5
(a)
0
x
5
10
15
20
-20 -20 -15 -10 -5
(b)
0
5
10
15
20
x
Figure 2. Projection of the attractor of the response Bonhöffer–van der Pol oscillator (green colour) on the drive Bonhöffer–van der Pol oscillator (red colour) with scaling factor of α = 2 via (a) active control; (b) the backstepping technique.
Pramana – J. Phys., Vol. 80, No. 5, May 2013
829
e1,e2,e3
K S Ojo, A N Njah and S T Ogunjo 20 15 10 5 0 -5 -10 -15 -20 -25 0
50
100
150
e1,e2,e3
200
250
300
350
200
250
300
350
time
(b) 20 15 10 5 0 -5 -10 -15 -20 -25 0
50
100
150
time
(a)
Figure 3. Error dynamics between two Bonhöffer–van der Pol oscillators with controller deactivated for 0 < t < 100 and activated for t ≥ 100 for a scaling factor α = 2 via (a) active control; (b) backstepping technique.
(0.5, −0.5, 0.5), a time step of 0.005, and fixing the parameter values of the system such that the systems exhibits double-scroll chaotic attractors as in figure 1 to ensure chaotic dynamics of the state variables. We solved systems (2) and (3) with the control function defined in (9) for λ = 1. The result displayed in figure 2a shows the projection of the response attractor on the drive system for α = 2.0. The two systems achieved projective synchronization as indicated by the convergence of the error state variables as soon as the controller is switched on for t ≥ 100 (figure 3a). The results obtained confirm the effectiveness of the single control function designed for the synchronization of a third-order double-scroll chaotic Bonhöffer–van der Pol oscillators via the active control technique.
4. Projective synchronization of chaos in an extended Bonhöffer–van der Pol oscillator via backstepping technique 4.1 Design of a single active controller Let the drive system of the Bonhöffer–van der Pol oscillator be as in (2) and the response Bonhöffer–van der Pol oscillator as x˙2 = −z 2 + Ax2 + tanh Bx2 + u(t), y˙2 = z 2 − δy2 , z˙ 2 = x2 − y2 , 830
Pramana – J. Phys., Vol. 80, No. 5, May 2013
(10)
Backstepping and modified active control in projective synchronization where u(t) is the control function to be determined. Subtracting (2) from (10) and rearrange the error dynamics, the error system can be written as e˙2 = e3 − δe2 , e˙3 = e1 − e2 , e˙1 = −e3 + Ae1 + tanh Bx2 − α tanh Bx1 + u(t),
(11)
where e1 = x2 − αx1 , e2 = y2 − αy1 and e3 = z 2 − αz 1 . The objective of this paper is to find a control function u(t) that can stabilize the error state in (11) at the origin. Firstly, we stabilize the first equation in (11) by regarding e3 as a controller, choosing a Lyapunov function V1 (e2 ) = 12 e22 and differentiating it with respect to time we have V˙1 = e2 e3 − δe22 .
(12) ˙ We estimate that the controller e3 = α1 (e2 ). Then, eq. (12) can be written as V1 = −δe22 + e2 α1 (e2 ). V˙1 is negative definite if the estimated function α1 (e2 ) = 0. The error W2 between e3 and α1 (e2 ) is W2 = e3 − α1 (e2 ).
(13)
Differentiating (13) with respect to time yields W˙ 2 = e1 − e2 .
(14)
regarding e1 as the controller Choose a Lyapunov function V2 (e2 , W2 ) = V1 (e3 ) + in (14) and assume that when e1 = α2 (e2 , W2 ) then, the time derivative of the Lyapunov function is (15) V˙2 = −δe22 + W2 α2 (e2 , W2 ) 1 W 2, 2 2
if α2 (e2 , W2 ) = −kW2 then V˙2 = −δe22 − kW22 < 0 is negative definite. Therefore, system (14) is asymptotically stable. The error W3 between e1 and α2 (e2 , W2 ) is W3 = e1 − α2 (e2 , W2 ) = e1 + kW2
(16)
e1 = W3 − kW2 .
(17)
or Substitution of (17) into the time derivative of (16) yields W˙ 3 = −e3 + Ae1 + tanh Bx2 − α tanh Bx1 + k(e1 − e2 ) + u(t).
(18)
and obtain its time We choose Lyapunov function V3 (e3 , W2 , W3 ) = V2 (e2 , W2 ) + derivative as V˙3 = −δe22 − kW22 + W3 [−e3 + (A + k)(W3 − W2 ) − ke2 + tanh Bx2 1 W2 2 3
− α tanh Bx1 + u(t)].
(19)
The control function is chosen as u(t) = e3 − (A + k)(W3 − kW2 ) − k(W3 − e2 ) − tanh Bx2 + α tanh Bx1
(20)
so that V˙3 = −δe22 − kW22 − kW32 < 0 which is negative definite, where k is a positive constant control parameter. It follows that all solutions of (11) converge to the manifold e1 = e2 = e3 = 0 as t → ∞ and hence systems (2) and (10) are globally synchronized. Thus, projective synchronization is achieved. Pramana – J. Phys., Vol. 80, No. 5, May 2013
831
K S Ojo, A N Njah and S T Ogunjo 4.2 Numerical simulation results To verify the effectiveness of the designed controller we used the fourth-order Runge– Kutta algorithm with initial conditions (x1 , y1 , z 1 ) = (1.0, −1.0, 1.0) and (x2 , y2 , z 2 ) = (0.5, −0.5, 0.5), with a time step of 0.005, fixing the parameter values of the system such that the system exhibits double-scroll chaotic attractors as in figure 1 to ensure chaotic dynamics of the state variables. We solve systems (2) and (10) with the control function defined in (20) for k = 1. The result is displayed in figure 2b which shows the projection of the response attractor on the drive system for α = 2. The two systems achieved projective synchronization as indicated by the convergence of the error state variables as soon as the controller is switched on for t ≥ 100 (figure 3b). The results obtained confirm the effectiveness of the single control function designed for the synchronization of a third-order double-scroll chaotic Bonhöffer–van der Pol oscillators via the backstepping technique.
5. Comparison of active control and backstepping techniques In order to make a detailed comparison between the two technique we obtained a relation between synchronization time and the control parameters. In our simulation, we varied the value of the control parameter k and λ between 0.1 and 10, and the result obtained
Synchronization time
180 160 140 120 100
0
2
4
(a)
6
8
10
6
8
10
k
Synchronization time
180 160 140 120 100
(b)
0
2
4
λ
Figure 4. Comparison of synchronization time between active control and backstepping technique with controller deactivated for 0 < t < 100 and activated for t ≥ 100 for a scaling factor α = 2 via (a) backstepping technique; (b) active control.
832
Pramana – J. Phys., Vol. 80, No. 5, May 2013
Backstepping and modified active control in projective synchronization for backstepping technique is shown in figure 4a, while, that for active control is shown in figure 4b. For clarity and effective comparison, the results in figures 4a and 4b are combined in figure 5. From the figures we notice that for active control, the least synchronization time is 111.8 which occurs at λ = 0.6 while, for backstepping technique the least synchronization time is 105.1 which occurs at k = 2.2 and 2.5. In order to obtain the fastest synchronization time (111.8) using the active control the choice of λ should be 0.6 while, for the backstepping technique k should be 2.2 or 2.5 in order to obtain fastest synchronization time (105.1). Moreover, it clear from figure 5 that the active control has a more stable synchronization time. Increase in the control parameters beyond the values quoted above amount to waste of energy. There is no significant difference in the time for the onset of synchronization. We also compared the result using the rate of error convergence e and the rate at which the ratio r/d tends to the predefined scaling factor (figures 6a and 6b) respectively, where e=
e12 + e22 + e32 ,
r=
x22 + y22 + z 22 ,
and d=
x12 + y12 + z 12 .
The results also show that the transient error dynamics via the backstepping technique converges faster than that of the active control technique.
180 Active control Backstepping
170
Synchronization time
160 150 140 130 120 110 100 0
2
4
6
8
10
k and λ
Figure 5. Synchronization time against the variation of the positive constant parameters with controller deactivated for 0 < t < 100 and activated for t ≥ 100 for a scaling factor α = 2 via (a) backstepping technique; (b) active control.
Pramana – J. Phys., Vol. 80, No. 5, May 2013
833
K S Ojo, A N Njah and S T Ogunjo
r/d
2.2 2 1.8 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0 90
95
100
(b)
105
110
115
105
110
115
time 30 25
e
20 15 10 5 0 90
(a)
95
100
time
Figure 6. (a) The time evolution of transient error dynamics; (b) the variation of the scaling factor via active control (red) and backstepping (green).
6. Conclusion We have investigated projective synchronization of double-scroll chaotic attractors of an extended Bonhöffer–van der Pol oscillator via backstepping and active control technique. In each synchronization scheme, we designed a single control function which achieved projective synchronization between two identical Bonhöffer–van der Pol oscillators evolving from different initial conditions. For the active control technique, the choice of the coefficient matrix of the error dynamics was chosen such that the number of control functions reduced from three to one, thereby, reducing the controller complexity in the design. The results show that the transient error dynamics convergence and synchronization time are achieved faster via the backstepping than via the active control technique. However, the control function obtained via the active control is simpler with a more stable synchronization time and hence, it is more suitable for practical implementation. Numerical simulations are presented to confirm the effectiveness of the analytical results. References [1] Z-Y Liu, C-J Liu, M-C Ho, Y-C Haung, T-F Hsu and I-M Jiang, Int. J. Bifurcat. Chaos 18(2), 3731 (2008) [2] L M Pecora and T L Carrol, Phys. Rev. Lett. 64, 821 (1990) [3] Y Li, J Zhang and Z Li, Int. J. Nonlinear Sci. 1(3), 131 (2006) [4] S Zhu and L Wu, Int. J. Mod. Phys. B 18, 2547 (2004) [5] Z Zhang and G Hu, Phys. Rev. E 62, 7882 (2000) [6] O Calvo, D R Chialvo, V M Eguituz, C R Mirasso and R Toral, Chaos 14, 7 (2004)
834
Pramana – J. Phys., Vol. 80, No. 5, May 2013
Backstepping and modified active control in projective synchronization [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36]
R Mainieri and J Rehacek, Phys. Rev. Lett. 82, 3042 (1999) G-H Li, Chaos, Solitons and Fractals 30, 71 (2005) H Chen and M Chen, Int. J. Nonlinear Sci. 2(3), 166 (2006) H-L Zhu and X-B Zhang, J. Information and Computing Sci. 4(1), 33 (2009) M Sun, L Tian, S Jian and J Xu, Chaos, Solitons and Fractals 32, 1725 (2007) G Cai, J Huang, L Tian and Q Wang, Int. J. Nonlinear Sci. 2(1), 1725 (2006) M Jiang and C-K Chen, Int. J. Bifurcat. Chaos 12, 1437 (2002) A N Njah and K S Ojo, Far East J. Dynamical Systems 11(2), 143 (2009) A N Njah and U E Vincent, J. Sound and Vibrations 319, 41 (2009) Y Lei, K-L Yung and Y Xu, J. Sound and Vibration 329, 973 (2010) Y Lei, W Xu, J Shen and T Fang, Chaos, Solitons and Fractals 28, 428 (2006) Z Zhang, Y Wang and Z Du, Appl. Math. Comput. 218, 6833 (2012) Z Zhang, H Shen and J Li, Appl. Math. Comput. 218, 4260 (2011) E W Bai and K E Lonngren, Chaos, Solitons and Fractals 8, 51 (1997) M C Ho and Y C Hung, Phys. Lett. 301, 424 (2002) Y Lie, W Xu, J Shen and F Fang, Chaos, Solitons and Fractals 28, 428 (2006) A N Njah, J. Sound and Vibrations 327, 322 (2009) M Krstic, I Kanellakopoulus and P O Kokotovic, Nonlinear and adaptive control design (John Wiley, New York, 1995) H Zheng, X Ma, M Li and J Zou, Chaos, Solitons and Fractals 26, 353 (2005) X Tian, J Zhang and Y Yang, Chaos, Solitons and Fractals 16, 37 (2003) S Mascolo, Backstepping design for design Lorenz chaos, Procceding 36th IEEE CDC (San Diego, CA) pp. 1500–1501 D Paydafar, F L Eldridge and J A Paydafar, The J. Physiol. 506, 515 (1998) R Fitzhugh, Biophys. J. 1, 445 (1961) S Rajasekar and M Lakshmanan, Physica D 67, 146 (1993) M Ramesh and S Nayaranan, Chaos, Solitons and Fractals 12, 2395 (2001) J Lu and J Chen, Int. J. Bifurcat. Chaos 16, 775 (2006) M E Yalcin, Cellular neural networks, multi-scroll chaos and synchronization: Theory, application and implication, Ph.D. thesis (Katholieke University Leuven, Belgium, 2004) M E Yalcin, J A Suykens and J Vandewalle, Cellular neural networks, multi-scroll chaos and synchronization (World Scientific, Singapore, 2005) S M Yu, Z G Ma, S S Qiu, S G Peng and Q H Lin, Chin. Phys. 13, 317 (2004) Y Nishiuchi, T Ueta and H Kawakami, Chaos, Solitons and Fractals 27, 941 (2006)
Pramana – J. Phys., Vol. 80, No. 5, May 2013
835
