# TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1

Save this PDF as:
Size: px
Start display at page:

Download "TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1"

## Transcription

1 COMMUNICATIONS ON Website: PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2 and M. Han 1 1 Department of Mathematics Shanghai Jiao Tong University, Shanghai, China Department of Applied Mathematics The University of Western Ontario, London, Ontario, Canada N6A 5B7 Abstract. In this paper, we report the existence of twelve small limit cycles in a planar system with 3rd-degree polynomial functions. The system has Z 2 - symmetry, with a saddle point, or a node, or a focus point at the origin, and two focus points which are symmetric about the origin. It is shown that such a Z 2 -equivariant vector field can have twelve small limit cycles. Fourteen or sixteen small limit cycles, as expected before, cannot not exist. 1. Introduction. One of the well-known mathematical problems is the second part of the 16th Hilbert s problem, which considers the maximal number of limit cycles for a general planar polynomial system. A recent survey article [5] (and more references therein) has comprehensively discussed this problem and reported the recent progress. To start discussion of the 16th Hilbert s problem, assume that the system under consideration is described by the following equations: ẋ = P n (x, y), ẏ = Q n (x, y), (1.1) where the dot denotes the differentiation with respect to time, and P n (x, y) and Q n (x, y) represent n th -degree polynomials of x and y. The problem is to find the maximum number of the limit cycles that system (1.1) may exhibit. In general, this is a very difficult question, in particular, for finding large (global) limit cycles. So far the best result reported in literature for cubic order systems is eleven limit cycles [5, 14]. If we restrict the problem to the neighborhood of isolated fixed points, then the question is reduced to studying degenerate Hopf bifurcations. In the past 50 years, many researchers have considered the local problem and many results have been obtained (e.g., see [1, 4, 6, 7, 8]). Recently, the following cubic order system: ẋ = a x + b y + a 30 x 3 + a 21 x 2 y + a 12 x y 2 + a 03 y 3, ẏ = b x + a y + b 30 x 3 + b 21 x 2 y + b 12 x y 2 + b 03 y 3, (1.2) has been investigated for the existence of limit cycles [2, 3, 10, 11, 12, 13]. Here, the Jacobian of the system evaluated at the equilibrium (x, y) = (0, 0) has eigenvalues a± b. Thus, the origin is either a saddle point when b > a, or a node if a > b. Further, let a 03 = b, b 03 = a and a 12 = a, (1.3) 2000 Mathematics Subject Classification. 34C05, 58F21. Key words and phrases. The 16th Hilbert problem, limit cycle, Hopf bifurcation, focus point. 515

2 516 P. YU AND M. HAN then system (1.2) is symmetric with the origin and has two Hopf-type focus points at (x, y) = (0, ±1). Under the above conditions, ten small limit cycles have been found from system (1.2) when the origin is a saddle point [2]. Recently, it has been shown that twelve small limit cycles exist for a special case of system (1.2) [13]. These results have been reported in several international conferences [3, 10, 11, 12]. Note that the third condition, a 12 = a given in Eq. (1.3), makes the eigenvalues of the system evaluated at the two focus points be a purely imaginary pair. Under the first two conditions of Eq. (1.3), the eigenvalues of the system evaluated at [ (a 12 a) ± (a 12 a) 2 8 (b 2 a 2 +b b 12 a a 12 ) (0, ±1) are given by λ ± = 1 2 from which the real part of the eigenvalues is obtained as v 0 = 1 2 (a 12 a) if (a 12 a) 2 < 8 (b 2 a 2 +b b 12 a a 12 ). The perturbation method based on multiple scales and Maple program developed in [9] was used to compute the focus values (equivalent to finding normal forms of differential equations), which was proved computationally efficient. Based on the focus values, parameter perturbation was carried out to show the existence of the exact small limit cycles [2, 13]. In this paper, we first summarize the results for the special case of twelve limit cycles, and then outline our study on the limit cycles of a general Z 2 -equivariant vector field. We shall show that the maximal number of small limit cycles of such a system is twelve. Fourteen or sixteen small limit cycles, as expected before, cannot exist. Now we briefly outline the result of twelve limit cycle for a special case [3, 10, 11, 12, 13]. Consider system (1.2) and assume b > a, which implies that the origin is a saddle point. In addition to the conditions given in Eq. (1.3), further suppose that the following condition: b 12 = a2 2 b b, (1.4) holds, then system (1.2) is symmetric with the origin and has two Hopf-type focus points at (0, ±1), with the eigenvalues ± i. Then, we have six free parameters: a, b, a 21, b 21, a 30 and b 30, which will be involved in the expressions of the focus values. Since the system is symmetric with the two focus points (0, ±1), we only need to consider one of them, say, (0, 1). Thus, introduce the transformation: x = 2 b u, y = a u v into system (1.2) to obtain ], u = v + 2 (a 21 b a 2 ) u a u v 3 2 v2 + 4 b (a 30 b + a 21 a)u 3 2 (a 21 b a 2 )u 2 v 2 a u v v3, v = u 4 (b 21 b 2 a 21 a b 2a b a 3 + a) u (2 a 2 2 b 2 + 1) u v 8 (b 30 b 3 a 30 a b 2 a 21 a 2 b + b 21 a b 2 a 2 b 2 + a a2 ) u (b 21 b 2 a 21 a b 2a b a 3 + a) u 2 v (2 a 2 b 2 + 1) u v 2, (1.5) whose Jacobian evaluated at the origin is in the real Jordan canonical form. Applying the Maple program [9] to system (1.5) yields the focus values, v i, expressed in terms of the parameters a, b, a 21, b 21, a 30 and b 30. Indeed, we have obtained the explicit symbolic expressions for the six focus values: v i (a, b, a 21, b 21, a 30, b 30 ), i = 1, 2,..., 6. For example, v 1 is given by

