Mathematics for Control Theory

Similar documents
MCE693/793: Analysis and Control of Nonlinear Systems

Lecture 2: Controllability of nonlinear systems

Mathematical Systems Theory: Advanced Course Exercise Session 5. 1 Accessibility of a nonlinear system

Chap. 3. Controlled Systems, Controllability

Controllability, Observability & Local Decompositions

Nonlinear Control Lecture 9: Feedback Linearization

Chap. 1. Some Differential Geometric Tools

A Sufficient Condition for Local Controllability

3.2 Frobenius Theorem

Lecture : Feedback Linearization

Output Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x)

Disturbance Decoupling Problem

MCE693/793: Analysis and Control of Nonlinear Systems

Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.

Exact Linearization Of Stochastic Dynamical Systems By State Space Coordinate Transformation And Feedback I g-linearization

Feedback Linearization

Topic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis

Symplectic and Poisson Manifolds

Introduction to Geometric Control

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

MCE693/793: Analysis and Control of Nonlinear Systems

Stabilization and Passivity-Based Control

Applied Nonlinear Control

Math 1060 Linear Algebra Homework Exercises 1 1. Find the complete solutions (if any!) to each of the following systems of simultaneous equations:

Distributions, Frobenious Theorem & Transformations

LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups

CONTROL OF THE NONHOLONOMIC INTEGRATOR

Nonlinear Control. Nonlinear Control Lecture # 25 State Feedback Stabilization

Nonlinear Control Systems

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1

fy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f))

MCE693/793: Analysis and Control of Nonlinear Systems

Complex line bundles. Chapter Connections of line bundle. Consider a complex line bundle L M. For any integer k N, let

A Do It Yourself Guide to Linear Algebra

10. Smooth Varieties. 82 Andreas Gathmann

APPLICATION OF INPUT STATE OF THE SYSTEM TRANSFORMATION FOR LINEARIZATION OF SELECTED ELECTRICAL CIRCUITS

Chapter 1 Numerical approximation of data : interpolation, least squares method

Observability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)

Math Lecture 4 Limit Laws

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.

LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS. Suba Thomas Harry G. Kwatny Bor-Chin Chang

Chapter 6. Differentially Flat Systems

Kevin James. MTHSC 3110 Section 4.3 Linear Independence in Vector Sp

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

Representations and Linear Actions

Applied Linear Algebra

Linear Algebra Massoud Malek

Zangwill s Global Convergence Theorem

Invariant Manifolds of Dynamical Systems and an application to Space Exploration

Math 480 The Vector Space of Differentiable Functions

The local structure of affine systems

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

Quotients of fully nonlinear control systems

Lyapunov Stability Theory

Vector Space Concepts

Nonholonomic Behavior in Robotic Systems

Gait Controllability for Legged Robots

Georgia Institute of Technology Nonlinear Controls Theory Primer ME 6402

REMARKS ON THE TIME-OPTIMAL CONTROL OF A CLASS OF HAMILTONIAN SYSTEMS. Eduardo D. Sontag. SYCON - Rutgers Center for Systems and Control

Vectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1. x 2. x =

Review of Some Concepts from Linear Algebra: Part 2

Robotics, Geometry and Control - A Preview

Linear Algebra Lecture Notes-I

Extremal Trajectories for Bounded Velocity Differential Drive Robots

Tangent bundles, vector fields

Vector Spaces. Addition : R n R n R n Scalar multiplication : R R n R n.

ENGINEERING MATH 1 Fall 2009 VECTOR SPACES

Vector Spaces and Linear Transformations

MATH 1902: Mathematics for the Physical Sciences I

AFFINE AND PROJECTIVE GEOMETRY, E. Rosado & S.L. Rueda 4. BASES AND DIMENSION

On the convergence of normal forms for analytic control systems

Chapter 3a Topics in differentiation. Problems in differentiation. Problems in differentiation. LC Abueg: mathematical economics

Geodesic Equivalence in sub-riemannian Geometry

