Intro. Computer Control Systems: F9

Size: px
Start display at page:

Download "Intro. Computer Control Systems: F9"

Transcription

1 Intro. Computer Control Systems: F9 State-feedback control and observers Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 21 dave.zachariah@it.uu.se

2 F8: Quiz! 2 / 21 dave.zachariah@it.uu.se

3 F8: Quiz! 1) For an observable system a the effect of all x(t) can be observed in y(t) b we have det O = 0 c we have stability 2 / 21 dave.zachariah@it.uu.se

4 F8: Quiz! 1) For an observable system a the effect of all x(t) can be observed in y(t) b we have det O = 0 c we have stability 2) If a state-space form of G(s) is a minimal realization, a A:s eigenvalues < G(s):s poles b A:s eigenvalues = G(s):s poles c there exist more compact state-space forms 2 / 21 dave.zachariah@it.uu.se

5 F8: Quiz! 1) For an observable system a the effect of all x(t) can be observed in y(t) b we have det O = 0 c we have stability 2) If a state-space form of G(s) is a minimal realization, a A:s eigenvalues < G(s):s poles b A:s eigenvalues = G(s):s poles c there exist more compact state-space forms 3) For a controllable system with state-feedback control a the closed-loop system is stable b the poles of the closed-loop system can be designed arbitrarily c no information about the system is required 2 / 21 dave.zachariah@it.uu.se

6 State-feedback control 3 / 21 dave.zachariah@it.uu.se

7 State-feedback control State-space form of linear time-invariant system ẋ = Ax + Bu y = Cx G(s) = C(sI A) 1 B u x y (si A) 1 B C 4 / 21 dave.zachariah@it.uu.se

8 State-feedback control Controller using state feedback u = Lx + l 0 r gives closed-loop system ẋ = (A BL)x + Bl 0 r y = Cx where r is the reference signal. 4 / 21 dave.zachariah@it.uu.se

9 State-feedback control Controler using state feedback gives closed-loop system u = Lx + l 0 r ẋ = (A BL)x + Bl 0 r y = Cx where r is the reference signal. rl 0 u x y (si A) 1 B C + L G c (s) = C(sI A + BL) 1 Bl 0 4 / 21 dave.zachariah@it.uu.se

10 Pole placement Rules of thumb for designing L Eigenvalues/poles given by det(si A + BL) = 0, which we can design Im{s} Re{s} 5 / 21 dave.zachariah@it.uu.se

11 Pole placement Rules of thumb for designing L Eigenvalues/poles given by det(si A + BL) = 0, which we can design Im{s} Re{s} Distance to the origin: Quick system but also sensitive to disturbances 5 / 21 dave.zachariah@it.uu.se

12 Estimating the states via simulation 6 / 21 dave.zachariah@it.uu.se

13 Estimating the states Via simulation Controller u = Lx + l 0 r requires states x which are often unknown. 7 / 21 dave.zachariah@it.uu.se

14 Estimating the states Via simulation Controller u = Lx + l 0 r requires states x which are often unknown. In practice, feedback using u = Lˆx + l 0 r where ˆx is an estimate of x. 7 / 21 dave.zachariah@it.uu.se

15 Estimating the states Via simulation Controller u = Lx + l 0 r requires states x which are often unknown. In practice, feedback using u = Lˆx + l 0 r where ˆx is an estimate of x. Naive idea: Estimate x by simulating the states where ˆx 0 is an initial guess. ˆx = Aˆx + Bu, ˆx(0) = ˆx 0 7 / 21 dave.zachariah@it.uu.se

16 Build intuition from simple systems State estimation via simulation Ex.: Damper y State-space form: [ ] [ ] ẋ(t) = x(t) + u(t), x(0) = x k/m 0 1/m 0 y(t) = [ 1 0 ] x(t) u 8 / 21 dave.zachariah@it.uu.se

17 Build intuition from simple systems State estimation via simulation Example using impulse u input u(t), output y(t) u(t) y(t) t [s] System with unknown initial state x 0 8 / 21 dave.zachariah@it.uu.se

