# Computer and Hardware Modeling of Periodically Forcedϕ6-Van der Pol Oscillator

Authors: A. O. Adelakun; A. N. Njah; O. I. Olusola; S. T. Wara
Journal: Active and Passive Electronic Components
 (2016)

---

## Abstract

Numerical simulation results for the dynamics ofϕ6-systems abound in the literature but their experimental results are yet to be known. This paper presents the chaotic dynamics of ϕ6-Van der Pol oscillator via electronic design, simulation, and hardware implementation. The results obtained are found to be in good agreement with numerical simulation results. The condition for stability of the fixed points is also computed and the method of multiple time scale is used to investigate the dynamical behaviour of the system. Therefore, the ϕ6-circuits which have rich dynamics and may have important applications in secure communications, random number generations, cryptography, and so forth have been practically implemented.


---

## Body

## 1. Introduction

In recent years the study of chaotic phenomena in self-excited oscillators with parametric and external-excitation has attracted much attention [1–6]. Among the 2-dimensional periodically forced oscillators, the most extensively studied examples are the Van der Pol oscillators [7, 8], Duffing oscillators [9], and Rayleigh oscillators [10]. Many works on self-excited, parametrically excited, and externally excited oscillators are also investigated in the literature. Recently, the rich dynamics of Van der Pol oscillator have been explored by theoretical analysis and numerical simulations [3, 11].Theϕ4-Van der Pol oscillator is a very significant classical model circuit that has been studied and even modified in some studies as reported by King et al. [12, 13], Fotsin et al. [14], and Fodjouong et al. [15]. However, with the advent of the study of chaotic motion by means of strange attractors, Poincare’ map, fractal dimension, it is necessary to seek for a better understanding of nonlinear systems with higher order nonlinear terms such as the ϕ6-Van der Pol oscillator.Aϕ6-Van der Pol oscillator or extended Van der Pol oscillator is the nonlinear oscillator with fifth-order nonlinear term. Its numerical simulation has been studied extensively by Siewe Siewe et al. [16, 17], Tchoukuegno et al. [18–20], and Yu and Li [21] and its synchronization behaviours for possible applications in information processing and secure communication by Njah et al. [22–24], and its application in the modeling of chemical reactions and brain activity by Chen and Dong [25] and Ott et al. [26]. In the past three decades, implementation of chaotic systems in power electronics has gone through intense development in many aspects of technology [27–29], including power devices, control methods, circuit design, computer-aided analysis, passive components, and packaging techniques. Like many areas of engineering, power electronics is mainly motivated by practical applications, and it often turns out that a particular circuit topology or system implementation has found widespread applications long before it is thoroughly analyzed and most of its subtleties are uncovered [30, 31].The aim of this paper is to design, simulate, and implement periodically forcedϕ6-Van der Pol oscillator by transforming the state equations into electronic circuit via analog components modeling. We also generate phase portraits for the triple-well and double-well attractors of the system. To the best of our knowledge, computer simulation with electronic software and hardware implementation of ϕ6-Van der Pol oscillators have not been studied. The rest of the paper is organized as follows: Brief description of the conditions for stability of fixed points for unperturbed system is discussed in the next section. Section 3 explains the numerical simulation of the oscillator for at least six physically relevant situations with corresponding attractors while Section 4 explains the analog simulations with electronic software (multisim) and hardware implementation of the system. Finally, Section 5 concludes the paper.


## 2. Analysis of the Unperturbed System

The general form of second-order differential equation of Van der Pol systems is given as(1)x¨=μ1-x2x˙-dVxdx+fcos⁡ωt,where μ<0 is the damping parameter and f and ω are amplitude and angular frequency, respectively, and its potential V(x) can be expressed in Taylor series as(2)Vx=12αx2+14βx4+16λx6,where α, β, and λ are constant parameters. The condition for stability of the fixed points of oscillator (1) for an unperturbed system can be assumed when μ=f=0. Thus, the system can therefore be expressed as a set of first-order differential equations of the form(3)x˙=y,y˙=-αx-βx3-λx5which corresponds to an integrable Hamiltonian system with the potential function in (2) and whose associated Hamiltonian function is(4)Hx,y=12y2+12αx2+14βx4+16λx6.From (3) and (4), we can compute for the fixed points and analyze their stabilities after some algebraic manipulations using the following conditions.(i)
Forβ2<4λα, the unperturbed system has only one fixed point (0,0), the centre, for which the system is single-well if α>0 and β>0 and single-hump potential if α<0 and β<0 as shown in Figures 1(a) and 1(b), respectively.(ii)
Forβ2=4λα with α<0, β>0, and λ>0, the unperturbed systems have three fixed (equilibrium) points, that is, one centre (0,0) and two saddles connected by two heteroclinic orbits (±x2,0), (5).(iii)
Forβ2>4λα with α>0, β<0, and λ>0, the unperturbed systems have five fixed points: one centre (0,0), two saddles connected by two heteroclinic orbits (±x2,0), and two saddle points connected to each other by one homoclinic orbit (±x1,0), (5) and (6).In previous works, Siewe Siewe et al. [16, 17] have shown that, for μ=0, system (4) is a Hamiltonian system with a pair of heteroclinic orbits defined as(5)xhett=±2x1sinhγ/2t-ξ+cosh⁡γt,yhett=±2γξ-1coshγ/2t2-ξ+coshγt3/2,and a pair of symmetric homoclinic trajectories connecting each unstable point to itself given by(6)xhomt=±2x1coshγ/2tξ+cosh⁡γt,yhomt=±2γξ-1sinhγ/2t2ξ+cosh⁡γt3/2,where(7)x1=-12λβ+Δ,x2=-12λβ-Δ,γ=x122λρ2-1,ξ=5-3ρ23ρ2-1,ρ2=β-Δβ+Δ,Δ=β2-4αλ.Figure 1
The smooth potential function,V(x), of the ϕ6-oscillators for (a) single-well: α=0.5 and β=0.175, (b) single-hump: α=-0.5 and β=-0.175, (c) double-well: α=-0.46, β=0.7, and λ=0.1, (d) double-hump: α=1.5, β=-0.7, and λ=-0.1, (e) triple-well: α=0.46, β=-0.7, and λ=0.1, and (f) triple-hump: α=-1, β=0.7, and λ=-0.1.
(a) (b) (c) (d) (e) (f)


