Abstract

This paper describes a new autonomous deterministic chaotic dynamical system having a single unstable saddle-spiral fixed point. A mathematical model originates in the fundamental structure of the class C amplifier. Evolution of robust strange attractors is conditioned by a bilateral nature of bipolar transistor with local polynomial or piecewise linear feedforward transconductance and high frequency of operation. Numerical analysis is supported by experimental verification and both results prove that chaos is neither a numerical artifact nor a long transient behaviour. Also, good accordance between theory and measurement has been observed.

1. Introduction

Heretofore, formation of the strange attractors was reported in many lumped electronic circuits, both isolated and driven. Chaotic waveforms have several fundamental properties: wideband and continuous frequency spectrum, sensitivity to small changes of initial conditions, dense state attractors with (unlike noise) non-integer metric dimension, and absence of closed-form analytic solution. Because of these features, especially presence of dominant harmonics in frequency spectrum, chaos can be observed in general models of all types of lumped analogue harmonic oscillators, regardless of type of active element used as energy source. For example, paper [1] describes a long-term unpredictable motion of the well-known Colpitts oscillator. Work [2] demonstrates existence of robust strange attractors unfolded during study of motion of the Hartley oscillator. Both oscillators require only single transistor and three accumulation elements for generation of the stable sinusoidal output voltages. Some interesting papers [3, 4] introduce a methodology on how to transform very popular Wien-bridge network configuration into chaotic oscillator. It turns out that simple mechanism dedicated to amplitude stabilization can cause both vector field stretching and folding required for the chaos evolution. In fact, the most famous chaotic electronic system, Chua’s oscillator, can be considered as connection of LC parallel resonant tank, timing network, and nonlinear resistor. This nonlinear element can be either polynomial [5] or piecewise linear [6]. Both types of nonlinearities are practically implemented using several commercially available devices and are the only active components in designed chaotic circuits. Study [7] proves that pair of coupled accumulation elements can be substituted by a single higher-order element such as the frequency dependent negative resistor. Single third-order immittance element can be used to create linear part of the vector field as it is demonstrated in paper [8]. Systematic design procedure leading to robust chaotic oscillators with the passive only nonlinear elements is provided in a couple of cookbooks [9, 10]. In these articles, nonlinear resistors have the ampere-voltage curve with single breakpoint and therefore serves as energy valve. Several papers utilize proposed approach to construct original structures of chaotic oscillators [11, 12].

Besides oscillators, unpredictable motion was detected in other naturally non-chaotic analogue building blocks. Recently, in the modern complex electronic systems, chance to observe robust chaos is strengthened because of high frequencies and amplitudes of useful processed signals. In such cases, parasitic properties (especially input and output impedances) of active elements necessarily introduce error terms in describing set of differential equations, lead to unexpected couplings between independent nodes, change values of working accumulation elements and increase mathematical order of circuit, cause additional functional nonlinearities or hysteresis effects on existing nonlinear transfer curves, etc. Let us provide a few examples of such systems. Presence of structurally stable strange attractors was confirmed in standard topologies of DC-DC converters, namely, buck [13], cuk [14], and boost [15, 16]. Very interesting reading about irregular behaviour in general class of DC-DC converters and its properties, including possible control mechanism, can be found in papers [17, 18]. Another transparent example of chaotic system based on switched power converter is described in study [19]. Impulsive chaos can be observed in switched capacitor circuits as well; see papers [20, 21] for details. Several research papers were dedicated to describing chaos in phase-locked loops: for details, consult several examples provided in papers [22, 23]. Recently, single-spiral as well as double-scroll attractors were detected within high frequency mathematical models of multistate static memory cells composed of the resonant tunnelling diodes; the number of these diodes determines total amount of fixed points. Chaos can be observed if ampere-voltage characteristic is considered piecewise linear [24] as well as cubic polynomial [25]. An interesting consequence coming from these studies is that there is a direct correlation between number of equilibria and complexity of generated chaotic attractors. Very complex vector field and geometric patterns of basins of attraction [26] suggest possible evolution of multiscroll chaotic attractors. Fingerprints of chaos can be observed in very simple networks that contain recently developed memory circuit elements: memristors [27, 28] or mem-elements in general [29, 30], both with integer and fractional order [31]. To end this, very brief, list of chaotic analogue building blocks, we can mention frequency filters as good candidates to exhibit chaotic solution. These filtering sections can be of different topologies, such as the state variable filters [32] or passive ladder filters with active nonlinear feedback. This active nonlinear feedback can be implemented by either two-terminal [5, 33, 34] or two-port [35] device. There, existence of chaos is conditioned by small system dissipation, i.e., filters with high quality factor. Speaking in terms of lumped chaotic circuit with minimum number of elements, it seems that driven R-L-diode circuit is the simplest one [36].

