ENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)
|
|
- Moris Cannon
- 5 years ago
- Views:
Transcription
1 ENGI Linear Approximation () Page 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 critical point at (0, 0) may be expressed as x = ax + by + P ( x, y) y = cx + d y + Q x, y () where a, b, c, d are all constants and P( x, y) Q( x, y) lim = 0 and lim ( xy, ) ( 0,0 ) ( xy, ) ( 0,0) x + y x + y = 0. Near the critical point (0, 0), this system may be approximated by the linear system x = ax + by y = cx + d y () Effect of Small Perturbations Small perturbations in the values of the coefficients a, b, c, d are reflected in small changes in the eigenvalues λ. If the eigenvalues are a pure imaginary pair, λ = ± jv, then the critical point is a centre. The effect of small changes in the coefficients will change the eigenvalues to the complex conjugate pair λ' = u' ± jv', where u' is small in magnitude and v' is close to v. The trajectories of the new system are likely to be spirals. The critical point will be an asymptotically stable focus if u' < 0, a stable (but not asymptotically stable) centre if u' = 0, but it will be an unstable focus if u' > 0. Therefore small perturbations in a linear system with pure imaginary eigenvalues are likely to result in radical changes in the trajectories (orbits) and may change the stable system into an unstable system.
2 ENGI Linear Approximation () Page 4. If the eigenvalues are a real equal pair, λ = λ, then a slight perturbation is likely to separate the roots into distinct values. If those values are still real, then the critical point remains a node. If the perturbed eigenvalues are a complex conjugate pair, then the nature of the trajectories will change into spirals and the critical point changes from a node into a focus. However, in both cases, an asymptotically stable critical point remains asymptotically stable after a small perturbation, while an unstable critical point remains unstable. In all other cases, a slight perturbation leaves the sign of the real part of both eigenvalues unchanged and affects neither the type of critical point nor the overall type of the orbits. These results are summarized in the two following theorems. Poincaré s Theorem: The singularities of the non-linear system () are identical to the singularities of the linear system (), except for the cases D < 0 and a + d = 0, which is a centre in the linear case, but may be a centre or a focus in the non-linear case; and D = 0, which is either a node or a focus in the non-linear case. Theorem on stability of the singularity at (0, 0): Linear approximation asymptotically stable unstable stable but not asymptotically stable Non-linear system asymptotically stable unstable any (unstable, stable or asymptotically stable)
3 ENGI Linear Approximation () Page 4.4 Example Perform a stability analysis on the system dx dy = x x + xy, = y xy 6y Find the critical points: dx dy x( x+ y) = 0 and y( x 6y) = 0 This generates four solutions: x = 0 and y = 0 x, y = 0,0 = = 0 x x + xy = 0 and y xy y = x = 0 and x 6y = 0 x, y = 0, x+ y = 0 and y = 0 x, y =,0 8 x+ y = 0 and x 6y = 0 x, y =, Linearize the system near each critical point Near the critical point (0, 0), it is obvious that the linear approximation to the system is dx dy = x, = y This is an uncoupled system of ODEs, whose general solution is quickly found (by direct integration or by the results on page 4.0) to be t t x t, y t = ce, c e The critical point (0, 0) is an unstable node. Also note that y = k x, so that the solutions are parabolas all sharing the same vertex at the origin, except for the lines x = 0 (from c = 0) and y = 0 (from c = 0). The orbits all begin at the origin and move away, (a hallmark of an unstable node). ( () ()) ( ) This diagram is valid for the linear approximation everywhere, but is valid only in the immediate neighbourhood of (0, 0) for the non-linear system.
4 ENGI Linear Approximation () Page 4.5 Example (continued) P x, y = x x + xy, Q x, y = y xy 6 y P P = x + y = x x y Q Q and = y = x y x y Near a critical point (a, b), the linear approximation to the system is P P x x ( ab, ) y ( ab, ) x a a+ b a x a = = y Q Q y b b a b y b x ( ab, ) y ( ab, ) Near the critical point ( 0, ), the linear approximation to the system is 4 x + 0 x 0 x = = y y y xz, = xy, leads to A simple change of variables to 4 0 x x = z z The matrix is triangular, allowing us to identify the eigenvalues immediately. The eigenvalues are of opposite sign (+4/ and ). Therefore the critical point at ( 0, ) is an unstable saddle point. From page 4.0, 4 4 ( 0 ) 0 D = a d + bc = + + = The eigenvector corresponding to eigenvalue λ = is 0 0 α a d D 0 = any non-zero multiple of β = =, c The eigenvector corresponding to eigenvalue λ = +4/ is 0 0 α + a d + D 0 = any non-zero multiple of β = = c
5 ENGI Linear Approximation () Page 4.6 Example (continued) The asymptotes are therefore x = 0 (corresponding to c = 0; inward because λ < 0 ) and x = 0( y ), (corresponding to c = 0; outward because λ > 0 ). This is sufficient to sketch the phase portrait, even without the general solution. The general solution near the critical point ( 0, ) is 4 4 ( (), ()) = ( (), () ) = 0, + t t t x t z t x t y t c e c e c e where c, c are arbitrary constants. All trajectories (except for c = 0) therefore come in from infinity near the asymptote x = 0 (where c = 0) and all trajectories (except for c = 0) return to infinity near the asymptote y = 0 x (where c = 0). The diagram is valid for the linear approximation everywhere, but is valid only in the immediate neighbourhood of ( 0, ) for the non-linear system.
6 ENGI Linear Approximation () Page 4.7 Example (continued) Near the critical point (, 0), the linear approximation to the system is x x x = = y 0 y 0 y Change variables to ( wy, ) = ( x, y). Then w w = y 0 y The matrix is triangular, with eigenvalues of opposite sign ( and +). Therefore the critical point at (, 0) is also a[n unstable] saddle point. D = a d + 4bc = + 0 = 4 a d D α = any non-zero multiple of β = =, 0 c 0 c = 0, so the alternative form is needed for the other eigenvector: b α = any non-zero multiple of β ( d a) + D = + = The general solution near the critical point (, 0) is ( () ()) () () t t t ( ) w t, y t = x t, y t = c e + c e, c e where c, c are arbitrary constants. All trajectories (except for c = 0) therefore come in from infinity near the asymptote y = 0 (where c = 0) and all trajectories (except for c = 0) return to infinity near the asymptote y = (x ) (where c = 0). Again, the diagram is valid for the linear approximation everywhere, but is valid only in the immediate neighbourhood of (, 0) for the non-linear system.
7 ENGI Linear Approximation () Page 4.8 Example (continued) Near the critical point ( 8 7, 7 ), the linear approximation to the system is x x x = = y y y Change variables to ( wz, ) = ( x 8, y ) w 8 8 w = z 7 6 z D = ( a d) + 4bc = + 4 = = D < 0 and (a + d) < 0 the critical point is an asymptotically stable focus λ 7 The eigenvalues, which are the solutions to det = 0, are λ ( a+ d) ± D 4 ± 8 7 λ = = = ± j Using the formula on page 4.0, u =, v= =, u d = + = ( wt, zt ) xt, yt = = 7 7 t t t t t t t e c cos sin c4 cos sin, ccos c sin Redefining the arbitrary constants as c5 = c and c6 = c 4, 7 7 the general solution near the critical point 8 ( 7 7) ( wt (), zt ()) = 8 ( xt (), yt () ) = 7 7, is t t t t t t e c5 cos + 7 sin + c6 7 cos + sin, c5cos + c6 sin where c 5, c 6 are arbitrary constants. t 7 All trajectories spiral in to the critical point [a phase portrait for the linear approximation is on the next page].
8 ENGI Linear Approximation () Page 4.9 Example (continued) Maple produces the following direction field plots: One can clearly see trajectories flowing outward from the unstable node at the origin in all directions. The natures of the other critical points are somewhat less obvious. Zooming in to the neighbourhood of one of the saddle points and to the neighbourhood of the stable focus produces the next pair of diagrams:
9 ENGI Linear Approximation () Page 4.40 Example (continued) Maple can also superimpose some trajectories on these phase portraits: Maple commands for this plot: > with(detools): > phaseportrait( [diff(x(t),t) = x(t) - x(t)^ + x(t)*y(t), diff(y(t),t) = *y(t) - x(t)*y(t) - 6*y(t)^], [x(t), y(t)], t=-0..0, [[x(0)=0., y(0)=0.0], [x(0)=0.8, y(0)=-0.0], [x(0)=., y(0)=0.0], [x(0)=0.05, y(0)=0.], [x(0)=0.05, y(0)=0.4], [x(0)=., y(0)=-0.0]], x=-0...4, y= , stepsize=.0, colour=red, linecolour=[blue, cyan, magenta, sienna, orange, black], title=`example Non-Linear Solution`);
10 ENGI Linear Approximation () Page 4.4 Example Perform a stability analysis on the system dx dy = y + x x y, = x + y x y Find the critical points: dx dy Clearly (x, y) = (0, 0) satisfies = = 0. dx dy = = x = 0 y = 0 dx dy = = y = 0 x= 0 If x 0 and y 0, then at any critical point ( ) ( ) 0 ( ) and ( ) y + x x y = x + y x y = y = x x y x = y x y x y y ( x y ) = = = y x x which has no solution for real (x, y). Therefore the only critical point is (0, 0). The linear approximation to the non-linear system is x x = y y D = a d + 4bc = + 4 = 4 < 0 (a + d) = > 0 the critical point is an unstable focus. Using the formula on page 4.0, a+ d a d D u = =, u d = = 0, v = = x t, y t = e t c sin t c cos t, c cost + c sin t ( () ()) ( 4 4 ) and c, c 4 are [real] arbitrary constants.
11 ENGI Linear Approximation () Page 4.4 Example (continued) Solution in the neighbourhood of the only critical point (0, 0): Now consider the distance r of any point (x, y) from the critical point (0, 0): dr dx dy r = x + y r = x + y [chain rule] From the original non-linear system: dr r = x( y + x( x y )) + y( x + y( x y )) = xy + x r xy + y r = r r ( r ) ( r ) dr < 0 > = r( r ) > 0 < Therefore solutions starting closer than one unit to the critical point spiral out, but solutions starting further away than one unit approach the critical point. A solution on the circle r = never changes its distance from the origin and stays on that circle, but is not stationary (because the only critical point is at (0, 0), not on that circle). Note that x + y = is a solution to the non-linear equation: dx dy = y + x( x y ) = y and = x + y( x y ) = x dy dy dx x = = dx y dy dy x But x + y = x + y = 0 = dx dx y No other solution may cross this solution. Therefore all non-trivial solutions approach the unit circle in the limit as t. The unit circle here is an example of a limit cycle.
12 ENGI Linear Approximation () Page 4.4 Example (continued) Further notes: Let u = r, then du = r dr and r dr = r ( r ) du = u( u) du = u u ( ) + du = ln u ln ( u) = t + C u u u u t + C t ln = t + C = e = ce u u u ce t uce t ce t t = + u = ce t ce + u =, where c = c u = r() t = t + ce + ce t If r(0) = r o, then t u ro ro ce = c = c = u ro ro Therefore r() t = r o + e ro t
13 ENGI Linear Approximation () Page 4.44 Example (continued) Also note that, from the polar coordinate system, ( xy, ) = ( rcos θ, r sinθ ) xy, = r cosθ rsin θθ, r sinθ + rcosθθ ( θ ) θ ( θ) θ θ ( θ) θ ( θ) xy yx = rcos r sin + rcos rcos rsin r cos + rsin rsinθ θ ( cos sin ) = r θ + θ θ = r θ But the non-linear system is x, y = y + x r, x + y r ( ) ( ) xy yx = x + xy r y xy r = x + y = r Therefore = = = + C r θ r θ θ t All paths move clockwise with constant angular speed. In Cartesian coordinates all orbits can therefore be described by x() t, y() t = cos t θo, sin t θo r o t + e ro ( ) where (x(0), y(0)) = r o (cos θ o, sin θ o )
14 ENGI Linear Approximation () Page 4.45 Example (A more challenging and tedious case, for reference only) Find and determine the nature of all critical points of the system dx x dy = e y, = y sinx () dx dy x = = 0 y = e and y = sinx x Critical points occur where the graphs of y = e and y = sin x intersect. The critical points are (0, 0) and (x i, y i ), where ( 4 π ) δ, ( 4 π ) ε ; (,,, ) x = i i + i = i i i Linearize for the critical point at (0, 0): near ( 0,0) ( 0,0) ( 0,0) 0 x P P e y x + y = e x + y x y Q Q y sin x x + y = ( cos 0) x + y x y near ( 0,0) ( 0,0) ( 0,0) Therefore the linear system that models the non-linear system () near (0, 0) is x x = y y Finding the eigenvalues: λ det ( A λi) = = ( + λ)( λ) ( )( ) = λ λ ( A I) det λ = 0 λ = λ = ± () The eigenvalues are real and of opposite sign. The critical point (0, 0) of () [and therefore also of ()] is an unstable saddle point.
15 ENGI Linear Approximation () Page 4.46 Example (continued) Using the results on page 4.0, D = ( a d) + 4bc = ( ) + 4( )( ) = 4+ 4 = 8 The eigenvectors are α ( a d) D = any non-zero multiple of β = c for λ = and α ( a d) + D + = any non-zero multiple of β = c for λ = + Choose a multiple of in both cases. The general solution of () is therefore () () () () ( t t t ) ( () ()) t xt, yt = c + e + c e, ce + ce y t lim = = = > 0 and t x t + y t + lim = = = < 0 t + x t The trajectories (except for c = 0) therefore come in from infinity along the asymptote y = x. The trajectories (except for c = 0) return to infinity along the asymptote y = + x. The phase space diagram for this solution (completely valid only for the linear system ()) is The phase space for the non-linear system () resembles this diagram only in the immediate neighbourhood of the critical point (0, 0).
16 ENGI Linear Approximation () Page 4.47 Example (continued) At other critical points (k, l), dx dy = = 0 l = sink e sink = 0 Linearizing (Taylor s series for P(x, y) about (x, y) = (k, l): x P P ( e y ) ( x k) + ( y l) = ( e )( x k) + ( )( y l near ( kl, ) x y near ( kl, ) ( kl) ( kl) ( kl),, Q Q y sin x ( x k) + y l = cos k x k + x y ( kl),, ( y l) Therefore the linear system that models the non-linear system () near (k, l) is x e x k = y cos k y l Finding the eigenvalues: e λ det ( A λi) = = ( e + λ)( λ) ( cos k)( ) cos k λ ( e k ) λ ( e cosk) 0 ( ) ( = λ + = λ = e ± e + 4 e + cosk ( e ) ( e ) ± + + 4cosk λ = Now k > 0 0< e < 0< e < Recall that k = xi = ( 4i π ) δi, ( 4i π ) + εi ; ( i =,,, ) Examining the right-hand critical point in each pair, k = ( 4i π ) + ε, i ( 0< εi ) cos k = cos 4i π + ε = cos ε π = cos π ε = sin ε ε > 0 ( ) Therefore e 4( e cosk) ( e ) ) i i i i i + + > > 0 and the two eigenvalues are real and of opposite sign. These critical points are therefore all unstable saddle points. ) ()
17 ENGI Linear Approximation () Page 4.48 Example (continued) Examining the left-hand critical point in each pair, k = 4i π δ, 0< δ i ( i ) ( δi) ( δi ) ( δi ) cos k = cos 4i π = cos π = cos + π = sinδ k But k is the solution to e sin k = 0 ( ) e = + sin k = + sin 4i π δi = + sin δ i π = sin δ π i + = cosδ i δ i e + cos k ( cosδi) sinδi δi < 0 ( k e ) 4( e cosk) ( + + < e e ± e + 4 e + cosk) But λ = Therefore the eigenvalues are a real distinct positive pair and the singularity is an unstable node. The locations and nature of the first five critical points are listed here. k (δ i is small) x y λ λ Type unstable saddle point unstable node unstable saddle point unstable node unstable saddle point Here is a Maple session to create direction field plots for the first three critical points of the non-linear system. > with(detools): > DEplot([diff(x(t),t) = -y(t) - + exp(-x(t)), diff(y(t),t) = y(t) - sin(x(t))], [x(t),y(t)], t=-.., x= , y= , title=`example Exact Solution`); > DEplot([diff(x(t),t) = -y(t) - + exp(-x(t)), diff(y(t),t) = y(t) - sin(x(t))], [x(t),y(t)], t=-.., x=4.5..5, y= , title=`example Exact Solution`); i
18 ENGI Linear Approximation () Page 4.49 Example (continued) Saddle Point at (0, 0) Unstable Node near (4.57, 0.99), Saddle Point near (4.84, 0.99) Solution curves can be traced by following the arrows (which at every location point in the direction of dy dx ).
ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.0 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential equations
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 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 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 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 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 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 informationMath 3301 Homework Set Points ( ) ( ) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, ( ) ( ) ( ) ( )
#7. ( pts) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, λ 5 λ 7 t t ce The general solution is then : 5 7 c c c x( 0) c c 9 9 c+ c c t 5t 7 e + e A sketch of
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 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 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 informationAutonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.
Autonomous equations Autonomous systems Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. i f i(x 1, x 2,..., x n ) for i 1,..., n As you
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 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 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 informationSection 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d
Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d July 6, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This
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 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 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 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 informationMATH 44041/64041 Applied Dynamical Systems, 2018
MATH 4441/6441 Applied Dynamical Systems, 218 Answers to Assessed Coursework 1 1 Method based on series expansion of the matrix exponential Let the coefficient matrix of α the linear system be A, then
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 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 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 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 informationAP Calculus (BC) Chapter 10 Test No Calculator Section. Name: Date: Period:
AP Calculus (BC) Chapter 10 Test No Calculator Section Name: Date: Period: Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The graph in the xy-plane represented
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
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
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 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 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 Chapter 9. u. (c) u(t) = 1 e t + c 2 e 3 t! c 1 e t 3c 2 e 3 t. (v) (a) u(t) = c 1 e t cos 3t + c 2 e t sin 3t. (b) du
Solutions hapter 9 dode 9 asic Solution Techniques 9 hoose one or more of the following differential equations, and then: (a) Solve the equation directly (b) Write down its phase plane equivalent, and
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 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 informationCreated by T. Madas LINE INTEGRALS. Created by T. Madas
LINE INTEGRALS LINE INTEGRALS IN 2 DIMENSIONAL CARTESIAN COORDINATES Question 1 Evaluate the integral ( x + 2y) dx, C where C is the path along the curve with equation y 2 = x + 1, from ( ) 0,1 to ( )
More informationFirst Midterm Exam Name: Practice Problems September 19, x = ax + sin x.
Math 54 Treibergs First Midterm Exam Name: Practice Problems September 9, 24 Consider the family of differential equations for the parameter a: (a Sketch the phase line when a x ax + sin x (b Use the graphs
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 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 informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
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 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 informationProject 6.3 Predator-Prey and Your Own Game Preserve
Project 6.3 Predator-Prey and Your Own Game Preserve The closed trajectories in the figure below represent periodic solutions of a typical predator-prey system, but provide no information as to the actual
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 informationODE, part 2. Dynamical systems, differential equations
ODE, part 2 Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Dynamical systems, differential equations Consider a system of n first order equations du dt = f(u, t),
More informationEven-Numbered Homework Solutions
-6 Even-Numbered Homework Solutions Suppose that the matric B has λ = + 5i as an eigenvalue with eigenvector Y 0 = solution to dy = BY Using Euler s formula, we can write the complex-valued solution Y
More informationDynamical 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 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 informationC3 papers June 2007 to 2008
physicsandmathstutor.com June 007 C3 papers June 007 to 008 1. Find the exact solutions to the equations (a) ln x + ln 3 = ln 6, (b) e x + 3e x = 4. *N6109A04* physicsandmathstutor.com June 007 x + 3 9+
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
More informationExam 1 Review SOLUTIONS
1. True or False (and give a short reason): Exam 1 Review SOLUTIONS (a) If the parametric curve x = f(t), y = g(t) satisfies g (1) = 0, then it has a horizontal tangent line when t = 1. FALSE: To make
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 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 informationENGI 9420 Lecture Notes 2 - Matrix Algebra Page Matrix operations can render the solution of a linear system much more efficient.
ENGI 940 Lecture Notes - Matrix Algebra Page.0. Matrix Algebra A linear system of m equations in n unknowns, a x + a x + + a x b (where the a ij and i n n a x + a x + + a x b n n a x + a x + + a x b m
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 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 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 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 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 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 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 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 informationMath 1272 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Exam ) This fraction appears in Problem 5 of the undated-? exam; a solution can be found in that solution set. (E) ) This integral appears in Problem 6 of the Fall 4 exam;
More informationLyapunov functions and stability problems
Lyapunov functions and stability problems Gunnar Söderbacka, Workshop Ghana, 29.5-10.5, 2013 1 Introduction In these notes we explain the power of Lyapunov functions in determining stability of equilibria
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 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 information1. Find and classify the extrema of h(x, y) = sin(x) sin(y) sin(x + y) on the square[0, π] [0, π]. (Keep in mind there is a boundary to check out).
. Find and classify the extrema of hx, y sinx siny sinx + y on the square[, π] [, π]. Keep in mind there is a boundary to check out. Solution: h x cos x sin y sinx + y + sin x sin y cosx + y h y sin x
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 3 2, 5 2 C) - 5 2
Test Review (chap 0) Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Solve the problem. ) Find the point on the curve x = sin t, y = cos t, -
More informationC3 A Booster Course. Workbook. 1. a) Sketch, on the same set of axis the graphs of y = x and y = 2x 3. (3) b) Hence, or otherwise, solve the equation
C3 A Booster Course Workbook 1. a) Sketch, on the same set of axis the graphs of y = x and y = 2x 3. b) Hence, or otherwise, solve the equation x = 2x 3 (3) (4) BlueStar Mathematics Workshops (2011) 1
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 informationENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations
ENGI 344 - Second Order Linear ODEs age -01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of
More informationMath 216 Final Exam 24 April, 2017
Math 216 Final Exam 24 April, 2017 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 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 informationSolved Examples. Given are two sets A {1, 2, -2, 3} and B = {1, 2, 3, 5}. Is the function f(x) = 2x - 1 defined from A to B?
Solved Examples Example 1: Given are two sets A {1, 2, -2, 3} and B = {1, 2, 3, 5}. Is the function f(x) = 2x - 1 defined from A to B? Solution : Out of all the ordered pairs, the ordered pairs which are
More informationFind the general solution of the system y = Ay, where
Math Homework # March, 9..3. Find the general solution of the system y = Ay, where 5 Answer: The matrix A has characteristic polynomial p(λ = λ + 7λ + = λ + 3(λ +. Hence the eigenvalues are λ = 3and λ
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 information10.1 Review of Parametric Equations
10.1 Review of Parametric Equations Recall that often, instead of representing a curve using just x and y (called a Cartesian equation), it is more convenient to define x and y using parametric equations
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 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 informationParametric Curves. Calculus 2 Lia Vas
Calculus Lia Vas Parametric Curves In the past, we mostly worked with curves in the form y = f(x). However, this format does not encompass all the curves one encounters in applications. For example, consider
More information36. Double Integration over Non-Rectangular Regions of Type II
36. Double Integration over Non-Rectangular Regions of Type II When establishing the bounds of a double integral, visualize an arrow initially in the positive x direction or the positive y direction. A
More informationENGI Parametric Vector Functions Page 5-01
ENGI 3425 5. Parametric Vector Functions Page 5-01 5. Parametric Vector Functions Contents: 5.1 Arc Length (Cartesian parametric and plane polar) 5.2 Surfaces of Revolution 5.3 Area under a Parametric
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 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 informationPhase Portraits of Nonlinear Differential Equations
ODE4-net.nb Phase Portraits of Nonlinear Differential Equations Nonlinear Differential Equations: x' = f (x, y) (1) Consider the system where f and g are functions of two y' = g(x, y) () variables x and
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 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 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 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 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 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 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 informationHomogeneous Constant Matrix Systems, Part II
4 Homogeneous Constant Matri Systems, Part II Let us now epand our discussions begun in the previous chapter, and consider homogeneous constant matri systems whose matrices either have comple eigenvalues
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 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 informationVector Field Topology. Ronald Peikert SciVis Vector Field Topology 8-1
Vector Field Topology Ronald Peikert SciVis 2007 - Vector Field Topology 8-1 Vector fields as ODEs What are conditions for existence and uniqueness of streamlines? For the initial value problem i x ( t)
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 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 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 informationSolutions to Final Exam Sample Problems, Math 246, Spring 2018
Solutions to Final Exam Sample Problems, Math 46, Spring 08 () Consider the differential equation dy dt = (9 y )y. (a) Find all of its stationary points and classify their stability. (b) Sketch its phase-line
More informationSECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes.
SECTION A 1. State the maximal domain and range of the function f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. 2. By evaluating f(0),
More information