3 12 LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM 517 v 1 = 1 2 b [ 2ab (4a 2 4b 2 +3)+3ba 30 +b 21 b (4b 2 4a 21 b 1) + 2a 21 a (2b 2 +2ba 21 4a 2 1) ]. (1.6) Since b 0, we can solve a 30 from the above equation to obtain a 30 = a 30 (b, a 21, a, b 21 ). With this a 30, executing the Maple program results in v 2, given by v 2 = 1 9 b2 (b N + 6 b b D b 30 ), (1.7) where b N and b D are expressed in terms of a, b, a 21, b 21. There are two cases in solving the equation v 2 = 0. Since b 0, so if b D 0 one can uniquely determine b 30 as b 30 = bn 6 b b D. This is generic case to be considered in the next section. When b D = 0, which is referred as a special case, b 30 is not used at this order, but b 21 can be determined, and then simplify v 2. It should be noted that since one more parameter is needed to solve the equation v 2 = 0, thus the remaining parameters b 30, a 21, a and b must be used to solve the four equations v i = 0, i = 2, 3, 4, 5. Hence, the special case gives finite possible solutions, while the generic case can have infinite possible solutions for solving the same number of equations. Now we can use the simplified v 2 to determine a, and then further simplify the expressions v 3, v 4 and v 5 to obtain v 3 = v 4 = v 5 = 2 a b (2 b 2 2 b a 21 1) 4 F F 1, 2 a b (2 b 2 2 b a 21 1) 7 F F 2, a b (2 b 2 2 b a 21 1) 10 F F 3, (1.8) where F i (i = 1, 2, 3) are i th -degree polynomials of b 30. Solving b 30 from F 1 = 0 and then eliminating a 21 from equations F 2 = 0 and F 3 = 0 results in a polynomial equation of b 2 : F 4 = b b b = 0, (1.9) which, in turn, gives a total of four solutions for the special case. We consider one of the solutions, given by the critical values: b = , a 21 = , b 30 = , a = , b 21 = , a 30 = , (1.10) and then b 12 = For these critical values, v The six parameters: a, b, a 21, b 21, a 30 and b 30 can then be perturbed from the critical values to find the exact six small amplitude limit cycles around the origin (u, v) = (0, 0). The results are summarized in the following theorem. The proof can be found in [13]. Theorem 1.1. Given the cubic system (1.2) which is assumed to have a saddle point at the origin and a pair of symmetric fine focus points at (x, y) = (0, 1) and (0, 1). Further suppose a 12 = b 03 = , a 03 = , b 12 = Then if b, a 21, b 30 (b, a 21 ), a(b, a 21 ), b 21 (b, a 21, a) and

4 518 P. YU AND M. HAN a 30 (b, a 21, a, b 21 ) are perturbed as b = b + ɛ 1, a 21 = a 21+ ɛ 2, b 30 = b 30 (b + ɛ 1, a 21 + ɛ 2) + ɛ 3, a = a(b + ɛ 1, a 21+ ɛ 2 ) ɛ 6, b 21 = b 21 (b + ɛ 1, a 21 + ɛ 2, a(b + ɛ 1, a 21 + ɛ 2) ɛ 6 ) ɛ 4, a 30 = a 30 (b + ɛ 1, a 21+ ɛ 2, a(b + ɛ 1, a 21+ ɛ 2 ) ɛ 6, b 21 (b + ɛ 1, a 21 + ɛ 2, a(b + ɛ 1, a 21 + ɛ 2) ɛ 6 ) ɛ 4 ) ɛ 5, (1.11) where 0 < ɛ 6 ɛ 5 ɛ 4 ɛ 3 (ɛ 2, ɛ 1 ) 1, system (1.2) has exact twelve small limit cycles. The notation (ɛ 2, ɛ 1 ) means that ɛ 2 and ɛ 1 are in the same order, with ɛ 2 = (δ + ɛ) ɛ 1 for some δ > 0 and some small ɛ > 0. The normal form of the perturbed system is then given by ṙ = r (v 0 + v 1 r 2 + v 2 r 4 + v 3 r 6 + v 4 r 8 + v 5 r 10 + v 6 r 12 ), (1.12) where 0 < v 0 v 1 v 2 v 3 v 4 v 5 v 6. A numerical example is presented in [11, 12, 13] to illustrate the theoretical result under the following perturbations: ɛ 1 = , ɛ 2 = , ɛ 3 = , ɛ 4 = , ɛ 5 = and ɛ 6 = The six positive roots for r are obtained as: r 1 = , r 2 = , r 3 = , r 4 = , r 5 = and r 6 = The phase portrait of the special case under the above perturbations is shown in Fig. 1, where the two boxes are used to mark the neighborhood of the fine focus points (0, ±1), one of which is enlarged, as shown in Fig. 2 to depict the six small limit cycles near (0, 1), three of them are stable and the other three are unstable. It should be pointed out that one cannot use numerical simulation to get the accurate small limit cycles in the neighborhood of a highly singular point, and has to use certain theoretical approach (like the one presented in this paper) to show the existence of the small limit cycles. In fact, the phase portrait of the system near the focus point (0, 1) is obtained from the normal form, rather than numerical simulations. 2. Twelve Limit Cycles in the Generic Case. In this section, we turn to the generic case in which b D 0 (see Eq. (1.7)), and thus b 30 is uniquely determined as b 30 = bn 6 b b D. For the generic case, there are six free coefficients: a, b, a 21, b 21, a 30 and b 30, which are involved in the expressions of the focus values. Therefore, unlike the special case, it is now possible to have v i = 0, i = 1, 2,..., 6, but v 7 0 This implies that it may be possible to have seven small limit cycles around each of the focus points (0, ±1). We need to find the expression of v 7 and solve the six polynomial equations. In particular, we have to solve the four very large coupled polynomial equations: v i = 0, i = 3, 4, 5, 6. This causes very high complexity in computation. Thus, we first need to simplify Eq. (1.5). To achieve this, let ā 21 = a 21 b a 2, ā 30 = a 30 b + a 21 a, b21 = b 21 b 2 a 21 a b 2a b a 3 + a, b30 = b 30 b 3 a 30 a b 2 a 21 a 2 b + b 21 a b 2 a 2 b 2 + a a2, (2.1)

5 12 LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM x x 1 Figure 1. The phase portrait of system (1.2) having 12 limit cycles, when the origin is a saddle point x x 1 Figure 2. The six limit cycles around the focus point (0, 1): - stable; unstable then Eq. (1.5) becomes u = v + 2 ā 21 u a u v 3 2 v2 + 4 b ā 30 u 3 2 ā 21 u 2 v 2 a u v v3, v = u 4 b 21 u 2 +2(2a 2 2b 2 +1)uv 8 b 30 u 3 + b 21 u 2 v (2a 2 b 2 +1)uv 2. (2.2) Thus, based on Eq. (2.2), we apply the Maple program [9] to find the focus values: v i = v i (a, b, ā 21, b 21, ā 30, b 30 ), i = 1, 2,..., 7. Then, letting v 1 = 0 yield the solution for ā 30 and setting v 2 = 0 results in b 30 = b 30 (a, b, ā 21, b 21 ) = b N 6 b D. The remaining focus values are 4 F F 1 v 3 = 9 [ a(8ā 21 1) 10 b 21 (aā 21 +2a 2 2b 2 + 1) + 8a(a 2 b 2 ) ] 2, v 4 = v 5 = v 6 = F F [ a(8ā 21 1) 10 b 21 (aā 21 +2a 2 2b 2 + 1) + 8a(a 2 b 2 ) ] 3, F F [ a(8ā 21 1) 10 b 21 (aā 21 +2a 2 2b 2 + 1) + 8a(a 2 b 2 ) ] 4, F F [ a(8ā 21 1) 10 b 21 (aā 21 +2a 2 2b 2 + 1) + 8a(a 2 b 2 ) ] 5, (2.3)

