Homoclinic Bifurcations in Planar Piecewise-Linear Systems
The problem of homoclinic bifurcations in planar continuous piecewise-linear systems with two zones is studied. This is accomplished by investigating the existence of homoclinic orbits in the systems. The systems with homoclinic orbits can be divided into two cases: the visible saddle-focus (or saddle-center) case and the case of twofold nodes with opposite stability. Necessary and sufficient conditions for the existence of homoclinic orbits are provided for further study of homoclinic bifurcations. Two kinds of homoclinic bifurcations are discussed: one is generically related to nondegenerate homoclinic orbits; the other is the discontinuity induced homoclinic bifurcations related to the boundary. The results show that at least two parameters are needed to unfold all possible homoclinc bifurcations in the systems.
Nonsmooth dynamical systems are naturally used to model many physical processes, such as impacting, friction, switching, and sliding systems. The study of the nonsmooth dynamical systems has attracted more and more attention in the recent decades. Piecewise-smooth systems, as an important branch of the nonsmooth dynamical systems, involve collision systems, Filippov systems, higher-order discontinuity systems, and so forth [1, 2]. In particular, the study of piecewise-linear systems is significantly important because it not only describes some processes such as circuits [3, 4], but also enables us to locally understand the bifurcation phenomena in the nonlinear systems .
This paper studies the planar piecewise-linear continuous vector fields with two zones. Without loss of generality, the considered plane is divided into two half-planes by the boundary coinciding with the vertical axis. The system is linear in each of the half-planes and continuous along the vertical axis. In 1998, Freire et al.  studied discontinuous bifurcations, proved that there exists at most one limit cycle, and proved that if the limit cycle exists then it is either attracting or repelling, which solved the problem proposed in 1991 by Lum and Chua . di Bernardo et al.  studied discontinuity induced bifurcations, including boundary equilibrium bifurcations and grazing bifurcations of limit cycles. Concerning the limit cycle problem, Simpson and Meiss  investigated the Hopf branch in this kind of systems and mentioned homoclinic loops in . One can refer to [8, 9] for the periodic orbit problems in piecewise-linear systems with multiple nonsmooth boundaries. Bifurcation of limit cycles by perturbing a piecewise-linear Hamiltonian system is considered in [10, 11]. For the study of bifurcations in piecewise-linear systems with some application problems, see also .
In this paper, we further study the existence problem of homoclinic orbits and homoclinic bifurcations in the above-mentioned planar piecewise-linear systems. Starting with transforming the system into a canonical form, we found that its homoclinic orbits exist only in two cases: one is the saddle-focus (or saddle-center) system, which is called nondegenerate homoclinic orbits; the other one has two nodes coinciding on the vertical axis which have opposite stability, which is called degenerate homoclinic system. For both cases, the necessary and sufficient conditions for the existence of homoclinic orbits are established and used for the study of the homoclinic bifurcation problem. For the nondegenerate case, we find that there are two kinds of homoclinic bifurcations: one is generic, that is, the limit equilibria of the homoclinic orbits digress from the vertical axis; the other is nongeneric, that is, the discontinuity induced homoclinic bifurcation. However, for the degenerate case, only discontinuity induced homoclinic bifurcation occurs. Finally, we point out that at least two parameters are needed to unfold all possible homoclinic bifurcations in the mentioned systems.
The rest of the paper is outlined as follows. In Section 2, we give a brief overview of the concerned system and transform the system into a canonical form. In Section 3, the existence conditions of homoclinic orbits for both nondegenerate and degenerate cases are provided. Then we study the nondegenerate homoclinic bifurcations and the discontinuity induced homoclinic bifurcations in Section 4. Conclusions are drawn in Section 5.
2. Canonical Forms
We start with a planar piecewise-linear system: where , and are matrices. By the requirement of continuity on the -axis, we may assume that the system (1) is of the form where Here we call the left (right) half-system of (2) and the region of the plane to be the left (right) zone. Suppose that . An equilibrium of (2), which is of course isolated, is called visible if it is an equilibrium of the left (right) half-system and (); otherwise, it is invisible. Let Then it is easy to check that if and only if the equilibria of (2) both coincide on the -axis. We mention that if the equilibrium of one half-system in (2) lies on the -axis, then by continuity it is also the equilibrium of the other half-system.
A homoclinic orbit has both -limit and -limit set consisting of one equilibrium, called its limit equilibrium, that with the orbit constitutes the homoclinic loop. A homoclinic orbit is called degenerate if its limit equilibrium lies on the -axis; otherwise, it is nondegenerate. The following examples for homoclinic orbits describe the two cases.
Example 1. Nondegenerate homoclinic orbits. Take , and in (3).(a)Set , and . Then (2) is a saddle-focus system with the nondegenerate homoclinic orbit as shown in Figure 1(a).(b)Set , and . Then (2) is a saddle-center system with the nondegenerate homoclinic orbit as shown in Figure 1(b).
From the above examples, we found that if a homoclinic orbit exists then some of the orbits will cross the -axis from left to right and some from right to left. That means that the sign of for the vector field (2) restricted on the -axis will vary with . So we should assume that . We mention in passing that the existence of a limit cycle for (2) also implies for the same reason. See also .
In order to study homoclinic orbits, it is convenient to simplify system (2) to a canonical form.
Lemma 3. For system (2), let . Then there exists a linear transformation such that (2) becomes the following canonical system: where and are, respectively, the trace and determinant of the matrixes , and is given by (5), that is, .
Remark 6. Write We see that, on the -axis, , the first component of the vector field (7), becomes . So any orbit of (7) crossing is from left to right, and crossing is from right to left. Moreover, for , the orbit crossing the origin is tangent to the -axis at the origin.
Remark 7. The eigenvalues at the equilibria (12) are We may then classify each half-system of (7) by and in the trace-determinant plane as shown in Figure 3. There are 5 types of linear systems: saddle (S), attracting node (AN), repelling node (RN), attracting focus (AF), and repelling focus (RF).
To end this section, we introduce two propositions .
Proposition 9. System (7) has at most one limit cycle or homoclinic orbit. It is either attracting or repelling.
Proposition 10. If system (7) has a limit cycle or a homoclinic loop, then .
3. Existence of Homoclinic Orbits
We focus on conditions for system (7) to have homoclinic orbits in this section. The homoclinic orbits can be divided into two classes: nondegenerate homoclinic orbits for visible saddle-focus (or saddle-center) systems and degenerate homoclinic orbits for twofold node systems.
3.1. Nondegenerate Case
Proof. Without loss of generality, we assume that the limit equilibium of is in the left zone, denoted by . Due to piecewise-linearity of (7), the orbit starting from should cross the positive -axis to enter the right zone and then back to the original . So is a visible saddle in the left zone. By the generalized Poincaré-Bendixon theorem , inside the homoclinic loop of there is one equilibrium , which is a visible focus (or center) in the right zone. So (7) is a saddle-focus (or saddle-center) system.
Remark 12. From Remark 5, it easy to see that the existence of nondegenerate homoclinic orbits implies . Moreover, the saddle of (7) as the limit equilibrium of is in the left zone for and in the right zone for .
The following theorem gives conditions for the existence of nondegenerate homoclinic orbits.
Theorem 13. Let system (7) be a saddle-focus (or saddle-center) system with a visible saddle in one zone and a visible focus (or center) in the other zone. Denote by the eigenvalues of the saddle half-system and by , the eigenvalues of the focus (or center) half-system. Write Then there exists a nondegenerate homoclinic orbit if and only if
Proof. We may assume that the saddle of (7) is in the left zone and the focus (or center) is in the right zone, denoted by and , respectively. So one has
By using (15), the eigenvalues and eigenvectors of the saddle half-system are, respectively,
and the eigenvalues of the focus (or center) half-system are , where
If the nondegenerate homoclinic orbit exists, will start with the unstable manifold of the saddle to intersect at a point, say, , then along an orbit of the right half focus (or center) system in the right zone to intersect at another point, say, , and will finally return along the stable manifold of in the left zone, as shown in Figure 4.
In order to find conditions for the existence of such homoclinic orbit , we consider the right zone governed by the focus (or center) half-system in (7). Firstly, we transform the right half-system by applying the coordinate transformation where Note that (22) is so (22) is -axis invariant, and the right half-system of (7) becomes The equilibrium is
By using the polar coordinate transformation Equation (25) then becomes Rewrite (28) as which has the solution where is a constant.
We will define a half-return map from to by such that for there exists an orbit from to within the right zone. In fact, we let the points and in the polar coordinate (27) be given by and , respectively. From (30), it follows that
As shown in Figure 5, one has where Then (32) becomes
In general, for any , let and be still given by and , respectively, as shown in Figure 5. We mention that and are determined uniquely by and due to . Write where the parameter . If and are fixed, we often omit the parameter . From the above discussion, we see that if and only if (31) holds. That means that We point out that, by using (34), formula (36) becomes
Next we consider the left zone governed by the saddle half-system in (7). If the nondegenerate homoclinic orbit exists, then, by using (20), it is not hard to check that the intersection points and are given by So and . We mention that, from (24), one has .
From the discussion of (37) and (38), it follows that the existence of the homoclinic orbit is equivalent to . Substituting for (39) by and , we have By using (26), (34), and (40), a further simplification yields which is just formula (16). The result is then obtained.
For the case , the conclusion is similar; just exchange and .
Remark 15. From , in (36), we find in the second term of (42) that On the other hand, in the first term of (42), one has where . From (43) and (44), we find that the existence of a nondegenerate homoclinic loop implies , as stated in Proposition 10.
Remark 16. In particular, we consider the saddle-center system, where . From (42), one has
it is easy to see that and are equivalent to and , respectively.
On the other hand, it is not hard to check that if and , then . In this case, all the trajectories inside the homoclinic loop are closed around the center.
Remark 17. It is easy to see that is a decreasing function of . From (38), , we have So and have the same sign.
3.2. Degenerate Case
For , the origin is the only equilibrium of both half-systems. If the system has twofold node with opposite stability, then, as pointed in Freire et al. , there exists a continuum of homoclinic loops, consisting of degenerate homoclinic orbits. Conversely, the existence of degenerate homoclinic orbits implies that the equilibrium is neither saddle nor focus (or nor center). So one has the following result.
Remark 19. For a continuum of homoclinic loops, see Example 2. In fact, by calculating the eigenvalue (19) with the corresponding eigenvectors (20), one may find that all the homoclinic loops lie between an unstable manifold and a stable manifold.
4. Homoclinic Bifurcations
Based on the study of homoclinic orbits in the last section, now the homoclinic bifurcations in (7) are discussed. We first provide the generic homoclinic bifurcation related to nondegenerate homoclinic orbits and then study the problem of discontinuity induced homoclinic bifurcations. The latter will be divided into two cases to study: nondegenerate and the degenerate.
4.1. Generic Homoclinic Bifurcations
In this subsection, we discuss bifurcations related to the nondegenerate homoclinic orbits, called generic homoclinic bifurcations. For simplicity, we still assume that the left half-system in (7) is saddle, so . Two kinds of homoclinic bifurcations are studied: for the saddle-focus system and for the saddle-center system.
For the saddle-focus system, . As we saw in Theorem 13 and Remark 15, the existence of a nondegenerate homoclinic loop depends on and . Varying the parameters related to slightly, the homoclinic loop may disappear.(1)When , the right zone is governed by an RF system. From Remark 17 and the analysis of the phase portraits, we find that, for , one has and there is no periodic orbits in the system and that, for , one has and there is an attracting limit cycle around the focus. (2)When , in the right zone is an AF system. For , there is no periodic orbits in the system. For , there is a repelling limit cycle around the focus.
We mention here that the system has at most one limit cycle by Proposition 9.
Now we discuss the saddle-center system case; that is, . (1)For vary and , there exists a limit cycle in the right zone tangent to the -axis. It is the unique limit cycle. (2)For and fixed vary , in the right zone is a focus system, no limit cycle in this case. (3)Varying , and , it turns into the saddle-focus system case.
Example 21. A bifurcation of saddle-center systems is shown in Figure 7, where the parameters are , and , and take for (a), for (b), and for (c).
In order to summarize the above conclusions we let , , and be fixed and consider the -plane, here and . Figure 8 is the homoclinic bifurcation diagram for the nondegenerate case. The line divides the plane into two parts: upper is for and bellow is for . All the trajectories are clockwise. From Figure 8, we can find that the parameters and unfold all possible generic homoclinic bifurcation behaviors.
4.2. Discontinuity Induced Homoclinic Bifurcations
In this subsection, we study the problem of the discontinuity induced homoclinic bifurcations (DIHBs) related to variation of boundary equilibria. Herein we take as a bifurcation parameter and discuss two kinds of DIHBs: the nondegenerate case and the degenerate case.
4.2.1. Nondegenerate Case
We still assume that the equilibrium of (7) in the left zone is a visible saddle for , and then the right zone has a visible focus (or center). It is obvious that whenever both of the visible equilibria will coincide on the nonsmooth boundary -axis and become invisible as , which is a nonsmooth fold boundary equilibrium bifurcation; see also . From the discussion in Section 3, the homoclinic orbits exist for and shrink to an equilibrium gradually as the scale parameter decreases from positive to then disappear as . This is the discontinuity induced homoclinic bifurcation in the nondegenerate case.
Example 22. Following Example 1(a), set , and in (7) and (8). Take near zero. Figures 9(a), 9(b), and 9(c) show a discontinuity induced homoclinic bifurcation of the system for , , and , respectively.
In Figure 9, we find that, for , the system is saddle-focus with a nondegenerate homoclinic orbit as shown in (a), that, for , two equilibria coincide at the origin and the homoclinic orbit becomes a broken line as in (b), and that for both equilibria become invisible and the homoclinic orbits disappear as in (c).
4.2.2. Degenerate Case
As mentioned in Section 3, under the conditions of Theorem 18 there exists a continuum of homoclinic loops when . Varying the parameter , the system (7) has a type of boundary equilibrium bifurcation, called persistent-node bifurcation. Here as the bifurcation parameter crosses the critical value , the visible node of one half-system in (7) becomes invisible, and the other invisible node becomes visible, while all homoclinic orbits disappear for . This is a discontinuity induced homoclinic bifurcation in the degenerate case.
Example 23. Following Example 2, for (7) and (8), set , and . We take near zero. Figures 10(a), 10(b), and 10(c) show the discontinuity induced homoclinic bifurcation for degenerate case for , , and , respectively.
We have studied in detail the homoclinic bifurcations in the planar piecewise-linear systems with two zones. Firstly, the system is transformed to a canonical form. Then both nondegenerate and degenerate homoclinic orbits are investigated. The necessary and sufficient conditions for the existence of both a nondegenerate homoclinic orbit and a continuum of degenerate homoclinic orbits are established for discussion of the homoclinic bifurcation problem. A generic homoclinic bifurcation diagram is shown as in Figure 8. We have found that there are two kinds of DIHBs, one is nondegenerate and the other is degenerate. Nonsmooth fold boundary equilibrium bifurcation occurs accompanying DIHB in nondegenerate case, and persistent-node boundary equilibrium bifurcation occurs accompanying DIHB in degenerate case. To sum up, at least two parameters are needed to unfold all possible homoclinic bifurcations.
This work was supported by the NNSF of China (Grant nos. 11072274 and 11272169).
M. di Bernardo, C. J. Budd, A. R. Champneys, and P. Kowalczyk, Piecewise-Smooth Dynamical Systems, vol. 163 of Applied Mathematical Sciences, Springer, London, UK, 2008.View at: MathSciNet
M. di Bernardo, D. J. Pagano, and E. Ponce, “Nonhyperbolic boundary equilibrium bifurcations in planar Filippov systems: a case study approach,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 18, no. 5, pp. 1377–1392, 2008.View at: Publisher Site | Google Scholar | MathSciNet
R. Lum and L. O. Chua, “Global properties of continuous piecewise linear vector fields—part I: simplest case in ℝ2,” International Journal of Circuit Theory and Applications, vol. 19, no. 3, pp. 251–307, 1991.View at: Google Scholar
D. Simpson, Bifurcations in Piece-Wise Continuous Systems, World Scientific, 2010.
R. I. Leine and H. Nijmeijer, Dynamics and Bifurcations of Non-Smooth Mechanical Systems, vol. 18 of Lecture Notes in Applied and Computational Mechanics, Springer, Berlin, Germany, 2004.View at: MathSciNet
L. Perko, Differential Equations and Dynamical Systems, vol. 7 of Texts in Applied Mathematics, Springer, New York, NY, USA, 3rd edition, 2001.View at: MathSciNet