Journal of Chaos

Volume 2013 (2013), Article ID 839038, 12 pages

http://dx.doi.org/10.1155/2013/839038

Research Article

Dynamical Properties and Finite-Time Hybrid Projective Synchronization Using Fractional Nonsingular Sliding Mode Surface in Fractional-Order Two-Stage Colpitts Oscillators

Romanic Kengne,1,2 Robert Tchitnga,1,2 Arnaud Nzeusseu Tchikankou,1,2 Aurelle Tchagna Kouanou,¹ and Anaclet Fomethe³

¹Laboratory of Electronics and of Signal Processing, Department of Physics, Faculty of Science, University of Dschang, P.O. Box 67, Dschang, Cameroon

²Research Group on Experimental and Applied Physics for Sustainable Development (EAPhySuD), P.O. Box 412, Dschang, Cameroon

³Laboratoire de Mécanique et de Modélisation des Systèmes, Département de Physique, Faculté des Sciences, Université de Dschang, B.P. 67, Dschang, Cameroon

Received 30 August 2013; Accepted 28 October 2013

Academic Editor: Uchechukwu E. Vincent

Copyright © 2013 Romanic Kengne et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The dynamics and robust finite-time hybrid projective synchronization of a fractional-order four-dimensional nonlinear system based on a two-stage Colpitts oscillator is investigated. The study of the fractional order stability of the equilibrium states of the system is carried out. The bifurcation diagram confirms the occurrence of Hopf bifurcation in the proposed system when the fractional-order passes a sequence of critical values; the Lyapunov exponent shows the different chaotic sequences of the system. Further, a fractional nonsingular terminal sliding surface and an appropriate robust fractional sliding mode control law are proposed for the finite-time hybrid projective synchronization of a fractional-order chaotic two-stage Colpitts oscillator by taking into account the effects of model uncertainties and the external disturbances. The fractional version of the Lyapunov stability is used to prove the finite-time existence of the sliding motion. Finally, some numerical simulations are presented to demonstrate the effectiveness and applicability of the proposed technique.

1. Introduction

Fractional calculus has an about 300-year-old history, but its applications to physics and engineering are rather recent [1]. Many systems are known to display fractional-order dynamics, such as viscoelastic systems, dielectric polarization, and electromagnetic waves [2–4], just to name some.

For some decades, there is a growing interest in investigating the chaotic behavior and dynamics of fractional-order dynamic systems; this can be understood as it has been found that fractional-order systems possess memory and display more sophisticated dynamics compared to their integral-order counterparts, something that is of great significance in secure communication [5–20]. It has been shown that several chaotic systems can remain chaotic when their models become fractional [5]. A three-dimensional fractional-order modified hybrid optical system is presented in [10] where it was shown that Hopf bifurcation occurs on the proposed system when the fractional order varies and passes a sequence of critical values. Despite these many examples the bifurcation of fractional-order nonlinear system has been studied using solely the Caputo derivative definition and limited to three-dimensional systems.

On the other hand, in the past two decades, a new direction of chaos research has emerged to address the more challenging problem of chaos synchronization due to its potential applications in laser physics, chemical reactions, secure communication, biomedicine, and so on [21–23]. The thrust of research within this area aims at achieving a master-slave synchronization between two chaotic systems by choosing various kinds of methods following the pioneering work of Pecora and Carroll [24]. In [11], based on the stability theory, a novel fractional-order controller is presented for the control and synchronization of the fractional-order Lorenz chaotic system via the fractional derivative. In [12], based on tracking control and the stability theory of nonlinear fractional-order systems, Zhou et al. propose the multidrive-one response synchronization technique which is simple and theoretically rigorous. In [13], the hybrid projective synchronization of different dimensional chaotic fractional-order systems is investigated based on the stability theory of linear fractional-order systems. Based on the fractional-order stability theory and control tracking, Zhou and Zhu investigate the function projective synchronization for the fractional-order chaotic systems [7]. In [25], based on the stability theory of the fractional-order system, the authors studies on projective synchronization of the new fractional-order chaotic system through designing the suitable nonlinear controller are investigated. All these examples perfectly clarify the importance of the fractional-order stability theory, but these works have a common drawback from a practical point of view: the lack of knowledge of analytical time of synchronization.