Searching for the most trivial algebraic description of continuous time chaotic flow was started by fundamental question raised in [37]. To date, this problem still represents up-to-date topic.

Existence of robust chaotic waveforms within very simple lumped electronic systems having one or two transistors is old but still an up-to-date problem. Some recent examples of transistor-based chaotic oscillators can be found in paper [38]. Circuits discussed here clearly belong to the simplest cases exhibiting chaotic dynamics because, apart from the transistor used, some comprise only three accumulation elements. A slightly more complicated chaotic network having two standard bipolar transistors is discussed in research work [39].

The main purpose of upcoming sections of this paper is to prove that chaos can be observed, under very specific conditions, in the fundamental cell of class C amplifier comprising single bipolar transistor. This fact has been already mentioned in a very recent publication [40]. Unlike this research, where both backward and forward transconductances of bipolar transistor two-port are considered as polynomial nonlinearity, this work assumes a case much closer to physical reality. In addition, two differently shaped strange attractors instead of one are presented and thoroughly analysed. Although describing mathematical model that gives basis to this work and corresponding circuit realization is not the simplest one, it still represents interesting chaotic dynamical system that belongs to class of systems having a single unstable equilibrium.

This paper is organized as follows. The next section provides a brief mathematical description of autonomous deterministic dynamical system to be analysed. The third section of this paper numerically quantifies chaotic motion produced by analysed system. The fourth section describes two flow-equivalent circuit designs based on class C amplifier cell with polynomial nonlinearity. The fifth section of this manuscript is dedicated to experimental verification of constructed chaotic oscillators. Finally, concluding remarks and future possible topics are revealed.

2. Mathematical Background

The principal structure of class C single transistor-based amplifier is given in Figure 1(a). For useful signals, this network can be redrawn as computational model shown in Figure 1(b). In its basic form, both circuits are linear and valid for small useful signals. Without input waveform, the resulting dynamical system can be described by the following set of ordinary differential equations:

Note that admittance parameters of bipolar transistor y11, y12, and y22 are considered constant. This represents common situation in many practical applications. Forward transconductance is supposed to be nonlinear, characterized by curve approximated by cubic polynomial of form . Symbolic form of Jacobi matrix isand characteristic polynomial calculated for single fixed point located at origin is

In further text, let us assume that biasing point forces transistor to work in a regime of constant current source; i.e., y22 = 0 S. Also, suppose that normalized values of accumulation elements are Cce = 0 F, Cbe = C = 1 F, L = 1 H and backward transconductance is linear; i.e., y12 = 1 S. Considering this, characteristic polynomial simplifies into

For the purpose of the upcoming optimization, namely, to reduce size of hyperspace of parameters that undergoes search, limit cases for sought parameters become important

3. Numerical Analysis

All numerical integration results presented in this paper were achieved by using four-order Runge–Kutta method with fixed step size. This routine was also used as core for calculation of the largest Lyapunov exponent (LLE) [41] together with Gram–Smith orthogonalization. LLE as the flow quantifier forms one fitness function for non-gradient optimization as admitted in papers [42, 43]. In fact, searching for chaotic solution of dynamical system is a multi-objective optimization problem. Firstly, chaotic behaviour should not be interchanged with long transient motion. Therefore, LLE needs to be calculated on the attractor. Also, slowly divergent motion needs to be removed from interesting solution. These potentially useful results are stored and finally verified by handmade analysis and attractor visualization. Boundedness of numerically integrated ω-limit set is verified via inside-sphere condition. Whole optimization routine was implemented in MATLAB. Since evaluation of individual objective functions is not mutually connected, searching process could be significantly accelerated by employing CUDA-based parallel processing. Nevertheless, searching within three-dimensional space of internal system parameters is still a very demanding operation.

Note that the optimization procedure described above and its results, i.e., numerical values of input admittance of transistor and shape of the cubic polynomial forward transconductance, do not have direct correlation with physical nature of one stage class C amplifier. However, admittance is a scalable network quantity (via impedance norm) and, for the large processed amplitudes of useful signal, saturation-type amplification property of bipolar transistor is mild hypothesis. Of course, normalized value of y12 = 1 S indicates strongly non-unilateral bipolar transistor stage.