## 3. Numerical Simulation

The general form of the periodically forced Van der Pol oscillator model is given by the second-order nonautonomous differential equation (1), where x is the state variable, dot(s) over x denotes the derivative with respect to time t, f is the amplitude of periodic forcing, ω is the angular frequency, μ>0 is the damping parameter, and V(x) is the potential. The potential for six physically relevant situations (single, double, and triple-well and hump) is shown in Figure 1. They correspond to different choices of the set of values for the parameter set α, β, and λ. For instance, the potential is double-well if (α< 0, β > 0, and λ > 0), double-hump if (α > 0, β < 0, and λ < 0), triple-well if (α > 0, β < 0, and λ > 0), and triple-hump if (α < 0, β > 0, and λ < 0).Van der Pol oscillator models describe various physical, electrical, mechanical, and engineering devices, for different potentialsV(x). The potentials V(x) are approximated by finite Taylor series for the ϕ6-chaotic oscillators as shown in (3), where α, β, and λ are constant parameters of the potential. The introduction of a new parameter λ in (2) yields ϕ6-Van der Pol equation (9) that can be expressed as(8)x¨=μ1-x2x˙-αx-βx3-λx5+fcos⁡ωt.Numerical simulation of (8) using fourth-order Runge-Kutta method gives the phase portraits presented in Figure 2 for double and triple scroll attractors. From the phase portraits it is obvious that the systems are chaotic for the parameter values used.Figure 2
Phase portrait of the chaotic attractor for theϕ6-Van der Pol oscillator with parameters values (a) μ=0.01, α=0.46, ω=0.86, f=9.0, β=1.0, and λ=0.1 for the double-well potentials and (b) μ=0.4, α=1.0, ω=3.14, f=9.0, β=-0.7, and λ=0.1, for the triple-well.
(a) (b)


## 4. Design of Analog Simulator

Circuit realization ofϕ6-Van der Pol system is shown in Figure 3 (schematic circuit) and Figure 4 (analog circuit) in which the nonlinear resistor (NR) can be modeled by using i-ϕ nonlinear characteristic given by the following relation:(9)iNR=dqdϕ·dϕdτ=αoϕ+βoϕ3+δoϕ5,where αo, βo, and δo are parameters of the nonlinear resistor gNR, ϕ is the flux over the resistor, q is the charge, C is the capacitance, L is the inductance, R is the resistance, V=Ecosωτ, and iNR is the current through the nonlinear resistor. Applying Kirchoff’s law to the circuit shown in Figure 3, we have(10)VR+VL+VC+VNR=Ecos⁡ωτ,iR+iL+iC+iNR=0.Then (10) becomes(11)d2ϕdt2-LC1-ϕ2dϕdt+RL2Cϕo+ϕo3+ϕo5=ELCcos⁡ωt.Taking μ=L/C, t=τ/LC, α=R/L2Cαo, β=R/L2Cβo, f=E/LC, and δ=R/L2Cδo, we can rewrite electronic equation (11) in form of flux as(12)d2ϕdt2-μ1-ϕ2dϕdt+αϕ+βϕ3+δϕ5=fcos⁡ωt.The electronic analog simulator for ϕ6-Van der Pol oscillator can be designed using the electronic multiplier, AJ633AN, which perform the multiplication operations. The analog circuit operates over a dynamic range of ±15D.C. The integrators are operational amplifiers, U5 and U6 (UA741CD) with feedback capacitors, and summations of the state variables are obtained using operational amplifiers U7 with the help of feedback resistors and multiple input resistors. Also, the simulator output showing the dynamics when varying the input components can be viewed directly on a software oscilloscope using an appropriate time scaling. The designed circuitry realization for the oscillator consists of two channels A and B which realize the two state variables x and x˙, respectively. The attractors were generated from ϕ6-Van der Pol circuit of Figure 4 by varying(13)λ=R9R8.The computer simulation of ϕ6-Van der Pol circuit on electronic software (multisim) (Figure 5(a)) and its hardware implementation on electronic breadboard (Figure 5(b)) both reproduced the numerical results in Figure 2 on the oscilloscope. Indeed, the practical implementation of the aforementioned system has application in technology such as in secure communication.Figure 3
Schematic circuit of Van der Pol model.Figure 4
Analog circuit ofϕ6-Van der Pol oscillator.Figure 5
Phase portrait ofϕ6-Van der Pol attractors for (a) computer simulation and (b) experimental implementation for (i) double-well at R8=50Ω and R9=200kΩ and (ii) triple-well at R8=5kΩ and R9=50kΩ.
(a) (b)


## 5. Conclusion

In this paper, we have used the analog electronic software (multisim) for the electronic circuit design and simulation ofϕ6-Van der Pol oscillator and have practically implemented the designed circuit on electronic breadboard using analog components. Our experimental results reproduced those obtained via numerical simulation and multisim simulation. Therefore the rich dynamics of extended Van der Pol oscillator which have wide application in science and technology have been practically implemented, for possible technological application such as in secure communication.



---

*Source: 3426713-2016-12-15.xml*