Studies on Rayleigh Equation
|
|
- Howard Short
- 6 years ago
- Views:
Transcription
1 Johan Matheus Tuwankotta Studies on Rayleigh Equation Intisari Dalam makalah ini dipelajari persamaan oscilator tak linear yaitu persamaan Rayleigh : x x = µ ( ( x ) ) x. Masalah kestabilan lokal di sekitar titik tetap dan hubungannya dengan paramater µ di bahas. Juga akan dibuktikan eksistensi selesaian periodik yang dalam hal ini juga merupakan limit cycle dengan menggunakan teorema Poincare-Bendixson. Selesaian periodik tersebut akan dihampiri dengan menggunakan metode Poincare-Linsted. Sebagai perbandingan juga di hitung selesaian numerik dan dibandingkan dengan selesaian asimtotik Abstract In this paper we consider a type of nonlinear oscillator known as Rayleigh equation, i.e. x x = µ ( ( x ) ) x. We study local the stability of the fixed point in the presence of positive parameter µ. We are also looking at the existence of the periodic solution which is also a limit cycle in this equation. Using Poincare-Bendixson s theorem we proof that the limit cycle exists. For small values of the parameter, we use asymptotic analysis to approximate the periodic solution. The method we apply in this paper is Poincare-Linsted method. We calculate also the numerical solution for the system and compare the result with the asymptotic. INTEGRAL, vol. 5 no., April 000
2 I. Introduction Consider a system of differential equations u = f t, u,µ.. () where u R, f is a vector valued function in R and µ is a parameter. We assume that solution exists and is unique with given a initial value. This system is also known as a planar system. The system depended on the parameter µ. For a fixed µ, let u µ (t) = ϕ µ (t) be a solution of the system (the subscript shows the dependency on the parameter). This vector valued function ϕ µ parameterize by t defines a flow in the two dimensional plane u = (x, y). The curve of points (x, y) generated by ϕ µ (t) when we let t run from - to is called the trajectory. A collection of trajectories of different initial value is calling the phase portrait of the system. A simple question that arises is; how to know the phase portrait of a system. The best way to have the phase portrait of a system is to calculate the general expression of the solution of the system. Unfortunately, this is not an easy thing to do especially when f is non-linear. The second best way is to calculate the integral of the system. An integral (in this case) is a real valued, two variables function F(x, y) such that d dt ( F( x, y)) ϕ = 0, which is the derivative of F with respect to time evaluated at a particular solution zero. Consequently, F(x, y) = C define a curve in the (x, y) plane such that the flow of system () will map the curve to it self. Such a curve is called invariant curve of the system. In elementary courses on differential equation such a curve is called implicit solution. (Unfortunately not every system has an integral). We remark that if we can calculate one of those (either the solution or the integral) we will have a global picture of the flow. It means that the analysis is valid for all time (t) and everywhere in the (x, y) plane. That is why these two methods are also known as global analysis technique. Naturally, when the global picture cannot be achieved, we can try to have a local picture. We can restrict our analysis in a small area around some interesting location. The question is, where the interesting location is. Some of the locations, which are considered to be interesting, are around a fixed point or around a periodic orbit. A fixed point (also known as a constant of motion point, or critical point, or singular point) is a point (x 0, y 0 ) in the plane such that f(x 0, y 0 ) = 0. A periodic orbit is a solution ϕ µ (t) where a real number T such that ϕ µ (t T) = ϕ µ (t) exists. Now supposed we have a periodic orbit that start at a particular point. After T time it will be back at that starting point. It means the trajectory will be a closed curve. The converse is also true if we assume that T is finite. The most common tool to understand locally the flow of the system is linearization. We expand f to its Taylor series around a particular solution, i.e. u = f( ϕµ )( u ϕµ )! () where the dots represent the higher order terms and represent the partial derivative with respect to spatial variables. And then we can use linear analysis to obtain the information about the flow. See [], [4] or [5] for example. Our goal in this paper is to analyze the non-linear oscillator known as the Rayleigh equation, i.e. x x = µ ( ( x ) ) x.(3) where µ is a positive parameter. We give the analysis in the neighborhood of INTEGRAL, vol. 5 no., April 000
3 the fixed point. We will also show that there exists a periodic orbit and that it is stable. This stable periodic orbit is also known as the stable limit cycle. Due to the fact that the analytic calculation of the limit cycle is complicated, we will construct an approximation of it s for a small parameter using Poincare-Linsted method. For a large parameter we can construct the approximation using boundary layer technique or singular perturbation technique (see [3]). These technique is mathematically non trivial so we will skip them. In this paper we will also construct the numerical comparison of the analysis. of the solutions being π-periodic. Thus, the origin (0,0) is a center point if µ = 0. In general this statement is not always true. II. Fixed Point Analysis To calculate the fixed point of the Rayleigh oscillator, we transform the system into a system of first order differential equations. This is done by setting x = x and y = x'. The Rayleigh equation then becomes x = y..(4) y = x µ ( y ) y It is easy to see that the fixed point of (4) is only (0,0). In general, for a system like in () the fixed point is also depended on µ. The linearized system of (3) written in matrix form is: ' x 0 x =...(5) ' y µ y From linear analysis we know that the stability of the fixed point (0,0) dependeds on the eigen-values of the linearized system (5). The information on the stability also gives information on the flow it-self around the point. In this case we have the eigen values µ ± µ 4 λ, =...(6) For µ =0 we see that the eigen-values are purely imaginary λ, = ± i. This is clear since if µ = 0, what we have is just a linear harmonic oscillator with all Figure : Phase portrait for µ = 0. The fixed point in this case is a center point and all solutions are periodic of period π. Obviously the eigen-value of a non linear system should be non linearly dependent on the system. Thus if we linearize locally, we restrict the domain so that the non- linear effect can be considered as a small perturbation. If the real part of the eigen value is nonzero then we can choose the domain small enough so that the non-linear perturbance will be small enough. It implies that the real part will still be the same in sign (see [4],[5] for details). This is not the case for zero real part of the eigen-value. In the case of all the real parts of the eigen-values is zero we have to consider a higher order term on the expanded system. In the case that there is at least one is non zero, we can use Center Manifold Theorem (see [] for details). This case does not arise on this problem so we omit it. For 0 < µ < we can write the eigen µ ± ω i values as λ, =, ω < µ. In INTEGRAL, vol. 5 no., April 000 3
4 this case, the flow around the origin is spiraling outward the fixed point. Thus, if we start on a point close to the origin the trajectory will move around the origin with increasing radius. It is interesting that most of the periodic orbits of the case µ=0, break up. In fact there is only one periodic orbit that survives and forms a limit cycle. The next case is if µ =. In this case we have only one eigen-value. This eigen value corresponds to single a eigen-vector k. Thus the origin will look like a source which is an improper node. The phase flow is flowing out from the origin along one direction which can be approximate by the eigenvector. Figure 3: Phase portrait for µ =. It represents a improper node where the flow is going outward along one direction. For µ >, we have two different real eigen-values. Thus we will have two linearly independent eigen-vectors. Since the eigen values are different, then the origin will also be improper node sources. The phase flow is flowing out of the origin along two direction only. The picture in figure () until (4) are drawn using Maple V. It is clear that the analysis of the origin coincides with the numerical result using Maple V. Figure : Phase portrait for µ = 0.5. The flow around the fixed point is spiraling outward the fixed point. Figure 4: Phase flow for µ =.5. III. Existence of the Periodic Solution To proof the existence of the periodic orbit, we will apply a fundamental 4 INTEGRAL, vol. 5 no., April 000
5 theorem of Poincare-Bendixson. This theorem will be stated without proof. Theorem: Poincare-Bendixson If D is a closed bounded region of the (x, y)- plane and a solution u(t) of a nonsingular system () is such that u(t) D, then the solution either a closed path, approaches a closed path, or approaches a fixed point. The theorem in other word says that if we have a closed and bounded region which the flow of the system in () is flowing into the region, then there is a periodic solution or a fixed point in the region. We will apply this theorem to a general type of oscillator with damping by constructing a region which is invariant to the flow. After that we will relate the Rayleigh equation to the general equation we have. Consider now an oscillator equation known as Lienard equation, i.e. x f ( x) x x = 0..(7) We define x F( x) = f ( s) ds 0 and assume it to be an odd function. Futher assume that for x > β (a real positive number) F is positive, goes to infinity for x goes to infinity, and monotonically increasing, while for 0 < x < α F is negative. We transform (7) to a first order system using transformation (x, x ) (x, y = x F(x)). Using this transformation we have x = y F( x)..(8) y = x. Now define a transformation to polar coordinate by R = ( x y ). Consequently we have R = -xf(x). It implies that for -α < x < α, R 0. Thus the flow on the circle domain with radius α is flowing out of the domain. What we need to construct is a domain where the flow is flowing into the domain. See figure (5). Before we start constructing the domain, we first make some important remark on the gradient field of the flow. From (8) we can calculate dy x dx = y F ( x. It ) implies that on the y axis the tangent is horizontal and on the curve y = F(x) the tangent is vertical. One can clarify easily that if we start at an initial value say A we will have a trajectory as in figure (5) (this can be easily done by considering the negative or positiveness of x or y at certain location). The reflection symmetry of the system gives the negative trajectory. Our purpose is now to proof that R(A) > R(D). This can be done by considering the line integral RD RA = dr. We split the ABCD curve ABCD into three segments AB, curve ABCD into three segments AB, BC, CD. We also have two expressions for dr, i.e. Figure 5: Constructing the invariant domain. xf x dr = dx or dr = F( x) dy. y F( x) INTEGRAL, vol. 5 no., April 000 5
6 We first consider segment AB and CD. If A (= (0, y A )) is high enough then we know that y-f(x) is large while -xf(x) is bounded. Thus for y A going to infinity, the integrals dr and dr go to AB CD zero. It also means that the integral is dominated by the line integral in the segment CD. In segment CD we have BC C dr = F( x) dy. Since we assume B that for x > β, F is positive and monotonically increasing, we see that the integral will be negative and monotonically decreasing (it is clear since the integration is from positive values to the negative values of y). Also it is continuosly dependent on y A. Thus if y A goes to infinity the integral goes to infinity. It means that we can choose y A large enough such that R(A) > R(D). In the same way we can show R(-A) < R(-D). Now consider the domain in figure (5) which is bounded by the circle with radius α, the two trajectories and the line segment connecting the trajectories. This domain defines an invariant closed and bounded domain. Moreover, the flow is going into the domain. Hence the Poincare- Bendixson s theorem applies. We know that the only fixed point of the Lienard s system is the origin. The conclusion is that we have at least one periodic solution. Moreover, since F is monotonically increasing and consequently the integral in the segment BC is monotonically decreasing, the condition that if α=β then leads to precisely one periodic solution. The question arises is how to apply this result to the Rayleigh equation. If we differentiate The Rayleigh equation, we have x µ ( 3( x ) ) x x = 0....(9) Define z = 3 x and transform (7) into z µ z z z = 0...(0) The solution of (9) is just a solution of ordinary Rayleigh equation with an additional constant. Thus by setting the constant to be zero, we have the original solution. Since (0) satisfies the assumption we now find that there exists at least one periodic solution. We can proof that in the case of (0), α=β. Thus we have a unique periodic solution. IV. Asymptotic Approximation The next aim is to approximate asymptotically the periodic solution of Rayleigh equation. We will apply the Poincare-Linsted Method to approximate the periodic solution. We will first give a short introduction to the method. Details of the method can be found in [4]. Consider a second order differential equation of the form x x = εf( x, x, ε) () It is easy to see that for ε = 0, all solutions are π-periodic. For ε 0 we assume that there exists a periodic solution with period T(ε) starting at initial values x(0) = a(ε) and x'(0) = 0. The fact that the period and the location or the periodic orbit are dependent on ε is natural. Define a time transformation θ = ωt such that in this new time variable, the period of the periodic orbit is π. We write ω - = - εη(ε). Obviously, the periodic solution is also dependent continuesly on ε. Thus we assume that we can write the periodic solution as x( θ) = x0 εx ε x!...() With the new time variable we can write () in the form of (the dot represents the derivative with respect to θ) 6 INTEGRAL, vol. 5 no., April 000
7 x" x"" x = ε ηx ( ) f x,, εη εη ε. If we write the right hand side as εgxx (, ", ε) and transform it to its integral form, then the periodicity condition of the solution leads to two equations, i.e. π F ( a, η) = cos( τ ) g x, x", ε, ηdτ = 0 F ( a, η) = sin( τ ) g x,", x ε, ηdτ = 0 0..(3) π 0 Implicit function theorem gives the conditions for the existence of a nontrivial solution of (3) in the neighborhood of ε = 0, provided ( F, F) 0.(4) a, η For the case that the existence condition is satisfied we will have a unique periodic solution. This is in agreement with the previous analysis (we have proven that the Rayleigh system has a unique periodic solution). Note that if (4) fail to hold, it does not mean that there is no periodic solution. In many cases such as in hamiltonian system, the periodic solution is not isolated so that apriory we know that the condition fails to hold. It is instructive to apply the method. Note that this is an asymptotic method so that it is valid for µ 0. Thus we take µ equal to small parameter ε. The calculations are rather routine and lengthy. We write the result of the calculation up to order 4. The calculation is done using Maple release V. The periodic solution is approximated by () where x 0 ( θ ) = x ( θ ) = x ( θ ) = x 3 ( θ ) = 78 x 4 ( θ) = cos( θ ) 6 cos( θ ) 3 sin( 3 θ ) 3 cos( 3 θ ) 88 3 cos( 5θ ) 3 sin( 7 θ) cos( θ) 3 cos( θ) 6 3 cos( 7θ) sin( θ ) sin( θ ) 3 cos( 5 θ ) 3 sin( 3θ ) 3 cos( 5θ) 3 cos( 9θ) and θ = ωt. Obviously this is a nontrivial thing to do but nevertheless, as we noted above it is instructive to do it. Furthermore, we also calculate the period, i.e T = π πε πε O( ε ) The location of the periodic solution is a = ε 743 O ε ( ε ) We will check on this result with a numerical integration of the system using MatLab 5.. V. Numerical Result We will now check the approximation above using a numerical integration. We plot both of the result in one picture so that we can immediately compare the result. We use only the first approximation x(θ) = x 0 (θ) and O(ε 5 ) approximation for ω. This has the advantage of a very long time-scale. It means that the approximation is valid for quite a long time and in this case until ε -4. The numerical integration is done using build in integrator in MatLab 5. for non-stiff system, i.e. ODE45. This integrator uses Runge-Kutta scheme for INTEGRAL, vol. 5 no., April 000 7
8 integration. We plot the result in figure (6), (7), (8) and (9). The curve plotted by symbol o represent the numerical solution and the line curve is for the asymptotic approximation. In figure (6) and (7) we take ε=0.05. Thus the approximation is still good until 6,000 second. For numerical integration this long time integration is not recommendable. We have to worry also about the numerical error. Thus we integrate with initial values x(0) = a and x (0) =0 for 0 second only. Obviously we can expect a very good result in the comparison. In figure (8) and (9), we increase the value of ε to 0.5. We integrate with the same initial values for 80 second. This is already longer than the timescale of validity of the approximation. Thus we expect to see an O(ε) deviation on the picture. In figure (8) the deviation is not very clear but in figure (9) it is. Figure 6: The comparison between numerical result and asymptotic approximation for ε=0.05. This picture represents the time evolution of the periodic solution. Figure 8: The comparison of the time evolution for ε=0.5. Figure 9: The comparison in the (x, x )- plane for ε=0.5. Figure 7: The comparison between numerical result and asymptotic approximation for ε=0.05 on the (x, x )-plane. To make the error even clearer we plot the error defined as the difference between the numerical solution and the asymptotic solution in figure (0) and (). 8 INTEGRAL, vol. 5 no., April 000
9 Beside the method we have used here there is also a simpler method to get the approximation. The method is called averaging method. This is a very natural method introduced intuitively by Lagrange et. al.. Unlike the Poincare- Linsted method, this method is more general since it can be applied to approximate non-periodic solutions. The disadvantage of the averaging method is that, to get a higher order approximation is non-trivial, while for the Poincare-Linsted method it is routine. Also in extending the timescale, using averaging method is rather difficult. For a higher dimension case, where everything is restricted, averaging is easier to apply. Figure 0: The error on x. Figure : The error on x. VI. Remarks and Acknowledgment We have shown an example of asymptotic approximation of a periodic solution of a planar system. This Poincare-Linsted method can also be extended to higher dimension system. We note that it is not trivial to do so. We would like to express our gratitude to Santi Goenarso for her diversified contribution. We also like to thank our student Maynerd Tambunan for providing a comparison work for the Maple calculation. References: [] Boyce, W.E., Di Prima, R.C., Elementary Differential Equations and Boundary Value Problems, John Wiley & Sons, Inc., New York et. al., 99. [] Carr, J., Applications of Center Manifold Theorem, Applied Mathematical Science 35, Springer- Verlag, New York, 98. [3] O Mailley Jr., R. E., Singular Perturbation Methods for Ordinary Differential Equations, Applied Mathematical Science 89, Springer- Verlag, New York, 99. [4] Verhulst, F., Nonlinear Differential Equations and Dynamical Systems nd ed., Springer-Verlag, Berlin, 996. [5] Wiggins, S., Introduction to Nonlinear Dynamical System, Text on Applied Mathematics, Springer-Verlag, New York, 990. [6] Guckenheimer, J., Holmes, P., Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields, Applied Mathematical Science 4, Springer-Verlag, New York, 983. Author Johan Matheus Tuwankotta is a lecturer at the Mathematics Department ITB. address: tuwankotta@math.uu.nl Received August 5, 999; revised September INTEGRAL, vol. 5 no., April 000 9
A Note on Singularly Perturbed System
Jurnal Matematika dan Sains Vol. 7 No. 1, April 00, hal 7-33 A Note on Singularly Perturbed System Hengki Tasman 1), Theo Tuwankotta ), and Wono Setya Budhi 3) Department of Mathematics, ITB E-mails: 1)
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationLectures on Periodic Orbits
Lectures on Periodic Orbits 11 February 2009 Most of the contents of these notes can found in any typical text on dynamical systems, most notably Strogatz [1994], Perko [2001] and Verhulst [1996]. Complete
More informationENGI Duffing s Equation Page 4.65
ENGI 940 4. - Duffing s Equation Page 4.65 4. Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More informationNonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53
More informationCanards at Folded Nodes
Canards at Folded Nodes John Guckenheimer and Radu Haiduc Mathematics Department, Ithaca, NY 14853 For Yulij Il yashenko with admiration and affection on the occasion of his 60th birthday March 18, 2003
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationMathematics Department, UIN Maulana Malik Ibrahim Malang, Malang, Indonesian;
P r o c e e d i n g I n t e r n a t i o n a l C o n f e r e n c e, 0 3, *, * * - * * The 4 th Green Technology Faculty of Science and Technology Islamic of University State Maulana Malik Ibrahim Malang
More informationDynamics of a mass-spring-pendulum system with vastly different frequencies
Dynamics of a mass-spring-pendulum system with vastly different frequencies Hiba Sheheitli, hs497@cornell.edu Richard H. Rand, rhr2@cornell.edu Cornell University, Ithaca, NY, USA Abstract. We investigate
More informationA Study of the Van der Pol Equation
A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple
More informationLecture - 11 Bendixson and Poincare Bendixson Criteria Van der Pol Oscillator
Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 11 Bendixson and Poincare Bendixson Criteria
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
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 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 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 information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
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 informationChapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics
Chapter 23 Predicting Chaos We have discussed methods for diagnosing chaos, but what about predicting the existence of chaos in a dynamical system. This is a much harder problem, and it seems that the
More informationPeriod-doubling cascades of a Silnikov equation
Period-doubling cascades of a Silnikov equation Keying Guan and Beiye Feng Science College, Beijing Jiaotong University, Email: keying.guan@gmail.com Institute of Applied Mathematics, Academia Sinica,
More informationThe Big Picture. Discuss Examples of unpredictability. Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986)
The Big Picture Discuss Examples of unpredictability Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986) Lecture 2: Natural Computation & Self-Organization, Physics 256A (Winter
More informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
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 information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
More information4 Second-Order Systems
4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
More informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x
More informationSection 9.3 Phase Plane Portraits (for Planar Systems)
Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Half-line solutions Unstable solution Stable
More informationWHAT IS A CHAOTIC ATTRACTOR?
WHAT IS A CHAOTIC ATTRACTOR? CLARK ROBINSON Abstract. Devaney gave a mathematical definition of the term chaos, which had earlier been introduced by Yorke. We discuss issues involved in choosing the properties
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More information27. Topological classification of complex linear foliations
27. Topological classification of complex linear foliations 545 H. Find the expression of the corresponding element [Γ ε ] H 1 (L ε, Z) through [Γ 1 ε], [Γ 2 ε], [δ ε ]. Problem 26.24. Prove that for any
More informationTWO ASPECTS OF A GENERALIZED FIBONACCI SEQUENCE
J. Indones. Math. Soc. Vol., No. (05, pp. 7. TWO ASPECTS OF A GENERALIZED FIBONACCI SEQUENCE J.M. Tuwankotta Analysis and Geometry Group, FMIPA, Institut Teknologi Bandung, Ganesha no. 0, Bandung, Indonesia
More informationSYNTHESIS OF SINEWAVE OSCILLATOR BASED ON THE MODIFIED VAN
ELECTRONICS 7 19 1 September, Sozopol, BULGARIA SYNTHESIS OF SINEWAVE OSCILLATOR BASED ON THE MODIFIED VAN DER POL EQUATION USING MELNIKOV THEORY Zhivko Dimitrov Georgiev 1, Todor Georgiev Todorov, Emil
More informationExponentially small splitting of separatrices of the pendulum: two different examples. Marcel Guardia, Carme Olivé, Tere M-Seara
Exponentially small splitting of separatrices of the pendulum: two different examples Marcel Guardia, Carme Olivé, Tere M-Seara 1 A fast periodic perturbation of the pendulum We consider a non-autonomous
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 information1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),
1.7. Stability and attractors. Consider the autonomous differential equation (7.1) ẋ = f(x), where f C r (lr d, lr d ), r 1. For notation, for any x lr d, c lr, we let B(x, c) = { ξ lr d : ξ x < c }. Suppose
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More informationENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)
ENGI 940 4.06 - Linear Approximation () Page 4. 4.06 Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationSample Solutions of Assignment 10 for MAT3270B
Sample Solutions of Assignment 1 for MAT327B 1. For the following ODEs, (a) determine all critical points; (b) find the corresponding linear system near each critical point; (c) find the eigenvalues of
More informationLecture 38. Almost Linear Systems
Math 245 - Mathematics of Physics and Engineering I Lecture 38. Almost Linear Systems April 20, 2012 Konstantin Zuev (USC) Math 245, Lecture 38 April 20, 2012 1 / 11 Agenda Stability Properties of Linear
More informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
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 informationBoulder School for Condensed Matter and Materials Physics. Laurette Tuckerman PMMH-ESPCI-CNRS
Boulder School for Condensed Matter and Materials Physics Laurette Tuckerman PMMH-ESPCI-CNRS laurette@pmmh.espci.fr Dynamical Systems: A Basic Primer 1 1 Basic bifurcations 1.1 Fixed points and linear
More informationObservations on the ponderomotive force
Observations on the ponderomotive force D.A. Burton a, R.A. Cairns b, B. Ersfeld c, A. Noble c, S. Yoffe c, and D.A. Jaroszynski c a University of Lancaster, Physics Department, Lancaster LA1 4YB, UK b
More informationMath 273 (51) - Final
Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationTHE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A 2 2 SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS
THE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS Maria P. Skhosana and Stephan V. Joubert, Tshwane University
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
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 informationCANARDS AT FOLDED NODES
MOSCOW MATHEMATICAL JOURNAL Volume 5, Number 1, January March 2005, Pages 91 103 CANARDS AT FOLDED NODES JOHN GUCKENHEIMER AND RADU HAIDUC For Yulij Ilyashenko with admiration and affection on the occasion
More informationNonlinear Oscillators: Free Response
20 Nonlinear Oscillators: Free Response Tools Used in Lab 20 Pendulums To the Instructor: This lab is just an introduction to the nonlinear phase portraits, but the connection between phase portraits and
More informationCalculus and Differential Equations II
MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,
More informationENERGY DECAY ESTIMATES FOR LIENARD S EQUATION WITH QUADRATIC VISCOUS FEEDBACK
Electronic Journal of Differential Equations, Vol. 00(00, No. 70, pp. 1 1. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp ENERGY DECAY ESTIMATES
More informationWIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY
WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY J.M. TUWANKOTTA Abstract. In this paper we present an analysis of a system of coupled oscillators suggested
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS
This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]
More informationCALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD
Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen
More informationAn Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems
An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The Hartman-Grobman and Poincaré-Bendixon Theorems
More informationSymmetry Properties of Confined Convective States
Symmetry Properties of Confined Convective States John Guckenheimer Cornell University 1 Introduction This paper is a commentary on the experimental observation observations of Bensimon et al. [1] of convection
More informationComplicated behavior of dynamical systems. Mathematical methods and computer experiments.
Complicated behavior of dynamical systems. Mathematical methods and computer experiments. Kuznetsov N.V. 1, Leonov G.A. 1, and Seledzhi S.M. 1 St.Petersburg State University Universitetsky pr. 28 198504
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More information(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,
2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationStability Analysis for ODEs
Stability Analysis for ODEs Marc R Roussel September 13, 2005 1 Linear stability analysis Equilibria are not always stable Since stable and unstable equilibria play quite different roles in the dynamics
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 informationQuestion: Total. Points:
MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make
More informationModule 2: Reflecting on One s Problems
MATH55 Module : Reflecting on One s Problems Main Math concepts: Translations, Reflections, Graphs of Equations, Symmetry Auxiliary ideas: Working with quadratics, Mobius maps, Calculus, Inverses I. Transformations
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 informationPeriodic Solutions Averaging Methods in Nonlinear Ordinary Differential Equations
Volume: 4 Issue: 7 23 3 Periodic Solutions Averaging Methods in Nonlinear Ordinary Differential Equations Dr. V. Ramadoss, M.Sc., M.Phil., Ph.D.,* & D.Sinduja** *Professor, Department of Mathematics, PRIST
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 informationHello everyone, Best, Josh
Hello everyone, As promised, the chart mentioned in class about what kind of critical points you get with different types of eigenvalues are included on the following pages (The pages are an ecerpt from
More informationTable of contents. d 2 y dx 2, As the equation is linear, these quantities can only be involved in the following manner:
M ath 0 1 E S 1 W inter 0 1 0 Last Updated: January, 01 0 Solving Second Order Linear ODEs Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections 4. 4. 7 and
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More information5.2.2 Planar Andronov-Hopf bifurcation
138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More information18.02 Multivariable Calculus Fall 2007
MIT OpenourseWare http://ocw.mit.edu 8.02 Multivariable alculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8.02 Lecture 8. hange of variables.
More informationPhase portraits in two dimensions
Phase portraits in two dimensions 8.3, Spring, 999 It [ is convenient to represent the solutions to an autonomous system x = f( x) (where x x = ) by means of a phase portrait. The x, y plane is called
More informationarxiv:physics/ v1 [physics.ed-ph] 24 May 2006
arxiv:physics/6521v1 [physics.ed-ph] 24 May 26 Comparing a current-carrying circular wire with polygons of equal perimeter: Magnetic field versus magnetic flux J P Silva and A J Silvestre Instituto Superior
More informationSample Solutions of Assignment 9 for MAT3270B
Sample Solutions of Assignment 9 for MAT370B. For the following ODEs, find the eigenvalues and eigenvectors, and classify the critical point 0,0 type and determine whether it is stable, asymptotically
More informationME 680- Spring Geometrical Analysis of 1-D Dynamical Systems
ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions
More informationAnalysis of the Takens-Bogdanov bifurcation on m parameterized vector fields
Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationMath 215/255 Final Exam (Dec 2005)
Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.
More informationQUARTERLY OF APPLIED MATHEMATICS
QUARTERLY OF APPLIED MATHEMATICS Volume LIV December 1996 Number 4 DECEMBER 1996, PAGES 601-607 ON EXISTENCE OF PERIODIC ORBITS FOR THE FITZHUGH NERVE SYSTEM By S. A. TRESKOV and E. P. VOLOKITIN Institute
More informationThere is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)
82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationDef. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)
Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November
More informationChapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationThe Big Picture. Python environment? Discuss Examples of unpredictability. Chaos, Scientific American (1986)
The Big Picture Python environment? Discuss Examples of unpredictability Email homework to me: chaos@cse.ucdavis.edu Chaos, Scientific American (1986) Odds, Stanislaw Lem, The New Yorker (1974) 1 Nonlinear
More informationOutline. Learning Objectives. References. Lecture 2: Second-order Systems
Outline Lecture 2: Second-order Systems! Techniques based on linear systems analysis! Phase-plane analysis! Example: Neanderthal / Early man competition! Hartman-Grobman theorem -- validity of linearizations!
More informationMath 308 Final Exam Practice Problems
Math 308 Final Exam Practice Problems This review should not be used as your sole source for preparation for the exam You should also re-work all examples given in lecture and all suggested homework problems
More information