# Existence of Center for Planar Differential Systems with Impulsive Perturbations

Authors: Dengguo Xu
Journal: Advances in Mathematical Physics
 (2015)

---

## 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.


---

## Body

## 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 Section2, 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 matrixM will be denoted by M
T and the set of positive integer numbers by N
+. The symbols sqrt
(
x
) and exp
⁡
(
x
) denote x and e
x, respectively.


## 2. Preliminaries

Consider the following planar differential systems with an impulsive perturbation:(1)
x
˙
=
f
x
,
x
∉
T
0
,
Δ
x
=
B
x
,
x
∈
T
0
,where x
=
[
x
1
,
x
2
]
T
∈
R
2; B is a known 2
×
2 constant matrix; f
(
x
)
:
R
2
→
R
2 is a real analytical vector function in the neighborhood of the origin such that f
(
0
)
=
0; and T
0 is a subset of R
2. The phase point of (1) moves between two consecutive intersections with the set T
0 along the trajectories of (1). When the solution to (1) meets the set T
0 at the moment τ, the point x
(
τ
) has a jump Δ
x
=
B
x, and x
(
τ
+
)
=
x
(
τ
)
+
B
x.The following assumptions will be imposed throughout the paper.(
H
1
)
SetT
0
=
⋃
i
=
1
p
‍
s
i
,
p
∈
N
+, where s
i are half-lines starting at the origin in the planar coordinate.
(
H
2
)
The Jacobian matrix off
(
x
) has a pair of conjugate complex roots at the origin. The roots are denoted as α
±
β
j, β
>
0, where j is complex unit.
(
H
3
)
SetB
=
b
11
b
12
b
21
b
22; we denote s
i
+
=
[
I
+
B
]
s
i
,
i
=
1,2
,
…
,
p. Let θ
i and γ
i be angles s
i and s
i
+ with respect to the positive half-line of the first coordinate axis, and 0
<
θ
1
<
γ
1
<
⋯
<
θ
p
<
γ
p
<
2
π.Next, one must make the transformationx
1
=
r
cos
⁡
θ, x
2
=
r
sin
θ. Then (1) is transformed, in polar coordinates, to the following differential system that possesses impulsive attributes:(2a)
r
˙
t
=
α
r
+
r
R
r
,
θ
,
θ
˙
t
=
β
+
Q
r
,
θ
θ
≠
θ
i
,
(2b)
Δ
r
=
k
θ
r
,
Δ
θ
=
Φ
θ
θ
=
θ
i
,where R
(
r
,
θ
) and Q
(
r
,
θ
) are at least first-order term in r with coefficients depending on sin
θ or cos
⁡
θ. When θ
=
θ
i, i
=
1,2
,
…
,
p, the point (
θ
i
,
r
i
) has jumps Δ
θ
i
=
θ
i
+
-
θ
i and Δ
r
i
=
r
i
+
-
r
i; and (3)
k
θ
=
cos
⁡
⁡
θ
+
b
11
cos
⁡
⁡
θ
+
b
12
sin
⁡
θ
2
+
sin
⁡
θ
+
b
21
cos
⁡
⁡
θ
+
b
22
sin
⁡
θ
2
1
/
2
-
1
,
Φ
θ
=
arctan
⁡
1
+
b
22
sin
⁡
θ
+
b
21
cos
⁡
⁡
θ
b
12
sin
⁡
θ
+
1
+
b
11
cos
⁡
⁡
θ
-
θ
.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 smallr
0
>
0, let r
(
θ
,
θ
0
,
r
0
) denote the solution to a system (2a) and (2b) which starts at the point (
θ
0
,
r
0
). Then the origin of system (1) is called the center if r
(
2
π
+
θ
,
θ
0
,
r
0
)
=
r
(
θ
,
θ
0
,
r
0
) and the focus if r
(
2
π
+
θ
,
θ
0
,
r
0
)
≠
r
(
θ
,
θ
0
,
r
0
). Moreover, the origin is called an asymptotically stable focus of the system (1) if r
(
2
π
+
θ
,
θ
0
,
r
0
)
<
r
(
θ
,
θ
0
,
r
0
) and an unstable focus of the system (1) if r
(
2
π
+
θ
,
θ
0
,
r
0
)
>
r
(
θ
,
θ
0
,
r
0
).


## 3. Linear Impulsive Differential Systems