In recent years some researchers have applied finite-time control and synchronization of fractional-order systems by means of the fractional-order Lyapunov stability theory. In [26], the authors design a novel fractional-order nonsingular sliding mode controller for robust synchronization problem of a class of fractional-order chaotic systems in the presence of model uncertainties and external disturbances, the stability of a novel fractional-order integral type sliding surface is proven, based on fractional-order Lyapunov stability theory, a robust sliding control law is derived to guarantee the occurrence of the sliding motion in finite time. In [27], using a nonsingular sliding mode surface, the authors study the finite-time synchronization problem of fractional-order chaotic/hyperchaotic systems in the presence of both model uncertainties and external disturbances. These last works will serve as our handbook in this paper.

It is important to seek bifurcations such as those of Hopf in systems because they are routes towards chaos in such systems. Concerning the system with fractional-order, very little work emphasizes the conditions leading to Hopf bifurcation. Moreover modified projective synchronization is a general case of the simple synchronization. But it had been shown that the finite-time synchronization has very great applications on a practical point of view. Finally, to the best of our knowledge there are few works interested in fractional-order hybrid projective synchronization in finite time and using a fractional nonsingular terminal sliding mode surface; the importance of this surface as that used in [27] is that it supports the robustness of the synchronization in the presence of the disturbances and uncertainties, and that it is easily stabilisable with zero and in finite time.

Motivated by the above discussion, at first in the present work, we propose to tackle the problem of bifurcation of a four-dimensional fractional-order nonlinear system. The two-stage well-studied Colpitts oscillator presented in [28] is a good candidate for the study, due to its broad band in frequency domain. Secondly we study the finite-time hybrid projective synchronization problem of fractional-order two-stage Colpitts oscillator in the presence of both model uncertainties and external disturbances. The modified nonsingular terminal sliding mode surface is introduced; its finite-time stability to zero is proved via the Lyapunov stability. So on basis of fractional-order Lyapunov stability theory, a robust control law is designed to force the trajectories of the synchronization error system onto the sliding surface within a finite time and remain on it forever. Numerical simulations demonstrate the applicability and efficiency of the nonlinear control law and verify the theoretical results of the paper.

The rest of this paper is organized as follows. In Section 2, the fractional-order system developed around a two-stage Colpitts oscillator is proposed and its dynamics studied. The numerical results of the dynamics are presented and discussed in Section 3, while the next section is devoted to the synchronization of two two-stage Colpitts oscillators. Finally, Section 5 concludes this work.

2. The Fractional-Order of a Two-Stage Colpitts Oscillator

2.1. Dynamics of the System

The proposed four-dimensional fractional-order system under study described by the set of (1) is obtained by modifying the integer-order two-stage Colpitts oscillator proposed in [28]:

Here, the parameters , , , and are positive reals, , and is the fractional order. According to [28], when , the system (1) exhibits chaotic behavior with the parameter values , , , and .

2.1.1. Stability Analysis

Theorem 1 (see [15, 16]). The following commensurate order system

with , , and , is asymptotically stable if and only if is satisfied for all eigenvalues of the matrix . Moreover, this system is stable if and only if is satisfied for all eigenvalues of with those critical eigenvalues that satisfy having geometric multiplicity of one.

Theorem 2 (see [17]). Consider the following linear fractional-order system:

where , , and , with and , .

Let be the lowest common multiple of the denominators ’s. The zero solution of system (3) is globally asymptotically stable in the Lyapunov sense if all roots ’s of the equation

satisfy .

2.1.2. Stability of the Equilibrium Points

In this section we proceed with commensurate order . Fractional order of the proposed two-stage Colpitts oscillator (1), when , has one equilibrium point, . The Jacobian matrix of system (1), evaluated at the equilibrium point, is

For around and to the equilibrium point , the equation with and becomes

Thus, for , and for , .

Therefore, on the basis of Theorem 2, system (1) is asymptotically stable at equilibrium point for and unstable for .

2.2. Hopf Bifurcation

One of the basic differences between the dynamical behavior of fractional-order systems and that of integer-order systems is that the limit set of a trajectory of integer-order system such as a limit cycle is solution for the system under consideration, while in the case of fractional-order systems, such a limit set of a trajectory may not be solution for this system [18]. In [19], the authors claimed that there are no periodic orbits in fractional-order systems, and in [20], an example is given where the solutions of the system are also not periodic but do converge to periodic signals, confirming in both cases what has been stipulated in [18].

