DOCSLIB.ORG

TWO DIMENSIONAL FLOWS Lecture 4: Linear and Nonlinear Systems

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
TWO DIMENSIONAL FLOWS Lecture 4: Linear and Nonlinear Systems TWO DIMENSIONAL FLOWS Lecture 4: Linear and Nonlinear Systems 4. Linear and Nonlinear Systems in 2D In higher dimensions, trajectories have more room to manoeuvre, and hence a wider range of behaviour is possible. 4.1 Linear systems: definitions and examples A 2-dimensional linear system has the form x˙ = ax + by y˙ = cx + dy where a, b, c, d are parameters. Equivalently, in vector notation x˙ = Ax (1) where a b x A = and x = (2) c d ! y ! The Linear property means that if x1 and x2 are solutions, then so is c1x1 + c2x2 for any c1 and c2. The solutions ofx ˙ = Ax can be visualized as trajectories moving on the (x,y) plane, or phase plane. 1 Example 4.1.1 mx¨ + kx = 0 i.e. the simple harmonic oscillator Fig. 4.1.1 The state of the system is characterized by x and v =x ˙ x˙ = v k v˙ = x −m i.e. for each (x,v) we obtain a vector (˙x, v˙) ⇒ vector field on the phase plane. 2 As for a 1-dimensional system, we imagine a fluid flowing steadily on the phase plane with a local velocity given by (˙x, v˙) = (v, ω2x). − Fig. 4.1.2 Trajectory is found by placing an imag- • inary particle or phase point at (x0,v0) and watching how it moves. (x,v) = (0, 0) is a fixed point: • static equilibrium! Trajectories form closed orbits around (0, 0): • oscillations! 3 The phase portrait looks like... Fig. 4.1.3 NB ω2x2 + v2 is constant on each ellipse. • This is simply the energy Example 4.1.2 x˙ a 0 x = y˙ 0 1 y ! − ! ! 4 The phase portraits for these uncoupled equa- tions are... Fig. 4.1.4 Solution is at x x0e = t y ! y0e− ! 5 Some terminology... x =0 is an attracting fixed point in Figs • ∗ (a) - (c) since x(t) x as t . → ∗ →∞ x = 0 is called Lyapunov Stable in Figs • ∗ (a) - (d) since all trajectories that start sufficiently close to x∗ remain close to it for all time. Fig. (d) shows that a fixed point can • be Lyapunov stable but not attracting ⇒ it is neutrally stable. It is also possible for a fixed point to be attracting but not Lyapunov stable! If a fixed point is both Lyapunov sta- • ble and attracting, we’ll call it stable, or sometimes asymptotically stable x is unstable in Fig. (e) because it is • ∗ neither attracting nor Lyapunov stable 6 4.2 Classification of Linear Systems Consider a general 2 2 matrix A such that × x˙ = Ax To solve: try x(t) = eλtv (v is a constant vector) λeλtv = eλtAv ⇒ Av = λv ⇒ Hence if we obtain the eigenvectors v and eigenvalues λ, we will have two independent a b solutions x(t). Recall that A = has c d ! eigenvalues λ1 and λ2, where τ + τ 2 4∆ τ τ 2 4∆ λ = − λ = − − 1 q 2 2 q 2 with τ = trace(A) = a + d ∆ = det(A) = ad bc − 7 Useful check when calculating eigenval- • ues: λ1 + λ2 = τ and λ1λ2 = ∆ x˙ 1 1 x Example 4.2.1 = y˙ 4 2 y ! − ! ! 1 λ1 = 2 with v1 = λ1 > 0 • ⇒ 1 ! hence solution grows 1 λ2 = 3 with v2 = λ2 < 0 • ⇒ − 4 ! hence solution decays − Fig. 4.2.1 8 straight line trajectories in Fig. 4.2.1 are • the eigenvectors v1 and v2 Example 4.2.2 Consider λ2 < λ1 < 0 Fig. 4.2.2 Both solutions decay exponentially! • 9 Example 4.2.3 What happens if λ1, λ2 are complex? Fixed point is either... Fig. 4.2.3 If λ , λ are purely imaginary, all solutions • 1 2 are periodic If λ = λ we get a star node or a degen- • 1 2 erate node 10 Classification of Fixed Points λ = 1(τ τ 2 4∆), where 1,2 2 ± − ∆= λ1λ2 andq τ = λ1 + λ2 Fig. 4.2.4 11 4.3 Phase Portraits Recallx ˙ = f(x), i.e. x˙1 = f1(x1,x2) x˙2 = f2(x1,x2) where x = (x1,x2) and f(x) = (f1(x), f2(x)) (not necessarily linear now). The trajectories x(t) wind their way through the phase plane. The entire phase plane is filled with trajec- tories! 4.4 Example of a phase portrait - Shows a sample of the qualitatively different trajectories 12 Fig. 4.4.1 x Fixed points A, B and C satisfy f( ∗) = 0 • and correspond to steady states or equi- libria Closed orbit D corresponds to periodic • solutions, i.e. x(t + T )= x(t) for all t for some T > 0 The existence and uniqueness theorem • given for 1-dimensional systems can be generalized to 2-dimensional systems ... fortunately different trajectories never intersect! ⇒ 13 4.5 Fixed points and Linearization This is the same idea as for 1-dimensional systems x˙ = f(x,y) y˙ = g(x,y) Suppose (x∗,y∗) is a fixed point. Expand around (x ,y ) using u = x x and v = y y . ∗ ∗ − ∗ − ∗ u˙ =x ˙ = f(x∗ + u,y∗ + v) ∂f ∂f 2 2 = f(x∗,y∗)+ u + v + O(u ,v ,uv) ∂x ∂y ∂f ∂f u + v ≃ ∂x ∂y Similarly ∂g ∂g v˙ u + v ≃ ∂x ∂y Hence a small disturbance around (x∗,y∗) evolves as ∂f ∂f u˙ u = ∂x ∂y v˙ ∂g ∂g v ! ∂x ∂y ! where the matrix is known as the Jacobian A matrix at (x∗,y∗), and is the multivariable equivalent of f ′(x∗) for 1-D systems. 14 Example 4.5.1 x˙ = x + x3 − y˙ = 2y − Fixed points occur wherex ˙ = 0 andy ˙ = 0 simultaneously. Hence x = 0 or x = 1 ± and y = 0 3 fixed points (0, 0), (1, 0) and ⇒ ( 1, 0) − Jacobian matrix A ∂x˙ ∂x˙ 2 ∂x ∂y 1 + 3x 0 A = = − ∂y˙ ∂y˙ 0 2 ∂x ∂y − ! 1 0 At (0, 0) A = − stable node 0 2 ⇒ − ! 2 0 At ( 1, 0) A = both are sad- ± 0 2 ! ⇒ dle points. − 15 Fig. 4.5.1 In general, we must obtain fixed points by solvingx ˙ = 0 andy ˙ = 0 simultaneously. e.g. x˙ = x(3 x 2y) − − y˙ = y(2 x y) − − yields fixed points (0, 0), (0, 2), (3, 0) and (1, 1) A In general, will not be diagonal at (x∗,y∗). Hence we must diagonalize A, i.e. find eigen- values λ1 and λ2 and eigenvectors v1 and v2 of A 16 Basically, we are doing the same here as be- fore for 2D linear systems, since we are treat- ing the nonlinear system as linear near (x∗,y∗). Knowledge of λ1 and λ2, and v1 and v2, en- ables us to sketch the phase portrait near (x∗,y∗). The fixed points can be classified according to their stability as follows: If Re(λ ) > 0 and Re(λ ) > 0 • 1 2 repeller (unstable node) ⇒ If Re(λ ) < 0 and Re(λ ) < 0 • 1 2 attractor (stable node) ⇒ If Re(λ ) > 0 but Re(λ ) < 0 (or vice • 1 2 versa) saddle ⇒ If λ and λ are both imaginary centre • 1 2 ⇒ 17 4.6 Example: Rabbits vs Sheep An example of the Lotka-Volterra model of competition between two species (e.g. rab- bits and sheep) grazing the same food supply (grass). Each species grows to its carrying capac- • ity in the absence of the other - logistic growth (rabbits faster...!) When species encounter each other, the • larger (sheep) has an advantage. Conflicts occur at a rate proportional to • the size of each population. Conflicts re- duce the growth rate of each species (but more for rabbits). A model encapsulating these properties could be (see above!) x˙ = x(3 x 2y) − − y˙ = y(2 x y) − − 18 Fixed points at 3 0 (0, 0) where A = λ = 3, 2 0 2 ! ⇒ 1 0 (0, 2) where A = − λ = 1, 2 2 2 ⇒ − − − − ! 1 2 (1, 1) where A = − − λ = 1 √2 1 1 ⇒ − ± − − ! 3 6 (3, 0) where A = − − λ = 3, 1 0 1 ⇒ − − − ! (0, 0): λ = 3, 2 unstable node (repeller) ⇒ λ = 2 v = (0, 1) “slow eigendirection ⇒ ′′ λ = 3 v = (1, 0) “fast eigendirection ⇒ ′′ General rule... Trajectories are tangential to the slow eigendirection (i.e. smallest λ ) at a node | | 19 (0, 2): λ = 1, 2 stable node (attrac- − − ⇒ tor) Once again... Trajectories are tangential • to the slow eigendirection at a node Here λ = 1 v = (1, 2) is the slow • − ⇒ − eigendirection. (1, 1): λ = 1 √2 saddle point − ± ⇒ 20 (3, 0): λ = 3, 1 stable node (attrac- − − ⇒ tor) Putting these together, the phase portrait becomes.... Fig. 4.6.1 NB: You don’t really need to calculate the eigenvectors to get the right shape! 21 Biological interpretation... In general, one species eventually drives • the other to extinction; which species even- tually dominates depends on initial pop- ulations x0 = (x0,y0) Basin of attraction of an attracting fixed • point x∗ defined as the set of initial con- ditions x0 such that x x as t . → ∗ →∞ In this case, basin boundary is the stable • manifold of the saddle node at (1, 1) Fig. 4.6.2 22 4.7 Conservative Systems Considerx ˙ = f(x). A conserved quantity of this system is a real-valued continuous func- tion E(x) that is constant on trajectories i.e. dE/dt = 0. Example 4.7.1 mz¨ = dV (z)/dz = F (z) − Take x = z and y =z ˙ ⇒ x˙ = y 1 y˙ = F (x) m 1 2 E(z)= 2mz˙ +V (z) is the total energy, which is constant 1 E(x) my2 + V (x) ⇒ ≡ 2 dE(x) = 0 dt since total energy is constant. 23 Example 6.5.2 θ¨+ sin θ = 0 e.g.
Load more

