# Bursting Oscillations in Shimizu-Morioka System with Slow-Varying Periodic Excitation

Authors: Xindong Ma; Shuqian Cao
Journal: Shock and Vibration
 (2018)

---

## Abstract

The coupling effect of two different frequency scales between the exciting frequency and the natural frequency of the Shimizu-Morioka system with slow-varying periodic excitation is investigated. First, based on the analysis of the equilibrium states, homoclinic bifurcation, fold bifurcation, and supercritical Hopf bifurcation are observed in the system under a certain parameter condition. When the exciting frequency is much smaller than the natural frequency, we can regard the periodic excitation as a slow-varying parameter. Second, complicated dynamic behaviors are analyzed when the slow-varying parameter passes through different bifurcation points, of which the mechanisms of four different bursting patterns, namely, symmetric “homoclinic/homoclinic” bursting oscillation, symmetric “fold/Hopf” bursting oscillation, symmetric “fold/fold” bursting oscillation, and symmetric “Hopf/Hopf” bursting oscillation via “fold/fold” hysteresis loop, are revealed with different values of the parameterb by means of the transformed phase portrait. Finally, we can find that the time interval between two symmetric adjacent spikes of bursting oscillations exhibits dependency on the periodic excitation frequency.


---

## Body

## 1. Introduction

Multiple time scales problems can be observed in many real systems, such as catalytic reactions in chemical systems [1, 2], electrophysiological experiments [3–6], and circuit systems [7, 8]. The coupling effect of these different fast-slow time scales generally results in a system that exhibits periodic motion characterized by a combination of relatively large amplitude and nearly harmonic small amplitude oscillations [9, 10], conventionally denoted by NK with N and K corresponding to the numbers of the large and small amplitude oscillations, respectively [11]. Variables that exhibit large amplitude oscillations generate a spiking state (SP) [12], whereas the systems may be in a quiescent state (QS) when all the variables exhibit small amplitude oscillations [13]. The dynamic behaviors connecting the spiking states and the quiescent states are called bursting oscillations. There are two types of bifurcations existing in one bursting oscillation, specifically one that transforms the system from spiking states to quiescent states and the other one that transforms the system from quiescent states to spiking states [14, 15]. Therefore, the bursting patterns not only depend on the patterns of the spiking states and quiescent states, but also depend on the patterns of the bifurcations that connect the two states [16].For a typical system with different coupled time scales, its dynamic behaviors can be described by a singularly perturbed system with two time scales of the following form [17]:(1)x˙=fx,uFastSpiking,u˙=εgx,uSlowModulation,where ε is a small parameter that satisfies 0<ε≪1 and represents different time scales. Vector x∈Rm models the dynamics of a relatively fast changing process, while vector u∈RL describes the relatively slowly changing quantity that modulates x. Characterization of the dynamic behaviors at different time scales has long intrigued scholars. Rinzel [18] introduced the fast-slow analysis method to examine the Chay-Keizer model and characterize dynamic behaviors at different time scales and the mechanism of bursting oscillations. Wang et al. [19] examined the synchronization of bursting through the introduction of finite delays. Kiss et al. [20] studied the bursting oscillations in electrochemical systems caused by changes in the geometrical structure. Li et al. [21] examined the mechanism of bursting oscillations in a nonsmooth generalized Chua’s circuit with two time scales. However, most of the current work focused on the single slow variable, the mechanism of which is relatively simple, and some work is concentrated on the autonomous system. However, the mechanism of bursting oscillations in nonautonomous systems with slow-varying periodic excitation requires further exploration.In the present work, we consider the Shimizu-Morioka system [22] with slow-varying excitation (SMSWSVE):(2)x˙=yy˙=x-ay-xz+Asin⁡ωtz˙=-bz+x2,where a∈R+, b∈R+, A is the forcing amplitude, and ω is the forcing frequency, which satisfies 0<ω≪1. When the forcing frequency ω is much smaller than the natural frequency ω0, we characterize the presented system as a typical fast-slow system involving two time scales. Inspired by Rinzel’s method, we present an analysis of bursting oscillation driven by slow-varying periodic excitation for the Shimizu-Morioka system. The purpose of this work is twofold. First, some results about bursting oscillations induced by the slow-varying periodic excitation are derived and how the excitation modulates the dynamics of bursting oscillations is drawn. Second, some important aspects of Rinzel’s method are highlighted in the forced Shimizu-Morioka system.The rest of this paper is organized as follows. In Section2, we summarize some of the results related to the stabilities and bifurcations of SMSWSVE, and based on this, we plot a two-parameter bifurcation set and obtain three SMSWSVE bifurcation diagrams following changes in parameter b, which exhibit typical dynamic evolution behaviors of SMSWSVE. In Section 3, we investigate the generation of bursting oscillations based on three cases of different values of the parameter b, of which four bursting patterns are obtained: symmetric “homoclinic/homoclinic” bursting oscillation for b=0.65; symmetric “fold/Hopf” bursting oscillation for b=1.2; symmetric “fold/fold” bursting oscillation; and symmetric “Hopf/Hopf” bursting oscillation via “fold/fold” hysteresis loop for b=2. In Section 4, we focus on the effects of the excitation amplitude and frequency on the bursting oscillation. Finally, Section 5 concludes the paper.


## 2. Stabilities and Bifurcations of the SMSWSVE

