Abstract and Applied Analysis

Volume 2014 (2014), Article ID 835482, 10 pages

http://dx.doi.org/10.1155/2014/835482

## Hopf Bifurcations of a Stochastic Fractional-Order Van der Pol System

^{1}State Key Laboratory for Strength and Vibration of Mechanical Structures, Xi’an Jiaotong University, Xi’an 710049, China^{2}School of Mathematics and Statistics, Tianshui Normal University, Tianshui 741001, China

Received 27 March 2014; Revised 21 May 2014; Accepted 21 May 2014; Published 9 June 2014

Academic Editor: Juan J. Trujillo

Copyright © 2014 Xiaojun Liu 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 Hopf bifurcation of a fractional-order Van der Pol (VDP for short) system with a random parameter is investigated. Firstly, the Chebyshev polynomial approximation is applied to study the stochastic fractional-order system. Based on the method, the stochastic system is reduced to the equivalent deterministic one, and then the responses of the stochastic system can be obtained by numerical methods. Then, according to the existence conditions of Hopf bifurcation, the critical parameter value of the bifurcation is obtained by theoretical analysis. Then, numerical simulations are carried out to verify the theoretical results.

#### 1. Introduction

Fractional calculus is a topic of more than 300 years old. The idea of fractional calculus has been known since the regular calculus. However, its development is very rapid just in recent several decades due to its application in physics, engineering, secure communications, and so on. Meanwhile, it has been found that many fractional-order nonlinear systems can demonstrate chaotic behavior, such as fractional-order Chua circuit [1], fractional-order Lorenz system [2], and fractional-order Chen system [3]. These examples and many other similar samples perfectly clarify the importance of consideration and analysis of dynamical systems with fractional-order models [4–6].

The VDP oscillator represents a nonlinear system with an interesting behavior that exhibits naturally in several applications. It has been used for study and design of many models including biological phenomena, such as the heartbeat, neurons, acoustic models, and radiation of mobile phones, and as a model of electrical oscillators. The deterministic VDP has very rich dynamical behaviors, such as saddle-node bifurcation, symmetry-breaking bifurcation, period-doubling bifurcation, Hopf bifurcation, and chaos. Recently, many researchers studied the bifurcation and periodic solutions of the VDP system in detail [7–11]. As research on fractional differential equations goes up, the VDP model with fractional order was proposed and investigated in [12, 13]. It was found that when the order is less than 1, the fractional-order VDP system also has rich dynamics similar to the corresponding integer-order one.

In the real world, uncertainty due to measuring, materials, manufacturing, and assembling is inevitable. These uncertainties usually can be modeled as random parameters with certain statistics. A stochastic system means a system with random parameters of given statistics, which is a mathematical model close to the real world. There are several methods to deal with the dynamical problems of stochastic systems. The first one is the Monte Carlo method [14], which is simple and universal but usually involves a great deal of computational effort. The second one is the stochastic finite element method [15], which involves the least computation but is usually restricted to systems with random variables of small perturbations only. And the third one is orthogonal polynomials [16–20], which is based on the expansion theory of orthogonal polynomials, without the limitation of small perturbations. There are a lot of references about the analysis of stochastic systems via orthogonal polynomials. For example, in [21, 22], the Chebyshev polynomial approximation was applied to analyze the bifurcation and chaos of stochastic Duffing system and stochastic Duffing-Van der Pol system. The period-doubling bifurcation for double-well stochastic Duffing system was investigated via Lagurre polynomial approximation in [23]. Recently, the synchronization for a stochastic fractional-order system with a random parameter via Laguerre polynomial was analyzed in [24]. However, the literature on the investigation of a stochastic fractional-order VDP system is few.

Inspired by the above discussion, in this paper, Hopf bifurcations of a fractional-order VDP system with a random parameter are studied. Firstly, the Chebyshev polynomial approximation method is applied to investigate the stochastic fractional-order system. Based on this method, the stochastic system is reduced to the equivalent deterministic one. Then, according to the existence conditions of Hopf bifurcation, the critical parameter value of the bifurcation is obtained by theoretical analysis. Meanwhile, numerical simulations are performed to verify the theoretical results. Besides, the Hopf bifurcation with the variation of derivative order is also observed by numerical computations.

