Dynamical Systems & Lyapunov Stability
|
|
- Erica Booker
- 5 years ago
- Views:
Transcription
1 Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University
2 Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence on parameters Invariant sets, nonwandering sets, limit sets Lyapunov Stability Autonomous systems Basic stability theorems Stable, unstable & center manifolds Control Lyapunov function
3 Basics of Nonlinear ODE s
4 Dynamical Systems d x ( t ) ( ) n = f x t, t, x R, t R non-autonomous dt d x ( ) n = f x t, R, R dt autonomous A solution on a time interval t, is a function [ t t ] n xt :[ t, t] R that satisfies the ode.
5 Vector Fields and Flow n We can visualize an individual solution as a graph xt : t R. For autonomous systems it is convenient to think of ( x) n as a vector field on R - f( x) assigns a vector to each point in R. As t varies, a solution xt traces a path through R n ( x) tangent to the field f. These curves are often called trajectories or orbits. The collection of all trajectories in R of the vector field f( x). n is called the flow f n
6 Auto at Constant Speed x V s = V V δ d dt ω = β f ( ωβ,, V, δ ) Body Frame y β u Y v θ m, J a b F l Space Frame X F r Notice the three equilibria.
7 Van der Pol x ( x = ).8 x x x x
8 Damped Pendulum x x = x x / sin x
9 Lipschitz Condition The existence and uniqueness of solutions depend on properties of the function f. In many applications f( xt, ) has continuous derivatives in x. We relax this - we require that f is Lipschitz in x. Def : f : R R is locally Lipschitz on an open subset D R if each point x D has a neighborhood U such that for some constant Note: n n n f x f x L x x C always are. L and all x U (continuous) functions need not be Lipschitz, C functions
10 The Lipschitz Condition A Lipschitz continuous function is limited in how fast it can change, A line joining any two points on the graph of this function will never have a slope steeper than its Lipschitz constant L, The mean value theorem can be used to prove that any differentiable function with bounded derivative is Lipschitz continuous, with the Lipschitz constant being the largest magnitude of the derivative.
11 Examples: Lipschitz =, [ 3, 7 ] L = 4 f x x x L = y f x = x +, x R x L = Fails 3..5 =.5, [ 5,5] f x x x =, [, 4] f x x x
12 Local Existence & Uniqueness Proposition (Local Existence and Uniqueness) Let f( xt, ) be piece-wise continuous in t and satisfy the Lipschitz condition (, ) (, ) f xt f yt L x y { n } for all xy, B= x R x x < r and all t [ t, t]. Then r there exists δ > such that the differential equation with initial condition (, ), x= f xt x t = x B r has a unique solution over [ t, t + δ ].
13 The Flow of a Vector Field x = f( x), x t = x x x, t this notation indicates 'the solution of the ode that passes through x at t = ' ( xt) More generally, let Ψ, denote the solution that passes through at. The function : n n x t = Ψ R R R satisfies Ψ ( xt, ) t ( (, )), (,) = f Ψ xt Ψ x = x Ψ is called the flow or flow function of the vector field f
14 Example: Flow of a Linear Vector Field ( xt, ) Ψ = = Ψ(, ) Ψ (, ) = t Example: At x Ax A xt xt e x + x cos t x sin t 3 x R, A=, Ψ ( xt, ) = xcos t xsin t t e x x x x.5
15 Invariant Set A set of points S R n is invariant with respect to the vector field f if trajectories beginning in S remain in S both forward and backward in time. Examples of invariant sets: any entire trajectory (equilibrium points, limit cycles) collections of entire trajectories
16 Example: Invariant Set x x x d x = x dt x 3 x x3 each of the three trajectories shown are invariant sets the x -x plane is an invariant set x.5
17 Limit Points & Sets A point ( t p) is called an ω-limit point of the trajectory Ψ, if there exists a sequence of time values t + such that q lim Ψ t, p = q t k ( t p) is said to be an α-limit point of Ψ, if there exists a sequence t q R k n of time values lim Ψ t, p = q k k t k such that The set of all ω-limit points of the trajectory through p is the ω-limit set, and the set of all α-limit points is the α-limit set. For an example see invariant set k
18 Introduction to Lyapunov Stability Analysis
19 Lyapunov Stability ε,, n x = f x f = f : D R locally Lipschitz The origin is δ a stable equilibrium point if for each ε >, there is a δ ε > such that x < x t < ε t > unstable if it is not stable, and asymptotically stable if δ can be chosen such that x δ x( t) < lim t δ x
20 Two Simple Results The origin is asymptotically stable only if it is isolated. The origin of a linear system x = Ax At is stable if and only if e N < t > It is asymptotically stable if and only if, in addition At e, t
21 Example: Non-isolated Equilibria x x = x x x x All points on the x axis are equilibrium points x x - -4
22 Positive Definite Functions n A function V : R R is said to be negative (semi-) definite if V ( x) radially unbounded if V ( x) as x T T V ( x) = x Qx Q = Q positive definite if V = and V x >, x positive semi-definite if V = and V x, x is positive (semi-) definite For a quadratic form:, the following are equivalent V x is positive definite the eigenvalues of Q are positive the principal minors of Q are positive
23 Lyapunov Stability Theorem V x is called a Lyapunov function relative to the flow of x = f x respect to the flow: if it is positive definite and nonincreasing with V ( x) V =, V x > for x V = f x x Theorem : If there exists a Lyapunov function on some neighborhood D is negative definite on of the origin, then the origin is stable. If D then it is asymptotically stable. D V = α V = β < α V = γ < β V
24 Example: Linear System T T V x = x Qx, Q = Q > where x = Ax T T T T T V = x Qx + x Qx = x QA + A Q x = x Px T QA A Q P + = Lyapunov Equation At P > Ais stable (Hurwitz) : e bounded t > So we can specify Q, compute P and test P. Or, specify P and solve Lyapunov equation for Q and test Q.
25 Example: Rotating Rigid Body xyz,, body axes; ω, ω, ω angular velocities in body coord's; x y z diag I, I, I I I I > inertia matrix x y z x y z Iz Iy ωx = ωω z y = aωω z y I x Ix I z ωy = ωω x z = bωω x z note abc,, > I y Iy Ix ωz = ωyωx = cωω y x I z z y x
26 Rigid Body, Cont d ( ωx ωy ωz) Note that a state,, is an equilibrium point if any two of the angular velocity components are zero, i.e., the ω, ω, ω axes are all equilibrium points. x y z Consider a point ω,,. Shift ω ω + ω. ω = aωω x z y b ω = ω + ω ω y x x z ω = c ω + ω ω z x x y x x x x z y x ω x ω =
27 Rigid Body, Cont d Energy does not work for ω. Obvious? So, how do we find Lyapunov function? We want V V,, =, x ( ωx, ωy, ωz) if ( ωx, ωy, ωz) (,, ) and ( ωx, ωy, ωz) > D V Lets look at all functions that satisfy V =,i.e., that satisfy the pde: V V V aωzωy + ( b( ωx + ωx) ωz) + c( ωx + ωx) ωy = ω ω ω x y z All solutions take the form: bω + bωω + aω cω + cωω + aω f, a a x x x y x x x z 8bcω x V ( ωx, ωy, ωz) = ca + bb + ( ca bb) = ω... x + cωy + bωz + h o t a
28 Rigid Body, Cont d Clearly, V =, V > on a neighborhood D of V = spin about x-axis is stability This is one approach to finding candidate Lyapunov functions The first order PDE usually has many solutions The method is connected to traditional first integral methods to the study of stability in mechanics Same method can be used to prove stability for spin about z-axis, but spin about y-axis is unstable why?
29 LaSalle Invariance Theorem Theorem : { n x R V x < c} n Suppose : is and let c denote a component of the region Suppose Ω is bounded and within Ω, V x. c V R R C Let E be the set of points within Ω where V x =, Let M be the largest invariant set within E. every sol'n beginning in Ω tends to M. as t. c c c Ω
30 Example: LaSalle s Theorem d dt x x, x = m kx m cx V V ( x, x ) = mx + kx V x x = cx (, ) M Ω c x E x
31 Lagrangian Systems (, ) L( xx, ) d L xx T = Q dt x x n x R generalized coordinates x = dx / dt generalized velocities ( ) ( ) T T( xx ) = xm ( x) x ( ) = ( ) + L R R L xx = T xx U x n : is the Lagrangian,,, kinetic energy:, total energy: V xx, T xx, U x
32 Lagrange-Poincare Systems (, ) L( px, ) d L px x = p, = Q dt p x T (, ) = (, ), (, ) = (, ) d (, ) =, = + L px T px U x T px pm x p pm x pm x p p dt p x T L px T L px T T pm x p M x M x p U x pm ( x) + p p + = Q x x x T T T T T x = p T T M x p p M x U x M ( x) p + p+ = Q T x x x
33 Example x x x T U =, =, = + x x x x x = cx ( + x ) x x V ( x) = +, V ( x) = cx for c> + x Notice that x the level sets are unbounded for V x = constant V x is not radially unbounded
34 Example, Cont d 4 4 x x t
35 First Integrals Definition : A first integral of the differential equation (, ) x = f xt ( xt) ( xt) is a scalar function ϕ, that is constant along trajectories, i.e., ( xt) ϕ, ϕ, ϕ ( xt, ) = f( x, t) + x t Observation : For simplicity, consider the autonomous case x = f x. Suppose ϕ x is a first integral ϕn ϕ ( x) and ϕ x,, x are arbitrary independent functions on a neighborhhod of the point x, i.e., det x ϕ ( x) n x= x Then we can define coordinate transformation x ( x) ( z) z, via ϕ z = ϕ ( x) z = f ( x) z = z constant x x= ϕ The problem has been reduced to solving n differential equations.
36 Noether s Theorem If the Lagrangian is invariant under a smooth parameter change of coordinates, h : M M, s R, then the Lagrangian system has a first integral Φ ( pq, ) L = dh s ( q) p ds s= s F X, v M M z θω, m x M + m m cosθ v T ( p, q) = [ v ω], cos U q mg θ m cosθ m ω = M + m m cosθ v m ωsinθ v m cosθ m ω + mωsinθ mvsinθ ω F + mg sinθ = X + s hs ( X, θ ) =, Φ ( pq, ) = ( M+ mv ) Momentum in X direction
37 Chetaev s Method Consider the system of equations x = f xt,, f, t = We wish to study the stability of the equilibrium point x =. ( xt) Obviously, if ϕ, is a first integral and it is also a positive definite function, then V xt, is not positive definite? = ϕ ( xt, ) establishes stability. But suppose ϕ ( xt, ) Suppose the system has kfirst integrals ϕ xt,,, ϕ xt, such that ϕ, t =. Chetaev suggested the construction of Lyapunov functions of the form: k k (, ) = αϕ (, ) + βϕ (, ) V xt xt xt i= i i i= i i k i
38 Chetaev Instability Theorem Let D be a neighborhood of the origin. Suppose there is a function V x : D R and a set D D such that ) V x is C on D, ) the origin belongs to the boundary of D, D, 3) V x >, V x > on D, 4) on the boundary Then the origin is unstable. of D inside D, i.e., on D DV, x = D r U x( t) D
39 Example, Rigid Body, Cont d ( ωy ) Consider the rigid body with spin about the y-axis (intermediate inertia), ω =,, Shifted equations: ω = aω ω + ω x z y y ω = bωω y x z ω = c ω + ω ω z y y x Attempts to prove stability fails. So, try to prove instability. V Consider ω, ω, ω = ωω x y z x z { ω ω ω ω ω ω } {( ω ω ω ) ω ω } Let B =,, + + < r and D =,, B >, > r x y z x y z x y z r x z z so that V > on D and V = on D r < ωy ωy + ωy > ( ωx ωy ωz) Br V = aω ω + ω + c ω + ω ω = ω + ω aω + cω z y y y y x y y z x We can take,, in which case V > on D instability B r D T x
40 Stability of Linear Systems - Summary Consider the linear system x = Ax Choose V x T T T V x = x A P + PA x : = x Qx a) if their exists a positive definite pair of matrices PQ, that satisfy = T x Px T A P + PA = Q (Lyapunov equation) the origin is asymptotically stable. b) if P has at least one negative eigenvalue and Q>, the origin is unstable. c) if the origin is asymptotically stable then for any Q >, there is a unique solution, P >, of the Lyapunov equation. sufficient condition necessary condition
41 Second Order Systems Consider the system T T T Mx + Cx + Kx =, M = M >, C = C >, K = K > (, ) T T E x x = x Mx + x Kx d E x x x T Mx x T Kx x T Cx Kx x T Kx dt (, ) [ ] d E ( x, x ) = dt = + = + + x T Cx Some interesting generalizations: T T ) C, ) C C, 3) K K The anti-symmetric terms correspond to circulatory forces (transfer conductances in power systems) they are nonconservative. The anti-symmetric terms correspond to gyroscope forces they are conservative.
42 Example V, δ -ja Assume uniform damping Assume e= Designate Gen as swing bus -jb -jd 4 -jc V 4, δ 4 Eliminate internal bus 4 e-jf V 3, δ 3 3 θ + γθ = P b sin θ b θ θ sin sin 3 θ + γθ = P + b sin θ θ b θ 3 θ = δ δ, θ = δ δ, P = P P, P = P P 3 3
43 Example Cont d This is a Lagrangian system with ( θ, θ ) = θ θ cos( θ ) cos( θ θ ) cos( θ ) U P P b b b T 3 3 ( ω, ω ) = ( ω + ω ), Q = [ γω γω ] To study stability choose total energy as Lyapunov function ( ω, ω ) ( θ, θ ) V = T + U V = γω γω T( ω ω) ( ω ω) U ( θ θ ) Note: T, = and, >, Equilibria corresponding to, a local minimum are stable. G. V. Arononvich and N. A. Kartvelishvili, "Application of Stability Theory to Static and Dynamic Stability Problems of Power Systems," presented at Second All-union Conference on Theoretical and Applied mechanics, Moscow, 965.
44 Example, Cont d Since π θ < π and π θ < π we should consider ( θ θ ) U, as a function on a torus U : T R P =, P =, b =, b =, b =
45 Example, Cont d P =.5, P =, b =, b =, b = 3 3
46 Example Cont d P = π /5, P =, b =, b =, b = 3 3
47 Example Cont d P = π /5, P =, b =, b =.5, b = 3 3
48 Control Lyapunov Function Consider the controlled system n x = f xu,, x D R containing x=, u R Find u x such that all trajectories beginning in D converge to x=. Definition : A control Lyapunov function (CLF) is a function V : D R with V =, V x >, x, such that V x D, u( x) V ( x) = f ( x, u( x) ) < x Theorem : (Artsteins Theorem) A differentiable CLF exists iff there exists a 'regular' feedback control u x. m
49 Example Consider the nonlinear mass-spring-damper system: v = q, m + q v + bv + k q + k q = u 3 v d q ( 3 ) dt v = k q kq bv u + m + q suppose the target state is v=, q=, Define r = v + αq, α >. A CLF candidate is V = r. This is positive definite wrt the target state 3 u bv kq kq V = rr = ( v + αq) v = v + αq κv, κ = > m( + q ) 3 u bv kq kq + α = κ = α + κ m( + q ) d q v closed loop dt v = κv αq 3 q v u bv kq kq m q q v
Topic # /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 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 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 information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
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 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 informationNonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points
Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence
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 informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
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 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 informationUsing Lyapunov Theory I
Lecture 33 Stability Analysis of Nonlinear Systems Using Lyapunov heory I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Outline Motivation Definitions
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 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 informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
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 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 informationStability in the sense of Lyapunov
CHAPTER 5 Stability in the sense of Lyapunov Stability is one of the most important properties characterizing a system s qualitative behavior. There are a number of stability concepts used in the study
More informationStability of Nonlinear Systems An Introduction
Stability of Nonlinear Systems An Introduction Michael Baldea Department of Chemical Engineering The University of Texas at Austin April 3, 2012 The Concept of Stability Consider the generic nonlinear
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 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 informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationModeling & Control of Hybrid Systems Chapter 4 Stability
Modeling & Control of Hybrid Systems Chapter 4 Stability Overview 1. Switched systems 2. Lyapunov theory for smooth and linear systems 3. Stability for any switching signal 4. Stability for given switching
More informationME 680- Spring Representation and Stability Concepts
ME 680- Spring 014 Representation and Stability Concepts 1 3. Representation and stability concepts 3.1 Continuous time systems: Consider systems of the form x F(x), x n (1) where F : U Vis a mapping U,V
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 informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
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 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 informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationNonlinear Control Lecture 5: Stability Analysis II
Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41
More informationControl Systems. Internal Stability - LTI systems. L. Lanari
Control Systems Internal Stability - LTI systems L. Lanari outline LTI systems: definitions conditions South stability criterion equilibrium points Nonlinear systems: equilibrium points examples stable
More informationThe Inverted Pendulum
Lab 1 The Inverted Pendulum Lab Objective: We will set up the LQR optimal control problem for the inverted pendulum and compute the solution numerically. Think back to your childhood days when, for entertainment
More informationQuestion 1: A particle starts at rest and moves along a cycloid whose equation is. 2ay y a
Stephen Martin PHYS 10 Homework #1 Question 1: A particle starts at rest and moves along a cycloid whose equation is [ ( ) a y x = ± a cos 1 + ] ay y a There is a gravitational field of strength g in the
More informationThere is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)
82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory
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 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 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 informationNonlinear Control Lecture 4: Stability Analysis I
Nonlinear Control Lecture 4: Stability Analysis I Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 4 1/70
More informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
More informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
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 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 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 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 informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
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 informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
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 informationNonlinear Dynamic Systems Homework 1
Nonlinear Dynamic Systems Homework 1 1. A particle of mass m is constrained to travel along the path shown in Figure 1, which is described by the following function yx = 5x + 1x 4, 1 where x is defined
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
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 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 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 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 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 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 informationNow I switch to nonlinear systems. In this chapter the main object of study will be
Chapter 4 Stability 4.1 Autonomous systems Now I switch to nonlinear systems. In this chapter the main object of study will be ẋ = f(x), x(t) X R k, f : X R k, (4.1) where f is supposed to be locally Lipschitz
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 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 informationSummary of topics relevant for the final. p. 1
Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a
More informationAppendix A Equations of Motion in the Configuration and State Spaces
Appendix A Equations of Motion in the Configuration and State Spaces A.1 Discrete Linear Systems A.1.1 Configuration Space Consider a system with a single degree of freedom and assume that the equation
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationRobotics. Control Theory. Marc Toussaint U Stuttgart
Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, Linear-Quadratic optimal control, Riccati equations (differential, algebraic, discrete-time), controllability,
More informationMEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction
MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/2/214 Outline Solving State Equations Variation
More informationTranscendental cases in stability problem. Hamiltonian systems
of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationRobust Control 5 Nominal Controller Design Continued
Robust Control 5 Nominal Controller Design Continued Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 4/14/2003 Outline he LQR Problem A Generalization to LQR Min-Max
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 informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationControl of Mobile Robots
Control of Mobile Robots Regulation and trajectory tracking Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Organization and
More information3. Fundamentals of Lyapunov Theory
Applied Nonlinear Control Nguyen an ien -.. Fundamentals of Lyapunov heory he objective of this chapter is to present Lyapunov stability theorem and illustrate its use in the analysis and the design of
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Spring Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs)
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs) is
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 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 informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationTranslational and Rotational Dynamics!
Translational and Rotational Dynamics Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Copyright 217 by Robert Stengel. All rights reserved. For educational use only.
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 informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More informationStabilization of a 3D Rigid Pendulum
25 American Control Conference June 8-, 25. Portland, OR, USA ThC5.6 Stabilization of a 3D Rigid Pendulum Nalin A. Chaturvedi, Fabio Bacconi, Amit K. Sanyal, Dennis Bernstein, N. Harris McClamroch Department
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 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 informationHybrid Systems Course Lyapunov stability
Hybrid Systems Course Lyapunov stability OUTLINE Focus: stability of an equilibrium point continuous systems decribed by ordinary differential equations (brief review) hybrid automata OUTLINE Focus: stability
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
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 informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
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 informationControllability, Observability & Local Decompositions
ontrollability, Observability & Local Decompositions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Lie Bracket Distributions ontrollability ontrollability Distributions
More informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
More informationPredictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado
Chapter 4 Dynamic Behavior Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado in Texas? Talk given by Edward Lorenz, December 972 meeting of the American Association for the
More informationControl of Robotic Manipulators
Control of Robotic Manipulators Set Point Control Technique 1: Joint PD Control Joint torque Joint position error Joint velocity error Why 0? Equivalent to adding a virtual spring and damper to the joints
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationDefinition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.
5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that
More informationNonlinear System Analysis
Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control
More informationPhysics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics
Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics Many interesting physics systems describe systems of particles on which many forces are acting. Some of these forces are immediately
More informationStability Implications of Bendixson s Criterion
Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,
More information