In order to analyze the influence of the slowly varying excitation, we can regard the external excitationAsin⁡(ωt) as a slow-varying parameter c. When the exciting frequency ω is much smaller than the natural frequency ω0, during an arbitrary period T, (T=2π/ω0), that is, t∈(t0,t0+T), where t0 is the initial timing point, the slow-varying parameter c may change between cA=Asin⁡ω0t0 and cB=Asin⁡ω0(t0+T), thereby implying that cA≈cB, such that the parameter c is kept almost constant during any arbitrary period T. Therefore, (2) can be considered as a generalized autonomous system, the equilibrium points of which can also be considered as generalized equilibrium points. The generalized autonomous system in (2) forms the fast subsystem, whereas c˙=Aωcos⁡ωt forms the slow subsystem.The equilibrium point of (2) can be expressed in the form of EQ(x,y,z)=(x0,0,x02/b) with respect to the periodic external excitation Asin⁡ωt as the parameter c, where x0 satisfies(3)x03-bx0-bc=0,the stability of which can be determined by the corresponding characteristic equation, which is defined as (4)λ3+a+bλ2+ab+x02b-1λ+3x02-b=0.Obviously, EQ is stable for a+b>0, 3x02-b>0, and (a+b)(ab+x02/b-1)-(3x02-b)>0. When the stability conditions are changed, different types of bifurcations may occur as follows.Fold Bifurcation. For 3x02-b=0, the three equilibrium points may meet together to form one degenerate equilibrium point, thereby generating eigenvalues for the characteristic equation that can be expressed as λ1=0, and λ2,3=-(a+b)±(a+b)2-4(ab+x02/b-1)/2, which imply the presence of fold bifurcation.Hopf Bifurcation. (a+b)(ab+x02/b-1)-(3x02-b)=0 defines the presence of a pair of pure imaginary eigenvalues, which implies the occurrence of supercritical Hopf bifurcation and thus results in periodic oscillation with a frequency of ΩH=ab+x02/b-1.Here we fix the parametera=1 and regard the parameters b and c as the control parameters. The two bifurcation sets of SMSWSVE are plotted in Figure 1, wherein the parameter plane (c,b) is divided into eight regions. According to our numerical calculations, region 8 exhibits an unstable equilibrium point E0 and two stable equilibria E±. When the parameters pass through the fold bifurcation curves LP1 and LP2 and enter regions 7 and 6, the numbers of the equilibria may decrease via fold bifurcation, thereby leaving two stable equilibria E+ in region 7 and E- in region 6. As the parameters pass through the Hopf bifurcation curves supH1 and supH2 and enter regions 5 and 4, the two stable equilibria, specifically E+ and E-, lose their stabilities via Hopf bifurcations and are left with two stable limit cycles, namely, LC1 in region 5 and LC2 in region 4. When the parameters pass through the fold bifurcation curves LP1 and LP2 and enter regions 3 and 2, the two stable limit cycles LC1 and LC2 are still observed, and the numbers of the equilibria may increase via fold bifurcation, specifically one stable equilibrium point and two unstable equilibrium points. As the parameters pass through the Hopf bifurcation curves supH1 and supH2 and enter region 1, the stable equilibrium point loses its stability, thereby producing one stable limit cycle via Hopf bifurcation. Therefore, three unstable equilibria and two stable limit cycles are observed in region 1.Figure 1
Two bifurcation sets of the system with respect to parametersb and c, for fixed a=1. BT1 and BT2 represent the two Bogdanov-Takens bifurcation points; LP1 and LP2 represent the twofold bifurcation curves; and supH1 and supH2 represent the two Hopf bifurcation curves.From the two bifurcation sets of SMSWSVE with respect to the parametersb and c, we can find that when the parameter b is taken at different values, the two types of curves corresponding to fold bifurcation and Hopf bifurcation exhibit different positions. For example, when the parameter b is smaller than 0667, no Hopf bifurcation is observed. When the parameter b is greater than 0667 but smaller than 1.420, the supH1 curve is observed on the left side of the LP1 curve and the supH2 curve is observed on the right side of the LP2 curve, thereby indicating that a sufficiently large amplitude A of the parameter c generates an increase (decrease) in the parameter c from point 0, wherein it first passes through the Hopf bifurcation points and then passes through the fold bifurcation points. When the parameter b is greater than 1.420, the supH1 curve is observed on the right side of the LP1 curve and the supH2 curve is observed on the left side of the LP2 curve, thereby indicating that the absolute value of the Hopf bifurcation points is greater than the absolute value of the fold bifurcation points.To comprehensively investigate the influence of the slow-varying periodic excitation, three cases with different values of the parameterb, that is, case A: b=0.65, case B: b=1.2, and case C: b=2.0, are considered, the bifurcations of which are presented in Figure 2, respectively.Figure 2
Bifurcation diagrams for (a)b=0.65; (b) b=1.2; and (c) b=2. The solid line represents the stable parts; the dashed line represents the unstable parts; ● represents the stable limit cycles; Hom1 and Hom2 represent the two homoclinic bifurcation points; LP1 and LP2 represent the twofold bifurcation points; supH1 and supH2 represent the two Hopf bifurcation points; LC1 and LC2 represent the two stable limit cycles; and E+, E-, and E0 represent the equilibrium branch of the generalized autonomous system.
(a) (b) (c)Forb=0.65 in Figure 2(a), twofold bifurcation points are observed at LP1 (0.310,-0.466) and LP2 (-0.310,0.466). And two homoclinic bifurcation points are observed at Hom1 (0.029,-0.034) and Hom2 (-0.029,0.034), at which the two stable limit cycles can be switched.Forb=1.2 in Figure 2(b), twofold bifurcation points are observed at LP1 (0.420,-0.640) and LP2 (-0.420,0.640), whereas two supercritical Hopf bifurcation points are observed at supH1 (0.203,1.187) and supH2 (-0.203,-1.187), at which the two limit cycles are produced via Hopf bifurcation.Forb=2.0 in Figure 2(c), twofold bifurcation points are observed at LP1 (0.539,-0.806) and LP2 (-0.539,0.806), at which the jumping phenomenon may occur. In addition, two supercritical Hopf bifurcation points are observed at supH1 (1.222,1.834) and supH2 (-1.222,-1.834), at which the two limit cycles are produced via Hopf bifurcation.Sincec=Asin⁡ωt, the change interval of the parameter c depends on its amplitude A. Therefore, the influence of the amplitude A on the system can be examined to characterize the dynamic behaviors of SMSWSVE when c passes through different bifurcation points.