In the present paper, we consider the final state of trajectory that appears at the Hopf bifurcation (after suppression of the transitory state). It is also not a periodic solution of the fractional-order system (1) but attracts nearby solutions.

Let us consider the following four-dimensional fractional-order commensurate system:

where , , and suppose that is an equilibrium point of this system. In the integer case , the stability of is related to the sign of , , where are the eigenvalues of the Jacobian matrix . If for all , then is locally asymptotically stable. If there exists an for which , then is unstable.

To undergo a Hopf bifurcation at the equilibrium point when , system (7) with must fulfill the following conditions:

(i)The Jacobian matrix must have two pairs of complex-conjugate eigenvalues and .(ii), with ;(iii), with ; and finally(iv).

In the fractional case, the stability of is related to the sign of , with . If for all , the is locally asymptotically stable. If there exists any for which , then, the equilibrium point is unstable. So, the function for fractional-order systems has a similar effect as the real part of eigenvalues in integer system. Therefore, we can extend the Hopf bifurcation condition to the fractional systems by replacing with as follows, compared with [10]:

(i),(ii).

2.3. Hopf Bifurcation versus the Parameter and the Fractional Order

In this subsection, we consider the parameter values for the search for the Hopf bifurcation around the equilibrium point .

Figure 1(a) depicts the solution of equation , while the black curve on Figure 1(b) recalls that for all , and for the blue curve in Figure 1(b) it can be noted that for all except for . We have ; thus, the proposed fractional-order Hopf bifurcation conditions are verified for all pair solutions of  , except for .

(a)

(b)

Figure 1: Critical values versus the fractional order . This curve depicts the couples of values for which the Hopf bifurcation occurs in the system.

3. Numerical Results

For numerical calculation of fractional-order derivatives, methods defined with (8) to (10) are usually used. The Grünwald-Letnikov (GL) method [14] is given in the following equation:

where indicates the integer part.

The Riemann-Liouville (RL) definition follows as

where is the gamma function. The Caputo definition of fractional derivatives can also be recalled as

Based on the fact that for a wide class of functions the three definitions—GL (8), RL (9), and Caputo’s (10)—are equivalent if , we can then use relation (11) derived from the GL definition (8). The new relation for the explicit numerical approximation of th derivative at the points has the following form:

where is the “memory length,” , with the time step of calculation and the binomial coefficients. For their calculation we can use for instance the following expression:

The binomial coefficients can also be expressed using a factorial. The gamma function can allow the generalization of the binomial coefficient to noninteger argument. Thus, relation (12) can be rewritten as follows:

3.1. Bifurcation and Chaos versus the Parameter

In this subsection, the dynamical behavior of system (1) is numerically investigated by means of bifurcation diagram and largest Lyapunov exponents, which measure the exponential rates of divergence or convergence of nearby trajectories in phase space. For taken as control parameter and the following other parameter values, fractional order , , and , the critical Hopf bifurcation value is localized at (see Figure 2(a)) and confirmed by the largest Lyapunov exponents presented in Figure 2(b). When , the equilibrium point is a locally asymptotically stable focus; the neighbors trajectories converge to . This is supported by the negative sign of the largest Lyapunov exponents. For , system (1) undergoes a Hopf bifurcation as mentioned above. The fixed point becomes unstable, and a period-one limit cycle appears. A period-two limit cycle follows for , leading to a new bifurcation at , as the system undergoes a period-four bifurcation. This bifurcations scenario continues through a period-height limit cycle for up to a critical value of corresponding to the appearance of a chaotic attractor. This chaotic behavior is confirmed by the existence of positive largest Lyapunov exponents. Figure 3 depicts the phase portraits presenting routes to chaos according to the abovementioned parameter values.

(a)

(b)

Figure 2: Bifurcation diagram expressing the (a) dynamics of the system variable and (b) largest Lyapunov exponent, both as a function of , with .

(a)

(b)

(c)

(d)

Figure 3: Phase portrait of system (1) for different values of , with ; (a) period-1 for , (b) period-2 for , (c) period-4 for , and (d) chaos for .