INTRODUCTION TO LIE ALGEBRAS. LECTURE 2.

Hamiltonian flows, cotangent lifts, and momentum maps

El péndulo invertido: un banco de pruebas para el control no lineal. XXV Jornadas de Automática

2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?

Cambridge University Press The Mathematics of Signal Processing Steven B. Damelin and Willard Miller Excerpt More information

Chapter 3. Vector spaces

An Introduction to Differential Flatness. Will Burgoyne May 7, 2007 ME598 Geometric Mechanics

CHAPTER 7 DIV, GRAD, AND CURL

2.3. VECTOR SPACES 25

Geometric Control Theory

1 Antiderivatives graphically and numerically

Adaptive backstepping for trajectory tracking of nonlinearly parameterized class of nonlinear systems

Linear Algebra and Robot Modeling

CHAPTER 1. Theory of Second Order. Linear ODE s

Linear Algebra. Linear Algebra. Chih-Wei Yi. Dept. of Computer Science National Chiao Tung University. November 12, 2008

LINEAR ALGEBRA W W L CHEN

MA 262, Fall 2017, Final Version 01(Green)

Global Formulation and Control of a Class of Nonholonomic Systems

(1) is an invertible sheaf on X, which is generated by the global sections

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

Pattern generation, topology, and non-holonomic systems

INPUT-STATE LINEARIZATION OF A ROTARY INVERTED PENDULUM

BIBO STABILITY AND ASYMPTOTIC STABILITY

Geometric Homogeneity and Configuration Controllability of Nonlinear Systems

Lecture 7. Lie brackets and integrability

Vector Spaces. (1) Every vector space V has a zero vector 0 V

i = f iα : φ i (U i ) ψ α (V α ) which satisfy 1 ) Df iα = Df jβ D(φ j φ 1 i ). (39)

Transcription:

Mathematics for Control Theory Geometric Concepts in Control Involutivity and Frobenius Theorem Exact Linearization Hanz Richter Mechanical Engineering Department Cleveland State University

Reading materials Reference: Harry Kwatny and Gilmer Blankenship [2000], Nonlinear Control and Analytical Mechanics: A Computational Approach, Birkhäuser, ISBN 978-0-8176-4147-4, Sections 3.3.4 and 3.3.5. www.pages.drexel.edu/~hgk22/courses/mem636_638/ncam.pdf Jean-Jacques Slotine and Weiping Li [1991], Applied Nonlinear Control, Prentice-Hall, ISBN: 0-13-040890-5, Chapter 6. 2 / 11

More insight from car parking example Given a point (x p,y p,φ,θ), and two instantaneous values of control u 1 and u 2, the four differential equations describing cart dynamics provide a tangent vector ẋ (an instantaneous direction of motion in 4 coordinates). We see that our choices of instantaneous values of u 1 and u 2 can be used to point the cart in various directions. Pointwise, the tangent vector is a simple linear combination of d(x(t)) and s(x(t)), with u 1 and u 2 as scalar weights. But our problem is not merely to point the cart. We want to reach new points by manipulating u 1 and u 2 continuously. A simple subclass of control actions is given by piecewise constant functions of time, with the restriction that u 1 u 2 = 0 (follow the direction of only one vector field at a time). 3 / 11

Insight from car parking example... If a control action (u 1,u 2 ) with u 1 u 2 = 0 which is constant in a time interval is used, the flow of the differential equation will transfer the initial point to a new one. At this new point, we can form a new linear combination to point the tangent vector in a desired direction. Clearly, an understanding of achievable directions and achievable trajectories is necessary. We used a sequence of forward and backward motions along the flows generated by two vector fields: K(t,x) = Φ t g Φ t f Φ t g Φ t f x The operator K(t,x) yielding the final point starting from x is called the commutator of f and g. 4 / 11