## 3. Dynamical Mechanisms of the Bursting Oscillations

The present section characterizes the dynamical behaviors of (2) for the three typical cases. In this article, SMSWSVE is a typical fast-slow system involving two time scales, which may generate bursting oscillations in the system. In the following analysis, the excitation frequency is fixed at ω=0.005.


### 3.1. Symmetric “Homoclinic/Homoclinic” Bursting in Case A

In case A forb=0.65, no stable equilibrium points are observed. However, one stable limit cycle LC1 is observed for c>-0.029, one stable limit cycle LC2 is observed for c<0.029, and two stable limit cycles are observed for -0.029<c<0.029 simultaneously. The stable limit cycle LC1 disappears via homoclinic bifurcation of the saddle point at Hom2 (-0.029,0.034), and the limit cycle LC2 disappears via homoclinic bifurcation of the saddle point at Hom1 (0.029,-0.034). The two stable limit cycles can be switched via two homoclinic bifurcation points, namely, Hom1 and Hom2.Figure3 presents the bursting oscillation for A=0.3, wherein the trajectory moves around two stable limit cycles, namely, LC1 and LC2, as presented in Figure 3(a), which are switched by the two homoclinic bifurcation points. In order to reveal the mechanism of this bursting oscillation, we turn to the help of the transformed phase portrait on the (c,x) plane, as presented in Figure 3(b). The trajectory, starting from the point A, moves along the attraction domain of the stable limit cycle LC1, and the movement exhibits a repetitive spiking state with large amplitude. As the parameter c decreases to −0.029, where the homoclinic bifurcation point Hom2 is located, LC1 disappears via homoclinic bifurcation. In addition, the trajectory jumps into the vector field of stable limit cycle LC2. Parameter c continues to decrease until it reaches its minimum value of c=-0.3, after which it gradually increases.Figure 3
Bursting oscillation forA=0.3. (a) Phase portrait on the (x,y) plane; (b) transformed phase portrait on the (c,x) plane with equilibrium branches; (c) time history.
(a) (b) (c)When the parameterc increases to 0.029, where the homoclinic bifurcation point Hom1 is located, LC2 disappears via saddle homoclinic bifurcation, and the trajectory jumps into the vector field of the stable limit cycle LC1. Parameter c continues to increase until it reaches its maximum value of c=0.3, after which it gradually decreases. The return of the trajectory to starting point A finished one periodic movement.Bursting oscillation is observed obviously in the time history in Figure3(c), wherein the up-state and the down-state are switched via the two homoclinic bifurcation points that form the bursting oscillation with the symmetric structure. Therefore, the bursting can be classified as the symmetric “homoclinic/homoclinic” type.


### 3.2. Symmetric “Fold/Hopf” Bursting in Case B

In case B forb=1.2, twofold bifurcation points LP1,2(±0.420,±0.640) and two Hopf bifurcation points supH1,2(±0.203,±1.187) are observed on the (c,x) plane, which result in the transition between the repetitive spiking and quiescent states.When the amplitudeA is taken at A=0.430, the slow variable c exhibits changes between −0.430 and 0.430, that is, c∈[-0.430,0.430], thereby implying the presence of the fold bifurcation points at c=±0.423 and the Hopf bifurcation points at c=±0.203. From the phase portrait on the (x,y) plane plotted in Figure 4(a), the attractors oscillate around two equilibrium points, as denoted by E±. In addition, a typical bursting oscillation is observed in the time history in Figure 4(b). The trajectory jumps between two spiking states SP±, which are symmetric to each other.Figure 4
Bursting oscillation forA=0.430. (a) Phase portrait on the (x,y) plane; (b) time history; (c) transformed phase portrait on the (c,x) plane with the bifurcation diagram.
(a) (b) (c)In order to reveal the mechanism of the bursting oscillation, we turn to the transformed portrait phase with the bifurcation diagram, as presented in Figure4(c). The trajectory starts from the point A with c=0.430, where c is taken at its positive extreme value. At this point, the attractor is categorized as the stable limit cycle LC1, thereby allowing the trajectory to move within the vector field of LC1, and repetitive spiking movement with large amplitude is observed. Parameter c then decreases, as indicated by the arrow. When the parameter c decreases to c=0.203, where the Hopf bifurcation point supH1 is observed, thereby resulting in the disappearance of stable limit cycle LC1 and the appearance of stable equilibrium point E+.Due to the slow passage effect [23], the trajectory will not move within the vector field of E+ immediately, it still moves in the vector field of LC1 for a short time, and then it turns into the vector field of E+. The repetitive spiking state movement is terminated, and the quiescent state movement is created. Parameter c then continues to decrease to c=-0.420, where the fold bifurcation point LP2 is located, thereby generating a trajectory jump into the vector field of LC2 and resulting in the termination of the quiescent state movement and the creation of repetitive spiking state movement. As the parameter c decreases to c=-0.430, it reaches its negative extreme value, after which the parameter c exhibits a gradual increase.When the parameterc increases to c=-0.203, where the other Hopf bifurcation point supH2 is observed, it results in the disappearance of the stable limit cycle LC2 and the appearance of the stable equilibrium point E-. The trajectory still moves in the vector filed of LC2 for a short time according to the slow passage effect, after which it turns into the vector field of E-. The Hopf bifurcation point supH2 terminates the repetitive spiking state movement and exhibits the quiescent state movement. When the parameter c increases to c=0.420, where the fold bifurcation point LP1 is observed, it results in the disappearance of the quiescent state movement and the appearance of the repetitive spiking state movement. In addition, this increase results in the loss of E-’s stability and the generation of the stable limit cycle LC1. Parameter c then continues to increase to c=0.430. The return of the trajectory back to the starting point A finished one periodic movement.The bursting oscillation can be classified as symmetric “fold/Hopf” bursting, since the repetitive spiking related to the limit cycle is generated by fold bifurcation and disappears by Hopf bifurcation.


