Mathematical Problems in Engineering
Volume 2020 (2020), Article ID 5670652, 10 pages
https://doi.org/10.1155/2020/5670652

Research Article

The Dynamics Behavior of Coupled Generalized van der Pol Oscillator with Distributed Order

Asma Al Themairi¹ and Ahmed Farghaly2,3

¹Department of Mathematical Sciences, College of Science, Princess Nourah Bint Abdulrahman University, Riyadh, Saudi Arabia
²Department of Mathematics, Faculty of Science, Assiut University, Assiut 71516, Egypt
³Department of Basic Science, College of Computer and Information Sciences, Majmaah University, Al-Majmaah 11952, Saudi Arabia

Correspondence should be addressed to Asma Al Themairi; aialthumairi@pnu.edu.sa

Received 3 April 2020; Revised 18 June 2020; Accepted 6 July 2020; Published 28 July 2020

Academic Editor: Isabel S. Jesus

Copyright © 2020 Asma Al Themairi and Ahmed Farghaly. 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.

In this paper, we presented different behaviors such as chaotic and hyperchaotic of the generalized van der Pol oscillator with distributed order. We introduced the parameter intervals of these behaviors by computing the Lyapunov exponents of the oscillator, which is a good test for classifying the dynamical systems’ solutions. The active control approach with the Laplace transform technique was used to realize the antisynchronization and control of the proposed oscillator. Finally, numerical investigations have been carried out on the dynamics of the proposed oscillator to verify the reliability of our analytical results.

1. Introduction

In 1920, van der Pol invented the van der Pol oscillator [1]. It describes the oscillation of a triode in an electrical circuit. It is a fundamental mathematical model, where it has many numerous applications and exciting features. This oscillator is used in designing many biological models such as heartbeats [2], designing physical models such as mobile and phone oscillators [3], and modeling of electrical systems [4]. Mathematically, there are many versions of the van der Pol oscillator like

where (1) was introduced in 1927 by Van der Pol and Van Der Mark [4], and (2) was submitted in 1991 by Kapitaniak and Steeb [5]. In this paper, we will introduce the sinusoidal forcing with amplitude and frequency as a fractional version of the generalized van der Pol of the form

where , , , , and are constants.

Fractional calculus plays an essential role in modern science. It is a different and distinct method for dealing with nonlinear systems along with the integer order. Fractional order models are adequate for the description of dynamical systems rather than integer order models. We can recognize, describe, and know dynamic phenomena such as chaos, hyperchaos, synchronization, and some other aspects of fractional order models faster and more accurately than those of the integer order of nonlinear systems. At present, the application of fractional calculus in most scientific fields has attracted much attention. So, the fractional calculus on the dynamical system was essential and exciting, which had been investigated recently by many researchers [6–10]. Here, in our paper, we used distributed order as a type of fractional calculus to study the dynamic behavior of nonlinear generalized van der Pol oscillator.

Distributed order calculus has been investigated for the first time as the extinction of fractional order calculus by Caputo [11]. Caputo et al. introduced useful properties for the distributed order calculus [12–15]. Fernández-Anaya et al. proposed the Lyapunov theorem and several features for distributed order nonlinear dynamic systems [16]. The distributed dynamic systems have various applications in engineering and physics [17–20]. Some of the solutions’ properties for the viscoelastic rod derivative of the distributed order have been studied in [17]. Tavazoei [18] provided numerous conclusions about the monotonicity of responses stage describing distributed order structures in irrational transfer functions. The distributed time of order for Schrödinger-form equation implemented using the local Galerkin discontinuous method [19]. In [20, 21], the writers were safely using complex distributed order structures. On the other hand, chaos and hyperchaotic solutions for distributed order dynamical systems are essential topics. Chen et al. [20] introduced a chaotic distributed order Lorenz system. Mahmoud et al. [22] presented the chaotic complex distributed order Lüand Chen systems.

