E209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions
|
|
- Natalie Patterson
- 6 years ago
- Views:
Transcription
1 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 the function V (x) = x 1 +x let s see what we can say about the stability of the origin. Since V (x) is a class-k function, by definition V (x) is a positive-definite and decrescent function. The Lie derivative is: V = (x 1 + x )(x 1) () Therefore, V only when x 1. Therefore, the equilibrium is only locally stable on G r = {x x 1}. (1)
2 Problem : Phase-Locked Loop. We are to analyze the phase-locked loop described by ÿ(t) + (a + bcos y(t))ẏ + csin y(t) = (3) Let x 1 = y and x = ẏ. We need to show that (,) is a stable equilibrium point if a b. The system is then ẋ 1 = x ẋ = (a + bcos x 1 )x csin x 1 Given the Lyapunov function candidate V (x 1,x ) = c(1 cos x 1 ) + x /, we need to start by showing that V (x,t) is positive definite and decrescent on G r. We assume that c > for the hint to be helpful. It is helpful to use the trig identity Thus, 1 cos x 1 = sin x 1 V = csin x 1 + x (5) For simplicity, choose G r = {x R n : x < π}. Thus, V monotonically increases with x from the origin until x 1 = ±π, when sin x 1 = sin π = 1, and the x 1 term is maximum with respect to x 1. To show that the function is decrescent on G r, we look at x ( csin 1 + x c x max, 1 ) (6) where the c condition in the max function arises from sin x 1 x 1 x 1 (7) and the 1 condition in the max function is for the case where x / dominates if c is too small. Therefore, φ( x ) = max( c, 1 ) x is a class-k function that bounds V from above, proving that V is decrescent on G r. To show that V is positive definite on G r, we use the fact that x ( csin 1 + x c min 1 ) π, x when x 1 (8) where the 1 condition in the min function is from the case where x c / dominates if c is too big. The π term is from solving for the coefficient necessary for x 1 to intercect csin x 1 at x 1 = π. That is, when csin π = kπ. The parabola can only intercept the square of the sin at one point between and π, and it does so from below, so it the condition in (8) is proven true. Thus, min ( c, 1 π ) x, another class-k function, bounds V from below on G r, proving that V is positive definite on G r. Now, we need to show that V is positive definite. V = V ẋ 1 + V ẋ x 1 x ( = csin x 1 cos x ) 1 (x ) + (x ) ( (a + bcos x 1 )x csin x 1 ) = cx sin x 1 x (a + bcos x 1) csin x 1 = x (a + bcos x 1) (4)
3 Clearly, x is always. The term (a + bcos x 1) is always, because a b, so a + bcos x 1 b(1 + cos x 1 ). Therefore, V, so the origin is a stable equilibrium point. An alternative method to prove stability, although it goes against the implication of the hint, is to find the Jacobian linearization of the system about the equilibrium. [ ] [ ] 1 1 Df (,) = = (9) bx sin x 1 ccos x 1 (a + bcos x 1 ) c (a + b) This ( leads to a characteristic equation of λ + (a + b)λ + c =. Solving for λ, we see that λ = 1 (a + b) ± ) (a + b) 4c. Thus, the real part of the eigenvalue is always less than zero, by the triangle inequality. There is the rather large exception to this stability claim, a = b =, which is allowed by the constraints, where the roots are purely imaginary. By Hartman-Grobman, we cannot tell the stability of the system from the linearization in this case. We can, however, tell from the Lyapunov analysis that we already did. If a = b =, then V =, so the solution trajectories would follow level curves of V encircling the origin. (,) 3
4 Problem 3. ÿ = u u = k(y) Using V (y,ẏ) = 1 ẏ + y o k(σ)dσ show that the origin is stable in the sense of Lyapunov if k() = and dk dy () >. Let x 1 = y and x = ẏ. V (x) = x ẋ + k(x 1 )ẋ 1 (1) V (x) = x ( k(x 1 )) + k(x 1 )x (11) V (x) = (1) Therefore, V (x). Setting ẋ 1 = ẋ =, yields (k 1 (),) and since k() =, the origin is an equilibrium point. Now we need to show that V (x) is positive-definite and decrescent. Since V (x) is time-invariant, it is decrescent. Since dk dx 1 () >, there is no deadband in k(x 1 ). Therefore, at least locally around (,) we can bound the function V (x) from below. Therefore, V (x) is LPDF, decrescent Lyapunov function. Therefore, the origin is locally stable. (13) 4
5 Problem 4: Control. We are to consider the system ẋ 1 = a(x x 1 ) ẋ = bx 1 x x 1 x 3 + u ẋ 3 = x 1 + x 1 x ax 3 where a > and b > are constants. We are asked if using Lyapunov theory, we can design a control law u(x) which globally stabilizes the origin. Let s choose a positive definite, decrescent candidate Lyapunov function, V (x), and see if we can define u(x) such that V (x). Let s try Then, we find the Lie derivative. V (x) = 1 x x + 1 x 3 (14) V (x) = x 1 ẋ 1 + x ẋ + x 3 ẋ 3 = x 1 a(x x 1 ) + x (bx 1 x x 1 x 3 + u) + x 3 (x 1 + x 1 x ax 3 ) = ax 1 x ax 1 + bx 1 x bx x 1 x x 3 + x u + x 1 x 3 + x 1 x x 3 ax 3 = (a + b)x 1 x + x 1 x 3 ax 1 bx ax 3 + x u In order to eliminate the terms that can be positive in the derivative, we could chose u(x) = x 1 (a + b) x 1x 3 x (15) However, this control law goes to infinity when x. Let s see what the real impact of the x 1 x 3 term is. Perhaps it is cancelled by the quadratic terms. Using the control law u(x) = x 1 (a + b), we would have V (x) = x 1 x 3 ax 1 bx ax 3 (16) Clearly, the x term is less than or equal to zero. The remainder of the terms can be simplified by seeing that we can rewrite x 1 x 3 in a quadratic form. (x 1 cx 3 ) = x 1 + cx 1 x 3 c x 3 x 1 x 3 = 1 c x 1 + c x 3 (x 1 cx 3 ) (17) where c is an arbitrary scaling constant. Substituting (17) back into (16) we obtain V (x) = 1 ( (x 1 x 3 ) a 1 ) ( x 1 a c ) x 3 (18) c Therefore, if a 1 c and a c, then V, and the origin is globally stabilized. To find the value of c that constrains a the least, we solve for c at the point where both inequalities transition from true to false at the same time. That leads to the system of equations a = 1 c a = c 5
6 The solution is c = and a = 1 1. Thus, if a, then the Lie derivative is less than or equal to zero, and the origin is globally stabilized. However, we do not have the option of picking a, therefore this candidate function does not work. Therefore, let s try a different candidate function. V (x) = 1 (αx 1 + βx + δx 3 ) (19) As long as α, β, and δ, the function V (x) is positive definite and decrescent. Taking the Lie derivative yields, V = αx 1 ẋ 1 + βx ẋ + δx 3 ẋ 3 V = aαx 1 βx aδx 3 + (δ β)x 1x x 3 + δx 1 x 3 + (bβ + aα)x 1 x + βu Let s choose β = δ, which will eliminate the term (δ β)x 1 x x 3. Also, let s choose u = 1 β (bβ + aα)x 1, which will eliminate one of the two positive terms. Substituting and simplifying yields, V = aαx 1 βx aβx 3 + βx 1 x 3 Now we need to show that V for some α and β. V can be written in the form where Q = V = x T Qx () aα.5β β.5β aβ (1) We now need to show that Q is positive-semidefinite so that V (x). In order for Q, the diagonal of Q needs to be positive, and the determinant needs to be greater than or equal. Therefore, a >, α >, and β >. The determinant of Q is (a αβ.5β ) () Setting the determinant equal to zero, and solving for β yields β = 8a α (3) Choosing α = 1, yeilds β = 8a. Therefore, with the function V (x) = 1 x 1 + 4a (x + x 3 ) which is positive definite and decrescent, and with u = (b + 1 8a )x 1, we have shown that V (x). Therefore, we have found a control input which makes the origin stable in the sense of Lyapunov. 6
7 Problem 5: RLC circuit with passive nonlinear resistor. The state equations of this RLC circuit are ẋ = y f(x) (4) ẏ = x (5) where x is the current, y is the capacitor voltage, and f(x) is the nonlinear resistor s effect. We can assume that f(x) is a continuous function, because a resistor is a physical device. We are to show that if the resistor is strictly passive, then the equations admit only one stable equilibrium. Let s start by finding the equilibria. ẏ = x = x = (6) ẋ = y f(x) = y = f(x) = f() (7) So the equilibrium is at (,f()). It is unique, because xf(x) = only when x =. Next, we use the stored energy function for the circuit, W = (x + y )/, as the Lyapunov function. It is obviously positive definite and decrescent. Then, V = Ẇ = xẋ + yẏ = x(y f(x) + y( x)) = xf(x) The resistor is passive, so xf(x) = only at x =, and is positive otherwise. Therefore, V, so the minimum of W is a stable equilibrium point. There is only one equilibrium point, so the one equilibrium point is stable. 7
8 Problem 6. We consider the system ẋ 1 = x + ǫx 1 (x 1 + x )sin (x 1 + x ) (8) ẋ = x 1 + ǫx (x 1 + x )sin (x 1 + x ) (9) We are to show that the linearization is inconclusive to determine the stability of the origin, and use the method of Lyapunov to study the stability of the origin between ǫ = 1 and ǫ = 1. First, notice that the origin is clearly an equilibrium point. The Jacobian linearization of the system is ] Df = [ f1 x 1 f 1 x f x 1 f x (3) where f 1 x 1 = ǫx 1(x 1 + x )cos (x 1 + x ) + ǫ(3x 1 + x )sin (x 1 + x ) f 1 x = ǫx 1 x (x 1 + x )cos (x 1 + x ) + ǫx 1x sin(x 1 + x ) 1 f x 1 = ǫx 1 x (x 1 + x )cos (x 1 + x ) + ǫx 1x sin(x 1 + x ) + 1 f x = ǫx (x 1 + x )cos (x 1 + x ) + ǫ(x 1 + 3x )sin (x 1 + x ) Evaluating the Jacobian linearization at the origin, it is simply [ 1 Df (,) = 1 ] (31) Thus, the eigenvalues of the linearized system are simply λ = ±j. Because the eigenvalues are purely imaginary, according to Hartman-Grobman, we cannot determine the stability of this equilibrium. Let us pick the simple and common Lyapunov function, V (x 1,x ) = 1 (x 1 + x ), and see how much information it gives us. The Lie derivative is V = x 1 x 1 + x x ( = x 1 x + ǫx 1 (x 1 + x )sin (x 1 + x )) ( + x x1 + ǫx (x 1 + x )sin (x 1 + x )) = ǫ(x 1 + x ) sin (x 1 + x ) The (x 1 + x ) term is always. Around the origin, the sin(x 1 + x ) term is for x < π. So, let us take G r = {x R n : x < π}. Then, in G r, V > if ǫ >, and V < if ǫ <. Therefore, if ǫ <, the origin is stable, and if ǫ >, the origin is unstable. If ǫ =, then clearly the system is a continuum of circular orbits, as can be seen by looking at the system equations. 8
9 Problem 7: Proving Stability of a Boundary control scheme. The governing equation is m tw(x,t) T xw(x,t) = (3) A vertical force u is applied at the free end of the string, at x = 1. The balance of forces in the vertical direction yields u(t) = T W(x,t) x (33) for all t. We are to show that if the boundary control u(t) = k W(x,t) t is applied, where k >, then Ė(t) for all t. We are given that E(t) = 1 1 ( ) W(x,t) m dx + 1 t (34) 1 ( ) W(x,t) T dx (35) x where the first term is the kinetic energy, and the second term is the strain energy. Taking the derivative, we obtain 1 ( ) ( W ) W 1 ( )( W ) W Ė(t) = m t t dx + T dx (36) x x t Using (3), we see that Thus, ( 1 Ė(t) = T ( W t W t )( W x = T W m x (37) ) dx + 1 ( )( W ) ) W dx x x t Now, we need to simlify this, so we use integration by parts, as given by the hint. To do so, the equations need to be of the form b a udv = uv b a b a vdu. We apply this form to the second integral above, by absorbing the second order derivatives into the differential term, and taking u = W x and v = W t. Then, ( 1 ( ) ( ) W W 1 ( ) ( )) W W Ė(t) = T d + d t x x t ( 1 ( ) ( ) W W = T d + W W 1 1 ( ) ( ) ) W W t x x t d t x = T W W 1 x t = T W W x t T W W x t The first line was found by using W x x = W d ( ) W x x = W dx for the first integral, and similar x calculus for the second integral. The integration limits throughout the problem are from x = to x = 1. The second line was found by applying integration by parts to the second integral in line 1. The third line was found by realizing that the two integrals in line two are equal and have opposite signs. 9 x= (38)
10 We use W(,t) = to see that W t x= =. In other words, the left end of the string is fixed in place. This eliminates the second term in the equation above, yielding ( Ė(t) = T W ) W x t Ė(t) = u(t) W t where the last step used (33). Now, if we define u(t) as in (34), with k >, we have ( ) W Ė(t) = k t (39) t 1
Solution 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 information2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ
1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory
MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A
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 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 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 informationE209A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions
E09A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions Michael Vitus Stanford University Winter 007 Problem : Planar phase portraits. Part a Figure : Problem a This phase portrait is correct.
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 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 7: Passivity
Nonlinear Control Lecture 7: Passivity Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 7 1/26 Passivity
More informationNonlinear systems. Lyapunov stability theory. G. Ferrari Trecate
Nonlinear systems Lyapunov stability theory G. Ferrari Trecate Dipartimento di Ingegneria Industriale e dell Informazione Università degli Studi di Pavia Advanced automation and control Ferrari Trecate
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 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 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 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 informationNonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems
Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin
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 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 informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
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 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 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 informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
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 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 informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
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 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 informationNonlinear differential equations - phase plane analysis
Nonlinear differential equations - phase plane analysis We consider the general first order differential equation for y(x Revision Q(x, y f(x, y dx P (x, y. ( Curves in the (x, y-plane which satisfy this
More 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 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 informationKing Fahd University of Petroleum and Minerals Prep-Year Math Program Math Term 161 Recitation (R1, R2)
Math 001 - Term 161 Recitation (R1, R) Question 1: How many rational and irrational numbers are possible between 0 and 1? (a) 1 (b) Finite (c) 0 (d) Infinite (e) Question : A will contain how many elements
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationSample Solutions of Assignment 10 for MAT3270B
Sample Solutions of Assignment 1 for MAT327B 1. For the following ODEs, (a) determine all critical points; (b) find the corresponding linear system near each critical point; (c) find the eigenvalues of
More 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 informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
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 information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
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 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 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 information16.30/31, Fall 2010 Recitation # 13
16.30/31, Fall 2010 Recitation # 13 Brandon Luders December 6, 2010 In this recitation, we tie the ideas of Lyapunov stability analysis (LSA) back to previous ways we have demonstrated stability - but
More informationMathematics (MEI) Advanced Subsidiary GCE Core 1 (4751) May 2010
Link to past paper on OCR website: http://www.mei.org.uk/files/papers/c110ju_ergh.pdf These solutions are for your personal use only. DO NOT photocopy or pass on to third parties. If you are a school or
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More information1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),
1.7. Stability and attractors. Consider the autonomous differential equation (7.1) ẋ = f(x), where f C r (lr d, lr d ), r 1. For notation, for any x lr d, c lr, we let B(x, c) = { ξ lr d : ξ x < c }. Suppose
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationOutput Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x)
Output Feedback and State Feedback EL2620 Nonlinear Control Lecture 10 Exact feedback linearization Input-output linearization Lyapunov-based control design methods ẋ = f(x,u) y = h(x) Output feedback:
More informationChapter 8B - Trigonometric Functions (the first part)
Fry Texas A&M University! Spring 2016! Math 150 Notes! Section 8B-I! Page 79 Chapter 8B - Trigonometric Functions (the first part) Recall from geometry that if 2 corresponding triangles have 2 angles of
More informationLyapunov Stability Analysis: Open Loop
Copyright F.L. Lewis 008 All rights reserved Updated: hursday, August 8, 008 Lyapunov Stability Analysis: Open Loop We know that the stability of linear time-invariant (LI) dynamical systems can be determined
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 informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationConsider an ideal pendulum as shown below. l θ is the angular acceleration θ is the angular velocity
1 Second Order Ordinary Differential Equations 1.1 The harmonic oscillator Consider an ideal pendulum as shown below. θ l Fr mg l θ is the angular acceleration θ is the angular velocity A point mass m
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 informationInvariant Manifolds of Dynamical Systems and an application to Space Exploration
Invariant Manifolds of Dynamical Systems and an application to Space Exploration Mateo Wirth January 13, 2014 1 Abstract In this paper we go over the basics of stable and unstable manifolds associated
More informationMEM 255 Introduction to Control Systems: Modeling & analyzing systems
MEM 55 Introduction to Control Systems: Modeling & analyzing systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline The Pendulum Micro-machined capacitive accelerometer
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 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 informationCalculus I Review Solutions
Calculus I Review Solutions. Compare and contrast the three Value Theorems of the course. When you would typically use each. The three value theorems are the Intermediate, Mean and Extreme value theorems.
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 informationThe Derivative. Appendix B. B.1 The Derivative of f. Mappings from IR to IR
Appendix B The Derivative B.1 The Derivative of f In this chapter, we give a short summary of the derivative. Specifically, we want to compare/contrast how the derivative appears for functions whose domain
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 informationI ml. g l. sin. l 0,,2,...,n A. ( t vs. ), and ( vs. ). d dt. sin l
Advanced Control Theory Homework #3 Student No.597 Name : Jinseong Kim. The simple pendulum dynamics is given by g sin L A. Derive above dynamic equation from the free body diagram. B. Find the equilibrium
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationEG4321/EG7040. Nonlinear Control. Dr. Matt Turner
EG4321/EG7040 Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear [System Analysis] and Control Dr. Matt
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 informationMathematics (MEI) Advanced Subsidiary GCE Core 1 (4751) June 2010
Link to past paper on OCR website: www.ocr.org.uk The above link takes you to OCR s website. From there you click QUALIFICATIONS, QUALIFICATIONS BY TYPE, AS/A LEVEL GCE, MATHEMATICS (MEI), VIEW ALL DOCUMENTS,
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationLinearization at equilibrium points
TMA4165 2007 Linearization at equilibrium points Harald Hanche-Olsen hanche@math.ntnu.no This note is about the behaviour of a nonlinear autonomous system ẋ = f (x) (where x(t) R n ) near an equilibrium
More informationComputer Problems for Methods of Solving Ordinary Differential Equations
Computer Problems for Methods of Solving Ordinary Differential Equations 1. Have a computer make a phase portrait for the system dx/dt = x + y, dy/dt = 2y. Clearly indicate critical points and separatrices.
More informationSOLUTIONS FOR PROBLEMS 1-30
. Answer: 5 Evaluate x x + 9 for x SOLUTIONS FOR PROBLEMS - 0 When substituting x in x be sure to do the exponent before the multiplication by to get (). + 9 5 + When multiplying ( ) so that ( 7) ( ).
More informationUpdated: January 16, 2016 Calculus II 7.4. Math 230. Calculus II. Brian Veitch Fall 2015 Northern Illinois University
Math 30 Calculus II Brian Veitch Fall 015 Northern Illinois University Integration of Rational Functions by Partial Fractions From algebra, we learned how to find common denominators so we can do something
More informationStability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming
arxiv:math/0603313v1 [math.oc 13 Mar 2006 Stability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming Erin M. Aylward 1 Pablo A. Parrilo 1 Jean-Jacques E. Slotine
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 informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
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 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 informationOrdinary Differential Equations
Ordinary Differential Equations (MA102 Mathematics II) Shyamashree Upadhyay IIT Guwahati Shyamashree Upadhyay ( IIT Guwahati ) Ordinary Differential Equations 1 / 10 Undetermined coefficients-annihilator
More information2t t dt.. So the distance is (t2 +6) 3/2
Math 8, Solutions to Review for the Final Exam Question : The distance is 5 t t + dt To work that out, integrate by parts with u t +, so that t dt du The integral is t t + dt u du u 3/ (t +) 3/ So the
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 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 informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
More informationSystems of Equations and Inequalities. College Algebra
Systems of Equations and Inequalities College Algebra System of Linear Equations There are three types of systems of linear equations in two variables, and three types of solutions. 1. An independent system
More informationRoots and Coefficients Polynomials Preliminary Maths Extension 1
Preliminary Maths Extension Question If, and are the roots of x 5x x 0, find the following. (d) (e) Question If p, q and r are the roots of x x x 4 0, evaluate the following. pq r pq qr rp p q q r r p
More informationCalculus of Variations
ECE 68 Midterm Exam Solution April 1, 8 1 Calculus of Variations This exam is open book and open notes You may consult additional references You may even discuss the problems (with anyone), but you must
More informationPartial Fraction Decomposition Honors Precalculus Mr. Velazquez Rm. 254
Partial Fraction Decomposition Honors Precalculus Mr. Velazquez Rm. 254 Adding and Subtracting Rational Expressions Recall that we can use multiplication and common denominators to write a sum or difference
More informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
More informationPractice Exam 1 Solutions
Practice Exam 1 Solutions 1a. Let S be the region bounded by y = x 3, y = 1, and x. Find the area of S. What is the volume of the solid obtained by rotating S about the line y = 1? Area A = Volume 1 1
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 informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationCalculus with business applications, Lehigh U, Lecture 03 notes Summer
Calculus with business applications, Lehigh U, Lecture 03 notes Summer 01 1 Lines and quadratics 1. Polynomials, functions in which each term is a positive integer power of the independent variable include
More informationLyapunov stability ORDINARY DIFFERENTIAL EQUATIONS
Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS An ordinary differential equation is a mathematical model of a continuous state continuous time system: X = < n state space f: < n! < n vector field (assigns
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 information1MA6 Partial Differentiation and Multiple Integrals: I
1MA6/1 1MA6 Partial Differentiation and Multiple Integrals: I Dr D W Murray Michaelmas Term 1994 1. Total differential. (a) State the conditions for the expression P (x, y)dx+q(x, y)dy to be the perfect
More informationMarch Algebra 2 Question 1. March Algebra 2 Question 1
March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question
More informationSTABILITY ANALYSIS OF DYNAMIC SYSTEMS
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY STABILITY ANALYSIS OF DYNAMIC SYSTEMS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e
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 information