Research Article | Open Access

# Existence of Center for Planar Differential Systems with Impulsive Perturbations

**Academic Editor:**Ming Mei

#### Abstract

We present a method that uses successor functions in ordinary differential systems to address the “center-focus” problem of a class of planar systems that have an impulsive perturbation. By deriving solution formulae for impulsive systems, several interesting criteria for distinguishing between the center and the focus of linear and nonlinear planar systems with state-dependent impulsions are established. The conditions describing the stability of the focus of the considered models are also given. The computing methods presented here are more convenient for determining the center of impulsive systems than those in the literature. Numerical examples are given to show the effectiveness of the theoretical results.

#### 1. Introduction

Many differential evolutionary processes that arise in physical, chemical, and biological phenomena are characterized by the fact that they experience an abrupt change of state at certain moments in time. These processes can be modeled using impulsive differential equations. It is known, for example, that many phenomena involve thresholds and exhibit impulsive attributes, including bursting rhythm models (medicine and biology), optimal control models (economics), pharmacokinetics, and frequency modulated systems [1–3]. The theory of impulsive differential systems is much richer than the corresponding theory of systems without impulsive attributes. In general, there are three kinds of impulsive differential systems: (i) systems with impulses at fixed intervals, (ii) systems containing impulses with variable timing, and (iii) autonomously impulsive systems. Over the past few decades, the majority of research has been concentrated on systems with impulses at fixed intervals [4–6], while the other two kinds of impulsive differential systems were relatively less studied due to the difficulties caused by state-dependent impulsion. However, there have been significant developments in the theory of determining the center-focus for differential dynamical systems without impulsive attributes; the reader can refer to the monographs [7, 8]. The center problem consists of distinguishing between a center and a focus at the origin of the system. In the case of a focus, the related problem is to analyze its stability and determine its highest possible order. Since the time of Poincaré [9], who initially defined the notion of the center for a real system of differential equations in the plane, the center-focus problem has been further extended by Liapunov [10], Bendixson [11], and Frommer [12]. Over the past couple decades, various kinds of techniques and algorithms have been developed, and extensive analyses have been consequently produced. Some results about center-focus problem for polynomial differential systems are presented in [13–16]. About the topic of discontinuous systems, there are also a lot of literatures; see, for example, [17–19].

There are many results regarding the bifurcation of impulsive differential systems, for example, [20–22]. However, to the best of our knowledge, the topic of the center-focus for systems with state-dependent impulsion has received less attention. In [23], the center-focus problem and the Hopf bifurcation of dynamical systems with impulses at variable intervals were studied by constructing a Poincaré map. An extension of these results to three-dimensional impulsive differential systems was presented in [24]. However, according to existing techniques, it is not easy to distinguish between the center and the focus in dynamical systems with a state-dependent impulsion. Motivated by this apparent lack of insight, in this paper we devoted our attention to a class of planar impulsive differential systems using successor functions. Our aim is to derive conditions for the existence of the center and focus and to analyze the stability of the focus in dynamical systems with a state-dependent impulsion. At the same time, a computational process will be presented. Due to the complexity caused by an impulse-jump in the systems, it is a challenge to obtain the solution of nonlinear impulsive systems using impulsive jumps of series solutions. In the case of linear impulsive systems, we extend the results in [23] to a more generalized scenario, and a convenient center-focus criterion is given. In the case of nonlinear impulsive systems, we propose a relatively simple method for distinguishing between the center and the focus of the systems, which is completely different from those in the existing literature.

The rest of this paper is organized as follows. In Section 2, the impulsive differential systems are formulated, and new definitions for the center and the focus for systems with a state-dependent impulsion are presented. The criteria for distinguishing between the center and the focus in linear and nonlinear impulsive systems are established by extension via the use of successor functions in Sections 3 and 4, respectively. Two numerical examples are given to demonstrate the effectiveness of the derived results in Section 5. Finally, our conclusions are drawn in Section 6.

The following conventions will be used throughout this paper. The transpose of matrix will be denoted by and the set of positive integer numbers by . The symbols and denote and , respectively.

#### 2. Preliminaries

Consider the following planar differential systems with an impulsive perturbation: where ; is a known constant matrix; is a real analytical vector function in the neighborhood of the origin such that ; and is a subset of . The phase point of (1) moves between two consecutive intersections with the set along the trajectories of (1). When the solution to (1) meets the set at the moment , the point has a jump , and .

The following assumptions will be imposed throughout the paper.Set , where are half-lines starting at the origin in the planar coordinate.The Jacobian matrix of has a pair of conjugate complex roots at the origin. The roots are denoted as , , where is complex unit.Set ; we denote . Let and be angles and with respect to the positive half-line of the first coordinate axis, and .