Commutator and Lie Bracket We saw that K(t,x)x = x iff [f,g] = 0. There s a more precise relationship (See Kwatny and Blankenship, Proposition 3.67): K(t,x) defines a continuous path in M and: d dt K(0+,x) = [f,g] x Examine the proof and pay close attention to Fig. 3.7. 5 / 11

Control by Exact Linearization Input to state linearization, also called exact linearization, if possible, provides a simpler way to specify control actions. For instance, a 3rd-order feedback linearizable system ẋ = f(x)+g(x)u can be written as: ż 1 = z 2 ż 2 = z 3 ż f = v where u = u(z,v) is a suitable input transformation and v is a new control. Then control design becomes trivial using v as the control. The state and control transformations are use to express the control law in original coordinates. 6 / 11

Complete Integrability Refer to Slotine and Li, Sect. 6.2 for the motivation of the definition below. Suppose the set {f i (x)},i = 1,2...m is linearly independent in some subset U R n. The set is said to be completely integrable in U if we can find n m scalar functions h : R n R satisfying: h i f j = 0 for 1 i n m, 1 j m, with the set of { h i (x)} linearly independent in U. Under this definition, we have m(n m) partial differential equations. If there s a solution, the set of vector fields is completely integrable. Exercise: Determine if the following set is completely integrable (can you find an subset of R n where the vector fields are linearly independent?) f(x) = [sin(x 1 ) 0 cos(x 2 )] T, g(x) = [cos(x 3 ) x 2 sin(x 2 )] T 7 / 11

Involutivity and Frobenius Theorem Suppose that for any pair of vector fields f i,f j from a linearly independent set {f 1,f 2,...f m }, the Lie bracket [f i,f j ] can always be expressed as a linear combination of the vector fields in the set (with possibly state-dependent coefficients). Then the set is involutive. The test for involutivity is just a pointwise rank test: for any i,j, [f i (x),f j (x)] must already be in the span of {f 1 (x),f 2 (x),...f m (x)} for all x in the subset where involutivity is to be verified. Example: Use a numerical or analytical procedure to test involutivity of the set {f,g} of the previous example in the same region where these fields are linearly independent. Frobenius Theorem: The linearly independent set {f 1,f 2,...f m } is completely integrable if anf only if it is involutive. 8 / 11

Exact Linearization of Single-Input Affine Systems Consider the system ẋ = f(x)+g(x)u where f and g are smooth vector fields. The system is input-to-state, or exactly linearizable if there is a diffeomorphic state transformation Φ : Ω R n with z = Φ(x) and an input transformation u = α(x)+β(x)v in some region Ω R n such that the new state z and the new input v are related by the linear, time-invariant state equation below: ż 1 = z 2 ż 2 = z 3. ż n = v Note that there s no output defined, and if this ambitious transformation is possible, there are no internal dynamics left. Among mechanical systems, robot dynamics are an important case of input-state linearizable systems with many inputs. 9 / 11

Necessary and Sufficient Condition for Exact Linearizability (Theorem 6.2 in Slotine and Li): The single-input, n-state nonlinear system ẋ = f(x)+g(x)u is exactly linearizable iff the following conditions hold in some region Ω R n : 1. The vector fields {g,ad f g,...ad f n 1g} are linearly independent in Ω 2. The set {g,ad f g,...ad f n 2g} is involutive in Ω. What happens with these conditions for linear systems (f(x) = Ax and g(x) = B)? The proof of this Theorem includes a method to construct the diffeomorphic state transformation and the input transformation. Suppose the system is exactly linearizable and let z(x) = [z 1 (x) z 2 (x)...z n (x)] T. Find z 1 first, to satisfy z 1 ad f ig = 0, i = 1,2..n 2 z 1 ad f n 1g 0 10 / 11

Recipe for Exact Linearizability... The state transformation is z(x) = [z 1 L f z 1...L f n 1z 1 ] T while the input transformation uses α(x) = Ln f z 1 L g L f n 1z 1 1 β(x) = L g L f n 1z 1 Examples: We solve Problems 6.3 and 6.7 from Slotine and Li. 11 / 11