The paper is organized as follows. In Section 2, the definitions for the fractional calculus and numerical algorithms are given. In Section 3, the Chebyshev polynomials are introduced as the basis of orthogonal polynomial approximation. In Section 4, the transformation of stochastic fractional-order VDP system into its equivalent deterministic one by Chebyshev polynomial approximation is shown. Section 5 is devoted to studying the Hopf bifurcation of the stochastic system. Finally, we summarize the results in Section 6.

#### 2. Fundamentals of Fractional Derivative

##### 2.1. Definition

Fractional calculus is a generalization of integration and differentiation to a noninteger-order integrodifferential operator which is defined by where is the derivative order which can be a complex number and is the real part of . The numbers and are the limits of the operator. There are many definitions for fractional derivatives. Three most frequently used ones are the Grunwald-Letnikov definition, the Riemann-Liouville, and the Caputo definitions.

The Grunwald-Letnikov definition (GL) derivative with fractional-order is described by where the symbol means the integer part.

The Riemann-Liouville (RL) definition of fractional derivatives is given by where is the gamma function.

The Caputo () fractional derivative is defined as follows:

It is well known that the initial conditions for the fractional differential equations with Caputo derivatives take the same form as for the integer-order ones, which is very suitable for practical problems [25]. Therefore, we will use the Caputo definition for the fractional derivatives in this paper.

##### 2.2. Numerical Algorithms

Obtaining numerical solutions of a fractional differential equation is not easily compared with that for ordinary differential equations. There are two approximation methods which can frequently be used to numerical computations on chaos and bifurcations with fractional differential equations. One is an improved version of Adams-Bashforth-Moulton algorithm based on the predictor-correctors scheme [26], which is a time-domain approach. The other is a method, known as frequency domain approximation [27], based on numerical analysis of fractional-order systems in the frequency domain.

Simulations for fractional-order systems using the time domain methods are complicated and due to long memory characteristics of these systems require a very long simulation time but on the other hand, it is more accurate [28]. Therefore, we employ the improved predictor-corrector algorithm for fractional-order differential equations in this paper.

In order to get the approximate solution of a fractional-order chaotic system by the improved predictor-corrector algorithm, the following equation is considered: where is just the value rounded up to the nearest integer and is the ordinary derivative of . Formula (5) is equivalent to the Volterra integral equation Now, for the sake of simplicity, we assume that we are working on a uniform grid with some integer and set . Using the standard quadrature techniques for the integral in (6), denote ; the integral is replaced by the trapezoidal quadrature formula at point where is the piecewise linear interpolant for with nodes and knots chosen at the , (). After some elementary calculations, the right hand side of (7) gives And if we use the product rectangle rule, the right hand of (7) can be written as where Then the predictor and corrector formulae for solving (6) are given, respectively, by The approximation accuracy of scheme (11) is .

#### 3. The Chebyshev Polynomials

In general, random parameters in some engineering structures are bounded in nature. So the probability density function (PDF for short) model for the bounded random variables is taken to be an arch-like PDF of the form [29]:

Figure 1 shows the arch-like PDF curve for bounded random variable . As the orthogonal polynomial basis for this kind of PDF, the second kind of Chebyshev polynomials is the only choice. The general expansion for the second kind of Chebyshev polynomials can be described as follows: Then, we have The recurrent formula for the Chebyshev polynomials of the second kind is The orthogonality of Chebyshev polynomial of the second kind can be expressed as Equation (16) represents a weighted orthogonal relationship. The weighting function is just the same as that for the PDF of the random parameter which obeys arch distribution, in (12), and the left-hand side of (16) can be regarded as the expectation of the product . Owing to the orthogonality of Chebyshev polynomial, any measurable function can be expanded as where . This expansion is called an orthogonal decomposition of random function , which is the theoretical base of orthogonal decomposition methods.

#### 4. The Stochastic Fractional-Order VDP System

In this paper, we consider the fractional-order VDP system with a random parameter. The system can be described by the following differential equations: where is the fractional order and , are state variables. , are deterministic parameters and is a random parameter, and it can be expressed as where and are the mean value and standard deviation of , respectively, and is regarded as the intensity of . is a random variable defined on .