### 3.3. Bursting Oscillation in Case C

In case C forb=2, twofold bifurcation points LP1,2(±0.539,±0.806) and two Hopf bifurcation points supH1,2(±1.222,±1.834) are observed on the (c,x) plane, resulting in the transition between the repetitive spiking and quiescent states. Two different types of bursting oscillations are observed when the amplitude A is taken at different values.


#### 3.3.1. Symmetric “Fold/Fold” Bursting Oscillation

When the amplitudeA is taken at A=1.0, the slow variable c changes between −1.0 and 1.0, thereby implying the presence of only twofold bifurcation points at c=±0.539. Under this condition, three equilibria are present: two stable foci denoted by E± and one unstable saddle point denoted by E0, as presented in Figure 5.Figure 5
Distribution of the equilibria forA=1.0.From the phase portrait on the(x,y) plane, as presented in Figure 6(a), the trajectory oscillates around two equilibrium points, namely, E±. A typical bursting oscillation is observed in the time history diagram in Figure 6(b). In addition, the trajectory jumps between two spiking states, specifically SP±, which are symmetrical to each other.Figure 6
Bursting oscillation forA=1.0. (a) Phase portrait on the (x,y) plane; (b) time history; (c) transformed phase portrait on the (c,x) plane with equilibrium branches.
(a) (b) (c)In order to reveal the mechanism of the bursting oscillation, we turn to the transformed phase portrait on the(c,x) plane with equilibrium branches, as presented in Figure 6(c). The trajectory starts from the point A with c=1.0, where the parameter c is taken at its positive extreme value. The trajectory moves within the vector field of the stable focus E+ that is associated with the quiescent state. When the parameter c decreases to c=-0.539, where the fold bifurcation point LP2 is observed, thereby resulting in a trajectory jump to the vector field of the other stable focus E- and ending the quiescent state. The trajectory exhibits high-frequency oscillation with a gradual decrease in the amplitude, which corresponds to the spiking state. As the parameter c decreases to its negative value, specifically c=-1, it may generate a gradual increase that corresponds to the quiescent state.When the parameterc increases to c=0.539, where the fold bifurcation point LP1 is observed, it results in a trajectory jump to the vector field of E+ and terminating the quiescent state. In addition, the trajectory is characterized as high-frequency oscillation with a gradual decreasing in the amplitude, which corresponds to the spiking state. As the parameter c increases to c=1, the trajectory returned back to its starting point A, thereby finishing one periodic movement.Bursting oscillation can be created since the system periodically switches between the two stable fociE+ and E- by twofold bifurcations and the two parts are symmetrical. So the bursting pattern can be classified as symmetric “fold-fold” type.


#### 3.3.2. Symmetric “SupHopf/SupHopf” Bursting Oscillation via “Fold/Fold” Hysteresis Loop