18 Build intuition from simple systems State estimation via simulation Naive estimate using perfect initial guess: ˆx = Aˆx + Bu, ˆx 0 = x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = / 21 dave.zachariah@it.uu.se

19 Build intuition from simple systems State estimation via simulation Naive estimate using perfect initial guess: ˆx = Aˆx + Bu, ˆx 0 = x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = 20 8 / 21 dave.zachariah@it.uu.se

20 Build intuition from simple systems State estimation via simulation Naive estimate using perfect initial guess: ˆx = Aˆx + Bu, ˆx 0 = x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = / 21 dave.zachariah@it.uu.se

21 Build intuition from simple systems State estimation via simulation Naive estimate using wrong initial guess: ˆx = Aˆx + Bu, ˆx 0 x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = / 21 dave.zachariah@it.uu.se

22 Build intuition from simple systems State estimation via simulation Naive estimate using wrong initial guess: ˆx = Aˆx + Bu, ˆx 0 x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = 20 9 / 21 dave.zachariah@it.uu.se

23 Build intuition from simple systems State estimation via simulation Naive estimate using wrong initial guess: ˆx = Aˆx + Bu, ˆx 0 x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = / 21 dave.zachariah@it.uu.se

24 Build intuition from simple systems State estimation via simulation x och ˆx correspond to different outputs: y = Cx versus ŷ = C ˆx output y(t) y(t) C ˆx(t) t [s] 9 / 21 dave.zachariah@it.uu.se

25 Estimating the states via observer 10 / 21 dave.zachariah@it.uu.se