It is well known that the stochastic function space, which is constituted with the responses of nonlinear dynamic systems with random parameters, has proved to be a Hilbert space with respect to convergence in the mean square. Therefore, the fractional-order system (18) with random parameters can be reduced into the equivalent deterministic system by using of the orthogonal polynomial expansion. It follows from the orthogonal polynomial approximation that the responses of system (18) can be expressed approximately by the following series: where represents the th Chebyshev polynomial and is the largest order of the polynomials. Then substituting (20) into the system (18), we can get and with the aid of the recurrent formulas of Chebyshev polynomial, all the nonlinear terms in the system can be further reduced into linear combination of as follows: where can be obtained with the help of the computer algebraic systems, such as Maple or MATLAB. Meanwhile, the term can be reduced to Substituting (22) and (23) into (21), then we can have Multiply both sides of (24) by , in sequence and take expectation with respect to . Then we can finally obtain the equivalent deterministic system. Remember that if in (24), then its left is strictly equivalent to the response of stochastic system (18). Otherwise, (24) is just approximately valid with a minimal residual error. Due to the requirement of computational precision, we take in the following numerical analysis. Therefore, we get the equivalent deterministic system approximately as We can get the numerical solutions , , of system (25) by the improved predictor-corrector algorithm. Therefore, the approximate random responses of the stochastic system (18) can be expressed as The ensemble mean responses (EMRs for short) of the stochastic fractional-order system are described as Especially for , namely, , this sample system is significant for reference. It is usually called a mean parameter system, of which the responses (SRMs for short) may be approximated as

#### 5. Hopf Bifurcations

A Hopf bifurcation occurs when a periodic solution or limit cycle, surrounding an equilibrium point, arises or goes away as a parameter varies. In this section, the Hopf bifurcation for the stochastic fractional-order VDP system will be studied with theoretical and numerical computations. Obviously, both the stochastic VDP and equivalent deterministic systems (18) and (25) only have one equilibrium point, namely, the origin (0,0). The Jacobian matrix for the linearized system evaluated at the origin is The corresponding characteristic equation is where are the coefficients of (30) and can be obtained by computing via Maple or MATLAB. According to the bifurcation theory of nonlinear dynamical systems and taking as a bifurcation parameter [22], we can obtain the existence conditions of Hopf bifurcation for the system. (1)The eigenvalues of the linearized system about the fixed point are .(2)The real and image parts of the eigenvalues satisfy , , , when .(3)The real parts of other eigenvalues are not equal to zeros.

For the high order of equivalent system (25), so it is hard to get the threshold value of bifurcation parameter directly. Therefore, we can get the threshold via the following theorem.

Theorem 1. *Let be the characteristic equation evaluated at the equilibrium point, the -dimensional Routh-Hurwitz determinant, and , () the coefficient of the characteristic equation. If there exists a threshold of bifurcation parameter which conforms to the following terms, then the first two existence conditions of Hopf bifurcation are satisfied:*(1)*, , ;*(2)*, ;*(3)*.*

According to Theorem 1, the following parameter relations which can make the equality are derived: Then, (31) is substituted into the 8-dimensional and 7-dimensional Routh-Hurwitz determinants , , and only the following parameters can make , satisfied: Meanwhile, the relations which can make all the coefficients satisfied () are We substitute the relations (33) into and find that is not equal to 0; therefore, according to Theorem 1, the first two existence conditions of Hopf bifurcation are satisfied. When (33) is substituted into the characteristic equation, all the eigenvalues of it can be obtained: Obviously, when the intensity of random parameter , the real parts of all the other eigenvalues are less than 0 except and . The characteristic equation does not have zero solution for . Therefore, all the existence conditions of Hopf bifurcation are satisfied. When , the Hopf bifurcation occurs around the equilibrium .

The parameters are taken as , , and , and the derivative order is . And initial conditions are and . According to the above analysis, we can get that the critical parameter value of the bifurcation is , which implies that if , the system will converge to the equilibrium; otherwise, it will converge to a limit cycle. The parameter values are chosen as and , respectively. Then we obtain the time trajectories and phase diagrams for the EMR and SRM of the stochastic system by simulation computations; see Figures 2 and 3. It is clearly seen that the system converges to the fixed point for and to a limit cycle for , which implies that the Hopf bifurcation occurs in the stochastic system (18). Meanwhile, the SRM and EMR of phase plots coincide well enough. It means that the Chebyshev polynomial approximation method is effective for the stochastic fractional-order system (18).