A continues to increase to the range of A>1.222, for example, A=2, thereby implying that the parameter c not only passes through the fold bifurcation points LP1 and LP2 at c=±0.539, but also passes through the supercritical Hopf bifurcation points supH1 and supH2 at c=±1.222, where the Hopf bifurcation occurs.According to the phase portrait on the(x,y) plane in Figure 7(a), the trajectory oscillates around two stable foci, namely, E±. In addition, a typical bursting oscillation is observed in the time history diagram in Figure 7(b), and the trajectory jumps between two spiking states, specifically SP±, both of which are symmetrical to each other. However, a great difference is observed between Figure 7(a) and Figure 6(a) because two limit cycles are observed in Figure 7(a).Figure 7
Bursting oscillation forA=2. (a) Phase portrait on the (x,y) plane; (b) time history; (c) overlap of the transformed phase portrait on the (c,x) plane with the system bifurcation diagram.
(a) (b) (c)In order to explore the mechanism of the bursting oscillation in Figure7, we turn to the overlap of the transformed phase portrait on the (c,x) plane with the system bifurcation diagram, as presented in Figure 7(c). We assume that the trajectory, starting from the point A with c=2, moves within the vector field of the limit cycle LC1, corresponding to the spiking state. When the parameter c decreases to c=1.222, where the Hopf bifurcation point supH1 is observed, thereby resulting in the disappearance of limit cycle LC1 and the presence of the stable focus E+. In addition, the Hopf bifurcation point supH1 resulted in the termination of the spiking state and the appearance of the quiescent state. As the parameter c decreases to c=-0.539, where the fold bifurcation point LP2 is observed, the trajectory jumps into the vector field of the stable focus E-, and the fold bifurcation point LP2 also results in the termination of the quiescent state and the appearance of the spiking state. When the parameter c decreases to c=-1.222, E- loses its stability via Hopf bifurcation, and a stable limit cycle LC2 is observed. As the parameter c decreases to an negative extreme value, specifically c=-2, it exhibits a gradual increase.When the parameterc increases to c=-1.222, LC2 disappears and E- is observed due to Hopf bifurcation. In addition, this increase results in the termination of the spiking state and the appearance of the quiescent state. As the parameter c increases to c=0.539, where the fold bifurcation point LP1 is observed, which results in the trajectory jumps to the vector field of E+, thereby terminating the quiescent state and generating the presence of the spiking state. When the parameter c exhibits an increase to c=1.222, where the Hopf bifurcation point supH1 is observed, E+ loses its stability and the stable limit cycle LC1 is produced via Hopf bifurcation. As the parameter c increases to c=2, the trajectory returns back to its starting point A, thereby finishing one periodic movement.From Figure7(c), a hysteresis loop produced by the equilibrium branches between twofold bifurcation points, specifically LP1,2, is observed. As a result, this bursting is classified as symmetric “Hopf/Hopf” type via the “fold/fold” hysteresis loop.


## 3.1. Symmetric “Homoclinic/Homoclinic” Bursting in Case A

In case A forb=0.65, no stable equilibrium points are observed. However, one stable limit cycle LC1 is observed for c>-0.029, one stable limit cycle LC2 is observed for c<0.029, and two stable limit cycles are observed for -0.029<c<0.029 simultaneously. The stable limit cycle LC1 disappears via homoclinic bifurcation of the saddle point at Hom2 (-0.029,0.034), and the limit cycle LC2 disappears via homoclinic bifurcation of the saddle point at Hom1 (0.029,-0.034). The two stable limit cycles can be switched via two homoclinic bifurcation points, namely, Hom1 and Hom2.Figure3 presents the bursting oscillation for A=0.3, wherein the trajectory moves around two stable limit cycles, namely, LC1 and LC2, as presented in Figure 3(a), which are switched by the two homoclinic bifurcation points. In order to reveal the mechanism of this bursting oscillation, we turn to the help of the transformed phase portrait on the (c,x) plane, as presented in Figure 3(b). The trajectory, starting from the point A, moves along the attraction domain of the stable limit cycle LC1, and the movement exhibits a repetitive spiking state with large amplitude. As the parameter c decreases to −0.029, where the homoclinic bifurcation point Hom2 is located, LC1 disappears via homoclinic bifurcation. In addition, the trajectory jumps into the vector field of stable limit cycle LC2. Parameter c continues to decrease until it reaches its minimum value of c=-0.3, after which it gradually increases.Figure 3
Bursting oscillation forA=0.3. (a) Phase portrait on the (x,y) plane; (b) transformed phase portrait on the (c,x) plane with equilibrium branches; (c) time history.
(a) (b) (c)When the parameterc increases to 0.029, where the homoclinic bifurcation point Hom1 is located, LC2 disappears via saddle homoclinic bifurcation, and the trajectory jumps into the vector field of the stable limit cycle LC1. Parameter c continues to increase until it reaches its maximum value of c=0.3, after which it gradually decreases. The return of the trajectory to starting point A finished one periodic movement.Bursting oscillation is observed obviously in the time history in Figure3(c), wherein the up-state and the down-state are switched via the two homoclinic bifurcation points that form the bursting oscillation with the symmetric structure. Therefore, the bursting can be classified as the symmetric “homoclinic/homoclinic” type.


## 3.2. Symmetric “Fold/Hopf” Bursting in Case B

