Table of Contents
Journal of Nonlinear Dynamics
Volume 2017, Article ID 5483956, 13 pages
Research Article

Complex Dynamics and Synchronization in a System of Magnetically Coupled Colpitts Oscillators

1Laboratory of Electronics and Signal Processing, Department of Physics, Faculty of Science, University of Dschang, P.O. Box 67, Dschang, Cameroon
2Department of Physics, Faculty of Science, University of Douala, P.O. Box 24157, Douala, Cameroon
3Laboratory of Mechanics and Modeling of Physical Systems, Faculty of Science, University of Dschang, P.O. Box 67, Dschang, Cameroon

Correspondence should be addressed to L. K. Kana; moc.oohay@lkgnagmek

Received 8 October 2016; Revised 13 February 2017; Accepted 23 February 2017; Published 10 April 2017

Academic Editor: Zhao Lu

Copyright © 2017 L. K. Kana 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.


We propose the use of a simple, cheap, and easy technique for the study of dynamic and synchronization of the coupled systems: effects of the magnetic coupling on the dynamics and of synchronization of two Colpitts oscillators (wireless interaction). We derive a smooth mathematical model to describe the dynamic system. The stability of the equilibrium states is investigated. The coupled system exhibits spectral characteristics such as chaos and hyperchaos in some parameter ranges of the coupling. The numerical exploration of the dynamics system reveals various bifurcations scenarios including period-doubling and interior crisis transitions to chaos. Moreover, various interesting dynamical phenomena such as transient chaos, coexistence of solution, and multistability (hysteresis) are observed when the magnetic coupling factor varies. Theoretical reasons for such phenomena are provided and experimentally confirmed with practical measurements in a wireless transfer.

1. Introduction

In electronics and in nonlinear dynamics, oscillators involving inductors are always modeled without any coupling from external magnetic field. However, in certain experimental condition, non-physically coupled systems may interact.

In this paper, we study the effect of the magnetic coupling on the behavior of a Colpitts coupled oscillator. Coupled chaotic systems have become a topic of active research in many scientific areas such as chemistry, biology, optics, and communications. In particular, the synchronization of coupled chaotic systems has received a lot of attention since the early 1990s [13]. This has been motivated by potential application of chaotic signals in such areas as signal encryption and communications [4].

The classical Colpitts oscillator is widely used in electronic devices and communication system as a source of sinusoidal waveform with low harmonic content [5]. This oscillator has been extensively used in the investigation of various dynamical phenomena. With a special parameters settings this oscillator exhibits a rich dynamic behavior: period doubling and chaotic behavior [57]. As a chaotic generator it can be used to generate noise-like broadband signals. The oscillator is very flexible as the fundamental frequency can be tuned from several kilohertz to several gigahertz, that is, to the microwave range.

Hyperchaos is characterized as chaotic dynamics with more than one positive Lyapunov exponent (LE). Hyperchaos has attracted increasing attention from the scientific and engineering communities [8, 9]. The interest devoted to hyperchaotic systems is motivated by their complex dynamics such as rich bifurcation, wider frequency bandwidths, and their potential application to information technologies, specifically to secure communication; hyperchaos has been firstly reported by Rössler [8], and the first circuit implementation of hyperchaos was realized by Matsumoto et al. [9]. Although the Colpitts oscillator is known to exhibit chaotic dynamics, it has also been shown to exhibit an hyperchaotic behavior when coupled in various configurations [10, 11].

In a complex system several attractors may coexist for a given set of system parameters. This coexistence is termed as multistability and has been found in almost all research areas of natural science, such as mechanics, electronics, biology, environmental science, and neuroscience [12, 13]. Multistable systems are characterized by a high degree of complexity in dynamical behavior due to the interaction among the coexisting attractors [14].