Synchronization of chaos has a critical part to play in dynamic systems. It has various applications in different fields [20, 21, 23–25] such as biology, physics, stable communication, and engineering. There are many methods of control to achieve synchronization between chaotic and hyperchaotic systems, such as linear feedback control [26], nonlinear feedback control [27], active control [28, 29], back-stepping design [30], tracking control [31], and adaptive control [24, 32]. These methods also used to hold synchronization between chaotic and hyperchaotic distributed order systems. Chen et al. [20] used the active control approach to synchronize the two chaotic Lorenz distributed order structures. Based on linear feedback control, Mahmoud et al. [22] presented the synchronization between chaotic complex distributed order Chen and Lü systems. By applyng the nonlinear feedback control and direct Lyapunov procedure, the synchronization among hyperchaotic complex distributed order van der Pol oscillators was investigated in [21].

System (3) can be constructed as a system of two differential equations of the first order when , ; then, we have

system (4) can be written with distributed order and complex version as

where , , .

The principal aims of this paper described as follows. (1) The hyperchaotic generalized van der Pol method has been introduced with complex parameter distributed order () of the form (4). (2) The dynamics of the system are analyzed, and we also evaluate the parameter intervals (, ) when the solution of this system is chaotic and hyperchaotic (3). The solution of system (3) is transformed into a periodic solution using linear feedback control. (4) A scheme to achieve antisynchronization between two generalized frameworks by van der Pol with distributed order hyperchaotic complex is stated. (5) The numerical simulations to test this theorem are presented.

The article is set out as follows. We are displaying some critical preliminaries in Section 2. Dynamics of the van der Pol generalized system of the distributed order complex and the parameter intervals at which chaotic and hyperchaotic solutions are calculated in Section 3. Using linear feedback control in Section 4, we manage its solution from chaotic and hyperchaotic to periodic. In Section 5, the antisynchronization between two identical systems of (3) is achieved through active control and transformation of the Laplace. The paper’s conclusion is presented in Section 6.

2. Preliminaries

The following section includes some definitions of the fractional order and the distributed order derivatives [16, 22, 33, 34], with useful remark and theorem that will be used later.

Definition 1. For any , the Caputo fractional derivative of function is given by

The Laplace transforms a fractional derivative of Caputo by

Definition 2. The distributed derivative of a continuous function is

where , , , , , .

Remark 1. The Laplace transform of the distributed derivative is given by:

where and .

Theorem 1. Let . If all poles of are in the open left-half complex plane, then .

3. Dynamics of the Complex Generalized van der Pol Oscillator with Distributed Order

The dynamics of the generalized van der Pol oscillator distributed order with complex parameters are studying in this section. We test the intervals of the parameters where there are chaotic and hyperchaotic approaches to the system.

The real form of system (5) can be written as

System (10) is seen as a generalization of several van der Pol oscillator variants (integer and fractional order). It has no fixed points because it is not an autonomous system. It is symmetric under the transformation .

Therefore, if is a solution of (10), then is also a solution of the same system. System (10) is also dissipative in the case of .

To show the solution’s behavior and to obtain the intervals for chaotic and hyperchaotic phenomena of system (10), we calculated the Lyapunov exponents, which are an excellent test to classify the solutions.

3.1. Lyapunov Exponents of System (10)

System (10) has four Lyapunov exponents , , which have been measured by Wolf algorithm [35]. The classification of signs of Lyapunov exponents are stated in Table 1 [36]:

Table 1: Signs from exponents of Lyapunov and the corresponding solution form.

Dynamics

−
−
−
−
Solution goes to equilibrium point

0
−
−
−
Periodic solution (limit cycles)

0
0
−
−
Quasiperiodic solution (2-torus)

0
0
0
−
Quasiperiodic solution (3-torus)

+
−
−
−
Chaotic behavior

+
+
−
−
Hyperchaotic behavior of order 2

+
+
+
−
Hyperchaotic behavior of order 3

Assuming that is an initial condition and , for the following cases, we determined Lyapunov exponents of system (10):

(i)Case 1: Fix  =  and vary such that and calculate Lyapunov exponents. Some results are shown in Table 2 which plotted in Figure 1.

Figure 1: Lyapunov exponents of system (10) for , and .

Figure 2: Lyapunov exponents of system (10) for , and .

Figure 3: Lyapunov exponents of system (10) for , and .

Table 2: Lyapunov exponent of system (10) for three cases.

