ME 406 Assignment #9 Solutions
|
|
- Valentine Wilkinson
- 5 years ago
- Views:
Transcription
1 ME 06 Assignment #9 Solutions sysid Mathematica 6.0., DynPac.0, ê8ê009 intreset; plotreset; imsize = 50; Problem The equations are x = -x, y = x +y, z = x + z. Let F = x + y +z. To show that the surface F = 0 is an invariant surface, we have show that points on the surface remain on the surface -- that is, df/dt evaluated on F = 0 is 0. We have df dt = x F x + y F y +z F z = -x(x) + (x +y) () + ( x + z) () = x + y +z, which is zero when evaluated on F = 0. Problem We define the system for DynPac. setstate@8x, y<d; setparm@8<d; slopevec = 9-x + y, x + y + x y=; sysname = "ProbNL"; eigval@80, 0<D 8-, < The eigenvalues are - and, hence neither has a zero real part, so the equilibrium is hyperbolic. It is in fact a saddle point, as DynPac tells us also.
2 Assignment#9sol.nb 0<D Abbreviations used in classifyd. L = linear, NL = nonlinear, R = repeated root. Z = one zero root, Z = two zero roots. This message printed once. unstable - saddle We now construct the linearized system. The matrix A of the linearized system is the derivative matrix evaluated at the equilibrium: A = dermatval@80, 0<D 88-, 0<, 8, << We save our nonlinear system and then define the linear system for DynPac. savesys@probnld; slopevec = A.statevec 8-x, x + y< sysname = "ProbLIN"; The manifolds of the linear system will be determined by the eigenvectors and eigenvalues. Eigensystem@AD 88-, <, 88-, <, 80, <<< Thus the stable manifold (associated with the eigenvalue -) is spanned by the vector {-,}. The unstable manifold (associated with the eigenvalue ) is spanned by {0,}. We plot these, using blue for stable, red for unstable. man = ParametricPlot@8 u, -u<, 8u, -, <, PlotStyle Ø RGBColor@0, 0, DD; man = ParametricPlot@80, u<, 8u, -, <, PlotStyle Ø RGBColor@, 0, 0DD;
3 Assignment#9sol.nb man<d ProbLIN It is actually just as easy to get these from DynPac by a numerical integration. We do that next. plrange = 88-.,.<, 8-.,.<<; asprat = ; axon = False; frameon = True; arrowflag = True; arrowvec = 8 ê <; rangeflag = True; ranger = plrange; t0 = 0.0; h = 0.0; nsteps = 800; eq = 80, 0<; eigstab = 8-, <; eigunstab = 80, <; eps = 0.0; initset = 8eq + eps * eigstab, eq - eps * eigstab, eq + eps * eigunstab, eq - eps * eigunstab< , 0.0<, 80.0, -0.0<, 80, 0.0<, 80, -0.0<< h = 8-0.0, -0.0, 0.0, 0.0<; setcolor@8blue, Blue, Red, Red<D;
4 Assignment#9sol.nb linport = portrait@initset, t0, h, nsteps,, D ProbLIN y Now we carry out the same construction for the full nonlinear system. restoresys@probnld; h = 8-0.0, -0.0, 0.0, 0.0<; x
5 Assignment#9sol.nb 5 nlport = portrait@initset, t0, h, nsteps,, D ProbNL y We combine these x labon = "Problem - Stable and Unstable Manifolds";
6 6 Assignment#9sol.nb nlportd Problem - Stable and Unstable Manifolds y x Problem We enter the matrix A. A = 880, ê, ê, ê, -7 ê <, 8 ê, ê, - ê, -5 ê, ê <, 8, - ê, - ê, - ê, 7 ê <, 85, 9 ê, - ê, -5 ê, 9 ê <, 85 ê, 5 ê,, -, 5 ê <<; MatrixForm@AD We get the eigenvalues and eigenvectors of A.
7 Assignment#9sol.nb 7 Eigensystem@AD :8-, -,,, <, ::-,,, 5, >, 80, 0, 0, 0, 0<, 80,, 0,, <, 80, 0, 0, 0, 0<, 80, 0, 0, 0, 0<>> We see that there is an eigenvalue - of multiplicity with one eigenvector, and an eigenvalue of multiplicity, also with one eigenvector. We will have to use generalized eigenvectors. We start with the eigenvalue -. Id = IdentityMatrix@5D 88, 0, 0, 0, 0<, 80,, 0, 0, 0<, 80, 0,, 0, 0<, 80, 0, 0,, 0<, 80, 0, 0, 0, << stabspan = NullSpace@MatrixPower@A + Id, DD 88, -, -, 0, <, 8-,,,, 0<< Thus the stable manifold is spanned by the two vectors stab = stabspan@@dd 8, -, -, 0, < stab = stabspan@@dd 8-,,,, 0< The original eigenvector must be in this space. We use this to verify our work so far. orig = 8- ê, ê, ê, 5 ê, <; Solve@C * stab + C * stab ã orig, 8C, C<D ::C Ø, C Ø 5 >> * stab + H5 ê L * stab :-,,, 5, > Everything checks. The most general initial condition which will decay to zero is the most general vector in this manifold, which is given by the expression below for arbitrary C and C. C * stab + C * stab 8 C - C, - C + C, - C + C, C, C< The unstable manifold is associated with the eigenvalue, and it is spanned by the generalized eigenvec-
8 8 Assignment#9sol.nb The unstable manifold is associated with the eigenvalue, and it is spanned by the generalized eigenvectors associated with that eigenvalue. unstabspan = NullSpace@MatrixPower@A - Id, DD 88-,, 0, 0, <, 8, 0, 0,, 0<, 8, -,, 0, 0<< Once again we check our work by verifying that the original eigenvector is in the space spanned by these. unstab = unstabspan@@dd 8-,, 0, 0, < unstab = unstabspan@@dd 8, 0, 0,, 0< unstab = unstabspan@@dd 8, -,, 0, 0< orig = 80,, 0,, <; Solve@C * unstab + C * unstab + C * unstab ã orig, 8C, C, C<D 88C Ø, C Ø, C Ø 0<< * unstab + * unstab 80,, 0,, < Problem We define the system for DynPac. setstate@8x, y, z<d; setparm@8<d; sysname = "Problem "; slopevec = 9-y + x z - x, x + y z + x y z, -z - x - y + z + SinAx E=; We find the eigenvalues at the origin. eigsys@80, 0, 0<D 88-, Â, -Â<, 880, 0, <, 8Â,, 0<, 8-Â,, 0<<< We have a complex conjugate pair, so the origin is a non-hyperbolic equilibrium, and its stability is not determined from the linearization. We see from the eigenvalues that there is a stable manifold of dimension and a center manifold of dimension. We now obtain a local approximation for the center manifold, which in turn will lead to a stability determination.
9 Assignment#9sol.nb 9 We see from the eigenvalues that the linear stable manifold is the z-axis, and the linear center manifold is the x-y plane. The nonlinear center manifold will have the form z = h(x, y) with h(0,0) = 0, h x (0,0) = 0, and h y (0,0) = 0. We differentiate this equation with respect to time to get () z = h x x + h y y. When the expressions for the time derivatives of x, y, and z are entered into (), with z replaced by h, we get the basic partial differential equation satisfied by h. We enter this equation into Mathematica. We define first a modification of the slope vector in which z is replaced by h. ClearAll@hD slopevec 9-x - y + x z, x + y z + x y z, -x - y - z + z + SinAx E= () dotvec = slopevec êê. z Ø h 9h x - x - y, x + h y + h x y, -h + h - x - y + SinAx E= Because we are going to keep only through cubic terms, we replace Sin[x ] with x. dotvec = dotvec ê. SinAx E Ø x 9h x - x - y, x + h y + h x y, -h + h - x + x - y = Now we expand h in a Taylor series, keeping terms through order. Note that the constant and linear terms are absent because of the conditions in (). h = a x + b x y + c y + d x + f x y + g x y + k y a x + d x + b x y + f x y + c y + g x y + k y Now we write z - h x x - h y y for Mathematica. hequat = dotvec@@dd - D@h, xd dotvec@@dd - D@h, yd dotvec@@dd -x - a x + x - d x - b x y - f x y - y - c y - g x y - k y + Ia x + d x + b x y + f x y + c y + g x y + k y M - I a x + d x + b y + 6 f x y + g y M I-x - y + x Ia x + d x + b x y + f x y + c y + g x y + k y MM - I b x + f x + c y + 6 g x y + k y M Ix + y Ia x + d x + b x y + f x y + c y + g x y + k y M + x y Ia x + d x + b x y + f x y + c y + g x y + k y MM This quantity must vanish. It is a power series in x and y, so terms of any given order must vanish separately.
10 0 Assignment#9sol.nb This quantity must vanish. It is a power series in x and y, so terms of any given order must vanish separately. We impose these conditions, starting with the quadratic terms. coeffxx = H ê L D@hequat, 8x, <D ê. Thread@8x, y< Ø 80, 0<D H- - a - bl coeffyy = H ê L D@hequat, 8y, <D ê. Thread@8x, y< Ø 80, 0<D H- + b - cl coeffxy = D@D@hequat, xd, yd ê. Thread@8x, y< Ø 80, 0<D a - b - c We solve these for a,b,c. quadans = Solve@8coeffxx ã 0, coeffxy ã 0, coeffyy ã 0<, 8a, b, c<d 88a Ø -, b Ø 0, c Ø -<< Thus the quadratic approximation to the center manifold is hquad = -x - y -x - y Now consider the orbits on the center manifold. We have xdot = slopevec@@dd ê. z Ø hquad -x - y + x I-x - y M ydot = slopevec@@dd ê. z Ø hquad x + y I-x - y M + x y I-x - y M Correct to third order, we have the following differential equations on the center manifold: dx dt = -y - x(x + y ) + higher order, dy dt = x -y(x + y ) + higher order. We convert to polar coordinates to get rr = -r + higher order. Equation () shows the stability of the equilibrium at the origin. () () Had there been no conclusion from this analysis, we would have had to return to the calculation of the center manifold and keep more terms.
ME 406 Assignment 6 Solutions
ME 406 Assignment 6 Solutions sysid Mathematica 6.0.3, DynPac 11.02, 3ê12ê2009 inreset; plotreset; Problem 1 ü Part a We look for equilibria. The first equation requires y = 0, and the second equation
More informationME 406 Trapping a Limit Cycle
orbtrap.nb 1 ME 406 Trapping a Limit Cycle sysid Mathematica 6.0.3, DynPac 11.01, 1ê13ê2009 intreset; plotreset; imsize = 250; ü Introduction In this notebook, we localize a limit cycle by orbit trapping.
More informationME 406 Bifurcations VII Subcritical Hopf Bifurcation
ME 406 Bifurcations VII Subcritical Hopf Bifurcation sysid Mathematica 4.1.2, DynPac 10.66, 3ê5ê2002 intreset; plotreset; 1. Introduction In this notebook, the seventh in a series of notebooks on bifurcations,
More informationME 406 Example of Stable and Unstable Manifolds
ME 46 Eample of Stable and Unstable Manifolds intreset; plotreset; ssid Mathematica 4..2, DnPac.66, 3ê8ê22 Introduction In this notebook, we stud an eample given in Differential Equations and Dnamical
More informationME 406 Further Results for the Project: Analysis of Solution Near the Spiral Equilibrium for a = 0
ME 406 Further Results for the Project: Analysis of Solution Near the Spiral Equilibrium for a = 0 1. Introduction sysid; plotreset; intreset; imsize = 250; Mathematica 6.0.3, DynPac 11.02, 4ê10ê2009 In
More informationME 406 Trapping a Limit Cycle
orbtrap.nb 1 ME 6 Trapping a Limit Ccle In[3]:= ssid Mathematica.1, DnPac 1.65, ê13ê In[33]:= In[3]:= intreset; plotreset; ü Introduction In this notebook, we localize a limit ccle b orbit trapping. The
More informationME 406 S-I-R Model of Epidemics Part 2 Vital Dynamics Included
ME 406 S-I-R Model of Epidemics Part 2 Vital Dynamics Included sysid Mathematica 6.0.3, DynPac 11.01, 1ê13ê9 1. Introduction Description of the Model In this notebook, we include births and deaths in the
More informationME 406 Predator-Prey Example Part 1 - Equilibrium
predpt1.nb 1 In[299]:= sysid Mathematica 6.0.3, DynPac 11.02, 2ê12ê2009 ME 406 Predator-Prey Example Part 1 - Equilibrium ü Introduction This is a simple model of an ecological system with three components:
More informationME 406 Orbital Stability
ME 406 Orbital Stability sysid Mathematica 4.., DynPac 0.66, 3êê00. Introduction In this notebook we develop some examples which illustrate the concept of orbital stability. We begin by reviewing the "standard"
More informationME 163 A Catalogue of Equilibria for Systems of Two Linear Autonomous First Order Equations
catalogued.nb 1 ME 163 A Catalogue of Equilibria for Sstems of Two Linear Autonomous First Order Equations Introduction In this notebook, we look at solutions of equations of the following form: d = a
More informationDYNAMICAL SYSTEMS Tutorial 17
map.nb DYNAMICAL SYSTEMS Tutorial 7 Iterated Maps sysid Mathematica 4..2, DynPac.67, 3ê9ê22 plotreset; intreset; ü Functions and Variables Used in This Tutorial asprat, bifurcmap, bifurc3dmap, bimap, borat,
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 informationHW3 - Due 02/06. Each answer must be mathematically justified. Don t forget your name. 1 2, A = 2 2
HW3 - Due 02/06 Each answer must be mathematically justified Don t forget your name Problem 1 Find a 2 2 matrix B such that B 3 = A, where A = 2 2 If A was diagonal, it would be easy: we would just take
More informationNonlinear Autonomous Dynamical systems of two dimensions. Part A
Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous
More informationMath 5490 November 5, 2014
Math 549 November 5, 214 Topics in Applied Mathematics: Introduction to the Mathematics of Climate Mondays and Wednesdays 2:3 3:45 http://www.math.umn.edu/~mcgehee/teaching/math549-214-2fall/ Streaming
More informationME 406 Using Eigenvector Methods with Mathematica to Solve Linear Autonomous Systems of First Order Differential Equations
ME 406 Using Eigenvector Methods with Mathematica to Solve Linear Autonomous Systems of First Order Differential Equations 1. Introduction In this notebook, we use the methods of linear algebra -- specifically
More informationStability of critical points in Linear Systems of Ordinary Differential Equations (SODE)
Stability of critical points in Linear Systems of Ordinary Differential Equations (SODE) In this chapter Mathematica will be used to study the stability of the critical points of twin equation linear SODE.
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 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 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 informationThe Higgins-Selkov oscillator
The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being
More 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 informationDynamical systems tutorial. Gregor Schöner, INI, RUB
Dynamical systems tutorial Gregor Schöner, INI, RUB Dynamical systems: Tutorial the word dynamics time-varying measures range of a quantity forces causing/accounting for movement => dynamical systems dynamical
More informationNon-Linear Problems. Almost Linear Systems Consider the system of equations: (1) Y ' = AY + ghyl
Non-Linear Problems Almost Linear Systems Consider the system of equations: (1) Y ' = AY + ghyl (here Y is a matrix). Suppose that the origin is an isolated equilibrium point (that is, there is some circle
More informationLimits and Continuity. 2 lim. x x x 3. lim x. lim. sinq. 5. Find the horizontal asymptote (s) of. Summer Packet AP Calculus BC Page 4
Limits and Continuity t+ 1. lim t - t + 4. lim x x x x + - 9-18 x-. lim x 0 4-x- x 4. sinq lim - q q 5. Find the horizontal asymptote (s) of 7x-18 f ( x) = x+ 8 Summer Packet AP Calculus BC Page 4 6. x
More informationProblem Set Number 5, j/2.036j MIT (Fall 2014)
Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:
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 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 informationCCNY 203 Fall 2011 Final Solutions
CCNY Fall Final Solutions a: Take cross products of the displacement vectors to get a normal to plane P: Cross@8,, < - 8,, -
More informationMath 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv
Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More 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 informationtutorial ii: One-parameter bifurcation analysis of equilibria with matcont
tutorial ii: One-parameter bifurcation analysis of equilibria with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht February 13, 2018 1 This session
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 informationTopic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis
Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability of
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationApril 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.
April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set
More informationCENTER MANIFOLD AND NORMAL FORM THEORIES
3 rd Sperlonga Summer School on Mechanics and Engineering Sciences 3-7 September 013 SPERLONGA CENTER MANIFOLD AND NORMAL FORM THEORIES ANGELO LUONGO 1 THE CENTER MANIFOLD METHOD Existence of an invariant
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More informationDynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting
Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where
More informationBifurcation of Fixed Points
Bifurcation of Fixed Points CDS140B Lecturer: Wang Sang Koon Winter, 2005 1 Introduction ẏ = g(y, λ). where y R n, λ R p. Suppose it has a fixed point at (y 0, λ 0 ), i.e., g(y 0, λ 0 ) = 0. Two Questions:
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 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 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 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 informationChapter 9 Global Nonlinear Techniques
Chapter 9 Global Nonlinear Techniques Consider nonlinear dynamical system 0 Nullcline X 0 = F (X) = B @ f 1 (X) f 2 (X). f n (X) x j nullcline = fx : f j (X) = 0g equilibrium solutions = intersection of
More 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 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 informationCHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION
CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION [Discussion on this chapter is based on our paper entitled Hopf-Bifurcation Ina Two Dimensional Nonlinear Differential Equation,
More informationBINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem
BINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem Michael Priolo Jerry Marsden, Ph.D., Mentor Shane Ross, Graduate Student, Co-Mentor Control and Dynamical Systems 107-81 Caltech,
More informationMathematical analysis of a third-order memristor-based Chua oscillators
Mathematical analysis of a third-order memristor-based Chua oscillators Vanessa Botta, Cristiane Néspoli, Marcelo Messias Depto de Matemática, Estatística e Computação, Faculdade de Ciências e Tecnologia,
More informationLecture 10 Polynomial interpolation
Lecture 10 Polynomial interpolation Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationProblem set 6 Math 207A, Fall 2011 Solutions. 1. A two-dimensional gradient system has the form
Problem set 6 Math 207A, Fall 2011 s 1 A two-dimensional gradient sstem has the form x t = W (x,, x t = W (x, where W (x, is a given function (a If W is a quadratic function W (x, = 1 2 ax2 + bx + 1 2
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 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 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 informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationApplication demonstration. BifTools. Maple Package for Bifurcation Analysis in Dynamical Systems
Application demonstration BifTools Maple Package for Bifurcation Analysis in Dynamical Systems Introduction Milen Borisov, Neli Dimitrova Department of Biomathematics Institute of Mathematics and Informatics
More informationME201/MTH281/ME400/CHE400 Newton's Law of Cooling
ME1/MTH281/ME0/CHE0 Newton's Law of Cooling 1. Introduction This notebook uses Mathematica to solve the problem presented in class, in section 5.1 of the notes. The problem is one of transient heat conduction
More informationOn the complex volume of hyperbolic knots Yoshiyuki Yokota Tokyo Metropolitan University
On the complex volume of hyperbolic knots Yoshiyuki Yokota Tokyo Metropolitan University Let M be a closed, oriented, hyperbolic 3-manifold. Then, the Chern-Simons invariant of M is defined by cs(m) =
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationMATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum
MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum Reconsider the following example from last week: dx dt = x y dy dt = x2 y. We were able to determine many qualitative features
More informationM3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011
M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics
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 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 informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationME3.6 Sheet 2 Answers
ME36 Sheet 2 Answers (i) We have dx/dt Ý6x 2xy Ý 8 fx,y,dy/dt y 2 Ý x 2 gx,y The critical points satisfy g y x and f Ý6x 2xy Ý 8 Substitute y x to get Ý6x 2x 2 Ý 8 x Ý 4x x 4 or Ý If instead we substitute
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
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 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 information1. Taylor Polynomials of Degree 1: Linear Approximation. Reread Example 1.
Math 114, Taylor Polynomials (Section 10.1) Name: Section: Read Section 10.1, focusing on pages 58-59. Take notes in your notebook, making sure to include words and phrases in italics and formulas in blue
More informationand A I L j T O S O L O LOWELL, MICHIGAN, THURSDAY, AUG. 1, 1935
Y D D Y 5 VD D Y D D - ( D K D ( > Kz j K D x j ; K D D K x z ( D K D ( ; ) ( ) V DY - j! ) : x x x j K D K j DY 95 Y-D Y D D j 5 4 V 89 j j j j 89 j 4998 9 j K 8 j V x j K x x 5 x x x j - K 4- -D K 4-
More informationSection 8.2 : Homogeneous Linear Systems
Section 8.2 : Homogeneous Linear Systems Review: Eigenvalues and Eigenvectors Let A be an n n matrix with constant real components a ij. An eigenvector of A is a nonzero n 1 column vector v such that Av
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 informationM.5 Modeling the Effect of Functional Responses
M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationSynopsis of Numerical Linear Algebra
Synopsis of Numerical Linear Algebra Eric de Sturler Department of Mathematics, Virginia Tech sturler@vt.edu http://www.math.vt.edu/people/sturler Iterative Methods for Linear Systems: Basics to Research
More informationPhase Plane Analysis
Phase Plane Analysis Phase plane analysis is one of the most important techniques for studying the behavior of nonlinear systems, since there is usually no analytical solution for a nonlinear system. Background
More informationCALCULUS MATH*2080 SAMPLE FINAL EXAM
CALCULUS MATH*28 SAMPLE FINAL EXAM Sample Final Exam Page of 2 Prof. R.Gentry Print Your Name Student No. SIGNATURE Mark This exam is worth 45% of your final grade. In Part I In Part II In part III In
More informationSolutions of Math 53 Midterm Exam I
Solutions of Math 53 Midterm Exam I Problem 1: (1) [8 points] Draw a direction field for the given differential equation y 0 = t + y. (2) [8 points] Based on the direction field, determine the behavior
More informationAutonomous Mechanical systems of 1 degree of freedom
Autonomous Mechanical systems of degree of freedom = f (, ), = y = y = f (, y) = f (, ) : position y : velocity f : force function Solution : = ( t; c, c), y = y( t; c, c ), ci = ci(, y) Integral of motion
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 information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationMath 312 Lecture Notes Linearization
Math 3 Lecture Notes Linearization Warren Weckesser Department of Mathematics Colgate University 3 March 005 These notes discuss linearization, in which a linear system is used to approximate the behavior
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
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 informationE209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions
E9A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions Michael Vitus Gabe Hoffmann Stanford University Winter 7 Problem 1 The governing equations are: ẋ 1 = x 1 + x 1 x ẋ = x + x 3 Using
More information. CALCULUS AB. Name: Class: Date:
Class: _ Date: _. CALCULUS AB SECTION I, Part A Time- 55 Minutes Number of questions -8 A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAMINATION Directions: Solve each of the following problems, using
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationVideo 6.1 Vijay Kumar and Ani Hsieh
Video 6.1 Vijay Kumar and Ani Hsieh Robo3x-1.6 1 In General Disturbance Input + - Input Controller + + System Output Robo3x-1.6 2 Learning Objectives for this Week State Space Notation Modeling in the
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 informationMA261-A Calculus III 2006 Fall Midterm 2 Solutions 11/8/2006 8:00AM ~9:15AM
MA6-A Calculus III 6 Fall Midterm Solutions /8/6 8:AM ~9:5AM. Find the it xy cos y (x;y)(;) 3x + y, if it exists, or show that the it does not exist. Assume that x. The it becomes (;y)(;) y cos y 3 + y
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More information6.3. Nonlinear Systems of Equations
G. NAGY ODE November,.. Nonlinear Systems of Equations Section Objective(s): Part One: Two-Dimensional Nonlinear Systems. ritical Points and Linearization. The Hartman-Grobman Theorem. Part Two: ompeting
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More information2. Second-order Linear Ordinary Differential Equations
Advanced Engineering Mathematics 2. Second-order Linear ODEs 1 2. Second-order Linear Ordinary Differential Equations 2.1 Homogeneous linear ODEs 2.2 Homogeneous linear ODEs with constant coefficients
More informationMath 216 First Midterm 19 October, 2017
Math 6 First Midterm 9 October, 7 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More information