Two sets of internal parameters leading to different shapes of strange attractors and associated non-integer Kaplan–Yorke dimension DKY have been discovered, namely,

Eigenvalues associated with fixed point for dynamical system with (6) areand for set (7) we get

Therefore, origin is saddle-node equilibrium point with stability index 1.

Figure 2 demonstrates shapes of typical strange attractors generated by dynamical system (1) with parameters (6). Numerical integration was calculated with final time 300 s and time step 10 ms. Because of vector field symmetry with respect to the horizontal plane z = 0, there are two mirrored strange attractors with the same geometric structure. Figure 2 also shows sensitivity to the small deviations in the initial conditions and vector field near fixed point  = (0, 0, 0)T. For this analysis, final time and time step were kept the same. Figure 3 is a visualization of LLE as two-dimensional function of shape of forward transconductance function . For values of input admittance y11 bigger than 0.6, chaos disappears, and class C amplifier produces various limit cycles. Order in circuit is restored thanks to the significant system dissipation. Of course, cases with y11 ≤ 0 are out of the question because of the physical essence of the input admittance. Note that hidden attractors including those with fractal geometric dimension are not excluded by investigations presented in this paper. Figure 4 shows distribution of dynamic energy over state space volume occupied by strange attractor. To visualize the local dynamics calculated at each point in 510 × 610 grid, numerical integration was up to final time 1 s with time step 10 ms. Value ranges for individual edges are uniform:  ∈ (–2.5, 2.5) and  ∈ (–3, 3).

Figure 5 provides information about shape of the strange attractor if dynamical system (1) has the internal parameters (7). It strongly resembles the so-called double hook attractor that can be generated by Chua’s circuit with three real eigenvalues associated with fixed point located at origin [44]. Unlike Chua’s circuit, there are two coexisting double hooks. Also note evolution of two neighbouring orbits plotted with respect to time: difference is visible after quite a long time; about 120 s lapses to phase shift and 170 s expires to observe exponential divergence.

3.1. Fractional-Order Case

Recently, discovery and presentation of new chaotic dynamical system are usually supported by numerical analysis of describing mathematical model under condition of fractionality; that is, one or all differential equations are of non-integer order. Such investigation could be performed on three levels. The first level is pure mathematical and involves discretization of ordinary [45] or partial [46] differential equations and performing iterations. This numerical approach is well suited for chaotic [4749] as well as hyperchaotic [50, 51] dynamics. The second level is to model fractional-order (FO) circuit element in the form of two-terminal or two-port using suitable circuit-oriented software such as OrCAD PSpice or Multisim. Both classes of FO elements (often called constant phase elements due to their fundamental properties keeping constant phase shift between response and driving force) could be considered either ideal or approximated in the frequency domain. To reach precise approximation of FO two-terminal fully passive circuit element, paper [52] could be consulted. Design procedure toward approximated FO two-port element is provided in a comprehensive study [53]. Since chaotic signal is wideband, the mentioned approximations should be valid, at least, within the same frequency range. To avoid laborious calculations, passive ladder FO capacitors of different math orders can be found in [54].

Because of practically oriented concluding remarks coming from this research study, the last method dedicated for modelling non-integer nature of accumulation elements associated with class C amplifier has been adopted. New circuit structure to be analysed is provided in Figure 1(c). Note that FO nature of accumulation elements is modelled by fully passive ladder networks having seven sections. This number turns to be enough, from the viewpoint of both required frequency range (five decades) and maximal allowed phase error (that is lower than 1 degree). To increase universality, let us denote mathematical order of the approximation associated with capacitor Cbe by k, sum of capacitors C + Cce by letter m, and math order of inductor by n. State vector changes into a more complex form

Behaviour of this lumped electronic system is uniquely determined by the following set of ordinary differential equations:where and stand for voltage across jth capacitor in approximation of Cbe and C + Cce FO capacitor, respectively. Similarly, symbol icj represents current flowing through jth inductor. Note that second, fourth, and sixth differential equation are presented in our system kth, mth, and nth times, respectively. Obviously, we are experiencing system of 3 + k + m + n first-order ordinary differential equations. Thus, instead of calculation of the binomial coefficients and discretization of differential equations, we can follow the standard way of numerical integration process. Some numerical results are graphically demonstrated via Figure 6. For mathematical orders of the approximated constant phase elements below 0.9, chaotic behaviour of class C amplifier disappears and for very low total orders of circuit only limit cycles were observed. Symbols α, β, and δ used in Figure 6 represent FO of first, second, and third differential equation, respectively. Let us discuss a specific FO case where orders of individual differential equations equal 0.9.