6 520 P. YU AND M. HAN where F i = F i (a, b, ā 21, b 21 ), i = 1, 2, 3, 4. Since F 1 is linear in ā 21, one can explicitly solve ā 21 from F 1 = 0 and then obtain the simplified three polynomial equations: Fi = Fi (a, b, b 21 ) = 0, i = 2, 3, 4. These three polynomial equations are coupled and have to be solved simultaneously. Eliminating b from the two equations F2 = 0 and F 3 = 0 yields the equation: F 5 = 160 b a b a2 b21 + a 3 = 0. (2.4) Similarly, eliminating b from another two equations F2 = 0 and F 4 = 0 gives the resultant equation: F 6 = 0, where the lengthy expression of F 6 is omitted. Solving F 5 = F 6 = 0 for a and b 21, gives possibility of finding solutions such that v i = 0, i = 1, 2,..., 6. It should be pointed out that the above variable elimination process do not miss any possible solutions. (It is possible to have extra solutions, and thus one has to verify all solutions using the original equations.) Finally, eliminating b 21 from the two equations: F 5 = F 6 = 0 yields F 7 = a( a a a a a a ) = 0. (2.5) It is obvious that Eq. (2.5) has only one real solution a = 0. But if a = 0, then b21 = 0 which, in turn, results in ā 21 = b 30 = ā 30 = 0, leading to a center. Thus, a must be non-zero. This implies that one cannot find possible non-zero values of a and b 21 such that F 5 = F 6 = 0. Therefore, there is no solution for F 2 = F 3 = F 4 = 0, and so no possible non-trivial solutions exist for a, b, ā 21, b 21, ā 30 and b 30 such that v i = 0, i = 1, 2,..., 6. Hence, fourteen small limit cycles are not possible even for the generic case when the origin is a saddle point, and the maximal number of the small limit cycles for the generic case is twelve. Next, we shall find all the solutions for the twelve limit cycles of the generic case. It has been shown in Section 1 that there exist four solutions for the special case. It should be noted that for each of the four solutions, there are a family of small limit cycles depending upon perturbations. For the generic case, on the other hand, since there is an extra free coefficient to be chosen, there are infinite solutions for the twelve limit cycles. Finding these solutions only requires F 2 = F 3 = 0 (F 4 0), which can be obtained from Eq. (2.4) in a parametric form. For example, numerically solving Eq. (2.4) for b 21 yields three solutions: b21 = a, a, a. Then, for each of these solutions, one finds two solutions for b 2. But for b 21 = a, one solution is negative. Hence, there exist a total of ten solutions for ( b 21, b). However, by checking F 2 = F 3 = 0, only four solutions are satisfied: two for the case when the origin is a saddle point and two for the case when the origin is a node. Note that for the special case discussed in the previous section, solutions exist only for the origin being a saddle point. In the following, we consider one of the solutions when the origin is a node. Let the critical values be: b = a, b 21 = a, ā 21 = a 2,

7 12 LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM x x 1 Figure 3. The phase portrait of system (1.2) having 12 limit cycles, when the origin is a node. b 30 = a2 ( a a 4 ) a 2, ā 30 = ( a2 ). (2.6) Then, similarly like Theorem 1.1 of the special case, we can prove the following theorem for the generic case. Theorem 2.1. Given the cubic system (1.2) which is assumed to have a saddle point or a node at the origin and a pair of symmetric fine focus points at (x, y) = 1+4 a2 (0, ±1). Further suppose a 12 = b 03 = a, a 03 = a, b 12 = 2 b b. Then, for an arbitrarily given a 0, if b, b 21, ā 21 (b, b 21 ), b 30 (b, ā 21, b 21 ) and ā 30 (b, ā 21, b 21 ) are perturbed as b = b + ɛ 1, b21 = b 21+ ɛ 2, ā 21 = ā 21 (b + ɛ 1, b 21 + ɛ 2) + ɛ 3, b30 = b 30 (b + ɛ 1, b 21+ ɛ 2, ā 21 (b + ɛ 1, b 21+ ɛ 2 ) + ɛ 3 ) + ɛ 4, ā 30 = ā 30 (b + ɛ 1, b 21 + ɛ 2, ā 21 (b + ɛ 1, b 21 + ɛ 2) + ɛ 3 ) + ɛ 5, a 12 = a + ɛ 6, (2.7) where 0 < ɛ 6 ɛ 5 ɛ 4 ɛ 3 ( ɛ 2, ɛ 1 ) 1, system (1.2) has exact twelve small limit cycles. The notation ( ɛ 2, ɛ 1 ) means that ɛ 2 and ɛ 1 are in the same order, with ɛ 2 = (δ + ɛ) ɛ 1 for some δ > 0 and some small ɛ > 0. (Note that here ɛ i can be either positive or negative.) To end this section, we present a numerical example for the case when the origin is a node. We choose a = 0.8 and take the following perturbations: ɛ 1 = , ɛ 2 = , ɛ 3 = , ɛ 4 = , ɛ 5 = , ɛ 6 = , under which the six positive roots for r are: r 1 = , r 2 = , r 3 = , r 4 = , r 5 = and r 6 = The phase portrait of this generic case when the origin is a node, under the above perturbations, is shown in Fig. 3, where the marked box near (0, 1) or (0, 1) contains six small limit cycles, like those shown in Fig. 2. But the stabilities of these limit cycles are in the reverse order of those given in Fig. 2.