Case 1: Fix  =  and vary

B

Dynamics

26
6.0041
3.1382
−6.1039
−6.7879
Hyperchaotic behavior of order 2

28
9.3759
3.7670
−4.5096
−6.2100
Hyperchaotic behavior of order 2

30
7.4983
4.1911
−5.3663
−7.3997
Hyperchaotic behavior of order 2

Case 2: Fix  =  and vary

Dynamics

21
8.0825
4.8029
−10.4143
−11.0868
Hyperchaotic behavior of order 2

25
7.4983
4.1911
−5.3663
−7.3997
Hyperchaotic behavior of order 2

27
8.7146
4.6502
−4.9561
−6.7303
Hyperchaotic behavior of order 2

Case 3: Fix  =  and vary

Dynamics

0.7
1.9130
−0.0303
−10.6506
−11.3279
Chaotic behavior

2
7.4983
4.1911
−5.3663
−7.3997
Hyperchaotic behavior of order 2

9
7.1457
3.4489
−6.1396
−7.6788
Hyperchaotic behavior of order 2

From Table 2, the solution of system (10) is hyperchaotic.

(ii)Case 2: Fix  =  and vary such that , we got the results of Lyapunov exponents shown in Table 2, plotted in Figure 2.

This implies our system (10) will have a hyperchaotic solution.

(iii)Case 3: Fix  =  and vary such that , we got Lyapunov exponents as shown in Table 2, plotted in Figure 3. Table 2 clearly shows that system (10) gives a chaotic and hyperchaotic solution of some intervals.

4. Supervision of Unusual Approaches of (10)

We applied the linear feedback control method to transform the system’s chaotic and hyperchaotic solution (10) to periodic one. System (10) after the introduction of the control functions became

where . The example , , , and was chosen to make control. As shown in Figures 4 and 5, the hyperchaotic solution converted to periodic solution.

(a)

(b)

(c)

(d)

Figure 4: The distributed order hyperchaotic complex generalized van der Pol system with control (11): (a) graph, (b) graph, (c) graph, and (d) graph.

(a)

(b)

Figure 5: The hyperchaotic complex distributed by order generalized van der Pol system with control (11): (a) space, and (b) space.

5. Chaos Synchronization of System (10)

In this section, using active control and Laplace transform, we present a theorem for achieving antisynchronization among the two same systems of (10). We consider (10) the master system, and the slave system takes the form

where is a suitable control function, which will be given in the last section.

Definition 3. The drive system (10) is said to achieve antisynchronization with the response system (12), if

where is vector of synchronization error, and is the matrix norm.
The system of error can be written as

Theorem 2. The antisynchronization between the master system (10) and the slave system (12) will be achieved if the control functions chose as follows:

where is the gain matrix.

Proof. Using the control functions (15), the error system (14) can be written as

By transforming system (16) by Laplace and applying Remark 1 in , , then we get:

By using Theorem 1 and Remark 1, we deduced that , .. Anti synchronization between the master system (10) and the slave system (12) can therefore be accomplished.
Numerically, if we take , and the initial values of the master system (10) and the slave system (12) are , respectively, and the same parameters of Figure 6, the antisynchronization between master system (10) and slave system (12) achieved as shown in Figures 7 and 8. Figure 7 shows the antisynchronization between the state variables in the master and slave systems. The errors of the synchronization approach to zero, as shown in Figure 8. It clear that there exists an agreement between numerical simulations and Theorem 2.

Figure 6: The hyperchaotic solution of system (10) in space.

(a)

(b)

(c)

(d)

Figure 7: Synchronization of the state variables of the master system (10) and the slave system (12): (a) and versus , (b) and versus , and (c) and versus , (d) and versus .

(a)

(b)

(c)

(d)

Figure 8: Synchronization error of the master system (10) and the slave system (12): (a) graph, (b) graph, (c) graph, and (d) graph.

6. Conclusions

