Limit Cycles and Invariant Curves in a Class of Switching Systems with Degree Four

<jats:p>In this paper, a class of switching systems which have an invariant conic <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1"><mml:msup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:mi>c</mml:mi><mml:msup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>c</mml:mi><mml:mo>∈</mml:mo><mml:mi>R</mml:mi></mml:math>, is investigated. Half attracting invariant conic <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2"><mml:msup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:mi>c</mml:mi><mml:msup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>c</mml:mi><mml:mo>∈</mml:mo><mml:mi>R</mml:mi></mml:math>, is found in switching systems. The coexistence of small-amplitude limit cycles, large amplitude limit cycles, and invariant algebraic curves under perturbations of the coefficients of the systems is proved.</jats:p>


Introduction
It is well known that the 16th problem stated in 1900 by D. Hilbert is considered to be the most difficult problem in the 23 problems; it is far from being solved.Over past three decades, there have been many good results about this problem.As far as the maximal number of small-amplitude limit cycles which are bifurcated from an elementary center or focus is concerned, the best known result obtained by Bautin in 1952 [1] is (2) = 3, where () denotes the maximal number of small-amplitude limit cycles around a singular point with  being the degree of polynomials in the system.For cubicdegree system, many good results have also been obtained.For example, a cubic system was constructed by Lloyd and Pearson [2] to show 9 limit cycles with the aid of purely symbolic computation.Moreover, Yu and Tian [3] proved that there can exist 12 limit cycles around an elementary center in a planar cubic-degree polynomial system.As far as we know this is the best result obtained so far for cubic-degree polynomial systems with all limit cycles around a single singular point.For  ≥ 4, because of the difficulty of computation of focal values, there are very few results.An example of a quartic system with 8 limit cycles bifurcating from a fine focus [4] was given by Huang et al.Theory of rotated equations and applications to a population model can be found in [5]; they gave a new method to solve the center problem.
Center and the coexistence of large and small-amplitude limit cycles problems are two closely related questions of the 16th problem.Algebraic trajectories play an important role in the dynamical behavior of polynomial systems, so it has been an interesting problem in planar polynomial systems.Over the past twenty years, many interesting results were got for quadratic systems; the authors in [17,18] proved that quadratic systems with a pair of straight lines or an invariant hyperbola, ellipse, can have no limit cycles other than the possible ellipse itself.Furthermore, if there is an invariant line, there can be no more than one limit cycle.The case of parabola was considered in [19].For cubic systems, there exist different classes of cubic systems in which there may coexist an invariant hyperbola or straight lines with limit cycles (see [20][21][22][23][24][25][26][27][28]).For a given family of real planar polynomial systems of ordinary differential equations depending on parameters, the problem of how to find the systems in the family which become time-reversible was solved in [29].
In modelling many practical problems in science and engineering, switching systems have been widely used recently.The richness of dynamical behavior found in switching systems covers almost all the phenomena discussed in general continuous systems.For example, the maximum number of limit cycles bifurcating from the periodic orbits of the quadratic isochronous centers of switching system was studied in [30].In [31], limit cycles in a class of continuous and switching cubic polynomial systems were investigated.Bifurcation of limit cycles in switching quadratic systems with two zones was considered in [30] .In [32,33], the authors considered nonsmooth Hopf bifurcation in switching systems.Limit cycles bifurcating from centers of discontinuous quadratic systems were studied by Chen and Du [34].Switching Bautin system was also investigated in [35]; they got 10 limit cycles for this system. 2 -equivariant cubic systems were also investigated, and 14 limit cycles were obtained in [36].Bifurcation of periodic orbits by perturbing high-dimensional piecewise smooth integrable systems was investigated in [37].Bifurcation theory for finitely smooth planar autonomous differential systems was considered in [38].All results obtained show that the dynamical behavior of switching systems is more complex than continuous system.About algebraic invariant curves, as far as we know, there are few papers to consider switching system with algebraic invariant curves.In this paper we are concerned with the limit cycle problem and the center problem for a class of degree four polynomial differential systems which have an invariant conic  2 +  2 = 1,  ∈ , and we prove the coexistence of large elliptic limit cycle that contains at least four small-amplitude limit cycles generated by Hopf bifurcations.
The rest of the paper is organized as follows.In the next section, we prove that the switching system (1) has an invariant conic  2 +  2 = 1,  ∈ , and there exists a large limit cycle in switching system (1); half attracting invariant conic  2 +  2 = 1,  ∈ , is found in switching systems.In Section 3, the first eight Lyapunov constants will be computed; bifurcation of limit cycles and center conditions of (1) are investigated.Section 4 is devoted to discuss the number of limit cycles with different parameter  of (1).At last, coexistence of invariant curve and limit cycles of ( 1) is drawn in Section 5.