Recommended publications
  • Ordinary Differential Equation 2
    Unit - II Linear Systems- Let us consider a system of first order differential equations of the form dx = F(,)x y dt (1) dy = G(,)x y dt Where t is an independent variable. And x & y are dependent variables. The system (1) is called a linear system if both F(x, y) and G(x, y) are linear in x and y. dx = a ()t x + b ()t y + f() t dt 1 1 1 Also system (1) can be written as (2) dy = a ()t x + b ()t y + f() t dt 2 2 2 ∀ = [ ] Where ai t)( , bi t)( and fi () t i 2,1 are continuous functions on a, b . Homogeneous and Non-Homogeneous Linear Systems- The system (2) is called a homogeneous linear system, if both f1 ( t )and f2 ( t ) are identically zero and if both f1 ( t )and f2 t)( are not equal to zero, then the system (2) is called a non-homogeneous linear system. x = x() t Solution- A pair of functions defined on [a, b]is said to be a solution of (2) if it satisfies y = y() t (2). dx = 4x − y ........ A dt Example- (3) dy = 2x + y ........ B dt dx d 2 x dx From A, y = 4x − putting in B we obtain − 5 + 6x = 0is a 2 nd order differential dt dt 2 dt x = e2t equation. The auxiliary equation is m2 − 5m + 6 = 0 ⇒ m = 3,2 so putting x = e2t in A, we x = e3t obtain y = 2e2t again putting x = e3t in A, we obtain y = e3t .
    [Show full text]
  • Phase Plane Diagrams of Difference Equations
    PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS Abstract. We will be determining qualitative features of a dis- crete dynamical system of homogeneous difference equations with constant coefficients. By creating phase plane diagrams of our system we can visualize these features, such as convergence, equi- librium points, and stability. 1. Introduction Continuous systems are often approximated as discrete processes, mean- ing that we look only at the solutions for positive integer inputs. Dif- ference equations are recurrence relations, and first order difference equations only depend on the previous value. Using difference equa- tions, we can model discrete dynamical systems. The observations we can determine, by analyzing phase plane diagrams of difference equa- tions, are if we are modeling decay or growth, convergence, stability, and equilibrium points. For this paper, we will only consider two dimensional systems. Suppose x(k) and y(k) are two functions of the positive integers that satisfy the following system of difference equations: x(k + 1) = ax(k) + by(k) y(k + 1) = cx(k) + dy(k): This system is more simply written in matrix from z(k + 1) = Az(k); x(k) a b where z(k) = is a column vector and A = . If z(0) = z y(k) c d 0 is the initial condition then the solution to the initial value prob- lem z(k + 1) = Az(k); z(0) = z0; 1 2 TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS is k z(k) = A z0: Each z(k) is represented as a point (x(k); y(k)) in the Euclidean plane.
    [Show full text]
  • Phase Plane Methods
    Chapter 10 Phase Plane Methods Contents 10.1 Planar Autonomous Systems . 680 10.2 Planar Constant Linear Systems . 694 10.3 Planar Almost Linear Systems . 705 10.4 Biological Models . 715 10.5 Mechanical Models . 730 Studied here are planar autonomous systems of differential equations. The topics: Planar Autonomous Systems: Phase Portraits, Stability. Planar Constant Linear Systems: Classification of isolated equilib- ria, Phase portraits. Planar Almost Linear Systems: Phase portraits, Nonlinear classi- fications of equilibria. Biological Models: Predator-prey models, Competition models, Survival of one species, Co-existence, Alligators, doomsday and extinction. Mechanical Models: Nonlinear spring-mass system, Soft and hard springs, Energy conservation, Phase plane and scenes. 680 Phase Plane Methods 10.1 Planar Autonomous Systems A set of two scalar differential equations of the form x0(t) = f(x(t); y(t)); (1) y0(t) = g(x(t); y(t)): is called a planar autonomous system. The term autonomous means self-governing, justified by the absence of the time variable t in the functions f(x; y), g(x; y). ! ! x(t) f(x; y) To obtain the vector form, let ~u(t) = , F~ (x; y) = y(t) g(x; y) and write (1) as the first order vector-matrix system d (2) ~u(t) = F~ (~u(t)): dt It is assumed that f, g are continuously differentiable in some region D in the xy-plane. This assumption makes F~ continuously differentiable in D and guarantees that Picard's existence-uniqueness theorem for initial d ~ value problems applies to the initial value problem dt ~u(t) = F (~u(t)), ~u(0) = ~u0.
    [Show full text]
  • Calculus and Differential Equations II
    Calculus and Differential Equations II MATH 250 B Linear systems of differential equations Linear systems of differential equations Calculus and Differential Equations II Second order autonomous linear systems We are mostly interested with2 × 2 first order autonomous systems of the form x0 = a x + b y y 0 = c x + d y where x and y are functions of t and a, b, c, and d are real constants. Such a system may be re-written in matrix form as d x x a b = M ; M = : dt y y c d The purpose of this section is to classify the dynamics of the solutions of the above system, in terms of the properties of the matrix M. Linear systems of differential equations Calculus and Differential Equations II Existence and uniqueness (general statement) Consider a linear system of the form dY = M(t)Y + F (t); dt where Y and F (t) are n × 1 column vectors, and M(t) is an n × n matrix whose entries may depend on t. Existence and uniqueness theorem: If the entries of the matrix M(t) and of the vector F (t) are continuous on some open interval I containing t0, then the initial value problem dY = M(t)Y + F (t); Y (t ) = Y dt 0 0 has a unique solution on I . In particular, this means that trajectories in the phase space do not cross. Linear systems of differential equations Calculus and Differential Equations II General solution The general solution to Y 0 = M(t)Y + F (t) reads Y (t) = C1 Y1(t) + C2 Y2(t) + ··· + Cn Yn(t) + Yp(t); = U(t) C + Yp(t); where 0 Yp(t) is a particular solution to Y = M(t)Y + F (t).
    [Show full text]
  • Phaser: an R Package for Phase Plane Analysis of Autonomous ODE Systems by Michael J
    CONTRIBUTED RESEARCH ARTICLES 43 phaseR: An R Package for Phase Plane Analysis of Autonomous ODE Systems by Michael J. Grayling Abstract When modelling physical systems, analysts will frequently be confronted by differential equations which cannot be solved analytically. In this instance, numerical integration will usually be the only way forward. However, for autonomous systems of ordinary differential equations (ODEs) in one or two dimensions, it is possible to employ an instructive qualitative analysis foregoing this requirement, using so-called phase plane methods. Moreover, this qualitative analysis can even prove to be highly useful for systems that can be solved analytically, or will be solved numerically anyway. The package phaseR allows the user to perform such phase plane analyses: determining the stability of any equilibrium points easily, and producing informative plots. Introduction Repeatedly, when a system of differential equations is written down, it cannot be solved analytically. This is particularly true in the non-linear case, which unfortunately habitually arises when modelling physical systems. As such, it is common that numerical integration is the only way for a modeller to analyse the properties of their system. Consequently, many software packages exist today to assist in this step. In R, for example, the package deSolve (Soetaert et al., 2010) deals with many classes of differential equation. It allows users to solve first-order stiff and non-stiff initial value problem ODEs, as well as stiff and non-stiff delay differential equations (DDEs), and differential algebraic equations (DAEs) up to index 3. Moreover, it can tackle partial differential equations (PDEs) with the assistance ReacTran (Soetaert and Meysman, 2012).
    [Show full text]
  • Phase-Plane Methods
    NPS ARCHIVE 1969 HSIUNG, Y. PHASE-PLANE METHODS Yun-Lan Hsiung DUDLEY KNOX LIBRARY NAVAL POSTGRADUATE SCHOOL MONTEREY, CA 93943-5101 United States Naval Postgraduate School THESIS PHASE-PLANE METHODS by Yun-Lan Hsiung October 1969 Tlnu document, kcud been app\ove.d fan. pubtic *e- leo*e and 6aJU; itb duVubution <u unlimUzd. T13S516 DUDLEY KNOX LIBRARY SCHOOL N M/AL POSTGRADUATE MONTEREY, CA 93943-5101 Phase -Plane Methods by Yun-Lan Hsiung Lieutenant Commander, Chinese Navy B. S., Chinese Naval Academy, 1954 Submitted in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE IN ELECTRICAL ENGINEERING from the NAVAL POSTGRADUATE SCHOOL October 1969 c \ \<o'=\ ABSTRACT The phase-plane method is a graphical method for linear and non-linear second-order systems. There are many techniques for obtaining the phase portrait which consists of a number of phase trajectories on the phase- plane. From the summary and discussion the isocline method is a most general and useful method. The only difficulty is in the labor required. Solution by digi- tal computer reduces this difficulty and makes the iso- cline method a much more useful method in non-linear systems analysis and design applications. LIBRARY IAYAL POSTGRADUATE SCHOOL MOHTJEBJSy, CALIF. 93940 TABLE OF CONTENTS Page I. INTRODUCTION 15 A. CONDITIONS FOR LINEARITY AND DEFINITION 15 OF NONLINEARITY B. CLASSES OF NONLINEARITY 15 C. NONLINEAR ANALYSIS 19 D. THE GRAPHICAL METHOD - THE PHASE-PLANE 20 (PHASE-PLANE ANALYSIS OF NONLINEAR SYSTEMS) II. SUMMARY OF USEFUL GRAPHICAL METHODS FOR 22 CONSTRUCTING PHASE TRAJECTORIES ON THE PHASE-PLANE A.
    [Show full text]
  • 7 Phase Planes
    F1.3YT2/YF3 67 7 Phase Planes 7.1 Introduction As we have discussed in detail in Section 2, a system of N ODEs can be expressed as x_ (t)=F(t; x(t)) where x1(t) . x : R RN ; x(t)= . (7.1) ! 0 1 x(t) BN C @ A F1(t; x1; ;xN) N+1 N .¢¢¢ F : R R ; F (t; x)= . (7.2) ! 0 1 F (t; x ; ;x ) BN 1 N C @ ¢¢¢ A A system is called autonomous if it is of the form x_ (t)=F(x(t)) (7.3) i.e., if there is no explicit t dependence on the right hand side. Example 7.1. (i) x_(t)=tx(t) is not autonomous. (ii) x_(t)=x2(t)is autonomous. (iii) The equation xÄ(t)+x _(t) + sin (x(t)) = 0 can be expressed as the autonomous system x_ 1(t)=x2(t) x_(t)= sin (x (t)) x (t): 2 ¡ 1 ¡ 2 where x(t)=x1(t). If we assume that F is su±ciently smooth, then Picard's theorem tells us that the solution to the initial value problem x_ (t)=F(x(t)); x(t0)=x0 (7.4) is unique. De¯nition 7.2. The set of points in x(t) RN satisfying (7.4) is called the trajectory 2 for (7.4) passing through x0. F1.3YT2/YF3 68 x2 x² x1 Physical Interpretation (N=2) When N = 2, we tend to use the notation (x; y) instead of (x1;x2). So we have a systems of autonomous equations of the form x_(t)=f(x(t);y(t)) (7.5) y_(t)=g(x(t);y(t)) An interpretation of these equations is that they specify the two components of the velocity of a particle moving in the x y plane.
    [Show full text]
  • 6 Phase Plots and Vector Field Plots Fall 2003
    6 Phase Plots and Vector Field Plots Fall 2003 Phase Plots When a particle moves along the x axis, we often represent the motion as a graph of the coordinate x or the velocity v = dx/dt as a function of time t. A different but very useful representation is to plot the instantaneous position x and velocity v of the particle as a point in a plane called the phase plane, with horizontal and vertical axes representing x(t) and v(t), respectively. Such a plot is called a phase plot. Each point in the x-v phase plane represents an instantaneous state of motion (position and velocity) of the system. As the motion progresses, the representative point (the phase point) traces out a path called the phase trajectory in the phase plane. A simple example is the phase plot for the undamped, undriven harmonic oscillator. The total energy E of the system is constant, and conservation of energy gives the equation 1 2 1 2 2 mv + 2 kx = E = constant. (1) The graph of this equation in the x-v plane (i.e., the phase plot) is an ellipse. As the motion evolves in time, the phase point moves around this ellipse, tracing out the phase plot once each cycle. It always moves in the clockwise sense; can you see why? For every periodic motion, the phase plot is a closed curve that is traced out once each cycle. When damping is present, the motion is not strictly periodic. The phase trajectory is no longer a closed curve but a spiral that curves into the origin as the motion dies down.
    [Show full text]
  • Appendix A. Phase Plane Analysis
    Appendix A. Phase Plane Analysis We discuss here, only very briefly, general autonomous second-order ordinary differen- tial equations of the form dx dy = f (x, y), = g(x, y). (A.1) dt dt We present the basic results which are required in the main text. There are many books which discuss phase plane analysis in varying depth, such as Jordan and Smith (1999) and Guckenheimer and Holmes (1983). A good, short and practical exposition of the qualitative theory of ordinary differential equation systems, including phase plane tech- niques, is given by Odell (1980). Phase curves or phase trajectories of (A.1) are solu- tions of dx f (x, y) = . (A.2) dy g(x, y) Through any point (x0, y0) there is a unique curve except at singular points (xs, ys) where f (xs, ys) = g(xs, ys) = 0. Let x → x −xs, y → y − ys; then, (0, 0) is a singular point of the transformed equation. Thus, without loss of generality we now consider (A.2) to have a singular point at the origin; that is, f (x, y) = g(x, y) = 0 ⇒ x = 0, y = 0. (A.3) If f and g are analytic near (0, 0) we can expand f and g in a Taylor series and, retaining only the linear terms, we get dx ax + by ab f f = , A = = x y + (A.4) dy cx dy cd gx gy (0,0) which defines the matrix A and the constants a, b, c and d. The linear form is equivalent to the system 502 Appendix A.
    [Show full text]
  • Phase Plane Analysis of Differential Equations
    Phase-plane Analysis of Ordinary November 1, 2017 Differential Equations Outline Phase-plane Analysis of Ordinary Differential Equations • Midterm exam two weeks from tonight covering ODEs and Laplace transforms Larry Caretto • Review last class Mechanical Engineering 501A • Introduction to phase-plane analysis Seminar in Engineering Analysis • Look at two simultaneous ODEs dy1/dt and dy2/dt plotted as y1 vs. y2 November 1, 2017 • Look at different “critical points” for different systems of equations 2 Review Laplace Transforms Review Partial Fractions • Use transform tables to transform terms • Method to convert fraction with several in differential equation for y(t) into an factors in denominator into sum of algebraic equation for Y(s) individual factors (in denominator) – Derivative transforms give initial conditions • Example is F(s) = 1/(s+a)(s+b) on y(t) and its derivatives • Write 1/(s+a)(s+b) = A/(s+a) + B/(s+b) • Manipulate Y(s) equation to sum of • Multiply by (s+a)(s+b) and equate individual terms, “Y(s) subterms”, that coefficients of like powers of s you can find in transform tables – 1 = A(s + b) + B(s + a) – Manipulation may require use of method of – A + B = 0 for s1 terms and 1 = bA + aB for partial fractions s0 terms 3 4 Review Partial Fractions II Review Partial Fraction Rules – A + B = 0 for s1 terms and 1 = bA + aB for • Repeated fractions for repeated factors 1 A A A A s0 terms n n1 2 1 n n n1 2 – Solving for A and B gives A = -B = 1/(b – a) (s a) (s a) (s a) (s a) s a • Complex factors (s + –i)(s +
    [Show full text]
  • GS. Graphing ODE Systems
    GS. Graphing ODE Systems 1. The phase plane. Up to now we have handled systems analytically, concentrating on a procedure for solving linear systems with constant coefficients. In this chapter, we consider methods for sketching graphs of the solutions. The emphasis is on the workd sketching. Computers do the work of drawing reasonably accurate graphs. Here we want to see how to get quick qualitative information about the graph, without having to actually calculate points on it. First some terminology. The sort of system for which we will be trying to sketch the solutions is one which can be written in the form ′ x = f(x,y) (1) ′ . y = g(x,y) Such a system is called autonomous, meaning the independent variable (which we under- stand to be t) does not appear explicitly on the right, though of course it lurks in the derivatives on the left. The system (1) is a first-order autonomous system; it is in standard form — the derivatives on the left, the functions on the right. A solution of such a system has the form (we write it two ways): x(t) x = x(t) (2) x(t) = , . y(t) y = y(t) It is a vector function of t, whose components satisfy the system (1) when they are substi- tuted in for x and y. In general, you learned in 18.02 and physics that such a vector function describes a motion in the xy-plane; the equations in (2) tell how the point (x,y) moves in the xy-plane as the time t varies.
    [Show full text]
  • Boyce/Diprima 9Th Ed, Ch 7.5: Homogeneous Linear Systems With
    Boyce/DiPrima 9th ed, Ch 7.5: Homogeneous Linear Systems with Constant Coefficients! Elementary Differential Equations and Boundary Value Problems, 9th edition, by William E. Boyce and Richard C. DiPrima, ©2009 by John Wiley & Sons, Inc. ! " We consider here a homogeneous system of n first order linear equations with constant, real coefficients: ! " This system can be written as x' = Ax, where Equilibrium Solutions ! " Note that if n = 1, then the system reduces to ! " Recall that x = 0 is the only equilibrium solution if a ≠ 0. ! " Further, x = 0 is an asymptotically stable solution if a < 0, since other solutions approach x = 0 in this case. ! " Also, x = 0 is an unstable solution if a > 0, since other solutions depart from x = 0 in this case. ! " For n > 1, equilibrium solutions are similarly found by solving Ax = 0. We assume detA ≠ 0, so that x = 0 is the only solution. Determining whether x = 0 is asymptotically stable or unstable is an important question here as well. Phase Plane ! " When n = 2, then the system reduces to ! " This case can be visualized in the x1x2-plane, which is called the phase plane. ! " In the phase plane, a direction field can be obtained by evaluating Ax at many points and plotting the resulting vectors, which will be tangent to solution vectors. ! " A plot that shows representative solution trajectories is called a phase portrait. ! " Examples of phase planes, directions fields and phase portraits will be given later in this section. Solving Homogeneous System ! " To construct a general solution to x' = Ax, assume a solution of the form x = ξert, where the exponent r and the constant vector ξ are to be determined.
    [Show full text]