In this work, we have investigated a new generalized van der Pol oscillator distributed order with a complex parameter (10). The literature is to be a generalization of several variants of van der Pol oscillator. We calculated Lyapunov exponents of that system for several parameter values, and we noticed that the system contains chaotic and hyperchaotic phenomena depicted in Figures 6, 9, and 10. We applied the linear feedback control method of system (10), and we could convert the chaotic and hyperchaotic solution to periodic one, as shown in Figures 4 and 5. We introduced Theorem 2 to fulfill the antisynchronization between two identical distributed order generalized van der Pol oscillators by using active control and Laplace transform method; the numerical simulations were consistent with the analytical study. The results are shown in Figures 7 and 8.

Figure 9: The hyperchaotic solution of system (10) in space.

Figure 10: The hyperchaotic solution of system (10) in space.

Data Availability

The authors declare that all data sources are original.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors are thankful to the deanship of scientific research Princess Nourah Bint Abdulrahman University, for supporting and funding the work through the research funding program under grant number (FRP144028).

References

B. Van der Pol, “Theory of amplitude of free and forced triode vibration,” Radio Rev, vol. 1, pp. 701–710, 1920.

L. Glass, Theory of Heart, Springer, Berlin, Germany, 1990.

B. Van der Pol, “On “relaxation-oscillations”,” The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, vol. 2, no. 11, pp. 978–992, 1926.

B. Van der Pol and J. Van Der Mark, “Frequency demultiplication,” Nature, vol. 120, no. 3019, pp. 363-364, 1927.

T. Kapitaniak and W.-H. Steeb, “Transition to hyperchaos in coupled generalized van der Pol equations,” Physics Letters A, vol. 152, no. 1-2, pp. 33–36, 1991.

M. Kavyanpoor and S. Shokrollahi, “Dynamic behaviors of a fractional order nonlinear oscillator,” Journal of King Saud University—Science, vol. 31, no. 1, pp. 14–20, 2019.

P. Pirmohabbati, A. h. Refahi Sheikhani, H. Saberi Najafi, and A. Abdolahzadeh Ziabari, “Numerical solution of full fractional Duffing equations with cubic-quintic-heptic nonlinearities,” AIMS Mathematics, vol. 5, no. 2, pp. 1621–1641, 2020.

S. Li, J. Niu, and X. Li, “Primary resonance of fractional-order Duffing-van der Pol oscillator by harmonic balance method,” Chinese Physics B, vol. 27, no. 12, Article ID 120502, 2018.

C. L. Ejikeme, “Solution to nonlinear duffing oscillator with fractional derivatives using homotopy analysis method (HAM),” Global Journal of Pure and Applied Mathematics, vol. 14, no. 10, pp. 1363–1388, 2018.

F. Wang and Y. Yang, “Fractional order barbalat’s lemma and its applications in the stability of fractional order nonlinear SystemsFractional order barbalat’s lemma and its applications in the stability of fractional order nonlinear systems,” Mathematical Modelling and Analysis, vol. 22, no. 4, pp. 503–513, 2017.

M. Caputo, Elasticita e dissipazione, Zanichelli, Bologna, Italy, 1969.

M. Caputo, “Mean fractional-order-derivatives differential equations and filters,” Annali dell’Universita di Ferrara, vol. 41, pp. 73–84, 1995.

M. Caputo, “Distributed order differential equations modelling dielectric induction and diffusion,” Fractional Calculus and Applied Analysis, vol. 4, no. 4, pp. 421–442, 2001.

R. L. Bagley and P. J. Torvik, “On the existence of the order domain and the solution of distributed order equations-Part I,” International Journal of Applied Mathematics, vol. 2, no. 7, pp. 865–882, 2000.

R. L. Bagley and P. J. Torvik, “On the existence of the order domain and the solution of distributed order equations-Part II,” International Journal of Applied Mathematics, vol. 2, no. 8, pp. 965–988, 2000.

G. Fernández-Anaya, G. Nava-Antonio, J. Jamous-Galante, R. Muñoz-Vega, and E. G. Hernández-Martínez, “Asymptotic stability of distributed order nonlinear dynamical systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 48, pp. 541–549, 2017.

T. M. Atanackovic, M. Budincevic, and S. Pilipovic, “On a fractional distributed-order oscillator,” Journal of Physics A: Mathematical and General, vol. 38, no. 30, pp. 6703–6713, 2005.

