Abstract

A class of polynomial differential systems with high-order nilpotent critical points are investigated in this paper. Those systems could be changed into systems with an element critical point. The center conditions and bifurcation of limit cycles could be obtained by classical methods. Finally, an example was given; with the help of computer algebra system MATHEMATICA, the first 5 Lyapunov constants are deduced. As a result, sufficient and necessary conditions in order to have a center are obtained. The fact that there exist 5 small amplitude limit cycles created from the high-order nilpotent critical point is also proved.

1. Introduction

In the qualitative theory of ordinary differential equations, bifurcation of limit cycles for planar polynomial differential systems which belong to second part of the Hilbert 16th problem is hot but intractable issue. Of course, this problem is far from being solved now; there are many famous works about this problem. Let be the maximum possible number of limit cycles bifurcating from infinity for analytic vector fields of degree . It was found that [1, 2], [3], [4].

When the critical point is a nilpotent critical point, let be the maximum possible number of limit cycles bifurcating from nilpotent critical points for analytic vector fields of degree . It was found that , , in [5], , in [6], and for a family of Kukles systems with 6 parameters in [7]. Recently, li and liu other found that , , [8, 9] employing the inverse integral factor method.

But when the order of nilpotent critical point is higher than three, it is very difficult to study the limit cycle because the inverse integral factor method could not be used to compute the singular values. There are few results about the bifurcations of limit cycles at a nilpotent critical point with high order.

In this paper, we will study the bifurcation of limit cycles and conditions of centers for a class of special system

Obviously, when , the system is

The origin is an element critical point, but it is a high order nilpotent critical point when .

This paper will be organized as follows. In Section 2, we state some preliminary knowledge given in [10] which is useful throughout the paper. In Section 3, we make some transformations to change system (2) into a system with an element singular. As an example, a special system is investigated. Using the linear recursive formulae in [10] to do direct computation, we obtain the first 5 Lyapunov constants and the sufficient and necessary conditions of center and the result that there exist 5 limit cycles in the neighborhood of the high-order nilpotent critical point are proved.

2. Preliminary Knowledge

The ideas of this section come from [10], where the center focus problem of critical points in the planar dynamical systems are studied. We first recall the related notions and results. For more details, please refer to [10].

In [10], the authors defined complex center and complex isochronous center for the following complex system: and gave two recursive algorithms to determine necessary conditions for a center and an isochronous center. We now restate the definitions and algorithms.

Lemma 1. For system (3), we can derive uniquely the following formal series: where , , such that

Definition 2. Let , , , . is called the singular point quantity of the origin of system (3) and is called the period constant of the origin of system (3).

Theorem 3. For system (3), the origin is a complex center if and only if , . The origin is a complex isochronous center if and only if , .

Theorem 4. For system (3), we can derive successively the terms of the following formal series: such that where , for all , , and for any integer , is determined by the following recursive formulae:

Theorem 5 (the constructive theorem of singular point quantities in [10]). A -order singular point quantity of system (3) at the origin can be represented as a linear combination of -order monomial Lie invariants and their antisymmetry forms, that is, where is a positive integer and is a rational number, and are -order monomial Lie invariants of system (3).

Theorem 6 (the extended symmetric principle in [10]). Let denote an elementary Lie invariant of system (3). If for all the symmetric condition is satisfied, then the origin of system (3) is a complex center. Namely, all singular point quantities of the origin are zero.

3. Simplification of System

In fact, system (1) is symmetric with axis when is even, so the origin is a center. Now we will consider system (1) when is odd.

By transformations

system (1) is changed into then by transformations system (11) could be changed into

Thus, the origin of system (13) is an element critical point. It could be investigated using the classical integral factor method.

Now, we consider the following system:

When , by those transformations, system (14) is changed into where

By transformation

system (15) is changed into where

After careful computation by using formula in (4), we have the following.

Theorem 7. For system (18), the first 5 Lyapunov constants at the origin are given by
When
When ,
When ,
When
In the above expression of , one has already let .

From Theorem 7, we obtain the following assertion.

Proposition 8. The first 5 Lyapunov constants at the origin of system (18) are zero if and only if one of the following conditions is satisfied:

Furthermore, we have the following.

Theorem 9. The origin of system (18) is a center if and only if the first 5 Lyapunov constants are zero; that is, one of the conditions in Proposition 8 is satisfied.

Proof. When one of conditions (25), (27), (28), (29), (30), (32), (33), and (34) holds, according to Theorems 6, we get all , .
When condition (26) holds, system (18) could be written as System (36) has an analytic first integral
When condition (31) holds, system (18) could be written as
System (38) has an analytic first integral
When condition (35) holds, system (18) could be written as
System (40) has an analytic first integral

Next, we will prove that when the critical point is a 5-order weak focus, the perturbed system of (15) can generate 5 limit cycles enclosing the origin of perturbation system (15).

Using the fact

we obtain the following.

Theorem 10. The origin of system (18) is a 5-order weak focus if and only if one of the following conditions is satisfied:

We next study the perturbed system of (15) as follows:

Theorem 11. If the origin of system (15) is a 5-order weak focus, for , making a small perturbation to the coefficients of system (15), then, for system (45), in a small neighborhood of the origin, there exist exactly 5 small amplitude limit cycles enclosing the origin .

Proof. It is easy to check that when condition (43) or (44) holds,
From the statement mentioned above, according to the classical theory of Bautin, there exist 5 limit cycles in a small enough neighborhood of the origin.

Acknowledgments

This research is partially supported by the National Nature Science Foundation of China (11201211, 61273012) and Project of Shandong Province Higher Educational Science and Technology Program (J12LI58).