Dynamical Systems & Lyapunov Stability

Size: px
Start display at page:

Download "Dynamical Systems & Lyapunov Stability"

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 # /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 information

Lyapunov Stability Theory

Lyapunov 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 information

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015

EN 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 information

1 Lyapunov theory of stability

1 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 information

MCE/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 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 information

Solution of Additional Exercises for Chapter 4

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 information

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Nonlinear 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 information

2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ

2 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 information

Nonlinear 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 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 information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/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 information

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Lecture 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 information

Using Lyapunov Theory I

Using 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 information

Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS

Lyapunov 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 information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 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 information

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Introduction 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 information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 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 information

Georgia Institute of Technology Nonlinear Controls Theory Primer ME 6402

Georgia 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 information

Stabilization and Passivity-Based Control

Stabilization 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 information

Stability in the sense of Lyapunov

Stability 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 information

Stability of Nonlinear Systems An Introduction

Stability 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 information

Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait

Prof. 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 information

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

ENGI 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 information

LECTURE 8: DYNAMICAL SYSTEMS 7

LECTURE 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 information

Modeling & Control of Hybrid Systems Chapter 4 Stability

Modeling & 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 information

ME 680- Spring Representation and Stability Concepts

ME 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 information

STABILITY. Phase portraits and local stability

STABILITY. 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 information

Lecture 9 Nonlinear Control Design

Lecture 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 information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

1. 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 information

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

Review: 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 information

In the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as

In 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 information

Nonlinear Control Lecture 5: Stability Analysis II

Nonlinear 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 information

Control Systems. Internal Stability - LTI systems. L. Lanari

Control 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 information

The Inverted Pendulum

The 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 information

Question 1: A particle starts at rest and moves along a cycloid whose equation is. 2ay y a

Question 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 information

There 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)

There 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 information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability 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 information

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

Nonlinear 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 information

16.30/31, Fall 2010 Recitation # 13

16.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 information

Nonlinear Control Lecture 4: Stability Analysis I

Nonlinear 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 information

Robust Control 2 Controllability, Observability & Transfer Functions

Robust 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 information

EE363 homework 7 solutions

EE363 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 information

Hybrid Systems - Lecture n. 3 Lyapunov stability

Hybrid 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 information

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems

Chapter #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 information

Intermediate Differential Equations. Stability and Bifurcation II. John A. Burns

Intermediate 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 information

3 Stability and Lyapunov Functions

3 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 information

2.10 Saddles, Nodes, Foci and Centers

2.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 information

Robot Control Basics CS 685

Robot 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 information

Nonlinear systems. Lyapunov stability theory. G. Ferrari Trecate

Nonlinear 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 information

UNIVERSITY OF MANITOBA

UNIVERSITY 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 information

MEM 255 Introduction to Control Systems: Modeling & analyzing systems

MEM 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 information

Nonlinear Dynamic Systems Homework 1

Nonlinear 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 information

THEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)

THEODORE 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 information

CDS 101 Precourse Phase Plane Analysis and Stability

CDS 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 information

Nonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points

Nonlinear 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 information

1 The pendulum equation

1 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 information

STABILITY ANALYSIS OF DYNAMIC SYSTEMS

STABILITY 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 information

Handout 2: Invariant Sets and Stability

Handout 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 information

Introduction to Nonlinear Control Lecture # 4 Passivity

Introduction 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 information

Now I switch to nonlinear systems. In this chapter the main object of study will be

Now 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 information

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

154 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 information

LMI Methods in Optimal and Robust Control

LMI 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 information

Summary of topics relevant for the final. p. 1

Summary 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 information

Appendix A Equations of Motion in the Configuration and State Spaces

Appendix 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.

ẋ = 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 information

Robotics. Control Theory. Marc Toussaint U Stuttgart

Robotics. 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 information

MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction

MEM 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 information

Transcendental cases in stability problem. Hamiltonian systems

Transcendental 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 information

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

Applied 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 information

Robust Control 5 Nominal Controller Design Continued

Robust 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 information

MCE693/793: Analysis and Control of Nonlinear Systems

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 information

for 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

for 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 information

Alberto Bressan. Department of Mathematics, Penn State University

Alberto 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 information

Control of Mobile Robots

Control 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 information

3. Fundamentals of Lyapunov Theory

3. 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 information

Mathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Spring Department of Mathematics

Mathematical 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 information

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67

ECEN 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 information

Mathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics

Mathematical 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 information

Chapter III. Stability of Linear Systems

Chapter 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 information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION 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 information

Chap. 3. Controlled Systems, Controllability

Chap. 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 information

Translational and Rotational Dynamics!

Translational 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 information

Invariant Manifolds of Dynamical Systems and an application to Space Exploration

Invariant 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 information

The Liapunov Method for Determining Stability (DRAFT)

The 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 information

Stabilization of a 3D Rigid Pendulum

Stabilization 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 information

A plane autonomous system is a pair of simultaneous first-order differential equations,

A 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 information

Half of Final Exam Name: Practice Problems October 28, 2014

Half 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 information

Hybrid Systems Course Lyapunov stability

Hybrid 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 information

Lecture 41: Highlights

Lecture 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 information

Math 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 # 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 information

Structural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.

Structural 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 information

Modeling and Analysis of Dynamic Systems

Modeling 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 information

Controllability, Observability & Local Decompositions

Controllability, 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 information

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.

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. 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 information

Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado

Predictability: 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 information

Control of Robotic Manipulators

Control 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 information

A conjecture on sustained oscillations for a closed-loop heat equation

A 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 information

Definition 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.

Definition 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 information

Nonlinear System Analysis

Nonlinear 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 information

Physics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics

Physics 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 information

Stability Implications of Bendixson s Criterion

Stability 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