Abstract

With the aid of computer algebra system Mathematica 8.0 and by the integral factor method, for a family of generalized nilpotent systems, we first compute the first several quasi-Lyapunov constants, by vanishing them and rigorous proof, and then we get sufficient and necessary conditions under which the systems admit analytic centers at the origin. In addition, we present that seven amplitude limit cycles can be created from the origin. As an example, we give a concrete system with seven limit cycles via parameter perturbations to illustrate our conclusion. An interesting phenomenon is that the exponent parameter controls the singular point type of the studied system. The main results generalize and improve the previously known results in Pan.

1. Introduction

A famous problem for the plane analytic systems of differential equations is under what conditions the local phase portrait at a critical point is topologically equivalent to the local phase portrait of the linear part of the system at . This problem has been solved by Poincaré and Bendixson for hyperbolic critical points and for elementary critical points, that is, for points having zero determinant and nonzero trace linear part. Another problem is to characterize the local phase portrait at an isolated critical point .

When the matrix of the linear part at the origin is not identically null but has its eigenvalues which are equal to zero, at this moment, the origin is a nilpotent critical point. An analytic system having an isolated nilpotent singularity at the origin, in some suitable coordinates, can be written as follows:

Suppose that the function satisfies , . Amelkin et al. proved (see, for instance, [1]) that the origin of system (1) is a monodromic critical point (i.e., a center or a focus) if and only if where is a positive integer. Andreev [2] shows what the behavior of the solutions in a neighborhood of the nilpotent critical point is, except if it is a center or a focus (nilpotent center problem). This last result can not distinguish between a focus and a center. Takens [3] and Bogdanov [4] find easy formal normal forms for the nilpotent critical points. Moussu [5] has found the normal form for a nilpotent center and asks if there exists an analytic normal form for the nilpotent centers. Takens [3] proves that for any analytic nilpotent center there exists an analytic change of variables such that the new system can be written as a system of the form (1) and it is a time-reversible nilpotent center. However, the aforementioned result is difficult to implement because, in general, in order to decide if an analytic system has a nilpotent center, we must know explicitly the analytic change of variables which writes system (1) in the Berthier-Moussu normal form. Writing systems in a convenient normal form, Teixeira and Yang [6] study the relationship between time-reversibility and the center-focus problem for elementary singular points and nilpotent singular points. Giné [7] develops a method which provides necessary conditions for obtaining a local analytic integral in a neighborhood of a generalized nilpotent singular point. García and Giné [8] give a necessary condition to have local analytic integrability in an analytic nilpotent center. Giacomini et al. [9] studied the centers of planar analytic vector fields which are a limit of linear-type centers.

In this paper, we consider the system of differential equations in the plane whose origin is a generalized nilpotent singular point The origin of system (3) is a third-order nilpotent singular point when , while it is a total degenerate singular point when .

In Section 2, we give some preliminary knowledge concerning the nilpotent critical point. In Section 3, we transform the origin into a third-order nilpotent singular point by a homeomorphism. Then, we compute the first several quasi-Lyapunov constants and derive the sufficient and necessary conditions for the origin to be an analytic center. In the last section, we prove that there exist seven small amplitude limit cycles in the neighborhood of the origin.

2. Computation of Quasi-Lyapunov Constants and Determination of Analytic Centers

In this section, we first introduce some definitions, notations, and symbols in order to make the paper compact and clear, followed by an algorithm to obtain the necessary conditions for the third-order nilpotent critical point of system (1) to be an analytic center, and then we present several methods to prove the sufficiency. More details are due to [10, 11].

If the origin of system (1) is a high-order critical point, it is called a th-order critical point when (2) is satisfied, and it could be broken into elementary critical points in the neighborhood of the origin in the complex plane. It is easy to show that the origin of system (1) is a third-order monodromic critical point if and only if and .

Without loss of generality, we assume that Otherwise, by letting , , and taking the transformation , , we obtain the desired form. Under (4), system (1) becomes the following real autonomous planar system:

In the following, we are going to classify the third-order nilpotent critical point.

Theorem 1. For system (5), there exists a series with nonzero convergence radius: such that, by the transformation system (5) is reduced to the following Liénard equation: where . In addition, the origin of system (5) is a center if and only if for all .

The following two definitions are taken from [10].

Definition 2. (1) If , then the origin of system (5) is called a third-order nilpotent critical point of zero-class.
(2) If and there exists a positive integer , such that , but , then the origin of system (5) is called a third-order nilpotent critical point of -class.
(3) If and for all positive integer , then the origin of system (5) is called a third-order nilpotent critical point of -class.

Definition 3. Let , be continuous and bounded functions with respect to and , , . If, for any integer , there exist continuous and bounded functions of and all : , such that then we say that and are equivalent, denoted by . If, for any integer , we have , we say that the sequences of functions and are equivalent, denoted by .

Remark 4. It is easy to see from Definition 3 that the following conclusions hold.(1)The equivalence relationship of two sequences of functions is self-reciprocal, symmetric, and transmissible.(2)If, for some integer , , then when , we have .(3)The relationship implies that .

The following three theorems were proved in [10].

Theorem 5. If the origin of system (5) is -class, then the origin of system (5) is a center if and only if there is an inverse integrating factor .
If the origin of system (5) is -class, then the origin of system (5) is a center if and only if for any natural number there exists an inverse integrating factor , where is a power series given by

Theorem 6. For any natural number and a given number sequence terms with the coefficients satisfying of the formal series (10) can be constructed successively such that that is, where .

Theorem 7. For , in (11) and (12), is uniquely determined by the recursive formula For , is uniquely determined by the recursive formula where In (15), one sets

It follows from Theorems 57 that if the origin of system (5) is a center with -class or -class, then, by choosing , such that

The following conclusion holds (see Theorem 3.4 of [10]).

Theorem 8. If the origin of system (5) is -class, then when the origin of system (5) is a center, in a neighborhood of the origin, system (5) has an analytic inverse integrating factor and an analytic first integral given by

According to [11], we have the following lemma.

Lemma 9. For system (5), if there exists a first integral which is the power series (20) in the neighborhood of the origin, then the origin of system (5) is a center.

By Theorem 1, if the origin of system (5) is a center, there exist analytic transformations in the neighborhood of the origin such that system (5) can be transformed into the Liénard equations:

The vector field defined in (21) is symmetrical with respect to the -axis. Further, by the transformation system (21) is reduced to In the -phase plane, the origin of system (23) is an elementary critical point (focus or center). Obviously, we have the following.

Lemma 10. If there exists a power series in , satisfying then can be written as a power series in , ; namely, .

Then, we have the following theorem.

Theorem 11. The origin of system (5) is an analytic center if and only if the origin of system (5) is a center and -class; namely, the origin of system (5) is a center and for any natural number , .

Corollary 12. If , the origin of system (5) is not an analytic center.

Theorem 13. The origin of system (5) is an analytic center if and only if there exists an analytic first integral in the neighborhood of the origin, which is the power series (23).

From Theorems 8 and 13, we further have the following.

Theorem 14. The origin of system (5) is an analytic center if and only if in the neighborhood of the origin of system (5) there exists an analytic inverse integrating factor

Similarly, by Theorems 5 and 11, we have the following.

Theorem 15. The origin of system (5) is an analytic center, if and only if, for any natural number , there exists an inverse integrating factor , where is the power series (10).

Moreover, Theorems 7 and 15 imply the following.

Theorem 16. The origin of system (5) is an analytic center, if and only if, for any natural number , there exists a power series satisfying in (12) for all .

Therefore, a new method of determining analytic nilpotent center for a given system has been obtained in Theorem 16.

Theorem 17. If the origin of system (5) is a nilpotent center and system (5) is symmetric with respect to the origin, namely, then the origin of system (5) is an analytic center.

Theorem 18. If system (5) is symmetric with respect to the -axis, then the origin of system (5) is an analytic center.