26 Estimating the states Correcting the state estimates Idea: Feedback the prediction error y C ˆx to correct ˆx Observer: an estimator with a correction term ˆx = Aˆx + Bu + K ( y C ˆx ), ˆx(0) = ˆx 0 }{{} correction 11 / 21 dave.zachariah@it.uu.se

27 Estimating the states Correcting the state estimates Idea: Feedback the prediction error y C ˆx to correct ˆx Observer: an estimator with a correction term Using matrix ˆx = Aˆx + Bu + K ( y C ˆx ), ˆx(0) = ˆx 0 }{{} correction we can design the estimator. k 1 k 2 K =. k n 11 / 21 dave.zachariah@it.uu.se

28 Estimating the states Correcting the state estimates Idea: Feedback the prediction error y C ˆx to correct ˆx Observer: an estimator with a correction term ˆx = Aˆx + Bu + K ( y C ˆx ), ˆx(0) = ˆx 0 }{{} correction u x y (si A) 1 B C ˆx Obs. 11 / 21 dave.zachariah@it.uu.se

29 Build intuition using simple systems State estimation using observer Estimation using observer: ˆx = Aˆx + Bu + K ( y C ˆx ), ˆx 0 x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = / 21 dave.zachariah@it.uu.se

30 Build intuition using simple systems State estimation using observer Estimation using observer: ˆx = Aˆx + Bu + K ( y C ˆx ), ˆx 0 x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = / 21 dave.zachariah@it.uu.se

31 Build intuition using simple systems State estimation using observer Estimation using observer: ˆx = Aˆx + Bu + K ( y C ˆx ), ˆx 0 x x 2 (velocity) [m/s] x 1 (position) [m] x versus ˆx at t = / 21 dave.zachariah@it.uu.se

32 State estimation Estimation error and observability Estimation error: x x ˆx [Board: derive evolution of estimation errors] 13 / 21 dave.zachariah@it.uu.se

33 State estimation Estimation error and observability Estimation error: x x ˆx [Board: derive evolution of estimation errors] Result Errors of observer described as system x(t) = e (A KC)t x(0) and therefore x(t) decays at a rate given by maximum Re{ s i } where s i are observer poles/eigenvalues of (A KC). 13 / 21 dave.zachariah@it.uu.se

34 State estimation Estimation error and observability Estimation error: x x ˆx Result 9.2 State-space form is observable (cf. det O 0) matrix K can be chosen such that x vanish arbitrarily quick 14 / 21 dave.zachariah@it.uu.se

35 State estimation Estimation error and observability Estimation error: x x ˆx Result 9.2 State-space form is observable (cf. det O 0) matrix K can be chosen such that x vanish arbitrarily quick K is solved by polynomial det(si A + KC) = 0 with desired roots in left halfplane 14 / 21 dave.zachariah@it.uu.se

36 State estimation Estimation error and observability Estimation error: x x ˆx Result 9.2 State-space form is observable (cf. det O 0) matrix K can be chosen such that x vanish arbitrarily quick K is solved by polynomial det(si A + KC) = 0 with desired roots in left halfplane Quick observer ˆx is however sensitive to measurement noise! 14 / 21 dave.zachariah@it.uu.se

37 Combining feedback with estimated states 15 / 21 dave.zachariah@it.uu.se

38 Feedback using estimated states Controller the Laplace domain rl 0 + u x y (si A) 1 B C L ˆx Obs. 16 / 21 dave.zachariah@it.uu.se

39 Feedback using estimated states Controller the Laplace domain rl 0 + u x y (si A) 1 B C L ˆx Obs. System and controller with observer: { { ẋ = Ax + Bu u and y = Cx ˆx = Lˆx + l 0 r = Aˆx + Bu + K(y C ˆx) 16 / 21 dave.zachariah@it.uu.se

40 Feedback using estimated states Controller the Laplace domain rl 0 + u x y (si A) 1 B C L ˆx Obs. Controller with observer: { U(s) = L L : X(s) + l 0 R(s) s X(s) = A X(s) + BU(s) + K ( Y (s) C X(s) ) [Board: solve for controller] 16 / 21 dave.zachariah@it.uu.se

41 Feedback using estimated states General linear feedback control General linear feedback form, ch.9.5 G&L Controller with observer can be written as where U(s) = F r (s)r(s) F y (s)y (s), F r (s) = (1 L(sI A + KC + BL) 1 B)l 0 F y (s) = L(sI A + KC + BL) 1 K r F r + u G y F y 17 / 21 dave.zachariah@it.uu.se

42 Resulting closed-loop system 18 / 21 dave.zachariah@it.uu.se

43 Feedback using estimated states Effect of estimation error rl 0 + u x y (si A) 1 B C L ˆx Obs. Study system and controller with observer: { ẋ = Ax + Bu and u = Lˆx + l 0 r y = Cx by substituting ˆx = x x [Board: derive the closed-loop system with estimation error x] 19 / 21 dave.zachariah@it.uu.se

44 Feedback using estimated states Effect of estimation error rl 0 + u x y (si A) 1 B C L ˆx Obs. Yields closed-loop system: ẋ = (A BL)x + y = Cx with additional error states effect of estimation error {}}{ BL x +Bl 0 r x = (A KC) x 0 [Board: write the closed-loop system in state-space form] 19 / 21 dave.zachariah@it.uu.se

45 Feedback using estimated states The closed-loop system with observer The closed-loop system with estimation error can be written as [ẋ ] [ ] [ ] [ ] A BL BL x B = + l x 0 A KC x 0 0 r }{{}}{{} Ã B y = [ C 0 ] }{{} C with extended state vector. [ x x This yields transfer function from r to y: ] G c (s) = C(sI Ã) 1 B 20 / 21 dave.zachariah@it.uu.se

46 Feedback using estimated states The closed-loop system with observer Closed-loop system transfer function, ch.9.5 G&L Insert matrices Ã, B and C yields G c (s) = C(sI Ã) 1 B = C(sI A + BL) 1 Bl 0 with same poles as if states were known and K is gone! r F r + u G y F y 20 / 21 dave.zachariah@it.uu.se

47 Summary and recap Rules of thumb for pole placement Estimation using observer Feedback using estimated states Closed-loop system with observer 21 / 21 dave.zachariah@it.uu.se

Intro. Computer Control Systems: F8

Intro. Computer Control Systems: F8 Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2

More information

POLE PLACEMENT. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 19

POLE PLACEMENT. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 19 POLE PLACEMENT Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 19 Outline 1 State Feedback 2 Observer 3 Observer Feedback 4 Reduced Order

More information

CONTROL DESIGN FOR SET POINT TRACKING

CONTROL DESIGN FOR SET POINT TRACKING Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded

More information

5. Observer-based Controller Design

5. Observer-based Controller Design EE635 - Control System Theory 5. Observer-based Controller Design Jitkomut Songsiri state feedback pole-placement design regulation and tracking state observer feedback observer design LQR and LQG 5-1

More information

Lec 6: State Feedback, Controllability, Integral Action

Lec 6: State Feedback, Controllability, Integral Action Lec 6: State Feedback, Controllability, Integral Action November 22, 2017 Lund University, Department of Automatic Control Controllability and Observability Example of Kalman decomposition 1 s 1 x 10 x

More information

Topic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback

Topic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall

More information

Control Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft

Control Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Control Systems Design, SC426 SC426 Fall 2, dr A Abate, DCSC, TU Delft Lecture 5 Controllable Canonical and Observable Canonical Forms Stabilization by State Feedback State Estimation, Observer Design

More information

6 OUTPUT FEEDBACK DESIGN

6 OUTPUT FEEDBACK DESIGN 6 OUTPUT FEEDBACK DESIGN When the whole sate vector is not available for feedback, i.e, we can measure only y = Cx. 6.1 Review of observer design Recall from the first class in linear systems that a simple

More information

ECE 388 Automatic Control

ECE 388 Automatic Control Controllability and State Feedback Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage:

More information

EL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)

EL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1) EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the

More information

EEE582 Homework Problems

EEE582 Homework Problems EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use

More information

Linear State Feedback Controller Design

Linear State Feedback Controller Design Assignment For EE5101 - Linear Systems Sem I AY2010/2011 Linear State Feedback Controller Design Phang Swee King A0033585A Email: king@nus.edu.sg NGS/ECE Dept. Faculty of Engineering National University

More information

Control Systems Design

Control Systems Design ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:

More information

ECEEN 5448 Fall 2011 Homework #4 Solutions

ECEEN 5448 Fall 2011 Homework #4 Solutions ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with

More information

EL2520 Control Theory and Practice

EL2520 Control Theory and Practice So far EL2520 Control Theory and Practice r Fr wu u G w z n Lecture 5: Multivariable systems -Fy Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden SISO control revisited: Signal

More information

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state

More information

State Feedback and State Estimators Linear System Theory and Design, Chapter 8.

State Feedback and State Estimators Linear System Theory and Design, Chapter 8. 1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 2 Homework 7: State Estimators (a) For the same system as discussed in previous slides, design another closed-loop state estimator,

More information

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7) EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and

More information

One-Sided Laplace Transform and Differential Equations

One-Sided Laplace Transform and Differential Equations One-Sided Laplace Transform and Differential Equations As in the dcrete-time case, the one-sided transform allows us to take initial conditions into account. Preliminaries The one-sided Laplace transform

More information

Pole placement control: state space and polynomial approaches Lecture 2

Pole placement control: state space and polynomial approaches Lecture 2 : state space and polynomial approaches Lecture 2 : a state O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename -based November 21, 2017 Outline : a state

More information

Identification Methods for Structural Systems

Identification Methods for Structural Systems Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from

More information

EE C128 / ME C134 Fall 2014 HW 9 Solutions. HW 9 Solutions. 10(s + 3) s(s + 2)(s + 5) G(s) =

EE C128 / ME C134 Fall 2014 HW 9 Solutions. HW 9 Solutions. 10(s + 3) s(s + 2)(s + 5) G(s) = 1. Pole Placement Given the following open-loop plant, HW 9 Solutions G(s) = 1(s + 3) s(s + 2)(s + 5) design the state-variable feedback controller u = Kx + r, where K = [k 1 k 2 k 3 ] is the feedback

More information

TRACKING AND DISTURBANCE REJECTION

TRACKING AND DISTURBANCE REJECTION TRACKING AND DISTURBANCE REJECTION Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 13 General objective: The output to track a reference

More information

Observability and state estimation

Observability and state estimation EE263 Autumn 2015 S Boyd and S Lall Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time observability

More information

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

Introduction to Computer Control Systems

Introduction to Computer Control Systems Introduction to Computer Control Systems Lecture 1: Introduction Dave Zachariah Div. Systems and Control, Dept. Information Technology, Uppsala University October 28, 2014 (UU/Info Technology/SysCon) Intro.

More information

Lecture 19 Observability and state estimation

Lecture 19 Observability and state estimation EE263 Autumn 2007-08 Stephen Boyd Lecture 19 Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time

More information

Control Systems. Frequency domain analysis. L. Lanari

Control Systems. Frequency domain analysis. L. Lanari Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic

More information

Solution via Laplace transform and matrix exponential

Solution via Laplace transform and matrix exponential EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative

More information

Final: Signal, Systems and Control (BME )

Final: Signal, Systems and Control (BME ) Final: Signal, Systems and Control (BME 580.) Instructor: René Vidal May 0th 007 HONOR SYSTEM: This examination is strictly individual. You are not allowed to talk, discuss, exchange solutions, etc., with

More information

State Feedback and State Estimators Linear System Theory and Design, Chapter 8.

State Feedback and State Estimators Linear System Theory and Design, Chapter 8. 1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 State Estimator In previous section, we have discussed the state feedback, based on the assumption that all state variables are

More information

Time Response of Systems

Time Response of Systems Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =

More information

Chapter 3. State Feedback - Pole Placement. Motivation

Chapter 3. State Feedback - Pole Placement. Motivation Chapter 3 State Feedback - Pole Placement Motivation Whereas classical control theory is based on output feedback, this course mainly deals with control system design by state feedback. This model-based

More information

SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER

SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER Exercise 54 Consider the system: ẍ aẋ bx u where u is the input and x the output signal (a): Determine a state space realization (b): Is the

More information

Control Systems. Laplace domain analysis

Control Systems. Laplace domain analysis Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output

More information

CDS Solutions to the Midterm Exam

CDS Solutions to the Midterm Exam CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2

More information

Zeros and zero dynamics

Zeros and zero dynamics CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)

More information

Synthesis via State Space Methods

Synthesis via State Space Methods Chapter 18 Synthesis via State Space Methods Here, we will give a state space interpretation to many of the results described earlier. In a sense, this will duplicate the earlier work. Our reason for doing

More information

ME 132, Fall 2015, Quiz # 2

ME 132, Fall 2015, Quiz # 2 ME 132, Fall 2015, Quiz # 2 # 1 # 2 # 3 # 4 # 5 # 6 Total NAME 14 10 8 6 14 8 60 Rules: 1. 2 sheets of notes allowed, 8.5 11 inches. Both sides can be used. 2. Calculator is allowed. Keep it in plain view

More information

1 (30 pts) Dominant Pole

1 (30 pts) Dominant Pole EECS C8/ME C34 Fall Problem Set 9 Solutions (3 pts) Dominant Pole For the following transfer function: Y (s) U(s) = (s + )(s + ) a) Give state space description of the system in parallel form (ẋ = Ax +

More information

Control Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft

Control Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Control Systems Design, SC4026 SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) and Observability Algebraic Tests (Kalman rank condition & Hautus test) A few

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

SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL

SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL 1. Optimal regulator with noisy measurement Consider the following system: ẋ = Ax + Bu + w, x(0) = x 0 where w(t) is white noise with Ew(t) = 0, and x 0 is a stochastic

More information

Controls Problems for Qualifying Exam - Spring 2014

Controls Problems for Qualifying Exam - Spring 2014 Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function

More information

There are none. Abstract for Gauranteed Margins for LQG Regulators, John Doyle, 1978 [Doy78].

There are none. Abstract for Gauranteed Margins for LQG Regulators, John Doyle, 1978 [Doy78]. Chapter 7 Output Feedback There are none. Abstract for Gauranteed Margins for LQG Regulators, John Doyle, 1978 [Doy78]. In the last chapter we considered the use of state feedback to modify the dynamics

More information

EE C128 / ME C134 Final Exam Fall 2014

EE C128 / ME C134 Final Exam Fall 2014 EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket

More information

University of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.

University of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA. University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is

More information

Chapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control

Chapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design

More information

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017

More information

ECE504: Lecture 9. D. Richard Brown III. Worcester Polytechnic Institute. 04-Nov-2008

ECE504: Lecture 9. D. Richard Brown III. Worcester Polytechnic Institute. 04-Nov-2008 ECE504: Lecture 9 D. Richard Brown III Worcester Polytechnic Institute 04-Nov-2008 Worcester Polytechnic Institute D. Richard Brown III 04-Nov-2008 1 / 38 Lecture 9 Major Topics ECE504: Lecture 9 We are

