EXISTENCE OF PERIODIC SOLUTIONS OF CERTAIN NON-LINEAR DIFFERENTIAL EQUATIONS by EUGENE BUTKOV A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF ARTS IN THE DEPARTMENT OF MATHEMATICS We accept this thesis as conforming to the standard required from candidates for the degree of MASTER OF ARTS. Members of the DEPARTMENT OF MATHEMATICS. THE UNIVERSITY OF BRITISH COLUMBIA April, 1956. ABSTRACT The theory of Poincare and Bendixson is applied to establish the existence of periodic solutions of the differential equation x + f(x, x)x + g(x) x = 0 One part of the work is concerned with those equations which can be considered as arising from small perturbations of other equations of the same type, already possessing periodic solutions. Two existence theorems are demonstrated and the stability and uniqueness of the periodic solutions is also discussed. The other part contains several, theorems stating sufficient conditions for existence of periodic solutions which cannot be treated by perturbation methods. ACKNOWLEDGEMENTS I am indebted to Dr. E. Leimanis of the Department of Mathematics at the University of British Columbia for directing and supervising my work. I would also like to express my gratitude to the National Research Council of Canada for considerable financial assistance. TABLE OF CONTENTS INTRODUCTION I CHAPTER ONE. PRELIMINARY THEOREMS ' 1 Section 1.1. The Phase Plane 1 Section 1 .2. Existence and General Behavior of the Solutions 2 Section 1 .3. Conservative Systems 4 Section 1 .4. The Theory of Poincare and Bendixson.... 5 Section 1.5. Non-conservative systems and their properties 7 CHAPTER TWO. SMALL PERTURBATIONS OF SYSTEMS WHOSE BEHAVIOUR IS KNOWN 10 Section 2 .1 . Perturbations of conservative systems ••• 10 Section 2 . 2 . Perturbations of systems with stable periodic solutions 25 CHAPTER THREE. SOME SPECIAL CASES OF SELF-EXCITED OSCILLATIONS 36 Section 3 . 1 . Preliminary remarks 36 Section 3 . 2 . Cases when f (x, v) is bounded 37 Section 3 . 3 . A case when f(x, v) is not bounded . . . . 64 BIBLIOGRAPHY : 71 INTRODUCTION The investigation carried out in the present work is concerned with the existence of periodic solutions of the differential equation where both variables t and x are assumed to be real. In the particular case when the functions f (x, v) and g(x) are constant, the above equation becomes a linear equation with constant coefficients and is readily integrable. The question of periodicity of a solution is then most readily settled by examining its analytical expression obtained by integration. A wide class of physical problems leads, however, to differential equations of type (A) which are not linear. The majority of such non-linear equations can be integrated only by approximation methods which have to be applied to each particular equation separately and which may require excessive time or labour. On the other hand the question of major physical interest is often just the eventual existence (and stability) of a periodic solution, while the actual shape of the resulting oscillation is of secondary importance. For this reason the question of the existence of periodic solutions has been investigated by a number of authors over the past thirty years. II One of the first investigators of the problem was B. van der Pol [ l l ] who considered the equation governing the oscillations in a circuit containing a vacuum triode. It was shown by Lienard [s"] that for al l values of |x equation (B) possesses a periodic solution such that al l non-periodic solutions approach, as t -»• » f the periodic one in an ttasymptotictt fashion. Similar results were later established for wider classes of non-linear equations of type (A) by N. Levinson and 0 . Smith [ 7 ] , Filippor [ 5 ] , A.V. Dragilev [ 4 ] , A. de Castro [ 2\ and others. All of these authors used in their investigations the Poincare-Bendixson theory of second order autonomons differential systems developed in the works of H. Poincare' [ l©] and I. Bendixson [ l ] . This theory forms also the basis of the present vrork. The particular techniques in demonstrating some of the theorems are, however, slightly different from those used by previous investigators. Chapter I gives a brief account of the well-known general properties of solutions of equation (A) and some basic theorems used later in the text. Chapters II and III contain several cases of sufficient conditions for the existence of non-trivial periodic solutions. Some of the results appear to be original. CHAPTER ONE. Preliminary Theorems 1.1. The Phase Plane Let us consider the differential equation This equation is equivalent to the following second order differential system inc. which, in turn, can be reduced to the first order differential equation (1 .1 .3) £ = - 4 C - . - 0 - * C « V f If a solution v = v(x) of ( 1 . 1 . 3 ) is known then the corresponding solution of (1 .1 .2) or (1 .1 .1) can be obtained by straight forward integration. Consequently, the problem is reduced to the investigation of solutions of (1.1.3 )• Let us. now consider the (x, v) - plane, known as the phase-plane. The solutions v=v(x) of ( 1 .1 .3) may be represented in the phase-plane by certain curves, called the trajectories (or characteristics) of the system (1 .1 .2) (or ( 1 . 1 . 1 ) ) . For a given point (x, v) the vector ( , 3£ )> as given by (1.1.2) and called the field-vector, has the direction of the tangent to the trajectory through this point. It is clear, therefore, that the direction of a trajectory through a given point is uniquely defined for a l l points in the phase plane except those where both and F^t vanish, which are called the singular points (or singularities). One such singularity is the origin. The others may occur on the x-axis at the points where g(x) = 0 • Points of the phase-plane which are not singular are called regular points. 1.2. Existence and General Behavior of the Solutions. It is well known (see, for instance [6] ) that the system (1.1.2) possesses a unique solution x «= x(t) satisfying given initial conditions x(t 0 ) « x Q , v(t c ) = v e provided the functions f (x, v) and g(x) are locally Lipschitzian; that is, they satisfy a Lipschitz condition in every bounded region in the phase-plane. This property of functions f(x, v) and g(x) will , consequently, be tacitly assumed throughout the present work. Trajectories of (1.1.2) in the phase plane possess some important properties which are constantly used in proofs of the theorems given in the subsequent two chapters. One such property is that two trajectories cannot intersect at a regular point. A rigorous proof of this fact can be found - 3 -elsewhere. It is also clear from (1.1.2) that the trajectories are directed curves, the positive direction being that of increase in t . In particular, a trajectory passing through a (regular) point M not on the x-axis must necessarily be crossing the vertical line through M from left to right i f M is above the x-axis and from right t© left if M is below the x-axis. At the regular points on the v-axis the slope of trajectories is vertical and these are the only points with this property. It is well known that periodic solutions of (1.1.2) correspond to closed trajectories (with no singularities on on them) and that each such closed trajectory, i f it exists, must contain at least--one singularity in its interior. Without loss of generality we may assume this singularity to be the origin because i f it were a point S = (a, 0) (a £ 0 ), the transformation x^ • x - a will reduce the equation (1.1.1) to an equation of the same type for x^ with the relevant singularity at the origin. From (1.1.2) it is clear that the origin itself (and, therefore, every other singularity) may be considered as a degenerated (closed) trajectory, corresponding to the trivial solution Since this solution is ever-present it is understood that any investigation of existence of periodic solutions actually involves non-trivial solutions only. - 4 -1.3• Conservative Systems. As mentioned in the introduction, the system (1.1.2) is generally not integrable explicitly. Apart from the linear and other particular cases, there exists an important exception to this rule, namely when f (x, v) = 0 . In this case the equation (1.1.3) assumes the form eta _ i' \ 'X (1.3.1) 5^ c = " * L ~ V 7F and is readily integrable, yielding (1.3.2) G C * ) - X - Con*t. where ^ G W = ] ^(^)-^% and X represents a parameter depending upon the init ial conditions, namely Solving (1.3.2) for* v = v(x) and substituting into the first of equations (1.1.2) we obtain the solution of ( l . l . l ) . The systems which have the property that f (x, v) = 0 are usually called conservative. It is to be noted that the question of existence of closed trajectories of (1.3.1) may be usually settled without actually solving the equation. For instance i f g(x) is positive at (0, 0) and, consequently, in a neighbourhood - 5 -R of the origin, then it is easy to see that there exist a neighbourhood of the origin such that al l trajectories passing through R^ are closed so that there is an infinity of periodic solutions. On the other hand, i f g(0)<0 and there are no other singular points i t can be shown that no periodic solutions are possible. Similar conclusions are also possible i f f(x, v) is an odd function of x, viz. at least in a neighbourhood of the origin. It was shown by E. McHarg [ ' 9 ] that if this is the case and i f g(0) > 0 then there exists a neighbourhood of the origin containing an infinity of closed trajectories, just as in the case of a conservative system. 1.4. The Theory of Poincare and Bendixson. Investigation of systems which do not fal l into the categories described in the preceding paragraph is considerably more difficult. The main tool used by practically all authors in this field is the so-called Poincare'-Bendixson theory of plane trajectories. In the particular question of existence of a periodic solution the following theorem is of utmost importance: Theorem 1.1 (Poincare-Bendixson) "If the arc of a trajectory of (1.1.2), corresponding to values of t such that t •> tQ (for some fixed tg) lies - 6 -entirely within a bounded closed region R which is free of singularities, then the trajectory is either a closed eurve itself, or it approaches spirally a closed trajectory of (1.1.2) lying within R.»» Proof of this theorem may be found, for instance, in [ l ] and will not be reproduced here. It is to be remarked that all existence theorems in this thesis are demonstrated by the application of Theorem 1.1. In addition to this we shall quote another helpful theorem. Theorem 1.2. (Dragilev) "Suppose the D.E. (1.4.1) ^ £ * 1 l » V * - 0 (where the functions f].(x, v) and g(x) are locally Lipschitzian) possesses a non-trivial periodic solution. Then the D.E. will also possess a non-trivial periodic solution provided (a) <?z (*iv) »s \oca\ly L ^ c l ^ **an (b) 4 * C o , o ) < O (c) M*>*> * $,C~,-fr) -for aH C« ,•«•). Moreover, i f and L2 are closed surves in the phase-plane representing the periodic solutions of (1.4*1) and ( 1 . 4 . 2 ) , respectively, and i f R is the interior of Li then L2 lies in E (closure of R). Proof of this theorem is given in Dragilev*s paper I k~\ as well as in a paper by Yelshin [ 1 2 ] who obtains the theorem as a special case of his "principle of comparison". 1.5. Non-conservative systems and their properties. general case when f(x, v) is an arbitrary (locally Lipschitzian) function of x and v . By putting f (x, v) » 0 we obtain an associated conservative system Although the general picture of the trajectories in the phase-plane may be quite different for ( 1 . 1 . 2 ) and ( 1 . 5 . 1 ) it is clear, for instance, that both systems have the same singular points. unless g(0) = 0 , the most important properties of these singularities also remain the same. In particular, it can be easily shown that if g(0) < 0 and if there are no singularities other .than the origin then there can exist no non-trivial periodic solution of ( 1 . 1 . 2 ) , just as is the case with ( 1 . 5 . 1 ) . Let us now consider the system (1 .1 .2) in the most (1 .5 .1) Moreover, it was already shown by Poincare that, - d -If we define the function \(x, v) by means of (1.5.2) * w h e r e GC*)- \ ^)«b<h o then X(x, v) will no longer remain constant along the trajectories of (1.1.2) but, instead of this, will possess the following property: Lemma 1.3. "If (xa , v p ) and (x t , v, ) are two points on a trajectory T of (1.1.2) then (1.5.3) X = - \ {O-.vVu-.a* where the integral on the r.h.s. is a curvilinear integral taken from (x0 , v0 ) to (x, , v, ) along the curve T Proof. Upon multiplication by vdx (1.1.3) becomes: Integrating from (x 0 , v s ) to (x, , v, ) along T we obtain *» From here (1.5.3) follows immediately by applying the definition (1.5.2). In conclusion, we shall state another lemma which will be often used in later text: - 9 -Lemma 1.4* "If g(0) > 0 and f ( 0 , 0) < 0 then no trajectory of,,(1.1.2) passing-through a regular point will approach the origin as t •* °° ." Proof. Consider the system of curves given by: ( 1 . 5 . 4 ) G C < X > ) = ^ ^ > 0 Because of the continuity of f(x, v) and g(x) there will exist a neighbourhood R1 of the.origin such that f (x, v) < 0 and g(x) •> 0 everywhere in R* . Then G(x) is non-negative in Rf and increasing as \x\ increases. From this, in turn, easily follows that ( 1 . 5 . 4 ) represents, for small values of Vj. a family of closed curves C\ encircling the origin such that i f \^ < ^2 then Cj^ lies in the interior of C^2 Consider now a trajectory of (1 .1 .2) passing through a regular point of Rf . From ( 1 . 5 . 2 ) and (1 .1 .2) it follows that along such a trajectory At 1* d t dt L AX It is seen that ^ > 0 within R1 so that the trajectories of ( 1 . 1 . 2 ) recede from the origin and the lemma is proven. - 10 -CHAPTER TWO. Small perturbations of systems whose behavior is known. 2.1. Perturbations of conservative systems. In this section we shall consider the systems which differ but very l i t t le from the conservative systems. It is convenient to represent such systems in the foim (2.1.1) fry where the parameter p. is assumed to be small. Putting u- <* 0 in equations (2.1.1) we obtain a conservative system which is usually called the generating system corresponding to (2.1.1): (2.1.2) We shall investigate the possibility of existence of a periodic solution of (2.1.1) around the singularity at the origin. For this purpose i t is natural to start with the consideration of system (2.1.2). We shall restrict our-selves to the case when g(0) •/> 0 • It was already mentioned before that if g(0) < 0 and i f the origin is the only singularity of the system then no (non-trivial) periodic solutions are possible. This will - 11 -be also true for the system (2.1.1), even i f p. is not small. Therefore we shall henceforward assume that g(0) > 0 . In this case there will exist a neighbourhood around the origin such that al l trajectories of (2.1.2) passing through it will be closed (see Fig. l ) . These trajectories will be given by GO*) - X where Q t * V \ <W) (2.1.3) and X - a variable non-negative parameter. If R^ is restricted in size the range of X will also be restricted. This can actually happen i f our system possesses other singularities beside the origin. To maintain full generality we shall restrict ourselves to an open region R around the origin defined as the interior of a closed trajectory of the type (2.1.3) corresponding to a particular value X • X1 X - V Suppose now that p. is no longer zero so that we are dealing with the system (2.1.1). If [i is small enough it is natural, however, to expect that the trajectories of (2.1.1) will differ but very l i t t le from the trajectories of (2.1.2). Fig. 1 This is indeed true and can be expressed as follows: - 12 -Lemma II.1. "Let [xi(t), vi(t)] and [x z(t), v 2(t)] be the solutions of ( 2 . 1 . 1 ) and ( 2 . 1 . 2 ) , respectively suth that where || (x, v)|| = )J x 2 + v 2 and both points [se.CO,"*,(.+,>)] and [/XxCto), ^ i C t © ^ ] belong to an open neighbourhood R of the origin. Then, as long as the points [xi(t) , v^t)] and [x2(t), V2(t)] remain within R the following inequality will hold for t > te-ll I>, Ct> C-0] -[**(t\* a(-fcVl|| £ where M • «u<p \ 40* ,>»•)• n*| and L is a constant Ox ,v) e R. Proof. The systems ( 2 . 1 . 2 ) and ( 2 . 1 . 1 ) can be respectively written in the vector notation as follows: -( 2 . 1 . 5 ) a ? = m ^ and ( 2 . 1 . 6 ) + where ^ « Z • O t t ^ G t ) ] * -^teyoc] Let <p = <p(t) and \Jr = ty(t) be the solution-vectors satisfying ( 2 . 1 . 5 ) and ( 2 . 1 . 6 ) respectively, such that - 13 -(fC±o) = I/ar.Ofe.), tr, t t 0 ) ] Y ^ o ) - [ * i G O / * t ( * 0 ] Then from ( 2 . 1 . 5 ) and ( 2 . 1 . 6 ) follows that -t q>(t) = qX-fco) + J R [ c f C « i ) ] - d s -to = ' Y C-to) + j { R fri*)] + & [Y(s)] J As so that ( 2 . 1 . 7 ) |yW-*?WH^ W^-W^-MI * j l qcY)- + jll*MII-<k "to t „ Let Since g(x) is Lipschitzian and f(x^ v)v is bounded by M we can write Consequently ( 2 . 1 . 7 ) becomes •t (2.1.d) 0(-t) s 6 + L • j ©(0 -^ s + M • M • Ct-^ to) to If we define ^ © t t ) = j e u V < U "to then the above inequality becomes © ' ( • 0 - L • © f t ) ^ e + M - M - C-t-to) - 14 -Multiplying by exp [ - L(t - to)] and integrating from tQ to t we obtain ® f t ) • [-Ui-to)] < £ [ \ - e^op[-Lft--to)]^ + so that © f t ) * i {e^[Ut-t.)] - |1[ + ^{e^p[Ut-to)]-l--MH(t-"U)]. Substituting this expression into (2.1.6) we get 9(t) < £ • -e^ > [LCt-to)] + ^ ^ [LCi-to)] - I ] and the lemma is proven. Let us now choose a point A(XQ, VQ) in R and consider a trajectory of (2.1.1) passing through i t . If JJ- = 0 this trajectory is closed and given by (2.1.9) ^ s v « = |/ 2 ^o- 2 G W W l l e r e A I Z X, / \ In view of Lemma II.1 it is obvious that for \i •£ 0 but sufficiently small our. trajectory will have, in general, a "spiral shape". By this we mean that it will encircle the origin and cross the radiusvector AO through A at some point C lying close to A . We are, of course, interested whether C will be closer to the origin than A or not. To answer this question we shall consider the following function: (2.1.10) WOO = ^ 4 [ * , * \ C » A ) ] - ^ & c , X M * - 15 -where the integral should be taken in clockwise direction around the closed curve defined by ( 2 . 1 . 9 ) . Function W(\) can be evaluated for al l values of \ such that 0 < \ < XQ . Its importance is seen from the following theorem. Theorem II.2. Let A = (x„ , v 0 ) 6 R and let W(\0) *> 0 where \ 0 = 1/2 v Q 2 + G(XQ) . Then a number u.Q •> 0 can be found such that for a l l values of \i satisfying 0 < JJ. < \l0 the trajectory T of (2.1.1) passing through A will make its successive crossing of radiusvector AO at a point C «* (x j^V )^ such that Proof. Let ^(t) and cp(t) be the vector-functions representing the solutions of (2.1.1) and ( 2 . 1 . 2 ) , respectively, which pass through point A at t = t 0 , viz. We know that <p(t) is periodic. Let its period be T . Then from Lemma II.1 follows that for t e E^o^o+T] we shall have (2.1.11) |K0O-<fCt)| - ^ { ^ p L T - l ] = I provided ty(t) remains within R . This inequality means that the point B defined by t ( t 0 + T) will be within a - 16 -c i r c l e of radius S centered at the point A (see F i g . 2 ) . By choosing \x s u f f i c i e n t l y small we can make 8 a r i b t r a r i l y sma l l . This w i l l also ensure that y ( t ) remains w i t h i n R for t € [ t „ , t 0 + T l so that (2.1.11) w i l l be ac tua l ly app l icab le . Let us continue ( i f necessary) the t ra jectory defined by ^ ( t ) beyond the point B u n t i l we in tersec t the rad ius -vector AO and c a l l the point of in t e r sec t ion C . We proceed to show that X(C) < \ (A) From Lemma 1.3 we know that (2.1.12) Q = - = . . Transforming t h i s l i n e - i n t e g r a l defined br cf (t) i n t 0 a n ordinary in t eg ra l <iefmed by w(t) with respect to t we get F i g . 2. to S i m i l a r l y , from (2.1.10) we obtain so that to - 17 -t«,*T to t. Since f(x, v) is Lipschitzian in R we have Therefore to "to Using the inequality (2.1.11) and also recalling that \ f (x, v)v | < M and we obtain ^ ^ T "t* 'to Finally, .vi is bounded in R-j_ so |vjj < b and we obtain where It follows that if n is chosen sufficiently small Q will have the same sign as u-W and, therefore, be positive, which means that \(B) < \(A) • - IB -To complete the proof we also note that c |X(c ) - X ( f t ) | = 1 j l * ^ C t y A a e l < 8 Now we have It follows that for sufficiently small n the quantity X(C) - X(A) will have the same sign as u-W . This yields X(C) < X(A) and the theorem is proven. From the proof of the theorem the following corollary is immediate. Corollary. Let A = (x0, v 0) 6 R , W(XQ) and C(xi, Vj.) be defined as in the theorem and let W(X0) < 0 • Then for sufficiently small n the inequality will hold. Using the results just developed we are now in a position to make some statements about the existence of closed trajectories of (2.1.1) for sufficiently small values of u- . It is quite clear that if such trajectories exist then they should be sought in the vicinity of such trajectories - 19 -Tx of (2.1.2) for which W(X) = 0 . Let us p lo t the funct ion W(X) v s . X (see F i g . 3 ) . Suppose W(X) has a zero for some pa r t i cu l a r value of X , say X = X* ^ 0 . Suppose t h i s zero i s i s o l a t e d and such that i n a neighbourhood of X* W ( X ) > 0 -for 1>Y W W < 0 [or \<\* Consider a t h i n r i n g around the o r i g i n bounded by the t r a j ec to r i e s f\*-£ a n c ^ H\*+£ of (2.1.2) corresponding to the values of X equal to I f e i s s u f f i c i e n t l y small we s h a l l have W(X* - e) < 0 and ¥(X*+ e) > 0 Let A]_ and A2 be the points of i n t e r sec t ion with the pos i t i ve v-semiaxis of the t r a j ec to r i e s (2.1.2) and X*-£ of F i g . 5 - 20 -By Theorem II.2 it follows that, for sufficiently small values of |J. , the trajectories T^ and T 2 of (2.1.1) passing through and A2 , respectively, will cross again the positive v-semiaxis at points and C 2 , respectively, such that lies below A-^ and C 2 lies above A 2 (see Fig. 5)» This establishes a region K (shaded on Fig. 5) bounded by arcs A^C^ and iC^ Cj? o f T l a n d T 2 » respectively, and by segments A-^ C^ and A 2 C 2 . It is clear that a trajectory of (2.1.2) passing at t » t 0 through a point inside K will remain inside K as t •-*• 00 and so it will satisfy the conditions of Theorem 1.1. It follows that for sufficiently small values of n the system (2.1.1) will possess a non-trivial periodic solution. The closed trajectory corresponding to this periodic solution will approach, as |i -> 0 , the closed trajectory Tx* of (2.1.2) corresponding to X = X* . Suppose now we are dealing with a zero X = X* ^ 0 of ¥(X) such that in a neighbourhood of X* W(V) < O W X>X* a n d W U ) > - 0 \or X<X* In this case we can construct, in the same fashion as before, a region K which is trajectory of (2.1.1) can enter. If we now make the change of independent variable - 21 -in the system (2.1.1) we can easily see that the new system, viz. (2.1.12) will have the same trajectories as (2.1.1) but their direction will change. Since the region K will satisfy the conditions of Theorem 1.1 for the system (2.1.12) it follows that the system (2.1.12) and, therefore, system (2.1.1) as well, will possess a non-trivial periodic solution i f \i is sufficiently small. We can now summarize our results in the following Theorem: Theorem II.3 For sufficiently small values of u. the system (2.1.1) will possess non-trivial periodic solutions reducing, as |i -*• 0 , to such closed trajectories Tx* (X* ^ 0) for which W(X*) • 0 and for which the quantity W(X) . (X - X*) does not change sign in some open neighbourhood of X* • The investigation conducted thus far has not considered the question of stability nor the related question of the number of periodic solutions arising in the vicinity of the curve ^ X* . These questions can be, at least partially, answered in the case when the function W(X) is c ontinuously different iabl e. 22 -In order to consider this problem we introduce the following definition of stability, sometimes also called the strong, or asymptotic stability. Definition. Aperiodic solution <p(t) of (1.1.2) is called stable i f e •> 0 can be found such that every solution t(t) of (1.1.2) satisfying I M - t o ) - q>cto)|| ^ £ will possess the following property: For every point cp(t*) of the periodic solution <p(t) a sequence ( t n | (jn = 0,1,2, ...) where t n -*• 0 0 as n -»•« can be found such that (2.1.13) i i m IKC-U) - if (t*)| = O It was already proved by Poincare that a non-trivial periodic solution of (2.1.1) will be stable i f the following condition is satisfied. o where x *= x(t) and v = v(t) represent the solution under consideration and T is its period. We shall now show that the function is nothing else than the integral 3 , as defined above, taken along a curve Indeed, if [x]_(X)f o] and [x2(X), o] are - 23 -the points of intersection of T*\ with the negative and positive x-seraiaxes, respectively, and if -v^O.V) ^ + ]j 2 X - Z G M ^ " C-xA) = - ^ 1X-1 G(*) then ( 2 . 1 . 1 0 ) can be written as: Tt(X) - i W1. 4 [ o c , ^ w ( - x , X ) ] • v ' O x ^ V eke Differentiating with respect to X and remembering that v+(x, X) and v~(x, X) vanish at x-j^ X) and x2(X) we obtain dX J I ^ c\\ ^X ) _ I f J ? - . * - , * f I a * J l d*' dX J dX ) But we have "Jj^ = v+ and "jx s • Substituting this in the above formula and changing back to the curvilinear integral we get Finally, i f we take t as the independent variable and put dx = v-dt we obtain At = 3 - 24 -T aw ax o where T is the period of our solution. Consider now again a closed trajectory f\* of ( 2 , 1 . 2 ) such that We have already seen that for sufficiently small values of u- there will exist a closed trajectory T of ( 2 . 1 . 1 ) within an arbitrarily thin ring around Because of the continuity of ^ it can be shown (we shall omit the detailed proof) that the value of 3 computed along T will have the same sign as "jj^ (X*), provided u. is sufficiently small. From this we may conclude that the resulting periodic solution of ( 2 . 1 . 1 ) will be stable. Since it is a well known theorem (see, for instance [ 7 ] ) that two adjacent closed trajectories cannot be both stable i t also follows that the curve T\* will generate a unique stable periodic solution of ( 2 . 1 . 1 ) . Remark. If cp(t) is a solution of an autonomous system then cp(t + t f c), where t* is arbitrary, is also a solution. In order to avoid reference to this fact every time we consider the uniqueness question, we shall henceforward assume that by saying that a solution cp(t) is unique we shall mean that there are no other solutions except those of the form cp(t + t*) - 25 -It is not hard to see that those zeros of W(\) for which -j^ < 0 will be similarly associated with unstable periodic solutions of (2.1.1) where the definition of an unstable periodic solution is the same as that of the stable one except that the sequence l^*^ is replaced by another one, {t n ] such that tQ-> - °o as n •+ » . It is to be remarked that the above theory is insufficient to treat the case when both W(\) and ^ vanish simultaneously for some value of X . 2 . 2 . Perturbations of systems with stable periodic solutions. In this section we shall discuss the case of small perturbations of non-conservative systems of type (1.1.). Suppose that the unperturbed (generating) system possesses a periodic solution corresponding to a closed trajectory _TL in the phase-plane. If we perturb this system by introducing an additional small term in the differential equation (1.1.1) then the question arises whether the perturbed system will also possess one or many periodic solutions, their respective trajectories reducing to J7. as the relative magnitude of the perturbation decreases to zero. This question may be partially answered by the following theorem: Theorem II.4. Consider the differential equation (2.2.1) ^ + * l « , , ) . » * < f c C « V ~ - 0 - 26 -where f(x, v) and g(x) are continuous, bounded and Lipschitzian functions in an open region R around the origin. Let there be no other singularities in R except the origin. If the equation (2.2.1) possesses a stable non-trivial periodic solution whose trajectory si lies in R then the differential equation ( 2 . 2 . 2 ) 3^ + [4C*i^^M^C*.*)]-^+ ^ V ' X . O M * 0 will also possess a non-trivial periodic solution provided h(x, v) is continuous, bounded and Lipschitzian in R and u. is sufficiently small. Proof. Let X(t) = I x © ^ * v o ( t ) J < w h e r e v o ( t ) * "35 ) be the stable periodic solution of (2.2.1) under consideration and SI - the closed trajectory associated with i t . Choose a point and consider the point A = (a, b) lying in the interior of SI. on the radiusvector through X(tQ) such that where e is the positive constant as defined for X(t 0) by the condition of stability (2.1.13). - 2? -*• ~x. Let q>(t) = [x 1 (t) , v^t)] and t(t) = [x 2(t), v 2(t)] (where ST± Ot) • ^ ^ 1-1,20 Fig. 6 be the solutions of (2.2.1) and (2.2.2), respectively, such that Let us follow the trajectory T, described in the phase-plane by cp(t) • Let A1 be the point where crosses the radiusvector through A after one complete revolution (see Fig. 6). It is obvious that A1 lies between X(t0) and A ; otherwise the stability condition (2.1.13) could not hold. As usual, we iset and note that \(fi')-Mr>) - p > O Consider now the trajectory ^ defined by ty(t). We shall show that after a complete revolution around the origin this trajectory shall cross the radiusvector through A at a point A" such that X(A") > X(A) • Let the constant T be defined by means of - 23 -Now, in exactly the same way as in Lemma II . l we can J deduce that for all t such that t 0 < t < t Q + T we shall have (2.2.3) . ||v«-q>C-b)U £ ^ r ^ p L T - l l ^ L L where M = sup h(x, v) and L is a constant. (x ,v )«R The inequality (2.2.3) means that the point Af" • ty(t0+T) will lie within a circle of radius S around A* -. In order to estimate \(An') - MA) we use Lemma 1.3-Putting dx = v»dt in (1.5.3) we transform the curvilinear integral into an ordinary one and obtain for our case: XCR") - M a ) » - ^ ( a ^ ) + M - ^ ^ - ^ - < t t t 0 Similarly, we can write to We shall also need the inequality *Cft"0 |X («'•) - XCR'01 = I j [ f O r ^ , y M^C*z^]v»Aoc( < ( N ^ H H - S x CR") • where W = sup M Now we can write | [ X ( f ^ - X U O ] - [ M f l V M n ) ] | ' -s l>CR ,0-ACfi"')h|C^^ ," >)->(^]-[^(B')-X^)]| * •to+T "to "to where K is the Lipschitz constant for f (x, v) in R . Using (2.2.3) once again we finally obtain l [ x c v ) - X ( n ) ] - [ x t B O - X { f t ) ] | £ where D is a quantity which remains bounded as u- •* 0 . It follows that for u- sufficiently small the quantity X(A") - X(A) will have the same sign as X(A') - \(A) and, therefore, be positive. It follows that Att is farther from the origin than A and, therefore, must lie between X(t0) and A , as stated. - 30 -Let us now choose a point C * (c, d) belonging to the exterior of Q. and on the radiusvector through Z(tQ) (see Fig. 7), such that •||xttoV(e,aMI An entirely analogous calculation shows that the trajectory 1"^ of (2.2.2) passing through C will , after a ful l revolution around the origin, cut the radiusvector through X(t 0) at a point C» closer to the origin than C . Fig. 7 It may be noted that in this case one has to ensure that the relevant trajectories remain within R . It is easy to see, however, that by choosing first 62 and then u. sufficiently small this can always be achieved. It is now clear that we have established a ring-shaped domain R1 (shaded area on Fib. 7) which the trajectories of (2.2.2) can enter (via the segments UU" and AAn) but cannot leave. The existence of a non-trivial periodic solution follows immediately from Theorem 1.1. 31 -The question which arises next is whether a periodic solution which exists under the conditions of Theorem II .4 will be unique and whether it will be stable. This question can be answered using the additional assumption that the stability integral T ( 2 . 2 . 4 ) 2 - J l * - 2 ^ +4<« .*V| -At . o has a positive value for the original (unperturbed) stable periodic solution of ( 2 . 2 . 1 ) (here T is the period of this solution). It was already mentioned that 3 "> 0 represents a sufficient condition for stability. It is not a necessary condition. One can actually construct differential equations having stable periodic solutions with 3 = 0 . However, such cases are rather artificial and in the equations arising from physical problems the condition 3 > 0 will hold for stable periodic solutions in practically al l cases. It may also be pointed out that the condition 3 > O may be often established without actually solving the equation and evaluating the integral (see, for instance, \_ 7 ] , the section on uniqueness theorems). With these remarks we shall now proceed to prove the following theorem. Theorem II.5 If —^— and ^ — exist and are bounded in R and i f the periodic solution X(t) of ( 2 . 2 . 1 ) is such - 32 -that the stability integral 3 associated with it is positive then the periodic solution of ( 2 . 2 . 2 ) described in Theorem II .4 is unique and stable. Proof. Let SL be the closed trajectory corresponding to X(t) • [ x 0(t), v 0(t)] and let T be its period. Similarly, let J l , and T^ = T^ (fx) . be the trajectory and the period for the periodic solution of ( 2 . 2 . 2 ) under consideration. From the proof of Theorem II .4 it is clear that given e1 > 0 and sufficiently small, we can find a value of [i such that will possess at least one point Q lying within the circle of radius e f centered at the point P - .^(.ioi.-tfot-to)"] , with -io(X*) + 0 , on SI Let X,C-0 - [•K.C'Oj'tf.Ctyi be the periodic solution of ( 2 . 2 . 2 ) satisfying and having, therefore, _QLi for its trajectory and T-^ for its period. At the time t = t 0 + T the vector function X(t) will return to the point P while X^(t) will be at some distance from Q where 0 as u- 0 . This readily follows from the considerations used in the existence proof. Since v 0(tQ) ^ 0 we can choose sufficiently small so that V]_(t0) £ 0 and then decrease u- again, i f necessary, to f u l f i l l p± < 1/2 • \ v ^ t ^ | . Let T x = T+t(uJ; then from - 33 -follows oeX-to^ -T.) where the integral is a curvilinear integral taken along _S1, in the direction of increasing t . ' Throughout the integration we have \<v\ -> \<\5,|-p, •> •> O so that \ T ( » 0 | < ; ; * It follows, that t(n) •* 0 as n •+ 0 .. In other words, T(uJ is continuous in u. in the neighbourhood of u. = 0 • Having established this, we are now in a position to show that the solution %(t) is stable. Let 3 and J, be the stability integrals for X(t) and X1(t) , respectively, viz. T o where fi(x, v) = f (x, v) + \i • h(x, v) . We shall show that -*• 0 as u. -»• 0 . Indeed, - 34 -T, T T T Using the fact that ^ , J* a n d h^x> v^ a r e bounded and that f(x, v) is Lipschitzian, we have < k M I x ^ - x w u + * ^ where and K>? are some constants independent of u. • Because of the continuity of = T (^u.) we can also write _ T, T, 1 ov T -At < where K3 is independent of |J, It is now clear that |o, -3 | < K/,, u- for some constant .K/,, and, therefore, - 35 -Since 3 •> 0 we conclude that for sufficiently-small u- we have 3, > 0 • We have therefore shown that a l l periodic solution of (2.2,2) generated by X(t) must lie" within a thin ring (whose width approaches zero as \x •* 0) around the trajectory SL of X(t) and that al l these solutions are stable. But it is a well-known theorem (see, for instance, [7*] ) that two closed stable trajectories having a common singularity and no other ones in their interiors cannot be adjacent, that is , the ring-shaped region bounded by these trajectories must contain at least one closed trajectory which is not stable. Since this is not possible in our case we conclude that the trajectory Sl{ is unique, and so is X,(t) . In conclusion, we shall remark that the theorems II.4 and II.5 will remain valid i f we replace "stable" by "unstable" and replace 3 > 0 by 3 < 0 . The proofs are completely analogous to those above except that a l l trajectoreis are followed in the direction of decreasing t . As an alternative, one may also introduce the new independent variable tT = - t and apply Theorems II .4 and II.5 directly to the new differential system. CHAPTER THREE: Some special cases of self-excited oscillations. 3.1. Preliminary remarks. In this chapter we shall state and prove several theorems establishing the existence of periodic solutions of differential equation (1.1.1) for some special cases, which cannot be treated as small perturbations of already known equations. The proofs of al l of these theorems are based on the same general principle. We proceed to show that al l trajectories which start at points sufficiently far from the origin have the general shape of spirals convolving inwards. At the same time the conditions f (0, 0) < 0 and g(0) > 0 permit us to use Lemma 1.4 to show that trajectories starting sufficiently close to the origin must recede from i t . This permits us to establish a closed region free of singularities which satisfied the conditions of Poincare-Bendixson theorem and ensures the existence of a periodic solution. At the same time we also establish that al l solutions are bounded. This is also an important property in applications. It should be noted that in a l l cases the conditions * imposed upon functions g(x) and f(x, v) are sufficient, but by no means necessary for the existence of periodic solutions. It is because of the rather restricted nature of the problem, as well as because of great variety of possible - 37 -behaviour patterns of non-linear systems that these conditions may appear to be fairly numerous or complicated. 3 * 2 . Cases when f(x, v) is bounded Theorem III.l Consider the differential equation ( 3 . 2 . 1 ) + {c~/0-* + «frC*y* = o v - ^ -Let the functions f(x, v) and g(x) be locally Lipschitzian and satisfy the following conditions (1) c^i-x) > O £or all « (2) _ o » bc, \-c\_>oo (3) 4Co,o^<o (4) <$C-x,t) -a CHvO) -for >-l0(W*o),all > N > M ^or ^ S - , all . Then the equation ( 3 . 2 . 1 ) possesses a non-trivial periodic solution. Proof. Define (5) q>0<0 It is clear that satisfies M -for «ff<-Zt cp(v) is also locally Lipschitzian and - 3 8 - -for a l l (x, v) . We shall prove the existence of a non-trivial periodic solution for the equation (3.2.2) ^ + cfCvV* +e % C*V*-0 tf= ~ Then the existence of a non-trivial periodic solution of (3.2^1) will follow from Theorem 1.2. Let L(a) » inf g(x) . Let us choose a •> 0 large enough to satisfy (6) LCaV a - 2.HI, > o (7) a ^ - H ) - ^ u a V a . 2 M W I -a and (8) r-t.*) = J ^ - x - A * > 2b4 G(-a) This is always possible in view of (l) and (2). Let us consider the point A(-a, -2b) and the trajectory T of (3.2.2) through i t . If we follow T in the direction of increasing t we can see that it must cross the x-axis at some point C such that x(c) < - a . To show this draw the straight line AB (see Fig. 8) with the slope The line AB intersects the x-axis at B with the abscissa CS) = - a -LCaVa - 2Mb - 39 -of T for x < - a , On the other hand the slope ^ and v < 0 s a t i s f i e s as fol lows from (1), (2) and (5) . From t h i s r e s u l t i t i s c lea r that T runs above AB and, therefore, cuts the x -ax i s at a point G(-a*, 0) with c^(fc) - a* < - a Let us now consider the behaviour of T i n the second quadrant The slope of T i s given now by (3.2.3) CV3C - 40 -This slope is greater than the slope of the curve P defined by (3.2.4) ^ = 0 = ^ Z G W - 2 G C - a*) 3C where * J e ^ V ^ o The curve T passes through point C , runs through the second quadrant with v (^x) increasing monotonieally and crosses the v-axis at [ 0, vi(0)^ . From -a* < -a follows that G(-a*) > G(-a) and then, by (8) and (3.1.4), that vx(0) > 2b . Since the slope of T is greater than the slope of P but bounded from above for -a* < x < 0 (except in the neighbourhood of the x-axis), it follows that T lies above P and must intersect the v-axis at a point E with v(E) > 2b . In the first quadrant the slope of T is bounded from above by M , therefore T must remain below the straight line ~ Hoc L x for some b 0 > 0 • This leaves ^ bounded away from zero and from (2) it follows that the slope of T will become and remain negative for sufficiently large x . Hence T will descend and cross the x-axis at a point denoted by G(c*, 0) (see Fig. g). Let F(c, 2b) be the point in the first quadrant where T intersects the horizontal line v = 2b for the first time. Without loss of generality we may assume that c satisfies the inequality (3.2.5) * C.N-M) > Vo Indeed, if this were not true we could find C]_ > 0 satisfying (3.2.5) and consider a trajectory T-^ through I?l(cl> 2b) and trace it back through the f irst , second and third quadrant, thus defining points D ,^ C-^ and A]_ corresponding to the points D, C and A on trajectory T. Clearly A^ would satisfy the conditions (6), (7) and (8). In this case we could apply al l our considerations to Ti instead of T . Let us now reflect the arc ACDEFG of T in the x-axis to obtain the "image" T« of T . Let A T , D», E» and F f be. the images of points A, D, E and F , respectively. Let us consider the point H(c*, -2b) and the trajectory S through i t . Let us follow S in the direction of increasing t , that is to the left of H . We shall proceed to show that S must intersect the v-axis at a point 3 which lies above E* . Let v, v1 , | ^ and ^ be the. ordinates and the slopes of S and T1 , respectively, within £o, c*]. 42 -We have (3.2.7) where the last equation follows easily from the fact that T 1 is an image of T . Therefore, Since both v and v* are negative, it follows that whenever v < v1 we have If, in addition, v < - 2b this becomes (3.2.9) ~ - ~ > - H - v ^ *>0 dot owe In particular, it follows that if S intersects TT within [o, c] , it must do so in the direction indicated by arrows on Fig. 3. This fact will be used later. Let us now return to the point H and follow the trajectory S to the left from i t . If S intersects the line x = c above the point F 1 then, in view of the preceding remark, i t will remain above T1 throughout the interval [o, c] and, consequently, will intersect the v-axis at a point 3 lying above E» . We shall therefore assume that v < v 1 at x • c . If S intersects T within [o, c] then our statement is proven, so we shall also make the assumption that v < v' < - 2b throughout the interval [o, c"| . In this case let x* be the least abscissa in [c, c*] for which v(x*) = - 2b . Then, since vf(x*) < 0 (3.2.10) -NS'C**^ - <* - vi**) - Z\> Now, between x = 0 and x = x"* we have v(x) < * 2b , as well as v < v* so that the inequality (3.2.9) is applicable. Integrative (3.2.9) between x • 0 and x = x* we obtain x* o hence _ > (M - M V c - > O as follows from (3.2.5) • This result contradicts the hypothesis that v < vT throughout the interval [o, c] . Consequently S must cross T f and the point 3 must lie above E' , as shown on Fig. g. Let us now return to the point . A and trace the trajectory T in the direction of decreasing t , that is, to the right. By arguments similar to those employed in tracing T through the second quadrant we can easily see that T must l ie below the line v = - 2b and since its slope is bounded from below, it must intersect the v-axis at some point K . We shall show that K must l ie below E* . This result follows immediately if the point DT lies above A since it is easy to check that the inequality (3.2.9) will also hold in - 44 -[-a, o] and this will ensure that T will remain below TT within this interval. We need, therefore, to consider only the case when D* lies below A • First, we shall derive an inequality for the distance DA' which is equal to the distance DfA . For this purpose consider the arcs CD and CA1 of T and TT respectively. Let v, v 1 , ^ and ^ be the ordinates and slopes of T and TT , rexpectively, within the interval [-a*, -a] . We can easily derive that (3.2.11) - - c f ^ ) - c ? C - ^ + ^ * V * - ( i . - i f ) Let x** be the greatest abscissa in the interval [-a*, -a] for which v = v' • Then we have v v f. inside the interval ^x** , -aj and, since both v and v f are positive and x < 0 , (3.2.11) becomes (3.2.12) ^ - ^ ' £ - t f M - C f C - v ' ) < M + M =ZM Therefore _ a But e * * * - - a* * -=c(s) = - a — -Consequently — 4t a PCV * 2 . M - — — ^ — — LCaVa - 1MW Since D1 A = US' we can write - 45 (3 .2 . 13 ) jyi\ 8 M l o Z L(a) - a - 2.H W We can now return to our considerations about the arc AK of T in the third quadrant. Let us denote again by v, vT , 3 ^ a n d 3<e t n e ordinates and slopes of T and T f , respectively, and assume that T remains above T 1 throughout the interval [-a, . Then-.the inequality ( 3 . 2 . ) will be again applicable, integrating ( 3 . 2 . 9 ) between x - - a and x • 0 we obtain But, by (3 .2 . 13 ) 8 M b 2 L(a.V a -so that as follows from ( 7 ) . This contradicts our hypothesis and shows that K must l ie below E* • Let us now consider the region RT bounded by the arc KCEG of T , by the vertical segment GH , by the arc H3 of S and by the vertical segment "3K . As outlined in the proof of Lemma 1 .4 we can construct, using the fact that 9 ( 0 ) < 0 , a neighbourhood of the origin which no trajectory can enter. Excluding this neighbourhood from R* we obtain a ring-shaped region R free of singularities which satisfies the oondition of Theorem 1.1 and the proof is complete. An additional remark can be made with regard to condition (2). This condition was used essentially only once, in order to shew that the trajectory T running through the first quadrant must finally attain a negative slope and cross the x-axis. If we replace (2) by a weaker condition (2») g(x) > k for \x \ sufficiently large then it can be shown that the argument will s t i l l hold provided that ,— where M is defined as in (4).. Theorem III.2 Consider the differential equation (3.2.14) ^ +[^C«,>r) + M .X(*) ]">T + ^('xV'DB-0 Let the functions h X and g be locally Lipschitzian and satisfy the following conditions: (1*) >0 {or &\\ -X. (2*) c^-oc) = Cfc(~) (|*) c^(-x) o» as 'ac-*oo (4*) ^ (0,0s) * O (5*) ^(oc^) > W >0 W \v\ > W , all -x j ^ C c^,-*) ^ - M (*4>o) jor W\<lo , ill so . - 47 -(6*) X(-*0=-X0=O (7*) | X W | < V. W all -x. Then the equation ( 3 . 2 . 1 4 ) possesses at least one non-trivial periodic solution for al l values of the parameter u. . Proof. Let us define N - M + CH + KJ) — K) - ( M + VO) for \ M > Zfe 4or - Z W f i ^ - t The function f(v) is locally Lipschitzian and satisfies the inequality for a l l (x, v) It follows from Theorem 1.2 that it will then suffice to prove the existence of a non-trivial periodic solution of the equat ion ( 3 . 2 . 1 5 ) ^ •+ 4 M.X(ac)].^+ = ° Let us choose a > 0 sufficiently large to satisfy (9*) (10*) L-a - ( M + M ^ > O •as *a a > - 48 -Let us consider the point P = (-a, -2b - p. K a) in the third quadrant of the phase-plane and the trajectory T* of (3.2.15) through it (see Fig. 9 ) . If we follow T* in the direction of decreasing t , c _r ^ r-. CO P, Zk 1 P • f t that is to the right, we can see, using (1*), ( 5 * ) , (S*) and (7 V ) that its slope, is bounded from above by p.K and from below by Fig. 9 where L» - sup g(x) . It follows that T* must l ie between x<a the straight lines PPi and PP2 with the slopes u-K and (-N-u-K-L1 • ) respectively and, consequently, must intersect the v-axis at a point denoted by A • It is also clear that v < - 2b for all points on the arc PA of T* • If we now follow 7 * from P in the direction of increasing t we see that ^~ is bounded from above by (M + u-K) and, therefore, T* must l ie above the straight line through P with that slope. This leaves the ratio V (which is non-negative in the third quadrant) bounded away from zero and then from (3*) it follows that, for \x\ sufficiently large, - 49 -w i l l become and remain negative. Hence T* must cross the x -ax i s at a point denoted by C = (-a*, 0) . We also notice that i f B i s the point i n t e r sec t ion of T* with the l i n e v = - 2b (the one wi th the largest abscissa i f there are more than one such points) then we have Let us now consider the point Q = (a, 2b + (j,Ka) i n the f i r s t quadrant and the t r a j ec to ry T** through i t (see F i g . 10) . By employing arguments analogous to those above we a r r ive at the conclusion that T** must cross the v-axis at a point D and that for a l l points on the arc DQ of T** we s h a l l have v i 2b , A l so , T * * must cross the x -ax i s at some point F = (a** , 0) and the hor izonta l l i n e v = 2b at some point E such that a. * * > rc ( E ) •> a v v F i g . 10 F i g . 11 - 50 -This done, we define c = max- (a*, a** ) and consider the point - (-c, 0) and the trajectory T^ through it (see Fig. 11). Let us follow T-L in the direction of decreasing "t through the third quadrant. It is quite obvious, in view of the preceding analysis, that Ti must cross the line v - - 2b at some point Bi = (-a^, -2b) with a, \sc(e-M>a and also must also must cross the v-axis at some point A^ . More-over, a l l points on the arc B , A| of Ti must satisfy v < - 2b . Indeed, i f c = a* this is so because T^ coincides with T* • If, however, c > a* phen B^ must l ie to the left of B and the slopes of all trajectories crossing the segment B^ B satisfy (on B^ B) the inequality (where (9*) was used again). From this immediately follows that Tj_ must remain below the line v = -2b between x = -a , and x = x(B) • Befdre we go further we shall find an upper bound for the quantity (c - a±) which will be needed later. Let us draw through B-JL a straight line B^D^ with the slope This line crosses the x-axis at a point Pi with the abscissa Co,") = - a , -n 51 -L-a - 2W CM + IUL!C) Since fin is an upper bound for the slope of arc 1^^ 1 °^ 1^ » ^is a r c a^s n o P ° i n t s below the segment DjB ,^ so can not lie to the left of • This yields (3.2.16) c - a , * W C o ^ - v a . l « — r Let us now reflect the arc . ciBx% o f ''"l i n t l i e v-axis to obtain the arc CZB» Ai of the "image-trajectory" T,1 , points C 2 « (c, 0) and B,T = (a,, -2b) being the images of and B^ , respectively. Let us- consider the trajectory T2 passing through C2 and follow it in the direction of increasing t , viz. through the fourth quadrant to the left. After T2 leaves the immediate neighbourhood of C 2 its slope becomes bounded from below and it follows that T2 must intersect the v-axis at some point F2 • We shall show that Fg lies above A2 To prove our assertion let v^ v* , 3 ^ and be the ordinates and the slopes of T2 and T.* , respectively. Then and = 40*') - M - X W 5V so that (3.2.17) Let E 2 be the point where T 2 intersects the vertical line x - ai . If E 2 is below B£ then we shall show that the distance B£ is bounded and estimate its bound. Assuming that E 2 is below B^ let x 2 be the least value of x in the interval [ a^, c] for which v(x) = v f(x) Between x - a± and x • x 2 we have v < v1 so that the third term on the r.h.s. of (3.2.17) is positive and we have Integrating between x = a^ and x = x 2 we get Remembering that v f(x 2) = v(x2) and x 2 < C and using (3.2.16) we obtain 8M102" V(a,) -^(a , ) < Z H Cc-d,Y £ Let us now consider the interval [ 0, a^ ] and make the hypothesis that T9 remains below TT throughout the * 1 interval. Then we have v < vT < 2b and (3.2.17) becomes Integrating between x = 0 and x = 3, , we obtain - 53 -or ; T ~ 2.Ma L-a - ZW Cvu nVO But the r.h.s. of this inequality is negative because of ( 9 * ) , so that v '(0) < v(0) . This, however, contradicts our hypothesis. It follows that T 2 must cross T£ in the interval [o, ai] • But then T 2 must remain (as t increases) above T£ because a l l trajectories crossing the arc A-jE* of T£ must do so in the direction shown by the arrows on Fig. 11; that is , they can only enter the region bounded by the arc A3BJ of T| and the lines x = 0 and v = -2b . This follows from (3.2.17) by putting v = v1 and remembering that f (v1) * f (v) = N on the arc A^Bj so that ^ - p = Z N > 0 drc doc This proves that F 2 must l i e above A-^ provided E 2 is below B£ . If, however, E 2 is above B^ or coincides with it , then T 2 must clearly remain above T£ throughout the interval [ 0 , 21^ ] (except, perhaps, at x = a^ ) because of the arguments just stated above. Now we can perform an exactly analogous construction for the upper half of the phase-plane by tracing T 2 from C 2 in the direction of decreasing t t i l l it intersects the - 54 -v-ax i s at A 2 , r e f l e c t i n g the arc C 2 A 2 of T 2 i n the v -ax i s to obta in the image-arc C - ^ of T£ (see F i g . 12), then continuing T-j_ beyond C-j_ and proving that i t s i n t e r -sect ion F-, with the v-axis must l i e below A 2 . F i g . 12 We can now see that the region R bounded by the arcs A}C]_Gi and A 2 C 2 F 2 of T]_ and T 2 , r espec t ive ly , and by the segments A 2 F ^ and A^F 2 has the property that no t ra jec to ry of (3.2.15) can leave i t as t °° . We can now conclude our proof by applying Lemma 1.4 and the Theorem 1.1 - 55 -in the same fashion as in the preceding Theorem. It can also be shown by considerations similar to those employed in the proof of Theorem III. l that if we replace the condition (3*) by the weaker condition $0=c) a-. W: for \x\ suff. large then the theorem will s t i l l hold for values of u. satisfying Theorem III.3 Consider the differential equation (3.2.18) ~ + [ * C « , * W M - c f C o c ^ ) ] - v ^ C o c V ^ - 0 ^ g . Let the functions f, q» and g be locally Lipschitzian and satisfy the/following conditions (1**) $t*aO > ° ^ 0 r 3 ^ ^ '•> (2**) (ft-*.) -*oo SI V*\ oo ( 3 * * ) J c o ^ - c o (4**) 4^,^) t K) >0 -for bc\ a , all -tf ; 4 0 , V) > - H (M > O) | or V*\ * a , all ^ . (5* *) (f (^, -^ - - ^(^V) ( 6 * * ) W C ^ i ^ M ^ L ^ o r s l 1 C"«>v) . (7**) lirvi t f ( " X i V ) -Vor -eacl* oc Then the differential equation (3.2.Id) will possess at least one non-trivial periodic solution for a l l values of \i . Proof. Let us choose v 0 > 0 such that ( a * * ) & = sup I^X-XJ-NO - c f C x ^ i M <. ~ for al l v ,^ V 2 such that \'v,\ * Wi\ avnd - a « -ac s a This is always possible in view of (7**). After this let us choose x 0 "> 0 such that (9**) KIf~ —. , n N ^ C » M » f l * V r tog KlC^o - a" ) — -e^, — — ( where Finally, increase v 0 , i f necessary, to satisfy the inequalities and ( 1 1 * * ) d, * W0=c>O - M ? c' x>^l <• Y for al l (x, v^ ) and {x^ v 2) such that 1^ ,1 ^ 'Vo Wil "5 *\*0 and - -are ^ =c s a - 57 -This increase in v 0 will not violate the inequalities (g**) and (9**) . The first assertion is evident while the second follows from the fact that the r.h.s. of (9* * ) is equal to 2 a 5 \ and, consequently, decreases as vQ. increases. Having done this, let us consider the point S* 8 5 (-x0,vQ) in the second quadrant of the phase plane. Let us draw through S* a straight line S*B (see Fib. 13) with the slope (3.2.19) m H + M. L ~ V where K2 • sup \ g(x)x| for -x 0 < x < a . The slope of a trajector of (3.2.18) is given by (3.2.20) ^ - - + C « . v ) -R s —— a —- a —— Fig. 13 A comparison between (3.2.19) and (3.2.20) shows that the slope of a trajectory passing through a point on S*B with x € [ -xQ, a] is less than m , so that al l trajectories cross S*B in this interval as it is indicated by arrows on Fig. 13. It is now clear that if we - 5* -consider a t ra jec tory T through B = (a, v 0 + m x c + ma } and trace i t to the l e f t from B then we s h a l l have v > v Q throughout the whole i n t e r v a l ^ - x 0 , a"] . Let us assume that T , when followed to the l e f t , crosses the x -ax is at a point P . The case where no such point ex is t s w i l l be discussed l a t e r . Let us now fo l low the t r a jec to ry T to the r igh t of B . We s h a l l show that i t w i l l cross the x - a x i s . Indeed, i t s slope i s given by (3.2.20) and we see that fo r x "> a and v 0 the f i r s t and the t h i r d terms on the r . h . s . are negative while the second term i s bounded by uL . Using the condi t ion (2*"*) and repeating the arguments employed i n the proofs of the preceding two theorems we conclude that T must cross the x-ax is at a point which we sha l l denote by C . Let us now r e f l e c t the arc PSABC of T i n the x -ax i s to obta in i t s image P S U ' B ' C (see F i g . 14), denoted by T f . We can now make the fo l lowing observation, a l l t r a j ec to r i e s of (3.2.IS) crossing T T between P and A f or between B ' and C do so i n the d i r e c t i o n indica ted by the arrows, v i z . They enter the reg ion R 0 enclosed F i g . 14 by T and T T . To see t h i s we compare the slope of a t ra jec tory given by (3.2.20) with the slope of T* given by - 59 -( 3 . 2 . 2 D j£ = K - = C > - ^ 4- H -q>c-=c,-v) - <%c*y f -Using the condi t ion ( 5 * * ) we see that the difference i s The r . h . s . of t h i s equal i ty i s greater than 2N > 0 for a l l points of T f between P and A' or between B 1 and C and our asser t ion fo l l ows . Let us now consider the t rajectory a, T* of (3.2.18) through the point B ' and fol low i t to the l e f t . At B f the t ra jec to ry TM must enter the reg ion RQ described before, but between B* and A' the t r a j ec to ry T" F i g . 15 may cross and recross T f any number of t imes. We s h a l l assume that T" f i n a l l y leaves R 0 at a point F (see F i g . 1 5 H The other p o s s i b i l i t y i s of no in teres t since i t would mean that T" never leaves R, as t -*• °° and the proof of the theorem w i l l f o l low immediately. Let doc and ^ ~ denote the ordinates and the slopes of TT and T" , r e spec t ive ly . Using (4 ) - 60 -and (3 * *) we get where K = sup | g(x)x \ for -a < x < a . Denoting v 1 - v n z 2M + ~ = o* , ~ \ = p the above inequality becomes <k^. > _ q- _ p«y Adding' pty to both sides and multiplying by exp px • dx we obtain Integrating between -a and x(F) and remembering that at x = x(F) \)r = 0 we get - ^ C - a > (-pa) ^ £ t * * ^ ^ - ? a ) " " P " * ^ ] V or y C- aN; ^ — [e<c^ > p {^O) *• a"} - l"] P Substituting the expressions for T and p and using x(F) + a < 2 a we finally obtain y eoqp — - 1 = £ - 61 -This inequality states, in fact, that the vertical distance ATAn (see Fig. 15) does not exceed e . Let us now follow T n beyond Att to the left. We shall show that T" must intersect the arc S'A' of T» . To prove' this we make the hypothesis that this is not so, that is T" remains below T* throughout the interval [-•ac©,-a] . As before, we compute the difference . . ^ -and apply the conditions ( 4 * * ) and (11* *) to derive the following inequality, valid in C'^^-a"] '• (3.2.22) **' _ d*v 2.N - - ~ h tv'-v") At x = -a we have (v1 - v") < e and, therefore drx. dot. v = N > O Z, 2. where the condition (10 ) was applied. Since J ~ (v* - vM) *> 0 at x - - a , the difference (vT - v") will start decreasing as we follow T" to the left. This in turn will increase ^ (vr - v") as follows from (3.2.22) and it is obvious that the inequality (3.2.23) jt^'-r*") - VO will hold throughout the interval j ^ - ^ * , - 3 ! • Integrating (3.2.23) from x = -xQ to x = - a and putting - 62 -v ' ( - a ) - v"(-a) < e we obtain v'(-ocO - v * £ - KH «o -a") Za£ (2M + H as fo l lows from ( 9 * * ) . This r e su l t contradicts our hypothesis I t fo l lows that there ex i s t s a point Q on the arc A'S* of TT where T" crosses T ' and re-enters the region R 0 . Let R be the region bounded by the arc PSABC of T , by arc CB» of T' , by arc B'Q of T" and, f i n a l l y , by arc QP of T f (see F i g . 16). From the considerations stated before i t i s c lear that no trajectory of (3.2.18) can leave R . We may now apply Lemma 1.4 and the Theorem 1.1 i n the usual fashion and the proof i s complete. I t remains to consider the case when the point P does not e x i s t . As before, we construct the image T 1 of our t ra jec tory T (see F i g . 17) and consider the t ra jec tory T" through B ' . Using the condi t ion (2* * ) and repeating the arguments already used before we can e a s i l y show that T" p -I' - 63 -V* F \ V yf T* ft" J L r must in tersec t the x -ax i s at some point P" . Let us fo l low T" beyond P" . Except for a neighbourhood of P" the slope of T" i s bounded and therefore T n must cross the v -ax i s and enter the f i r s t quadrant. From F i g . 17 there on we use the same argument as for the t ra jec tory T to show that T" w i l l cross the x-ax is at some point F" . I t i s easy to see that F" must l i e to the l e f t of C since otherwise T" w u l d have had to cross T which i s impossible. The reg ion RT bounded by the arc B T A"P T , F" of T" , by the segment F"C and by the arc CB' of T' i s equivalent to R as far as our purpose i s concerned and the theorem fol lows i n t h i s case as w e l l . As before, we also have the fo l lowing c o r o l l a r y : Coro l l a ry . I f the condi t ion (2**) i s replaced by c^t-x) *. It for | x | suff . large then the theorem w i l l hold for the values of \x s a t i s fy ing L - 64 -3.3 A case when f (x, v) is not bounded Theorem III.4 Consider the differential equation (3.3.D ^ + {Coc.vV-* + «%c~Voc - O ~ Let the functions f and g be locally Lipschitzian and satisfy the following conditions: (a) <%L~) -*0 (b) <frC-«) = c^L-x) (c) ^Co,o)<0 (d) ^C-x.,^) M > o except, perhaps, in some closed neighbourhood R of the origin and in the region for which 0 <~ <- -rr\ where W is a positive constant (e) there exists cc,>0 and b ^ = O such that Then the equation (3.3.1) possesses a non-trivial periodic solution. Proof. With x^ and b defined as in (e), let us denote rr 5- = a 1 and choose -Xo > ^ large enough to satisfy N T T the inequality (3.3.2) 2 - oc,2" + Wl) >0 - 65 -Let us consider the curves defined by (3.3.3) where and X i s a non-negative parameter. These curves are t r a j ec to r i e s of the d i f f e r e n t i a l equation and, because of the condit ions (a) and (b) they are closed curves around the o r i g i n , symmetric with respect to both coordinate axes. They w i l l be c a l l e d "X-curves". I t i s easy to see that we can always choose x£ > 0 large enough so that the region R w i l l l i e e n t i r e l y w i t h i n the i n t e r i o r of the X-curve given by Let us now denote x„ = max ( x 1 . x") and consider the point o o o A = ( -x f t ) 0) i n the phase-NT plane (see F i g . I S ) . Let us draw through A a X-curve L and a t ra jectory T of D' = ( -x 0 , m x 0 ) and l e t H « ( x 0 , 0 ) . We s h a l l fo l low point where L in tersec ts the ray v - - m x and l e t ( 3 . 3 . 1 ) . Let B be the F i g . IS - 66 -the trajectory T in the direction of increasing t and show that it must cross the v-axis and then descend and cross the x-axis at some point P lying between the origin and the point H . As we start from A , the trajectory T must remain to the right of the vertical line - DA and so must cross the ray v = - inn x at some point C . If C is closer to the origin than B (or coincides with B) then our assertion follows immediately since on the arc HI of L we have f (x, v) 2- N > 0 and it is easy to see that T cannot leave the interior of L via the arc BH. We shall, therefore, assume that C lies between D and B . Then, clearly X(C) < \(D) . We also have Let us follow the trajectory T beyond C . Because of the conditions (a) and (d) it easily follows that T must cross the v-axis, then descend and cross the x-axis at some point P . From an argument employed before it is clear that if T crosses the arc BH of L then P must l ie X(O) - A ( r V ) - £ [<v(D)"]Z + G(~o^ - Gteo) = and, therefore, (3.3.4) Vie) - MR) <, ^ l x 0 l - 67 -to the left of H and our assertion will be proven. To demonstrate that this will indeed be the case we shall make the contrary hypothesis and assume that T lies above the arc BH of L and crosses the line x = xQ at some point F lying above H . This hypothesis is illustrated on Fig. 18. Let E be the intersection of T with the v-axis. From Lemma 1.3 and the condition (d) i t readily follows that (3.3.5) M E ) < : M c ) Let v^ = v-^ (x) and v 2 s v2(x) represent the equations of curves T and L^ respectively, in the interval [ 0, x 0] • Using our hypothesis we have (3.3.6) > N • j •»•, • doc - M • J ^ • <^ - N -Si o o where si is the area under the curve L in the first quadrant. We are now going to show that To see this let us consider the ellipse L 1 defined by _ (3.3.S) -vr-s = ^ ( * ) = a* U -ac.*- « l - 63 -The are under L» in the first quadrant is equal to fa 1 z - 4 - • ceo ; hence, i f we show that v 2 •a v^ for 0 < x < , the inequality (3.3.7) will follow: Indeed, we have ^ = 2XO0 - 2 • GCoc) = 2- GC-xo) - 2- GC*) Together with (3.3.S) this yields 0 We have to show that 8 > 0 for a l l x e £ 0, xQJ . If x € [ x j^ xQ^ this follows immediately from the condition (e) which can be rewritten as 0^=0 - a* > ^ ^ O for it t-ac, If, however, x fe [ 0, x^ "] then we can write •ac -ar, - 69 -where the conditions (e) and (3•3.2), as well as the fact that XQ > x£ , were used. The validity of (3.3.7) is, therefore, established. Now using (3.3.4), (3-3.5), (3.3.6) and (3.3.7) we obtain Since X(A) - \(H) it follows that \(F) < \(H) . But, according to our hypothesis, F lies above H on the same vertical, so that X(F) >• \(H). The contradiction proves that the hypothesis is false and, consequently, that T must intersect the arc BH of L . It follows that the point P lies to the left of H . If we now consider the trajectory T* of (3.3*1) passing through H then we can show, in the same fashion as above, that T' must cross the v-axis and then ascend and cross the x-axis at a point Q lying to the right of A (see Fig. 19). The region bounded by the arcs AP and HQ of T and T* , respectively, and the segments AQ and PH possesses the property that no trajectory can leave it as t -*• 0 0 • The existence of a non-trivial periodic solution follows readily now by application of Lemma 1*4 and the Theorem 1.1. BIBLIOGRAPHY I. Bendixson, Acta. Math. 24, (1901), 1-88. A. de Castro, Boll. Un. Mat. Italiana III, 8, (1953). 26-29. A. de Castro, Boll. Un. Mat. Italiana III, 9, (1954), 369-372. A.V. Dragilev, Prikl. Mat. Mekh. 16, (1952), 85-88. A.F. Filippov, Mat. Sbornik 30(72), (1952), 171-180. S. Lefschetz, Lectures on Diff. Equations, Princeton Univ. Press 1948. N. Levinson and O.K. Smith, Duke Math. J. 9, (1942), 382-403. A. Lienard, Rev. gen. d f Electric it e* 23, (1928), 901-946. E. McHarg, J . London Math. Soc. 22(86), (1947), 83-85. H. Poincare, Oeuvres I, pp. 3-222. B. van der Pol, Phil. Mag. II, 7, (1926) 97^-992. M.I. Yel'shin, Mat. Sbornik 34 (1954), 323-330.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Existence of periodic solutions of certain non-linear...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Existence of periodic solutions of certain non-linear differential equations Butkov, Eugene 1956
pdf
Page Metadata
Item Metadata
Title | Existence of periodic solutions of certain non-linear differential equations |
Creator |
Butkov, Eugene |
Publisher | University of British Columbia |
Date Issued | 1956 |
Description | The theory of Poincaré and Bendixson is applied to establish the existence of periodic solutions of the differential equation ẍ + f(x, ẋ)ẋ + g(x) x = 0 One part of the work is concerned with those equations which can be considered as arising from small perturbations of other equations of the same type, already possessing periodic solutions. Two existence theorems are demonstrated and the stability and uniqueness of the periodic solutions is also discussed. The other part contains several, theorems stating sufficient conditions for existence of periodic solutions which cannot be treated by perturbation methods. |
Subject |
Differential equations |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2012-02-02 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0080654 |
URI | http://hdl.handle.net/2429/40444 |
Degree |
Master of Arts - MA |
Program |
Mathematics |
Affiliation |
Science, Faculty of Mathematics, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
AggregatedSourceRepository | DSpace |
Download
- Media
- 831-UBC_1956_A8 B9 E9.pdf [ 4.75MB ]
- Metadata
- JSON: 831-1.0080654.json
- JSON-LD: 831-1.0080654-ld.json
- RDF/XML (Pretty): 831-1.0080654-rdf.xml
- RDF/JSON: 831-1.0080654-rdf.json
- Turtle: 831-1.0080654-turtle.txt
- N-Triples: 831-1.0080654-rdf-ntriples.txt
- Original Record: 831-1.0080654-source.json
- Full Text
- 831-1.0080654-fulltext.txt
- Citation
- 831-1.0080654.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
https://iiif.library.ubc.ca/presentation/dsp.831.1-0080654/manifest