8 522 P. YU AND M. HAN 3. Other Possible Cases for the Z 2 -equivariant Vector Field. In the previous two sections, we have discussed both special and generic cases which exhibit twelve limit cycles. It has been noted that for the special case, the origin can only be a saddle point, while for the generic case, the origin can be either a saddle point or a node. Further, we have shown that the maximal number of small limit cycles is twelve. Existence of fourteen small limit cycles is not possible. Recall Eq. (1.4) that b 12 was chosen so that the frequency of the limit cycles is 1. If we let b 12 be free, that is, the frequency can be chosen arbitrarily, then we have seven coefficients a, b, a 21, b 21, a 30, b 30 and b 12, free to be chosen in the expressions of the focus values. This seems it might be possible to set v i = 0, i = 1, 2,..., 7, but v 8 0, implying that eight (in total sixteen for the whole system) limit cycles may exist. Also we want to ask: for any Z 2 -equivariant vector field, what is the maximal number of small limit cycles? In this section, we first show that for an arbitrary Z 2 -equivariant vector field, the origin can be a focus point, or a saddle point or a node. The first case has not been considered in the previous two sections. In the following, we shall show that for any of these cases, the maximal number of small limit cycles is twelve. First, we take the standard complex form of the Z 2 -equivariant vector field given by [5]: ż = F 2 (z, z), z = F 2 (z, z), (3.1) where F 2 (z, z) = P (w 1, w 2 ) + i Q(w 1, w 2 ), w 1 = 1 2 (z + z), w 2 = 1 2 i (z z), and F 2 (z, z) = (A 0 + A 1 z 2 + A 2 z 4 ) z + (A 3 + A 4 z 2 + A 5 z 4 ) z + (A 6 + A 7 z 2 ) z 3 + (A 8 + A 9 z 2 ) z 3, up to 5 th order. Here, A j = a j + i b j are complex. For a consistent comparison with system (1.2), we transform Eq (3.1) into the real form and truncated to the 3rd order: ẇ 1 = (a 0 + a 3 ) w 1 (b 0 b 3 ) w 2 + (a 1 + a 4 + a 6 + a 8 ) w1 3 (b 1 b b 6 3 b 8 ) w 2 1 w 2 + (a 1 + a 4 3 a 6 3 a 8 ) w 1 w 2 2 (b 1 b 4 b 6 + b 8 ) w 3 2, ẇ 2 = (b 0 + b 3 ) w 1 + (a 0 a 3 ) w 2 + (b 1 + b 4 + b 6 + b 8 ) w (a 1 a a 6 3 a 8 ) x 2 w 2 + (b 1 + b 4 3 b 6 3 b 8 ) w 1 w (a 1 a 4 a 6 + a 8 ) w 3 2. (3.2) The two eigenvalues of the Jacobian of Eq. (3.2) evaluated at (w 1, w 2 ) = (0, 0) are a 0 ± a b2 3 b2 0. Thus, the origin (0, 0) is a saddle point or a node when a b 2 3 b 2 0 0, while a focus point or a center if a b 2 3 b 2 0 < 0. When a b2 3 b2 0 = 0, the origin is either a node or a double zero singular point. First we consider a b 2 3 b for which the origin is a saddle point or a node. For this case, we introduce the following linear transformations: [ ] [ ] [ ] w1 0 b 0 b 3 = w 2 x a b2 3, (3.3) b2 0 a 3 y into Eq. (3.2), and then let a = a 0, b = ± a b2 3 b2 0, and rename the coefficients of the third order terms in the resulting equations as a 30, a 21, a 12, a 03, b 30, b 21, b 12 and b 03, which finally yields Eq. (1.2). It is easy to see that system (1.2) is invariant under the rotation about the origin with angle π. To consider small limit

9 12 LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM 523 cycles, one needs fine focus point to generate degenerate Hopf bifurcation. Due to the Z 2 -symmetry, it is naturally to best have two symmetric focus points about the origin. Then, if one focus point has N limit cycles, the whole system would have 2N limit cycles, though requiring the two symmetric focus points be Hopf type gives some constraints on the coefficients. If consider one Hopf type focus point at the origin, it may have couple of more free coefficients to choose, but the total number of the limit cycles would be less than that of the case of the two symmetric focus points. This will be discussed next for the case: a b2 3 b2 0 < 0. In addition, noticing that the two eigenvectors of the saddle point (the origin) are (1, 1) and ( 1, 1), so the two focus points can be assumed on the y-axis (or x-axis). Further, by a proper scaling, we can assume that the two focus points are located at (0, ±1), which leads to the conditions a 03 = b and b 03 = a, and making the two focus point be Hopf type requires a 12 = a. These are exact three conditions given by Eq (1.3). Before considering the case a b 2 3 b 2 0 < 0, we want to show that the condition 1+4 a2 given in Eq. (1.4): b 12 = 2 b b, which normalizes the frequency of limit cycles to be 1, is not necessary. In other words, if we assume that the frequency can be selected arbitrarily, then b 12 is free to choose. However, it can been shown that this extra free coefficient does not increase possibility of extra limit cycles. That is, one of the coefficients b 12, a and b (not necessarily b 12 ), can be chosen arbitrarily, which does not affect the solutions of limit cycles. To prove the above statement, let ω be an arbitrary frequency of the limit cycles, then introducing the following parameter scaling and the time scaling: a ω a, b ω b, a 30 ω a 30, a 21 ω a 21, b 30 ω b 30, b 21 ω b 21, τ = ω t, (3.4) into system (1.2) results in the following equations: du dτ = v + 2 ā 21 u a u v 3 2 v2 + 4 b ā 30 u 3 2 ā 21 u 2 v 2 a u v v3, dv dτ = u 4 b 21 u 2 +2(2a 2 2b 2 +1)uv 8 b 30 u 3 +4 b 21 u 2 v (2a 2 b 2 +1)uv 2, (3.5) in which a shifting of the focus point (0, 1) to (0, 0) has been used. The above equation is exactly in the same form of Eq. (2.2), except that the time now is τ. The coefficients ā 21, b 21, ā 30 and b 30 are given by Eq. (2.1). Therefore, based on the above discussion and the results obtained in the previous two sections, we have shown that for the case a 2 3 +b2 3 b2 0 0, the maximal number of small limit cycles that system (1.2) can have is twelve. Next, we turn to the case when a b2 3 b2 0 < 0. For this case, the origin is either a focus point (a 0 0) or a center (a 0 = 0). Similarly, let a = a 0 and b = (a b2 3 b2 0 ) > 0, then the eigenvalues of the Jacobian of the system evaluated at the origin are a ± b i. Applying a linear transformation into Eq. (3.2) and renaming the coefficients of the third order terms in the resulting equations as a 30, a 21, a 12, b 30, b 21, b 12 and b 30 produces ẋ = a x + b y + a 30 x 3 + a 21 x 2 y + a 12 x y 2 + a 03 y 3, ẏ = b x + a y + b 30 x 3 + b 21 x 2 y + b 12 x y 2 + b 03 y 3. (3.6)