More information

1 Continuous-time Systems

1 Continuous-time Systems Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output

More information

State will have dimension 5. One possible choice is given by y and its derivatives up to y (4)

State will have dimension 5. One possible choice is given by y and its derivatives up to y (4) A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To

More information

Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science : MULTIVARIABLE CONTROL SYSTEMS by A.

Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science : MULTIVARIABLE CONTROL SYSTEMS by A. Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Q-Parameterization 1 This lecture introduces the so-called

More information

DESIGN OF OBSERVERS FOR SYSTEMS WITH SLOW AND FAST MODES

DESIGN OF OBSERVERS FOR SYSTEMS WITH SLOW AND FAST MODES DESIGN OF OBSERVERS FOR SYSTEMS WITH SLOW AND FAST MODES by HEONJONG YOO A thesis submitted to the Graduate School-New Brunswick Rutgers, The State University of New Jersey In partial fulfillment of the

More information

Control Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft

Control Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) vs Observability Algebraic Tests (Kalman rank condition & Hautus test) A few

More information

Perspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems

Perspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems ECE 3640 Lecture State-Space Analysis Objective: systems To learn about state-space analysis for continuous and discrete-time Perspective Transfer functions provide only an input/output perspective of

More information

Model reduction via tangential interpolation

Model reduction via tangential interpolation Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time