3.2. Bifurcation and Chaos versus the Fractional Order

The fractional order is taken as control parameter, while is fixed at . The critical Hopf bifurcation value is localized at , using the previously proposed conditions. The resulting bifurcation diagram (Figure 4(a)) for the second variable of the set of (1) is plotted as a function of the fractional order and corresponding largest Lyapunov exponents in Figure 4(b).

(a)

(b)

Figure 4: Bifurcation diagram expressing the (a) dynamics of the system variable and (b) largest Lyapunov exponent, both as a function of , with .

When , the equilibrium point is a locally asymptotically stable focus confirmed by the negative sign of the largest Lyapunov exponents; the neighbors trajectories converge to this origin. For , system (1) undergoes a Hopf bifurcation as mentioned above. The fixed point becomes unstable, and a period-1 limit cycle appears for . As the fractional order parameter nears the value , a new bifurcation occurs for period-2 limit cycle. This is followed by a period-4 limit cycle at . This bifurcation scenario continues up to a critical value where a chaotic attractor appears, sustained by the existence of positive largest Lyapunov exponents. For a periodic steady state, all spikes in the power spectrum are harmonically related to the fundamental, whereas a broadband noise like power spectrum is associated with a chaotic steady state. The periodicity of the attractor (i.e., total number of frequencies in a wave) is deduced by counting the number of spikes located at the left-hand side of the highest spike (the latter is included). Indeed, we have obtained the complete scenarios to chaos presented in Figure 5. Specifically, the following scenario was observed when monitoring the control parameter: fixed point behavior period-1 period-2 period-4 chaos.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5: Phase portrait of system (1) for different values of , with and corresponding power spectra: ((a) and (e)) period-1 for , ((b) and (f)) period-2 for , ((c) and (g)) period-4 for , and ((d) and (h)) chaos for .

4. Finite-Time Hybrid Projective Synchronization of Two Fractional-Order Two-Stage Colpitts Oscillators

4.1. Analytic Results

This section is devoted to the finite-time hybrid projective synchronization of the drive and response commensurate fractional order of a two-stage Colpitts system using a robust fractional nonsingular terminal sliding mode controller, for . The drive system is defined as follows:

The drive system (14) can be written as well as

Accordingly, the response system takes the following form:

where represents the nonlinear controllers and and represent unknown model uncertainty and external disturbances of the system. By subtracting (15) from (16) and setting

the following set of equations defining the errors is obtained:

Theorem 3 (see [27]). Let be an equilibrium point for the nonautonomous fractional-order system , where satisfies the Lipschitz condition with Lipschitz constant and . Assume that there exists a Lyapunov function satisfying

where , , , and are positive constants. Then the equilibrium point of the system is Mittag-Leffler stable.

Assumption 4. The uncertainty terms and external disturbance are bounded by and , , where and are positive constants.

Lemma 5. Assume , , and are real numbers; then the following inequality holds: .

4.1.1. Main Results

In this subsection a modified nonsingular terminal sliding surface is proposed as

where are the sliding surface parameters to be introduced later and are the scaling factor to content in synchronization error. The nonsingular terminal sliding surface as defined in (20) present a major advantage on that proposed in [27], namely its dependence on scaling factor no matter the initial conditions which, allows to improve the synchronization time between the driven system and response system. For the existence of the sliding mode it is necessary and sufficient that and [27]. Therefore, the dynamics of the proposed nonsingular terminal sliding mode can be obtained as

Theorem 6. The system (20) is finite time stable and its trajectories converge to the equilibrium in a finite-time, , determined by

where and .

Proof. Consider the candidate following positive definite Lyapunov function:

The time derivative of the Lyapunov function along the trajectories of (21) is

Using Lemma 5, one can obtain

Hence, according to Theorem 3, the error system (21) will converge to zero asymptotically. In order to show that the sliding motion occurs in finite time, we can obtain the convergence time as follows.

From inequality (25) we have

What leads us to

Taking integral of both sides of (27) from 0 to and letting , we have

Therefore, the state trajectories of the error system (21) will converge to in the finite time . This completes the proof.

A control law which forces the error trajectories to go onto the sliding surface within a finite time and remain on it forever is designed as follows:

, where are the sliding surface parameters to be introduced later, are the scaling factor to content in synchronization error, and and are positive constants.

Theorem 7. If the error system (18) is controlled with control law (29), then the states of the system will move toward the sliding surface and will approach the sliding surface in a finite time, , given by

Proof. Choosing a Lyapunov function in the form of

and taking the time derivative, one has

Inserting (20) into (32), we have

Inserting (18) into (33), we have

Using expression into (34), we obtain

On the basis of Assumption 4, we get

what leads to

Introducing (29) into (37), we have

Using Lemma 5, one can obtain

where .

Hence, according to Theorem 3, the error system (18) will converge to asymptotically.

In order to show that the sliding motion occurs in finite time, we can obtain the convergence time as follows.

From inequality (39), we have

what leads us to

Taking integral of both sides of (41) from 0 to and letting , we have

Remark 8. According to Theorems 6 and 7, the sliding mode control law (29) and the sliding surface (20) can make the response system (16) reach the drive system (15) in the finite-time .

4.2. Numerical Results

The drive and response system (14) and (16) are numerically integrated with the parameter values , , , , and ; the uncertainties and external disturbance in (16) are selected by

Control laws (29) use the nonsingular surface defined at (20) and the rest of functions of (29) are chosen as:

and the parameters of (29) are selected as follows: , , , , and . The initial conditions are given as for the driver and for the response. The analytics time of synchronization is , whereas in Figure 6 we show that the numerical time of synchronization is . These results confirm our theoretical result. Figure 7 displays for the parameters selects as follows: , , , and are given in the key. The initial conditions are as follows: for the master and for the slave. Figure 7 shows us that when is far from one, the synchronization time turns fast to zero; this decreasing numerical time of synchronization justifies the decreasing analytical time of synchronization which is, respectively, , , , and for , , , and .

(a)

(b)

(c)

Figure 6: (a) Time evolution of the sliding mode surface, showing the (b) synchronization process between the two coupled chaotic oscillators starting from different initial condition and (c) synchronization errors obtained for , , , , , , , , , and .

Figure 7: Time evolution of .

5. Conclusion

In this paper the dynamics and synchronization of a proposed four-dimensional fractional-order two-stage Colpitts oscillator have been investigated using analytical and numerical methods. The analytic method proved the existence of the Hopf bifurcation as well as the beach of the control parameter for which the system is stable. On the basis of fractional Lyapunov stability theory we determined with success the conditions under which the synchronization of two systems is achieved. For numerical simulation we used the Grünwald-Letnikov method, the largest Lyapunov exponents, and the bifurcation diagrams to show the period-doubling bifurcation routes to chaos as well as the Hopf bifurcation. The numerical analysis validates the conditions of Hopf bifurcation. For the finite-time hybrid projective synchronization the numerical investigation validates also the analytic conditions which achieve synchronization. Numerical simulations have been used to show the effectiveness of the proposed synchronization techniques.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, New Jersey, NJ, USA, 2001.

R. C. Koeller, “Application of fractional calculus to the theory of viscoelasticity,” Journal of Applied Mechanics, vol. 51, no. 2, pp. 299–307, 1984.

H. H. Sun, A. A. Abdelwahab, and B. Onaral, “Linear approximation of transfer function with a pole of fractional order,” IEEE Transactions on Automatic Control, vol. 29, no. 5, pp. 441–444, 1984.

O. Heaviside, Electromagnetic Theory, Chelsea, New York, NY, USA, 1971.

I. Grigorenko and E. Grigorenko, “Chaotic dynamics of the fractional Lorenz system,” Physical Review Letters, vol. 91, no. 3, Article ID 034101, 2003.

K. Sun and J. C. Sprott, “Bifurcations of fractional-order diffusionless lorenz system,” Electronic Journal of Theoretical Physics, vol. 6, no. 22, pp. 123–134, 2009.

P. Zhou and W. Zhu, “Function projective synchronization for fractional-order chaotic systems,” Nonlinear Analysis, vol. 12, no. 2, pp. 811–816, 2011.

D. Chen, Y. Q. Chen, and H. Sheng, “Fractional variational optical flow model for motion estimation,” in Proceedings of the 4th IFAC Workshop Fractional Differentiation and Its Application, pp. 18–20, Badajoz, Spain, October 2010.