10 524 P. YU AND M. HAN Then via the scalings given in Eq. (3.4), one can shift the point (0, 1) to (0, 0) to obtain du dτ = v + 2 ā 21 u a u v 3 2 v2 + 4 b ā 30 u 3 2 ā 21 u 2 v 2 a u v v3, dv dτ = u 4 b 21 u 2 +2(2a 2 +2b 2 +1)uv 8 b 30 u 3 +4 b 21 u 2 v (2a 2 +b 2 +1)uv 2. (3.7) where ā 21 and ā 30 are given in Eq. (2.1), while b 21 and b 30 become b21 = b 21 b 2 a 21 a b + 2a b a 3 + a, b30 = b 30 b 3 a 30 a b 2 a 21 a 2 b + b 21 a b 2 + a 2 b 2 + a a2. (3.8) Based on Eq. (3.7), we can follow the same procedures, described in Sections 1 (for the special case) and 2 (for the generic case), to consider the case when the origin is a focus point. In fact, all the formulas for this case are similar to the case when the origin is a saddle point or a node. Therefore, for this case, we can conclude that the maximal number of small limit cycles is also twelve and they indeed exist. One may think that in addition to the two symmetric fine focus points, making the origin be a center may increase the number of small limit cycles. This requires a = 0, for which one can easily show that all the three fixed points (the origin and the two symmetric points) are centers. Thus, this does not increase the number of limit cycles. Finally, let us consider the origin being a focus point when a = 0 (i.e., a 0 = 0), and investigate the maximal number of small limit cycles which can exist around the origin. Under the condition a = 0, Eq. (3.6) has eight free coefficients a ij s and b ij s, while b = ω > 0 can be chosen arbitrarily, which does not affect the number of limit cycles. In general, eight free coefficients may yield nine small limit cycles. But a careful consideration (following a similar procedure as that given in the previous sections) shows that the maximal number of small limit cycles around the origin is five, much less than twelve for the case of two symmetric focus points. The case a = b = 0, giving rise to a double zero singular point, obviously cannot have twelve small limit cycles. Summarizing the results obtained in this section, we have the following theorem. Theorem 3.1. For the cubic system (1.2), suppose it has property of Z 2 -symmetry, then the maximal number of small limit cycles that the system can exhibit is twelve, i.e., H 2 (3) = 12. Here, as usual, H(n) denotes the Hilbert number of a planar system with nth-degree polynomials. The subscript 2 indicates the Z 2 -equivariant vector field. Finally, a numerical example is given to show the twelve small limit cycles around the two symmetrical focus points when the origin is a focus point. Taking again a = 0.8, and using the following perturbations: ɛ 1 = , ɛ 2 = , ɛ 3 = , ɛ 4 = , ɛ 5 = and ɛ 6 = , yields the six positive roots for r : r 1 = , r 2 = , r 3 = , r 4 = , r 5 = and r 6 = The phase portrait of this case, when the origin is a focus point, is shown in Fig. 4, where the areas near the points (0, ±1) marked by the two boxes, contain the twelve small limit cycles.

11 12 LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM x x 1 Figure 4. The phase portrait of system (1.2) having 12 limit cycles, when the origin is a focus point. 4. Conclusion. A cubic order, planar system with Z 2 -symmetry is studied for the existence of small limit cycles. With the aid of normal form computation for degenerate Hopf bifurcations, it has been shown that such a system can have maximal twelve small limit cycles. The system has two symmetric fine focus points on the y-axis (or x-axis), while the origin can be a saddle point, or a node, or a focus point. For each case, solutions can be classified as special case or generic case. For the special case, the solutions for twelve limit cycles are finite, while for the generic case, it can have infinite solutions. Furthermore, it has been shown that one coefficient (related to the frequency of limit cycles) can be chosen arbitrarily, which does not affect the number of limit cycles. Fourteen or more small limit cycles are not possible for Z 2 -equivariant vector fields. Acknowledgments. This work was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC No. R2686A02) and the National Natural Science Foundation of China (NNSFC No ). REFERENCES [1] N.N. Bautin, On the number of limit cycles which appear with the variation of coefficients from an equilibrium position of focus or center type, Mat. Sbornik (N.S.), 30 (1952), [2] M. Han Y. Lin, and P. Yu, A study on the existence of limit cycles of a planar system with 3rd-degree polynomials, Int. J. Bifurcation & Chaos, 14 (2004), [3] M. Han and P. Yu, Existence of twelve small limit cycles in a cubic order planar system, Shanghai International Symposium on Nonlinear Science and Applications, Shanghai, November 9-13, [4] I.S. Kukles, Necessary and sufficient conditions for the existence of center, Dokl. Akad. Nauk, 42 (1944), [5] J. Li, Hilbert s 16th problem and bifurcations of planar polynomial vector fields, Int. J. Bifurcations & Chaos, 13 (2003), [6] J. Li and Z. Liu, Bifurcation set and limit cycles forming compound eyes in a perturbed Hamiltonian system, Publications Mathematiques, 35 (1991), [7] Y. Liu and J. Li, On the singularity values of complex autonomous differential systems, Sci. in China (Series A), 3 (1989), (in Chinese). [8] K.E. Malkin, Criteria for center of a differential equation, Volg. Matem. Sbornik, 2 (1964), [9] P. Yu, Computation of normal forms via a perturbation technique, J. Sound & Vib., 211 (1998),

12 526 P. YU AND M. HAN [10] P. Yu, Limit cycles in 3rd-order planar system, International Congress of Mathematicians, Beijing, China, August 20-28, [11] P. Yu, Twelve limit cycles in a cubic case of the 16th Hilbert problem, Workshop on Bifurcation Theory and Spatio-Temporal Pattern Formation in PDE, Toronto, Canada, December 11-13, [12] P. Yu, Study of limit cycles in planar systems, International Workshop on Bifurcation Theory and Applications, Shanghai, China, May 23-26, [13] P. Yu and M. Han, Twelve limit cycles in a cubic case of the 16th Hilbert problem, Int. J. Bifurcation & Chaos (to appear). [14] T. Zhang, H. Zang and M. Han, Bifurcations of limit cycles in a cubic system, Chaos, Solitons and Fractals, 20(2004), Received June 2003; revised June address: address:

### THE CENTER-FOCUS PROBLEM AND BIFURCATION OF LIMIT CYCLES IN A CLASS OF 7TH-DEGREE POLYNOMIAL SYSTEMS

Journal of Applied Analysis and Computation Volume 6, Number 3, August 016, 17 6 Website:http://jaac-online.com/ DOI:10.1194/01605 THE CENTER-FOCUS PROBLEM AND BIFURCATION OF LIMIT CYCLES IN A CLASS OF

### A matching pursuit technique for computing the simplest normal forms of vector fields

Journal of Symbolic Computation 35 (2003) 59 65 www.elsevier.com/locate/jsc A matching pursuit technique for computing the simplest normal forms of vector fields Pei Yu, Yuan Yuan Department of Applied

### B5.6 Nonlinear Systems

B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre

### APPPHYS217 Tuesday 25 May 2010

APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag

### Models Involving Interactions between Predator and Prey Populations

Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate

### Chapter 2 Hopf Bifurcation and Normal Form Computation

Chapter 2 Hopf Bifurcation and Normal Form Computation In this chapter, we discuss the computation of normal forms. First we present a general approach which combines center manifold theory with computation

### Stability of Dynamical systems

Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u

### Commun Nonlinear Sci Numer Simulat

Commun Nonlinear Sci Numer Simulat 14 (2009) 4041 4056 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Symbolic computation

### CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS

International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems

### FOUR LIMIT CYCLES FROM PERTURBING QUADRATIC INTEGRABLE SYSTEMS BY QUADRATIC POLYNOMIALS*

International Journal of Bifurcation and Chaos, Vol 22, No 22) 25254 28 pages) c World Scientific Publishing Company DOI: 42/S28274252549 FOUR LIMIT CYCLES FROM PERTURBING QUADRATIC INTEGRABLE SYSTEMS

### Dynamical analysis and circuit simulation of a new three-dimensional chaotic system

Dynamical analysis and circuit simulation of a new three-dimensional chaotic system Wang Ai-Yuan( 王爱元 ) a)b) and Ling Zhi-Hao( 凌志浩 ) a) a) Department of Automation, East China University of Science and

### 7 Planar systems of linear ODE

7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

### USING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS

(Preprint) AAS 12-629 USING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS John E. Hurtado and Andrew J. Sinclair INTRODUCTION Although the solutions with respect to time are commonly sought

### B5.6 Nonlinear Systems

B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A

### Chaos, Solitons & Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena

Caos, Solitons & Fractals 5 (0) 77 79 Contents lists available at SciVerse ScienceDirect Caos, Solitons & Fractals Nonlinear Science, and Nonequilibrium and Complex Penomena journal omepage: www.elsevier.com/locate/caos

### CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 14, 2019, at 08 30 12 30 Johanneberg Kristian

### Additive resonances of a controlled van der Pol-Duffing oscillator

Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University

### Constructing a chaotic system with any number of equilibria

Nonlinear Dyn (2013) 71:429 436 DOI 10.1007/s11071-012-0669-7 ORIGINAL PAPER Constructing a chaotic system with any number of equilibria Xiong Wang Guanrong Chen Received: 9 June 2012 / Accepted: 29 October

### 2 Lecture 2: Amplitude equations and Hopf bifurcations

Lecture : Amplitude equations and Hopf bifurcations This lecture completes the brief discussion of steady-state bifurcations by discussing vector fields that describe the dynamics near a bifurcation. From

### AN EFFICIENT METHOD FOR COMPUTING THE SIMPLEST NORMAL FORMS OF VECTOR FIELDS

International Journal of Bifurcation and Chaos, Vol. 3, No. (2003) 9 46 c World Scientific Publishing Company AN EFFICIENT METHOD FOR COMPUTING THE SIMPLEST NORMAL FORMS OF VECTOR FIELDS P. YU and Y. YUAN

### Boulder School for Condensed Matter and Materials Physics. Laurette Tuckerman PMMH-ESPCI-CNRS

Boulder School for Condensed Matter and Materials Physics Laurette Tuckerman PMMH-ESPCI-CNRS laurette@pmmh.espci.fr Dynamical Systems: A Basic Primer 1 1 Basic bifurcations 1.1 Fixed points and linear

### Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.

Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the

### 5.2.2 Planar Andronov-Hopf bifurcation

138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω

### Lotka Volterra Predator-Prey Model with a Predating Scavenger

Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics

### LIMIT CYCLES AND INVARIANT PARABOLAS FOR AN EXTENDED KUKLES SYSTEM

Miskolc Mathematical Notes HU e-issn 787-243 Vol. 5 (204), No., pp. 29 225 LIMIT CYCLES AND INVARIANT PARABOLAS FOR AN EXTENDED KUKLES SYSTEM BÉLA SZABÓ AND IVÁN SZÁNTÓ Received 3 October, 203 Abstract.

### 2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

### DYNAMICAL SYSTEMS WITH A CODIMENSION-ONE INVARIANT MANIFOLD: THE UNFOLDINGS AND ITS BIFURCATIONS

International Journal of Bifurcation and Chaos c World Scientific Publishing Company DYNAMICAL SYSTEMS WITH A CODIMENSION-ONE INVARIANT MANIFOLD: THE UNFOLDINGS AND ITS BIFURCATIONS KIE VAN IVANKY SAPUTRA

### Problem Set Number 5, j/2.036j MIT (Fall 2014)

Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:

### COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS

International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 9 No. III (September, 2015), pp. 197-210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL

### Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional Dispersive Long Wave Equation

Commun. Theor. Phys. (Beijing, China) 43 (005) pp. 975 98 c International Academic Publishers Vol. 43, No. 6, June 15, 005 Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional

### Bifurcations of Traveling Wave Solutions for a Generalized Camassa-Holm Equation

Computational and Applied Mathematics Journal 2017; 3(6): 52-59 http://www.aascit.org/journal/camj ISSN: 2381-1218 (Print); ISSN: 2381-1226 (Online) Bifurcations of Traveling Wave Solutions for a Generalized

### CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 08, 2018, at 08 30 12 30 Johanneberg Kristian

### CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,

### Computation of the simplest normal forms with perturbation parameters based on Lie transform and rescaling

Journal of Computational and Applied Mathematics 144 (2002) 359 373 wwwelseviercom/locate/cam Computation of the simplest normal forms with perturbation parameters based on Lie transform and rescaling

### ANALYSIS AND CONTROLLING OF HOPF BIFURCATION FOR CHAOTIC VAN DER POL-DUFFING SYSTEM. China

Mathematical and Computational Applications, Vol. 9, No., pp. 84-9, 4 ANALYSIS AND CONTROLLING OF HOPF BIFURCATION FOR CHAOTIC VAN DER POL-DUFFING SYSTEM Ping Cai,, Jia-Shi Tang, Zhen-Bo Li College of

### Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples

### Dynamics at infinity and a Hopf bifurcation arising in a quadratic system with coexisting attractors

Pramana J. Phys. 8) 9: https://doi.org/.7/s43-7-55-x Indian Academy of Sciences Dynamics at infinity and a Hopf bifurcation arising in a quadratic system with coexisting attractors ZHEN WANG,,,IRENEMOROZ

### Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x

### DYNAMICS OF THREE COUPLED VAN DER POL OSCILLATORS WITH APPLICATION TO CIRCADIAN RHYTHMS

Proceedings of IDETC/CIE 2005 ASME 2005 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference September 24-28, 2005, Long Beach, California USA DETC2005-84017

### BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs

BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits: equilibria cycles connecting orbits compact invariant manifolds strange

### A new four-dimensional chaotic system

Chin. Phys. B Vol. 19 No. 12 2010) 120510 A new four-imensional chaotic system Chen Yong ) a)b) an Yang Yun-Qing ) a) a) Shanghai Key Laboratory of Trustworthy Computing East China Normal University Shanghai

### An improved lower bound on the number of limit cycles bifurcating from a quintic Hamiltonian planar vector field under quintic perturbation

An improved lower bound on the number of limit cycles bifurcating from a quintic Hamiltonian planar vector field under quintic perturbation Tomas Johnson Department of Mathematics Uppsala University Box

### 8.1 Bifurcations of Equilibria

1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

### A plane autonomous system is a pair of simultaneous first-order differential equations,

Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

### Hello everyone, Best, Josh

Hello everyone, As promised, the chart mentioned in class about what kind of critical points you get with different types of eigenvalues are included on the following pages (The pages are an ecerpt from

### Characterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddle-node equilibrium point of type 0

Anais do CNMAC v.2 ISSN 1984-82X Characterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddle-node equilibrium point of type Fabíolo M. Amaral Departamento

### Anti-synchronization of a new hyperchaotic system via small-gain theorem

Anti-synchronization of a new hyperchaotic system via small-gain theorem Xiao Jian( ) College of Mathematics and Statistics, Chongqing University, Chongqing 400044, China (Received 8 February 2010; revised

### tutorial ii: One-parameter bifurcation analysis of equilibria with matcont

tutorial ii: One-parameter bifurcation analysis of equilibria with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht February 13, 2018 1 This session

### STABILITY OF 2D SWITCHED LINEAR SYSTEMS WITH CONIC SWITCHING. 1. Introduction

STABILITY OF D SWITCHED LINEAR SYSTEMS WITH CONIC SWITCHING H LENS, M E BROUCKE, AND A ARAPOSTATHIS 1 Introduction The obective of this paper is to explore necessary and sufficient conditions for stability

### Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.

Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the

### Slow Manifold of a Neuronal Bursting Model

Slow Manifold of a Neuronal Bursting Model Jean-Marc Ginoux 1 and Bruno Rossetto 2 1 PROTEE Laboratory, Université du Sud, B.P. 2132, 83957, La Garde Cedex, France, ginoux@univ-tln.fr 2 PROTEE Laboratory,

### Bifurcations of Travelling Wave Solutions for the B(m,n) Equation

American Journal of Computational Mathematics 4 4 4-8 Published Online March 4 in SciRes. http://www.scirp.org/journal/ajcm http://dx.doi.org/.436/jasmi.4.4 Bifurcations of Travelling Wave Solutions for

### Dynamical behaviour of a controlled vibro-impact system

Vol 17 No 7, July 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(07)/2446-05 Chinese Physics B and IOP Publishing Ltd Dynamical behaviour of a controlled vibro-impact system Wang Liang( ), Xu Wei( ), and

### New interpretation of homotopy perturbation method

From the SelectedWorks of Ji-Huan He 26 New interpretation of homotopy perturbation method Ji-Huan He, Donghua University Available at: https://works.bepress.com/ji_huan_he/3/ International Journal of

### 4 Second-Order Systems

4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization

### Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems

University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September

### ZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS

This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]

### Classification of Phase Portraits at Equilibria for u (t) = f( u(t))

Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the

### Example of a Blue Sky Catastrophe

PUB:[SXG.TEMP]TRANS2913EL.PS 16-OCT-2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of

### Effects of Atomic Coherence and Injected Classical Field on Chaotic Dynamics of Non-degenerate Cascade Two-Photon Lasers

Commun. Theor. Phys. Beijing China) 48 2007) pp. 288 294 c International Academic Publishers Vol. 48 No. 2 August 15 2007 Effects of Atomic Coherence and Injected Classical Field on Chaotic Dynamics of

### Elements of Applied Bifurcation Theory

Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Third Edition With 251 Illustrations Springer Introduction to Dynamical Systems 1 1.1 Definition of a dynamical system 1 1.1.1 State space 1 1.1.2

### Kinematics of fluid motion

Chapter 4 Kinematics of fluid motion 4.1 Elementary flow patterns Recall the discussion of flow patterns in Chapter 1. The equations for particle paths in a three-dimensional, steady fluid flow are dx

### Dynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited

Dynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited arxiv:1705.03100v1 [math.ds] 8 May 017 Mark Gluzman Center for Applied Mathematics Cornell University and Richard Rand Dept.

### University of Utrecht

University of Utrecht Bachelor thesis The bifurcation diagram of the second nontrivial normal form of an axially symmetric perturbation of the isotropic harmonic oscillator Author: T. Welker Supervisor:

### 7 Two-dimensional bifurcations

7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed

### ROBUST STABLE NONLINEAR CONTROL AND DESIGN OF A CSTR IN A LARGE OPERATING RANGE. Johannes Gerhard, Martin Mönnigmann, Wolfgang Marquardt

ROBUST STABLE NONLINEAR CONTROL AND DESIGN OF A CSTR IN A LARGE OPERATING RANGE Johannes Gerhard, Martin Mönnigmann, Wolfgang Marquardt Lehrstuhl für Prozesstechnik, RWTH Aachen Turmstr. 46, D-5264 Aachen,

### Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in Complex Networks

International Journal of Difference Equations ISSN 0973-5321, Volume 4, Number 1, pp. 155 163 (2009) http://campus.mst.edu/ijde Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in

### Research Article Adaptive Control of Chaos in Chua s Circuit

Mathematical Problems in Engineering Volume 2011, Article ID 620946, 14 pages doi:10.1155/2011/620946 Research Article Adaptive Control of Chaos in Chua s Circuit Weiping Guo and Diantong Liu Institute

### Autonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.

Autonomous equations Autonomous systems Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. i f i(x 1, x 2,..., x n ) for i 1,..., n As you

### DIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS

WORLD SCIENTIFIC SERIES ON NONLINEAR SCIENCE Series Editor: Leon O. Chua Series A Vol. 66 DIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS Jean-Marc Ginoux Université du Sud, France World Scientific

### Math 4200, Problem set 3

Math, Problem set 3 Solutions September, 13 Problem 1. ẍ = ω x. Solution. Following the general theory of conservative systems with one degree of freedom let us define the kinetic energy T and potential

### ! 1.1 Definitions and Terminology

! 1.1 Definitions and Terminology 1. Introduction: At times, mathematics aims to describe a physical phenomenon (take the population of bacteria in a petri dish for example). We want to find a function