Invariant Curve and Large Limit Cycle of (1)
In this section, we will prove that the switching system (1) has an invariant conic  2 +  2 = 1,  ∈ , and there exists a large limit cycle in switching system (1).

Bifurcation of Limit Cycle and Center Conditions of (1)
First of all, it is easy to know that the origin of upper half system and lower half system is a fine focus if  1 ̸ = 0, so we let  =  1 = 0 in order to consider the center conditions and the number of small limit cycles.With the aid of symbolic computation, we obtain the following result.Theorem 3.For system ( ), the first eight Lyapunov constants at the origin are given by with two cases: (I)  = ± √ 2/2. 4 =  (64 + 21)  2  288 (4 + ) , where + 39259691311104000 24 ) . ( Note that in computing the above expressions   ,  = 2, . . ., 8, and The following proposition follows directly from Theorem 3.

Proposition 4. e first eight Lyapunov constants at the origin of system ( ) become zero if and only if one of the following conditions is satisfied:
ey are also the center conditions of system ( ).
Proof.When the conditions in (11) hold, system (1) can be brought to Obviously, the system is symmetric with the -axis, and so the origin is a center of system (14).When the conditions in (12) hold, system (1) can be rewritten as which is symmetric with the -axis, and so the origin is a center of system (15).When the conditions in (13) hold, system (1) becomes a continuous system By elementary integration, the above system in Ω = (, ) |  2 +  2 < 1 is topologically equivalent to the system which has the first analytic integral Remark .The phase plane of system ( 16) can be drawn by Maple; see Figure 2.
As far as limit cycles are concerned, it follows from Theorem 3 that at most 8 limit cycles can bifurcate from the origin of system (1).We have the following theorem.Theorem 6.If the origin of system ( ) is a 8th-order weak focus, then for 0 <  1 ,  2 ≪ 1, 8 small-amplitude limit cycles can bifurcate from the origin of the perturbed system ( ).
Proof.When the origin of system ( 1) is a 8th-order weak focus, the conditions should be satisfied.Furthermore, one has the following: When  = −1.9427,When  = −0.763201, ( So it implies that 8 small-amplitude limit cycles can bifurcate from the origin of the perturbed system (1).

Number of Limit Cycles with Different
Parameter  of (1) In this section, we devote to discuss the number of limit cycles with different parameter  of (1).The following theorem could be concluded from Theorem 3.
Theorem 7. e number of limit cycles with different parameter  of ( ) can be shown in the Table .Proof.Let  1 =  2 =  3 = 0, it is easy to obtain that If  = 0, it is easy to check that the origin is a three-order weak focus.Furthermore, if  ̸ = 0, when  = √ 2/2, the Lyapunov constants in case 1 yield that the origin is a seventh order weak focus.

Coexistence of Invariant Curve and Limit
Cycles of (1) From above discussion, we study the coexistence of invariant curve and limit cycles of (1), by perturbation method of small parameters, the following conclusions could be got easily; for example, when  = −2.09067,there exist eight small limit cycles at least and  2 − 2.09067 2 = 1 is an invariant algebraic curve.The distribution of limit cycle can be drawn in Figure 3.When  = − √ 2/2, there exist seven small limit cycles at least and  2 − ( √ 2/2) 2 = 1 is an invariant algebraic curve.The distribution of limit cycle can be drawn in Figure 4.
When  = 0, there exist three small limit cycles at least and  = 1 and  = −1 are two invariant lines.The distribution of limit cycle can be drawn in Figure 5.
When  = 0.581824, there exist eight small limit cycles and a large limit cycle  2 + 0.581824 2 = 1 at the same time, namely, nine limit cycles in total for this system.The distribution of limit cycle can be drawn in Figure 6.
When  = √ 2/2, there exist seven small limit cycles at least and there is a large limit cycle  2 + ( √ 2/2) 2 = 1 at the same time.The distribution of limit cycles can be drawn in Figure 7.
When  = 1, there exist six small limit cycles at least and there is a large limit cycle  2 +  2 = 1 at the same time.The distribution of limit cycles can be drawn in Figure 8.

Conclusion
In this paper, a class of switching systems is investigated; the coexistence of small limit cycles and algebraic an invariant curve is proves.An interesting phenomenon that the algebraic invariant curve  2 + 2 = 1,  > 0, can be half attracting is found.

Table 1 :
The maximum number of small limit cycles around the origin for different parameter c.