Consider the following linear impulsive differential system:(4)
x
˙
=
A
x
,
x
∉
T
0
,
Δ
x
=
B
x
,
x
∈
T
0
,where A
=
α
-
β
β
α
,
β
>
0, and x
=
[
x
1
,
x
2
]
T. T
0 and the matrix B are the same as those in (
H
1
) and (
H
3
).Using polar coordinates,x
1
=
r
cos
⁡
θ
,
x
2
=
r
sin
θ, the system (4) reduces to(5a)
r
˙
t
=
α
r
,
θ
˙
t
=
β
θ
≠
θ
i
,
(5b)
Δ
r
=
k
θ
r
,
Δ
θ
=
Φ
θ
θ
=
θ
i
,where k
(
θ
) and Φ
(
θ
) are the same as those in (2a) and (2b).Theorem 2.
Letθ
p
+
1
=
2
π
+
θ
0. For θ
∈
θ
i
-
1
,
θ
i
,
i
=
1,2
,
…
,
p
+
1, the solution to system (5a) and (5b) with an initial value of (
θ
0
,
r
0
) can be formulated as (6)
r
θ
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
α
β
θ
-
θ
0
-
∑
j
=
1
i
-
1
Φ
θ
j
+
∑
l
=
1
i
-
1
Δ
r
l
exp
⁡
⁡
α
β
θ
-
θ
l
-
∑
j
=
l
i
-
1
Φ
θ
j
,where θ
0
<
θ
1 and (7)
Δ
r
l
=
r
0
k
θ
l
·
∏
i
=
1
l
-
1
1
+
k
θ
i
exp
⁡
⁡
α
β
θ
l
-
θ
0
-
∑
i
=
1
l
-
1
Φ
θ
i
,
l
=
1,2
,
…
,
i
-
1
.Proof.
Wheni
=
1
,
θ
∈
(
θ
0
,
θ
1
], the solution to system (5a) and (5b) is (8)
r
θ
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
α
t
,
θ
=
θ
0
+
β
t
.It follows that (9)
r
θ
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
α
β
θ
-
θ
0
.This solution intersects with θ
=
θ
1 at point (
θ
1
,
r
1
), where (10)
r
1
=
r
θ
1
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
α
β
θ
1
-
θ
0
.From (5b), we have (11)
θ
1
+
=
θ
1
+
Φ
θ
1
,
r
1
+
=
r
1
+
Δ
r
1
=
r
0
exp
⁡
⁡
α
β
θ
1
-
θ
0
+
Δ
r
1
,
Δ
r
1
=
k
θ
1
r
1
=
r
0
k
θ
1
exp
⁡
⁡
α
β
θ
1
-
θ
0
.So, when i
=
2, θ
∈
(
θ
1
,
θ
2
]. The solution to system (5a) and (5b) is given by (12)
r
θ
,
θ
0
,
r
0
=
r
1
+
exp
⁡
⁡
α
t
=
r
0
exp
⁡
⁡
α
β
θ
1
-
θ
0
+
Δ
r
1
exp
⁡
⁡
α
t
,
θ
=
θ
1
+
+
β
t
=
θ
1
+
Φ
θ
1
+
β
t
,from which we obtain (13)
r
θ
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
α
β
θ
-
θ
0
-
Φ
θ
1
+
Δ
r
1
exp
⁡
⁡
α
β
θ
-
θ
1
-
Φ
θ
1
,which implies that (6) holds when i
=
1 and i
=
2.
Now, suppose that (6) is true when i
=
m. Namely, as θ
∈
(
θ
m
-
1
,
θ
m
], (14)
r
θ
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
α
β
θ
-
θ
0
-
∑
j
=
1
m
-
1
Φ
θ
j
+
∑
l
=
1
m
-
1
Δ
r
l
exp
⁡
⁡
α
β
θ
-
θ
l
-
∑
j
=
l
m
-
1
Φ
θ
j
.This solution intersects with θ
=
θ
m at point (
θ
m
,
r
m
), where (15)
r
m
=
r
θ
m
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
α
β
θ
m
-
θ
0
-
∑
j
=
1
m
-
1
Φ
θ
j
+
∑
l
=
1
m
-
1
Δ
r
l
exp
⁡
⁡
α
β
θ
m
-
θ
l
-
∑
j
=
l
m
-
1
Φ
θ
j
;it gives that from (5b) (16)
θ
m
+
=
θ
m
+
Φ
θ
m
,
r
m
+
=
r
m
+
Δ
r
m
.We can conclude that when θ
∈
(
θ
m
,
θ
m
+
1
], the solution to system (5a) and (5b) is (17)
r
θ
,
θ
0
,
r
0
=
r
m
+
exp
⁡
⁡
α
t
,
θ
=
θ
m
+
+
β
t
;it follows that (18)
r
θ
,
θ
0
,
r
0
=
r
m
+
exp
⁡
⁡
α
β
θ
-
θ
m
+
=
r
m
+
Δ
r
m
exp
⁡
⁡
α
β
θ
-
θ
m
-
Φ
θ
m
=
r
0
exp
⁡
⁡
α
β
θ
-
θ
0
-
∑
j
=
1
m
Φ
θ
j
+
∑
l
=
1
m
Δ
r
l
exp
⁡
⁡
α
β
θ
-
θ
l
-
∑
j
=
l
m
Φ
θ
j
,where (19)
Δ
r
l
=
k
θ
l
r
l
=
r
0
k
θ
l
∏
i
=
1
l
-
1
1
+
k
θ
i
·
exp
⁡
⁡
α
β
θ
l
-
θ
0
-
∑
i
=
1
l
-
1
Φ
θ
i
,
l
=
1,2
,
…
,
i
-
1
.This implies that (6) holds when i
=
m
+
1. Using the reasoning of mathematical induction, we can immediately conclude that (6) is true for any θ
∈
θ
i
-
1
,
θ
i
,
i
=
1,2
,
…
,
p
+
1. The proof is complete.According to Definition1 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.
Denoter
(
θ
,
θ
0
,
r
0
) as the solution to system (5a) and (5b) with an initial value (
θ
0
,
r
0
), and let r
p
+
1
=
r
(
θ
p
+
1
,
θ
0
,
r
0
). Then the origin of system (4) is(i)
a center ifr
p
+
1
=
r
0,
(ii)
an asymptotically stable focus ifr
p
+
1
<
r
0,
(iii)
an unstable focus ifr
p
+
1
>
r
0.Remark 4.
Whenα
≠
0, 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 α
≠
0. 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θ
0
=
0 and θ
∈
(
θ
p
,
θ
p
+
1
], the solution (6) yields(20)
r
θ
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
α
β
θ
-
∑
j
=
1
p
Φ
θ
j
+
∑
l
=
1
p
Δ
r
l
exp
⁡
⁡
α
β
θ
-
θ
l
-
∑
j
=
l
p
Δ
θ
j
=
r
0
exp
⁡
⁡
α
β
θ
-
∑
j
=
1
p
Φ
θ
j
+
∑
l
=
1
p
r
0
k
θ
l
∏
i
=
1
l
-
1
1
+
k
θ
i
×
exp
⁡
⁡
α
β
θ
l
-
∑
i
=
1
l
-
1
Φ
θ
i
×
exp
⁡
⁡
α
β
θ
-
θ
l
-
∑
j
=
l
p
Φ
θ
j
=
r
0
exp
⁡
⁡
α
β
θ
-
∑
j
=
1
p
Φ
θ
j
+
∑
l
=
1
p
r
0
k
θ
l
∏
i
=
1
l
-
1
1
+
k
θ
i
×
exp
⁡
⁡
α
β
θ
-
∑
j
=
1
p
Φ
θ
j
=
r
0
exp
⁡
⁡
α
β
θ
-
∑
j
=
1
p
Φ
θ
j
×
1
+
∑
l
=
1
p
k
θ
l
∏
i
=
1
l
-
1
1
+
k
θ
i
=
r
0
exp
⁡
⁡
α
β
θ
-
∑
j
=
1
p
Φ
θ
j
∏
i
=
1
l
1
+
k
θ
i
,which 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:(21)
x
˙
=
A
x
+
f
x
,
x
∉
T
0
,
Δ
x
=
B
x
,
x
∈
T
0
,where A
=
α
-
β
β
α, β
>
0; x
=
[
x
1
,
x
2
]
T; and f
(
x
)
=
f
1
(
x
1
,
x
2
)
f
2
(
x
1
,
x
2
) is a real analytical vector function in the neighborhood of the origin such that f
(
x
)
=
o
(
x
) and f
i
(
0,0
)
=
0
,
i
=
1,2. T
0 and B are the same as in (4).Let us subject (21) to the transformation x
1
=
r
cos
⁡
θ, x
2
=
r
sin
θ. Then (21) is transformed to(22a)
r
˙
t
=
α
r
+
r
R
r
,
θ
,
θ
˙
t
=
β
+
Q
r
,
θ
θ
≠
θ
i
,
(22b)
Δ
r
=
k
θ
r
,
Δ
θ
=
Φ
θ
θ
=
θ
i
,where k
(
θ
) and Φ
(
θ
) are the same as in (5a) and (5b). Next, (23)
R
r
,
θ
=
1
r
f
1
r
cos
⁡
⁡
θ
,
r
sin
⁡
θ
cos
⁡
⁡
θ
+
1
r
f
2
r
cos
⁡
⁡
θ
,
r
sin
⁡
θ
sin
⁡
θ
,
Q
r
,
θ
=
1
r
f
2
r
cos
⁡
⁡
θ
,
r
sin
⁡
θ
cos
⁡
⁡
θ
-
1
r
f
1
r
cos
⁡
⁡
θ
,
r
sin
⁡
θ
sin
⁡
θ
.It is clear that (22a) yields the following form: (24)
d
r
d
θ
=
α
r
+
r
R
r
,
θ
β
+
Q
r
,
θ
.Using the Taylor expansion, (24) can be expressed as follows: (25)
d
r
d
θ
=
R
1
θ
r
+
R
2
θ
r
2
+
R
3
θ
r
3
+
⋯
,where, as is shown in [25], R
1
(
θ
)
=
α
/
β and R
i
θ
,
i
=
2,3
,
…, is a polynomial in sin
⁡
θ and cos
⁡
θ. Taking into account the continuity of the solution with respect to its initial values and the fact that r
=
0 is a solution to system (25), the solution r
(
θ
,
θ
0
,
r
0
) which starts at point (
θ
0
,
r
0
) exists for θ
∈
(
θ
0
,
θ
1
] when r
=
r
0 is a sufficiently small constant. This solution is analytical in r
0; thus r
(
θ
,
θ
0
,
r
0
) follows the Taylor expansion about r
0: (26)
r
θ
,
θ
0
,
r
0
=
r
11
θ
r
0
+
r
12
θ
r
0
2
+
r
13
θ
r
0
3
+
⋯
;it follows from the initial condition r
(
θ
0
,
θ
0
,
r
0
)
=
r
0 that (27)
r
11
θ
0
=
1
,
r
12
θ
0
=
r
13
θ
0
=
⋯
=
0
.Thus we can conclude that when θ
∈
(
θ
0
,
θ
1
], the solution to a nonlinear impulsive differential system (22a) and (22b) that starts at point (
θ
0
,
r
0
) is given by (26). The solution r
(
θ
,
θ
0
,
r
0
) intersects θ
=
θ
1 at point (
θ
1
,
r
1
), where (28)
r
1
=
r
11
θ
1
r
0
+
r
12
θ
1
r
0
2
+
r
13
θ
1
r
0
3
+
⋯
.From (22b), we have (29)
θ
1
+
=
θ
1
+
Φ
θ
1
,
r
1
+
=
r
1
+
Δ
r
1
=
1
+
k
θ
1
r
1
,using the same reasoning as before, when θ
∈
(
θ
1
,
θ
2
], the solution to system (22a) and (22b), which starts at point (
θ
1
+
,
r
1
+
), can be expressed by the following Taylor series:(30)
r
θ
,
θ
0
,
r
0
=
r
21
θ
r
1
+
+
r
22
θ
r
1
+
2
+
r
23
θ
r
1
+
3
+
⋯
=
r
21
θ
1
+
k
θ
1
r
1
+
r
22
θ
1
+
k
θ
1
2
r
1
2
+
r
23
θ
1
+
k
θ
1
3
r
1
3
+
⋯
=
∑
i
=
1
∞
r
2
i
θ
1
+
k
θ
1
i
r
1
i
.Using the initial condition r
(
θ
1
+
,
θ
0
,
r
0
)
=
r
1
+, one finds that r
2
i
θ
,
i
=
1,2
,
…, satisfys the following conditions: (31)
r
21
θ
1
+
=
1
,
r
22
θ
1
+
=
r
23
θ
1
+
=
⋯
=
0
.The solution (30) intersects with θ
=
θ
2 at point (
θ
2
,
r
2
), where (32)
r
2
=
∑
i
=
1
∞
r
2
i
θ
2
1
+
k
θ
1
i
r
1
i
.From (22b), we have (33)
r
2
+
=
r
2
+
Δ
r
2
=
1
+
k
θ
2
r
2
,
θ
2
+
=
θ
2
+
Δ
θ
2
.When θ
∈
(
θ
2
,
θ
3
], the solution to system (22a) and (22b) starting at point (
θ
2
+
,
r
2
+
) has the following form: (34)
r
θ
,
θ
0
,
r
0
=
r
31
θ
r
2
+
+
r
32
θ
r
2
+
2
+
r
33
θ
r
2
+
3
+
⋯
=
r
31
θ
1
+
k
θ
2
r
2
+
r
32
θ
1
+
k
θ
2
2
r
2
2
+
r
33
θ
1
+
k
θ
2
3
r
2
3
+
⋯
=
∑
i
=
1
∞
r
3
i
θ
1
+
k
θ
2
i
r
2
i
.And r
3
i
(
θ
)
(
i
=
1,2
,
…
) satisfies the following conditions: (35)
r
31
θ
2
+
=
1
,
r
32
θ
2
+
=
r
33
θ
2
+
=
⋯
=
0
.When reasoning via mathematical induction, we can immediately conclude that when θ
∈
θ
j
-
1
,
θ
j
,
j
=
1,2
,
…
,
p
+
1, θ
p
+
1
=
2
π
+
θ
0, the solution to system (22a) and (22b) is given by (36)
r
θ
,
θ
0
,
r
0
=
∑
i
=
1
∞
r
j
i
θ
1
+
k
θ
j
-
1
i
r
j
-
1
i
,where r
1
=
∑
i
=
1
∞
r
1
i
(
θ
1
)
r
0
i and (37)
r
j
-
1
=
∑
i
=
1
∞
r
j
-
1
,
i
θ
j
-
1
1
+
k
θ
j
-
2
i
r
j
-
2
i
,
j
=
3,4
,
…
.And r
j
i
θ
,
i
=
1,2
,
…, satisfies the following conditions: (38)
r
j
1
θ
j
-
1
+
=
1
,
r
j
2
θ
j
-
1
+
=
r
j
3
θ
j
-
1
+
=
⋯
=
0
.Whenθ
∈
θ
j
-
1
,
θ
j
,
j
=
1,2
,
…
,
p
+
1, substituting the solution r
(
θ
,
θ
0
,
r
0
) into (25) and equating coefficient of r
j
-
1
m
,
m
=
1,2
,
…, on both sides in (25), r
j
i
θ
,
i
=
1,2
,
…, can be obtained. This is deduced as follows.Substituting (26) into (25) yields the following equation: (39)
r
11
′
θ
r
0
+
r
12
′
θ
r
0
2
+
r
13
′
θ
r
0
3
+
⋯
=
R
1
θ
r
11
θ
r
0
+
r
12
θ
r
0
2
+
⋯
+
R
2
θ
r
11
θ
r
0
+
r
12
θ
r
0
2
+
⋯
2
+
R
3
θ
r
11
θ
r
0
+
r
12
θ
r
0
2
+
⋯
3
+
⋯
.Equating the coefficient of r
0
m
,
m
=
1,2
,
…, on both sides of (39), one obtains (40)
r
11
′
θ
=
R
1
θ
r
11
θ
,
r
12
′
θ
=
R
1
θ
r
12
θ
+
R
2
θ
r
11
2
θ
,
r
13
′
θ
=
R
1
θ
r
13
θ
+
2
R
2
θ
r
11
θ
r
12
θ
+
R
3
θ
r
11
3
θ
,
⋮By solving these ordinary differential equations with the initial conditions in (27), it follows that (41)
r
11
θ
=
exp
⁡
⁡
α
β
θ
-
θ
0
,
r
12
θ
=
exp
⁡
⁡
α
β
θ
∫
θ
0
θ
R
2
θ
exp
⁡
⁡
α
β
θ
-
2
θ
0
d
θ
,
r
13
θ
=
exp
⁡
⁡
α
β
θ
∫
θ
0
θ
2
R
2
θ
exp
⁡
⁡
α
β
θ
-
3
θ
0
·
∫
θ
0
θ
R
2
θ
exp
⁡
⁡
α
β
θ
d
θ
+
R
3
θ
·
exp
⁡
⁡
α
β
2
θ
-
3
θ
0
d
θ
,
⋮Substituting (30) into (25) yields (42)
r
21
′
θ
1
+
k
θ
1
r
1
+
r
22
′
θ
1
+
k
θ
1
2
r
1
2
+
r
23
′
θ
1
+
k
θ
1
3
r
1
3
+
⋯
=
R
1
θ
·
r
21
θ
1
+
k
θ
1
r
1
+
r
22
θ
1
+
k
θ
1
2
r
1
2
+
r
23
θ
1
+
k
θ
1
3
r
1
3
+
⋯
+
R
2
θ
r
21
θ
1
+
k
θ
1
r
1
+
r
22
θ
1
+
k
θ
1
2
r
1
2
+
r
23
θ
1
+
k
θ
1
3
r
1
3
+
⋯
2
+
R
3
θ
r
21
θ
1
+
k
θ
1
r
1
+
r
22
θ
1
+
k
θ
1
2
r
1
2
+
r
23
θ
1
+
k
θ
1
3
r
1
3
+
⋯
3
+
⋯
.Equating the coefficient of r
1
m
,
m
=
1,2
,
…, 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): (43)
r
21
θ
=
exp
⁡
⁡
α
β
θ
-
θ
0
,
r
22
θ
=
exp
⁡
⁡
α
β
θ
∫
θ
0
θ
R
2
θ
exp
⁡
⁡
α
β
θ
-
2
θ
0
d
θ
,
r
23
θ
=
exp
⁡
⁡
α
β
θ
∫
θ
0
θ
2
R
2
θ
exp
⁡
⁡
α
β
θ
-
3
θ
0
·
∫
θ
0
θ
R
2
θ
exp
⁡
⁡
α
β
θ
d
θ
+
R
3
θ
·
exp
⁡
⁡
α
β
2
θ
-
3
θ
0
d
θ
,
⋮When reasoning via mathematical induction, we can conclude that when θ
∈
θ
j
-
1
,
θ
j
,
j
=
1,2
,
…
,
p
+
1, (44)
r
j
1
θ
=
exp
⁡
⁡
α
β
θ
-
θ
0
,
r
j
2
θ
=
exp
⁡
⁡
α
β
θ
∫
θ
0
θ
R
2
θ
exp
⁡
⁡
α
β
θ
-
2
θ
0
d
θ
,
r
j
3
θ
=
exp
⁡
⁡
α
β
θ
∫
θ
0
θ
2
R
2
θ
exp
⁡
⁡
α
β
θ
-
3
θ
0
·
∫
θ
0
θ
R
2
θ
exp
⁡
⁡
α
β
θ
d
θ
+
R
3
θ
exp
⁡
⁡
α
β
2
θ
-
3
θ
0
d
θ
,
⋮The following theorem yields the solution to system (22a) and (22b).Theorem 6.
Letθ
p
+
1
=
2
π
+
θ
0. For θ
∈
θ
j
-
1
,
θ
j, j
=
1,2
,
…
,
p
+
1, the solution to system (22a) and (22b), which starts at the point (
θ
0
,
r
0
), is given by (45)
r
θ
,
θ
0
,
r
0
=
∑
i
=
1
∞
r
j
i
θ
1
+
k
θ
j
-
1
i
r
j
-
1
i
,where (46)
k
θ
0
=
0
,
r
1
=
∑
i
=
1
∞
r
1
i
θ
1
r
0
i
,
r
j
-
1
=
∑
i
=
1
∞
r
j
-
1
,
i
θ
j
-
1
1
+
k
θ
j
-
2
i
r
j
-
2
i
,
j
=
3,4
,
…
,
p
+
1
.We define the successor function as(47)
ρ
r
0
=
r
2
π
+
θ
0
,
θ
0
,
r
0
-
r
0
.Thus, according to Definition 1, the following theorem holds.Theorem 7.
Assume thatr
0 is sufficiently small, and r
(
θ
,
θ
0
,
r
0
) is the solution to system (22a) and (22b) that starts at point (
θ
0
,
r
0
). The origin is the center of system (21) if ρ
(
r
0
)
=
0. Otherwise, the origin is either an asymptotically stable focus of system (21) under the condition ρ
(
r
0
)
<
0 or an unstable focus under the condition ρ
(
r
0
)
>
0.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:(48a)
x
˙
1
=
0.1
x
1
-
x
2
,
x
˙
2
=
x
1
+
0.1
x
2
x
∉
l
,
(48b)
Δ
x
1
=
3
2
b
0
-
1
x
1
-
1
2
b
0
x
2
,
Δ
x
2
=
1
2
b
0
x
1
+
3
2
b
0
-
1
x
2
x
∈
l
,where x
=
[
x
1
,
x
2
]
T, b
0
>
0 is a real parameter, and the line l is given by equation x
2
=
(
3
/
3
)
x
1
,
x
1
>
0. Therefore θ
1
=
π
/
6, θ
0
<
θ
1
<
2
π and b
11
=
(
3
/
2
)
b
0
-
1, b
12
=
-
(
1
/
2
)
b
0, b
21
=
(
1
/
2
)
b
0, and b
22
=
(
3
/
2
)
b
0
-
1. Using polar coordinates, x
1
=
r
cos
⁡
θ, x
2
=
r
sin
θ, (48a) and (48b) reduce to(49a)
r
˙
t
=
0.1
r
,
θ
˙
t
=
1
θ
≠
θ
i
,
(49b)
Δ
r
=
k
θ
r
,
Δ
θ
=
Φ
θ
θ
=
θ
i
,where (50)
k
θ
=
cos
⁡
⁡
θ
+
b
11
cos
⁡
⁡
θ
+
b
12
sin
⁡
θ
2
+
sin
⁡
θ
+
b
21
cos
⁡
⁡
θ
+
b
22
sin
⁡
θ
2
1
/
2
-
1
=
b
0
-
1
,
Φ
θ
=
arctan
⁡
1
+
b
22
sin
⁡
θ
+
b
21
cos
⁡
⁡
θ
b
12
sin
⁡
θ
+
1
+
b
11
cos
⁡
⁡
θ
-
θ
=
arctan
⁡
3
sin
⁡
θ
+
cos
⁡
⁡
θ
-
sin
⁡
θ
+
3
cos
⁡
⁡
θ
-
θ
.When θ
∈
(
θ
0
,
θ
1
], the solution to (49a) and (49b), which starts at point (
θ
0
,
r
0
), is (51)
r
θ
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
0.1
t
,
θ
=
θ
0
+
t
;that is, (52)
r
θ
,
θ
0
,
r
0
=
r
0
exp
⁡
⁡
0.1
θ
-
θ
0
.This solution intersects with θ
=
θ
1 at point (
θ
1
,
r
1
), where (53)
r
1
=
r
0
exp
⁡
⁡
0.1
θ
1
-
θ
0
.Considering (54)
θ
1
+
=
θ
1
+
Δ
θ
1
=
π
3
,
r
1
+
=
r
1
+
Δ
r
1
=
1
+
k
θ
1
r
1
=
r
0
b
0
exp
⁡
⁡
0.1
θ
1
-
θ
0
;thus, when θ
∈
(
θ
1
,
2
π
), the solution to (49a) and (49b) is (55)
r
θ
,
θ
0
,
r
0
=
r
1
+
exp
⁡
⁡
0.1
t
,
θ
=
θ
1
+
+
t
.Thus, when θ
∈
(
θ
1
,
2
π
], the solution, in polar coordinates, to system (49a) and (49b) that starts at the point (
θ
0
,
r
0
) is given by (56)
r
θ
,
θ
0
,
r
0
=
r
1
+
exp
⁡
0.1
θ
-
θ
1
+
=
r
0
b
0
e
0.1
θ
-
0.1
θ
0
-
π
/
60
.When θ
=
2
π
+
θ
0, r
(
θ
,
θ
0
,
r
0
)
=
r
0
b
0
e
11
π
/
60. Thus, the following proposition follows.Proposition 10.
The origin is the center (focus) of (48a) and (48b) if b
0
=
e
-
11
π
/
60
(
b
0
≠
e
-
11
π
/
60
). Moreover, the origin is a(n) stable (unstable) focus of (48a) and (48b) if b
0
<
e
-
11
π
/
60
(
b
0
>
e
-
11
π
/
60
).Numerical simulations are presented in Figure1. The initial states of system (48a) and (48b) are x
1
(
0
)
=
5, x
2
(
0
)
=
0. Here, we find a center for b
0
=
e
-
11
π
/
60
≈
0.5622 and an unstable focus for b
0
=
e
-
π
/
12
≈
0.7697; see Figures 1(a) and 1(b). For a stable focus, the control parameter is chosen as b
0
=
e
-
π
/
4
≈
0.4559; see Figure 1(c).Figure 1
The center and focus of system (48a) and (48b). (a) The center, where b
0
=
0.5622. (b) An unstable focus, where b
0
=
0.7697. (c) A stable focus, where b
0
=
0.4559.
(a)
(b)
(c)Example 11.
Consider the following nonlinear impulsive differential system:(57a)
x
˙
1
=
-
x
2
,
x
˙
2
=
x
1
+
x
1
2
-
x
2
2
x
∉
l
,
(57b)
Δ
x
1
=
3
2
b
0
-
1
x
1
-
1
2
b
0
x
2
,
Δ
x
2
=
1
2
b
0
x
1
+
3
2
b
0
-
1
x
2
x
∈
l
,where x
=
[
x
1
,
x
2
]
T and the line l is given by equation x
2
=
(
3
/
3
)
x
1
,
x
1
>
0. Thus θ
1
=
π
/
6, 0
=
θ
0
<
θ
1
<
2
π. Applying the transformation in polar coordinates, (57b) is just transformed as (49b) and (57a) yield (58)
d
r
d
θ
=
r
2
sin
⁡
θ
cos
⁡
⁡
2
θ
1
+
r
cos
⁡
⁡
θ
cos
⁡
⁡
2
θ
;using the Taylor expansion, we derive the following differential equation: (59)
d
r
d
θ
=
sin
⁡
θ
cos
⁡
⁡
2
θ
r
2
-
sin
⁡
θ
cos
⁡
⁡
θ
cos
⁡
2
⁡
2
θ
r
3
+
sin
⁡
θ
cos
⁡
2
⁡
θ
cos
⁡
3
⁡
2
θ
r
4
+
⋯
;it gives that (60)
R
1
θ
=
0
,
R
2
θ
=
sin
⁡
θ
cos
⁡
⁡
2
θ
,
R
3
θ
=
-
sin
⁡
θ
cos
⁡
⁡
θ
cos
⁡
2
⁡
2
θ
,
R
4
θ
=
sin
⁡
θ
cos
⁡
2
⁡
θ
cos
⁡
3
⁡
2
θ
,
⋮When θ
∈
[
θ
0
,
θ
1
], the solution to system (59), which starts at the point (
θ
0
,
r
0
), is (61)
r
θ
,
θ
0
,
r
0
=
r
11
θ
r
0
+
r
12
θ
r
0
2
+
r
13
θ
r
0
3
+
⋯
;it follows from (41) that (62)
r
11
θ
=
1
,
r
12
θ
=
1
2
cos
⁡
⁡
θ
-
1
6
cos
⁡
⁡
3
θ
-
1
3
,
r
13
θ
=
1
18
cos
⁡
⁡
θ
20
cos
⁡
5
θ
-
42
cos
⁡
3
θ
+
8
cos
⁡
2
θ
+
27
cos
⁡
θ
-
12
-
1
18
⋮And we have(63)
θ
1
=
π
6
,
k
θ
1
=
b
0
-
1
,
Δ
θ
1
=
π
6
,
θ
1
+
=
π
3
,
r
1
+
=
b
0
r
1
.According to Theorem 6, when θ
∈
(
θ
1
,
2
π
], the solution to system (59) with the impulsive perturbation (49b), which starts at the point (
θ
0
,
r
0
), is (64)
r
θ
,
θ
0
,
r
0
=
∑
i
=
1
∞
r
2
i
θ
b
0
i
r
1
i
=
r
21
θ
b
0
r
1
+
r
22
θ
·
b
0
2
r
1
2
+
r
23
θ
b
0
3
r
1
3
+
⋯
=
r
21
θ
·
b
0
r
11
θ
1
r
0
+
r
12
θ
1
r
0
2
+
r
13
θ
1
r
0
3
+
⋯
+
r
22
θ
·
b
0
2
r
11
θ
1
r
0
+
r
12
θ
1
r
0
2
+
r
13
θ
1
r
0
3
+
⋯
2
+
r
23
θ
·
b
0
3
r
11
θ
1
r
0
+
r
12
θ
1
r
0
2
+
r
13
θ
1
r
0
3
+
⋯
3
+
⋯
.So, we have from (61) and (63) that(65)
r
2
π
,
0
,
r
0
-
r
0,0
,
r
0
=
r
21
2
π
b
0
r
11
θ
1
-
r
11
0
r
0
+
r
21
2
π
b
0
r
12
θ
1
+
r
22
2
π
b
0
2
r
11
2
θ
1
-
r
12
0
r
0
2
+
r
21
2
π
b
0
r
13
θ
1
+
2
r
22
2
π
b
0
2
r
11
θ
1
r
12
θ
1
+
r
23
2
π
b
0
3
r
11
3
θ
1
-
r
13
0
r
0
3
+
⋯
=
b
0
-
1
r
0
+
3
4
-
1
3
b
0
r
0
2
+
9
32
-
3
6
-
1
18
b
0
r
0
3
+
⋯
.Then, the following proposition follows.Proposition 12.
Assume thatr
0 is sufficiently small. If b
0
>
1, the origin is an unstable focus of system (57a) and (57b). If b
0
<
1, the origin is a stable focus of system (57a) and (57b). Otherwise, the coefficient of r
0 is zero and the coefficient of r
0
2 is positive in (65) when b
0
=
1; thus the origin is an unstable focus of system (57a) and (57b).Numerical simulations are displayed in Figure2. The initial states of system (57a) and (57b) are x
1
(
0
)
=
0.2
,
x
2
(
0
)
=
0. Here, we find a stable focus for b
0
=
0.5 and an unstable focus for b
0
=
2.1; see Figures 2(a) and 2(b). We find that the origin of the system (57a) and (57b) is still an unstable focus when b
0
=
1; see Figure 2(c).Figure 2
The focus of system (57a) and (57b). (a) A stable focus, where b
0
=
0.5. (b) An unstable focus, where b
0
=
2.1. (c) An unstable focus, where b
0
=
1.
(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.



---

*Source: 479480-2015-05-05.xml*