In case B forb=1.2, twofold bifurcation points LP1,2(±0.420,±0.640) and two Hopf bifurcation points supH1,2(±0.203,±1.187) are observed on the (c,x) plane, which result in the transition between the repetitive spiking and quiescent states.When the amplitudeA is taken at A=0.430, the slow variable c exhibits changes between −0.430 and 0.430, that is, c∈[-0.430,0.430], thereby implying the presence of the fold bifurcation points at c=±0.423 and the Hopf bifurcation points at c=±0.203. From the phase portrait on the (x,y) plane plotted in Figure 4(a), the attractors oscillate around two equilibrium points, as denoted by E±. In addition, a typical bursting oscillation is observed in the time history in Figure 4(b). The trajectory jumps between two spiking states SP±, which are symmetric to each other.Figure 4
Bursting oscillation forA=0.430. (a) Phase portrait on the (x,y) plane; (b) time history; (c) transformed phase portrait on the (c,x) plane with the bifurcation diagram.
(a) (b) (c)In order to reveal the mechanism of the bursting oscillation, we turn to the transformed portrait phase with the bifurcation diagram, as presented in Figure4(c). The trajectory starts from the point A with c=0.430, where c is taken at its positive extreme value. At this point, the attractor is categorized as the stable limit cycle LC1, thereby allowing the trajectory to move within the vector field of LC1, and repetitive spiking movement with large amplitude is observed. Parameter c then decreases, as indicated by the arrow. When the parameter c decreases to c=0.203, where the Hopf bifurcation point supH1 is observed, thereby resulting in the disappearance of stable limit cycle LC1 and the appearance of stable equilibrium point E+.Due to the slow passage effect [23], the trajectory will not move within the vector field of E+ immediately, it still moves in the vector field of LC1 for a short time, and then it turns into the vector field of E+. The repetitive spiking state movement is terminated, and the quiescent state movement is created. Parameter c then continues to decrease to c=-0.420, where the fold bifurcation point LP2 is located, thereby generating a trajectory jump into the vector field of LC2 and resulting in the termination of the quiescent state movement and the creation of repetitive spiking state movement. As the parameter c decreases to c=-0.430, it reaches its negative extreme value, after which the parameter c exhibits a gradual increase.When the parameterc increases to c=-0.203, where the other Hopf bifurcation point supH2 is observed, it results in the disappearance of the stable limit cycle LC2 and the appearance of the stable equilibrium point E-. The trajectory still moves in the vector filed of LC2 for a short time according to the slow passage effect, after which it turns into the vector field of E-. The Hopf bifurcation point supH2 terminates the repetitive spiking state movement and exhibits the quiescent state movement. When the parameter c increases to c=0.420, where the fold bifurcation point LP1 is observed, it results in the disappearance of the quiescent state movement and the appearance of the repetitive spiking state movement. In addition, this increase results in the loss of E-’s stability and the generation of the stable limit cycle LC1. Parameter c then continues to increase to c=0.430. The return of the trajectory back to the starting point A finished one periodic movement.The bursting oscillation can be classified as symmetric “fold/Hopf” bursting, since the repetitive spiking related to the limit cycle is generated by fold bifurcation and disappears by Hopf bifurcation.


## 3.3. Bursting Oscillation in Case C

In case C forb=2, twofold bifurcation points LP1,2(±0.539,±0.806) and two Hopf bifurcation points supH1,2(±1.222,±1.834) are observed on the (c,x) plane, resulting in the transition between the repetitive spiking and quiescent states. Two different types of bursting oscillations are observed when the amplitude A is taken at different values.


### 3.3.1. Symmetric “Fold/Fold” Bursting Oscillation

When the amplitudeA is taken at A=1.0, the slow variable c changes between −1.0 and 1.0, thereby implying the presence of only twofold bifurcation points at c=±0.539. Under this condition, three equilibria are present: two stable foci denoted by E± and one unstable saddle point denoted by E0, as presented in Figure 5.Figure 5
Distribution of the equilibria forA=1.0.From the phase portrait on the(x,y) plane, as presented in Figure 6(a), the trajectory oscillates around two equilibrium points, namely, E±. A typical bursting oscillation is observed in the time history diagram in Figure 6(b). In addition, the trajectory jumps between two spiking states, specifically SP±, which are symmetrical to each other.Figure 6
Bursting oscillation forA=1.0. (a) Phase portrait on the (x,y) plane; (b) time history; (c) transformed phase portrait on the (c,x) plane with equilibrium branches.
(a) (b) (c)In order to reveal the mechanism of the bursting oscillation, we turn to the transformed phase portrait on the(c,x) plane with equilibrium branches, as presented in Figure 6(c). The trajectory starts from the point A with c=1.0, where the parameter c is taken at its positive extreme value. The trajectory moves within the vector field of the stable focus E+ that is associated with the quiescent state. When the parameter c decreases to c=-0.539, where the fold bifurcation point LP2 is observed, thereby resulting in a trajectory jump to the vector field of the other stable focus E- and ending the quiescent state. The trajectory exhibits high-frequency oscillation with a gradual decrease in the amplitude, which corresponds to the spiking state. As the parameter c decreases to its negative value, specifically c=-1, it may generate a gradual increase that corresponds to the quiescent state.When the parameterc increases to c=0.539, where the fold bifurcation point LP1 is observed, it results in a trajectory jump to the vector field of E+ and terminating the quiescent state. In addition, the trajectory is characterized as high-frequency oscillation with a gradual decreasing in the amplitude, which corresponds to the spiking state. As the parameter c increases to c=1, the trajectory returned back to its starting point A, thereby finishing one periodic movement.Bursting oscillation can be created since the system periodically switches between the two stable fociE+ and E- by twofold bifurcations and the two parts are symmetrical. So the bursting pattern can be classified as symmetric “fold-fold” type.


### 3.3.2. Symmetric “SupHopf/SupHopf” Bursting Oscillation via “Fold/Fold” Hysteresis Loop

