# Hopf Bifurcation of Compound Stochastic van der Pol System

Authors: Shaojuan Ma; Qianling Zhang
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

Hopf bifurcation analysis for compound stochastic van der Pol system with a bound random parameter and Gaussian white noise is investigated in this paper. By the Karhunen-Loeve (K-L) expansion and the orthogonal polynomial approximation, the equivalent deterministic van der Pol system can be deduced. Based on the bifurcation theory of nonlinear deterministic system, the critical value of bifurcation parameter is obtained and the influence of random strengthδ and noise intensity σ on stochastic Hopf bifurcation in compound stochastic system is discussed. At last we found that increased δ can relocate the critical value of bifurcation parameter forward while increased σ makes it backward and the influence of δ is more sensitive than σ. The results are verified by numerical simulations.


---

## Body

## 1. Introduction

In recent years, the research of dynamics behaviors have been extensively concentrated on the stochastic nonlinear dynamical systems, whose internal uncertainty of system forming random parameter and external actions of system forming random excitation are modeled as different noise processes.There are some results about these different kinds of stochastic systems. The first kind of stochastic dynamical systems are the systems driven by stochastic parameter excitation or external noise excitation idealized as Gaussian white noise processes. Dynamics behaviors for stochastic van der Pol system and its coupled systems with parametric and external forces have been investigated by the stochastic averaging method [1–6], such as response analysis, Hopf bifurcation, and chaos. The existence of Hopf bifurcation about stochastic van der Pol system treating driving term as a randomly fluctuating quantity is obtained [7, 8]. For the nonlinear system with external Gaussian white noise, various properties of the response are presented [9–11]. With the cumulant-neglect closure method, the statistical moments in a Duffing system responding to random excitation are established [12]. The possibility for onset of chaotic motion in the Duffing oscillator under parameter excitation of bound noise is studied by the stochastic Melnikov method [13]. The second kind of stochastic systems are the systems with only internal random parameter considered as time-invariant and that served as random fields or as random variables. The orthogonal series expansion using a vector of random variables to represent the random field is widely exploited and applied in the stochastic structural systems [14–17]. Li and Fang et al. [18–21] have developed the orthogonal polynomial method to approximate the dynamics response of linear random structure by a series of orthogonal basis functions. Owing to the representation of random parameters as random fields greatly increases the size and complexity of the analysis, Ma et al. [22–25] have applied the method of orthogonal polynomial approximation to study various dynamics behaviors of nonlinear systems with random parameter idealized as random variables and have verified this method as feasible and effective.However, there also exist other stochastic systems with both random fields and external noise processes. By transforming from stochastic Duffing system to the deterministic structural Duffing system with external noise excitation, [26] has obtained the statistical moments and the probability density function of response to study the reliability analysis of structure. The hybrid polynomial chaos and pseudoexcitation method have been proposed to evaluate the random vibration response in the linear structure system with uncertain parameters under external random excitation [27]. With the combination of the orthogonal polynomial approximation and the K-L expansion, the sequential orthogonal decomposition is developed to discuss the response analysis of ground motion by Li [28, 29]. The precise integration method is used to obtain the K-L decomposition of the nonstationary filtered white noise random excitation; then based on the linear superposition principle, each element of the K-L vector is added to each equation of the extended order system by the form of one to one and a small amount of K-L vector used to compute response of the extended order system can obtain a good result [30].In this paper, we name the stochastic system with both random parameter treating as random variables and external noise processes as compound (or double) stochastic dynamical system. For random variables conforming to arch-like probability density function, a bound random parameter van der Pol system with external Gaussian white noise is considered, and the response of the system is presented by two sets of orthogonal polynomial bases in random function space under the condition of convergence in mean square. Compared with [26, 27, 30], the random response of this paper for the stochastic van der Pol system is directly expressed as the double series form of deterministic response using the known orthogonal polynomial functions by means of the sequential orthogonal decomposition. So this paper can be viewed as an extensive and supplementary version of the related literature. On the other hand, there are a few results to explore the stochastic bifurcation in this kind of stochastic system. So the stochastic Hopf bifurcation of van der Pol system with a bound random parameter and Gaussian white noise is explored in detail in this paper.This paper is organized as follows. In Section2, a new van der Pol system with a bound random parameter and Gaussian white noise is simplified into its equivalent deterministic one. The critical value of Hopf bifurcation parameter is given by the bifurcation theory of deterministic system and the effect of random strength and noise intensity on the critical value of bifurcation parameter is analyzed in Section 3. In the same section the results are verified by numerical simulations. Section 4 concludes the influence of random parameter and Gaussian white noise on Hopf bifurcation in stochastic van der Pol system.


## 2. The Equivalent Deterministic One for Compound Stochastic van der Pol System

Excited by internal random parameter and Gaussian white noise, compound stochastic van der Pol system is established in this paper:(1)
x
¨
+
α
x
2
-
β
¯
x
˙
+
x
=
σ
F
t
,where α and β
¯ are deterministic parameter and random parameter, respectively. F
(
t
) is Gaussian white noise with zero mean and its intensity is σ.Random parameterβ
¯ can be defined as(2)
β
¯
=
β
+
δ
u
,where u is random variable obeying arch-like probability density function defined on [
-
1,1
] and β and δ are all the statistical parameters of β
-.Using eigenfunctions as the orthogonal bases, stochastic noise processF
(
t
) can be expressed as the denumerable series form called the K-L expansion [31–33], so the K-L expansion of F
(
t
) is (3)
F
t
=
∑
n
=
0
N
ξ
n
λ
n
f
n
t
,where λ
n and f
n
(
t
) are the eigenvalues and eigenfunctions of the covariance function of F
(
t
), and eigenfunctions f
n
(
t
) form a set of complete orthogonal bases in complete Hilbert space and satisfy ∫
0
∞
f
n
(
t
)
f
m
(
t
)
d
t
=
δ
n
m, in which δ
n
m is the Kronecker function. ξ
n are independent orthogonal random variables obeying normal distribution and have the properties of E
(
ξ
n
)
=
0, E
(
ξ
n
2
)
=
1.Eigenvalues and eigenfunctions of the K-L expansion for Gaussian white noise with zero mean are given [34], namely, λ
n
=
1
/
[
(
n
+
1
/
2
)
2
π
2
] and f
n
(
t
)
=
2
(
(
n
+
1
/
2
)
π
t
).The only fixed point of system (1) is A
σ
F
(
t
)
,
0 when σ
≠
0 and can be further transformed into the original point O
(
0,0
). Through the coordinate transformation, system (1) is rewritten as coupled first-order system (here is still noted by x):(4)
x
˙
=
y
,
y
˙
=
-
x
+
β
¯
-
α
σ
2
F
2
t
y
-
α
x
+
2
σ
F
t
x
y
.By the sequential orthogonal decomposition, the response of system (4) can be expressed in the following double series form:(5)
x
ξ
,
u
,
t
=
∑
i
=
0
N
∑
j
=
0
M
H
i
ξ
T
j
u
x
i
j
t
,
y
ξ
,
u
,
t
=
∑
i
=
0
N
∑
j
=
0
M
H
i
ξ
T
j
u
y
i
j
t
,where x
i
j
(
t
)
=
∫
Ω
ξ
∫
Ω
u
x
ξ
,
u
,
t
H
i
ξ
T
j
u
ω
ξ
θ
u
d
ξ
d
u and y
i
j
(
t
)
=
∫
Ω
ξ
∫
Ω
u
y
(
ξ
,
u
,
t
)
H
i
(
ξ
)
T
j
(
u
)
ω
(
ξ
)
θ
(
u
)
d
ξ
d
u, in which ω
(
ξ
) and θ
(
u
) are normal probability density function and arch-like probability density function, respectively. Corresponding to these probability density functions, Hermite orthogonal polynomial H
i
(
ξ
) and Chebyshev orthogonal polynomial of the second kind T
j
(
u
) can be taken as orthogonal bases. N and M are the truncation order of the series we have taken. In addition, it is assumed that random variables ξ and u are mutually independent.Substituting (2), (3), and (5) into (4), equivalent deterministic van der Pol system can be given when N
→
∞ and M
→
∞. So as to make the following numerical simulations effective and convenient, we take N
=
M
=
1 and approximately obtain the equivalent deterministic system of compound stochastic van der Pol system:(6)
x
˙
00
t
=
y
00
t
,
y
˙
00
t
=
-
x
00
t
+
β
y
00
t
+
δ
2
y
01
t
-
α
M
00
t
-
α
σ
2
λ
0
+
λ
1
y
00
t
,
x
.
01
t
=
y
01
t
,
y
˙
01
t
=
-
x
01
t
+
β
y
01
t
+
δ
2
y
00
t
-
α
M
01
t
-
α
σ
2
λ
0
+
λ
1
y
01
t
,
x
˙
10
t
=
y
10
t
,
y
˙
10
t
=
-
x
10
t
+
β
y
10
t
+
δ
2
y
11
t
-
α
M
10
t
-
α
σ
2
λ
0
+
λ
1
y
10
t
,
x
˙
11
t
=
y
11
t
,
y
˙
11
t
=
-
x
11
t
+
β
y
11
t
+
δ
2
y
10
t
-
α
M
11
t
-
α
σ
2
λ
0
+
λ
1
y
11
t
,where M
i
j
(
t
)
(
i
=
0,1
,
j
=
0,1
) notes the coefficient of H
i
ξ
T
j
u in the linear combination.The deterministic parameter for system (6) can be taken as α
=
1.0 and initial values are (7)
x
00
0
=
0.3
,
y
00
0
=
0.2
,
x
i
j
t
=
y
i
j
t
=
0
i
=
0,1
,
j
=
0,1
.


