Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
|
|
- Horatio Webb
- 5 years ago
- Views:
Transcription
1 Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the phase diagram associated to We just need to find the solutions to Xy which are y ± In order to classify the fixed points we have to either look at the sign of the rhs of in the various regions between fixed points or use the criterium coming from linearisation which means looking at the sign of the derivative of the rhs of at each of the fixed points Following this last criterium we find that So we find that: X y y 4 yy 3y + 4 X 4 > therefore y is a repellor X 8 > therefore y is an attractor X 8 < therefore y is also an attractor Hence the phase diagram is b Solve for each of the following initial conditions: y 3 y indicating the range of values of t for which each solution is defined The equation can be solved by separating variables in the usual way: t y dy y dt t t t y yy y + y 8 + y 4 y y 8 log y 4y y 4y with y yt the initial condition This can be solved to [ mark] yt ± e 4t t e 8t t + 4 y The ± sign will be + for y > and for y < in order to satisfy the initial condition Let us now solve for each initial condition: for y 3 we have yt e 4t e 8t 9 [ marks] This solution is well defined only when e 8t 9 > This gives the condition t > 8 log Thus the range of values for which the solution is defined is t 8 log 9
2 For y the solution becomes [4 marks] e4t yt e 8t + 3 where we had to introduce a minus sign to guarantee that the initial condition is fulfilled In this case the denominator is never zero and e 8t + 3 can never become negative so the solution is valid for all values of t The interval of definition is t [4 marks] c Provide a sketch of your solutions against the variable t clearly showing the behaviour of the various functions as t ± In order to carry out the sketch it is important to study the behaviour of both solutions as t ± The solution yt e 4t y 3 e 8t 9 is clearly always positive for those values of t for which it is defined and as t it tends to the value which is the fixed point we identified before This was to be expected as the point is an attractor In this case it does not make sense to consider t as the solution is only defined for t 8 log 9 Finally since the denominator vanishes exactly at t 8 log 9 the function should diverge as it gets closer to that value Figure : The red figure corresponds to the initial condition y 3 As we can see it diverges at a small negative value of t which was computed above The blue graph is the solution with y The green lines are the fixed points y ± Concerning the other solution e4t yt y e 8t + 3 it is clearly a negative function for all values of t As t it tends to the value - which is the fixed point identified before We should have expected that since - is an attractor As t the solution tends to the value which is also compatible with our fixed point analysis as zero was identified as a repellor A plot of the two solutions is given above d Write down the linearized version of equation about a generic fixed point y a At a fixed point Xa and therefore a Taylor expansion of Xy to first order in y gives dy dt X ay a 3a + 4y a [ marks] Page
3 [ mark] e Use the result of d to find the solution of the linearized equation about the smallest of the fixed points of if the initial condition is y 3 We already found in a that X 8 therefore the linearized equation is dz dt 8z for z y + The solution to this equation is z Ae 8t where A is a constant which will be fixed by the initial condition The final answer is: yt e 8t [Seen] Consider the following system of linear differential equations: dx dt where α β γ c c are real constants ẋ αx + γx + c dx dt ẋ βx + αx + c a Write the equations in the standard form ẋ Ax + b Write down the vector a which is a fixed point of ẋ α γ x c + β α c The fixed point is the solution to the vector equation that is ẋ Ax + b x A α γ b α βγ β α x αc γc c α βγ βc c +αc α βγ [ mark] [4 marks] b Compute the eigenvalues and eigenvectors of the matrix A The eigenvalues of A are obtained by finding the zeroes of the characteristic polynomial of A: the solutions to this equation are A λi α λ βγ α βγ + λ αλ and λ α + 4α 4α βγ λ α 4α 4α βγ α + βγ α βγ The eigenvectors corresponding to these eigenvalues are obtained by solving: α γ a α ± a βγ β α b b In components this gives two equations for each eigenvector: aα + bγ α ± βγa aβ + bα α ± βγb Page 3
4 simplifying bγ ± βγa aβ ± βγb The two equations are not independent and setting a this gives b ± β eigenvectors would be: E β γ E c Define the vector z in terms of which the equation can be re-written in the form ż Az We just need to change variables to z x a where a is the fixed point obtained in section a d Indicate the conditions if any that α β and γ must satisfy so that the fixed point is: an improper node a focus a star node β γ γ Therefore the [ marks] [ mark] An improper node arises when the eigenvalues of the matrix A are equal and real In our case the eigenvalues can only be equal if βγ which corresponds to β or γ [ mark] A focus occurs when the eigenvalues are complex conjugated to each other This will happen only if β > and γ < or viceversa [ mark] Finally a star node occurs when the matrix A is proportional to the identity matrix This will only happen if β γ αc γc α βγ βc +αc α βγ [ mark] e If α β and γ classify the fixed point Find the values of c and c so that the fixed point is a We found the fixed point in section a c c a c +c f This is solved by c 3 and c For α β and γ show that the matrix A is already in its Jordan Normal form For these values of α β γ and the values of c c obtained in e find the general solution of equations For the values of α β γ given in e the eigenvalues of A are λ + i and λ i The Jordan form is therefore: which coincides with the matrix A itself J [ mark] Page 4
5 We will use the change of variables of section c so that the equation becomes of the form ż Az where and z x a a from section e The equation in z can be written in components as: ż z + z ż z + z and can be solved by using the standard change of variables z r cos θ z r sin θ This leads to the equations where A is a constant and ṙ r r Ae t θ θ t + B where B is an integration constant Therefore the general solutions are z Ae t cos t + B z Ae t sin t + B Therefore for the original variables x x the general solutions are x + Ae t cos t + B x + Ae t sin t + B or in vector form x + Ae t cos t + B + Ae t sin t + B [ marks] 3 [Seen] Consider the following system of first order nonlinear differential equations: ẋ x x x ẋ x x x 3 a Show that the fixed points of the system are either at the origin or on a line In order to find the fixed points we have to solve the simultaneous equations: x x x x x x which give solutions x x and x x which is indeed a line b Calculate the Jacobian matrix for a general point x x and deduce that the fixed point at the origin is simple whereas the fixed points on the line are not The Jacobian matrix is obtained as usual by computing the derivatives A xx X x X X x X x x x x x x x x with X x x x x x and X x x x x x Page
6 Evaluating this matrix at the origin we obtain A which is the identity matrix The determinant of this matrix is obviously and therefore nonvanishing So the origin is a simple fixed point For any point α α lying on the line x x we find that the Jacobian matrix is α α A α α + α + α which has determinant deta therefore the line of fixed points x x contains only non-simple fixed points c By writing down the linearization about the fixed point at the origin classify its nature Solve the linearized equation We have already almost answered this question above when we computed the Jacobian matrix at the origin The linearized equation about the origin is: x A x x which x and since A x is the identity matrix the fixed point at the origin is an unstable star node The equation is very easy to solve In components we have that which is solved by ẋ x ẋ x x C e t x C e t with C C arbitrary constants Since x /x C /C phase space trajectories near the fixed point will be straight lines going by the origin d By dividing the two equations in 3 obtain a new equation of the form dx dx Xx x Find general solution to this equation Explain why your solution is what one would expect given the nature of the fixed point at the origin The resulting equation is dx x x x dx x x x x x which only is well defined when x x and x The equation can be solved by separation of variables dx dx x x which gives or ln x ln x + C x e C x where C is an integration constant Therefore in general phase space trajectories are straight lines going by the origin Page 6
7 This is what one expects from the fact that the origin is an unstable star node which is precisely characterized by trajectories which are straight lines with various slopes depending on the initial condition e Use the results of c and d to sketch the phase space diagram associated to 3 Pay special attention to the behaviour of phase space trajectories near the line of fixed points found in a Does the direction of the trajectories change at this line? Why? The phase space diagram is given in the figure below Figure : The phase space diagram of equation 3 The blue line is the line of fixed points x x The phase space trajectories change direction exactly at the line x x this is because the sign of the derivatives ẋ and ẋ changes at this line Below the line that is for x < x we have that x x > therefore ẋ x x x > ẋ x x x > if both x x > This is compatible with trajectories moving away from the origin in the direction of increasing x and x However for x > x that is above the line of fixed points we have exactly the opposite behaviour that is ẋ x x x < ẋ x x x < if both x x > Therefore here the trajectories have to move towards the origin at time grows that is in the direction of decreasing values of x and x 4 [Unseen] Consider the following second order linear differential equation: where ẍ d x dt and ẋ dx dt ẍ 4ẋ x 4 a Consider now the following two first order linear differential equations ẋ ax + bx + c ẋ dx + ex + f Page 7
8 where x x above Find the values of the constants a b c d e and f so that the system is equivalent to 4 By defining x x and x ẋ ẋ we can rewrite 4 as: ẋ x ẋ 4x + x so a c f and b d e 4 b Using the result of a write equation 4 in the standard matrix form ẋ Ax + b ẋ x ẋ 4 x with b in this case [ mark] c Find the eigenvalues and eigenvectors of the matrix A Hence construct the Jordan normal form of A and the matrix P which relates A to its Jordan form J as A P JP The eigenvalues of A are obtained by finding the zeroes of the characteristic polynomial of A: A λi λ4 λ λ λ The eigenvectors corresponding to these eigenvalues are obtained by solving: a a 4 b b In components this gives two equations: b a a + 4b b which are not independent Taking a we find that the eigenvector associated to λ is E and for λ we obtain the equation 4 In components this gives two equations: a b a b b a a + 4b b which are not independent Taking a we find that the second eigenvector is E The Jordan normal form of A is just the diagonal matrix J and the matrix P is built just by setting its columns to the eigenvectors above: P Page 8
9 d Find the fixed point of the system of equations and classify its nature In this case the fixed point is just at the origin [ mark] Because the eigenvalues are real and one eigenvalue is negative and the other positive this would be a saddle e Find the general solution of the system of equations ẏ Jy Hence find the general solution for x In this case the canonical system of equations is simply: ẏ y ẏ y which is solved by y C e t and y C e t where C and C are arbitrary constants Now we need to take the solutions from part e and multiply them by the matrix P to obtain x: y y + y x P y C e t + C e t x y y + y C e t + C e t x f This can be rewritten as x C e t E + C e t E Sketch the phase diagram both in the y y plane and in the x x plane Figure 3: The y y diagram is a typical saddle As expected the x x diagram would be a rotated and slightly stretched version of the first diagram where the symmetry axes are now the directions determined by the eigenvectors E E [Unseen] Consider the following first order nonlinear differential equations [4 marks] ẋ x x x ẋ + x x x 6 which can be interpreted as modelling the populations of two competing species Page 9
10 a Find the four fixed points of the system To find the fixed points we have to solve the simultaneous equations: x x x + x x x An obvious solution is x x For x we can also have x and for x we can also have x Finally if both x and x are not zero there is also the solution x x 6 is also a fixed point Therefore we have the four points: 6 b Find the general form of the Jacobian matrix of the system and evaluate it at each fixed point Classify the fixed points according to their linearization The Jacobian matrix is obtained by computing the derivatives of the rhs of the equations 6 wrt x and x Defining X x x x / x x and X x x + x / x x we obtain A X x X x X X x x At the various fixed points it becomes: A x x x x x + x which is diagonal with real eigenvalues of different signs This means that the origin is a saddle A [ mark] In this case we need to find the eigenvalues of A in order to classify the fixed point They are the solution to the equation λ λ that is λ and λ Therefore as before this would also be a saddle A 3 The eigenvalues of A are the solutions to the equation 3 λ λ that is λ and λ 3 Therefore both eigenvalues are real and positive which corresponds to an unstable node Finally 3 A 6// 6 The eigenvalues of A 6// are the solutions to the equation 3 λ λ + λ + 4λ + 3 that is λ ± i These are complex conjugated eigenvalues with negative real part which would describe a stable focus Page
11 c For only those fixed points lying in the region x and x solve the linearized equations The only fixed points lying on the indicated region are and 6/ / 4 The first fixed point is a saddle and the linearized equation is already in the canonical form Therefore the general solution to this linearized equation is: where C C are arbitrary constants x C e t + C e t [ mark] The fixed point 6/ / 4 is a stable focus In order to solve the linearized equation about the fixed point we need to carry out the standard transformations First we change variables to z x a where a is the fixed point so that we bring the fixed point to the origin Second we rewrite the resulting linearized equation ż Az in the canonical form and solve this equation: ẏ Jy where J is the Jordan form of A and z P y where P is a nonsingular matrix that relates A and J as A P JP [ mark] In order to construct the matrix P we need the eigenvectors of A The eigenvalues λ ± i mean that the Jordan form of the matrix A is J The eigenvectors of A are obtained by solving the equations 3 6 a a b ± i b which can be written as 3a 6b ± i a a b ± i b As usual the two equations for each eigenvalue are not independent b 6 ± i 6 therefore the eigenvectors and eigenvalues are and E E 6 i i 6 λ + i Taking a we obtain λ i The matrix P that relates A to its Jacobian form is obtained by taking the real and imaginary parts of the first eigenvector and setting them as columns of J P 6 6 and we can check that indeed A P JP Page
12 [ mark] The equation ẏ Jy can be solved in the usual way for J of the form above see section g of question The solution takes the form y C e t cos t + C + C e t sin t + C d where C C are constants In order to get the solutions for the variable z we need to multiply the vector above by the matrix P and finally to get x we have to add the fixed point x a + P y 6 6 In components: + C e t cos t + C P + C e t cos t + C 6 + C e t sin t + C P + C e t sin t + C x 6 + C e t cos t + C x C t e cos 6 t + C C e t sin 6 t + C 6 In the same region of the x x plane considered in section c and using the solutions to the linearized equations sketch the phase diagram associated to the nonlinear system of equations You may employ the vector field diagram below as a guiding tool in order to produce your sketch Figure 4: The sketch here is not entirely accurate as I have ignored the presence of the other fixed points at - and Specially - would affect the shape of the trajectories on the lower right hand corner of the picture as one can see from the vector field diagram The students were not expected to notice this Page
Dynamical Systems Solutions to Exercises
Dynamical Systems Part 5-6 Dr G Bowtell Dynamical Systems Solutions to Exercises. Figure : Phase diagrams for i, ii and iii respectively. Only fixed point is at the origin since the equations are linear
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 informationMath 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016
Math 4B Notes Written by Victoria Kala vtkala@math.ucsb.edu SH 6432u Office Hours: T 2:45 :45pm Last updated 7/24/206 Classification of Differential Equations The order of a differential equation is the
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 informationSystems of Linear ODEs
P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here
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 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 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 informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationCalculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y
Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationNonlinear differential equations - phase plane analysis
Nonlinear differential equations - phase plane analysis We consider the general first order differential equation for y(x Revision Q(x, y f(x, y dx P (x, y. ( Curves in the (x, y-plane which satisfy this
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 informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationWeek Quadratic forms. Principal axes theorem. Text reference: this material corresponds to parts of sections 5.5, 8.2,
Math 051 W008 Margo Kondratieva Week 10-11 Quadratic forms Principal axes theorem Text reference: this material corresponds to parts of sections 55, 8, 83 89 Section 41 Motivation and introduction Consider
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 informationSolutions of Spring 2008 Final Exam
Solutions of Spring 008 Final Exam 1. (a) The isocline for slope 0 is the pair of straight lines y = ±x. The direction field along these lines is flat. The isocline for slope is the hyperbola on the left
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
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 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 Higgins-Selkov oscillator
The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being
More 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 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 informationPractice Final Exam Solutions for Calculus II, Math 1502, December 5, 2013
Practice Final Exam Solutions for Calculus II, Math 5, December 5, 3 Name: Section: Name of TA: This test is to be taken without calculators and notes of any sorts. The allowed time is hours and 5 minutes.
More informationEngg. Math. I. Unit-I. Differential Calculus
Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle
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 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 informationBifurcation Analysis of Non-linear Differential Equations
Bifurcation Analysis of Non-linear Differential Equations Caitlin McCann 0064570 Supervisor: Dr. Vasiev September 01 - May 013 Contents 1 Introduction 3 Definitions 4 3 Ordinary Differential Equations
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 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 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 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 informationModelling and Mathematical Methods in Process and Chemical Engineering
Modelling and Mathematical Methods in Process and Chemical Engineering Solution Series 3 1. Population dynamics: Gendercide The system admits two steady states The Jacobi matrix is ẋ = (1 p)xy k 1 x ẏ
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationMath 113 Winter 2005 Key
Name Student Number Section Number Instructor Math Winter 005 Key Departmental Final Exam Instructions: The time limit is hours. Problem consists of short answer questions. Problems through are multiple
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
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 information2t t dt.. So the distance is (t2 +6) 3/2
Math 8, Solutions to Review for the Final Exam Question : The distance is 5 t t + dt To work that out, integrate by parts with u t +, so that t dt du The integral is t t + dt u du u 3/ (t +) 3/ So the
More information+ i. cos(t) + 2 sin(t) + c 2.
MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,
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 informationVectors, matrices, eigenvalues and eigenvectors
Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a
More informationMathematics 104 Fall Term 2006 Solutions to Final Exam. sin(ln t) dt = e x sin(x) dx.
Mathematics 14 Fall Term 26 Solutions to Final Exam 1. Evaluate sin(ln t) dt. Solution. We first make the substitution t = e x, for which dt = e x. This gives sin(ln t) dt = e x sin(x). To evaluate the
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 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 informationState will have dimension 5. One possible choice is given by y and its derivatives up to y (4)
A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To
More informationMath 216 First Midterm 18 October, 2018
Math 16 First Midterm 18 October, 018 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 informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationMA 527 first midterm review problems Hopefully final version as of October 2nd
MA 57 first midterm review problems Hopefully final version as of October nd The first midterm will be on Wednesday, October 4th, from 8 to 9 pm, in MTHW 0. It will cover all the material from the classes
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 informationDynamical Systems. August 13, 2013
Dynamical Systems Joshua Wilde, revised by Isabel Tecu, Takeshi Suzuki and María José Boccardi August 13, 2013 Dynamical Systems are systems, described by one or more equations, that evolve over time.
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
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 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 informationALGEBRAIC GEOMETRY HOMEWORK 3
ALGEBRAIC GEOMETRY HOMEWORK 3 (1) Consider the curve Y 2 = X 2 (X + 1). (a) Sketch the curve. (b) Determine the singular point P on C. (c) For all lines through P, determine the intersection multiplicity
More informationSolutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011
Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 8, Solutions to Problems 5 are omitted since they involve topics no longer covered on the Comprehensive
More informationLinear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form
Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form d x x 1 = A, where A = dt y y a11 a 12 a 21 a 22 Here the entries of the coefficient matrix
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 informationTHE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK. Summer Examination 2009.
OLLSCOIL NA héireann, CORCAIGH THE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK Summer Examination 2009 First Engineering MA008 Calculus and Linear Algebra
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 informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline
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 informationMathematics 1EM/1ES/1FM/1FS Notes, weeks 18-23
2 MATRICES Mathematics EM/ES/FM/FS Notes, weeks 8-2 Carl Dettmann, version May 2, 22 2 Matrices 2 Basic concepts See: AJ Sadler, DWS Thorning, Understanding Pure Mathematics, pp 59ff In mathematics, a
More informationPractice Final Exam Solutions
Important Notice: To prepare for the final exam, study past exams and practice exams, and homeworks, quizzes, and worksheets, not just this practice final. A topic not being on the practice final does
More informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationFINAL EXAM SOLUTIONS, MATH 123
FINAL EXAM SOLUTIONS, MATH 23. Find the eigenvalues of the matrix ( 9 4 3 ) So λ = or 6. = λ 9 4 3 λ = ( λ)( 3 λ) + 36 = λ 2 7λ + 6 = (λ 6)(λ ) 2. Compute the matrix inverse: ( ) 3 3 = 3 4 ( 4/3 ) 3. Let
More informationMath 142, Final Exam. 12/7/10.
Math 4, Final Exam. /7/0. No notes, calculator, or text. There are 00 points total. Partial credit may be given. Write your full name in the upper right corner of page. Number the pages in the upper right
More informationReview Problems for Exam 2
Review Problems for Exam 2 This is a list of problems to help you review the material which will be covered in the final. Go over the problem carefully. Keep in mind that I am going to put some problems
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 informationE209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions
E9A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions Michael Vitus Gabe Hoffmann Stanford University Winter 7 Problem 1 The governing equations are: ẋ 1 = x 1 + x 1 x ẋ = x + x 3 Using
More informationMath 310 Introduction to Ordinary Differential Equations Final Examination August 9, Instructor: John Stockie
Make sure this exam has 15 pages. Math 310 Introduction to Ordinary Differential Equations inal Examination August 9, 2006 Instructor: John Stockie Name: (Please Print) Student Number: Special Instructions
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 informationHomogeneous Constant Matrix Systems, Part II
4 Homogeneous Constant Matrix Systems, Part II Let us now expand our discussions begun in the previous chapter, and consider homogeneous constant matrix systems whose matrices either have complex eigenvalues
More informationMath 5BI: Problem Set 6 Gradient dynamical systems
Math 5BI: Problem Set 6 Gradient dynamical systems April 25, 2007 Recall that if f(x) = f(x 1, x 2,..., x n ) is a smooth function of n variables, the gradient of f is the vector field f(x) = ( f)(x 1,
More informationExam Two. Phu Vu. test Two. Take home group test April 13 ~ April 18. Your group alias: Your group members: Student name
Exam Two Take home group test April 3 ~ April 8 Your group alias: Your group members: (leave it blank if you work alone on this test) Your test score Problem Score Total page /7 Problem : (chapter 9, applications
More informationMathQuest: Differential Equations
MathQuest: Differential Equations Geometry of Systems 1. The differential equation d Y dt = A Y has two straight line solutions corresponding to [ ] [ ] 1 1 eigenvectors v 1 = and v 2 2 = that are shown
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 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 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 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 informationUNIVERSITY OF SOUTHAMPTON. A foreign language dictionary (paper version) is permitted provided it contains no notes, additions or annotations.
UNIVERSITY OF SOUTHAMPTON MATH055W SEMESTER EXAMINATION 03/4 MATHEMATICS FOR ELECTRONIC & ELECTRICAL ENGINEERING Duration: 0 min Solutions Only University approved calculators may be used. A foreign language
More informationPractice Final Exam Solutions
Important Notice: To prepare for the final exam, one should study the past exams and practice midterms (and homeworks, quizzes, and worksheets), not just this practice final. A topic not being on the practice
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 informationExam Basics. midterm. 1 There will be 9 questions. 2 The first 3 are on pre-midterm material. 3 The next 1 is a mix of old and new material.
Exam Basics 1 There will be 9 questions. 2 The first 3 are on pre-midterm material. 3 The next 1 is a mix of old and new material. 4 The last 5 questions will be on new material since the midterm. 5 60
More informationStudy Guide/Practice Exam 3
Study Guide/Practice Exam 3 This study guide/practice exam covers only the material since exam. The final exam, however, is cumulative so you should be sure to thoroughly study earlier material. The distribution
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 informationB. Differential Equations A differential equation is an equation of the form
B Differential Equations A differential equation is an equation of the form ( n) F[ t; x(, xʹ (, x ʹ ʹ (, x ( ; α] = 0 dx d x ( n) d x where x ʹ ( =, x ʹ ʹ ( =,, x ( = n A differential equation describes
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 informationAutonomous Systems and Stability
LECTURE 8 Autonomous Systems and Stability An autonomous system is a system of ordinary differential equations of the form 1 1 ( 1 ) 2 2 ( 1 ). ( 1 ) or, in vector notation, x 0 F (x) That is to say, an
More informationMath Exam III - Spring
Math 3 - Exam III - Spring 8 This exam contains 5 multiple choice questions and hand graded questions. The multiple choice questions are worth 5 points each and the hand graded questions are worth a total
More information4 Partial Differentiation
4 Partial Differentiation Many equations in engineering, physics and mathematics tie together more than two variables. For example Ohm s Law (V = IR) and the equation for an ideal gas, PV = nrt, which
More information6.3. Nonlinear Systems of Equations
G. NAGY ODE November,.. Nonlinear Systems of Equations Section Objective(s): Part One: Two-Dimensional Nonlinear Systems. ritical Points and Linearization. The Hartman-Grobman Theorem. Part Two: ompeting
More informationYou may use a calculator, but you must show all your work in order to receive credit.
Math 2410-010/015 Exam II April 7 th, 2017 Name: Instructions: Key Answer each question to the best of your ability. All answers must be written clearly. Be sure to erase or cross out any work that you
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
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 informationFirst and Second Order ODEs
Civil Engineering 2 Mathematics Autumn 211 M. Ottobre First and Second Order ODEs Warning: all the handouts that I will provide during the course are in no way exhaustive, they are just short recaps. Notation
More informationMath 113 Winter 2005 Departmental Final Exam
Name Student Number Section Number Instructor Math Winter 2005 Departmental Final Exam Instructions: The time limit is hours. Problem consists of short answer questions. Problems 2 through are multiple
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 informationFurther Mathematical Methods (Linear Algebra)
Further Mathematical Methods (Linear Algebra) Solutions For The 2 Examination Question (a) For a non-empty subset W of V to be a subspace of V we require that for all vectors x y W and all scalars α R:
More information