Half of Final Exam Name: Practice Problems October 28, 2014
|
|
- Adelia Underwood
- 5 years ago
- Views:
Transcription
1 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 will be comprehensive.. Consider this enhanced predator prey system for the populations of two species x, y. Find the nullclines and use this information to sketch the phase space. Indicate the equilibrium points and the direction of flow. Determine the stability types of the rest points and find the stable and unstable spaces of the linearized flows at these points. Sketch enough trajectories to describe the global flow. x = x( + y) y = y( x y) x = if x = or y =. y = if y = or x + y =. Thus both equations are satisfied at the equilibrium points P = (, ), P 2 = (, ) and P = (2, ). We note that x > if y < and x < if y >. Similarly y > if x + y < and y < if x + y >. Thus in the four regions NE, NW, SW and SE cut by y = and x + y = going clockwise around P, the direction of flow is SE, NE, NW ande SW, respectively. The Jacobean is At P we have DF (x, y) = y y DF (, ) = x. 2y x ( ) so that the eigenvalues and eigenvectors are λ =, V = (, ); λ 2 =, V 2 = (, ). This is a saddle with E s and E u the horizontal and vertical axes, resp. At P 2 we have DF (, ) = ( ) 2 so that the eigenvalues and eigenvectors are λ = 2, V = (5, ); λ 2 =, V 2 = (, ). This is a saddle with E s the vertical axis and E u = span{v }. At P we have DF (2, ) = ( ) 2. The trace is negative and the determinant is positive, so that the eigenvalues are both negative or have negative real parts. Indeed the eigenvalues are λ = 2 ( + 7i). This is a stable spiral. All of this information is given in Fig.. The computed phase plane generated by D- XplorMath is given in Fig. 2.
2 Figure : Nullclines, Equilibrium Points, Directions and Stable/Unstable Spaces. Figure 2: Trajectories of the predator prey system of Problem.. 2
3 2. Consider the folowing model of reaction diffusion equation system. Determine if it is a gradient or a Hamiltonian system. Use this fact to sketch the global flow. Identify the stability type of the equilibrium points. x = 2(y x) + x( x 2 ) y = 2(y x) + y( y 2 ) We see that f y (x, y) = 2 and g x (x, y) = 2 are equal so there is a potential function F (x, y) = 2 x2 + 4 x4 2xy + 2 y2 + 4 x4 and the system is a gradient system ( ) x = F (x, y). y The critical points satisfy = f(x, y) = F x = 2y x x = g(x, y) = F y = 2x y y Thus the equilibrium points are the intersections of the curves y = 2 (x + x ) and y = 2 (y + y ), which are the points P = (, ), P 2 = (, ) and P = (, ). Computing the Hessian, ( ) ( ) Fxx F xy x 2 H(x, y) = = 2 y 2 At P, F yx F yy H(, ) = ( ) 2 2 whose eigenvalues are λ = 2 and λ 2 = 6, hence P 2 and P are relative minima. At P 2 and P, H(, ) = H(, ) = ( ) whose eigenvalues are λ = and λ 2 =, hence is a saddle. It follows that P is a saddle for the flow and both P 2 and P are stable nodes. The flow is perpendicular to the level surfaces. The D-XplorMath plot of level sets and trajectories is given in Fig. 4.. Suppose that f(x, t) C ( R 2, R ) and that (, x ) R 2. Assume that f satisfies a local Lipschitz Condition: for every compact set K R 2 there is L R so that f(t, x) f(t, y) L x y whenever (t, x), (t, y) K. Then there is an ɛ > and a function x(t) C ([, + ɛ), R 2 ) which solves the initial value problem x = f(t, x) () x( ) = x. This is the short time existence theorem for non-autonomous equations. Note that we only assume continuity with respect to t so that we cannot consider time as an independent variable in a system of one higher dimension and quote the result we proved in class for autonomous equations. However, we may follow the proof we gave verbatim.
4 Figure : Flow of gradient system and level curves of the potential. For any r > let R = [, + r] [x r, x + r] be a compact compact box and L its corresponding Lipschitz constant. Let { Let ɛ = min r, M } r,. Let L + M = sup f(t, x). (t,x) R X = {z(t) C([, + ɛ]) : z(t) x r whenever t [, + ɛ]} The space X is a closed convex subset of the complete metric space C([, + ɛ]) under the sup norm. The sup norm for w C([, + ɛ]) is given by w = sup w(s). s [,+ɛ] Thus the subspace (X, ) is also a complete metric space. We will show that there is a unique function z X that solves the integral equation x(t) = x + f(s, x(s)) ds for all t [, + ɛ] (2) Define the nonlinear operator T : X X for u X and t [, + ɛ] by the formula (T u)(t) = x + f(s, u(s)) ds. 4
5 Note that T u is continuous since it is the indefinite integral of a continuous function s f(s, u(s)). To see that is satisfies the bound, we have for any for u X and t [, + ɛ] (T u)(t) x = f(s, u(s)) ds f(s, u(s)) ds Thus we have shown that T : X X. M ds = M(t ) Mɛ M r M r. Now we shall show that the operator is a contraction, namely, for some θ (, ) we have for all u, v X T u T v θ u v. () In this situation, we define the Picard Iterates recursively: x (t) = x and x n+ (t) = (T x n )(t) for n N. If the mapping T is a contraction, then the Picard Iterates {x n (t)} are dominated by a geometric sequence, thus form a Uniformly Cauchy sequence on [, + ɛ], thus converges uniformly to a unique continuous function in the complete metric space (X, ) satisfying x = T x, which is the integral equation (2). This fact from analysis is called the Contraction Mapping Theorem or the Banach Fixed Point Theorem. For u, v X and t [, + ɛ] we have by the Lipschitz bound, (T u)(t) (T v)(t) = f(s, u(s)) f(s, v(s)) ds f(s, u(s)) f(s, v(s)) ds L u(s) v(s) ds L u v ds = (t ) u v L u v. L + Taking supremum over t [, + ɛ] gives () with θ = L L +. It remains to argue that the solution of the integral equation (2) also solves the initial value problem (). If x(t) solves (), then by the Fundamental Theorem of Calculus, for t [, + ɛ], x(t) x = dx t (s) ds = f(s, x(s)) ds dt which is the integral equation (2). On the other hand, if x(t) X solves the integral equation (2), then it is differentiable since by the Fundamental Theorem of Calculus, x(x) is the integral of a continuous function. Moreover, its derivative is given by the integrand Also, at the initial point dx (t) = f(t, x(t)). dt so the initial condition holds as well. x( ) = x + f(s, x(s)) ds = x +, 5
6 4. Suppose that f(x, t) C ( R 2, R ), (, x ) R 2 and a >. Assume that f satisfies a Lipschitz Condition: there is L R so that f(t, x) f(t, y) L x y whenever (t, x), (t, y) [, + a] R. Show that the soliution of the initial value problem () is unique. Suppose that both u(t) and v(t) satisfy the initial value problem (). Using ther integral equation, we find u(t) v(t) = f(s, u(s)) f(s, v(s)) ds f(s, u(s)) f(s, v(s)) ds L u(s) v(s)) ds Thus if we define the function F (t) = L u(s) v(s)) ds we find fromt the Fundamental Theorem of Calculus that F (t) = u(t) v(t)) LF (t). It follows that ( e L(t t) F ) = e L(t ) (F LF ). Thus the positive function satisfies F (t) F ( ) = for all t. Thus F (t) = or u(t) = v(t) for all t. 5. Generalize the method of Picard Approximation to accommodate non-autonomous equtions. Find the first few iterates of the initial value problem and compare to the Maclaurin series for the solution. x = 2tx x() =. The generalization is to use the non-autonomous integral equation. anything in X. Then recursively devine for n N, In this example, let x (t) =. Then x (t) = x x 2 (t) = x x (t) = x x 4 (t) = x 2sx (s) ds = 2 2sx (s) ds = 2 2sx 2 (s) ds = 2 2sx (s) ds = 2 x n+ (t) = x + s ds = t 2 ; f(s, x(s)) ds. s ( s 2) ds = t t4 ; Put x (t) = x or s ( s ) s4 ds = t t4 6 t6 ; s ( s s4 6 ) s6 ds = t t4 + 6 t6 24 t8 ; 6
7 The Maclaurin series solution is found by assuming the solutin is a convergent power series and plugging into the ODE to find equations for the coefficients. Assuming we get and x (t) = x(t) = + a n t n i= na n t n = a + 2a 2 t + a t 2 + 4a 4 t 4 + n= 2tx = 2t which results in the equations 2a n t n+ = 2t 2a t 2 2a 2 t 2a t 4 n= a =, a =, a 2 =, na n = 2a n 2, for n. This results in a 2k = for all odd terms and a 2k = 2 2k a 2(k ) whose solution is a 4 = 2, a 6 = 6, a 8 = 24 and so on. In general, for k N we have a 2k = ( )k. k! In other words x(t) = e t2. 6. Suppose that both f(t, x), g(x, t) C ( R 2, R ), (, x ) R 2 and a > satisfies the Lipschitz Condition: there is L R so that f(t, x) f(t, y) L x y whenever (t, x), (t, y) [, + a] R. Show that the soliution of the initial value problem () depends continuously on the equation. That is, assume that for some η > the two functions satisfy f(t, x) g(t, x) η for all (t, x) [, + a] R. Suppose we look at two solutions of the different equations { x = f(t, x), x( ) = x. { y = g(t, y) y( ) = x. Then x(t) y(t) η L ( ) e L(t t) for all t [, + a]. 7
8 As usual, we will try to apply Gronwall s Inequality to the difference of solutions. For t [, + a], x(t) y(t) = x + f(s, x(s)) ds x g(s, y(s)) ds = f(s, x(s)) g(s, y(s)) ds f(s, x(s)) g(s, y(s)) ds f(s, x(s)) g(s, x(s)) + g(s, x(s)) g(s, y(s)) ds f(s, x(s)) g(s, x(s)) + g(s, x(s)) g(s, y(s)) ds η + L x(s) y(s) ds = η(t ) + L x(s) y(s) ds We apply Gronwall s Inequality (Problem 5) with f(t) = η(t ), u = x(t) y(t) and v = L. For t [, + a], x(t) y(t) η(t ) + ηl e L(t s) (s ) ds = η ( ) e L(t t). L 7. Prove Gronwall s Inequality. If f, u, v : [, + a] R are continuous, u, v and satisfy u(t) f(t) + u(s) v(s) ds for all t then Let Then ( ) u(t) f(t) + exp v(u) du f(s)v(s) ds. s F (t) = u(s) v(s) ds. F = uv v(f + F ) so ( e ) t v(s) t ds v(s) F = e t ds (F vf ) e v(s) t ds vf which implies Multiplying by the integrating factor yields e t t v(r) t dr F (t) e s v(r) t dr v(s) f(s) ds. u(t) f(t) + e s v(r) dr v(s) f(s) ds. 8
9 8. Analyze the stability properties of each rest point. At each rest point, determine the dimensions of the stable and unstable manifolds. Compute a basis for the stable and unstable manifolds for the associated linearized system. We recast the equation as a system. x 7x + x + e x = x = y y = z z = ( ex ) y + 7 z The equilibrium point satisfies y =, z = so = e x so x =. The Jacobean at (,, ) is A = DF () = 7 The characteristic polynomial is ( ) 7 p(λ) = ( λ)( λ) λ + ( λ) = λ + 7 λ2 λ ( = (λ ) λ 2 4 λ ) The eigenvalues are λ =, λ 2 = and λ = 2 7. Thus the matrix is hyperbolic with two positive and one negative eigenvalue. The origin is of saddle type, thus unstable. It has a local one dimensional stable manifold and a two dimensional local unstable manifold. We solve for eigenvectors. = (A λ I)V = = (A λ 2 I)V 2 = = (A λ I)V = It follows that V and V 2 span the unstable space E u of the linearied equation y = Ay and V spans the stable space E s. 9
10 9. Find the stable and unstable manifolds. x = x y = y + x 2 z = z + y 2 The only equilibrium point is the origin. There the linearized system is A = DF () = diag(,, ) so that the stable and unstable spaces are E s = {(x, y, z) : z = } and E u = {(x, y, z) : x = y = }. To find the stable manifold, we seek solutions that are tangent to E s at the origin. The system may be solved by back substitution. Indeed, x = x means x(t) = x e t. Hence y = y + x 2 = y + x 2 e 2t so y(t) = ( y + x 2 ) e t x 2 e 2t Also z = z + y 2 = z + ( y + x) 2 2 ( e 2t 2x 2 y + x) 2 e t + x 4 e 4t so ( z(t) = z + ) ( y + x 2 2 ( ) 2 x2 y + x 2 ) + 5 x4 e t ( y + x 2 ) 2 e 2t + ( 2 x2 y + x 2 ) e t 5 x4 e 4t The stable manifold tends to zero as t so the e t coefficient must vanish, namely, points of W s satisfy z = ( y + x 2 2 ( ) + 2 x2 y + x 2 ) 5 x4 = y2 6 x2 y x4. The points in the unstable manifold W u tend to zero as t, which means that the coefficients of the e t and e 2t, e t and e 4t must vanish, namely x = y =.. Classify the stability types of all the equilibrium points of the system. x = 4y + 2xy 8 y = 4y 2 x 2 At the equilibrium, the second equation = y says x = ±2y. Using x = 2y in = x says = 2y +2y 2 4 or y =, 2, giving two rest points P = ( 2, ) and P 2 = (4, 2). Using x = 2y says y 2 + y + 2 = which has no real roots. The Jacobean is At P, ( ) 2y 4 + 2x A = DF = 2x 8y A = ( ) whose eigenvalues are λ = 5 ± i 2i, thus by the linearized stability, P is a stable node node.
11 At P 2, A = ( 4 4 ) 8 6 whose eigenvalues are λ = 2,, thus by the linearized stability, P 2 is an unstable node.. Show that there is a nontrivial periodic solution. x = y y = x + ( x 2 2y 2 )y Observe that the origin is the only rest point. Consider the function Then along the flow, in polar coordinates, V (x, y) = 2 x2 + 2 y2. V (x, y) = ( x 2 2y 2 )y 2 Thus inside the circle x 2 +y 2 2 we have x2 2y 2 2(x 2 +y 2 ) so that V and the flow is not inward. Outside the circle x 2 + y 2 we have x 2 2y 2 x 2 y 2 so that V and the flow is not outward. It follows that the annulus 2 r2 is a closed, forward invariant set without rest points. Hence it must contain its limit cycles. Indeed, for any point in the annulus, its omega limit set ω(x) must be a nonempty, closed limit set in the annulus which contains no equilibrium point since the origin is not in the annulus. By the Poincaré Bendixson Theorem, ω(x) is a closed orbit. 2. Consider the prameterized system. Find the bifurcation values µ = µ and describe the nature of the bifurcation at the µ. x = µx y + x x 2 + y 2 y = x + µy + y x 2 + y 2 Write the system in polar coordinates. With x = r cos θ and y = r sin θ we get so r cos θ r sin(θ)θ = µr cos θ r sin θ + r 2 cos θ r sin θ + r cos(θ)θ = r cos θ + µr sin θ + r 2 sin θ θ = r = µr + r 2 If µ then r > all r > so that the origin is an unstable node. If µ < then at the radius r = µ, there is a limit cycle C. In this case if < r < µ we have r < and for r > µ, r >. In other words, the cycle C is unstable: nearby trajectories depart from C as t. At the same time, the origin is a stable rest point. Thus this system undergoes a Hopf bifurcation. As µ passes through zero from negative to positive, an unstable limit cycle C and a stable rest point zero collapse to form an unstable rest point zero.
12 . Determine the stability type of the zero solution. x + (x ) + x = The corresponding system is x = y y = x y Consider the Lyapunov function V (x, y) = 2 x2 + 2 y2. We have V = xx + yy = y 4. Since v, the origin is stable. In fact, in any closed ball P = B(, δ), the zero set is Z = {(x, y) : V (x, y) = } = {(x, y) : y = }. In this set there are no invariant subsets other than {(, )} because, if y = and x then y so the trajectory leaves Z. By LaSalle s Principle, (, ) is asymptotically stable and P is in the basin of its attraction. But as δ > is arbitrary, every point (x, y ) R 2 is attracted to the origin. 2
Problem 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 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 informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
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 informationMain topics and some repetition exercises for the course MMG511/MVE161 ODE and mathematical modeling in year Main topics in the course:
Main topics and some repetition exercises for the course MMG5/MVE6 ODE and mathematical modeling in year 04. Main topics in the course:. Banach fixed point principle. Picard- Lindelöf theorem. Lipschitz
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
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 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 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 informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
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 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 informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More informationProblem List MATH 5173 Spring, 2014
Problem List MATH 5173 Spring, 2014 The notation p/n means the problem with number n on page p of Perko. 1. 5/3 [Due Wednesday, January 15] 2. 6/5 and describe the relationship of the phase portraits [Due
More informationChapter 9 Global Nonlinear Techniques
Chapter 9 Global Nonlinear Techniques Consider nonlinear dynamical system 0 Nullcline X 0 = F (X) = B @ f 1 (X) f 2 (X). f n (X) x j nullcline = fx : f j (X) = 0g equilibrium solutions = intersection of
More informationSummary of topics relevant for the final. p. 1
Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
More informationApplied Dynamical Systems
Applied Dynamical Systems Recommended Reading: (1) Morris W. Hirsch, Stephen Smale, and Robert L. Devaney. Differential equations, dynamical systems, and an introduction to chaos. Elsevier/Academic Press,
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationSolutions to Final Exam Sample Problems, Math 246, Spring 2011
Solutions to Final Exam Sample Problems, Math 246, Spring 2 () Consider the differential equation dy dt = (9 y2 )y 2 (a) Identify its equilibrium (stationary) points and classify their stability (b) Sketch
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationNonlinear Autonomous Dynamical systems of two dimensions. Part A
Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More informationMA5510 Ordinary Differential Equation, Fall, 2014
MA551 Ordinary Differential Equation, Fall, 214 9/3/214 Linear systems of ordinary differential equations: Consider the first order linear equation for some constant a. The solution is given by such that
More informationMATH 32A: MIDTERM 2 REVIEW. sin 2 u du z(t) = sin 2 t + cos 2 2
MATH 3A: MIDTERM REVIEW JOE HUGHES 1. Curvature 1. Consider the curve r(t) = x(t), y(t), z(t), where x(t) = t Find the curvature κ(t). 0 cos(u) sin(u) du y(t) = Solution: The formula for curvature is t
More informationEcon 204 Differential Equations. 1 Existence and Uniqueness of Solutions
Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More 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 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 informationMath 234 Final Exam (with answers) Spring 2017
Math 234 Final Exam (with answers) pring 217 1. onsider the points A = (1, 2, 3), B = (1, 2, 2), and = (2, 1, 4). (a) [6 points] Find the area of the triangle formed by A, B, and. olution: One way to solve
More information2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.
Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following
More informationd(x n, x) d(x n, x nk ) + d(x nk, x) where we chose any fixed k > N
Problem 1. Let f : A R R have the property that for every x A, there exists ɛ > 0 such that f(t) > ɛ if t (x ɛ, x + ɛ) A. If the set A is compact, prove there exists c > 0 such that f(x) > c for all x
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
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 informationAPPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.
APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
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 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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More 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 informationMathematical Analysis Outline. William G. Faris
Mathematical Analysis Outline William G. Faris January 8, 2007 2 Chapter 1 Metric spaces and continuous maps 1.1 Metric spaces A metric space is a set X together with a real distance function (x, x ) d(x,
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 informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More 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 informationME 680- Spring Representation and Stability Concepts
ME 680- Spring 014 Representation and Stability Concepts 1 3. Representation and stability concepts 3.1 Continuous time systems: Consider systems of the form x F(x), x n (1) where F : U Vis a mapping U,V
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationEconomics 204 Fall 2011 Problem Set 2 Suggested Solutions
Economics 24 Fall 211 Problem Set 2 Suggested Solutions 1. Determine whether the following sets are open, closed, both or neither under the topology induced by the usual metric. (Hint: think about limit
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
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 informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
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 informationNonlinear Dynamical Systems Lecture - 01
Nonlinear Dynamical Systems Lecture - 01 Alexandre Nolasco de Carvalho August 08, 2017 Presentation Course contents Aims and purpose of the course Bibliography Motivation To explain what is a dynamical
More informationMath 10C - Fall Final Exam
Math 1C - Fall 217 - Final Exam Problem 1. Consider the function f(x, y) = 1 x 2 (y 1) 2. (i) Draw the level curve through the point P (1, 2). Find the gradient of f at the point P and draw the gradient
More informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More informationMAT 211 Final Exam. Spring Jennings. Show your work!
MAT 211 Final Exam. pring 215. Jennings. how your work! Hessian D = f xx f yy (f xy ) 2 (for optimization). Polar coordinates x = r cos(θ), y = r sin(θ), da = r dr dθ. ylindrical coordinates x = r cos(θ),
More informationMath 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv
Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationStable and Unstable Manifolds for Planar Dynamical Systems
Stable and Unstable Manifolds for Planar Dynamical Systems H.L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 184 October 7, 211 Key words and phrases: Saddle point, stable and
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 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 informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationMATH 32A: MIDTERM 1 REVIEW. 1. Vectors. v v = 1 22
MATH 3A: MIDTERM 1 REVIEW JOE HUGHES 1. Let v = 3,, 3. a. Find e v. Solution: v = 9 + 4 + 9 =, so 1. Vectors e v = 1 v v = 1 3,, 3 b. Find the vectors parallel to v which lie on the sphere of radius two
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 informationVolumes of Solids of Revolution Lecture #6 a
Volumes of Solids of Revolution Lecture #6 a Sphereoid Parabaloid Hyperboloid Whateveroid Volumes Calculating 3-D Space an Object Occupies Take a cross-sectional slice. Compute the area of the slice. Multiply
More informationExistence and uniqueness of solutions for nonlinear ODEs
Chapter 4 Existence and uniqueness of solutions for nonlinear ODEs In this chapter we consider the existence and uniqueness of solutions for the initial value problem for general nonlinear ODEs. Recall
More informationThird In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix
Third In-Class Exam Solutions Math 26, Professor David Levermore Thursday, December 2009 ) [6] Given that 2 is an eigenvalue of the matrix A 2, 0 find all the eigenvectors of A associated with 2. Solution.
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationNow I switch to nonlinear systems. In this chapter the main object of study will be
Chapter 4 Stability 4.1 Autonomous systems Now I switch to nonlinear systems. In this chapter the main object of study will be ẋ = f(x), x(t) X R k, f : X R k, (4.1) where f is supposed to be locally Lipschitz
More informationSolutions of Math 53 Midterm Exam I
Solutions of Math 53 Midterm Exam I Problem 1: (1) [8 points] Draw a direction field for the given differential equation y 0 = t + y. (2) [8 points] Based on the direction field, determine the behavior
More informationDo not write below here. Question Score Question Score Question Score
MATH-2240 Friday, May 4, 2012, FINAL EXAMINATION 8:00AM-12:00NOON Your Instructor: Your Name: 1. Do not open this exam until you are told to do so. 2. This exam has 30 problems and 18 pages including this
More informationMath 216 Final Exam 14 December, 2012
Math 216 Final Exam 14 December, 2012 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationFIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material
More informationMATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation
MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering one-dimensional systems, we saw that subtle changes in parameter
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1
ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate
More informationMAS331: Metric Spaces Problems on Chapter 1
MAS331: Metric Spaces Problems on Chapter 1 1. In R 3, find d 1 ((3, 1, 4), (2, 7, 1)), d 2 ((3, 1, 4), (2, 7, 1)) and d ((3, 1, 4), (2, 7, 1)). 2. In R 4, show that d 1 ((4, 4, 4, 6), (0, 0, 0, 0)) =
More informationLecture Notes for Math 524
Lecture Notes for Math 524 Dr Michael Y Li October 19, 2009 These notes are based on the lecture notes of Professor James S Muldowney, the books of Hale, Copple, Coddington and Levinson, and Perko They
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationDYNAMICAL SYSTEMS AND NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS
DYNAMICAL SYSTEMS AND NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS J. Banasiak School of Mathematical Sciences University of KwaZulu-Natal, Durban, South Africa 2 Contents 1 Solvability of ordinary differential
More informationDepartment of Mathematics IIT Guwahati
Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,
More informationSection 5. Graphing Systems
Section 5. Graphing Systems 5A. The Phase Plane 5A-1. Find the critical points of each of the following non-linear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A-2. Write each
More information= e t sin 2t. s 2 2s + 5 (s 1) Solution: Using the derivative of LT formula we have
Math 090 Midterm Exam Spring 07 S o l u t i o n s. Results of this problem will be used in other problems. Therefore do all calculations carefully and double check them. Find the inverse Laplace transform
More informationNon-Linear Dynamics Homework Solutions Week 6
Non-Linear Dynamics Homework Solutions Week 6 Chris Small March 6, 2007 Please email me at smachr09@evergreen.edu with any questions or concerns reguarding these solutions. 6.8.3 Locate annd find the index
More informationMath 118, Fall 2014 Final Exam
Math 8, Fall 4 Final Exam True or false Please circle your choice; no explanation is necessary True There is a linear transformation T such that T e ) = e and T e ) = e Solution Since T is linear, if T
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 informationMath 216 First Midterm 19 October, 2017
Math 6 First Midterm 9 October, 7 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationMath 350 Solutions for Final Exam Page 1. Problem 1. (10 points) (a) Compute the line integral. F ds C. z dx + y dy + x dz C
Math 35 Solutions for Final Exam Page Problem. ( points) (a) ompute the line integral F ds for the path c(t) = (t 2, t 3, t) with t and the vector field F (x, y, z) = xi + zj + xk. (b) ompute the line
More informationImplicit Functions, Curves and Surfaces
Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then
More informationApproximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle
Malaya J. Mat. 4(1)(216) 8-18 Approximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle B. C. Dhage a,, S. B. Dhage a and S. K. Ntouyas b,c, a Kasubai,
More informationMath 361: Homework 1 Solutions
January 3, 4 Math 36: Homework Solutions. We say that two norms and on a vector space V are equivalent or comparable if the topology they define on V are the same, i.e., for any sequence of vectors {x
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 information