J. C. Gutiérrez-Vega and C. López-Mariscal, “Nondiffracting vortex beams with continuous orbital angular momentum order dependence,” Journal of Optics A, vol. 10, no. 1, Article ID 015009, pp. 1–8, 2008.

M.-S. Abdelouahab, N.-E. Hamri, and J. Wang, “Hopf bifurcation and chaos in fractional-order modified hybrid optical system,” Nonlinear Dynamics, vol. 69, no. 1-2, pp. 275–284, 2012.

P. Zhou, “Control and synchronization of the fractional-order Lorenz chaotic system via fractional-order derivative,” Mathematical Problems in Engineering, vol. 2012, Article ID 214169, 14 pages, 2012.

P. Zhou, R. Ding, and Y. X. Cao, “Multi drive-one response synchronzation for fractional-order chaotic systems,” Nonlinear Dynamics, vol. 70, no. 2, pp. 1263–1271, 2012.

S. Wang, Y. Yu, and M. Diao, “Hybrid projective synchronization of chaotic fractional order systems with different dimensions,” Physica A, vol. 389, no. 21, pp. 4981–4988, 2010.

S. H. Hosseinnia, R. Ghaderi, A. Ranjbar N., M. Mahmoudian, and S. Momani, “Sliding mode synchronization of an uncertain fractional order chaotic system,” Computers and Mathematics with Applications, vol. 59, no. 5, pp. 1637–1643, 2010.

E. Ahmed, A. M. A. El-Sayed, and H. A. A. El-Saka, “Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models,” Journal of Mathematical Analysis and Applications, vol. 325, no. 1, pp. 542–553, 2007.

D. Matignon, “Stability results for fractional differential equations with applications to control processing,” IEEE-SMC Computational Engineering in Systems Applications, vol. 2, pp. 963–968, 1996.

A. Razminia, V. J. Majd, and D. Baleanu, “Chaotic incommensurate fractional order Rössler system: active control and synchronization,” Advances in Difference Equations, vol. 2011, p. 15, 2011.

M. S. Tavazoei, M. Haeri, M. Attari, S. Bolouki, and M. Siami, “More details on analysis of fractional-order Van der Pol oscillator,” JVC/Journal of Vibration and Control, vol. 15, no. 6, pp. 803–819, 2009.

M. S. Tavazoei and M. Haeri, “A proof for non existence of periodic solutions in time invariant fractional order systems,” Automatica, vol. 45, no. 8, pp. 1886–1890, 2009.

M. S. Tavazoei, “A note on fractional-order derivatives of periodic functions,” Automatica, vol. 46, no. 5, pp. 945–948, 2010.

L. O. Chua and M. Itah, “Chaos synchronization in Chua's circuits,” Journal of Circuits, Systems and Computers, vol. 3, no. 1, pp. 93–108, 1993.

L. O. Chua, T. Yang, and G. Q. Zhong, “Adaptive synchronization of Chua's oscillator,” International Journal of Bifurcation and Chaos, vol. 6, no. 1, pp. 189–201, 1996.

G. R. Chen and X. Dong, From Chaos to Order, World Scientific, Singapore, 1998.

L. M. Pecora and T. L. Carroll, “Synchronization in chaotic systems,” Physical Review Letters, vol. 64, no. 8, pp. 821–824, 1990.

X. Wu and H. Wang, “A new chaotic system with fractional order and its projective synchronization,” Nonlinear Dynamics, vol. 61, no. 3, pp. 407–417, 2010.

M. P. Aghababa, “Robust stabilization and synchronization of a class of fractional-order chaotic systems via a novel fractional sliding mode controller,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 6, pp. 2670–2681, 2012.

M. P. Aghababa, “finite-time chaos control and synchronization of fractional-order nonautonomous chaotic (hyperchaotic) systems using fractional non-singular terminal sliding mode technique,” Nonlinear Dynamics, vol. 69, no. 1-2, pp. 247–261, 2012.

J. Kengne, J. C. Chedjou, V. A. Fono, and K. Kyamakya, “On the analysis of bipolar transistor based chaotic circuits: case of a two-stage colpitts oscillator,” Nonlinear Dynamics, vol. 67, no. 2, pp. 1247–1260, 2012.