Generally speaking, chaos synchronization can be referred to as behavior in which two or more coupled systems exhibit identical chaotic oscillations. For such dynamical systems, the loss of synchronization between the subsystems is strictly related to the Lyapunov exponents (LEs) of the global system. When the subsystems lose synchrony, a transition from chaotic to hyperchaotic behavior takes place; that is, at least two LEs become positive [11]. Several approaches to chaotic synchronization have been proposed in the literature, for example, adaptive control, active-backstepping design, active control, backstepping, sliding mode control, and self-adaptive backstepping [14]. In most of these methods, the controllers and the design approach are often very complex and may be difficult to achieve in practice.

In this paper, a new method is introduced. It has been experimentally and has been numerically characterized, leading to a simple, cheap, and easy technique for the study of dynamic and synchronization of Colpitts coupled oscillator. When compared with other circuit techniques applied to the chaotic oscillator [11], our approach requires a simple control circuitry (based on the magnetic coupling: wireless coupling). A previous work in this domain dealt with chaos control and synchronization of a system. It consisted of chaotic modified Van der Pol-Duffing oscillator, which is able to generate low frequency chaotic signal [15].

The idea underlying our approach is also applied to chaotic Colpitts oscillator. It demonstrates how it is possible to control and synchronize a system in a simple and effective way, with high frequency oscillations. This coupling gives a possibility of obtaining a rich and complex dynamic system. To the best of our knowledge, the magnetic coupling of two chaotic Colpitts oscillators has never been investigated.

Our aims in this work are the following: (a) extend the study of magnetic coupling to higher frequencies chaotic systems (they are more useful in practical secure communication systems), (b) investigate the conditions leading to practical synchronization of the coupled Colpitts circuits, and (c) present an electronic implementation of the magnetically coupled systems. The rest of the paper is organized as follows: In Section 2, the circuit of the magnetically coupled system is presented and an appropriate model is derived to describe the dynamical behavior of the global oscillator. Section 3 deals with the complex dynamical behavior of the proposed circuit. In Section 4, based on the Lyapunov theory factor, we demonstrate the effectiveness of the proposed controller. In Section 5, numerical and experimental simulation results are provided to illustrate the effectiveness and feasibility of the control and synchronization technique. Finally, conclusions are given in Section 6.

2. Circuit Description of the Model and Differential Equations

2.1. The Model

The Colpitts oscillator shown in Figure 1(a) contains a bipolar junction transistor 2N3904 (BJT), resistors, a feedback capacitor , and a resonant circuit. The resonant circuit consists of an inductor and two capacitors and . The model of the transistor is shown in Figure 1(b), with a nonlinear element (voltage-controlled) and a linear current-controlled current source. This oscillator has been used in electronic devices and communication systems, due to the fact that it can exhibit chaos [5].

Figure 1: The electrical circuit of the system consisting of two mutually coupled Colpitts oscillators (a) and the circuit diagram of a bipolar junction transistor (BJT) (b). The circuit elements are 2N3904;  ;  ;  ;  ;  ;  ;  ; and .
2.2. State Equations

Let us consider the two magnetic core coils of Figure 1(a) and the magnetic coupling defined in [15]. The following set of differential equations are obtained by applying Kirchhoff’s laws to the circuit shown in Figure 1(a):where ; ,  , , and is the current gain (transistor gain) of the BJT. The relationship between the transistor gain and common base forward short-circuit gain of the transistor is defined as , is the mutual inductance ; if it is less than , then the coupling magnetic ( is the length of the coil and is the distance between the windings) [15] is also .

In our model (Figure 1), the values of electronic circuit components are  V, , ,  , , , and .

For the numerical investigation, an appropriate rescaling is introduced with the corresponding change of variables: ;  ,  ,  ,  ; one can rewrite (1) to obtain the following dimensionless set of differential equations:where

2.3. Equilibrium Points and Their Stability

We consider the normalized system (2). The equilibrium points of this system can be found by solving the system of algebraic equations:

We obtain the equilibrium points    with where

The equilibrium point does not depend on the control parameter . By linearizing system (2) around this point, we obtain the following Jacobian matrix:

The stability of the equilibrium point can be analyzed by solving the characteristic equation , where is the identity matrix. Then, the eigenvalues are obtained by solving the following characteristic equation: with where