## 3. Bifurcation Parameter and Numerical Simulations

According to the Hopf bifurcation theory of deterministic nonlinear system [35, 36], the critical value of bifurcation parameter for the equivalent deterministic system (6) is obtained as follows:(8)
β
c
=
80
σ
2
-
9
π
2
δ
18
π
2
.By the following numerical simulations, the influence of random strengthδ and noise intensity σ on bifurcation parameter β
c is discussed in this section. Figure 1 clearly shows the theoretical relationship among δ, σ, and β
c. β
c decreases as δ increases while β
c increases as σ increases. In other words, δ can induce the generation of Hopf bifurcation but σ delays it.Figure 1
The function diagram between the bifurcation parameter and two kinds of random intensity.System (4) is the deterministic van der Pol system when we take δ
=
σ
=
0. Meanwhile, the bifurcation parameter β equals to β
c
=
0 and the system shows limit cycles at the fixed point O
(
0,0
). Hopf bifurcation phase diagrams and time history diagrams are plotted in Figure 2 when we, respectively, take β less than and greater than β
c.Figure 2
The Hopf bifurcation phase diagrams and time history diagrams of deterministic van der Pol system.
(a)
β
<
β
c
(b)
β
>
β
c
(c)
β
<
β
c
(d)
β
>
β
cThe stochastic van der Pol system with a bound random parameter is the case ofδ
≠
0 and σ
=
0. Random variable u continually changes for the random response; we can take its ensemble mean response to replace it. Taking δ
=
0.1, we obtain β
c
=
-
0.05. When β
c
>
β
=
-
0.11, it takes more time to gradually converge at the fixed point for the trajectory of ensemble mean response than the deterministic system shown in Figures 3(a) and 3(c). The trajectory of ensemble mean response is gradually convergent closed curve when β
c
<
β
=
-
0.045 and time used for converging is shorter than the deterministic one shown in Figures 3(b) and 3(d). Compared with the deterministic van der Pol system, the critical value of bifurcation parameter is smaller and Hopf bifurcation occurred in advance.Figure 3
The bifurcation phase diagrams and time history diagrams of a bound random parameter van der Pol system.
(a)
δ
=
0.1, β
<
β
c
(b)
δ
=
0.1, β
>
β
c
(c)
δ
=
0.1, β
<
β
c
(d)
δ
=
0.1, β
>
β
cIt is also stochastic van der Pol system withδ
=
0 and σ
≠
0. We take the ensemble mean response about random variable ξ to substitute the random response. Taking σ
=
0.1, we obtain β
c
=
0.0045. When β
c
>
β
=
-
0.085, the trajectory is gradually convergent at the fixed point described in Figures 4(a) and 4(c). The trajectory is gradually convergent closed curve when β
c
<
β
=
0.005 also described in Figures 4(b) and 4(d). Compared with the deterministic van der Pol system, the noise intensity σ of stochastic van der Pol system makes the critical value of bifurcation parameter larger and Hopf bifurcation occurred later.Figure 4
The bifurcation phase diagrams and time history diagrams of van der Pol system with Gaussian white noise.
(a)
σ
=
0.1 and β
<
β
c
(b)
σ
=
0.1 and β
>
β
c
(c)
σ
=
0.1 and β
<
β
c
(d)
σ
=
0.1 and β
>
β
cThe compound stochastic van der Pol system isδ
≠
0 and σ
≠
0. We still take the ensemble mean response about random variables u and ξ to stand for the random response (5) and give the bifurcation phase diagram and time history diagram when δ
=
σ
=
0.1 and β
<
β
c (see Figure 5). Comparing Figure 5(b) with Figure 2(c), it is found that increased δ postpones time of convergence and promotes the generation of Hopf bifurcation at the fixed point, and increased σ shortens time of convergence and postpones the generation of Hopf bifurcation. Meanwhile, it is obvious that the range for the larger value of the ensemble mean response shown in Figure 5(b) is the same as the result shown in Figure 3(c) and is larger than the result described in Figure 4(c).Figure 5
The bifurcation phase diagrams and time history diagrams of compound stochastic van der Pol system.
(a)
δ
=
σ
=
0.1 and β
<
β
c
(b)
δ
=
σ
=
0.1 and β
<
β
cFrom the above numerical simulations, it is discovered thatδ and σ have the different influence on the critical value of Hopf bifurcation parameter in the compound stochastic van der Pol system. As we increase the random parameter strength δ, the critical value of bifurcation parameter is smaller than the deterministic system. While we increase the noise intensity σ, the critical value of bifurcation parameter is larger than the deterministic system. Simultaneously, the influence of δ on the critical value of bifurcation parameter is more sensitive than σ taking the same value.


## 4. Conclusion

In this paper, we deal with Hopf bifurcation of compound stochastic van der Pol system with random parameter and Gaussian white noise. It shows that the method of sequential orthogonal decomposition combining the orthogonal polynomial approximation with the K-L expansion is effective to deduce the stochastic van der Pol system into its equivalent deterministic system; then the theory of Hopf bifurcation in nonlinear deterministic system is applied to the equivalent one and thereby Hopf bifurcation of stochastic van der Pol system under random parameter and Gaussian white noise is analyzed. In addition, it is found that these two kinds of random intensity have the different effects on the critical value of Hopf bifurcation parameter for compound stochastic van der Pol system and the influence ofδ is more outstanding than σ in the same value. Finally, we verify the results by numerical simulations.The compound stochastic van der Pol system established in this paper is also available for other nonlinear stochastic dynamical systems. And various dynamics behaviors in such systems as an extensive and supplementary version of the related literature should have been studied deeply, such as different bifurcations, stability, chaotic motion, and numerical analysis of response. It is especially worth taking into account the convergence study of response prediction because of a series with finite terms.



---

*Source: 4246916-2016-06-06.xml*