When the others parameters are fixed and , we also found that Hopf bifurcation occurs with the variation of the derivative order. The stochastic system converges to the equilibrium when and to a limit cycle for , which means that the Hopf bifurcation occurs when the derivative order ; see Figure 4, from which it is obtained that the derivative order can also be taken as a bifurcation parameter for fractional-order system.

The simulation time is taken as , one integration step of time is taken as , derivative order , and divide the time interval into 8 equal parts. The time trajectory and phase diagram for are shown in Figure 5, from which it appears that the converges to a fixed point as . However, the time trajectory is zoomed in when belongs to the interval ; see Figure 6(a), from which we can see that the trajectory gradually decreases but does not converge to a point directly. Secondly, the time trajectory is zoomed in when ; see Figure 6(b). Meanwhile, we also obtain the trajectories for , , and which are shown in Figures 6(c)–6(e), respectively. Comparing these figures, it is seen that the amplitude of decreases as the time increases. Based on these observations, we can suggest that the EMRs of the stochastic system (18) have an extremely long transient process while converging to the equilibrium as .

#### 6. Conclusions

In this paper, the Hopf bifurcation of a fractional-order VDP system with a random parameter is analyzed. An equivalent deterministic system is obtained with the approximation principle of Chebyshev orthogonal polynomials. Based on bifurcation theory of nonlinear dynamic systems, the critical parameter value of the bifurcation is obtained by theoretical analysis. Meanwhile, numerical simulations are performed to verify the theoretical results. Besides, the derivative order can also be taken as a bifurcation parameter, and the Hopf bifurcation is observed with the variation of the order.

It is worth mentioning that the orthogonal polynomial approximation is an effective method for integer-order stochastic systems, but the universality of the method for fractional-order ones still remains an open question. Due to the absences of analytical methods, the threshold value of bifurcation when the derivative order is taken as a bifurcation parameter cannot be obtained by theoretical analysis but only the numerical experiments.

#### Conflict of Interests

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

#### Acknowledgment

This work is supported by the National Natural Science Foundation of China (NSFC) under Grant nos. 11332008, 11172223, and 11172224.

#### References

- C. P. Li, W. H. Deng, and D. Xu, “Chaos synchronization of the Chua system with a fractional order,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 360, no. 2, pp. 171–185, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - Z. M. . Ge and C. Y. Ou, “Chaos in a fractional order modified Duffing system,”
*Chaos, Solitons & Fractals*, vol. 34, no. 2, pp. 262–291, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. G. Lu and G. Chen, “A note on the fractional-order Chen system,”
*Chaos, Solitons and Fractals*, vol. 27, no. 3, pp. 685–688, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - H. Sun, W. Chen, and Y. Chen, “Variable-order fractional differential operators in anomalous diffusion modeling,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 388, no. 21, pp. 4586–4592, 2009. View at Publisher · View at Google Scholar · View at Scopus - Y. Liang and W. Chen, “A survey on computing Lévy stable distributions and a new MATLAB toolbox,”
*Signal Processing*, vol. 93, no. 1, pp. 242–251, 2013. View at Publisher · View at Google Scholar · View at Scopus - Y. J. Liang and W. Chen, “Reliability analysis for sluice gate anti-sliding stability using Lévy stable distributions,”
*Signal Processing*, 2014. View at Publisher · View at Google Scholar - V. Venkatasubramanian, “Singularity induced bifurcation and the van der Pol oscillator,”
*IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications*, vol. 41, no. 11, pp. 765–769, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - R. Mettin, U. Parlitz, and W. Lauterborn, “Bifurcation structure of the driven van der Pol oscillator,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 3, no. 6, pp. 1529–1555, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - U. Parlitz and W. Lauterborn, “Period-doubling cascades and devils staircases of the driven van der Pol oscillator,”
*Physical Review A*, vol. 36, no. 3, pp. 1428–1434, 1987. View at Publisher · View at Google Scholar · View at Scopus - A. Buonomo, “The periodic solution of van der Pol's equation,”
*SIAM Journal on Applied Mathematics*, vol. 59, no. 1, pp. 156–171, 1998. View at Publisher · View at Google Scholar · View at MathSciNet - J.-X. Xu and J. Jiang, “The global bifurcation characteristics of the forced van der Pol oscillator,”
*Chaos, Solitons and Fractals*, vol. 7, no. 1, pp. 3–19, 1996. View at Publisher · View at Google Scholar · View at MathSciNet - R. S. Barbosa, J. A. T. MacHado, B. M. Vinagre, and A. J. Caldern, “Analysis of the Van der Pol oscillator containing derivatives of fractional order,”
*Journal of Vibration and Control*, vol. 13, no. 9-10, pp. 1291–1301, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - J.-H. Chen and W.-C. Chen, “Chaotic dynamics of the fractionally damped van der Pol equation,”
*Chaos, Solitons and Fractals*, vol. 35, no. 1, pp. 188–198, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. Shinozuka, “Probability modeling of concrete structures,”
*ASCE Journal of the Engineering Mechanics Division*, vol. 98, no. 6, pp. 1433–1451, 1972. View at Google Scholar · View at Scopus - R. Ghamem and P. Spans,
*Stochastic Finite Element: A Spectral Approach*, Springer, Berlin, Germany, 1991. - C. L. Pettit and P. S. Beran, “Spectral and multiresolution Wiener expansions of oscillatory stochastic processes,”
*Journal of Sound and Vibration*, vol. 294, no. 4, pp. 752–779, 2006. View at Publisher · View at Google Scholar · View at Scopus - O. P. Le Maître, O. M. Knio, H. N. Najm, and R. G. Ghanem, “Uncertainty propagation using Wiener-Haar expansions,”
*Journal of Computational Physics*, vol. 197, no. 1, pp. 28–57, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. P. Le Maître, H. N. Najm, R. G. Ghanem, and O. M. Knio, “Multi-resolution analysis of Wiener-type uncertainty propagation schemes,”
*Journal of Computational Physics*, vol. 197, no. 2, pp. 502–531, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Xiu and G. E. Karniadakis, “Modeling uncertainty in steady state diffusion problems via generalized polynomial chaos,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 191, no. 43, pp. 4927–4948, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Wan and G. E. Karniadakis, “An adaptive multi-element generalized polynomial chaos method for stochastic differential equations,”
*Journal of Computational Physics*, vol. 209, no. 2, pp. 617–642, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. L. Leng, C. L. Wu, X. P. Ma, G. Meng, and T. Fang, “Bifurcation and chaos analysis of stochastic Duffing system under harmonic excitations,”
*Nonlinear Dynamics*, vol. 42, no. 2, pp. 185–198, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S.-J. Ma, W. Xu, W. Li, and T. Fang, “Analysis of stochastic bifurcation and chaos in stochastic Duffing-van der Pol system via Chebyshev polynomial approximation,”
*Chinese Physics*, vol. 15, no. 6, article 017, pp. 1231–1238, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. J. Ma, W. Xu, and T. Fang, “Analysis of period-doubling bifurcation in double-well stochastic Duffing system via Laguerre polynomial approximation,”
*Nonlinear Dynamics*, vol. 52, no. 3, pp. 289–299, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S.-J. Ma, Q. Shen, and J. Hou, “Modified projective synchronization of stochastic fractional order chaotic systems with uncertain parameters,”
*Nonlinear Dynamics*, vol. 73, no. 1-2, pp. 93–100, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Liu, L. Hong, and L. Yang, “Fractional-order complex T system: bifurcations, chaos control, and synchronization,”
*Nonlinear Dynamics*, vol. 75, no. 3, pp. 589–602, 2014. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Diethelm, N. J. Ford, and A. D. Freed, “A predictor-corrector approach for the numerical solution of fractional differential equations,”
*Nonlinear Dynamics*, vol. 29, no. 1–4, pp. 3–22, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Charef, H. H. Sun, Y.-Y. Tsao, and B. Onaral, “Fractal system as represented by singularity function,”
*IEEE Transactions on Automatic Control*, vol. 37, no. 9, pp. 1465–1470, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. S. Tavazoei and M. Haeri, “Limitations of frequency domain approximation for detecting chaos in fractional order systems,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 69, no. 4, pp. 1299–1320, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. Borwein and T. Erdélyi,
*Polynomials and Polynomial Inequalities*, Springer, New York, NY, USA, 1995. View at Publisher · View at Google Scholar · View at MathSciNet