Using the Routh-Hurwitz criteria on system (8), it appears that its trivial fixed point is stable if the following condition is satisfied: where(i)If , and , condition (11) can be solved choosing as the control parameter and keeping the other parameters constant; for this condition is not satisfied and the equilibria point is unstable saddle focus. Figure 2(a) illustrated such assertion, where we plot in the complex plane ; the imaginary part of eigenvalues as a function of their real part. We observe that, for , there exists the real part of the corresponding eigenvalues which is positive.(ii)If , , condition (11) is satiated for all values for . In Figure 2, we show the eigenvalues in the complex plane ; . Equation (8) is solved using the Newton-Raphson algorithm for the range of parameter : the real parts of the eigenvalues are all negative; then the equilibrium point is a stable spiral point. Let us consider the coefficient as the control parameter, since it can vary through experiment with the distance between the magnetic core coils. Applying the condition for the occurrence of the Hopf bifurcation: in (8) there exists a pair of purely imaginary complex eigenvalues and , where is called the critical value for the onset of the bifurcation. In the case of , the critical values of the parameter control are obtained by solving the equation:with

Figure 2: Representation of the eigenvalues, solution of the characteristic equation (11), in the complex plane (Re(, Im()) when (a) and when (b).

The critical value can be obtained by making recourse to numerical methods. We have one critical value ; the corresponding eigenvalues are ; ; ; ; ; ; ; and .

3. Dynamical Behavior of the Mutually Coupled Colpitts Systems

3.1. Bifurcation Analysis

As we mentioned above, one of the main advantages of this circuit for experimental applications is the richness of its dynamics. Here, the types of behaviors are identified using two indicators. The first indicator is the bifurcation diagram, the second being the graph of Lyapunov exponent’s spectra [7]. In order to define different types of behaviors in our model, system (4) is solved numerically through the fourth-order Runge-Kutta integration algorithm. For each set of parameters used in this work, the time step is always and the simulations are done with variables and constant parameters in extended mode [16, 17]. For each case, system (2) is integrated for a sufficiently long time and the transient is discarded with the time step being always . According to [18], the dynamics of the nonlinear system (2) can be classified in terms of the Lyapunov exponents as follows:(1)For an equilibrium point, (2)For a limit cycle (periodic orbits), (3)For 2-torus (quasiperiodic orbits), .(4)For 3-torus, (5)For chaotic orbits, , and .(6)For hyperchaotic orbits, , and .

The evolution process of the system is analyzed precisely by the means of the Lyapunov exponents spectrum, bifurcation diagrams, and phase portraits when varying the parameter in tiny step 0 < < 0.9; the other parameters are kept constant (the values defined in Section 2). The graph of the Lyapunov exponent spectrum and its corresponding bifurcation diagrams are depicted, in Figure 4. The bifurcation diagram is obtained by plotting local maxima of the states ; ; ; and . It can be seen that the bifurcation diagram coincides well with the spectrum of the Lyapunov exponents. In Figures 3 and 4, the following scenarios emerge when monitoring the control parameter : hyperchaos quasiperiodic hyperchaos chaos. Figure 5 shows some numerical phase portraits, obtained by Matlab simulation (1) and PSpice (2) analysis, respectively, when . Figure 5(a) shows a quasiperiodic attractor (). Figures 5(b) and 5(c) show, respectively, hyperchaos and chaos ( and ). Practically, the distance between two magnetically coupled Colpitts systems or the position or type of magnetic coupling core could have drastic effects on the oscillators dynamics, even if the number of oscillators is increased. This is an interesting behavior in the sense of control of the collective dynamics of magnetically coupled nonlinear oscillators [15].

Figure 3: Bifurcation diagram of magnetically coupled Colpitts systems (2) showing the coordinates , , , and . In terms of the control parameter , the other system parameters are ,  ,  ,  ,  ,  ,  , and .
Figure 4: Lyapunov exponents spectrum of magnetically coupled Colpitts systems (2) versus the control parameter . The rest of the parameters are those of Figure 3.
Figure 5: Qualitative comparison of numerical phase forms and experimental (PSpice) phase portraits of the oscillator, obtained for same specific values of the control parameter : (a) a quasiperiodic attractor for ; hyperchaotic attractor for ; and chaotic attractor for
3.2. Multistability or Coexistence of Attractors

Generally, multistable systems are characterized by a high degree of complexity in dynamical behavior due to the interaction among the coexisting attractors. First of all, the dynamics of a multistable system are extremely sensitive to initial conditions. Due to the coexistence of different attractors and complex fractal basin boundary structures very small perturbations of the initial state may influence the final attractor. Second, the qualitative behavior of the system often changes under the variation of the parameters system. Similarly, attractors exist only in small intervals of the parameters system. A slight change in a control parameter may cause a rapid change in the number and type of coexisting attractors. Third, multistable systems are extremely sensitive to noise. Noise may cause a popping process between various attractors. As an example of the system which possesses multistability, we consider the system described by (2) and consider as a control parameter. In Figure 6 we present the bifurcation diagrams of system (2), with the values of magnetic coupling on the horizontal axes. In the blue region of Figure 6 the values of have been calculated for the increasing values of while in the red region for the decreasing values of , period-3 and chaos coexist for ; chaos and period-5 coexist for ; period-3 and chaos coexist for , , and ; period-3 (increases) and period-5 coexist for ; and period-3 and period-6 coexist for .

Figure 6: Bifurcation diagrams for (blue) increasing and (red) decreasing of magnetic coupling.

For instance, In Figure 7 we present the coexistence of two different attractors, which can be obtained under the initial conditions , ,  ,  ;  ;  ,  ;  .

Figure 7: Coexistence of two different attractors, 3-period and hyperchaos, for and initial conditions .
3.3. Complex Transient Chaos and Hyperchaos

The appearance of chaos and hyperchaos on finite time scales is known as transient chaos and transient hyperchaos. In our system, the appearance of periodic or chaos or hyperchaos motion strongly depends not only on the circuit parameters but also on the initial conditions. When the magnetic coupling of the stimulus is selected as , a phenomenon of transient chaos is generated in the periodically coupled Colpitts oscillator (the trajectories of the system have a transition from transient chaotic to steady periodic behaviors with time evolutions). Figure 8(b) displays the time-domain wave form of variable , where chaotic attractor is located in the time interval and then a limit cycle with period-3 forms after ; Figure 8(b) depicts two phase portraits of chaos and period-3 attractors in different phase planes.

Figure 8: Transition from chaotic to period doubling: time-domain waveform of variable (a), the phase portrait in the time interval (b), and the phase portrait in the time interval (c) for initial conditions ;  ;  ;  ;  ;  ;  ;   when .

The transient hyperchaos phenomenon exists when : Figure 8 shows the time-domain waveforms and the phase portrait of the hyperchaotic attractor combined with period-3, where hyperchaotic attractor is located in the time interval and then a limit cycle with period-3 forms after . Figure 9(a) displays the time-domain wave form of variable , where hyperchaotic attractor is located in the time interval (Figure 9(b)), and then a limit cycle with period-3 forms after (Figure 9(c)).

Figure 9: Transition from hyperchaos to period doubling: time-domain waveform of variable (a), the phase portrait in the time interval (b), and the phase portrait in the time interval (c) for initial conditions ; ; ; ; ; ; ; when .

4. Synchronization of Two Modified Colpitts Oscillators

4.1. Case of Two Identical Oscillators

Although the coupling scheme implemented here is bidirectional, we need to define a master (drive) and a slave (response) system. We consider here that the drive system is given by the following set of coupled differential equations:

According to our aim, we define the controlled response system as the following set of differential equations:where is the control parameter. We define the error system as the difference between the signals from the master and the slave system as ,  ,  , and . By considering the time derivatives of the error signals together with (15) and (16), we obtain the error dynamics system: where .

It is clear that the synchronization problem is replaced by the equivalent problem of stabilizing (17) using a suitable choice of the control law .

Proposition 1. The synchronization problem for the magnetically coupled Colpitts system is to achieve the asymptote of the zero solution of the error system (17) in the sense that , as where .

Proof. Consider a Lyapunov function candidate in the following form:If the following condition holds with , the time derivative of function (18) with respect to the system of errors equations (17) satisfieswith ; ; and .
Since is positive definite function and is a negative function, it follows that the origin equilibrium point of system (17) is asymptotically stable. Therefore, from (20) we can easily show that is integrable with respect to time, for any initial condition as . This result concludes the proof.

4.2. Case of Nonidentical Oscillator (Mismatched Oscillators)

In this section, we consider two oscillators and with the different magnetic core coils, coupled via magnetic coupling (Figure 1). In this case the overall system is given by

A three-parameter diagram showing regions of synchronization (blue) in the () space computed, with the other parameter values stated in Section 2, is presented in Figure 10.

Figure 10: Three parameters of system , showing regions of synchronization (blue) and unsynchronization (brown) with the initial conditions ; ; ; ; ; ; ; .

5. Numerical and Experimental Simulations of the Controlled Synchronization

5.1. Numerical Study

In this section, we present numerical simulation results to verify the effectiveness of the controller. In all cases, we select the parameters , , , , , , , and such that the chaotic state in Figure 3 is maintained and use the following initial conditions for the drive-response system: , .

The graphs of Figure 11 illustrate the time evolutions of the synchronization errors. According to these graphs, one observes that the synchronization is reached, as the errors amplitudes are decreasing with time towards a small limit close to zero. Since Figure 11 presents the performance of the controller, the error dynamics are found to converge to the zero solution as , implying that the synchronization between systems (15) and (16) has been achieved.

Figure 11: Graphs of time variation of the synchronization errors , , , and .

Figure 12 presents the bifurcation diagrams, the synchronization boundaries of the system according to the coupling parameter .

Figure 12: Bifurcation diagram for the in-phase synchronization, respectively, , , , and .
5.2. Experimental Study

The complex behavior of the system (Figure 1) has been experimentally investigated. Various bifurcation and the synchronization state are examined, when the distance between the magnetic core coils is monitored, while keeping all values of electronic components constant. It is important to mention that the analog voltages obtained from our circuit are directly the voltages or of (1).

Figure 13 shows the photograph of the experimental Colpitts oscillators. This circuit is carried out using bipolar junction transistors (2N3904) and we set the following values: , , , and .

Figure 13: Experimental circuit.

Some examples are reported in Figure 14, in particular the synchronization between the oscillators obtained for .

Figure 14: In-phase synchronization of the magnetically coupled chaotic circuits (a) ; (b) ; and (c) . versus ; versus ; versus ; (4) versus .

6. Conclusion

The aim of this paper was to introduce a new method which is a simple, cheap, and easy technique for the study of dynamic and synchronization of Colpitts coupled oscillator. This method has gone through experiment and has been numerically characterized. The influence of the magnetic field via a magnetic coupling between the two oscillators has been presented. The stability of the fixed points has been investigated and the analytical formulas describing the stability of the system have been established. The bifurcation analysis of the phenomena leading to the creation of multistability of the system has been studied. Numerical simulations and electronic circuit implementation of the system show that the complex dynamics (transient chaos and transient hyperchaos) of the circuit are heavily dependent on the initial state of this system and on the magnetic coupling factor . We have shown that the magnetic coupling may be used to generate the hyperchaotic behavior. This can also drive two identical chaotic oscillators in a synchronized state as well as achieving the global stabilization of the system to its regular dynamics [15]. The dynamics of the coupled system have been investigated and Lyapunov stability theory has been applied to prove that, under some conditions, the drive-response system can achieve practical synchronization. The electronic circuit implementation has been realized to demonstrate the effectiveness of the controller. In future works, the following aspect will be considered: (i) the case of two nonidentical Colpitts oscillators and (ii) the collective dynamics of several magnetically coupled Colpitts systems.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


  1. E.-W. Bai and K. E. Lonngren, “Synchronization of two Lorenz systems using active control,” Chaos, Solitons and Fractals, vol. 8, no. 1, pp. 51–58, 1997. View at Publisher · View at Google Scholar · View at Scopus
  2. H. B. Fotsin and P. Woafo, “Adaptive synchronization of a modified and uncertain chaotic van der Pol-Duffing oscillator based on parameter identification,” Chaos, Solitons and Fractals, vol. 24, no. 5, pp. 1363–1371, 2005. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  3. U. E. Vincent, R. K. Odunaike, J. A. Laoye, and A. A. Gbindinninuola, “Adaptive backstepping control and synchronization of a modified and chaotic Van der Pol-Duffing oscillator,” Journal of Control Theory and Applications, vol. 9, no. 2, pp. 273–277, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  4. C. Li, X. Liao, and K.-W. Wong, “Lag synchronization of hyperchaos with application to secure communications,” Chaos, Solitons & Fractals, vol. 23, no. 1, pp. 183–193, 2005. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  5. M. P. Kennedy, “Chaos in the colpitts oscillator,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 41, no. 11, pp. 771–774, 1994. View at Publisher · View at Google Scholar · View at Scopus
  6. G. M. Maggio, O. De Feo, and M. P. Kennedy, “Nonlinear analysis of the Colpitts oscillator and applications to design,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 46, no. 9, pp. 1118–1130, 1999. View at Publisher · View at Google Scholar · View at Scopus
  7. J. Kengne, J. C. Chedjou, G. Kenne, and K. Kyamakya, “On the analysis of semiconductor diode-based chaotic,” Nonlinear Dynamics, vol. 1071, pp. 1301–1309, 2014. View at Google Scholar
  8. O. E. Rössler, “An equation for hyperchaos,” Physics Letters A, vol. 71, no. 2-3, pp. 155–157, 1979. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  9. T. Matsumoto, L. O. Chua, and K. Kobayashi, “Hyper chaos: laboratory experiment and numerical confirmation,” IEEE Transactions on Circuits and Systems, vol. 33, no. 11, pp. 1143–1147, 1986. View at Publisher · View at Google Scholar · View at MathSciNet
  10. H. Zhang, X.-K. Ma, Y. Yang, and C.-D. Xu, “Generalized synchronization of hyperchaos and chaos using active backstepping design,” Chinese Physics, vol. 14, no. 1, pp. 86–94, 2005. View at Publisher · View at Google Scholar · View at Scopus
  11. A. Čenys, A. Tamaševičius, A. Baziliauskas, R. Krivickas, and E. Lindberg, “Hyperchaos in coupled Colpitts oscillators,” Chaos, Solitons and Fractals, vol. 17, no. 2-3, pp. 349–353, 2003. View at Publisher · View at Google Scholar · View at Scopus
  12. M. D. Shrimali, A. Prasad, R. Ramaswamy, and U. Feudel, “The nature of attractor basins in multistable systems,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 18, no. 6, pp. 1675–1688, 2008. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  13. S. Kim, S. H. Park, and C. S. Ryu, “Multistability in coupled oscillator systems with time delay,” Physical Review Letters, vol. 79, pp. 1616–1619, 1997. View at Google Scholar
  14. U. Feudel, “Complex dynamics in multistable systems,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 18, no. 6, pp. 1607–1626, 2008. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  15. L. K. Kana, A. Fomethe, H. B. Fotsin, and P. H. Louodop, “A magnetic coupling based strategy for synchronization of a system consisting of chaotic modified van der pol dung oscillators,” Journal of Nonlinear Dynamics, 2014. View at Google Scholar
  16. G. P. King and S. T. Gaito, “Bistable chaos. I. Unfolding the cusp,” Physical Review A, vol. 46, no. 6, pp. 3092–3099, 1992. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  17. P. Ruiz, J. M. Gutiérrez, and J. Güémez, “Experimental mastering of nonlinear dynamics in circuits by sporadic pulses,” Chaos, Solitons and Fractals, vol. 36, no. 3, pp. 635–645, 2008. View at Publisher · View at Google Scholar · View at Scopus
  18. A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, “Determining Lyapunov exponents from a time series,” Physica D: Nonlinear Phenomena, vol. 16, no. 3, pp. 285–317, 1985. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus