Nonlinear Control Theory. Lecture 9
|
|
- Paul Eaton
- 6 years ago
- Views:
Transcription
1 Nonlinear Control Theory Lecture 9 Periodic Perturbations Averaging Singular Perturbations Khalil Chapter (9, 10) , 11
2 Today: Two Time-scales Averaging ẋ = ǫf(t,x,ǫ) The state x moves slowly compared to f. Singular perturbations ẋ = f(t,x,z,ǫ) ǫż = (t,x,z,ǫ) The state x moves slowly compared to z.
3 Example: Vibrating Pendulum I m θ l asinωt x=lsin(θ), y=lcos(θ) asin(ωt)
4 Newton s law in tangential direction m(l θ aω 2 sinωtsinθ) = m sinθ k(l θ+aωcosωtsinθ) (incl. viscous friction in joint) Let ǫ=a/l,τ= ωt, α= ω 0 l/ωa, and β=k/mω 0 x 1 = θ x 2 = ǫ 1 (dθ/dτ)+cosτsinθ f 1 (τ,x) = x 2 cosτsinx 1 f 2 (τ,x) = α βx 2 α 2 sinx 1 the state equation is given by +x 2 cosτcosx 1 cos 2 τsinx 1 cosx 1 dx dτ = ǫf(τ,x)
5 Averaging Assumptions Consider the system ẋ = ǫf(t,x,ǫ), x(0)=x 0 wherefand its derivatives up to second order are continuous and bounded. Letx av be defined by the equations ẋ av = ǫf av (x av ), x av (0)=x 0 1 f av (x) = lim T T T 0 f(τ,x,0)dτ
6 Example: Vibrating Pendulum II The averaged system ẋ = ǫf av (x) [ ] x = ǫ 2 α βx 2 α 2 sinx sin2x 1 has f av x (π,0) = [ ] 0 1 α α β which is Hurwitz for0< α<1/ 2, β>0. Can this be used for rigorous conclusions?
7 Periodic Averaging Theorem Letfbe periodic intwith periodt. Let x = 0 be an exponentially stable equilibrium of ẋ av = ǫf(x av ). If x 0 is sufficiently small, then x(t,ǫ) = x av (t,ǫ)+o(ǫ) for allt [0, ] Furthermore, for sufficiently small ǫ>0, the equation ẋ = ǫ f(t, x, ǫ) has a unique exponentially stable periodic solution of period T in an O(ǫ) neighborhood of x = 0.
8 General Averaging Theorem Under certain conditions on the convergence of 1 T f av (x) = lim f(τ,x,0)dτ T T 0 there exists a C > 0 such that for sufficiently small ǫ > 0 for allt [0,1/ǫ]. x(t,ǫ) x av (t,ǫ) < Cǫ
9 Example: Vibrating Pendulum III The Jacobian of the averaged system is Hurwitz for 0< α<1/ 2, β>0. For a/l sufficiently small and ω > 2ω 0 l/a the unstable pendulum equilibrium(θ, θ)=(π,0) is therefore stabilized by the vibrations.
10 Periodic Perturbation Theorem Consider ẋ = f(x)+ǫ (t,x,ǫ) wheref,, f/ x and / x are continuous and bounded. Let be periodic intwith periodt. Let x = 0 be an exponentially stable equilibrium point for ǫ = 0. Then, for sufficiently small ǫ>0, there is a unique periodic solution which is exponentially stable. x(t, ǫ) = O(ǫ)
11 Proof ideas of Periodic Perturbation Theorem Let φ(t,x 0,ǫ) be the solution of ẋ = f(x)+ǫ (t,x,ǫ), x(0)=x 0 Exponential stability of x = 0 for ǫ = 0, plus bounds on the magnitude of, shows existence of a bounded solution x for small ǫ>0. The implicit function theorem shows solvability of x = φ(t,0,x,ǫ) for small ǫ. This gives periodicity of x. Putz=x x. Exponential stability ofx=0for ǫ=0gives exponential stability ofz=0for small ǫ>0.
12 Proof idea of Averaging Theorem For small ǫ>0defineuandy by Then u(t,x) = t 0 x = y+ ǫu(t,y) [f(τ,x,0) f av (x)]dτ ẋ=ẏ+ ǫ u(t,y) + ǫ u(t,y) ẏ t y [ I+ ǫ u ] ẏ = ǫf(t,y+ ǫu,ǫ) ǫ u y t (t,y) = ǫf av (y)+ǫ 2 p(t,y,ǫ) Withs=ǫt, dy ( s ) = f av (y)+ǫq ds ǫ,y,ǫ which has a unique and exponentially stable periodic solution for small ǫ. This gives the desired result.
13 Application: Second Order Oscillators For the second order system introduce y = rsin φ ÿ+ω 2 y = ǫ (y,ẏ) (1) ẏ/ω = rcosφ (rsin φ,ωrcosφ)cos φ f(φ,r,ǫ) = ω 2 (ǫ/r) (rsin φ,ωrcosφ)sin φ f av (r) = 1 2π = Then (1) is equivalent to 2π 1 2πω 2 dr dφ 0 2π 0 f(φ,r,0)dφ (rsin φ,ωrcosφ)cos φdφ = ǫf(φ,r,ǫ) and the periodic averaging theorem may be applied.
14 Illustration: Van der Pol Oscillator I C i C L i L i + V resistive element linear osc-part For an ordinary resistance we will get a damped oscillation. For a negative resistance/admittance chosen as we get i=h(v) = ( V+ 1 3 V3 ) }{{} ivesvanderpol eq. i C +i L +i=0, i=h(v) CL d2 V dt 2 +V+Lh (V) dv dt =0
15 Example: Van der Pol Oscillator I The vacuum tube circuit equation (a k a thevanderpolequation) gives f av (r) = 1 2π The averaged system ÿ+y = ǫẏ(1 y 2 ) 2π 0 = 1 2 r 1 8 r3 dr dφ rcos φ(1 r 2 sin 2 φ)cos φdφ ( ) 1 = ǫ 2 r 1 8 r3 has equilibriar=0,r=2with df av dr = 1 r=2 so small ǫ give a stable limit cycle, which is close to circular with radiusr=2.
16 Singular Perturbations Consider equations of the form ẋ = f(t,x,z,ǫ), x(0)=x 0 ǫż = (t,x,z,ǫ) z(0)=z 0 For small ǫ>0, the first equation describes theslowdynamics, while the second equation defines thefastdynamics. The main idea will be to approximatexwith the solution of the reduced problem x = f(t,x,h(t, x),0) x(0)=x 0 whereh(t, x) is defined by the equation 0 = (t,x,h(t,x),0)
17 Example: DC Motor I i R L u EMK=kω ω J J dω dt L di dt = ki = kω Ri+u Withx= ω,z=iand ǫ=lk 2 /JR 2 we get ẋ = z ǫż = x z+u
18 Linear Singular Perturbation Theorem Let the matrixa 22 have nonzero eigenvalues γ 1,...,γ m and let λ 1,...,λ n be the eigenvalues ofa 0 =A 11 A 12 A 1 22 A 21. Then, δ>0 ǫ 0 >0such that the eigenvalues α 1,...,α n+m of the matrix [ ] A11 A 12 A 21 /ǫ A 22 /ǫ satisfy the bounds for0<ǫ<ǫ 0. λ i α i < δ, i=1,...,n γ i n ǫα i < δ, i=n+1,...,n+m
19 Proof A 22 is invertible, so it follows from the implicit function theorem that for sufficiently small ǫ the Riccati equation ǫa 11 P ǫ +A 12 ǫp ǫ A 21 P ǫ P ǫ A 22 = 0 has a unique solutionp ǫ =A 12 A O(ǫ). The desired result now follows from the similarity transformation [ ][ ] [ ] I ǫpǫ A11 A 12 I ǫpǫ 0 I A 21 /ǫ A 22 /ǫ 0 I = [ ][ ] I ǫpǫ A11 A 12 + ǫa 11 P ǫ 0 I A 21 /ǫ A 22 /ǫ+a 21 P ǫ [ A0 +O(ǫ) 0 = A 22 /ǫ+o(1) ]
20 Example: DC Motor II In the example we have ẋ = z ǫż = x z+u [ ] A11 A 12 A 21 A 22 = [ ] A 11 A 12 A 1 22 A 21 = 1 so stability of the DC motor model for small is verified. ǫ = Lk2 JR 2 See Khalil for example where reduced system is stable but fast dynamics unstable.
21 The Boundary-Layer System For fixed(t, x) the boundary layer system dŷ dτ = (t,x,ŷ+h(t,x),0), ŷ(0)=z 0 h(0,x 0 ) describes the fast dynamics, disregarding variations in the slow variables t, x.
22 Tikhonov s Theorem Consider a singular perturbation problem with f,,h, / x C 1. Assume that the reduced problem has a unique bounded solution x on[0, T] and that the equilibrium ŷ=0of the boundary layer problem is exponentially stable uniformly in(t,x). Then uniformly fort [0,T]. x(t,ǫ) = x(t)+o(ǫ) z(t,ǫ) = h(t, x(t))+ŷ(t/ǫ)+o(ǫ)
23 Example: High Gain Feedback u + + ψ( ) k 1 /s u p ẋ p =... y=.. y k 2 Closed loop system Reduced model ẋ p = Ax p +Bu p 1 p = ψ(u u p k 2 Cx p ) k 1 u ẋ p = (A Bk 2 C)x p +Bu
24 Proof ideas of Tikhonov s Theorem Replacefand withfandgthat are identical for x <r, but nicer for largex. For small ǫ,g(t,x,y,ǫ) is close tog(t,x,y,0). y-bound for G(,,, 0)-equation y-bound forg-equation x,y-bound forf,g-equations For small ǫ > 0, the x, y-solutions of the F, G-equations will satisfy x <r. Hence, they also solve thef, -equations
25 The Slow Manifold For small ǫ>0, the system ẋ = f(x,z) ǫż = (x,z) has the invariant manifold z = H(x,ǫ) It can often be computed approximately by Taylor expansion H(x,ǫ) = H 0 (x)+ǫh 1 (x)+ǫ 2 H 2 (x)+ whereh 0 satisfies 0 = (x,h 0 )
26 z z The Fast Manifold ǫ=0.001 x = x + z epsilon = z = 1/epsilon atan(1 z x) x x = x + z z = 1/epsilon atan(1 z x) ǫ=0.1 epsilon = x
27 Example: Van der Pol Oscillator III Consider With we have the system with slow manifold d 2 v ds 2 µ(1 v2 ) dv ds +v = 0 x = 1 dv µ ds +v 1 3 v3 z = v t = s/µ ǫ = 1/µ 2 ẋ = z ǫż = x+z 1 3 z3 x = z 1 3 z3
28 z z Illustration: Van der Pol III Phase plot for van der Pol example ǫ=0.001 x = z z = 1/epsilon ( x + z 1/3 z 3 epsilon = ) x x = z z = 1/epsilon ( x + z 1/3 z 3 ) ǫ=0.1 epsilon = x
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 information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More 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 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 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 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 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 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 informationDisturbance Decoupling Problem
DISC Systems and Control Theory of Nonlinear Systems, 21 1 Disturbance Decoupling Problem Lecture 4 Nonlinear Dynamical Control Systems, Chapter 7 The disturbance decoupling problem is a typical example
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
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 informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
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 Lecture 1: Introduction
Nonlinear Control Lecture 1: Introduction Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 1 1/15 Motivation
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 informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
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 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 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 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 informationLinear Systems Theory
ME 3253 Linear Systems Theory Review Class Overview and Introduction 1. How to build dynamic system model for physical system? 2. How to analyze the dynamic system? -- Time domain -- Frequency domain (Laplace
More informationThe Higgins-Selkov oscillator
The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being
More 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 informationnonlinear oscillators. method of averaging
Physics 4 Spring 7 nonlinear oscillators. method of averaging lecture notes, spring semester 7 http://www.phys.uconn.edu/ rozman/courses/p4_7s/ Last modified: April, 7 Oscillator with nonlinear friction
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 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 informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle
High-Gain Observers in Nonlinear Feedback Control Lecture # 2 Separation Principle High-Gain ObserversinNonlinear Feedback ControlLecture # 2Separation Principle p. 1/4 The Class of Systems ẋ = Ax + Bφ(x,
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 informationSolutions. Problems of Chapter 1
Solutions Problems of Chapter 1 1.1 A real square matrix A IR n n is invertible if and only if its determinant is non zero. The determinant of A is a polynomial in the entries a ij of A, whose set of zeros
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 information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More 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 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 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 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 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 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 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 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 informationNonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2009 Farzaneh Abdollahi Nonlinear Control Lecture 2 1/68
More informationLecture Notes. W. Tang. School of Mathematical & Statistical Sciences Arizona State University
Applied Math / Dynamical Systems Workshop Lecture Notes W. Tang School of Mathematical & Statistical Sciences Arizona State University May. 23 rd -25 th, 2011 Tang (ASU) AMDS Workshop Notes May. 23 rd
More informationSystem Control Engineering 0
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/
More 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 informationNonlinear Control Lecture 9: Feedback Linearization
Nonlinear Control Lecture 9: Feedback Linearization Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 9 1/75
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 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 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 informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
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 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 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 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 informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More informationPhys 7221 Homework # 8
Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with
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 informationChapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx
Chapter 1 Lecture Notes Chapter 1 Oscillatory Motion Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx When the mass is released, the spring will pull
More informationE209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions
E9A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions Michael Vitus Gabe Hoffmann Stanford University Winter 7 Problem 1 The governing equations are: ẋ 1 = x 1 + x 1 x ẋ = x + x 3 Using
More 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 informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationAnalysis of Dynamical Systems
2 YFX1520 Nonlinear phase portrait Mathematical Modelling and Nonlinear Dynamics Coursework Assignments 1 ϕ (t) 0-1 -2-6 -4-2 0 2 4 6 ϕ(t) Dmitri Kartofelev, PhD 2018 Variant 1 Part 1: Liénard type equation
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation
High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y
More information10 Transfer Matrix Models
MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important
More informationA Study of the Van der Pol Equation
A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
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 informationNonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
More informationStability and Control of dc Micro-grids
Stability and Control of dc Micro-grids Alexis Kwasinski Thank you to Mr. Chimaobi N. Onwuchekwa (who has been working on boundary controllers) May, 011 1 Alexis Kwasinski, 011 Overview Introduction Constant-power-load
More informationPHY 3221 Fall Homework Problems. Instructor: Yoonseok Lee. Submit only HW s. EX s are additional problems that I encourage you to work on.
PHY 3221 Fall 2012 Homework Problems Instructor: Yoonseok Lee Submit only HW s. EX s are additional problems that I encourage you to work on. Week 1: August 22-24, Due August 27 (nothing to submit) EX:
More informationLecture 7. Please note. Additional tutorial. Please note that there is no lecture on Tuesday, 15 November 2011.
Lecture 7 3 Ordinary differential equations (ODEs) (continued) 6 Linear equations of second order 7 Systems of differential equations Please note Please note that there is no lecture on Tuesday, 15 November
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 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 informationLab 6a: Pole Placement for the Inverted Pendulum
Lab 6a: Pole Placement for the Inverted Pendulum Idiot. Above her head was the only stable place in the cosmos, the only refuge from the damnation of the Panta Rei, and she guessed it was the Pendulum
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 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 informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationModule Two: Differential Calculus(continued) synopsis of results and problems (student copy)
Module Two: Differential Calculus(continued) synopsis of results and problems (student copy) Srikanth K S 1 Syllabus Taylor s and Maclaurin s theorems for function of one variable(statement only)- problems.
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: April 06, 2018, at 14 00 18 00 Johanneberg Kristian
More informationPhysics 235 Chapter 4. Chapter 4 Non-Linear Oscillations and Chaos
Chapter 4 Non-Linear Oscillations and Chaos Non-Linear Differential Equations Up to now we have considered differential equations with terms that are proportional to the acceleration, the velocity, and
More informationControl Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli
Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationNumerics and Control of PDEs Lecture 1. IFCAM IISc Bangalore
1/1 Numerics and Control of PDEs Lecture 1 IFCAM IISc Bangalore July 22 August 2, 2013 Introduction to feedback stabilization Stabilizability of F.D.S. Mythily R., Praveen C., Jean-Pierre R. 2/1 Q1. Controllability.
More informationSummer 2017 MATH Solution to Exercise 5
Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a
More informationNov : Lecture 22: Differential Operators, Harmonic Oscillators
14 Nov. 3 005: Lecture : Differential Operators, Harmonic Oscillators Reading: Kreyszig Sections:.4 (pp:81 83),.5 (pp:83 89),.8 (pp:101 03) Differential Operators The idea of a function as something that
More informationNPTEL Online Course: Control Engineering
NPTEL Online Course: Control Engineering Ramkrishna Pasumarthy Assignment-11 : s 1. Consider a system described by state space model [ ] [ 0 1 1 x + u 5 1 2] y = [ 1 2 ] x What is the transfer function
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 informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter
More informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
More informationChapter 4: Partial differentiation
Chapter 4: Partial differentiation It is generally the case that derivatives are introduced in terms of functions of a single variable. For example, y = f (x), then dy dx = df dx = f. However, most of
More informationFirst Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields
HT017 First Year Physics: Prelims CP1 Classical Mechanics: Prof Neville Harnew Problem Set III : Projectiles, rocket motion and motion in E & B fields Questions 1-10 are standard examples Questions 11-1
More informationSystem Control Engineering
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Feedback Systems: An Introduction for Scientists and Engineers Author: Karl J. Astrom and Richard M. Murray Text:
More informationSingular perturbation methods for slow fast dynamics
Nonlinear Dyn (2007) 50:747 753 DOI 10.1007/s11071-007-9236-z ORIGINAL ARTICLE Singular perturbation methods for slow fast dynamics Ferdinand Verhulst Received: 2 February 2006 / Accepted: 4 April 2006
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More information