A continues to increase to the range of A>1.222, for example, A=2, thereby implying that the parameter c not only passes through the fold bifurcation points LP1 and LP2 at c=±0.539, but also passes through the supercritical Hopf bifurcation points supH1 and supH2 at c=±1.222, where the Hopf bifurcation occurs.According to the phase portrait on the(x,y) plane in Figure 7(a), the trajectory oscillates around two stable foci, namely, E±. In addition, a typical bursting oscillation is observed in the time history diagram in Figure 7(b), and the trajectory jumps between two spiking states, specifically SP±, both of which are symmetrical to each other. However, a great difference is observed between Figure 7(a) and Figure 6(a) because two limit cycles are observed in Figure 7(a).Figure 7
Bursting oscillation forA=2. (a) Phase portrait on the (x,y) plane; (b) time history; (c) overlap of the transformed phase portrait on the (c,x) plane with the system bifurcation diagram.
(a) (b) (c)In order to explore the mechanism of the bursting oscillation in Figure7, we turn to the overlap of the transformed phase portrait on the (c,x) plane with the system bifurcation diagram, as presented in Figure 7(c). We assume that the trajectory, starting from the point A with c=2, moves within the vector field of the limit cycle LC1, corresponding to the spiking state. When the parameter c decreases to c=1.222, where the Hopf bifurcation point supH1 is observed, thereby resulting in the disappearance of limit cycle LC1 and the presence of the stable focus E+. In addition, the Hopf bifurcation point supH1 resulted in the termination of the spiking state and the appearance of the quiescent state. As the parameter c decreases to c=-0.539, where the fold bifurcation point LP2 is observed, the trajectory jumps into the vector field of the stable focus E-, and the fold bifurcation point LP2 also results in the termination of the quiescent state and the appearance of the spiking state. When the parameter c decreases to c=-1.222, E- loses its stability via Hopf bifurcation, and a stable limit cycle LC2 is observed. As the parameter c decreases to an negative extreme value, specifically c=-2, it exhibits a gradual increase.When the parameterc increases to c=-1.222, LC2 disappears and E- is observed due to Hopf bifurcation. In addition, this increase results in the termination of the spiking state and the appearance of the quiescent state. As the parameter c increases to c=0.539, where the fold bifurcation point LP1 is observed, which results in the trajectory jumps to the vector field of E+, thereby terminating the quiescent state and generating the presence of the spiking state. When the parameter c exhibits an increase to c=1.222, where the Hopf bifurcation point supH1 is observed, E+ loses its stability and the stable limit cycle LC1 is produced via Hopf bifurcation. As the parameter c increases to c=2, the trajectory returns back to its starting point A, thereby finishing one periodic movement.From Figure7(c), a hysteresis loop produced by the equilibrium branches between twofold bifurcation points, specifically LP1,2, is observed. As a result, this bursting is classified as symmetric “Hopf/Hopf” type via the “fold/fold” hysteresis loop.


## 3.3.1. Symmetric “Fold/Fold” Bursting Oscillation

When the amplitudeA is taken at A=1.0, the slow variable c changes between −1.0 and 1.0, thereby implying the presence of only twofold bifurcation points at c=±0.539. Under this condition, three equilibria are present: two stable foci denoted by E± and one unstable saddle point denoted by E0, as presented in Figure 5.Figure 5
Distribution of the equilibria forA=1.0.From the phase portrait on the(x,y) plane, as presented in Figure 6(a), the trajectory oscillates around two equilibrium points, namely, E±. A typical bursting oscillation is observed in the time history diagram in Figure 6(b). In addition, the trajectory jumps between two spiking states, specifically SP±, which are symmetrical to each other.Figure 6
Bursting oscillation forA=1.0. (a) Phase portrait on the (x,y) plane; (b) time history; (c) transformed phase portrait on the (c,x) plane with equilibrium branches.
(a) (b) (c)In order to reveal the mechanism of the bursting oscillation, we turn to the transformed phase portrait on the(c,x) plane with equilibrium branches, as presented in Figure 6(c). The trajectory starts from the point A with c=1.0, where the parameter c is taken at its positive extreme value. The trajectory moves within the vector field of the stable focus E+ that is associated with the quiescent state. When the parameter c decreases to c=-0.539, where the fold bifurcation point LP2 is observed, thereby resulting in a trajectory jump to the vector field of the other stable focus E- and ending the quiescent state. The trajectory exhibits high-frequency oscillation with a gradual decrease in the amplitude, which corresponds to the spiking state. As the parameter c decreases to its negative value, specifically c=-1, it may generate a gradual increase that corresponds to the quiescent state.When the parameterc increases to c=0.539, where the fold bifurcation point LP1 is observed, it results in a trajectory jump to the vector field of E+ and terminating the quiescent state. In addition, the trajectory is characterized as high-frequency oscillation with a gradual decreasing in the amplitude, which corresponds to the spiking state. As the parameter c increases to c=1, the trajectory returned back to its starting point A, thereby finishing one periodic movement.Bursting oscillation can be created since the system periodically switches between the two stable fociE+ and E- by twofold bifurcations and the two parts are symmetrical. So the bursting pattern can be classified as symmetric “fold-fold” type.


## 3.3.2. Symmetric “SupHopf/SupHopf” Bursting Oscillation via “Fold/Fold” Hysteresis Loop

