Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
|
|
- Noel McDonald
- 5 years ago
- Views:
Transcription
1 Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points Linearisation Limit cycles Bifurcations Stability Exercises Matlab Based Exercises for EEE Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /0
2 Nonlinear Dynamics. Linear Combination Until now we have seen only linear (autonomous) systems which unfortunately can only fully describe a small class of real systems. In general a nonlinear system is given by: x t f x t,u t () with x,u, f n q n The main property of a linear system (as we have seen) is that if you have solutions x and x then also their linear combination is a solution. For example, if we have: xt 3 x and x 3t e and x 3t 7e then x x x e e t t 4 is also a solution as: d 0e 4e 3t 3t dt 30e 4e e 3t 3t 3t 3 0e 4e 30e 4e e 3t 3t 3t 3t 3t However this is not the case in nonlinear systems: x t x and a solution: x t C as x and t C x A 3 is not a solution: t C x t C t C but Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /0
3 3 x and t C x t C t C t C x. Equilibrium points Before we see more differences between linear and nonlinear systems it is important to understand the concept of the equilibrium or rest or singular point which is defined as the point where xt 0 () For first order linear systems we have that 0 0 x t A t x t B t u t x t A t B t u t (3) For autonomous systems: 0 state spaces were drawn with respect to the origin: x t and for that reason in Chapter 3 all the Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 3/0
4 One of the biggest differences between linear and nonlinear systems is that a nonlinear system may have more than one s. Example.: Find the s of The s are: x ' x x x ' x x x' x x 0 x x x x 0 0 x' x x x x x x,x, x,x, So we have Eps at (,) and at (-,-). Example.: Find the s of The s are: x ' x x 3x x 0x x ' x x 4x x x 3xx 0x 0 x x 4x 0 From the nd expression we have that x 0 or x 4 From the st expression we have for x 0 0 x i.e. x,x, and for 00 xx 4 we have that either x 0 which implies that x 0 (as before) or x 3x 0 0 which is a nd order polynomial and can easily be solved as: syms x x, f=x-x, f=x^+x^-, [xs xs]=solve(f,f) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 4/0
5 37 x 5 x,x 5, 4 49 x, x x,x, 4 So the system has the following s: x,x,, x,x,, x,x, In linear systems we were using the term stable/unstable system but effectively we will referring to the equilibrium point. Now that we have multiple equilibriums we need to determine the stability of each. It has to be stated here that for nonlinear systems it does not make sense to talk about stable/unstable systems, but for stable/unstable s as it is possible to have systems with multiple s and some of them to be stable and other to be unstable. Also, the stability of linear systems was determined by the eigenvalues of the state matrix. Hence, now we need a state matrix for each fixed point. 3. Linearisation In order to determine if each of the s of a nonlinear system is stable or not we have to take a local picture of each. This local picture is called a linearization. I.e. we will describe each in a neighbourhood around it. The main tool for that is the Taylor series. Remember that the Taylor series around a point x 0 is defined as: 3 x x x x x x 3 f x 0 f x 0 f x 0 f x f x x! x! x 3! xx xx xx Notice how the command solve returns the results. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 5/0
6 Example.3: Approximate the function y=sin(x) with respect to x around x=0, using several expansion orders. The TS of sin(x) around x=0 is: sin x x x x 3! 5! 7! 9! x x Ox Now let s try to approximate sin(x): Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 6/0
7 3 So by increasing the order of the Taylor Series we get a better approximation. But let s get back to the first term y=x, obviously it is not a good approximation but if we focus close to the origin then: 3 x=-*pi:0.0:*pi; plot(x,sin(x)), hold on, plot(x,x), plot(x,x-x.^3/6), plot(x,x-x.^3/6+x.^5/0), plot(x,xx.^3/6+x.^5/0-x.^7/5040), plot(x,x-x.^3/6+x.^5/0-x.^7/5040+x.^9/36880) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 7/0
8 So we see that close to the origin effectively the two traces coincide, i.e. our approximation is good close to the point of expansion. This is exactly how we will study a nonlinear system and the stability of their s. By approximating them locally using a Taylor series expansion and by keeping only the linear term: Assume that we have a nonlinear system x x f x f x x 3 f3 x x t f x and a at x=x, then the TS expansion is: f x t xt f x xt x HOT xt xx By keeping only the linear term we have: f x t xt f x xt x xt xx (4) As 0 f x : f x t xt xt x xt xx dx Obviously we also have that: 0 and hence: dt Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 8/0
9 dt xt d x t x f x t xx x t x We can define now a new variable: x t x t x and we have that: dx dt f x t x t xx x Also, f x t xt xx f x t f x t f x t x t x t xn t f x t x t fn xt fn x t x t xn t xx A (this matrix is called the Jacobian of f) and hence: dx dt Ax (5) Which is a linear equation that we know how to solve and how to determine its stability. Similarly if xt f x,u then: dx f x t,u t f x t,u t x dt x t u t xx,u u xx,u u u (6) with u u u Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 9/0
10 dx which will give us: Ax dt Bu which is the form for a linear state space model. Example.4: Find the s and determine the stability of its s x x For the s we set x x x 0 0. Now the Jacobian is: x f x A x x x And when this is evaluated at the s: A Hence, the state space models are: x x for x x x for x This implies that the eigenvalue of the st is - and hence this is stable, while the eigenvalue of the nd is + and hence that is unstable. So we expect the following D state space: - + It is interesting to note here, that if we start at the left of the -, I will diverge to -. While otherwise I will converge to +: Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 0/0
11 4 Example.5: Find the s and determine their stability for x ' x x x ' x x We have seen that the s are: x,x, x,x, So we have: x x x x x ' x x x x A x ' x x x x x x x x x x Now the Jacobian evaluated at the s: 4 hold on, for x0=-:0.5:.5, sim('non_d.slx'), plot(x(:,),x(:,)), end Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /0
12 3 7i A eigs which implies that (,) is an unstable focus. A 7 eigs which implies that (,) is a saddle. Now in order to understand what happens we have to draw the eigenvectors of the saddle (the eigenvectors of the focus do not give us anything as we have seen). The saddle has eigenvalues at.56 and -.56 with the corresponding eigenvectors being: e and e 0. 8 Hence, the state space is: 5 5 syms x x, f=x-x, f=x^+x^-, [xs xs]=solve(f,f), J=jacobian([f f], [x x]); A=subs(J,{x,x},{xs() xs()}), A=subs(J,{x,x},{xs() xs()}), eigsa=eval(eig(a)) eigsa=eval(eig(a)) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /0
13 Chapter 4 4. Limit cycles Apart from the previous difference of having multiple equilibria in nonlinear systems, another feature of these systems is that it is possible to have a persistence periodic motion, a limit cycle. Assume the Van der Pol system x x x x 0. The time and the state space responses are: Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 3/0
14 5. Bifurcations Another difference between linear and nonlinear systems is the phenomenon of bifurcation, which (very briefly) is the creation (or destruction) of equilibria by changing some (one in our case) parameters of the system. This can be the inductance in an RL circuit, the mass of a mass-spring system, the force applied on a bridge, the temperature in a water tank So now, we study systems of the form x f x, p, where p is a vector that contains the system s parameters that may change and effectively the algebraic equation f x, p 0 has different number of solutions for different parameters in p. Example.6: Assume 6 the system x r x for various values of r. The equilibria are: x r x x r So if r>0 there are no equilibrium points 6 Example taken from Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 4/0
15 While if r<0: x r x x r i.e. equilibria which for various values of r are: hold on for r=-0:0.:0 plot(r,sqrt(-r),'k.') plot(r,-sqrt(-r),'k.') end The linearised system closed to these points is: r rx 0 x x x t e x rx x x t0 e r The response for r=-: Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 5/0
16 6. Stability In linear systems we called a system stable if the orbit was not diverging to infinity (exponentially or with oscillations). In nonlinear systems, we have multiple equilibria and we must check the stability locally around each equilibrium point. But we need a more general definition that will allow us to focus on global stability. Assume that we have a general nonlinear ODE and a specific solution x(t). Obviously as this will depend on the initial conditions we can denote it as x t,t 0. Now we want to test the stability of the general orbit φ (note that the can be one such orbit). To do that, we add a perturbation (say at t=t 0 ) and then we observe the perturbed orbit. If it stays close to the nominal then it is called stable. Effectively this means that for each required Δ we can always find a suitable δ. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 6/0
17 x t,t 0 0 x t,t 0 0 x t,t t 0 0 As we are focused on equilibria we have the following definition: x t,t t 0 0 An equilibrium point x is stable if for each R>0 there exists an r>0 such that E E, then x x t r 0 x x t R E One of the most important theorems that are used to study the stability of nonlinear systems is called: Lyapunov Stability Theorem Let function V of the states x, such as: Vx 0 Vx 0 x lim V x then the equilibrium point is globally asymptotic stable. The difficulty here is how to find the function V, so for the purposes of this module, this function V will be given. Example.7: Determine the stability of the origin of d x x y xy y dt x y x y V x, y x y as a candidate, by using Lyapunov Function and by using local linearization: Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 7/0
18 V x, y x y 0, V x y xx yy x x y xy y x y x y x xy x y 4xy y x y x xy y 4x y x y 4x y 0 x y xy y xy F A x y x y xy x r A 0 r 0 r 4 r r 3 0 r complex with negative real part, so stable 7 7 syms x x, f=-x+x-x*x^; f=-*x-x-x^*x; [xs xs]=eval(solve(f,f)); J=jacobian([f f], [x x]); A=subs(J,{x,x},{0 0}); eigsa=eval(eig(a)) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 8/0
19 7. Exercises ) For the system d u u v dt v, 4 u v i. Find the fixed points and classify them according to their stability. ii. Make a draft sketch showing the fixed points and a few trajectories. ) For the system point at the origin, using: i. V x, y x y d x x y xy y, prove that the system has a stable equilibrium dt x y x y. ii. Local linearization. 3) For the Lorenz equations: x a y x, y bx y xz, z xy cz where abc,, 0 are constants: i. Find the equilibrium point(s). ii. Show that if b, then the origin is the only equilibrium point, and that there three equilibrium points if b. Discuss the stability of the origin for b< and for b>. It is given that for b< we have that b> we have that a a 4ab. 4) Determine the stability of the equilibrium points of Sketch a draft bifurcation diagram. 5) Assume an input matrix a a 4ab while for dx dt B 3 and four state matrices: A, A, A 3, A rx x for r<0 and r>0. i. Find the equilibrium points of each system. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 9/0
20 ii. iii. Find the eigenvalues and hence describe the response of each system starting from x0 0 T a) and b) from the equilibrium point. Create the Simulink model for the above systems and hence crosscheck your previous answers. 6) For a given function f x cos x i. Find the st, nd, 3rd and 4th order TS expansion around x=0. ii. Using Matlab crosscheck that the fourth order is indeed a better approximation that the st order one. 7) A nonlinear system is given by i. Find the equilibrium points. ii. iii. iv. Determine the stability of each point. Create a draft sketch of the state space. x x y x xy, y x y Create the Simulink model of system and hence crosscheck your answers. v. For each, simulate the linearized approximation of the nonlinear system and compare the responses. 8) For the system x r x : i. Describe its response based on the values of r. ii. Create a Simulink model that will crosscheck your answer. 9) Repeat question 8 for x a x, y y 8. Matlab Based Exercises for EEE8086 Crosscheck your analysis in section 7 using Matlab. This must include mfiles and Simulink models. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 0/0
MCE693/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 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 informationSection 5. Graphing Systems
Section 5. Graphing Systems 5A. The Phase Plane 5A-1. Find the critical points of each of the following non-linear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A-2. Write each
More 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 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 informationStability lectures. Stability of Linear Systems. Stability of Linear Systems. Stability of Continuous Systems. EECE 571M/491M, Spring 2008 Lecture 5
EECE 571M/491M, Spring 2008 Lecture 5 Stability of Continuous Systems http://courses.ece.ubc.ca/491m moishi@ece.ubc.ca Dr. Meeko Oishi Electrical and Computer Engineering University of British Columbia,
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 informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
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 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 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 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 informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
More informationModule 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs
Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs 3 Analytical Solutions of Linear ODE-IVPs Before developing numerical
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 informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 533 Homeworks Spring 07 Updated: Saturday, April 08, 07 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. Some homework assignments refer to the textbooks: Slotine
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 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 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 informationHandout 2: Invariant Sets and Stability
Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More 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 informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationMath 273 (51) - Final
Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More 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 informationQuestion: Total. Points:
MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make
More information1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy
.. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution
More informationENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)
ENGI 940 4.06 - Linear Approximation () Page 4. 4.06 Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with
More informationSolution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)
Math 380 Practice Final Solutions This is longer than the actual exam, which will be 8 to 0 questions (some might be multiple choice). You are allowed up to two sheets of notes (both sides) and a calculator,
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 informationQualitative Analysis of Tumor-Immune ODE System
of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 15, 2002 This work was supported in part by a grant from the W.M. Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1. Simplified System
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 informationNAME: MA Sample Final Exam. Record all your answers on the answer sheet provided. The answer sheet is the only thing that will be graded.
NAME: MA 300 Sample Final Exam PUID: INSTRUCTIONS There are 5 problems on 4 pages. Record all your answers on the answer sheet provided. The answer sheet is the only thing that will be graded. No books
More informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
More informationSystem Control Engineering 0
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 08, 2018, at 08 30 12 30 Johanneberg Kristian
More informationMath 216 First Midterm 18 October, 2018
Math 16 First Midterm 18 October, 018 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
More 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 Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
More informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x
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 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 informationMath 216 Final Exam 24 April, 2017
Math 216 Final Exam 24 April, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 5 Homewors Fall 04 Updated: Sunday, October 6, 04 Some homewor assignments refer to the textboos: Slotine and Li, etc. For full credit, show all wor. Some problems require hand calculations. In those
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More 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 informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More informationALGEBRAIC GEOMETRY HOMEWORK 3
ALGEBRAIC GEOMETRY HOMEWORK 3 (1) Consider the curve Y 2 = X 2 (X + 1). (a) Sketch the curve. (b) Determine the singular point P on C. (c) For all lines through P, determine the intersection multiplicity
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More 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 informationDef. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)
Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November
More 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 2009 Farzaneh Abdollahi Nonlinear Control Lecture 2 1/68
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 informationWorkshop on Theoretical Ecology and Global Change March 2009
2022-3 Workshop on Theoretical Ecology and Global Change 2-18 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological
More informationStability Analysis for ODEs
Stability Analysis for ODEs Marc R Roussel September 13, 2005 1 Linear stability analysis Equilibria are not always stable Since stable and unstable equilibria play quite different roles in the dynamics
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
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 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 informationVideo 8.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar
Video 8.1 Vijay Kumar 1 Definitions State State equations Equilibrium 2 Stability Stable Unstable Neutrally (Critically) Stable 3 Stability Translate the origin to x e x(t) =0 is stable (Lyapunov stable)
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
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 informationNonconstant Coefficients
Chapter 7 Nonconstant Coefficients We return to second-order linear ODEs, but with nonconstant coefficients. That is, we consider (7.1) y + p(t)y + q(t)y = 0, with not both p(t) and q(t) constant. The
More informationSolutions to Math 53 First Exam April 20, 2010
Solutions to Math 53 First Exam April 0, 00. (5 points) Match the direction fields below with their differential equations. Also indicate which two equations do not have matches. No justification is necessary.
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 informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More information1 The relation between a second order linear ode and a system of two rst order linear odes
Math 1280 Spring, 2010 1 The relation between a second order linear ode and a system of two rst order linear odes In Chapter 3 of the text you learn to solve some second order linear ode's, such as x 00
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
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 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 informationJune 2011 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations
June 20 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations The topics covered in this exam can be found in An introduction to differential equations
More informationMATH 415, WEEKS 14 & 15: 1 Recurrence Relations / Difference Equations
MATH 415, WEEKS 14 & 15: Recurrence Relations / Difference Equations 1 Recurrence Relations / Difference Equations In many applications, the systems are updated in discrete jumps rather than continuous
More informationCalculus 2502A - Advanced Calculus I Fall : Local minima and maxima
Calculus 50A - Advanced Calculus I Fall 014 14.7: Local minima and maxima Martin Frankland November 17, 014 In these notes, we discuss the problem of finding the local minima and maxima of a function.
More informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More information4 Partial Differentiation
4 Partial Differentiation Many equations in engineering, physics and mathematics tie together more than two variables. For example Ohm s Law (V = IR) and the equation for an ideal gas, PV = nrt, which
More 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 informationIntermediate Differential Equations. Stability and Bifurcation II. John A. Burns
Intermediate Differential Equations Stability and Bifurcation II John A. Burns Center for Optimal Design And Control Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute and
More informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More informationNonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points
Nonlinear Control Lecture # 2 Stability of Equilibrium Points Basic Concepts ẋ = f(x) f is locally Lipschitz over a domain D R n Suppose x D is an equilibrium point; that is, f( x) = 0 Characterize and
More informationMath 341 Fall 2006 Final Exam December 12th, Name:
Math 341 Fall 2006 Final Exam December 12th, 2006 Name: You may use a calculator, your note card and something to write with. You must attach your notecard to the exam when you turn it in. You cannot use
More informationChapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis
Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs... 3. Second Order ODEs... 7 3. General Maerial...
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov
More informationLecture - 11 Bendixson and Poincare Bendixson Criteria Van der Pol Oscillator
Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 11 Bendixson and Poincare Bendixson Criteria
More informationGeorgia Institute of Technology Nonlinear Controls Theory Primer ME 6402
Georgia Institute of Technology Nonlinear Controls Theory Primer ME 640 Ajeya Karajgikar April 6, 011 Definition Stability (Lyapunov): The equilibrium state x = 0 is said to be stable if, for any R > 0,
More informationMath 216 Final Exam 14 December, 2017
Math 216 Final Exam 14 December, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More 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 informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More information+ i. cos(t) + 2 sin(t) + c 2.
MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,
More informationChapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More 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 information