Remark 19. If system (5) is symmetric with respect to the -axis, then the origin of system (5) may not be an analytic center. For example, system (21) is symmetric with respect to the -axis, but the origin is an analytic center if and only if for all .

Eventually, by Theorem 11, we have the following.

Theorem 20. The origin of system (5) is an analytic center if and only if system (5) can be changed into by the analytic transformation (7).

3. Analytic Center Conditions

In this section we will derive conditions for the origin of system (3) to be an analytic center.

After the change and renaming with , system (3) takes the form

According to Theorem 7, we have the following.

Lemma 21. Assume that is a natural number. One can derive a power series (10) for system (29) under which (12) is satisfied, where In addition, for any natural numbers , , is given by the following recursive formula: and, for any natural number , is given by the following recursive formula:

Applying Lemma 21 and computing with Mathematica, we have

From , the analytic center problem of system (29) can be broken down into three cases: (1) , (2) , and (3) .

Case 1. Consider .
In this case, further calculation gives the following.

Theorem 22. For system (29), the first three quasi-Lyapunov constants at the origin are given by In the above expression of , one has already let , .

From Theorem 22, we obtain the following assertion.

Theorem 23. For system (29), all the quasi-Lyapunov constants at the origin are zero if and only if the first three quasi-Lyapunov constants at the origin are zero; that is, one of the following two conditions holds: Relevantly, both conditions are the analytic center conditions of the origin.

Proof. When condition (35) is satisfied, system (29) can be brought to whose vector field is symmetric with respect to the origin.
When condition (36) is satisfied, system (29) can be brought to which is Hamiltonian and possesses the analytic first integral

Case 2. Consider .
In this case, yields that Under this condition, system (29) becomes whose vector field is symmetric with respect to axis . Thus, system (29) has an analytic center at the origin when condition (40) is satisfied.
Obviously, condition (35) is a special case of condition (40).

Case 3. Consder .
In this case, we have the following.

Proposition 24. For system (29), one can determine successively the terms of the formal series , such that where is the th quasi-Lyapunov constant at the origin of system (29), .

Further calculation gives the following.

Theorem 25. For system (29), the first seven quasi-Lyapunov constants at the origin are given by In the above expression of , we have already let , .

, in expression (43) vanish if and only if one of conditions (36) and (40) holds.

Therefore, we see from Cases 13 the following.

Corollary 26. The origin of system (29) ((3)) is an analytic center if and only if condition (36) or (40) holds.

Actually, when condition (36) holds, system (3) goes over to which is Hamiltonian and possesses the analytic first integral

When condition (40) holds, system (3) goes over to whose vector field is symmetric with respect to axis .

4. Multiple Bifurcation of Limit Cycles

In this section, we are going to establish the conditions for to be at most a weak focus of order seven, and, more, we will prove that the perturbed system of (29) can generate seven limit cycles enclosing an elementary node at the origin.

From the fact , , we obtain the following.

Theorem 27. For system (29), the origin is a weak focus of order seven if and only if

Consider the perturbed system of (29):

In order to get seven limit cycles, we only need to show that, when condition (47) holds, the Jacobian of the function group with respect to the parameter group does not equal zero. A direct calculation gives that

The above discussions indicate the following.

Theorem 28. If the origin of system (29) ((3)) is a weak focus of order seven, for , making a small perturbation to the coefficient group , then, for system (48), in a small neighborhood of the origin, there exist exactly seven small amplitude limit cycles enclosing the origin , which is an elementary node.

Example 29. Take where is an arbitrary nonzero real constant.
Straightforward computations by using Theorem 25 give the first seven quasi-Lyapunov constants at the origin of system (48): Then, for , system (48) has seven limit cycles in a small neighborhood of the origin, where , .

Conflict of Interests

The authors declare that they have no conflict of interests.

Acknowledgments

This work is supported in part by the National Nature Science Foundation of China (11101126 and 11261010) and the Scientific Research Foundation for Doctoral Scholars of HAUST (09001524).