A continues to increase to the range of A>1.222, for example, A=2, thereby implying that the parameter c not only passes through the fold bifurcation points LP1 and LP2 at c=±0.539, but also passes through the supercritical Hopf bifurcation points supH1 and supH2 at c=±1.222, where the Hopf bifurcation occurs.According to the phase portrait on the(x,y) plane in Figure 7(a), the trajectory oscillates around two stable foci, namely, E±. In addition, a typical bursting oscillation is observed in the time history diagram in Figure 7(b), and the trajectory jumps between two spiking states, specifically SP±, both of which are symmetrical to each other. However, a great difference is observed between Figure 7(a) and Figure 6(a) because two limit cycles are observed in Figure 7(a).Figure 7
Bursting oscillation forA=2. (a) Phase portrait on the (x,y) plane; (b) time history; (c) overlap of the transformed phase portrait on the (c,x) plane with the system bifurcation diagram.
(a) (b) (c)In order to explore the mechanism of the bursting oscillation in Figure7, we turn to the overlap of the transformed phase portrait on the (c,x) plane with the system bifurcation diagram, as presented in Figure 7(c). We assume that the trajectory, starting from the point A with c=2, moves within the vector field of the limit cycle LC1, corresponding to the spiking state. When the parameter c decreases to c=1.222, where the Hopf bifurcation point supH1 is observed, thereby resulting in the disappearance of limit cycle LC1 and the presence of the stable focus E+. In addition, the Hopf bifurcation point supH1 resulted in the termination of the spiking state and the appearance of the quiescent state. As the parameter c decreases to c=-0.539, where the fold bifurcation point LP2 is observed, the trajectory jumps into the vector field of the stable focus E-, and the fold bifurcation point LP2 also results in the termination of the quiescent state and the appearance of the spiking state. When the parameter c decreases to c=-1.222, E- loses its stability via Hopf bifurcation, and a stable limit cycle LC2 is observed. As the parameter c decreases to an negative extreme value, specifically c=-2, it exhibits a gradual increase.When the parameterc increases to c=-1.222, LC2 disappears and E- is observed due to Hopf bifurcation. In addition, this increase results in the termination of the spiking state and the appearance of the quiescent state. As the parameter c increases to c=0.539, where the fold bifurcation point LP1 is observed, which results in the trajectory jumps to the vector field of E+, thereby terminating the quiescent state and generating the presence of the spiking state. When the parameter c exhibits an increase to c=1.222, where the Hopf bifurcation point supH1 is observed, E+ loses its stability and the stable limit cycle LC1 is produced via Hopf bifurcation. As the parameter c increases to c=2, the trajectory returns back to its starting point A, thereby finishing one periodic movement.From Figure7(c), a hysteresis loop produced by the equilibrium branches between twofold bifurcation points, specifically LP1,2, is observed. As a result, this bursting is classified as symmetric “Hopf/Hopf” type via the “fold/fold” hysteresis loop.


## 4. Effects of the Excitation Amplitude and Frequency on the Bursting Dynamics

We have investigated the dynamical mechanisms of bursting oscillations. In this section, we focus on the effects of excitation amplitude and excitation frequency on the bursting dynamics.First, the effect of the excitation amplitude on the bursting dynamics is considered. Let us take the symmetric “fold/fold” bursting oscillation in case C as an example. According to the phase portrait related toc in Figures 6 and 7, bursting oscillation is created when the slow excitation c passes through twofold bifurcation values cLP1 and cLP2 periodically. Since c=Asin⁡ωt, if the excitation A satisfies A<cLP, the bursting oscillations are not created because the slow excitation cannot pass through the twofold bifurcation values, while if A>cLP, they may be created.Second, the effect of the excitation frequency on the bursting dynamics is considered. The three bursting patterns are symmetrical in this article. Therefore, the time interval between two adjacent spikes of bursting in one oscillation period is half a period, which is the minimum time intervalΔT=1/2T=π/ω between two adjacent spikes, where T is the period of one movement. Take the waveforms in Figure 7 as an example, when ω1=0.0025, ω2=0.005, and ω3=0.01, respectively, the corresponding time intervals between two adjacent spikes are measured, as presented in Figure 8, of which the numerical results agreed well with the analytical ones.Figure 8
Numerical results of the time interval between two adjacent spikes of bursting oscillations corresponding to Figure7 for the slow excitation frequency for ω1=0.0025, ω2=0.005, and ω3=0.01, respectively. (a) ΔT1≈1256; (b) ΔT2≈628; and (c) ΔT3≈314.
(a) (b) (c)


## 5. Conclusions

The Shimizu-Morioka system may exhibit bursting oscillations with different waveforms under certain external forcing conditions. This article theoretically explores the mechanisms of different bursting patterns in the presence of two different frequency scales between the exciting frequency and the natural frequency. First, we regard the slow excitationAsin⁡ωt as a slow-varying parameter c. Second, according to the two bifurcation sets of (c,b), three cases are considered: namely, b=0.65, b=1.2, and b=2. Third, we investigate different bursting oscillations of different cases of b: for b=0.65, one type of bursting oscillation, namely, symmetric “homoclinic/homoclinic” bursting, is created when c slow-varying passes through two homoclinic bifurcation points; for b=1.2, symmetric “fold/Hopf” is created when c slow-varying passes through twofold bifurcation points and two Hopf bifurcation points; for b=2, two types of bursting oscillations, namely, symmetric “fold/fold” bursting when c slow-varying passes through only two fold bifurcation points, and symmetric “Hopf/Hopf” bursting via “fold/fold” hysteresis loop is created when c slow-varying passes through twofold bifurcation points and two Hopf bifurcation points. In a word, these four bursting oscillations are created due to the coexistence of different attractors. Fourth, we introduce two analysis strategies, specifically the transformed phase portrait and the overlap of the phase portrait related to c with the bifurcation diagram, to explore the bursting oscillation mechanism in terms of the bifurcations and equilibria structures. Last of all, we find that the bursting oscillation can be easily controlled by the amplitude and frequency of the slow excitation.



---

*Source: 5207910-2018-03-14.xml*