Dulac-Cherkas function and its application to estimation of limit cycles number for some classes of planar autonomous systems
|
|
- Godfrey West
- 6 years ago
- Views:
Transcription
1 Dulac-Cherkas function and its application to estimation of limit cycles number for some classes of planar autonomous systems Castro Urdiales / September 12-16, 2011 Alexander Grin Yanka Kupala State University of Grodno a joint work with Leonid A. Cherkas Belarusian State University of Informatics and Radioelectronics Klaus R. Schneider Weierstrass Institute for Applied Analysis and Stochastics 1
2 Abstract The talk deals with the estimation of limit cycles number and their localization for real planar autonomous system dx dt = P(x,y), dy dt = Q(x,y), Ω R2, P,Q C 1 (Ω). (1) It is based on the method suggested by L.A. Cherkas (1997) with the use of a Dulac function in the form B = Ψ(x,y) 1/k, 0 k R, Ψ(x,y) C 1 (Ω), satisfying the inequality Φ kψdiv X + Ψ x P + Ψ Q > 0(< 0), (x,y) Ω, X = (P,Q). (2) y In this case Ψ is called a Dulac-Cherkas (DC) function for system (1) in the region Ω. The main subject is limit cycles of systems (1). Main goal is an application DC function to estimation of limit cycles number for some classes of planar autonomous systems. 2
3 Outline 1. DC function, its relationship to Dulac criterion and main properties. 2. Numerical aspects for the construction of DC function. 3. Application of DC function to some classes of system (1). 4. An extension of DC function to system (1) with cylindrical phase space. 5. Conclusions 3
4 DC function, its relationship to Dulac criterion and main properties The Dulac criterion for the nonexistence of a limit cycle: Theorem 0.1 Let D Ω be a simply-connected region. If there exists a function B C 1 (D) such that div(bx) does not change sign in D and is not identically zero on any subregion of D, then (1) has no closed orbit lying entirely in D. Remark 0.2 The proof of Theorem 0.1 is based on Green s Theorem. For B 1, Theorem 0.1 represents the well-known criterion of Bendixson. Definition 0.3 A function B C 1 (D) such that div(bx) has the same sign in some connected region D and is not identically zero on any subregion of D is called a Dulac function to system (1) in D. 4
5 In other words: A function B(x,y) C 1 (Ω) such that the expression F(x,y) = div(bx) = B x P + B y Q + B( P x + Q y ) does not change sign in Ω and such that any curve F(x,y) = 0 does not contain a limit cycle, is called a Dulac function for system (1). The existence of Dulac function can be used to estimate the number of limit cycles of system (1) in some regions. Namely, system (1) has at most p 1 limit cycles in a p connected region Ω if a Dulac function exists. The condition imposed by the criterion on the function F(x,y) implies that the curve B(x,y) = 0 is transversal to the vector field X. Consequently, the limit cycles do not cross it, and their analysis can be performed separately in each of the regions B > 0 and B < 0. However, the criterion itself provides no method for the construction of the function B and for the localization of limit cycles lying in the region Ω. The most widespread forms of B were x a y b, a,b R and e xa y b. 5
6 There is a connection between a DC function and a Dulac function. Lemma 0.4 Let Ω D be connected, let Ψ be a DC function of system (1) in Ω. Then B := Ψ 1/k is a Dulac function in each subregion of Ω where Ψ is positive or negative. The proof is based on the fact that in a subregion of Ω with Ψ > 0 we have div ( Ψ 1/k X) = sign Ψ 1 k Ψ 1 1 Ψ k x P + sign Ψ 1 k Ψ 1 1 Ψ k y Q + Ψ 1 k div X = sign Ψ 1 k Ψ 1 k 1 Φ. (3) 6
7 The curve W defined by plays an important role. W := {(x,y) Ω : Ψ(x,y) = 0} Lemma 0.5 Any trajectory of system (1) meeting the curve W intersects W transversally. Lemma 0.6 The curve W does not contain any equilibrium of system (1). Lemma 0.7 Let W 1 and W 2 be two different smooth local open branches of the curve W such that W 1 W 2 is not connected, that is, W 1 W 2 is empty. Then W 1 and W 2 do not meet. Lemma 0.8 The curve W decomposes the region D in subregions on which Ψ is definite and where the transition from one subregion to an adjacent subregion is connected with a sign change of Ψ. Theorem 0.9 Let Ψ be a DC function of system (1) in D Ω. Then any limit cycle of system (1) which is entirely located in D does not intersect the curve W. [L.A. Cherkas, 1997]; [A. Grin, K. Schneider, 2007] 7
8 The following theorem represents a generalization of Dulac s criterion on the nonexistence of a limit cycle. Theorem 0.10 Let Ω be a connected region and let Ψ be a DC function for system (1) in Ω. Furthermore, we suppose that W decomposes Ω in simply connected subregions Ω i,i = 1,...,l. Then system (1) has no limit cycle in Ω. We note that in Theorem 0.10 the sign of k in the Dulac function plays no role. This is due to the property that the subregions are simply connected and Ψ does not vanish in these regions. The following generalization of Dulac s criterion admits that Ψ can change its sign in simply connected region, but in that case, the sign of k is essential. 8
9 Theorem 0.11 Let Ω be a simply-connected region and let Ψ be a DC function for system (1) in Ω, where k is positive. Then system (1) has no limit cycle in Ω. Idea of the proof: We assume Φ > 0 in Ω. Then we suppose that the decomposition of Ω by W contains a multiply connected region Ω i and that Ω i contains a limit cycle Γ of system (1) surrounding in its interior a closed component W 0 of W. We use the fact that dψ is positive on W dt 0 and consider the region Ψ > 0 or Ψ < 0. In the case Ψ > 0 from (2) integration along the limit cycle Γ yields Γ div f dt = T 0 1 ( kψ Φ dψ dt ) T dt = 0 Φ kψ dt > 0, hence, Γ is orbitally unstable. By the same way we obtain that in any equilibrium point in the interior of Γ the inequality divf = Φ dψ dt kψ = Φ kψ > 0 holds. It leads to the contradiction which proves that there is no limit cycle in the domain Ψ > 0. The case Ψ < 0 can be treated analogously. We get that any trajectory starting near Γ has to cross W 0 but this contradicts the fact that dψ dt > 0. 9
10 The proof of Theorem 0.11 does not work in case k < 0. The following result include the case k < 0. Theorem 0.12 Let Ω D be a connected subset of R 2, and let Ψ be a DC function for system (1) in Ω. Then in each p-connected subdomain Ω of Ω with Ω W Ω, that means any point of the boundary of Ω belongs to W or to Ω, the number of limit cycles is at most p 1 and any limit cycle in Ω is stable (unstable) if k sign (ΨΦ) < 0 (k sign (ΨΦ) > 0). The proof of this theorem follows the same line as in the proof of Theorem Theorem 0.13 Let system (1) has an unique equilibrium point an antisaddle O, divx(o) 0, in a symply-connected region Ω and there is a function Ψ satisfying the condition (2). If in Ω the equation Ψ(x,y) = 0 defines a nest of q embedded ovals then in each of the q 1 doubly-connected subregions Ω i bounded by neighboring ovals system (1) has exactly one limit cycle; and in the Ω in the whole it has at most q limit cycles. 10
11 Remark 0.14 The substantial difference of such approach to estimate the number of limit cycles from the classical one is that we must not localize limit cycles because this localization follows from the topological analysis of the curve Ψ(x,y) = 0. Furthermore, the inequality (2) is linear in Ψ and is considered not on a part of Ω but on the entire region. This approach was exploited by us to estimate the number of limit cycles for some classes of polynomial systems with the help of construction of polynomial DC functions. This approach was exploited also by Gasull A., Giacomini H. A new criterion for controlling the number of limit cycles of some generalized Liénard equations (2002), Gasull A., Giacomini H. Upper bounds for the number of limit cycles through linear differential equations (2006), Gasull A., Giacomini H., Llibre J. New criteria for the existence and non-existence of limit cycles in Liénard differential systems (2008), Gasull A., Giacomini H. Upper bounds for the number of limit cycles of some planar polynomial differential systems (2008). 11
12 Numerical aspects for the construction of DC function We look for a function Ψ(x,y) in the form of a linear combination of given basic functions Ψ = Ψ(x,y,C) = n C j Ψ j (x,y), C j = const, C = (C 1,..., C n ). (4) j=1 Then Φ = Φ(x,y,C) = n C j Φ j (x,y), (5) j=1 and the existence of a linear combination (5) satisfying inequality (2) for a given value of k in the case of a closed bounded domain Ω is equivalent to the inequality max min C j 1 (x,y) Ω Φ(x, y, C) > 0. (6) 12
13 The linear dependence of the function Φ on the additional variables C j, j = 1,...,n, permits one to find the maximin (6) from the corresponding linear programming problem: to find max L, L R, with the restrictions on a grid chosen in the domain Ω. n C j Φ j (x,y) L, C j 1, (x,y) Ω (7) j=1 To solve the problem (7) we can apply e.g. simplex methods and construct the function Ψ in the form of polynomials or in the form of splines. 13
14 Application of DC function to some classes of system (1) If Ω = Ω x = {(x,y) : x [x 1,x N0 ],y R} we construct the function Ψ in the form Ψ(x,y) = n Ψ j (x)y j, Ψ j C 1 (R), (8) j=0 for systems with l 1. dx dt = y, dy l dt = h j (x)y j, j=0 Φ(x,y) h j C 0 (R), m Φ i (x)y i, (9) i=0 where Φ i (x) is a function of the known coefficient functions h 0 (x),...,h l (x), of the unknown coefficient functions Ψ 0 (x),..., Ψ n (x), of their first derivatives Ψ 0(x),..., Ψ n(x), and of k. The highest power m of y in (29) is m = max{n + 1,n + l 1}. To determine the functions Ψ j (x),j = 0,...,n, and the real number k we reduce Φ(x,y) to the following cases 1. Φ(x,y) = Φ 0 (x), Φ i (x) = 0 for i = 1,...,m. (10) 2. Φ(x,y) = Φ 0 (x) + Φ 1 (x)y + Φ 2 (x)y 2, 3. Φ(x,y) = Φ 0 (x) + Φ 2 (x)y 2 + Φ 4 (x)y 4, for l = 1 and l = 2 the relations (10) represent a system of n + 1 linear differential equations to determine the n + 1 functions Ψ j,j = 0,...,n. 14
15 In case l = 1 we have Liènard system and the system (10) can be solved successively by simple quadratures, starting with Ψ n. ẋ = y, ẏ = h 0 (x) + h 1 (x)y, (11) In case l = 2 the system (10) can also be integrated by solving inhomogeneous linear differential equations, starting with Ψ n. We note that the functions Ψ j depend on the parameter k, but we get no restriction on k in the process of solving this system. To fulfill the condition (2) we have to choose k and the integration constants appropriately. x l [x 1,x N0 ], i = 1,N 0. L max, n C i Φi (x l ) L 0, C 1, (12) i=0 15
16 dx dt = y F(x), dy dt = x, (13) F(x) = Tñ+1 (x) = 1 2ñ 1(2xT ñ(x) Tñ 1 (x)), T 0 (x) = 1, T 1 (x) = x, ñ 1, (14) Example We put k = 1 and ñ = 11 in (14) for system (13). If N 0 = 89, 1 x 1. In this case the highest degree of Φ(x) is 100, min = There exist 5 limit cycles. 16
17 In case l = 3 (Kukles system) the first equation of the system (10) is an algebraic equation which determines the constant k uniquely as k = n. The remaining equations represent a system of n + 1 linear differential equations. Its general solution depends on n integration constants which can be used to try to fulfill the relations (2). If this is not possible we have to look for corresponding conditions on the functions h i. The paper [Cherkas L.A., Grin A., 2010] contains two algorithms to construct Φ(x,y) > 0: for odd y p all Φ p (x) = 0, for even y p all Φ p (x) 0 and Φ 0 (x) > 0. 17
18 For Kukles system with P 0 (x) C 1 we look for ẋ = yp 0 (x), ẏ = h 0 (x) + h 1 (x)y + h 2 (x)y 2 + h 3 y 3, (15) n Ψ = Ψ i (x)y n i,n N, (16) i=1 To fulfil Φ(x,y) > 0 we require Φ w (x) < ε for odd y w, ε sufficiently small and Φ v (x) 0 for even y v and Φ 0 (x) > 0. In this case we take all Ψ i (x) in the form Ψ i (x) = m i j=0 C ijx j, C ij R, m i N and solve the linear programming problem L max, m m C j Φ vj (x l ) L > 0, C j Φ wj (x l ) ε < 0, j=0 j=0 the vector C j consists of coefficients C ij from all Ψ i (x) and has dimension m = m m n + n. Example For system (15) with P 0 (x) = 1 + x 2 /6, h 0 (x) = x(1 + x 2 ), h 1 (x) = 1 x 2, h 2 (x) = (x 1)/100, h 3 = 1 we constructed Ψ in the form (16) by using k = 1, n = 3, m 1 = 6, m 2 = 7, m 3 = 8, ε = , N 0 = 100, [ 1.5; 1.5]. For the solution (C,L ) the equation Ψ = 0 defines unique oval and polynomial Φ(x,y) > 0 on the whole plane. It allows to prove the uniqueness of limit cycle globally. 18
19 In case l 4 system (31) consist of n + 1 linear differential equations and l 2 algebraic equations to determine k and the functions Ψ 0,..., Ψ n. Thus, this system has generically no solution. But using reduction to Φ(x,y) = Φ 0 (x) we can construct classes of mentioned systems with DC function in the form Ψ = Ψ 0 (x) + Ψ 1 (x)y + Ψ 2 (x)y 2 that means existence at most one limit cycle if Φ 0 (x) > 0 for all real values of x. [Cherkas L.A., Grin A., Schneider K.R. 2011] 19
20 We have tested the derived approaches to classes of system (1) where P(x,y) and Q(x,y) are periodic functions in the variable x, including the case of the Abel equation and the generalized pendulum equation dy dx = q 0(x) + q 1 (x)y q w (x)y w, (17) dy dx = q 0(x) + q 1 (x)y q w (x)y w. (18) y 20
21 An extension of mentioned approaches to systems with cylindrical phase space We consider systems of two autonomous differential equations on the cylinder Z := S 1 R under the assumption dϕ dt = P(ϕ,y), dy dt = Q(ϕ,y), (19) (A 1 ). P and Q belong to the class C 1 (Ω,R), and are 2π-periodic in the first variable. We denote by Ω the finite region on Z bounded by two closed curves γ 1 and γ 2 which surround the cylinder Z and do not intersect. Two kinds of limit cycles of (19) in Ω: A limit cycle of the first kind is contractible to a point in Ω: it requires the existence of an equilibrium of (19) in Ω. A limit cycle of the second kind surrounds the cylinder Z: it does not need the existence of any equilibrium in Ω. 21
22 Let f be the vector field defined on Ω by system (19), let D be a subregion of Ω. Definition 0.15 A function B C 1 (D,R) with the properties (i). B(ϕ,y) = B(ϕ + 2π,y) (ϕ,y) D, (ii). div(bf) := (BP) + (BQ) = (gradb,f) + Bdivf 0 ( 0) in D, ϕ y where div(bf) vanishes only on a subset of D of measure zero, is called a Dulac function of system (19) in D. Theorem 0.16 Let B be a Dulac function of (19) in the region D Ω. If the boundary D of D is connected and contractible to a point, then (19) has no limit cycle in D. In case that D consists of two closed curves in Ω which do not intersect and which are not contractible to a point then (19) has no limit cycle of the first kind in D and at most one limit cycle of the second kind of (19) in D. 22
23 Definition 0.17 A function Ψ C 1 (D,R) with the properties (i). Ψ(ϕ,y) = Ψ(ϕ + 2π,y) (ϕ,y) D, (ii). The set W := {(ϕ,y) D : Ψ(ϕ,y) = 0} has measure zero. (iii). There is a real number k 0 such that Φ(ϕ,y) := (gradψ,f) + kψdivf 0 ( 0) in D, (20) where the set V := {(ϕ,y) D : Φ(ϕ,y) = 0} has the properties (a). V has measure zero, (b). If Γ is a limit cycle of (19), then it holds Γ V Γ. (iv). (gradψ,f) W 0 (21) is called a Dulac-Cherkas function of (19) in D. Remark 0.18 If the inequalities in (20) hold strictly, then also condition (21) is valid. 23
24 We assume that the bounded region D Ω is formed by the curves 1 and 2 surrounding the cylinder Z. Theorem 0.19 Let Ψ be a DC function of (19) in D. Then it holds: (i). If the set W is empty, then system (19) has at most on limit cycle of the second kind in D. (ii). If the set W contains at least two branches connecting the curves 1 and 2, then system (19) has no limit cycle of the second kind in D. (iii) If the curve W consists in D of s closed branches (ovals) W 1,W 2,...,W s surrounding the cylinder Z and if D contains no equilibrium of (19), then system (19) has at least s 1 but not more than s + 1 limit cycle of the second kind in D. Theorem 0.20 Let Ψ be a DC function of (19) in D. Then the curve W can contain s 2 ovals surrounding the cylinder only if k is negative. From (20) we get such that we have Hence, it holds divf = Φ dψ dt kψ Γ divfds = 1 k Γ Φ Ψ ds. sign divfds = signk. Γ Remark 0.21 In case that k is negative we have divfds < 0. Thus, we can conclude that Γ is a simple limit cycle, that means Γ Γ is hyperbolic. 24
25 The construction of a DC function in a closed finite region Ω Z can be reduced to a linear programming problem L max, n c j φ j (ϕ p,y p ) L 0, c 1, (22) j=1 on a grid of nodes (ϕ p,y p ),p = 1,...,N 0, in the region Ω. As base functions on the cylinder Z we use the functions y i cos(lϕ) and y i sin(lϕ) where i and l are positive integer, such that As an example we consider the system n+1 [ m ( Ψ(ϕ,y) := y i 1 ail cos((l 1)ϕ) + b il sin((l 1)ϕ) )]. (23) dϕ dt = y2 + 1, i=1 dy dt l=1 = (y sinϕ)(y sinϕ)(y 0.3sinϕ). (24) We estimate the number of limit cycles of the second kind for system (24) in the region Ω : π ϕ π, 2π y 2π, using the ansatz (23) with m = 4,n = 4 and k = 0.5 in the corresponding function Φ. We cover Ω with a uniform grid containing 900 nodes. 25
26 Analytical construction of a DC function for a class of systems We consider the vector field X l (ϕ,y) on the cylinder Z defined by the system dϕ dt = y, dy l dt = h j (ϕ)y j, l 1 (25) j=0 under the assumption that the functions h 0,...,h l are continuous and 2π-periodic, and that For a DC function Ψ(ϕ,y) of (25) in Z we make the ansatz h l (ϕ) 0. (26) Ψ(ϕ,y) = n Ψ j (ϕ)y j, n 1 (27) j=0 assuming that the functions Ψ 0,..., Ψ n are continuously differentiable and 2π-periodic. Ψ n (ϕ) 0, (28) We describe an algorithm to determine the functions Ψ j (ϕ) in (27) and the constant k such that the corresponding function Φ(ϕ,y) does not depend on y. For the sequel we represent Φ(ϕ,y) in the form Φ(ϕ,y) m Φ i (ϕ)y i, (29) i=0 the highest power m of y equals m = max{n + 1,n + l 1}. (30) 26
27 Our goal is to determine Ψ j (ϕ),j = 0,...,n, and the number k in such a way that Φ i (ϕ) 0 for i = 1,...,m. (31) Then it holds Φ(ϕ,y) Φ 0 (ϕ) Ψ 1 (ϕ)h 0 (ϕ) + kψ 0 (ϕ)h 1 (ϕ). (32) If we additionally require Φ 0 (ϕ) 0 ( 0) for 0 ϕ 2π (33) and if Φ 0 (ϕ) vanishes only at finitely many points of ϕ, then Ψ is a DC function of (25) in Z. For l 3 system (31) consists of l 2 algebraic equations and n + 1 linear differential equations 0 (n + lk)h l (ϕ)ψ n (ϕ), 0 (k(l 1) + n)h l 1 (ϕ)ψ n (ϕ) + (n 1 + lk)h l (ϕ)ψ n 1 (ϕ), 0 (n 1 + k(l 1))h l 1 (ϕ)ψ n 1 (ϕ) + (n + k(l 2))h l 2 (ϕ)ψ n (ϕ) + (n 2 + lk)h l (ϕ)ψ n 2,... 0 Ψ n(ϕ) + nh 2 (ϕ)ψ n (ϕ) + (n 1)h 3 (ϕ)ψ n 1 (ϕ) (n l)h l+2 (ϕ)ψ n l (ϕ) + 2kh 2 (ϕ)ψ n (ϕ) + 3kh 3 (ϕ)ψ n 1 (ϕ) klh l (ϕ)ψ n l+2 (ϕ),... 0 Ψ 1(ϕ) + (1 + 2k)h 2 (ϕ)ψ 1 (ϕ) + 3kh 3 (ϕ)ψ 0 (ϕ) + (2 + k)h 1 (ϕ)ψ 2 (ϕ) + 3h 0 (ϕ)ψ 3 (ϕ), 0 Ψ 0(ϕ) + 2kh 2 (ϕ)ψ 0 (ϕ) + (k + 1)h 1 (ϕ)ψ 1 (ϕ) + 2h 0 (ϕ)ψ 2 (ϕ). (34) 27
28 Thus, for l 4 from (30) we get m = n + l 1 and system (34) has generically no solution. The algebraic equation Φ m (ϕ) 0 determines the constant k uniquely as k = n l. We consider the system (25) in the case l = 5 dϕ dt = y, dy dt = h 0(ϕ) + h 1 (ϕ)y + h 2 (ϕ)y 2 + h 3 (ϕ)y 3 + h 4 (ϕ)y 4 + h 5 (ϕ)y 5. Our aim is to determine the functions h i, 0 i 5, by means of a DC function Ψ such that (35) has a given number of limit cycles surrounding the cylinder. We use the ansatz Ψ(ϕ,y) = Ψ 0 (ϕ) + Ψ 1 (ϕ)y + Ψ 2 (ϕ)y 2, (36) (35) Ψ 2 (ϕ) 0 for 0 ϕ 2π. (37) Φ(ϕ,y) Φ 0 (ϕ) = Ψ 1 (ϕ)h 0 (ϕ) + kψ 0 (ϕ)h 1 (ϕ). (38) k = 2 5, (39) h 2 (ϕ) = h 3 (ϕ) = h 4 (ϕ) = 5Ψ 1(ϕ)h 5 (ϕ), (40) 2Ψ 2 (ϕ) 5 ( 4Ψ 2 (ϕ) 2Ψ 0 (ϕ) + 3Ψ2 1(ϕ) Ψ 2 (ϕ) ) h 5 (ϕ), (41) 5 [( 9 ) 6Ψ 2 (ϕ) 2 Ψ 0(ϕ) + 3Ψ2 1(ϕ) Ψ1 (ϕ)h 5 (ϕ) ] Ψ 8Ψ 2 (ϕ) Ψ 2 (ϕ) 2(ϕ), (42) 28
29 5 [ h 1 (ϕ) = 8Ψ 2 (ϕ) + Ψ 1(ϕ) + Ψ 1(ϕ)Ψ 2(ϕ) 6Ψ 2 (ϕ) ( 3Ψ 2 0(ϕ) + 3Ψ 0(ϕ)Ψ 2 1(ϕ) 2Ψ 2 (ϕ) Ψ4 1(ϕ) 16Ψ 2 2(ϕ) ) h5 (ϕ) ], Ψ 2 (ϕ) (43) 1 [ h 0 (ϕ) = Ψ 2Ψ 2 (ϕ) 0(ϕ) + 3Ψ 1(ϕ)Ψ ( 1(ϕ) Ψ 2(ϕ) Ψ 2 1 (ϕ) 8Ψ 2 (ϕ) Ψ 2 (ϕ) 16Ψ 2 (ϕ) + 2Ψ ) 0(ϕ) 3 ( Ψ2 0(ϕ) 5Ψ 0(ϕ)Ψ 2 ) 1(ϕ) + 3Ψ4 1(ϕ) Ψ1 (ϕ)h 5 (ϕ) ]. 16Ψ 2 (ϕ) 128Ψ 2 2(ϕ) Ψ 2 2(ϕ) The 2π-periodic continuous function h 5 can be chosen arbitrarily. (44) Φ 0 (ϕ) = 3h ( 5(ϕ) Ψ 3 4Ψ 0(ϕ) + 3Ψ2 0(ϕ)Ψ 2 1(ϕ) 2 2(ϕ) 4Ψ 2 (ϕ) 1 ( + Ψ 1 (ϕ)ψ 2Ψ 2 (ϕ) 0(ϕ) + Ψ 1(ϕ)Ψ 0 (ϕ) 2 ). Ψ 2(ϕ)Ψ 3 1(ϕ) 16Ψ 2 2(ϕ) 3Ψ 2(ϕ)Ψ 1 (ϕ)ψ 0 (ϕ) 4Ψ 2 (ϕ) 3Ψ 0(ϕ)Ψ 4 1(ϕ) 16Ψ 2 2(ϕ) + 3Ψ 1(ϕ)Ψ 2 1(ϕ) 8Ψ 2 (ϕ) ) + Ψ6 1(ϕ) 64Ψ 3 2(ϕ) (45) If we choose Ψ 2 (ϕ) = 1 2, Ψ 1(ϕ) = 1, Ψ 0 (ϕ) = cosϕ 10, (46) then the equation Ψ(ϕ,y) = 0 describes two ovals W 1 := {(ϕ,y) Z : y = cos ϕ}, W 2 := {(ϕ,y) Z : y = cos ϕ}. 29
30 We set h 5 (ϕ) 1 (47) and get h 0 (ϕ) h 1 (ϕ) cos ϕ 2 h 2 (ϕ) cos ϕ cos2 ϕ 2 + sin ϕ, (48) + 15 cos2 ϕ, (49) cos ϕ, (50) h 3 (ϕ) cosϕ, (51) Φ 0 (ϕ) cosϕ 4 h 4 (ϕ) 5. (52) cos2 ϕ 2 + sinϕ 3 cos 3 ϕ Theorem 0.22 System (35) with h i defined by (47) - (52) has a unique limit cycle of second kind on the cylinder Z which is hyperbolic and unstable and located between W 1 and W 2. Now we use the ansatz Ψ(ϕ,y) = Ψ 0 (ϕ) + Ψ 1 (ϕ)y + Ψ 2 (ϕ)y 2 + Ψ 3 (ϕ)y 3 + Ψ 4 (ϕ)y 4, Ψ 4 (ϕ) 0. (53) Φ(ϕ,y) Φ 2 (ϕ)y 2 + Φ 1 (ϕ)y + Φ 0 (ϕ). (54) 30
31 If we introduce the function Ψ : Z R by and setting h 5 = 10 we obtain Ψ(ϕ,y) := y4 2 + ( cosϕ 10 2) y , (55) dϕ dt = y, dy dt = sin ϕ 20 + ( cos ϕ 2 cos2 ϕ 16 (56) ) ( 5 cos ϕ) y + 50 y 3 10y 5 2 ( 87 cosϕ Φ(ϕ,y) cos2 ϕ ( + sin ϕ cos ϕ sin ϕ cos3 ϕ ) y ) y cos ϕ + cos2 ϕ The set W consists of four branches surrounding the cylinder W 1 := {(ϕ,y) Z : y = 2 cos ϕ 10 + cos2 ϕ 40 cos ϕ }, 10 W 2 := {(ϕ,y) Z : y = 2 cos ϕ 10 cos2 ϕ 40 cos ϕ }, 10 W 3 := {(ϕ,y) Z : y = 2 cosϕ 10 cos2 ϕ 40 cos ϕ }, 10 W 4 := {(ϕ,y) Z : y = 2 cos ϕ 10 + cos2 ϕ 40 cos ϕ }. 10 (57) Theorem 0.23 On the cylinder Z system (56) has four limit cycles of second kind which are hyperbolic. 31
32 Conclusions 1. l 1 dx dt = p j (x)y j, j=0 dy l dt = h j (x)y j, h j C 0 (R), p j C 1 (R), j=0 2. dx dt = y εf(x), dy dt = g(x), (58) dx dt = y, dy dt = x + ε(h 0(x) + h 1 (x)y + h 2 (x)y 2 + h 3 (x)y 3 ), (59) 32
33 Thank you! 33
Weierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 198-8 Construction of generalized pendulum equations with prescribed maximum
More informationAbsence of limit cycles for Kukles-type systems
AVANZA. Vol. iv. Fm - Iit, Uacj (014) 5 13. Isbn: 978-607-94-5-3 Absence of limit cycles for Kukles-type systems Osuna, Osvaldo, Rodríguez-Ceballos, Joel, Villaseñor, Gabriel Abstract In this note, we
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationNonlinear Autonomous Dynamical systems of two dimensions. Part A
Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationProblem List MATH 5173 Spring, 2014
Problem List MATH 5173 Spring, 2014 The notation p/n means the problem with number n on page p of Perko. 1. 5/3 [Due Wednesday, January 15] 2. 6/5 and describe the relationship of the phase portraits [Due
More informationProblem 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:
More informationUniqueness of Limit Cycles For Liénard Differential Equations of Degree Four
Uniqueness of Limit Cycles For Liénard Differential Equations of Degree Four Chengzhi Li and Jaume Llibre (Peking University and Universitat Autònoma de Barcelona) Salou, Spain, October 1-5, 2012 Congratulations
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
More informationLecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps
Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit
More informationNew approach to study the van der Pol equation for large damping
Electronic Journal of Qualitative Theor of Differential Equations 2018, No. 8, 1 10; https://doi.org/10.1422/ejqtde.2018.1.8 www.math.u-szeged.hu/ejqtde/ New approach to stud the van der Pol equation for
More informationCenters 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
More informationSynchronization, Chaos, and the Dynamics of Coupled Oscillators. Supplemental 1. Winter Zachary Adams Undergraduate in Mathematics and Biology
Synchronization, Chaos, and the Dynamics of Coupled Oscillators Supplemental 1 Winter 2017 Zachary Adams Undergraduate in Mathematics and Biology Outline: The shift map is discussed, and a rigorous proof
More informationCOEXISTENCE OF ALGEBRAIC AND NON-ALGEBRAIC LIMIT CYCLES FOR QUINTIC POLYNOMIAL DIFFERENTIAL SYSTEMS
Electronic Journal of Differential Equations, Vol. 217 (217), No. 71, pp. 1 7. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu COEXISTENCE OF ALGEBRAIC AND NON-ALGEBRAIC LIMIT
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationfy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f))
1. Basic algebra of vector fields Let V be a finite dimensional vector space over R. Recall that V = {L : V R} is defined to be the set of all linear maps to R. V is isomorphic to V, but there is no canonical
More informationStability Implications of Bendixson s Criterion
Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationNotes for Lecture 10
February 2, 26 Notes for Lecture Introduction to grid-based methods for solving Poisson and Laplace Equations Finite difference methods The basis for most grid-based finite difference methods is the Taylor
More informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More information1. < 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
More informationALGEBRAIC GROUPS. Disclaimer: There are millions of errors in these notes!
ALGEBRAIC GROUPS Disclaimer: There are millions of errors in these notes! 1. Some algebraic geometry The subject of algebraic groups depends on the interaction between algebraic geometry and group theory.
More informationAn extended complete Chebyshev system of 3 Abelian integrals related to a non-algebraic Hamiltonian system
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 4, 2018, pp. 438-447 An extended complete Chebyshev system of 3 Abelian integrals related to a non-algebraic Hamiltonian
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationPERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND
PERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND JOSEP M. OLM, XAVIER ROS-OTON, AND TERE M. SEARA Abstract. The study of periodic solutions with constant sign in the Abel equation
More informationSequences. Chapter 3. n + 1 3n + 2 sin n n. 3. lim (ln(n + 1) ln n) 1. lim. 2. lim. 4. lim (1 + n)1/n. Answers: 1. 1/3; 2. 0; 3. 0; 4. 1.
Chapter 3 Sequences Both the main elements of calculus (differentiation and integration) require the notion of a limit. Sequences will play a central role when we work with limits. Definition 3.. A Sequence
More informationnonlinear oscillators. method of averaging
Physics 4 Spring 7 nonlinear oscillators. method of averaging lecture notes, spring semester 7 http://www.phys.uconn.edu/ rozman/courses/p4_7s/ Last modified: April, 7 Oscillator with nonlinear friction
More informationExercises from other sources REAL NUMBERS 2,...,
Exercises from other sources REAL NUMBERS 1. Find the supremum and infimum of the following sets: a) {1, b) c) 12, 13, 14, }, { 1 3, 4 9, 13 27, 40 } 81,, { 2, 2 + 2, 2 + 2 + } 2,..., d) {n N : n 2 < 10},
More informationElementary ODE Review
Elementary ODE Review First Order ODEs First Order Equations Ordinary differential equations of the fm y F(x, y) () are called first der dinary differential equations. There are a variety of techniques
More informationCommutative Banach algebras 79
8. Commutative Banach algebras In this chapter, we analyze commutative Banach algebras in greater detail. So we always assume that xy = yx for all x, y A here. Definition 8.1. Let A be a (commutative)
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationMASTERS EXAMINATION IN MATHEMATICS SOLUTIONS
MASTERS EXAMINATION IN MATHEMATICS PURE MATHEMATICS OPTION SPRING 010 SOLUTIONS Algebra A1. Let F be a finite field. Prove that F [x] contains infinitely many prime ideals. Solution: The ring F [x] of
More informationA note on linear differential equations with periodic coefficients.
A note on linear differential equations with periodic coefficients. Maite Grau (1) and Daniel Peralta-Salas (2) (1) Departament de Matemàtica. Universitat de Lleida. Avda. Jaume II, 69. 251 Lleida, Spain.
More information4.1 Primary Decompositions
4.1 Primary Decompositions generalization of factorization of an integer as a product of prime powers. unique factorization of ideals in a large class of rings. In Z, a prime number p gives rise to a prime
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationProblem: A class of dynamical systems characterized by a fast divergence of the orbits. A paradigmatic example: the Arnold cat.
À È Ê ÇÄÁ Ë ËÌ ÅË Problem: A class of dynamical systems characterized by a fast divergence of the orbits A paradigmatic example: the Arnold cat. The closure of a homoclinic orbit. The shadowing lemma.
More informationA BENDIXSON DULAC THEOREM FOR SOME PIECEWISE SYSTEMS
A BENDIXSON DULAC THEOREM FOR SOME PIECEWISE SYSTEMS LEONARDO P. C. DA CRUZ AND JOAN TORREGROSA Abstract. The Bendixson Dulac Theorem provides a criterion to find upper bounds for the number of limit cycles
More informationThe Higgins-Selkov oscillator
The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being
More information(1 + 2y)y = x. ( x. The right-hand side is a standard integral, so in the end we have the implicit solution. y(x) + y 2 (x) = x2 2 +C.
Midterm 1 33B-1 015 October 1 Find the exact solution of the initial value problem. Indicate the interval of existence. y = x, y( 1) = 0. 1 + y Solution. We observe that the equation is separable, and
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationTHE NUMBER OF POLYNOMIAL SOLUTIONS OF POLYNOMIAL RICCATI EQUATIONS
This is a preprint of: The number of polynomial solutions of polynomial Riccati equations, Armengol Gasull, Joan Torregrosa, Xiang Zhang, J. Differential Equations, vol. 261, 5071 5093, 2016. DOI: [10.1016/j.jde.2016.07.019]
More informationDYNAMICAL SYSTEMS
0.42 DYNAMICAL SYSTEMS Week Lecture Notes. What is a dynamical system? Probably the best way to begin this discussion is with arguably a most general and yet least helpful statement: Definition. A dynamical
More informationHilbert s 16. problem
Hilbert s 16. problem Aarhus University Aarhus 24 th March, 2017 Hilbert s Problem Does there exist (minimal) natural numbers H n, for n = 2, 3,..., such that any planar autonomous system of ODE s dx dt
More informationLectures on Dynamical Systems. Anatoly Neishtadt
Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances
More informationProblem Set 5: Solutions Math 201A: Fall 2016
Problem Set 5: s Math 21A: Fall 216 Problem 1. Define f : [1, ) [1, ) by f(x) = x + 1/x. Show that f(x) f(y) < x y for all x, y [1, ) with x y, but f has no fixed point. Why doesn t this example contradict
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationMULTIDIMENSIONAL SINGULAR λ-lemma
Electronic Journal of Differential Equations, Vol. 23(23), No. 38, pp.1 9. ISSN: 172-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) MULTIDIMENSIONAL
More informationQualitative Classification of Singular Points
QUALITATIVE THEORY OF DYNAMICAL SYSTEMS 6, 87 167 (2005) ARTICLE NO. 94 Qualitative Classification of Singular Points Qibao Jiang Department of Mathematics and Mechanics, Southeast University, Nanjing
More informationCourse 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra
Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra D. R. Wilkins Contents 3 Topics in Commutative Algebra 2 3.1 Rings and Fields......................... 2 3.2 Ideals...............................
More information2.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
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationLYAPUNOV STABILITY OF CLOSED SETS IN IMPULSIVE SEMIDYNAMICAL SYSTEMS
Electronic Journal of Differential Equations, Vol. 2010(2010, No. 78, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu LYAPUNOV STABILITY
More informationLIMIT 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.
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More information4th Preparation Sheet - Solutions
Prof. Dr. Rainer Dahlhaus Probability Theory Summer term 017 4th Preparation Sheet - Solutions Remark: Throughout the exercise sheet we use the two equivalent definitions of separability of a metric space
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
More informationAlternate Locations of Equilibrium Points and Poles in Complex Rational Differential Equations
International Mathematical Forum, Vol. 9, 2014, no. 35, 1725-1739 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.410170 Alternate Locations of Equilibrium Points and Poles in Complex
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationMethod of Lagrange Multipliers
Method of Lagrange Multipliers A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram September 2013 Lagrange multiplier method is a technique
More informationPattern generation, topology, and non-holonomic systems
Systems & Control Letters ( www.elsevier.com/locate/sysconle Pattern generation, topology, and non-holonomic systems Abdol-Reza Mansouri Division of Engineering and Applied Sciences, Harvard University,
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
More informationL. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS
Rend. Sem. Mat. Univ. Pol. Torino Vol. 57, 1999) L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Abstract. We use an abstract framework to obtain a multilevel decomposition of a variety
More informationTransverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization)
Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization) Anton Shiriaev 1,2, Leonid Freidovich 1, Sergey Gusev 3 1 Department
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationMCE693/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
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More information2. Metric Spaces. 2.1 Definitions etc.
2. Metric Spaces 2.1 Definitions etc. The procedure in Section for regarding R as a topological space may be generalized to many other sets in which there is some kind of distance (formally, sets with
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More information2) Let X be a compact space. Prove that the space C(X) of continuous real-valued functions is a complete metric space.
University of Bergen General Functional Analysis Problems with solutions 6 ) Prove that is unique in any normed space. Solution of ) Let us suppose that there are 2 zeros and 2. Then = + 2 = 2 + = 2. 2)
More informationAlgebras of singular integral operators with kernels controlled by multiple norms
Algebras of singular integral operators with kernels controlled by multiple norms Alexander Nagel Conference in Harmonic Analysis in Honor of Michael Christ This is a report on joint work with Fulvio Ricci,
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationEntrance 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
More informationTopological properties
CHAPTER 4 Topological properties 1. Connectedness Definitions and examples Basic properties Connected components Connected versus path connected, again 2. Compactness Definition and first examples Topological
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationLecture D16-2D Rigid Body Kinematics
J. Peraire 16.07 Dynamics Fall 2004 Version 1.2 Lecture D16-2D Rigid Body Kinematics In this lecture, we will start from the general relative motion concepts introduced in lectures D11 and D12, and then
More informationStokes Theorem. MATH 311, Calculus III. J. Robert Buchanan. Summer Department of Mathematics. J. Robert Buchanan Stokes Theorem
tokes Theorem MATH 311, alculus III J. Robert Buchanan Department of Mathematics ummer 2011 Background (1 of 2) Recall: Green s Theorem, M(x, y) dx + N(x, y) dy = R ( N x M ) da y where is a piecewise
More informationGeneralized Liénard Equations, Cyclicity and Hopf Takens Bifurcations
QUALITATIVE THEORY OF DYNAMICAL SYSTEMS 5, 195 222 (2004) ARTICLE NO. 80 Generalized Liénard Equations, Cyclicity and Hopf Takens Bifurcations Magdalena Caubergh Universiteit Hasselt, Campus Diepenbeek,
More informationTranscendental cases in stability problem. Hamiltonian systems
of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial
More informationD-MATH Algebra I HS18 Prof. Rahul Pandharipande. Solution 1. Arithmetic, Zorn s Lemma.
D-MATH Algebra I HS18 Prof. Rahul Pandharipande Solution 1 Arithmetic, Zorn s Lemma. 1. (a) Using the Euclidean division, determine gcd(160, 399). (b) Find m 0, n 0 Z such that gcd(160, 399) = 160m 0 +
More informationControllability of the linear 1D wave equation with inner moving for
Controllability of the linear D wave equation with inner moving forces ARNAUD MÜNCH Université Blaise Pascal - Clermont-Ferrand - France Toulouse, May 7, 4 joint work with CARLOS CASTRO (Madrid) and NICOLAE
More informationCharacterizing planar polynomial vector fields with an elementary first integral
Characterizing planar polynomial vector fields with an elementary first integral Sebastian Walcher (Joint work with Jaume Llibre and Chara Pantazi) Lleida, September 2016 The topic Ultimate goal: Understand
More informationCompetitive and Cooperative Differential Equations
Chapter 3 Competitive and Cooperative Differential Equations 0. Introduction This chapter and the next one focus on ordinary differential equations in IR n. A natural partial ordering on IR n is generated
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More information1.3 Vertex Degrees. Vertex Degree for Undirected Graphs: Let G be an undirected. Vertex Degree for Digraphs: Let D be a digraph and y V (D).
1.3. VERTEX DEGREES 11 1.3 Vertex Degrees Vertex Degree for Undirected Graphs: Let G be an undirected graph and x V (G). The degree d G (x) of x in G: the number of edges incident with x, each loop counting
More informationMATH 189, MATHEMATICAL METHODS IN CLASSICAL AND QUANTUM MECHANICS. HOMEWORK 2. DUE WEDNESDAY OCTOBER 8TH,
MATH 189, MATHEMATICAL METHODS IN CLASSICAL AND QUANTUM MECHANICS. HOMEWORK 2. DUE WEDNESDAY OCTOBER 8TH, 2014 HTTP://MATH.BERKELEY.EDU/~LIBLAND/MATH-189/HOMEWORK.HTML 1. Lagrange Multipliers (Easy) Exercise
More informationHandlebody Decomposition of a Manifold
Handlebody Decomposition of a Manifold Mahuya Datta Statistics and Mathematics Unit Indian Statistical Institute, Kolkata mahuya@isical.ac.in January 12, 2012 contents Introduction What is a handlebody
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
More information