On the periodic solutions of the Michelson continuous and discontinuous piecewise linear differential system

Applying new results from the averaging theory for continuous and discontinuous differential systems, we study the periodic solutions of two distinct versions of the Michelson differential system: a Michelson continuous piecewise linear differential system and a Michelson discontinuous piecewise linear differential system. The tools here used can be applied to general nonsmooth differential systems.

the periodic solutions of continuous and discontinuous piecewise linear differential systems using new results on the averaging theory developed for these nonsmooth dynamical systems.
The applications of these new results on averaging theory for studying the periodic solutions of nonsmooth dynamical systems required convenient changes of variable for writing the initial differential system in the normal form of the averaging theory. We shall do these applications to a version of a continuous and a discontinuous piecewise linear differential systems coming from the Michelson differential system. Such piecewise linear differential systems also have been studied by other authors, see Carmona et al. (2012Carmona et al. ( , 2014Carmona et al. ( , 2010Carmona et al. ( , 2015Carmona et al. ( , 2008. We must mention that the ideas and tools used for doing these applications of the averaging theory for studying the periodic solutions of these two Michelson nonsmooth dynamical systems can be applied to general nonsmooth dynamical systems. The Michelson differential system is given bẏ with (x, y, z) ∈ R 3 and the parameter c ≥ 0. The dot denotes derivative with respect to an independent variable t, usually called the time. This system is due to Michelson (1986) for studying the traveling solutions of the Kuramoto-Sivashinsky equation. It also arises in the analysis of the unfolding of the nilpotent singularity of codimension three, see Dumortier et al. (2001) and Freire et al. (2002). This system has been largely investigated from the dynamical point of view. Michelson (1986) proved that if c > 0 is sufficiently large, then system (1) has a unique bounded solution which is a transversal heteroclinic orbit connecting the two finite singularities (− √ 2c, 0, 0) and ( √ 2c, 0, 0). When c decreases, there will appear a cocoon bifurcation (see Kokubu et al. 2007;-T Lau 1992;Michelson 1986). In Llibre and Zhang (2011), there is an analytical proof of the existence of a zero-Hopf bifurcation for system (1).
In Carmona et al. (2008Carmona et al. ( , 2010, the authors consider a continuous piecewise linear version of Michelson differential system changing the non linear function x 2 in (1) by the piecewise linear function |x|. For such system, they proved that some dynamical aspects of the Michelson system remain as the existence of a reversible T-point heteroclinic cycle.
Doing the change of variable (x, y, z, c) → (2ε X, 2εY, 2ε Z , 2εd) with d ≥ 0 and ε > 0 sufficiently small to the Michelson differential system (1), followed by the change of the function X 2 → |X |, and denoting again X , Y , and Z by x, y, and z, we obtain the systeṁ that we call the Michelson continuous piecewise linear differential system. We note that this system is reversible, because it is invariant under the change of variables (x, y, z, t) → (−x, y, −z, −t).
Many problems in physics, economics, biology, and applied areas are modeled by discontinuous differential systems, but there exist only few analytical techniques for studying their periodic solutions. In Llibre et al. (2015), the authors extended the averaging theory to the discontinuous differential systems. An improvement of this result for a much bigger class of discontinuous differential systems is given in Llibre and Novaes (2016).
If in the continuous Michelson differential system (1), we change the continuous function |x| by the discontinuous one |x| + sign(x), where we obtain the Michelson discontinuous piecewise linear differential system given bẏ In this paper, first, we study analytically the periodic solutions of the Michelson continuous and discontinuous piecewise linear differential systems. Thus, our first main result is the following.
Theorem 1 For all d > 0 and ε = ε(d) > 0 sufficiently small, the Michelson continuous piecewise linear differential system (2) has a periodic solution of the form:

Moreover, this periodic solution is linearly stable.
We recall that a periodic solution is asymptotically stable if all the eigenvalues corresponding to the fixed point of the Poincaré map associated with this solution have negative real part, then this periodic solution is locally asymptotically stable. If one of the eigenvalues has positive real part, the periodic solution is unstable. If all the eigenvalues have zero real parts, then we say that the periodic solution is linearly stable; in this case, the linear stability does not provide any information on the kind of stability that the periodic solution has when we take into account the nonlinear terms, for more details, see Theorem 11.6 of Verhulst (2000).
Note that the periodic orbit obtained in our Theorem 1 is not reversible, because their dominant terms (−d 2 π , d 2 π sin t , d 2 π cos t) are not invariant under the change of variables (x, y, z, t) → (−x, y, −z, −t). Therefore, this periodic orbit has no relation with the periodic orbit of the noose bifurcation studied by the Michelson continuous piecewise linear differential system (1.2) of Carmona et al. (2015) which is reversible. We also note that our Michelson continuous piecewise linear differential system (2) and system (1.2) of Carmona et al. (2015) do not coincide.
Theorem 1 is proved in Sect. 3. Its proof uses an extension of the classical averaging theory for smooth differential systems to continuous differential systems given in Buica et al. (2007) or Llibre et al. (2014).
Theorem 2 For ε > 0 sufficiently small, the Michelson discontinuous piecewise linear differential system (3) satisfies the following statements.
(a) If (−1 + 2d 2 )π < 0, then system (3) has two periodic solutions (x(t, ε), r (t, ε), θ(t, ε)) of the form: where and a takes the value of the two unique zeros of the function (3) has a periodic solution of the form (4) given by the unique zero of the function g(a) in the interval (−1, 1).

Preliminary
For proving Theorems 1 and 2, we apply two recent results from the averaging theory: one for the continuous piecewise linear differential systems and the other for the discontinuous piecewise linear differential systems. In this section, we present these results and some remarks necessary for their applications.

Continuous piecewise linear differential systems
From Theorem B of Llibre et al. (2014), we get the following result adapted to the next system (5), see also Theorems 11.5 and 11.6 of Verhulst (2000).
Theorem 3 Consider the following system: where D ⊂ R n is an open subset, ε is a small parameter, the functions F i : R × D → R n for i = 0, 1 and R : R × D × (−ε 0 , ε 0 ) → R n are T periodic in the variable t, and for each t ∈ R, the functions F 0 (t, .) ∈ C 1 , F 1 (t, .) ∈ C 0 , D x F 0 , and R ∈ C 0 are locally Lipschitz in the second variable. We denote by x(t, z, ε) the solution of system (5), such that x(0, z, ε) = z. Assume that there exists an open and bounded subset V of D with its closure V ⊂ D, such that for each z ∈ V , the solution x(t, z, 0) is T periodic. We denote by M z (t) the fundamental matrix solution of the variational equation: and det(D z F (a)) = 0, then for ε > 0 sufficiently small, system (5) has a T -periodic solution x(t, a ε , ε), such that a ε → a as ε → 0. Moreover, the linear stability type of the periodic solution x(t, a ε , ε) is given by the eigenvalues of the matrix D z F (a).
Note that the stability of the periodic solutions of system (5), when it is applied to the Michelson continuous piecewise linear differential system, can be obtained from the stability of a differential system associated with it. In fact given the continuous system (1), consider a band of amplitude ε > 0 around the plane x = 0 and a differentiable extension of the continuous system (1) to this band. Studying the limit of this extended differentiable system when ε → 0, we conclude that the linear stability of system (1) is given by the eigenvalues of D z F (a).

Discontinuous piecewise linear differential systems
Let D ⊂ R n be an open subset and h : R × D → R a C 1 function having 0 as regular value. Consider F 1 , F 2 : R × D → R n continuous functions and = h −1 (0). We define the Filippov's system aṡ where The manifold is divided in the closure of two disjoint regions, namely, Crossing region ( c ) and Sliding region ( s ): Consider the differential system associated with system (6): where χ + , χ − are the characteristic functions defined as Systems (6) and (7) do not coincide in , but applying the Filippov's convention for the solutions of systems (6) and (7) (see Filippov 1988) passing through a point (t, x) ∈ , we have that these solutions do not depend on the value of F(t, x), so the solutions are the same. Let P be the space formed by the periodic solutions of (7). If dim P = dim D = d, then the following result follows directly from Theorem B of Llibre and Novaes (2016).
Theorem 4 Consider the differential system: with F 1 i ∈ C 1 , for i = 0, 1, and R 1 , R 2 are continuous functions which are Lipschitz in the second variable, and all these functions are T -periodic functions in the variable t ∈ R.
For z ∈ D and ε > 0 sufficiently small, denote by x(t, z, ε) the solution of system (8), such that x(0, z, ε) = z. 1 (s, x(s, z, 0))ds where x(s, z, 0) is a periodic solution of (8) with ε = 0, such that x(0, z, 0) = z, and M(s, z) is the fundamental matrix of the variational systemẏ = D x F 0 (t, x(t, z, 0))y associated with the unperturbed system evaluated on the periodic solution x(s, z, 0), such that M(0, z) = I d. Moreover, we assume the following hypotheses.

Define the averaged function
(H − ) There exists an open bounded subset C ⊂ D, such that, for ε sufficiently small, every orbit starting in C reaches the set of discontinuity only at its crossing region. (H + ) For a ∈ C with F (a) = 0, there exists a neighborhood U ⊂ C of a, such that F (z) = 0, for all z ∈Ū \{a} and det(D z F (a)) = 0.
Then, for ε > 0 sufficiently small, there exists a T -periodic solution x(t, a ε , ε) of (8), such that a ε → a as ε → 0. Moreover, the linear stability of the periodic solution x(t, ε) is given by the eigenvalues of the matrix D z F (a).
The same arguments for computing the kind of stability of the obtained periodic solution from the eigenvalues of the Jacobian matrix for the continuous piecewise differential systems also work for the discontinuous piecewise differential systems.

Proof of Theorem 1
Doing to the Michelson continuous piecewise linear differential system (2) the change to cylindrical coordinates x = x, y = r sin θ and z = r cos θ , the system becomeṡ Taking θ as the new independent variable, we can write the previous differential system as The unperturbed system is For each (x 0 , r 0 ), the solution (x(θ, which is 2π periodic for all (x 0 , r 0 ) = (0, 0). When r 0 = 0, we have a straight line of equilibrium points. Now, note that the function F 0 (θ, (x, r )) = (r sin θ, 0) is C ∞ and in particular C 1 and that the function F 1 (θ, (x, r )) = (2d 2 − |x|)(sin 2 θ, cos θ) is C 0 , and both are Lipschitz. Therefore, the differential system (9) satisfies the assumptions of Theorem 3. Then, by Theorem 3, we need to calculate the averaged function: is the fundamental matrix of the variational differential system associated with system (11) evaluated on the periodic solution (x 0 + r 0 (1 − cos θ), r 0 ), such that M(0) is the identity matrix. Therefore, we have Note that g(θ ) = g(−θ) and g(θ ) are 2π periodic. Therefore To calculate this integral, we need to study the zeros of the function G(θ ) = x 0 +r 0 (1−cos θ).
As G(θ ) = 0 if and only if θ = ± arccos( x 0 +r 0 r 0 ) and the function arccos(x) takes real values when x ∈ [−1, 1], we have to consider the following three cases.
Case 3 Analogously to Case 1, we have Finally, coming back to system (2), we find the periodic solution: This concludes the proof of Theorem 1.

Proof of Theorem 2
Doing the change to cylindrical coordinates x = x, y = r sin θ , and z = r cos θ , the Michelson discontinuous piecewise linear differential system becomeṡ Now, taking as new independent variable, the angle θ we get the system: where the prime denotes the derivative with respect to θ . This differential system satisfies the assumptions of Theorem 4, so we shall apply it for finding some of its periodic solutions.
As in the study of the continuous differential system, we separate the calculation of the averaged function corresponding to the discontinuous system (13) in the same three cases that appear in the proof of Theorem 1.
Therefore, by continuity of the function g(X ) in the intervals (−∞, 0) and (0, ∞), it follows that g(X ) has one zero in the interval (−∞, 0) if (−1 + 2d 2 )π < 0 and that g(X ) always has one zero in the interval (0, ∞). Moreover, since the derivative g (X ) > 0 in those two intervals, such zeros are the unique zeros of the function g(X ).
This concludes the proof of Theorem 2.