By adopting numerical values provided in paper [54] for mathematical order 0.9, frequency range of constant phase element can be downscaled from original interval 10 Hz up to 1 MHz. To analyse normalized circuit, pseudo-capacitance of the constant phase element will be 1 s0.1/F (FO capacitor) and pseudo-inductance approximately equals 1 s0.1/F (FO inductor). For the first case, resulting numerical values of resistors and capacitors will be

Of course, these values hold for second FO capacitor as well; i.e., index a can be replaced by b. FO inductor can be modelled by values

LLE can be utilized to obtain piecewise linear approximation of a smooth polynomial vector field [55]. In our case, transconductance could be substituted by the three-segment odd-symmetrical function with parameters

For the first discovered set of “chaotic” parameters (6), piecewise linear approximation can be adopted withwhile for the second set of parameters (7) leading to evolution of robust chaos optimal piecewise linear approximation is

The corresponding strange attractor is provided in Figure 7 together with the plot of dynamic energy distribution over three horizontal planes of state space. Note the high degree of attractor similarity with the results provided in Figure 5.

An interesting question that could be raised here is the possibility to generate the so-called multi-scroll [56] or multi-grid [57] strange attractors. This kind of complex behaviour is associated with specific geometry of vector field. State space can be divided into “periodic” formation of the affine segments with the same properties: existence of single fixed point with saddle-spiral local geometry with stability index 1. In each mentioned segment, state trajectory is attracted toward eigen-plane and spirals away to another segment. Famous double-scroll attractor observed in Chua’s oscillator can be considered as ancestor for this kind of motion and special type of multi-spiral attractor. However, eigenvalues associated with parameter sets (15) areand for set (16) we obtain

These eigenvalues do not lead to vector field configuration promising from the viewpoint of generation of multi-spiral attractor. Also, numerical investigation does not reveal such kind of dynamical motion. In outer regions of the state space, i.e., if  < – or  > , dynamical system exhibits strong dissipation along all directions and state orbit is pushed toward inner segment as demonstrated by means of Figure 8.

Figure 9 demonstrates complex nature of basins of attraction calculated for dynamical system (1) with parameter group (6). Individual colours distinguish between three ω-limit sets: chaotic attractor (red), solution close to limit cycle (green), and fixed point (black). Due to vector field symmetry, only horizontal state space slices having z ≥ 0 are provided. Figure 10 represents the same plots but related to parameter set (7).

By introducing constant term to polynomial function , fixed point located at the origin turns into virtual. For small values of parameter c, this equilibrium with less dynamical system still generates strange attractors like the original. Thus, there is the overlap of discovered chaotic dynamical system into quite rare cases without equilibrium.

4. Circuit Design

To meet the basic criteria for presentation of new chaotic system [58], structural stability of desired strange attractor needs to be proven. Because of limited capability of human eye to catch short time transient movement, practical construction and laboratory verification of chaotic oscillator is the best way how to demonstrate such stability.

If a set of ordinary differential equations and numerical values of internal parameters is known, the problem of circuit synthesis belongs to simple and straightforward task with multiple correct solutions. If minimization of circuits elements is not design priority, we can follow the concept of analogue computers. Only three building blocks are required: inverting integrators, differential amplifiers, and two-ports with a nonlinear transfer characteristic. Final network can work in either voltage [59, 60] or current mode [61]. Under specific circumstances, resulting circuit can be simplified by substitution of passive two-ports [62]. More complex nonlinear transfer curves can be implemented by digital subcircuits [63]. The method leading to full on-chip implementation of chaotic oscillator can be found in an example in paper [64]. Circuity realization of chaotic oscillator based on model (1) is shown by means of Figure 11(a). It contains two voltage-feedback operational amplifiers TL082, current feedback operational amplifier with the frequency compensation node (TZ) AD844, and two four quadrant analogue multipliers AD633 with versatile voltage transfer function  = K⋅( – ) ( – ). All state variables can be easily measured as voltages at output of inverting integrators. Circuit is fed by symmetrical ±15 V voltage supply. Dynamical behaviour of this network is described by the following set of ordinary differential equations:where K = 0.1 represents internally trimmed transfer constant of multiplier.

