|
|
- Eunice Griffin
- 5 years ago
- Views:
Transcription
1 Part 1: Overview of Ordinary Dierential Equations 1
2
3 Chapter 1 Basic Concepts and Problems 1.1 Problems Leading to Ordinary Dierential Equations Many scientic and engineering problems are modeled by systems of ordinary dierential equations (ODEs). Some examples [] follow. 1. Mechanical Vibrations. Let y(t) denote the displacement at time t of a block of mass m that is connected to a spring of stiness k, a damper of resistance c, and an oscillator f(t) (Figure 1.1.1). If the system is released from position y 0 velocity y 0 then its subsequent motion satises 0 with a m d y dt + c dy dt + ky = f(t) t>0 y(0) = y 0 dy(0) dt = y 0 0 : This is an example of an initial value problem (IVP) for a second-order linear ODE. The variable y(t) is called: The variable t is called: The ODE is second order because: It's linear because: This is an IVP because:
4 k m c f(t) y(t) Figure 1.1.1: Motion of a spring-mass-damper system.. Ecology. Consider a population of predators and prey living in an ecological niche. The predators survive by eating the prey and the prey exist on an independent source of food. A classical situation involves foxes and rabbits. Let P (t) and p(t), respectively, denote the populations of predators and prey at time t. Given the initial populations P 0 and p 0 of predators and prey, their subsequent populations satisfy the Lotka-Volterra equations dp dt = p(a ; P ) t>0 p(0) = p 0 dp dt = P (;c + p) t>0 P (0) = P 0 where a, c,, and are positive constants corresponding to the prey's natural growth rate, the predator's natural death rate, the prey's death rate upon coming into contact with predators, and the predator's growth rate upon coming into contact with prey. This example involves an IVP for a system of two rst-order nonlinear ODEs. The ODEs are nonlinear because: The system is rst order because:. Column Buckling. The lateral displacement y(t) at position t of a clamped-hinged bar of length l that is subjected to a load P (Figure 1.1.) may be approximated by the Euler-Bernoulli equations d 4 y dt 4 with the boundary conditions y(0) = dy(0) dt + d y =0 0 <t<l dt =0 y(l) = d y(l) dt =0: 4
5 The parameter = P=EI, where EI is the exural rigidity ofthe bar. This is an example of a boundary value problem (BVP) for a fourth-order linear ODE. This is a boundary value problem because: Observe that y(t) = 0 is a solution of this problem for all values of and l. A more interesting problem is to determine y and those values of and l for which non-trivial (y(t) 6= 0) solutions exist. As such, this BVP is also a dierential eigenvalue problem. Nontrivial solutions y(t) are called eigen functions and their corresponding values of are eigenvalues. P t y(t) P Figure 1.1.: Buckling of an elastic column. 4. Pendelum Oscillations ([1], Section 1.). The position (x(t) y(t)) at time t of a particle of mass m oscillating on a pendelum of length l (Figure 1.1.) is m d x dt = ;T sin = ; T l y x md = mg ; T cos = mg ; T y t>0 dt l where g is the acceleration of gravity, T (t) is the tension in the string, and (t) is the angle of the pendelum relative to the vertical at time t (Figure 1.1.). These equations are, however, insucient to guarantee that the particle stays on the string. To ensure that this is so, we must supplement the ODEs by the algebraic constraint x + y = l : The two second-order dierential equations and the constraint comprise a system of ve dierential algebraic equations (DAEs) for the unknowns x(t), y(t), and T (t). Initial conditions specify x(0), dx(0)=dt, y(0), and dy(0)=dt, but not T (0). 5
6 Of course, for this problem, it's easy to eliminate the constraint byintroducing the change of variables x = l sin, y = l cos. This would reduce the DAEs to the second-order ODE d = ; g sin : dt l Such simplications would not be possible with more dicult systems. θ y T x mg Figure 1.1.: Oscillations of a simple Pendelum. IVPs will comprise our initial study (Part of these notes). We'll take up BVPs next (Part ) and conclude with a study of DAEs (Part 4). In almost all cases, it will suce to develop and analyze numerical methods for IVPs and BVPs that are written as rst-order vector systems of ODEs in the explicit form y 0 (t) =f(t y) t>0 (1.1.1a) where y(t) = 6 4 y 1 (t) y (t). y m (t) f(t y) = f 1 (t y 1 ::: y m ) f (t y 1 ::: y m ). f m (t y 1 ::: y m ) 5 (1.1.1b) and y 0 := dy=dt. Data for IVPs would is specied as y(0) = y 0 = 6 4 y 10 y 0. 5 : (1.1.1c) y m0 6
7 Boundary data is more complex. The most general boundary conditions specify a nonlinear relationship between y(0) and y(l) oftheform g 1 (y(0) y(l)) g (y(0) y(l)) g(y(0) y(l)) = 6 = 0: (1.1.1d) 4. 5 g m (y(0) y(l)) In many cases the extension of a scalar to a vector IVP will be obvious and we can study methods for a scalar IVP y 0 (t) =f(t y) t>0 y(0) = y 0 (1.1.) Example Higher-order ODEs can be written as rst-order systems. Consider the m th order equation z (m) = g(t z z 0 ::: z (m;1) ) where z (i) := d i z=dt i, and introduce the new variables y 1 = z y = z 0 y = z 00 ::: y m = z (m;1) : Then, we have a system in the form (1.1.1) with y 0 1 = y y 0 = y ::: y 0 m = g(t y 1 y ::: y m ): For an IVP, that data would prescribe as z (i) (0) = c i i =0 1 ::: m; 1 with the understanding that z (0) = z. Written in terms of y, we have y 1 (0) = c 0 y (0) = c 1 ::: y m (0) = c m;1 : With DAEs and sometimes with IVPs and BVPs, it will be necessary to consider implicit dierential systems F(t y y 0 )=0: (1.1.) If the 0 is nonsingular then, by the implicit function theorem, (1.1.) is equivalent to (1.1.1) however, it may be natural to solve (1.1.) in its implicit form to maintain, e.g., sparsity. DAEs involve systems 0 is singular.
8 Example Let us write the oscillating pendelum problem as a rst-order system by letting y 1 = x y = x 0 y = y y 4 = y 0 : Then, we have y 0 1 = y y 0 = ; T ml y 1 y 0 = y 4 y 0 4 = g ; T ml y y 5 = T with the constraint y 1 + y = l : This can be written in the form (1.1.) with y 0 ; y 1 y 0 + y 1y 5 F(t y y 0 ml )= 6 y 0 ; y 4 4 y 0 ; g + y y 5 4 ml y + y ; l 1 5 = 0: The Jacobian is clearly singular. F y = Problems 5 1. Letting y 1 = p and y = P, write the predator-prey model dp dt = p(a ; P ) t>0 p(0) = p 0 dp dt = P (;c + p) t>0 P (0) = P 0 in the vector form (1.1.1).. Write the mechanical vibration problem m d y dt in the form (1.1.1a-c) + c dy dt + ky = f(t) t>0 y(0) = y 0 8 dy(0) dt = y 0 0
9 . Write the column buckling problem d 4 y dt 4 + d y =0 0 <t<l y(0) = dy(0) y(l) =0 y(l) = d =0 dt dt dt in the form (1.1.1a,b,d) 1. Existence, Uniqueness, and Stability of IVPs Before considering numerical procedures for (1.1.1, 1.1.), let us review some results from ODE theory that ensure the existence of unique solutions which depend continuously on the initial data. For simplicity, letusfocus on a scalar IVP having the form (1.1.). Denition A function f(t y) satises a Lipschitz condition in a domain D if there exists a non-negative constant L such that jf(t y) ; f(t z)j Ljy ; zj 8(t y) (t z) D: (1..4) Satisfying a Lipschitz condition guarantees that solutions y(t) of (1.1.) are unique as expressed by the following theorem. Theorem Let f(t y) be continuous and satisfy a Lipschitz condition on f(t y) j 0 t T ;1 <y<1g: Then the IVP (1.1.) has a unique continuously dierentiable solution on [0 T] for all y 0 (;1 1). Proof. The proof appears in most books on ODE theory, e.g., []. Remark 1. A Lipschitz conditions guarantees unique solutions and is not needed for esistence. For example, the IVP y 0 =y = y(0) = 0 has the two solution y(t) = 0 and y(t) = t. This f(t y) = y = Lipschitz condition. does not satisfy a 9
10 Remark. Theorem 1..1 applies when f satisfys a Lipschitz condition on a compact, rather than an unbounded, domain D as long as the solution y(t) remains in D []. In addition to existence and uniqueness, we will want to know something about the stability of solutions of the IVP. In particular, we will usually be interested in the sensitivity of the solution to small changes in the data. Perturbations arise naturally in numerical computation due to discretization and round o errors. A formal study of sensitivity would lead us to the following notion of a well-posed system. Denition 1... The IVP (1.1.) is well-posed if there exists positive constants k and ^ such that, for any ^, the perturbed IVP z 0 = f(t z)+(t) z(0) = y (1..5) satises jy(t) ; z(t)j k (1..6) whenever j 0 j <and j(t)j <for t [0 T]. Again, a Lipschitz condition ensures that we are dealing with a well-posed IVP. Theorem 1... If f(t y) satises a Lipschitz condition on f(t y) j 0 t T ;1 <y<1g then y 0 = f(t y) is well posed on [0 T] with respect to all initial data. Proof. Let (t) =z(t) ; y(t) and subtract (1.1.) from (1..5) to obtain 0 (t) =f(t z) ; f(t y)+(t) (0) = 0 : Taking an absolute value and using the Lipschitz condition (1..4) j 0 (t)j Lj(t)j + j(t)j j(0)j = 0 : 10
11 We easily see that j(t)j 0 j 0 (t)j provided that j(t)j 0 exists and it is fairly easy to show that j(t)j 0 exists. Additionally, by assumption, max j(t)j < 0tT max j 0 j < 0tT so j(t)j 0 Lj(t)j + j(0)j <: Multiply by the integrating factor e ;Lt to obtain (e ;Lt j(t)j) 0 e ;Lt j(0)j <: Integrating j(t)j Since the denominator is smallest when t = 0 jz(t) ; y(t)j thus, y 0 = f(t y) iswell posed on [0 T]. L[(L +1)e Lt ; 1] : L = k 8t [0 T] The notion of awell-posed system is related to the more common notion of stability as indicated by the following denition. Denition 1... Consider the dierential equation y 0 = f(t y) and (without loss of generality) let the origin (0 0) be anequilibrium point, i.e., f(0 0) = 0. Then, the origin is: stable if a perturbation of the initial condition jy(0)j <grows no larger than for subsequent times, i.e., if jy(t)j <for t>0. asymptotically stable if it is stable and jy(0)j <implies that lim t!1 jy(t)j =0. unstable if it is not stable. Remark. This denition could, like Denition 1.., also involve perturbations of f(t y). We have omitted these for simplicity. 11
12 An autonomous system is one where f(t y) does not explicitly depend on t, i.e., f(t y) =f(y). If (0 0) is an equilibrium point then, in this case, f(0) = 0. Expanding the solution in a Taylor's series in y, we have y 0 (t) =f(y) =f(0) + f y (0)y(t)+O(y ): Since f(0) = 0, y 0 (t) =f y (0)y(t)+O(y ) where f y =@y. Recall that a function g(y) iso(y p ) if there exists a constant C > 0 such that jg(y)j Cy p as y! 0. Letting = f y (0), we see that the stability of an autonomous system is related to that of the simple linear IVP y 0 = y t>0 y(0) = : (1..) Additional details on autonomous and non-autonomous systems appear in Birkho and Rota [] or Hairer et al. [5], Section 1.1. The solution of the linear IVP (1..) is y(t) =e t hence, (1..) is stable when Re() 0, asymptotically stable when Re() < 0, and unstable when Re() > 0. These conclusions remain true for the original nonlinear autonomous problem when Re() 6= 0 and y is small enough for the O(y ) term to be negligible relative toy. This cannot happen in the stable case (Re() = 0) hence, it requires a more careful analysis. Remark 4. The analysis of non-autonomous systems is similar [1]. Analyzing the stability of autonomous vector systems y 0 = f(y(t)) (1..8) 1
13 is more complicated. Once again, assume that y = 0 is an equilibrium point and expand f(0) in a series y 0 = f(0)+f y (0)y + O(kyk )=f y (0)y + O(kyk ): We need a brief digression for a few denitions. Denition The Jacobian matrix of a vector-valued function f(y) with respect to y is the matrix f = @ym Denition The norm of a vector y is a scalar kyk such that 1. kyk 0 and kyk =0if and only if y = 0,. kyk = jjkyk for any scalar, and. kx + yk kxk + kyk. We'll identify specic vector norms when they are needed. 5 : (1..9) For the moment, let's return to our stability analysis and let A = f y (0). Then, the stability of the nonlinear autonomous system (1..8) is related to that of the linear system The solution of this system is y 0 = Ay y(0) = y 0 ky 0 k: (1..10a) y(t) =e At y 0 (1..10b) where the matrix exponential is dened by the series Often, A can be diagonalized as e At = I + ta + t + t + : (1..10c)! A! A 1 T ;1 AT = = m 5 (1..11) 1
14 where 1 ::: m are the eigenvalues of A and the columns of T are the corresponding eigenvectors. In this case, we may easily verify that the solution of (1..10) is y(t) =Te t T ;1 y 0 : Thus, (1..10) is stable if all of the eigenvalues have non-positive real parts, asymptotically stable if all of the eigenvalues have negative real parts, and unstable otherwise. Unfortunately, not all matrices are diagonalizable. However, A can always be reduced to the Jordan canonical form 1 T ;1 AT = 6 (1..1a) l where each Jordan block has the form i = 6 4 i 1 i (1..1b) i The dimension of the Jordan block i corresponds to the multiplicity ofthe eigenvalue i. Thus, if i is simple, the block is a scalar. With this, it is relatively easy to show [] that y = 0 is stable when either Re( i ) < 0orRe( i )=0and i is simple, i =1 ::: m, asymptotically stable when Re( i ) < 0, i =1 ::: m, and unstable otherwise. Example Consider the predator-prey model of Section 1.1 y 0 1 = y 1 (a ; y ) y 0 = y (;c + y 1 ) where a,, c, and are positive constants. This autonomous problem has two equilibrium points: y1 y = c= a= :
15 The point [0 0] T corresponds to extinction of predators and prey, whereas [c= a=] T implies a co-existence of both species. The stability of[0 0] T the linear system 0 y1 a 0 = y 0 ;c y1 y : may be analyzed by using Since the system matrix is diagonal, its eigenvalues are 1 = a, = ;c. Thus, 1 > 0 and [0 0] T is unstable. In order to analyze the stability of the second equilibrium point, we introduce the change of variables z 1 = y 1 ; c z = y ; a to move the equilibrium point to the origin. Substituting this transformation into the predator-prey equations gives z 0 1 =(z 1 + c )[a ; (z + a )] = ;z ( c + z 1) z 0 =(z + a )[;c + (z 1 + c )] = z 1( a + z ): The linearized system is now 0 z1 = z 0 ;c= a= 0 z1 z : The eigenvalues of this system are 1 = p i c=a. Since the eigenvalues are imaginary, the linear system is stable however, the stability of the nonlinear system is undetermined without additional analysis. Problems 1. =@y is continuous on a closed convex domain D. Show that f(t y) satises a Lipschitz condition on D with (Hint: use the Mean Value y) L = max j j: (t Are the following IVPs well posed ([4], p. )? Explain. y 0 = p 1 ; y y(0) = 0 y 0 = p y ; 1 y(0) = : 15
16 16
17 Bibliography [1] U.M. Ascher and L.R. Petzold. Computer Methods for Ordinary Dierential Equations and Dierential-Algebraic Equations. SIAM, Philadelphia, [] G. Birkho and G.-C. Rota. Ordinary Dierential Equations. John Wiley and Sons, New York, third edition, 198. [] W.E. Boyce and R.C. DiPrima. Elementary Dierential Equations and Boundary Value Problems. John Wiley and Sons, New York, third edition, 19. [4] C.W. Gear. Numerical Initial Value Problems in Ordinary Dierential Equations. Prentice Hall, Englewood Clis, 191. [5] E. Hairer, S.P. Norsett, and G. Wanner. Solving Ordinary Dierential Equations I: Nonsti Problems. Springer-Verlag, Berlin, second edition,
Math 266: Ordinary Differential Equations
Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationMATH 4B Differential Equations, Fall 2016 Final Exam Study Guide
MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide GENERAL INFORMATION AND FINAL EXAM RULES The exam will have a duration of 3 hours. No extra time will be given. Failing to submit your solutions
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 informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We
More information! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»
Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable
More information3.3. SYSTEMS OF ODES 1. y 0 " 2y" y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0.
.. SYSTEMS OF ODES. Systems of ODEs MATH 94 FALL 98 PRELIM # 94FA8PQ.tex.. a) Convert the third order dierential equation into a rst oder system _x = A x, with y " y" y + y = x = @ x x x b) The equation
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
More informationThe Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 2 Dianne P. O Leary c 2008 The Plan
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.
More informationLocal Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid)
Biology Medicine & Natural Product Chemistry ISSN: 089-6514 Volume 5 Number 1 016 Pages: 9-14 DOI: 10.1441/biomedich.016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationLab #2 - Two Degrees-of-Freedom Oscillator
Lab #2 - Two Degrees-of-Freedom Oscillator Last Updated: March 0, 2007 INTRODUCTION The system illustrated in Figure has two degrees-of-freedom. This means that two is the minimum number of coordinates
More informationNumerical Integration of Equations of Motion
GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.
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 informationA brief introduction to ordinary differential equations
Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)
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 informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationCoupled differential equations
Coupled differential equations Example: dy 1 dy 11 1 1 1 1 1 a y a y b x a y a y b x Consider the case with b1 b d y1 a11 a1 y1 y a1 ay dy y y e y 4/1/18 One way to address this sort of problem, is to
More informationChapter 13 Internal (Lyapunov) Stability 13.1 Introduction We have already seen some examples of both stable and unstable systems. The objective of th
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More information2nd-Order Linear Equations
4 2nd-Order Linear Equations 4.1 Linear Independence of Functions In linear algebra the notion of linear independence arises frequently in the context of vector spaces. If V is a vector space over the
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 informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More 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 informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationChapter 4. Systems of ODEs. Phase Plane. Qualitative Methods
Chapter 4 Systems of ODEs. Phase Plane. Qualitative Methods Contents 4.0 Basics of Matrices and Vectors 4.1 Systems of ODEs as Models 4.2 Basic Theory of Systems of ODEs 4.3 Constant-Coefficient Systems.
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 informationCoordinate Curves for Trajectories
43 The material on linearizations and Jacobian matrices developed in the last chapter certainly expanded our ability to deal with nonlinear systems of differential equations Unfortunately, those tools
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
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 information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationLecture 2: Review of Prerequisites. Table of contents
Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
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 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 informationInteractions. Yuan Gao. Spring Applied Mathematics University of Washington
Interactions Yuan Gao Applied Mathematics University of Washington yuangao@uw.edu Spring 2015 1 / 27 Nonlinear System Consider the following coupled ODEs: dx = f (x, y). dt dy = g(x, y). dt In general,
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 informationOrdinary Differential Equations
Ordinary Differential Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro September 19, 2014 1 / 55 Motivation
More informationLecture 6. Eigen-analysis
Lecture 6 Eigen-analysis University of British Columbia, Vancouver Yue-Xian Li March 7 6 Definition of eigenvectors and eigenvalues Def: Any n n matrix A defines a LT, A : R n R n A vector v and a scalar
More information(Refer Slide Time: 00:32)
Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra
More informationSolution of Linear State-space Systems
Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state
More informationMathematics for chemical engineers. Numerical solution of ordinary differential equations
Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species
More informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationA note on the uniform perturbation index 1
Rostock. Math. Kolloq. 52, 33 46 (998) Subject Classification (AMS) 65L5 M. Arnold A note on the uniform perturbation index ABSTRACT. For a given differential-algebraic equation (DAE) the perturbation
More informationLAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE
Page 1 of 15 LAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE Friday, December 14 th Course and No. 2007 MATH 2066 EL Date................................... Cours et no........................... Total no.
More informationMIDTERM REVIEW AND SAMPLE EXAM. Contents
MIDTERM REVIEW AND SAMPLE EXAM Abstract These notes outline the material for the upcoming exam Note that the review is divided into the two main topics we have covered thus far, namely, ordinary differential
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 informationMath 2a Prac Lectures on Differential Equations
Math 2a Prac Lectures on Differential Equations Prof. Dinakar Ramakrishnan 272 Sloan, 253-37 Caltech Office Hours: Fridays 4 5 PM Based on notes taken in class by Stephanie Laga, with a few added comments
More informationThe Exponential of a Matrix
The Exponential of a Matrix 5-8- The solution to the exponential growth equation dx dt kx is given by x c e kt It is natural to ask whether you can solve a constant coefficient linear system x A x in a
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationA(t) L T *(t) t/t t/t
he Real Story behind Floquet-Lyapunov Factorizations RAYMOND J. SPIERI Department of Mathematics and Statistics Acadia University Wolfville, NS BP X CANADA raymond.spiteri@acadiau.ca http//ace.acadiau.ca/~rspiteri
More informationEcon Lecture 14. Outline
Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationOur goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
More informationLecture Notes in Mathematics. Arkansas Tech University Department of Mathematics
Lecture Notes in Mathematics Arkansas Tech University Department of Mathematics Introductory Notes in Ordinary Differential Equations for Physical Sciences and Engineering Marcel B. Finan c All Rights
More informationThe Fundamental Theorem of Calculus: Suppose f continuous on [a, b]. 1.) If G(x) = x. f(t)dt = F (b) F (a) where F is any antiderivative
1 Calulus pre-requisites you must know. Derivative = slope of tangent line = rate. Integral = area between curve and x-axis (where area can be negative). The Fundamental Theorem of Calculus: Suppose f
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationSecond-Order Linear ODEs
Second-Order Linear ODEs A second order ODE is called linear if it can be written as y + p(t)y + q(t)y = r(t). (0.1) It is called homogeneous if r(t) = 0, and nonhomogeneous otherwise. We shall assume
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More 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 informationDifferential Equations and Linear Algebra Supplementary Notes. Simon J.A. Malham. Department of Mathematics, Heriot-Watt University
Differential Equations and Linear Algebra Supplementary Notes Simon J.A. Malham Department of Mathematics, Heriot-Watt University Contents Chapter 1. Linear algebraic equations 5 1.1. Gaussian elimination
More informationSolutions to the Review Questions
Solutions to the Review Questions Short Answer/True or False. True or False, and explain: (a) If y = y + 2t, then 0 = y + 2t is an equilibrium solution. False: (a) Equilibrium solutions are only defined
More informationLecture I Introduction, concept of solutions, application
S. Ghorai Lecture I Introduction, concept of solutions, application Definition. A differential equation (DE) is a relation that contains a finite set of functions and their derivatives with respect to
More informationFinal Exam Review. Review of Systems of ODE. Differential Equations Lia Vas. 1. Find all the equilibrium points of the following systems.
Differential Equations Lia Vas Review of Systems of ODE Final Exam Review 1. Find all the equilibrium points of the following systems. (a) dx = x x xy (b) dx = x x xy = 0.5y y 0.5xy = 0.5y 0.5y 0.5xy.
More informationEcon 204 Differential Equations. 1 Existence and Uniqueness of Solutions
Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation
More 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 informationMATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES
MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES J. WONG (FALL 2017) What did we cover this week? Basic definitions: DEs, linear operators, homogeneous (linear) ODEs. Solution techniques for some classes
More informationAMATH 351 Mar 15, 2013 FINAL REVIEW. Instructor: Jiri Najemnik
AMATH 351 Mar 15, 013 FINAL REVIEW Instructor: Jiri Najemni ABOUT GRADES Scores I have so far will be posted on the website today sorted by the student number HW4 & Exam will be added early next wee Let
More informationSufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations
/. Inst. Maths Applies (1977) 19, 377-383 Sufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations W. M. GETZ AND D. H. JACOBSON National Research Institute for Mathematical
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationMAT 22B - Lecture Notes
MAT 22B - Lecture Notes 4 September 205 Solving Systems of ODE Last time we talked a bit about how systems of ODE arise and why they are nice for visualization. Now we'll talk about the basics of how to
More informationNotes for Numerical Analysis Math 5466 by S. Adjerid Virginia Polytechnic Institute and State University (A Rough Draft) Contents Numerical Methods for ODEs 5. Introduction............................
More information4. Higher Order Linear DEs
4. Higher Order Linear DEs Department of Mathematics & Statistics ASU Outline of Chapter 4 1 General Theory of nth Order Linear Equations 2 Homogeneous Equations with Constant Coecients 3 The Method of
More informationy0 = F (t0)+c implies C = y0 F (t0) Integral = area between curve and x-axis (where I.e., f(t)dt = F (b) F (a) wheref is any antiderivative 2.
Calulus pre-requisites you must know. Derivative = slope of tangent line = rate. Integral = area between curve and x-axis (where area can be negative). The Fundamental Theorem of Calculus: Suppose f continuous
More informationOn The Numerical Solution of Differential-Algebraic Equations(DAES) with Index-3 by Pade Approximation
Appl. Math. Inf. Sci. Lett., No. 2, 7-23 (203) 7 Applied Mathematics & Information Sciences Letters An International Journal On The Numerical Solution of Differential-Algebraic Equations(DAES) with Index-3
More informationTHETA-LOGISTIC PREDATOR PREY
THETA-LOGISTIC PREDATOR PREY What are the assumptions of this model? 1.) Functional responses are non-linear. Functional response refers to a change in the rate of exploitation of prey by an individual
More informationLecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems
Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Shawn D. Ryan Spring 2012 Last Time: We finished Chapter 9: Nonlinear Differential Equations and Stability. Now
More informationIntroduction to Ordinary Differential Equations with Mathematica
ALFRED GRAY MICHAEL MEZZINO MARKA. PINSKY Introduction to Ordinary Differential Equations with Mathematica An Integrated Multimedia Approach %JmT} Web-Enhanced Includes CD-ROM TABLE OF CONTENTS Preface
More informationBoyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields
Boyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields Elementary Differential Equations and Boundary Value Problems, 11 th edition, by William E. Boyce, Richard C. DiPrima,
More informationThe Unifying View on Ordinary Differential Equations and Automatic Differentiation, yet with a Gap to Fill Alexander Gofen.
The Unifying View on Ordinary Differential Equations and Automatic Differentiation, yet with a Gap to Fill Alexander Gofen Conventions Differentiation is understood as for holomorphic functions in the
More informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
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 informationTopic 6: Projected Dynamical Systems
Topic 6: Projected Dynamical Systems John F. Smith Memorial Professor and Director Virtual Center for Supernetworks Isenberg School of Management University of Massachusetts Amherst, Massachusetts 01003
More informationNonlinear System Analysis
Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control
More informationLecture 19: Solving linear ODEs + separable techniques for nonlinear ODE s
Lecture 19: Solving linear ODEs + separable techniques for nonlinear ODE s Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth
More informationNon-Dimensionalization of Differential Equations. September 9, 2014
Non-Dimensionalization of Differential Equations September 9, 2014 1 Non-dimensionalization The most elementary idea of dimensional analysis is this: in any equation, the dimensions of the right hand side
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationMath Assignment 2
Math 2280 - Assignment 2 Dylan Zwick Spring 2014 Section 1.5-1, 15, 21, 29, 38, 42 Section 1.6-1, 3, 13, 16, 22, 26, 31, 36, 56 Section 2.1-1, 8, 11, 16, 29 Section 2.2-1, 10, 21, 23, 24 1 Section 1.5
More informationSMA 208: Ordinary differential equations I
SMA 208: Ordinary differential equations I First Order differential equations Lecturer: Dr. Philip Ngare (Contacts: pngare@uonbi.ac.ke, Tue 12-2 PM) School of Mathematics, University of Nairobi Feb 26,
More informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
More informationApplications of eigenvalues/eigenvectors
Example 1. Applications of eigenvalues/eigenvectors Predator-Prey models. Foxes and Rabbits share a large forest, with the foxes eating rabbits and rabbits eating the abundant vegetation. The sizes of
More informationInitial value problems for ordinary differential equations
Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP
More informationM A : Ordinary Differential Equations
M A 2 0 5 1: Ordinary Differential Equations Essential Class Notes & Graphics D 19 * 2018-2019 Sections D07 D11 & D14 1 1. INTRODUCTION CLASS 1 ODE: Course s Overarching Functions An introduction to the
More information