Next, one must make the transformation , . Then (1) is transformed, in polar coordinates, to the following differential system that possesses impulsive attributes:where and are at least first-order term in with coefficients depending on or . When , , the point has jumps and ; and

Similar to the concepts of the center and focus in the ordinary differential systems [7, 25, 26], the definition of the center and focus for impulsive differential systems (1) is given as follows.

*Definition 1. *For any sufficiently small , let denote the solution to a system (2a) and (2b) which starts at the point . Then the origin of system (1) is called the center if and the focus if . Moreover, the origin is called an asymptotically stable focus of the system (1) if and an unstable focus of the system (1) if .

#### 3. Linear Impulsive Differential Systems

Consider the following linear impulsive differential system: where , and . and the matrix are the same as those in and .

Using polar coordinates, , the system (4) reduces towhere and are the same as those in (2a) and (2b).

Theorem 2. *Let . For , the solution to system (5a) and (5b) with an initial value of can be formulated as **where and *

*Proof. *When , the solution to system (5a) and (5b) is It follows that This solution intersects with at point , where From (5b), we have So, when , . The solution to system (5a) and (5b) is given by from which we obtain which implies that (6) holds when and .

Now, suppose that (6) is true when . Namely, as , This solution intersects with at point , where it gives that from (5b) We can conclude that when , the solution to system (5a) and (5b) is it follows that where This implies that (6) holds when . Using the reasoning of mathematical induction, we can immediately conclude that (6) is true for any . The proof is complete.

According to Definition 1 and Theorem 2, the following theorem holds, and it can distinguish between the center and the focus of a linear impulsive differential system.

Theorem 3. *Denote as the solution to system (5a) and (5b) with an initial value , and let . Then the origin of system (4) is*(i)*a center if ,*(ii)*an asymptotically stable focus if ,*(iii)*an unstable focus if .*

*Remark 4. *When , the origin is only a focus of the differential part in system (4). However, according to Theorem 3, the origin may become a center of the impulsive differential system (4) when . Whether the focus is stable or not has no dependency on the sign of . Therefore, the conclusion in Theorem 3 is actually an extension of the theory of centers and foci in ordinary differential systems.

*Remark 5. *When and , the solution (6) yieldswhich is the solution of the nonperturbed system in [23]. Therefore, Theorem 2 of this paper extends the existing results to a more generalized case.

#### 4. Nonlinear Impulsive Differential Systems

In this section, successor functions are used to study the center-focus problem of a class of nonlinear impulsive differential systems. We will consider the following nonlinear impulsive differential systems:where , ; ; and is a real analytical vector function in the neighborhood of the origin such that and . and are the same as in (4).

Let us subject (21) to the transformation , . Then (21) is transformed towhere and are the same as in (5a) and (5b). Next, It is clear that (22a) yields the following form: Using the Taylor expansion, (24) can be expressed as follows: where, as is shown in [25], and , is a polynomial in and . Taking into account the continuity of the solution with respect to its initial values and the fact that is a solution to system (25), the solution which starts at point exists for when is a sufficiently small constant. This solution is analytical in ; thus follows the Taylor expansion about : it follows from the initial condition that Thus we can conclude that when , the solution to a nonlinear impulsive differential system (22a) and (22b) that starts at point is given by (26). The solution intersects at point , where From (22b), we have using the same reasoning as before, when , the solution to system (22a) and (22b), which starts at point , can be expressed by the following Taylor series:Using the initial condition , one finds that , satisfys the following conditions: The solution (30) intersects with at point , where From (22b), we have When , the solution to system (22a) and (22b) starting at point has the following form: And satisfies the following conditions: When reasoning via mathematical induction, we can immediately conclude that when , , the solution to system (22a) and (22b) is given by where and And , satisfies the following conditions:

When , substituting the solution into (25) and equating coefficient of , on both sides in (25), , can be obtained. This is deduced as follows.

Substituting (26) into (25) yields the following equation: Equating the coefficient of , on both sides of (39), one obtains By solving these ordinary differential equations with the initial conditions in (27), it follows that Substituting (30) into (25) yields Equating the coefficient of , on both sides in (42) and solving the ordinary differential equations using the initial conditions in (31), one obtains the following result, which is similar to that obtained in (41): When reasoning via mathematical induction, we can conclude that when ,

The following theorem yields the solution to system (22a) and (22b).

Theorem 6. *Let . For , , the solution to system (22a) and (22b), which starts at the point , is given by **where *