Nonideal and parasitic properties of used active devices clearly limit the frequency spectrum of generated signal. This could be especially important in the case of oscillators where output waveforms (state variables) are the broad-band noise-like chaotic signals. Fortunately, time constant of designed chaotic oscillator can be shortened and resulting waveforms judged, both using suitable circuit simulator or laboratory measurement. Because of better correspondence between theoretical expectations and circuit analysis, OrCAD PSpice program has been utilized for this purpose. Individual active elements were modelled using third level of abstraction and influential parasitic properties swept parametrically. From the viewpoint of designed circuit, the most problematic part is AD633 exhibiting heavily distorted input-output characteristics for time constant of chaotic oscillator shorter than 1 μs. Further interesting reading about this topic can be found in [65].

Figure 11(b) recalls duality principle and construct circuit that is dual to Figure 1(b). It comprises two current-controlled voltage sources, but only z21 (i1) is nonlinear. This circuit is described by the same differential equations as (1) but after considering dual variables (voltages/currents) and parameters (admittances/impedances). In detail,where state vector becomes x = (i1, i2, )T. Direct circuitry realization of the state description of class C amplifier (1) can be found in Figure 11(c). The corresponding set of the ordinary differential equations arewhere is transconductance of operational transconductance amplifier LM13700 precisely adjusted by external DC current. Because the desired strange attractors are located inside a cube with edges ±3 V, dynamical range of integrated circuit LM13700 does not present a problem. The whole chaotic circuit can be supplied by symmetrical ±15 V. To turn around output current of LM13700, input voltage nodes should be swapped.

5. Experimental Verification

The chaotic oscillator proposed in the previous section was constructed using breadboard; see Figure 12. For experimental verification, time constant was set to τ = 10 μs. This fact leads to choosing all capacitors 10 nF and, using basic calculations, values of remaining passive circuit components: R1 = 50 kΩ, R2 = 10 kΩ, R3 = 3.9 kΩ, R4 = 10 kΩ, R5 = 127 Ω, and R6 = 10 kΩ. To trace evolution of chaos, two resistors are considered variable, namely, R3 (parameter b) and R5 (parameter a). Figure 12 also shows realization of printed circuit board of two-port described by admittance parameters with fully adjustable transconductances; i.e., shapes of both and can be adjusted independently. Individual measurement results are provided by means of Figures 1316.

6. Conclusion

Starting with a very simple mathematical model of class C amplifier with single bipolar transistor, two geometrically distinct strange attractors have been discovered. There are two conditions for chaos evolution: non-unilaterality and nonlinear forward transconductance of transistor. Numerical analysis includes calculation of LLE, energy distributions, and basins of attraction.

Based on this finding, a robust chaos generator is designed and experimentally verified within this paper. Circuitry realization follows analogue computer concept and utilizes cheap and commercially available active elements. Several different shapes of the self-excited strange attractors were observed during experimental verification.

The problem solved in this work opens new topics for future research. For example, existence of chaos in common building blocks dedicated for analogue signal processing such as multi-stage amplifiers, cascodes, current mirrors, and differential amplifiers could be proven via the same fully numerical approach as proposed in this paper. Also, staying focused on topology of class C amplifier, “chaotic” configurations of internal parameters that are closely related to standard operational regime could be localized. Particular attention should be paid to the promising cases where backward transconductance is very low. Finally, input signal as a driving force acts as an additional degree of freedom that can increase the overall probability of chaotic solution.

Data Availability

There are two groups of results (data outputs) presented in this manuscript. 1. Those coming from numerical routines used for mathematical analysis of investigated dynamical system (integration of state trajectories, sensitivity properties, calculation of the largest Lyapunov exponent, and visualization of basins of attraction). 2. Experimental verification of designed chaotic oscillator. This type of output can be re-simulated by using OrCAD PSpice circuit simulator. In both cases, source files, schematics, and simulation profiles can be forwarded to candidates based on the e-mail request submitted to the author of this paper, that is, [email protected]. Please do not hesitate to contact the author also in the case of some questions, suggestions, and comments. Details about input parameters for numerical procedures can be found in the main body of the paper.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author would like to thank Dr. Tomas Gotthans for his long-time support, much useful advice, and recommendations about the concept of this paper. The research described in this paper was financed by Grant Agency of Czech Republic under project no. 19-22248S. For research, the infrastructure of the SIX Center was used.