### Centers of projective vector fields of spatial quasi-homogeneous systems with weight (m, m, n) and degree 2 on the sphere

Electronic Journal of Qualitative Theory of Differential Equations 2016 No. 103 1 26; doi: 10.14232/ejqtde.2016.1.103 http://www.math.u-szeged.hu/ejqtde/ Centers of projective vector fields of spatial

### 0 as an eigenvalue. degenerate

Math 1 Topics since the third exam Chapter 9: Non-linear Sstems of equations x1: Tpical Phase Portraits The structure of the solutions to a linear, constant coefficient, sstem of differential equations

### HIGHER ORDER CRITERION FOR THE NONEXISTENCE OF FORMAL FIRST INTEGRAL FOR NONLINEAR SYSTEMS. 1. Introduction

Electronic Journal of Differential Equations, Vol. 217 (217), No. 274, pp. 1 11. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu HIGHER ORDER CRITERION FOR THE NONEXISTENCE

### A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term

ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science

### 18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in

806 Problem Set 8 - Solutions Due Wednesday, 4 November 2007 at 4 pm in 2-06 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state

### An Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson Equation

Commun. Theor. Phys. (Beijing, China) 50 (008) pp. 309 314 c Chinese Physical Society Vol. 50, No., August 15, 008 An Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson

### EE222 - Spring 16 - Lecture 2 Notes 1

EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued

### COMPUTATION OF SIMPLEST NORMAL FORMS OF DIFFERENTIAL EQUATIONS ASSOCIATED WITH A DOUBLE-ZERO EIGENVALUE

International Journal of Bifurcation and Chaos, Vol, No 5 (2) 37 33 c World Scientific Publishing Company COMPUTATION OF SIMPLEST NORMAL FORMS OF DIFFERENTIAL EQUATIONS ASSOCIATED WITH A DOUBLE-ZERO EIGENVALUE

### GEOMETRIC CONFIGURATIONS OF SINGULARITIES FOR QUADRATIC DIFFERENTIAL SYSTEMS WITH TOTAL FINITE MULTIPLICITY m f = 2

Electronic Journal of Differential Equations, Vol. 04 (04), No. 59, pp. 79. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GEOMETRIC CONFIGURATIONS

### GEOMETRIC CONFIGURATIONS OF SINGULARITIES FOR QUADRATIC DIFFERENTIAL SYSTEMS WITH TOTAL FINITE MULTIPLICITY m f = 2

This is a preprint of: Geometric configurations of singularities for quadratic differential systems with total finite multiplicity m f = 2, Joan Carles Artés, Jaume Llibre, Dana Schlomiuk, Nicolae Vulpe,

### Synchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time Systems

Commun. Theor. Phys. (Beijing, China) 48 (2007) pp. 871 876 c International Academic Publishers Vol. 48, No. 5, November 15, 2007 Synchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time

### The Hopf-van der Pol System: Failure of a Homotopy Method

DOI.7/s259--9-5 ORIGINAL RESEARCH The Hopf-van der Pol System: Failure of a Homotopy Method H. G. E. Meijer T. Kalmár-Nagy Foundation for Scientific Research and Technological Innovation 2 Abstract The

### Constructing Chaotic Systems with Total Amplitude Control

International Journal of Bifurcation and Chaos, Vol. 25, No. 10 (2015) 1530025 (14 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127415300256 Constructing Chaotic Systems with Total Amplitude

### Stability Computations for Nilpotent Hopf Bifurcations in Coupled Cell Systems

Stability Computations for Nilpotent Hopf Bifurcations in Coupled Cell Systems M. Golubitsky M. Krupa April 18, 006 Abstract Vanderbauwhede and van Gils, Krupa, and Langford studied unfoldings of bifurcations

### Dynamical systems tutorial. Gregor Schöner, INI, RUB

Dynamical systems tutorial Gregor Schöner, INI, RUB Dynamical systems: Tutorial the word dynamics time-varying measures range of a quantity forces causing/accounting for movement => dynamical systems dynamical

### Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.

Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may

### LIMIT CYCLES FOR A CLASS OF ELEVENTH Z 12 -EQUIVARIANT SYSTEMS WITHOUT INFINITE CRITICAL POINTS

LIMIT CYCLES FOR A CLASS OF ELEVENTH Z 1 -EQUIVARIANT SYSTEMS WITHOUT INFINITE CRITICAL POINTS ADRIAN C. MURZA Communicated by Vasile Brînzănescu We analyze the complex dynamics of a family of Z 1-equivariant

### MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear

### Chaos Control of the Chaotic Symmetric Gyroscope System

48 Chaos Control of the Chaotic Symmetric Gyroscope System * Barış CEVHER, Yılmaz UYAROĞLU and 3 Selçuk EMIROĞLU,,3 Faculty of Engineering, Department of Electrical and Electronics Engineering Sakarya

### 1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

### Bifurcations at in a model for 1:4 resonance

Ergod. Th. & Dynam. Sys. (1997), 17, 899 931 Printed in the United Kingdom c 1997 Cambridge University Press Bifurcations at in a model for 1:4 resonance BERND KRAUSKOPF Vakgroep Wiskunde, Rijksuniversiteit

### Elements of Applied Bifurcation Theory

Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Second Edition With 251 Illustrations Springer Preface to the Second Edition Preface to the First Edition vii ix 1 Introduction to Dynamical Systems

### The 3D restricted three-body problem under angular velocity variation. K. E. Papadakis

A&A 425, 11 1142 (2004) DOI: 10.1051/0004-661:20041216 c ESO 2004 Astronomy & Astrophysics The D restricted three-body problem under angular velocity variation K. E. Papadakis Department of Engineering

### Invariant manifolds of the Bonhoeffer-van der Pol oscillator

Invariant manifolds of the Bonhoeffer-van der Pol oscillator R. Benítez 1, V. J. Bolós 2 1 Dpto. Matemáticas, Centro Universitario de Plasencia, Universidad de Extremadura. Avda. Virgen del Puerto 2. 10600,

### SUB-MANIFOLD AND TRAVELING WAVE SOLUTIONS OF ITO S 5TH-ORDER MKDV EQUATION

Journal of Applied Analysis and Computation Volume 7, Number 4, November 07, 47 430 Website:http://jaac-online.com/ DOI:0.94/0706 SUB-MANIFOLD AND TRAVELING WAVE SOLUTIONS OF ITO S 5TH-ORDER MKDV EQUATION

### A Comparison of Two Predator-Prey Models with Holling s Type I Functional Response

A Comparison of Two Predator-Prey Models with Holling s Type I Functional Response ** Joint work with Mark Kot at the University of Washington ** Mathematical Biosciences 1 (8) 161-179 Presented by Gunog