M. S. Tavazoei, “Fractional/distributed-order systems and irrational transfer functions with monotonic step responses,” Journal of Vibration and Control, vol. 20, no. 11, pp. 1697–1706, 2014.

T. Aboelenen, “Local discontinuous Galerkin method for distributed-order time and space-fractional convection-diffusion and Schrödinger-type equations,” Nonlinear Dynamics, vol. 92, no. 2, pp. 395–413, 2018.

J. Chen, C. Li, and X. Yang, “Chaos synchronization of the distributed-order Lorenz system via active control and applications in chaotic masking,” International Journal of Bifurcation and Chaos, vol. 28, no. 10, Article ID 1850121, 2018.

G. M. Mahmoud, A. A. Farghaly, T. M. Abed-Elhameed, S. A. Aly, and A. A. Arafa, “Dynamics of distributed-order hyperchaotic complex van der Pol oscillators and their synchronization and control,” The European Physical Journal Plus, vol. 135, no. 1, 2019.

G. M. Mahmoud, T. Aboelenen, T. M. Abed-Elhameed, and A. A. Farghaly, “Generalized Wright stability for distributed fractional-order nonlinear dynamical systems and their synchronization,” Nonlinear Dynamics, vol. 97, no. 1, pp. 413–429, 2019.

G. M. Mahmoud, E. E. Mahmoud, and A. A. Arafa, “On projective synchronization of hyperchaotic complex nonlinear systems based on passive theory for secure communications,” Physica Scripta, vol. 87, no. 5, Article ID 055002, 2013.

G. M. Mahmoud, A. A. Farghaly, T. M. Abed-Elhameed, and M. M. Darwish, “Adaptive dual synchronization of chaotic (hyperchaotic) complex systems with uncertain parameters and its application in image encryption,” Acta Physica Polonica B, vol. 49, no. 11, p. 1923, 2018.

T. Carletti, R. Serra, I. Poli, M. Villani, and A. Filisetti, “Sufficient conditions for emergent synchronization in protocell models,” Journal of Theoretical Biology, vol. 254, no. 4, pp. 741–751, 2008.

M. Rafikov and J. M. Balthazar, “On control and synchronization in chaotic and hyperchaotic systems via linear feedback control,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 7, pp. 1246–1255, 2008.

G. M. Mahmoud, T. M Abed-Elhameed, and A. A. Farghaly, “Double compound combination synchronization among eight n-dimensional chaotic systems,” Chinese Physics B, vol. 27, no. 8, Article ID 080502, 2018.

G. M. Mahmoud, T. M. Abed-Elhameed, and M. E. Ahmed, “Generalization of combination-combination synchronization of chaotic n-dimensional fractional-order dynamical systems,” Nonlinear Dynamics, vol. 83, no. 4, pp. 1885–1893, 2016.

G. M. Mahmoud, M. E. Ahmed, and T. M. Abed-Elhameed, “Active control technique of fractional-order chaotic complex systems,” The European Physical Journal Plus, vol. 131, p. 200, 2016.

S. Bowong and F. M. Moukam Kakmeni, “Synchronization of uncertain chaotic systems via backstepping approach,” Chaos, Solitons & Fractals, vol. 21, no. 4, pp. 999–1011, 2004.

G. M. Mahmoud, M. E. Ahmed, and T. M. Abed-Elhameed, “On fractional-order hyperchaotic complex systems and their generalized function projective combination synchronization,” Optik, vol. 130, pp. 398–406, 2017.

Z. Gao, Y. Wang, L. Zhang, Y. Huang, and W. Wang, “The dynamic behaviors of nodes driving the structural balance for complex dynamical networks via adaptive decentralized control,” International Journal of Modern Physics B, vol. 32, no. 24, Article ID 1850267, 2018.

K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, Wiley, Hoboken, NY, USA, 1993.

D. G. Duffy, Transform Methods for Solving Partial Differential Equations, Chapman and Hall/CRC, London, UK, 2004.

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.

G. M. Mahmoud, M. A. Al-Kashif, and S. A. Aly, “Basic properties and chaotic synchronization of complex Lorenz system,” International Journal of Modern Physics C, vol. 18, no. 2, pp. 253–265, 2007.