More information

LMIs for Observability and Observer Design

LMIs for Observability and Observer Design LMIs for Observability and Observer Design Matthew M. Peet Arizona State University Lecture 06: LMIs for Observability and Observer Design Observability Consider a system with no input: ẋ(t) = Ax(t), x(0)

More information

Topic # Feedback Control

Topic # Feedback Control Topic #5 6.3 Feedback Control State-Space Systems Full-state Feedback Control How do we change the poles of the state-space system? Or,evenifwecanchangethepolelocations. Where do we put the poles? Linear

More information

Contents lecture 5. Automatic Control III. Summary of lecture 4 (II/II) Summary of lecture 4 (I/II) u y F r. Lecture 5 H 2 and H loop shaping

Contents lecture 5. Automatic Control III. Summary of lecture 4 (II/II) Summary of lecture 4 (I/II) u y F r. Lecture 5 H 2 and H loop shaping Contents lecture 5 Automatic Control III Lecture 5 H 2 and H loop shaping Thomas Schön Division of Systems and Control Department of Information Technology Uppsala University. Email: thomas.schon@it.uu.se,

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

EC Control Engineering Quiz II IIT Madras

EC Control Engineering Quiz II IIT Madras EC34 - Control Engineering Quiz II IIT Madras Linear algebra Find the eigenvalues and eigenvectors of A, A, A and A + 4I Find the eigenvalues and eigenvectors of the following matrices: (a) cos θ sin θ

More information

BALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS

BALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007

More information

Systems Engineering/Process Control L4

Systems Engineering/Process Control L4 1 / 24 Systems Engineering/Process Control L4 Input-output models Laplace transform Transfer functions Block diagram algebra Reading: Systems Engineering and Process Control: 4.1 4.4 2 / 24 Laplace transform

More information

Pole Placement Design

Pole Placement Design Department of Automatic Control LTH, Lund University 1 Introduction 2 Simple Examples 3 Polynomial Design 4 State Space Design 5 Robustness and Design Rules 6 Model Reduction 7 Oscillatory Systems 8 Summary

More information

D(s) G(s) A control system design definition

D(s) G(s) A control system design definition R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure

More information

Linear Systems. Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output

Linear Systems. Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output Linear Systems Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output Our interest is in dynamic systems Dynamic system means a system with memory of course including

More information

Full State Feedback for State Space Approach

Full State Feedback for State Space Approach Full State Feedback for State Space Approach State Space Equations Using Cramer s rule it can be shown that the characteristic equation of the system is : det[ si A] 0 Roots (for s) of the resulting polynomial

More information

Modelling and Control of Dynamic Systems. Stability of Linear Systems. Sven Laur University of Tartu

Modelling and Control of Dynamic Systems. Stability of Linear Systems. Sven Laur University of Tartu Modelling and Control of Dynamic Systems Stability of Linear Systems Sven Laur University of Tartu Motivating Example Naive open-loop control r[k] Controller Ĉ[z] u[k] ε 1 [k] System Ĝ[z] y[k] ε 2 [k]

More information

2E1252 Control Theory and Practice

2E1252 Control Theory and Practice 2E1252 Control Theory and Practice Lectre 11: Actator satration and anti wind-p Learning aims After this lectre, yo shold nderstand how satration can case controller states to wind p know how to modify

More information

Topics in control Tracking and regulation A. Astolfi

Topics in control Tracking and regulation A. Astolfi Topics in control Tracking and regulation A. Astolfi Contents 1 Introduction 1 2 The full information regulator problem 3 3 The FBI equations 5 4 The error feedback regulator problem 5 5 The internal model

More information

Course Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10)

Course Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10) Amme 35 : System Dynamics Control State Space Design Course Outline Week Date Content Assignment Notes 1 1 Mar Introduction 2 8 Mar Frequency Domain Modelling 3 15 Mar Transient Performance and the s-plane

More information

Advanced Control Theory

Advanced Control Theory State Feedback Control Design chibum@seoultech.ac.kr Outline State feedback control design Benefits of CCF 2 Conceptual steps in controller design We begin by considering the regulation problem the task

More information

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries . AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace

More information

Design of Positive Linear Observers for Positive Systems via Coordinates Transformation and Positive Realization

Design of Positive Linear Observers for Positive Systems via Coordinates Transformation and Positive Realization Design of Positive Linear Observers for Positive Systems via Coordinates Transformation and Positive Realization Juhoon Back 1 Alessandro Astolfi 1 1 Department of Electrical and Electronic Engineering

More information

Chapter 9 Observers, Model-based Controllers 9. Introduction In here we deal with the general case where only a subset of the states, or linear combin

Chapter 9 Observers, Model-based Controllers 9. Introduction In here we deal with the general case where only a subset of the states, or linear combin Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 9 Observers,

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

LECTURE NOTES ON CONTROL

LECTURE NOTES ON CONTROL Department of Control for Transportation and Vehicle Systems Faculty of Transportation Engineering and Vehicle Engineering Budapest University of Technology and Economics Tamás Tettamanti PhD., Qiong Lu

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

Transfer Functions. Chapter Introduction. 6.2 The Transfer Function

Transfer Functions. Chapter Introduction. 6.2 The Transfer Function Chapter 6 Transfer Functions As a matter of idle curiosity, I once counted to find out what the order of the set of equations in an amplifier I had just designed would have been, if I had worked with the

More information

= m(0) + 4e 2 ( 3e 2 ) 2e 2, 1 (2k + k 2 ) dt. m(0) = u + R 1 B T P x 2 R dt. u + R 1 B T P y 2 R dt +

= m(0) + 4e 2 ( 3e 2 ) 2e 2, 1 (2k + k 2 ) dt. m(0) = u + R 1 B T P x 2 R dt. u + R 1 B T P y 2 R dt + ECE 553, Spring 8 Posted: May nd, 8 Problem Set #7 Solution Solutions: 1. The optimal controller is still the one given in the solution to the Problem 6 in Homework #5: u (x, t) = p(t)x k(t), t. The minimum

More information

ECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77

ECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77 1/77 ECEN 605 LINEAR SYSTEMS Lecture 7 Solution of State Equations Solution of State Space Equations Recall from the previous Lecture note, for a system: ẋ(t) = A x(t) + B u(t) y(t) = C x(t) + D u(t),

More information

EECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.

EECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators. Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 5-8 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3

More information

Internal Model Principle

Internal Model Principle Internal Model Principle If the reference signal, or disturbance d(t) satisfy some differential equation: e.g. d n d dt n d(t)+γ n d d 1 dt n d 1 d(t)+...γ d 1 dt d(t)+γ 0d(t) =0 d n d 1 then, taking Laplace

More information

State Variable Analysis of Linear Dynamical Systems

State Variable Analysis of Linear Dynamical Systems Chapter 6 State Variable Analysis of Linear Dynamical Systems 6 Preliminaries In state variable approach, a system is represented completely by a set of differential equations that govern the evolution

More information

Introduction to Modern Control MT 2016

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

EEE 184: Introduction to feedback systems

EEE 184: Introduction to feedback systems EEE 84: Introduction to feedback systems Summary 6 8 8 x 7 7 6 Level() 6 5 4 4 5 5 time(s) 4 6 8 Time (seconds) Fig.. Illustration of BIBO stability: stable system (the input is a unit step) Fig.. step)

More information

Section 5 Dynamics and Control of DC-DC Converters

Section 5 Dynamics and Control of DC-DC Converters Section 5 Dynamics and ontrol of D-D onverters 5.2. Recap on State-Space Theory x Ax Bu () (2) yxdu u v d ; y v x2 sx () s Ax() s Bu() s ignoring x (0) (3) ( si A) X( s) Bu( s) (4) X s si A BU s () ( )

More information

Linearization problem. The simplest example

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

Linear dynamical systems with inputs & outputs

Linear dynamical systems with inputs & outputs EE263 Autumn 215 S. Boyd and S. Lall Linear dynamical systems with inputs & outputs inputs & outputs: interpretations transfer function impulse and step responses examples 1 Inputs & outputs recall continuous-time

More information

EECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.

EECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators. Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total

More information

State estimation and the Kalman filter

State estimation and the Kalman filter State estimation and the Kalman filter PhD, David Di Ruscio Telemark university college Department of Technology Systems and Control Engineering N-3914 Porsgrunn, Norway Fax: +47 35 57 52 50 Tel: +47 35

More information

Exam in Automatic Control II Reglerteknik II 5hp (1RT495)

Exam in Automatic Control II Reglerteknik II 5hp (1RT495) Exam in Automatic Control II Reglerteknik II 5hp (1RT495) Date: August 4, 018 Venue: Bergsbrunnagatan 15 sal Responsible teacher: Hans Rosth. Aiding material: Calculator, mathematical handbooks, textbooks

More information

Control engineering sample exam paper - Model answers

Control engineering sample exam paper - Model answers Question Control engineering sample exam paper - Model answers a) By a direct computation we obtain x() =, x(2) =, x(3) =, x(4) = = x(). This trajectory is sketched in Figure (left). Note that A 2 = I

More information

Power System Control

Power System Control Power System Control Basic Control Engineering Prof. Wonhee Kim School of Energy Systems Engineering, Chung-Ang University 2 Contents Why feedback? System Modeling in Frequency Domain System Modeling in

More information

Time Response Analysis (Part II)

Time Response Analysis (Part II) Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary

More information

9.5 The Transfer Function

9.5 The Transfer Function Lecture Notes on Control Systems/D. Ghose/2012 0 9.5 The Transfer Function Consider the n-th order linear, time-invariant dynamical system. dy a 0 y + a 1 dt + a d 2 y 2 dt + + a d n y 2 n dt b du 0u +

More information

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics

More information

Steady State Kalman Filter

Steady State Kalman Filter Steady State Kalman Filter Infinite Horizon LQ Control: ẋ = Ax + Bu R positive definite, Q = Q T 2Q 1 2. (A, B) stabilizable, (A, Q 1 2) detectable. Solve for the positive (semi-) definite P in the ARE:

More information

Observability. It was the property in Lyapunov stability which allowed us to resolve that

Observability. It was the property in Lyapunov stability which allowed us to resolve that Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property

More information