Probability density and stochastic stability for the coupled Van der Pol oscillator system

Transcription

1 STATISTICS RESEARCH ARTICLE Probability density and stochastic stability for the coupled Van der Pol oscillator system Shenghong Li, and Quanxin Zhu * Received 9 Ma7 Accepted Januar8 First Published 5 Januar8 *Corresponding author Quanxin Zhu, School of Mathematical Sciences and Institute of Finance and Statistics, Nanjing Normal university, Nanjing 3, China zqx@6.com Reviewing editor Zudi Lu, University of Southampton, UK Additional information is available at the end of the article Abstract In this paper, the probability density and almost sure asymptotic stability of the coupled Van der Pol oscillator system under the noise excitations are investigated. Through the stochastic averaging method and slow changing process theorem, averaged Fokker Planck Kolmogorov equation and exact solution of the dynamical system are obtained. Especially, the marginal density functions of the system excited by the additive white noises are derived. Then, the effects of coupled parameters and noise parameters on the marginal density functions are discussed through the numerical figures. In addition, the almost sure asymptotic stability under the parametric excitation by means of the maximal Lyapunov exponent is studied, and the stable demarcation points about noise intensity are presented. Subjects Statistics &; Probability; Probability; Stochastic Models & Processes Keywords stochastic averaging method; FPK equation; almost sure asymptotic stability; maximal Lyapunov exponent. Introduction A great deal of attention has been paid to the investigation of stochastic dynamical system excited by the noises since Itô theory was constructed (Itô, 95. Lots of stochastic methods were established, and one of which was the stochastic averaging method. The classic stochastic averaging ABOUT THE AUTHORS Quanxin Zhu is a professor of Nanjing Normal University, and he has obtained the Alexander von Humboldt Foundation of Germany. Zhu is an associate editor of MPE, JAM, TJMAA, JAAM, RAM and he is a reviewer of Mathematical Reviews and Zentralblatt-Math. Also, Zhu is a senior member of the IEEE and he is the lead guest editor of Recent Developments on the Stability and Control of Stochastic Systems in MPE. Zhu has obtained Annual Chinese and the award of one Hundred The Most Inuential International Academic Paper. Moreover, he has been one of most cited Chinese researchers in 46, Elsevier. Zhu is a reviewer of more than 5 other journals and he is the author or coauthor of more than journal papers. His research interests include random processes, stochastic control, stochastic differential equations, stochastic partial differential equations, stochastic stability, non-linear systems, Markovian jump systems, stochastic complex networks and finance mathematics. PUBLIC INTEREST STATEMENT The stochastic stability is an important topic of non-linear dynamical field. At present, we mainly focus on the study of the stochastic differential equation, stochastic dynamics and stability. Through the invariant measure, the stochastic P-bifurcation behaviours are analysed. By means of the Lyapunov exponent, Stochastic D-bifurcation, stability and stable boundaries are discussed. 8 The Author(s. This open access article is distributed under a Creative Commons Attribution (CC-BY 4. license. Page of

2 method has constituted a rigorous theoretical system and all the result has been concluded in (Lin & Cai, 995. Then, Zhu promoted the theories for the Hamilton system excited by the white noise, bounded noise and wide band noise (Zhu, Huang, & Suzuki, a, ; Zhu, Huang, & Yang, 997. Stochastic averaging method is a combination of stochastic averaging principle and Fokker Planck Kolmogorov (FPK equation. Its merits are that after applying stochastic averaging method, the FPK equation is often greatly simplified. Moreover, the dimension of obtained non-linear system is less half than that of original one, which decreases the difficulty of solving FPK equation. The exact stationary solution of FPK equation were firstly generalized to some non-linear damping systems and multiple degree of freedom non-linear systems by Caughey (Caughey & Ma, 98a, 98b; Caughey & Payne, 967, then, was extended to some general non-linear dissipated and partially integrable Hamilton systems under external and parametric excitation of Gaussian white noises (Zhu, Cai, & Lin, 99, 99; Zhu, Huang, & Suzuki, b. Lately, the stochastic averaging method was improved to apply a non-linear oscillator system (Zhao, 5, particularly, extended to the system excited by the fractional Gaussian noise excitation (Deng & Zhu, 5. But so far, the exact solution of high-dimensional FPK equation is still a very difficult problem. The almost sure asymptotic stability of a random system can be determined according to the sign of the maximal Lyapunov exponent. If the value of the maximal Lyapunov exponent is less than zero, the system is stable, otherwise, unstable (Arnold, 998. A method for exact evaluation of the maximal Lyapunov exponent in a linear Itô stochastic differential equation was first created by Khasminskii (967. Subsequently, many results for two-dimensional systems excited by the white noises or real ones were presented by some researchers (Arnold & Papanicolaou, 986; Mitchell & Kozin, 974; Nishoka, 976. By applying the stochastic averaging method, Ariaratnam and Xie (99 obtained the asymptotic expansions of maximal Lyapunov exponent for two coupled oscillators excited by a real noise. Through a perturbation method, a more general four-dimensional linear stochastic system was investigated by Doyle and Namchchivaya (994. Then, the maximal Lyapunov exponents for a co-dimension two bifurcation system under a white noise and a bounded noise excitation were given by Li and Liu ( and Liu and Liew (4. However, there have been no any results about the research of almost sure stability for above four-dimensional dynamical systems up to date. In practice, due to the role of various modes, the dynamical systems mostly present as coupled multi-degree of freedom oscillators whose dynamical behaviours are also very complex, have aroused the wide attention of many scholars. Ariaratnam and Xie (99, 994, Doyle and Namachchivaya (994, Namachchivaya, Vedula, and Onu ( and Liu and Zhu (3, 4, 5 researched this kind of the system. In this paper, the probability response and stochastic stability for coupled Van der Pol oscillator system excited by the white noises are investigated. In Section, the formulation of the problem is presented through the stochastic averaging method. The marginal probability densities of the system excited by the additive white noises are obtained, and the effects of system parameters on probability density are discussed in Section 3. In Section 4, the analytical expression of the maximal Lyapunov exponent is derived and the almost sure stability is given. The matter of the article is summarized in Section 5.. Model and formulation Consider the following coupled Van der Pol oscillator system subjected to the noise excitations. The motion equation of the system is of the form ẍ + x + ε(α + α x + α ẋ = ε(xξ (t+ξ (t ÿ + y + ε(β + β x + β ẏ = ε(yξ 3 (t+ξ 4 (t ( where ω i (i =, is an irreducible natural frequency, and α j and β j (j =,, are coupled parameters. ξ i (t(i =,, 3, 4 is independently Gaussian white noises with zero mean and the intensity D i (D i. ɛ is a very small positive parameter. Thus, Equation ( defines a quasi-hamiltonian Page of

3 system with light damping. Let = x, ẋ =, = y, ẏ =, and introduce the Wong-Zakai correction terms, thus, Equation ( is rewritten as the Itô stochastic differential equations, i.e. d = dt, d = [ ε(α + α x + α y ]dt + ε( D db (t+ D db (t, d = dt, d = [ ε(β + β x + β y ]dt + ε( db (t+ D 4 db (t, ( where B i (t(i =,, 3, 4 is a unit Wiener process. Obviously, the Hamiltonian system governed by Equation ( is quasi-periodic and integrable as ɛ = according to Zhu et al. (997, and is supposed non-resonant. Thus, the corresponding Hamiltonian function of the Equation ( is H = H + H, H =( x + x, H =(ω y + y. (3 It can be seen from the above equations that the model ( is a quasi-integrable Hamilton system. Via the following transformation z =( x + x / =( y + y /, θ = arctan( / ω, θ = arctan( / ω and according to Itô s formula, the Itô stochastic differential equations about z i, θ i (i =, are obtained (4 dz i = εμ i dt + ε (σ i db (t+σ i db (t, dθ i =( + ε θ i dt ε ( σ i db (t+ σ i db (t, i =,, (5 where μ = (α + α + α y x ω +(D x + D ω, σ = D, σ = D, μ = (β + β + β y y ω +(D 3 y + D ω, σ = 4, σ = ( D 4, θ = (α + α + α y z +(D x + D z, σ = D z, σ = D z, (6 ( θ = (β + β + β y z +(D 3 y + D 4 z, σ = z, σ = D z. In Equation (5, z i (i =, is slowly changing process and θ i (i =, is fast changing process. According to Khasminskii s theorem (968, z i (i =, weakly converges to a two-dimensional diffusion process as ε in the time interval of the order ɛ. If the two-dimensional limit process is still marked as z i (i =,, the averaged Itô stochastic differential equation for the process is in the following dz i = m i (Zdt + σ ik (ZdB i (t i =,, (7 and the corresponding FPK equation is p t = [a z i (Zp]+ [b i= i z i= j= i z ij (Zp], j where a i (Z =m i (Z, b ij (Z = σ ik (Z σ jk (Z, Z =(z T, i, j =,, and k= (8 a = α z α z α z z + D z ω + D, a = β z β z β z z + z ω + D 4 b = D z ω + D z, b = b =, b = z ω + D 4 z. (9 Page 3 of

4 3. The probability density of the coupled system Because of the complexity of Equation (8, it is almost impossible that the exact solution can be deduced. Consequently, we only consider the coupled system is additively excited by the noises. When there is only external excitation, i.e. D = =, one can calculate the exact stationary solution of Equation (8. Under the compatibility condition of α D = β D 4 = γ, the exact stationary solution of Equation (8 is derived as p(z =C exp{ φ(z }, ( where C is the normalized constant, ϕ(z is a probability potential and its expression is as follows φ(z =(4α z + α z 4D +(4β z + β z 4D 4 + γz z ( Through integrating over z from to + about Equation (, the marginal density function of z is obtained. Then substituting z =( x + x, the stationary joint probability density of and is given as { } p(, =C πβ D 4 exp (D 4 γz + β D 4 β (α z 4α ω z 4D. ( Similarly, the stationary joint probability density of and is also obtained { } p(, =C πd α exp (D γz + α D α (β z + 4β ω z 4D 4. (3 where z =( y + y. In order to intuitively observe the variety of the density function, Figures 6 are displayed. Figures give the change of the joint density function of and with the increase of parameters α and β, respectively. It can be seen from Figure that the top of joint probability density is like a crater shape, which is stochastic analogue of limit cycle and the most probability value lies on the circle. Simultaneously, both the height and the diameter of crater decline with the increase of parameter α. In Figure, the joint probability density appears the shape of mountain peak, and whose height falls with the increased β. The varieties of joint density function of and with the increase of parameters α and β are presented in Figures 3 4, respectively. Since the system ( is coupled, the curved surfaces of joint probability densities are similar to Figures. In detail, Figure 3 is alike Figures and 4 is close to Figure, which is because the status of α and β is identical in the coupled oscillator system. The effects of noise intensities on the joint probability density function are described in Figures 5 6. Figures 5 6 show that the altitude of peaks and crater raises as the intensities of noises change from D =. to D =., D 4 = Almost sure asymptotic stability When there is only parametric excitation in the system (, thus, D = D 4 =. In this section, we discuss the almost sure asymptotic stability of system (. According to the theory of stochastic dynamics, the invariant manifold of non-linear dynamical system is tangent with the invariant subspace of its corresponding linear system in equilibrium point, and so almost sure asymptotic stability of the two systems is consistent, which is central conclusion of Oseledec multiplicative ergodic theorem in the theory of stochastic dynamics. Therefore, in order to investigate the stability of coupled Van der Page 4 of

5 Figure. The joint density function of and =, =., α =.8, α =., β =.8, β = β =., D =., (a α =.5, (b α =.6, (c α =.7, (d α =.8. p(, 4 3 (a p(, (b 4 3 (c 4 3 (d p(, p(, Figure. The joint density function of and =, =., α = α =.8, α =., β =.8, β =., (a (b D =., (a β =.6, (b β =.7, (c β =.8, (d β =.9. p(,.5 p(,.5 (c (d p(,.5 p(,.5 Page 5 of

6 Figure 3. The joint density function of and =, =., α =.8, α =., β =.8, β = β =., D =., (a α =.5, (b α =.6, (c α =.7, (d α =.8. p(,.5.5 (a p(,.5 (b (c (d p(,.5 p(,.5 Pol oscillator system in the vicinity of equilibrium point, we only need to study the stability of its corresponding linear system. Consider Itô stochastic linearized differential equation of Equation (7 / dz = (α + D z dt +( D / z db (t / (4 dz = (β + z dt +( / z db (t Introduce the transformation ρ = [ ln(z + z ]/, φ i = z i /(z + z, i =, (5 Obviously, φ i [, ], i =,. The stochastic differential equations for ρ and ϕ can be obtained by Itô lemma and the expression ϕ = ϕ dρ = q ( φ dt +Σ ( φ db (t dφ = m ( φ dt + σ ( φ db (t (6 where q(φ =(c β D φ φ +[c α D ( φ 3 ω]( φ, Σ(φ = / D φ +( / ( φ, m(φ =(c φ ( φ D φ ( φ ω + φ ( φ, σ(φ =c φ ( φ, and c = β α + D, c = α β + ω, c ± = D. (7 Page 6 of

7 Figure 4. The joint density function of and =, =., α = α =.8, α =., β =.8, β =., (a.6 (b D =., (a β =.6, (b β =.7, (c β =.8, (d β =.9. p(, p(, (c (d p(,.3. p(,.3... Figure 5. The joint density function of and =, =., α = α =.8, α =., β =.8, β =., (a (b β =.6, (a D =., (b D =., D 4 =.44. p(,.5 p(,.5 It can be seen from Equation (6 the process ϕ is independent of variable ρ, and so ϕ is a diffusion process. Thus, the corresponding FPK equation of (5 is as follows p t = [m(φ φ p]+ [σ(φ φ p]. (8 There are two first kind singular boundaries ϕ =,, and no other singular point exists at < ϕ <. Via a simple calculation, we can easily obtain the diffusion index of Equation (6 is and the drift Page 7 of

8 Figure 6. The joint density function of and =, =., α = α =.8, α =., β =.8, β =., β =.8, (a D =.8, D 4 =.96, (b D =.. p(, (a p(,.5.5 (b index is. ϕ =, are both repulsively natural boundaries under the conditions c > c,c > c, thus, ϕ (t is ergodic in the interval (, and the stationary probability density of Equation (7 exists. Through the direct integral, one obtains p(φ =C φ (c c ( φ (c c, φ (, where C = Γ (4 D c Γ [(c c c ]Γ [(c + c c ]. (9 According to Khasminskii theory (Deng & Zhu, 5, the maximal Lyapunov exponent of system (4 is λ = q(φ p(φ dφ = [(c β (c c +(c α (c + c ] Γ(4 D c c Γ(c + c [c D (c c ω + c (c + c ω ] Γ(4 D c c 4 Γ((c +c c ( where c, c, c ± are shown in Equation (7. Figure 7. The plan of the maximal Lyapunov exponent and noise intensityd =, =., α., β =., = λ D = D Page 8 of

9 Figure 8. The plan of the maximal Lyapunov exponent and noise intensity =, =., α., β =., D = λ -.5 = The plan of maximal Lyapunov exponent and noise intensities D and is given in Figures 7 8. The maximal Lyapunov exponents rise with the increased D and, and are less than zero when D <.5585 and <.3997, at this moments, the system ( is almost sure asymptotic stable. Moreover, the maximal Lyapunov exponent increases more quickly in Figure 8 than in Figure 7 as noise intensities are increased. 5. Conclusion In this paper, through applying the stochastic averaging method, the averaged FPK equation of coupled Van der Pol oscillator system under the noise excitation is obtained. Moreover, its exact marginal densities are deduced under the compatibility conditions. Then, the effects of the parameters in the system ( on the marginal densities are discussed. From the Figures 6, it can be seen that the curved surfaces of density functions decline with the increase of parameters α and β, however, the influences of the noise intensities D and are contrary to them. Finally, the almost sure asymptotic stability of the system ( is researched through the maximal Lyapunov exponent, and whose analytical expression is derived. Furthermore, the stable boundary point about noise strengths D and are presented in Figures 7 8. Funding This work was jointly supported by the Alexander von Humboldt Foundation of Germany (Fellowship CHN/6339, the National Natural Science Foundation of China (63748, 698, the Natural Science Foundation of Jiangsu Province (BK655, Qing Lan Project of Jiangsu, the postdoctoral foundation of Jiangsu Province (569B, and the Priority Academic Program Development of Jiangsu Higher Education Institutions. Author details Shenghong Li, lishenhong@just.edu.cn Quanxin Zhu zqx@6.com ORCID ID http//orcid.org/795-6x School of Mathematical Sciences and Institute of Finance and Statistics, Nanjing Normal University, Nanjing 3, China. School of Science, Jiangsu University of Science and Technology, Zhenjiang 3, China. Citation information Cite this article as Probability density and stochastic stability for the coupled Van der Pol oscillator system, Shenghong Li & Quanxin Zhu, Cogent Mathematics & Statistics (8, References Ariaratnam, S. T., & Xie, W. C. (994. Almost-sure stochastic stability of coupled non-linear oscillators. International Journal of Non-Linear Mechanics, 9, https//doi.org/.6/-746( Page 9 of

10 Ariaratnam, S. T., & Xie, W. C. (99. Lyapunov exponents and stochastic stability of coupled linear systems under real noise excitation. ASME Journal of Applied Mechanics, 59, https//doi.org/.5/ Arnold, L. (998. Random dynamical systems. Berlin Springer. https//doi.org/.7/ Arnold, L., & Papanicolaou, G. (986. Asympototic analysis of the Lyapunov exponents and rotation numbers of the random oscillator and applications, SIAM. Journal of Applied Mathematics, 46, Caughey, T. K., & Ma, F. (98a. The exact steady-state solution of a class of non-linear stochastic systems. International Journal of Non-Linear Mechanics, 7, https//doi.org/.6/-746(893- Caughey, T. K., & Ma, F. (98b. The steady-state response of a class of dynamical systems to stochastic excitation. Journal of Applied Mechanics, 49, https//doi.org/.5/ Caughey, T. K., & Payne, H. J. (967. On the response of a class of self-excited oscillators to stochastic excitation. International Journal of Non-Linear Mechanics,, 5 5. https//doi.org/.6/-746(679-8 Deng, M. L., & Zhu, W. Q. (5. Stochastic averaging of quasinon-integrable Hamiltonian systems under fractional Gaussian noise excitation. Nonlinear Dynamics, 83, 5 7. Doyle, M. M., & Namachchivaya, N. S (994. Almost-sure asymptotic stability of a general four-dimensional system driven by real noise. Journal of Statistical Physics, 75, https//doi.org/.7/bf8687 Itô, K. (95. On stochastic differential equations memoris. American Mathmatical society, 4, 5. Khasminskii, R. Z. (967. Necessary and sufficient conditions for the asymptotic stability of linear stochastic systems. Theory of probability and application,, https//doi.org/.37/9 Khasminskii, R. Z. (968. On the averaging principle for Itô stochastic differential equations. Kibernetka, 3, Li, S. H., & Liu, X. B. (. The maximal Lyapunov exponent of a co-dimension two-bifurcation system excited by a bounded noise. Acta Mechanica Sinica, 8, https//doi.org/.7/s49--3-y Lin, Y. K., & Cai, G. Q. (995. Probabilistic structural dynamics. Advanced theory and applications. New York, NY McGraw-Hill. Liu, X. B., & Liew, K. M. (4. The maximal Lyapunov exponent for a three-dimension stochastic system, ASME. Journal of Applied Mechanics, 7, Liu, W. Y., & Zhu, W. Q. (3. Stochastic stability of quasi non-integrable Hamiltonian systems under parametric excitations of Gaussian and Poisson white noises. Probabilistic Engineering Mechanics, 3, Liu, W. Y., & Zhu, W. Q. (4. Stochastic stability of quasiintegrable and resonant Hamiltonian systems under parametric excitations of combined Gaussian and Poisson white noises. International Journal of Non-Linear Mechanics, 77, 5 6. https//doi.org/.6/j.ijnonlinmec Liu, W. Y., & Zhu, W. Q. (5. Lyapunov function method for analyzing stability of quasi-hamiltonian systems under combined Gaussian and Poisson white noise excitations. Nonlinear Dynamics, 8, 5. Mitchell, P. R., & Kozin, F. (974. Sample stability of second order linear differential equations with wide band noise coefficients, SIAM. Journal of Applied Mathematics, 7, Namachchivaya, N. S., Vedula, L., & Onu, K. (. Stochastic stability of coupled oscillators with one asymptotically stable and one critical mode. Journal of Applied Mechanics, 78, Nishoka, K. (976. On the stability of two-dimensional linear stochastic systems. Kodai Mathematical Seminar Reports, 7, 3. Zhao, Y. P. (5. An improved energy envelope stochastic averaging method and its application to a nonlinear oscillator. Journal of Vibration and Control,,. Zhu, W. Q., Cai, G. Q., & Lin, Y. K. (99. On exact stationary solutions of stochastically perturbed Hamiltonian systems. Probabilistic Engineering Mechanics, 5, https//doi.org/.6/66-89(99-8 Zhu, W. Q., Cai, G. Q., & Lin, Y. K. (99. Stochastically perturbed haniltonian systems nonlinear stochastic mechanics, proceedings of IUTAM symposium (pp Berlin Springer-Verlag. Zhu, W. Q., Huang, Z. L., & Suzuki, Y. (a. Response and stability of strongly nonlinear oscillators under wide band random excitation. International Journal of Non-Linear Mechanics, 36, https//doi.org/.6/s-746(93-7 Zhu, W. Q., Huang, Z. L., & Suzuki, Y. (b. Exact stationary solutions of stochastically excited and dissipated partially integrable Hamiltonian systems. International Journal of Non-Linear Mechanics, 36, https//doi.org/.6/s-746( Zhu, W. Q., Huang, Z. L., & Suzuki, Y. (. Stochastic averaging and Lyapunov exponent of quasi partially integrable Hamiltonian systems. International Journal of Non-Linear Mechanics, 37, https//doi.org/.6/s-746(8-x Zhu, W. Q., Huang, Z. L., & Yang, Y. Q. (997. Stochastic averaging of quasi-integrable Hamiltonian systems. Journal of Applied Mechanics, 64, https//doi.org/.5/ The Author(s. This open access article is distributed under a Creative Commons Attribution (CC-BY 4. license. You are free to Share copy and redistribute the material in any medium or format Adapt remix, transform, and build upon the material for any purpose, even commercially. The licensor cannot revoke these freedoms as long as you follow the license terms. Under the following terms Attribution You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use. No additional restrictions You may not apply legal terms or technological measures that legally restrict others from doing anything the license permits. Page of