CDS 101: Lecture 5.1 Controllability and State Space Feedback

Similar documents
CDS 101: Lecture 5.1 Reachability and State Space Feedback

CDS 101: Lecture 5.1 Reachability and State Space Feedback

CDS 101: Lecture 5-1 Reachability and State Space Feedback

CDS 101: Lecture 8.2 Tools for PID & Loop Shaping

CS537. Numerical Analysis and Computing

CS321. Numerical Analysis and Computing

Castiel, Supernatural, Season 6, Episode 18

Notes on iteration and Newton s method. Iteration

CDS 101: Lecture 4.1 Linear Systems

Lyapunov Stability Analysis for Feedback Control Design

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.

Introduction to Signals and Systems, Part V: Lecture Summary

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

State Space Representation

Time-Domain Representations of LTI Systems

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

EE Control Systems

Chapter 2 Feedback Control Theory Continued

TAYLOR AND MACLAURIN SERIES

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

Answer: 1(A); 2(C); 3(A); 4(D); 5(B); 6(A); 7(C); 8(C); 9(A); 10(A); 11(A); 12(C); 13(C)

Chapter 4 : Laplace Transform

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

Taylor Polynomials and Approximations - Classwork

Appendix: The Laplace Transform

Math 312 Lecture Notes One Dimensional Maps

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

Chapter 7: The z-transform. Chih-Wei Liu

Introduction to Algorithms

Signals and Systems. Problem Set: From Continuous-Time to Discrete-Time

Continuous Models for Eigenvalue Problems

Chapter 4. Fourier Series

Dynamic Response of Linear Systems

Representing Functions as Power Series. 3 n ...

Algebra II Notes Unit Seven: Powers, Roots, and Radicals

( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.

Introduction to Algorithms

Crash course part 2. Frequency compensation

You may work in pairs or purely individually for this assignment.

Discrete-Time Signals and Systems. Discrete-Time Signals and Systems. Signal Symmetry. Elementary Discrete-Time Signals.

Lecture 7: Polar representation of complex numbers

SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES

ADVANCED DIGITAL SIGNAL PROCESSING

( ) = p and P( i = b) = q.

Discrete-Time System Properties. Discrete-Time System Properties. Terminology: Implication. Terminology: Equivalence. Reference: Section 2.

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

Maximum and Minimum Values

Applications in Linear Algebra and Uses of Technology

Example 2. Find the upper bound for the remainder for the approximation from Example 1.

LINEARIZATION OF NONLINEAR EQUATIONS By Dominick Andrisani. dg x. ( ) ( ) dx

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 4 Solutions [Numerical Methods]

2.004 Dynamics and Control II Spring 2008

State space systems analysis

Olli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5

Exponential and Trigonometric Functions Lesson #1

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

Mon Apr Second derivative test, and maybe another conic diagonalization example. Announcements: Warm-up Exercise:

DIGITAL SIGNAL PROCESSING LECTURE 3

September 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1

Lab(8) controller design using root locus

CALIFORNIA INSTITUTE OF TECHNOLOGY. Control and Dynamical Systems CDS 270. Eugene Lavretsky Spring 2007

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

Lesson 10: Limits and Continuity

Some Variants of Newton's Method with Fifth-Order and Fourth-Order Convergence for Solving Nonlinear Equations

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. XII - Feedback Linearization of Nonlinear Systems - Alberto Isidori and Claudio De Persis

Statistical and Mathematical Methods DS-GA 1002 December 8, Sample Final Problems Solutions

Transfer Function Analysis

Chapter 6 Principles of Data Reduction

Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

8. Applications To Linear Differential Equations

Taylor Series (BC Only)

Lecture 13. Graphical representation of the frequency response. Luca Ferrarini - Basic Automatic Control 1

Math E-21b Spring 2018 Homework #2

Student s Printed Name:

Mon Feb matrix inverses. Announcements: Warm-up Exercise:

Chapter 7 z-transform

Calculus 2 - D. Yuen Final Exam Review (Version 11/22/2017. Please report any possible typos.)

1.225J J (ESD 205) Transportation Flow Systems

mx bx kx F t. dt IR I LI V t, Q LQ RQ V t,

MATH 372: Difference Equations

EE 4343 Lab#4 PID Control Design of Rigid Bodies

FIR Filter Design: Part I

Lecture 11. Solution of Nonlinear Equations - III

NBHM QUESTION 2007 Section 1 : Algebra Q1. Let G be a group of order n. Which of the following conditions imply that G is abelian?

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

Bayesian Methods: Introduction to Multi-parameter Models

ME 539, Fall 2008: Learning-Based Control

µ and π p i.e. Point Estimation x And, more generally, the population proportion is approximately equal to a sample proportion

Solution of EECS 315 Final Examination F09

ECE 308 Discrete-Time Signals and Systems

Continuous Functions

Run-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE

Precalculus MATH Sections 3.1, 3.2, 3.3. Exponential, Logistic and Logarithmic Functions

2C09 Design for seismic and climate changes

a for a 1 1 matrix. a b a b 2 2 matrix: We define det ad bc 3 3 matrix: We define a a a a a a a a a a a a a a a a a a

Abstract Vector Spaces. Abstract Vector Spaces

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

Transcription:

CDS, Lecture 5. CDS : Lecture 5. Cotrollability ad State Space Feedback Richard M. Murray 8 October Goals: Deie cotrollability o a cotrol system Give tests or cotrollability o liear systems ad apply to eamples Describe the desig o state eedback cotrollers or liear systems Readig: Packard, Poola ad Horowitz, Dyamic Systems ad Feedback, Sectio 5 Optioal: Friedlad, Sectios 5.-5.4 (cotrollability oly Optioal: A. D. Lewis, A Mathematical Itroductio to Feedback Cotrol, Chapter (available o web page Lecture Review 4.: rom Liear Last Systems Week u = A + Bu y = C+ ( = y + - 5 + - 5.5-5 - 5 -.5 5-5 t ( = At A( t τ ( + ( τ τ + ( τ = yt Ce Ce Bu d t Properties o liear systems Liearity with respect to iitial coditio ad iputs Stability characterized by eigevalues May applicatios ad tools available Provide local descriptio or oliear systems 8 Oct CDS /8/

CDS, Lecture 5. Cotrol Desig Cocepts System descriptio: sigle iput, sigle output oliear system (MIMO also OK = (, u, ( give y = h(, u Stability: stabilize the system aroud a equilibrium poit Give equilibrium poit e, id cotrol law u=α( such that lim t ( = or all ( t e Cotrollability: steer the system betwee two poits Give,, id a iput u(t such that = u t t = T = (, ( takes ( ( Trackig: track a give output trajectory Give y d (t, id u=α(,t such that ( lim yt ( y( t = or all ( t d y(t y d (t t 8 Oct CDS 3 Cotrollability o Liear Systems = A + Bu y = C+, ( give De A liear system is cotrollable i or ay, ad ay time T > there eists a iput u:[,t] such that the solutio o the dyamics startig rom (= ad applyig iput u(t gives (T=. Remarks I the deiitio, ad do ot have to be equilibrium poits we do t ecessarily stay at ater time T. Cotrollability is deied i terms o states does t deped o output Ca characterize cotrollability by lookig at the geeral solutio o a liear system: T AT A( T τ ( T = e + e Bu( τ dτ τ = I itegral is ull rak, the we ca id a iput to achieve ay desired ial state. 8 Oct CDS 4 /8/

CDS, Lecture 5. Tests or Cotrollability Thm A liear system is cotrollable i ad oly i the cotrollability matri is ull rak. = A + Bu y = C+, ( give B AB A B A B T ( = AT A( T τ + ( τ τ = T e e Bu dτ Remarks Very simple test to apply. I MATLAB, use ctrb(a,b I this test is satisied, we say the pair (A,B is cotrollable Some isight ito the proo ca be see by epadig the matri epoetial = + ( + ( + + ( + B ABT ( τ ABT ( τ A BT ( τ AT ( τ e B I A T τ A T τ A T τ B = + + + + + 8 Oct CDS 5 Eample #: Liearized pedulum o a cart θ m Questio: ca we locally cotrol the positio o the cart by proper choice o iput? Approach: look at the liearizatio aroud the upright positio (good approimatio to the ull dyamics i θ remais small M 8 Oct CDS 6 /8/ 3

CDS, Lecture 5. Eample #, co t: Liearized pedulum o a cart θ m M Cotrollability matri B AB A B A 3 B Full rak as log as costats are such that colums ad 3 are ot multiples o each other cotrollable as log as g(m+m ca steer liearizatio betwee poits by proper choice o iput 8 Oct CDS 7 Costructive Cotrollability Give that system is cotrollable, how do we id iput to traser betwee states? Simple case: chai o itegrators 3 = = = = u Fid curve (t such that ( ( = ( ( = ( T ( T ( ( T Choose iput as ut ( = ( t ( Cotrollable caoical orm I cotrollable, ca show there eists a liear chage o coordiates such that z = z+ u a a a z z z = ( z ut ( = z ( t + ( az ( + + az 8 Oct CDS 8 /8/ 4

CDS, Lecture 5. Cotrol Desig Cocepts System descriptio: sigle iput, sigle output oliear system (MIMO also OK = (, u, ( give y = h(, u Stability: stabilize the system aroud a equilibrium poit Give equilibrium poit e, id cotrol law u=α( such that lim t ( = e or all ( t Cotrollability: steer the system betwee two poits Give,, id a iput u(t such that = (, u( t takes ( t = ( t = Trackig: track a give output trajectory Give y r (t, id u=α(,t such that ( lim yt ( y( t = or all ( t d 8 Oct CDS 9 y(t y d (t t State space cotroller desig or liear systems = A + Bu y = C+, ( give T ( = AT A( T τ + ( τ τ = T e e Bu dτ Goal: id a liear cotrol law u=k such that the closed loop system = A + BK = ( A + BK is stable at e =. Remarks Stability determied by eigevalues use K to make eigevalues o (A+BK stable Ca also lik eigevalues to perormace (eg, iitial coditio respose Questio: whe ca we place the eigevalues ayplace that we wat? Thm The eigevalues o (A+BK ca be set to arbitrary values i ad oly i the pair (A,B is cotrollable. MATLAB: K = place(a, B, eigs 8 Oct CDS /8/ 5

CDS, Lecture 5. Natural dyamics = b a b r = a d b Eample #: Predator prey Cotrolled dyamics: modulate ood supply = b ( + u a b r = a d b 8 6 Q: ca we move rom some iitial populatio o oes ad rabbits to a speciied oe i time T by modulatio i the ood supply? Q: ca we stabilize the populatio aroud the curret equilibrum poit Approach: try to aswer this questio locally, aroud the equilibrium poit 5 5 8 Oct CDS 4 8 6 4 Eample #: Problem setup Equilibrium poit calculatio = br ( + u a b = a d b e = (78.76 94.46 Liearizatio Compute liearizatio aroud equil. poit, e : A= B = u ( e, ue ( e, ue Redeie local variables: z=- e, v=u-u e % Compute the equil poit % predprey.m cotais dyamics = ilie('predprey(,'; eq = solve(, [5,5]; % Compute liearizatio A = [ br - a*eq( - *b*eq(, -a*eq(; a*eq(, -d + a*eq( - *b*eq( ]; B = [br*eq(; ]; d z br a, e b, e a, e z b r, e v dt z = a, e d a, e b, e z + + Cotrollable? YES, i b r,a (check [B AB] ca locally steer to ay poit 8 Oct CDS /8/ 6

CDS, Lecture 5. Eample #: Stabilizatio via eigevalue assigmet d z br a, e b, e a, e z b r, e v dt z = a, e d a, e b +, e z + Cotrol desig: v = Kz = K( e u = u + v Place poles at stable values Choose λ=-, - K = place(a, B, [-; -]; e Modiy dyamics to iclude cotrol = b ( K( a b a d b r e = 8 6 4 8 6 4 5 5 8 Oct CDS 3 How to Assig Eigevalues Eigevalue locatio determies time respose For each eigevalue λ i =σ i + jω i, get cotributio o the orm ( ω ω σt y ( t = e asi( t + bcos( t i Imag Ais - Pole-zero map ω p Amplitude To: Y(..5..5 Step Respose From: U( T π/ω p Repeated eigevalues ca give additioal terms o the orm t k e σ +jω - -5 - -5 5 Real Ais 5 5 5 3 Time (sec. Eigevalue assigemet chages the dyamics o the system Illustrates oe o the mai priciples o eedback Ca be used to make ustable poits stable, icrease the speed o the respose, etc Cautio: eigevalue assigmet aects magitude o iput required The amout o actuator eort required to chage the dyamics rom the atural (ope loop dyamics to the desired (closed loop dyamics ca be large We will lear more sophisticated ways to deal with this tradeo later i the course 8 Oct CDS 4 /8/ 7

CDS, Lecture 5. State Space Feedback Desig Tools Eigevalue assigmet (also called pole placemet Choose eigevalue locatios to correspod to desired dyamics Optimal cotrol (CDS b Choose the eedback that miimizes a cost uctio: ( T T J = Q+ urudt u = K K = lqr(a, B, Q, R Key advatage is the eplicit tradeo betwee state error ad iput magitude Ote very heard to relate the weights (Q, R to the desired system behavior Noliear systems (CDS I liearizatio is cotrollable, ca show that liear eedback provides local stability More geerally, ca look at oliear state eedback, u = α( Optimal cotrol problem ca t be solved i closed orm Alterative: Lyapuov based desig ( itegrator backsteppig, Lypuov redesig, etc 8 Oct CDS 5 Summary: Cotrollability ad State Space Feedback = A + Bu y = C+ B AB A B A B u = u + K( e e Key cocepts Cotrollability: id u s.t. Cotrollability rak test or liear systems State eedback to assig eigevalues 8 6 4 8 6 4 5 5 8 6 4 8 6 4 5 5 8 Oct CDS 6 /8/ 8