We define the successor function as Thus, according to Definition 1, the following theorem holds.

Theorem 7. *Assume that is sufficiently small, and is the solution to system (22a) and (22b) that starts at point . The origin is the center of system (21) if . Otherwise, the origin is either an asymptotically stable focus of system (21) under the condition or an unstable focus under the condition .*

*Remark 8. *In order to obtain the solution to nonlinear impulsive differential systems, we use the impulsive jumps of the series solutions in differential systems. This method is fundamentally different from that in [23], and our method plays an important role in distinguishing between the center and the focus of impulsive differential systems.

#### 5. Numerical Examples

In this section, two numerical examples are given to illustrate the theoretical results presented in the previous sections.

*Example 9. *Consider the following linear impulsive differential system:where , is a real parameter, and the line is given by equation . Therefore , and , , , and . Using polar coordinates, , , (48a) and (48b) reduce towhere When , the solution to (49a) and (49b), which starts at point , is that is, This solution intersects with at point , where Considering thus, when , the solution to (49a) and (49b) is Thus, when , the solution, in polar coordinates, to system (49a) and (49b) that starts at the point is given by When , . Thus, the following proposition follows.

Proposition 10. *The origin is the center (focus) of (48a) and (48b) if . Moreover, the origin is a(n) stable (unstable) focus of (48a) and (48b) if .*

Numerical simulations are presented in Figure 1. The initial states of system (48a) and (48b) are , . Here, we find a center for and an unstable focus for ; see Figures 1(a) and 1(b). For a stable focus, the control parameter is chosen as ; see Figure 1(c).

**(a)**

**(b)**

**(c)**

*Example 11. *Consider the following nonlinear impulsive differential system:where and the line is given by equation . Thus , . Applying the transformation in polar coordinates, (57b) is just transformed as (49b) and (57a) yield using the Taylor expansion, we derive the following differential equation: it gives that When , the solution to system (59), which starts at the point , is it follows from (41) that And we haveAccording to Theorem 6, when , the solution to system (59) with the impulsive perturbation (49b), which starts at the point , is So, we have from (61) and (63) thatThen, the following proposition follows.

Proposition 12. *Assume that is sufficiently small. If , the origin is an unstable focus of system (57a) and (57b). If , the origin is a stable focus of system (57a) and (57b). Otherwise, the coefficient of is zero and the coefficient of is positive in (65) when ; thus the origin is an unstable focus of system (57a) and (57b).*

Numerical simulations are displayed in Figure 2. The initial states of system (57a) and (57b) are . Here, we find a stable focus for and an unstable focus for ; see Figures 2(a) and 2(b). We find that the origin of the system (57a) and (57b) is still an unstable focus when ; see Figure 2(c).

**(a)**

**(b)**

**(c)**

#### 6. Conclusion

In this paper, we studied the center-focus problem of planar differential systems that have a state-dependent impulsion. A criterion for distinguishing between the center and the focus in linear and nonlinear impulsive differential systems was established using successor functions. Our results are potentially useful if applied to the theory of impulsive differential systems. An extension of these results to three-dimensional impulsive differential systems is an interesting topic that will be considered in future research.

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The author would like to thank the editors and reviewers for their helpful suggestions. This research is partially supported by the National Natural Science Foundation of China (61463002) and the Yunnan provincial Natural Science Foundation of China (2012FB175).

#### References

- V. Lakshmikantham, D. D. Bainov, and P. S. Simeonov,
*Theory of Impulsive Differential Equations*, vol. 6 of*Series in Modern Applied Mathematics*, World Scientific Publishing, Singapore, 1989. View at: Publisher Site | MathSciNet - D. D. Bainov and P. S. Simeonov,
*Stability Theory of Differential Equations with Impulse Effects: Theory and Applications*, Ellis Horwood, Chichester, UK, 1989. - S. G. Pandit and S. G. Deo,
*Differential Systems Involving Impulses*, Springer, New York, NY, USA, 1982. View at: MathSciNet - V. D. Milman and A. D. Myshkis, “On the stability of motion in the presence of impulses,”
*Siberian Mathematical Journal*, vol. 1, no. 2, pp. 233–237, 1960. View at: Google Scholar | MathSciNet - J. Q. Lu, D. W. Ho, and J. D. Cao, “A unified synchronization criterion for impulsive dynamical networks,”
*Automatica*, vol. 46, no. 7, pp. 1215–1221, 2010. View at: Publisher Site | Google Scholar | MathSciNet - X. N. Liu and L. S. Chen, “Complex dynamics of holling type {II} lotka-volterra predator-prey system with impulsive perturbations on the predator,”
*Chaos, Solitons and Fractals*, vol. 16, no. 2, pp. 311–320, 2003. View at: Publisher Site | Google Scholar | MathSciNet - G. Sansone and R. Conti,
*Non-linear Differential Equations*, Pergamon Press, New York, NY, USA, 1964. View at: MathSciNet - L. Perko,
*Differential Equations and Dynamical Systems*, Springer, New York, NY, USA, 2000. - H. Poincaré, “Mémoire sur les courbes définies par une équation différentielle (II),”
*Journal de Mathématiques Pures et Appliquées*, vol. 8, pp. 251–296, 1882. View at: Google Scholar - M. A. Liapunov,
*Problème Général de la Stabilité du Mouvement*, Princeton University Press, Princeton, NJ, USA, 1947. View at: MathSciNet - I. Bendixson, “Sur les courbes définies par des équations différentielles,”
*Acta Mathematica*, vol. 24, no. 1, pp. 1–88, 1901. View at: Publisher Site | Google Scholar - M. Frommer, “Die integralkurven einer gewöhnlichen differentialgleichung erster ordnung in der umgebung rationaler unbestimmtheitsstellen,”
*Mathematische Annalen*, vol. 99, no. 1, pp. 222–272, 1928. View at: Publisher Site | Google Scholar | MathSciNet - D. M. Wang, “A class of cubic differential system with 6-tuple focus,”
*Journal of Differential Equations*, vol. 87, no. 2, pp. 305–315, 1990. View at: Publisher Site | Google Scholar | MathSciNet - L. S. Chen, Z. Y. Lu, and D. M. Wang, “A class of cubic systems with two centers or two foci,”
*Journal of Mathematical Analysis and Applications*, vol. 242, no. 2, pp. 154–163, 2000. View at: Publisher Site | Google Scholar | MathSciNet - Y. R. Liu and J. B. Li, “Bifurcations of limit cycles and center problem for a class of cubic nilpotent system,”
*International Journal of Bifurcation and Chaos*, vol. 20, no. 8, pp. 2579–2584, 2010. View at: Publisher Site | Google Scholar | MathSciNet - A. P. Sadovskii and T. V. Shcheglova, “Center conditions for a polynomial differential system,”
*Differential Equations*, vol. 49, no. 2, pp. 151–165, 2013. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - H. Giacomini, J. Giné, and J. Llibre, “The problem of distinguishing between a center and a focus for nilpotent and degenerate analytic systems,”
*Journal of Differential Equations*, vol. 227, no. 2, pp. 406–426, 2006. View at: Publisher Site | Google Scholar | MathSciNet - A. Gasull, J. Llibre, V. Mañosas, and F. Mañosas, “The focus-centre problem for a type of degenerate system,”
*Nonlinearity*, vol. 13, no. 3, pp. 699–729, 2000. View at: Publisher Site | Google Scholar | MathSciNet - A. Gasull and J. Torregrosa, “Center-focus problem for discontinuous planar differential equations,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 13, no. 7, pp. 1755–1765, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L. Wang, L. Chen, and J. J. Nieto, “The dynamics of an epidemic model for pest control with impulsive effect,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 3, pp. 1374–1386, 2010. View at: Publisher Site | Google Scholar | MathSciNet - L. N. Qian, Q. S. Lu, Q. G. Meng, and Z. S. Feng, “Dynamical behaviors of a prey-predator system with impulsive control,”
*Journal of Mathematical Analysis and Applications*, vol. 363, no. 1, pp. 345–356, 2010. View at: Publisher Site | Google Scholar | MathSciNet - G.-R. Jiang and Q.-G. Yang, “Periodic solutions and flip bifurcation in a linear impulsive system,”
*Chinese Physics B*, vol. 17, no. 11, pp. 4114–4122, 2008. View at: Publisher Site | Google Scholar - M. U. Akhmet, “Perturbations and Hopf bifurcation of the planar discontinuous dynamical system,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 60, no. 1, pp. 163–178, 2005. View at: Publisher Site | Google Scholar | MathSciNet - M. U. Akhmet and M. Turan, “Bifurcation of three-dimensional discontinuous cycles,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 71, no. 12, pp. e2090–e2102, 2009. View at: Publisher Site | Google Scholar | MathSciNet - F. Dumortier, “Singularities of vector fields on the plane,”
*Journal of Differential Equations*, vol. 23, no. 1, pp. 53–106, 1977. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. D. Bruno,
*Local Methods in Nonlinear Differential Equations*, Springer, Berlin, Germany, 1989.

#### Copyright

Copyright © 2015 Dengguo Xu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.