Lifted approach to ILC/Repetitive Control

Size: px
Start display at page:

Download "Lifted approach to ILC/Repetitive Control"

Transcription

1 Lifted approach to ILC/Repetitive Control Okko H. Bosgra Maarten Steinbuch TUD Delft Centre for Systems and Control TU/e Control System Technology Dutch Institute of Systems and Control DISC winter semester 2003/2004

2 Contents Formulation lifted approach Internal model principle Design of ILC without Q-filter using LQ theory Analysis of existing design results Analysis of Repetitive Control Properties, convergence, performance

3 classical ILC structure 1 L(z) Q(z) f kn C(z) z N f k P (z) y ref ε k y k d Update in trial space: f kn = Q(z)(f k L(z)ε k )

4 classical ILC structure 2 y ref, d assumed to be periodic inputs, period length N ε k f k d y ref C(z) P (z) y k control input f k, disturbance input y ref d, output ε k input/output relationship: ε k = (I P C) 1 (y ref d) (I P C) 1 P fk

5 System formulation 1 y ref d S(z) f k P (z) ε k Define P (z) : = (I P C) 1 P 1 S(z) : = (I P C)

6 System formulation 2 Better to represent these relations in the time-domain: (A p, B p, C p ) minimal realization for P (z), (A c, B c, C c, D c ) minimal realization for C(z): P (z) = C p (zi A p ) 1 B p C(z) = C c (zi A c ) 1 B c D c i.e. P (z) strictly proper, C(z) proper. State space representation for joint feedback system: x k1 = Ax k Bf k N r k where r k = y ref d ε k = Cx k r k

7 System formulation 3 ( ) ( ) ( ) A = C = ( A p B p D c C p B p C c B c C p A c ) C p 0 N = B p D c B c B = B p 0 r k f k A, B, C, N ε k

8 System formulation 4 Split up of system: r k periodic, r k asymptotically periodic r k A, N, C f k A, B, C r k ε k Simplification: consider r k as being periodic

9 System formulation 5 r k f k A, B, C ε k Servo problem ILC: x 0 = 0, x N = 0, x 2N = 0, x 3N = 0,... typical pattern of r k : pick and place r k 0 N 2N 3N time

10 System formulation 6 Periodic disturbance suppression Repetitive Control typical pattern of r k : r k 0 N 2N 3N time e.g. drive control with nonuniform torque pattern x 0 = 0, x N, x 2N, x 3N,... follow from RC algorithm Thus ILC and RC require slightly different formulations

11 System formulation in trial domain 1 Causal LTI system (A, B, C) yields impulse response h k R m l h k = CA k 1 B k = 1, 2, 3,... = 0 k 0 The system has m outputs, l inputs Let j = 0, 1, 2,... denote the trial number and define y j = y Nj y Nj1 y Nj2. f j = f Nj f Nj1 f Nj2. r = r 0 r 1 r 2. ε j = ε Nj ε Nj1 ε Nj2. y NjN 1 f NjN 1 r N 1 ε NjN 1

12 System formulation in trial domain 2 The lifted system representation in trial domain is now: x NjN = F x Nj Gf j y j = Hx Nj Jf j ε j = y j r where h h 1 h h 2 h 1 h J = h N h 2 h 1 h 0 ( ) G = A N 1 B A N 2 B... AB B C CA H = CA 2. F = A N CA N 1

13 System formulation in trial domain 3 time invariance J is Toeplitz causality J is lower triangular Toeplitz matrix J represents convolution (A, B, C) minimal (F, G, H) minimal In trial domain, r is constant disturbance vector For ILC, x k = 0 for k = Nj, j = 0, 1, 2,... system model reduces to ε j = Jf j r

14 Internal model principle 1 Asymptotic rejection of constant disturbance provided that 1. Error-driven dynamic disturbance model is added to controller 2. Each controller disturbance mode can propagate throught plant, i.e. is not cancelled against transmission zero 3. Dynamics of plant and controller are asymptotically stabilised by feedback Vector-valued constant disturbance r generated in trial domain by integrator having initial value ic = r z 1 r ic

15 Internal model principle 2 Adding disturbance model in feedback loop with gain L: r f j1 f j x j1 x j y j ε j z 1 G z 1 H F J L Asymptotic rejection: ε j 0 for j while r = constant 0 Then y j must compensate r y j r for j

16 Internal model principle 3 Only possible if rank of steady-state gain of (F, G, H, J) is Nm and if f j has at least dimension Nm Closed-loop state-space model in trial domain: ( ) ( ) ( ) ( ) f j1 I Nl LJ LH f j L = r G F 0 x j1 ε j = ( ) ( ) f j J H r Closed-loop system matrix: ( ) ( ) ( I Nl LJ LH I Nl 0 = G F G F x j x j I Nl 0 ) ( ) L J H

17 Internal model principle 4 System is controllable if controllability matrix has rank N l n: [ ] I Nl I Nl I Nl I Nl... rank = 0 G G F G G F G F 2 G... [ ] I Nl rank = n Nl 0 G F G F 2 G... System is observable if observability matrix has rank N l n: J H I Nm 0 [ ] J HG HF rank J HG HF G HF 2 = rank I Nm H J H I Nm H F H G F I n.... This requires both matrices to have rank at least Nl n:

18 Internal model principle 5 I Nm 0 I Nm 0 0 H I rank Nm H I Nm H F H = rank 0 F H = Nm n 0 F 2 H.... This requires m l, i.e. number of outputs number of inputs [ ] [ ] J H J H(F I n ) 1 rank = rank = G F I n G I n [ ] J H(I n F ) 1 G 0 rank = nrank (J H(zI n F ) 1 G) 0 I z=1 n Nm Nl plant steady-state gain matrix must have rank Nl

19 Internal model principle 6 If (F, G, H) is of order zero (ILC case) then the Nm Nl matrix J must have rank Nl, i.e. its columns must be linearly independent If (F, G, H) is of order greater than zero (RC case) then we may assume F = A N 0 Then the requirement is rank(j HG) = Nl where C h N h 3 h 2 h 1 CA. ( ) h3 h 2 HG = CA 2 A N 1 B... AB B = h CA N 1 h 2N 1 Contribution to rank from upper right influence of previous trial h N

20 Internal model principle 7 Result: Suppose L asymptotically stabilizes system. Then disturbance r is asymptotically rejected if m l. Proof: L asymptotically stabilizes all dynamics ( ) λ I Nl LJ LH i G F < 1 i = 1, 2,..., Nl n which implies ( ) LJ LH rank = Nl n G I n F i.e. this (Nl n) (Nl n) matrix is invertible For asymptotically stable transfer function matrix P (z), the steady state gain matrix is P (z) z=1

21 Internal model principle 8 Steady state gain matrix between r and ε j : ( ) [ ( )] 1 ( ) I Nl LJ LH L I Nm J H I Nln = G F 0 ( ) ( ) 1 ( ) LJ LH L = I Nm J H =: Ω G I n F 0 The rank of this steady state gain matrix Ω follows from I Nm J H I Nm Γ 1 Γ 2 rank L LJ LH = rank L I Nl 0 0 G I n F 0 0 I n where ( ) ( ) 1 LJ LH Γ = J H G I n F

22 The rank of this expression now equals Ω 0 0 rank 0 I Nl 0 = rank Ω Nl n 0 0 I n Internal model principle 9 From its structure, the rank of the first matrix is I Nm J H rank L LJ LH Nm n 0 G I n F so that rank Ω Nm Nl There is asymptotic rejection of r if this rank is zero, i.e. for m l, or number of inputs number of outputs

23 Internal model principle: conclusion 10 J has dimension Nm Nl. Feedback L requires: system controllable: always system observable if m l and rank(j HG) = Nl i.e. HG can contribute in making J full rank disturbance rejection if m l Thus requirements: m = l and (J HG) square and full rank If not satisfied: not all modes asymptotically stable r not fully compensated

24 Iterative learning control 1 For ILC, the initial state in each trial is zero, which applies in machine operations like pick-and-place tasks. Thus the system is r f j1 f j y j ε j z 1 J L

25 Iterative learning control 2 Controller contains all dynamics in system no Q-filter thus uncompromised convergence J square, not necessarily of full rank convergence if eigenvalues of I LJ smaller than 1 faster convergence for smaller eigenvalues of I LJ L can be time-varying, non-causal, i.e. no reason to restrict L to be Toeplitz additive noise disturbances may act on J If J full rank, feedback system is first order with state feedback good robustness properties attainable

26 Iterative learning control 3 System relations: f j1 = (I Nl LJ)f j Lr (1) ε j = Jf j r (2) Asymptotically stable if λ i (I Nl LJ) < 1 i = 1, 2,..., Nl This implies that LJ is non-singular, so that non-singularity of LJ is necessary in order that the system is asymptotically stable. Transfer function matrix betreen r and ε j (sensitivity function) ε j = [I Nl J(zI Nl I Nl LJ) 1 L]r For z = 1 we have steady state gain I Nl J( LJ) 1 L

27 Iterative learning control 4 Thus LJ non-singular and steady-state gain zero requires J and L to be square and non-singular If J singular: stabilization and full rejection of r not possible Two reasons for loss of rank for J: strictly proper system, or additional delays non-minimum phase zeroes in i/o behaviour of system If Markov parameters h 0, h 1, h 2,..., h d 1 are zero, then rank J is at most Nm d In that case: no full stabilization, no full rejection

28 ILC example r = im J = im 1 0 = im J = J is output matrix in system with unit system matrix third mode unobservable feedback has no effect on pole location of third mode

29 ILC example 2.1 H(z) = z 2 z(z 1) = 2 z 3 z 1 [ ] zi A B = C zero at z = impulse response h = z z observability matrix H =

30 ILC example J = By definition of zeros there exists x 0 and exponential input signal containing zero as exponential factor, such that resulting output is zero: ( ) =

31 ILC example 2.3 or J f j H x 0 = 0 Then upper bound on smallest non-zero singular value of J is given by σ(j) < H x 0 f j For nonminimum-phase zero f j will grow for larger dimensions N forcing σ(j) to be small Thus delays and non-minimum phase zeroes create loss of rank for J

32 ILC singular values 1 If J is square but rank-deficient, the singular value decomposition provides the rank as the dimension of Σ 1 : ( ) ( ) ( ) Σ 1 0 V1 T J = U 1 U 2 0 Σ 2 V2 T where Σ 2 is zero. Then decompose f, the input to J, as where Jf = Jf 1 Jf 2 f 1 imv 1, f 2 imv 2 Thus Jf = U 1 Σ 1 V T 1 f 1 as V T 1 f 2 = 0, V T 2 f 1 = 0

33 and replace f 1 by V 1 f, i.e dimension reduction of state ILC singular values 2 r u j f j1 f j y j ε j z 1 JV 1 L All poles can now be stabilised or assigned

34 ILC pole assignment f j1 = (I LJV 1 )f k Lr f 0 = 0 ε j = JV 1 f j r If or then LJV 1 = α I LU 1 Σ 1 V T 1 V 1 = α I L = ασ 1 1 U T 1 If α = 1 then dead-beat. The output y follows r only in subspace im J =im U 1

35 Alternative for pole assignment: LQ optimal control f j1 = f j u j ILC LQ optimal control 1 y j = J V 1 f j LQ criterion: Cr = = (yk T Q y k u T k R u k ) k=1 k=1 (f T k V T 1 J T QJV 1 f k u T k R u k ) Useful choice: Q = I, R = βi

36 ILC LQ optimal control 2 Solution to LQ-optimal control problem with unit system and input matrices: u j = (βi P ) 1 P f j where P is stabilizing solution of algebraic Riccati equation or P = P V T 1 J T JV 1 P (βi P ) 1 P 0 = Σ 2 1 P (βi P ) 1 P Thus P is diagonal, with entries p i on diagonal, and Σ 1 = diag (σ i ): ( ) p i = 1 2 σ2 i 1 1 4β σ 2 σi 2 i β

37 ILC LQ optimal control 3 u j = (βi Σ 2 1 βi) 1 (Σ 2 1 βi) f j = LJV 1 f j Solving for L gives the result L = (2βI Σ 2 1) 1 (Σ 1 βσ 1 )U T 1

Iterative Learning Control Analysis and Design I

Iterative Learning Control Analysis and Design I Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations

More information

Video 6.1 Vijay Kumar and Ani Hsieh

Video 6.1 Vijay Kumar and Ani Hsieh Video 6.1 Vijay Kumar and Ani Hsieh Robo3x-1.6 1 In General Disturbance Input + - Input Controller + + System Output Robo3x-1.6 2 Learning Objectives for this Week State Space Notation Modeling in the

More information

Analysis of Discrete-Time Systems

Analysis of Discrete-Time Systems TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time

More information

Analysis of Discrete-Time Systems

Analysis of Discrete-Time Systems TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time

More information

Design strategies for iterative learning control based on optimal control

Design strategies for iterative learning control based on optimal control Selected Topics in Signals, Systems and Control Vol. 2, September 2 Design strategies for iterative learning control based on optimal control Rob Tousain, Eduard van der Meché and Okko Bosgra Mechanical

More information

Here represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.

Here represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities. 19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that

More information

Iterative Learning Control (ILC)

Iterative Learning Control (ILC) Department of Automatic Control LTH, Lund University ILC ILC - the main idea Time Domain ILC approaches Stability Analysis Example: The Milk Race Frequency Domain ILC Example: Marine Vibrator Material:

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

High Precision Control of Ball Screw Driven Stage Using Repetitive Control with Sharp Roll-off Learning Filter

High Precision Control of Ball Screw Driven Stage Using Repetitive Control with Sharp Roll-off Learning Filter High Precision Control of Ball Screw Driven Stage Using Repetitive Control with Sharp Roll-off Learning Filter Tadashi Takemura and Hiroshi Fujimoto The University of Tokyo --, Kashiwanoha, Kashiwa, Chiba,

More information

Part 1: Introduction to the Algebraic Approach to ILC

Part 1: Introduction to the Algebraic Approach to ILC IEEE ICMA 2006 Tutorial Workshop: Control Algebraic Analysis and Optimal Design Presenters: Contributor: Kevin L. Moore Colorado School of Mines YangQuan Chen Utah State University Hyo-Sung Ahn ETRI, Korea

More information

Chapter 13 Digital Control

Chapter 13 Digital Control Chapter 13 Digital Control Chapter 12 was concerned with building models for systems acting under digital control. We next turn to the question of control itself. Topics to be covered include: why one

More information

On Spectral Factorization and Riccati Equations for Time-Varying Systems in Discrete Time

On Spectral Factorization and Riccati Equations for Time-Varying Systems in Discrete Time On Spectral Factorization and Riccati Equations for Time-Varying Systems in Discrete Time Alle-Jan van der Veen and Michel Verhaegen Delft University of Technology Department of Electrical Engineering

More information

Design Methods for Control Systems

Design Methods for Control Systems Design Methods for Control Systems Maarten Steinbuch TU/e Gjerrit Meinsma UT Dutch Institute of Systems and Control Winter term 2002-2003 Schedule November 25 MSt December 2 MSt Homework # 1 December 9

More information

Contents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31

Contents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31 Contents Preamble xiii Linear Systems I Basic Concepts 1 I System Representation 3 1 State-Space Linear Systems 5 1.1 State-Space Linear Systems 5 1.2 Block Diagrams 7 1.3 Exercises 11 2 Linearization

More information

CBE507 LECTURE III Controller Design Using State-space Methods. Professor Dae Ryook Yang

CBE507 LECTURE III Controller Design Using State-space Methods. Professor Dae Ryook Yang CBE507 LECTURE III Controller Design Using State-space Methods Professor Dae Ryook Yang Fall 2013 Dept. of Chemical and Biological Engineering Korea University Korea University III -1 Overview States What

More information

magnitude [db] phase [deg] frequency [Hz] feedforward motor load -

magnitude [db] phase [deg] frequency [Hz] feedforward motor load - ITERATIVE LEARNING CONTROL OF INDUSTRIAL MOTION SYSTEMS Maarten Steinbuch and René van de Molengraft Eindhoven University of Technology, Faculty of Mechanical Engineering, Systems and Control Group, P.O.

More information

APPROXIMATE REALIZATION OF VALVE DYNAMICS WITH TIME DELAY

APPROXIMATE REALIZATION OF VALVE DYNAMICS WITH TIME DELAY APPROXIMATE REALIZATION OF VALVE DYNAMICS WITH TIME DELAY Jan van Helvoirt,,1 Okko Bosgra, Bram de Jager Maarten Steinbuch Control Systems Technology Group, Mechanical Engineering Department, Technische

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

System Identification Using a Retrospective Correction Filter for Adaptive Feedback Model Updating

System Identification Using a Retrospective Correction Filter for Adaptive Feedback Model Updating 9 American Control Conference Hyatt Regency Riverfront, St Louis, MO, USA June 1-1, 9 FrA13 System Identification Using a Retrospective Correction Filter for Adaptive Feedback Model Updating M A Santillo

More information

Chapter 7 Interconnected Systems and Feedback: Well-Posedness, Stability, and Performance 7. Introduction Feedback control is a powerful approach to o

Chapter 7 Interconnected Systems and Feedback: Well-Posedness, Stability, and Performance 7. Introduction Feedback control is a powerful approach to o 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 7 Interconnected

More information

Iterative Learning Control for Tailor Rolled Blanks Zhengyang Lu DCT

Iterative Learning Control for Tailor Rolled Blanks Zhengyang Lu DCT Iterative Learning Control for Tailor Rolled Blans Zhengyang Lu DCT 2007.066 Master Internship Report Supervisors: Prof. Maarten Steinbuch (TU/e) Camile Hol (Corus) Eindhoven University of Technology Department

More information

Statistical and Adaptive Signal Processing

Statistical and Adaptive Signal Processing r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory

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

Lecture 2. Linear Systems

Lecture 2. Linear Systems Lecture 2. Linear Systems Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering April 6, 2015 1 / 31 Logistics hw1 due this Wed, Apr 8 hw due every Wed in class, or my mailbox on

More information

Design of iterative learning control algorithms using a repetitive process setting and the generalized KYP lemma

Design of iterative learning control algorithms using a repetitive process setting and the generalized KYP lemma Design of iterative learning control algorithms using a repetitive process setting and the generalized KYP lemma 22 th July 2015, Dalian, China Wojciech Paszke Institute of Control and Computation Engineering,

More information

An LQ R weight selection approach to the discrete generalized H 2 control problem

An LQ R weight selection approach to the discrete generalized H 2 control problem INT. J. CONTROL, 1998, VOL. 71, NO. 1, 93± 11 An LQ R weight selection approach to the discrete generalized H 2 control problem D. A. WILSON², M. A. NEKOUI² and G. D. HALIKIAS² It is known that a generalized

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

Process Modelling, Identification, and Control

Process Modelling, Identification, and Control Jan Mikles Miroslav Fikar 2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. Process Modelling, Identification, and

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

Subspace Identification Methods

Subspace Identification Methods Czech Technical University in Prague Faculty of Electrical Engineering Department of Control Engineering Subspace Identification Methods Technical Report Ing. Pavel Trnka Supervisor: Prof. Ing. Vladimír

More information

21.4. Engineering Applications of z-transforms. Introduction. Prerequisites. Learning Outcomes

21.4. Engineering Applications of z-transforms. Introduction. Prerequisites. Learning Outcomes Engineering Applications of z-transforms 21.4 Introduction In this Section we shall apply the basic theory of z-transforms to help us to obtain the response or output sequence for a discrete system. This

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

Kalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and

Kalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1

More information

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII CONTENTS VOLUME VII Control of Linear Multivariable Systems 1 Katsuhisa Furuta,Tokyo Denki University, School of Science and Engineering, Ishizaka, Hatoyama, Saitama, Japan 1. Linear Multivariable Systems

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 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:

More information

EL2520 Control Theory and Practice

EL2520 Control Theory and Practice EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller

More information

RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK

RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK TRNKA PAVEL AND HAVLENA VLADIMÍR Dept of Control Engineering, Czech Technical University, Technická 2, 166 27 Praha, Czech Republic mail:

More information

Problem Set 5 Solutions 1

Problem Set 5 Solutions 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel

More information

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors /88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix

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

Lecture 4 and 5 Controllability and Observability: Kalman decompositions

Lecture 4 and 5 Controllability and Observability: Kalman decompositions 1 Lecture 4 and 5 Controllability and Observability: Kalman decompositions Spring 2013 - EE 194, Advanced Control (Prof. Khan) January 30 (Wed.) and Feb. 04 (Mon.), 2013 I. OBSERVABILITY OF DT LTI SYSTEMS

More information

4.0 Update Algorithms For Linear Closed-Loop Systems

4.0 Update Algorithms For Linear Closed-Loop Systems 4. Update Algorithms For Linear Closed-Loop Systems A controller design methodology has been developed that combines an adaptive finite impulse response (FIR) filter with feedback. FIR filters are used

More information

Lecture 7 LQG Design. Linear Quadratic Gaussian regulator Control-estimation duality SRL for optimal estimator Example of LQG design for MIMO plant

Lecture 7 LQG Design. Linear Quadratic Gaussian regulator Control-estimation duality SRL for optimal estimator Example of LQG design for MIMO plant L7: Lecture 7 LQG Design Linear Quadratic Gaussian regulator Control-estimation duality SRL for optimal estimator Example of LQG design for IO plant LQG regulator L7:2 If the process and measurement noises

More information

Jordan blocks. Defn. Let λ F, n Z +. The size n Jordan block with e-value λ is the n n upper triangular matrix. J n (λ) =

Jordan blocks. Defn. Let λ F, n Z +. The size n Jordan block with e-value λ is the n n upper triangular matrix. J n (λ) = Jordan blocks Aim lecture: Even over F = C, endomorphisms cannot always be represented by a diagonal matrix. We give Jordan s answer, to what is the best form of the representing matrix. Defn Let λ F,

More information

Cross Directional Control

Cross Directional Control Cross Directional Control Graham C. Goodwin Day 4: Lecture 4 16th September 2004 International Summer School Grenoble, France 1. Introduction In this lecture we describe a practical application of receding

More information

AUTOMATIC CONTROL COMMUNICATION SYSTEMS LINKÖPING

AUTOMATIC CONTROL COMMUNICATION SYSTEMS LINKÖPING "!# $ %'&)(+* &-,.% /03254-687:9@?A?AB54 C DFEHG)IJ237KI#L BM>A>@ION B5P Q ER0EH?@EHBM4.B3PTSU;V68BMWX2368ERY@BMI Q 7K[25>@6AWX7\4)6]B3PT^_IH7\Y\6A>AEHYK25I#^_4`MER47K7\>AER4` a EH4GbN

More information

Analysis and Synthesis of Single-Input Single-Output Control Systems

Analysis and Synthesis of Single-Input Single-Output Control Systems Lino Guzzella Analysis and Synthesis of Single-Input Single-Output Control Systems l+kja» \Uja>)W2(ja»\ um Contents 1 Definitions and Problem Formulations 1 1.1 Introduction 1 1.2 Definitions 1 1.2.1 Systems

More information

Let the plant and controller be described as:-

Let the plant and controller be described as:- Summary of Fundamental Limitations in Feedback Design (LTI SISO Systems) From Chapter 6 of A FIRST GRADUATE COURSE IN FEEDBACK CONTROL By J. S. Freudenberg (Winter 2008) Prepared by: Hammad Munawar (Institute

More information

Control System Design

Control System Design ELEC4410 Control System Design Lecture 19: Feedback from Estimated States and Discrete-Time Control Design Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science

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

Nonlinear Control Lecture 9: Feedback Linearization

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

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

Andrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain)

Andrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain) 1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Linear systems (frequency domain) 2 Motivations Consider an LTI system Thanks to the Lagrange s formula we can compute the motion of

More information

Efficient Equalization for Wireless Communications in Hostile Environments

Efficient Equalization for Wireless Communications in Hostile Environments Efficient Equalization for Wireless Communications in Hostile Environments Thomas Strohmer Department of Mathematics University of California, Davis, USA strohmer@math.ucdavis.edu http://math.ucdavis.edu/

More information

Notes on Determinants and Matrix Inverse

Notes on Determinants and Matrix Inverse Notes on Determinants and Matrix Inverse University of British Columbia, Vancouver Yue-Xian Li March 17, 2015 1 1 Definition of determinant Determinant is a scalar that measures the magnitude or size of

More information

AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Introduction to Automatic Control & Linear systems (time domain)

AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Introduction to Automatic Control & Linear systems (time domain) 1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Introduction to Automatic Control & Linear systems (time domain) 2 What is automatic control? From Wikipedia Control theory is an interdisciplinary

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

A Comparison of Fixed Point Designs And Time-varying Observers for Scheduling Repetitive Controllers

A Comparison of Fixed Point Designs And Time-varying Observers for Scheduling Repetitive Controllers Proceedings of the 46th IEEE Conference on Decision and Control New Orleans, LA, USA, Dec. 12-14, 27 A Comparison of Fixed Point Designs And Time-varying Observers for Scheduling Repetitive Controllers

More information

1 Linear Algebra Problems

1 Linear Algebra Problems Linear Algebra Problems. Let A be the conjugate transpose of the complex matrix A; i.e., A = A t : A is said to be Hermitian if A = A; real symmetric if A is real and A t = A; skew-hermitian if A = A and

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

Properties of Open-Loop Controllers

Properties of Open-Loop Controllers Properties of Open-Loop Controllers Sven Laur University of Tarty 1 Basics of Open-Loop Controller Design Two most common tasks in controller design is regulation and signal tracking. Regulating controllers

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

Every real system has uncertainties, which include system parametric uncertainties, unmodeled dynamics

Every real system has uncertainties, which include system parametric uncertainties, unmodeled dynamics Sensitivity Analysis of Disturbance Accommodating Control with Kalman Filter Estimation Jemin George and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY, 14-44 The design

More information

only nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr

only nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr The discrete algebraic Riccati equation and linear matrix inequality nton. Stoorvogel y Department of Mathematics and Computing Science Eindhoven Univ. of Technology P.O. ox 53, 56 M Eindhoven The Netherlands

More information

Finite-Time Behavior of Inner Systems

Finite-Time Behavior of Inner Systems 1134 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 7, JULY 2003 Finite-Time Behavior of Inner Systems Jobert H. A. Ludlage, Member, IEEE, Siep Weiland, Anton A. Stoorvogel, Senior Member, IEEE,

More information

João P. Hespanha. January 16, 2009

João P. Hespanha. January 16, 2009 LINEAR SYSTEMS THEORY João P. Hespanha January 16, 2009 Disclaimer: This is a draft and probably contains a few typos. Comments and information about typos are welcome. Please contact the author at hespanha@ece.ucsb.edu.

More information

Remark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.

Remark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero. Sec 5 Eigenvectors and Eigenvalues In this chapter, vector means column vector Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called

More information

Learn2Control Laboratory

Learn2Control Laboratory Learn2Control Laboratory Version 3.2 Summer Term 2014 1 This Script is for use in the scope of the Process Control lab. It is in no way claimed to be in any scientific way complete or unique. Errors should

More information

The Essentials of Linear State-Space Systems

The Essentials of Linear State-Space Systems :or-' The Essentials of Linear State-Space Systems J. Dwight Aplevich GIFT OF THE ASIA FOUNDATION NOT FOR RE-SALE John Wiley & Sons, Inc New York Chichester Weinheim OAI HOC OUOC GIA HA N^l TRUNGTAMTHANCTINTHUVIIN

More information

Algebra C Numerical Linear Algebra Sample Exam Problems

Algebra C Numerical Linear Algebra Sample Exam Problems Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric

More information

Background Mathematics (2/2) 1. David Barber

Background Mathematics (2/2) 1. David Barber Background Mathematics (2/2) 1 David Barber University College London Modified by Samson Cheung (sccheung@ieee.org) 1 These slides accompany the book Bayesian Reasoning and Machine Learning. The book and

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

A Comparative Study on Automatic Flight Control for small UAV

A Comparative Study on Automatic Flight Control for small UAV Proceedings of the 5 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'18) Niagara Falls, Canada June 7 9, 18 Paper No. 13 DOI: 1.11159/cdsr18.13 A Comparative Study on Automatic

More information

Problem Set 4 Solution 1

Problem Set 4 Solution 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 4 Solution Problem 4. For the SISO feedback

More information

Introduction to Signals and Systems Lecture #4 - Input-output Representation of LTI Systems Guillaume Drion Academic year

Introduction to Signals and Systems Lecture #4 - Input-output Representation of LTI Systems Guillaume Drion Academic year Introduction to Signals and Systems Lecture #4 - Input-output Representation of LTI Systems Guillaume Drion Academic year 2017-2018 1 Outline Systems modeling: input/output approach of LTI systems. Convolution

More information

Structured State Space Realizations for SLS Distributed Controllers

Structured State Space Realizations for SLS Distributed Controllers Structured State Space Realizations for SLS Distributed Controllers James Anderson and Nikolai Matni Abstract In recent work the system level synthesis (SLS) paradigm has been shown to provide a truly

More information

Automatic Control Systems theory overview (discrete time systems)

Automatic Control Systems theory overview (discrete time systems) Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations

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

Digital Control Systems State Feedback Control

Digital Control Systems State Feedback Control Digital Control Systems State Feedback Control Illustrating the Effects of Closed-Loop Eigenvalue Location and Control Saturation for a Stable Open-Loop System Continuous-Time System Gs () Y() s 1 = =

More information

Geometric Control Theory

Geometric Control Theory 1 Geometric Control Theory Lecture notes by Xiaoming Hu and Anders Lindquist in collaboration with Jorge Mari and Janne Sand 2012 Optimization and Systems Theory Royal institute of technology SE-100 44

More information

Jerk derivative feedforward control for motion systems

Jerk derivative feedforward control for motion systems Jerk derivative feedforward control for motion systems Matthijs Boerlage Rob Tousain Maarten Steinbuch Abstract This work discusses reference trajectory relevant model based feedforward design. For motion

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

Iterative learning control (ILC) is based on the notion

Iterative learning control (ILC) is based on the notion Iterative learning control (ILC) is based on the notion that the performance of a system that executes the same task multiple times can be improved by learning from previous executions (trials, iterations,

More information

Chapter 15 - Solved Problems

Chapter 15 - Solved Problems Chapter 5 - Solved Problems Solved Problem 5.. Contributed by - Alvaro Liendo, Universidad Tecnica Federico Santa Maria, Consider a plant having a nominal model given by G o (s) = s + 2 The aim of the

More information

Optimal Polynomial Control for Discrete-Time Systems

Optimal Polynomial Control for Discrete-Time Systems 1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should

More information

Discrete and continuous dynamic systems

Discrete and continuous dynamic systems Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty

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

10/8/2015. Control Design. Pole-placement by state-space methods. Process to be controlled. State controller

10/8/2015. Control Design. Pole-placement by state-space methods. Process to be controlled. State controller Pole-placement by state-space methods Control Design To be considered in controller design * Compensate the effect of load disturbances * Reduce the effect of measurement noise * Setpoint following (target

More information

Classifying and building DLMs

Classifying and building DLMs Chapter 3 Classifying and building DLMs In this chapter we will consider a special form of DLM for which the transition matrix F k and the observation matrix H k are constant in time, ie, F k = F and H

More information

Robotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007

Robotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007 Robotics & Automation Lecture 25 Dynamics of Constrained Systems, Dynamic Control John T. Wen April 26, 2007 Last Time Order N Forward Dynamics (3-sweep algorithm) Factorization perspective: causal-anticausal

More information

Lecture 13: Internal Model Principle and Repetitive Control

Lecture 13: Internal Model Principle and Repetitive Control ME 233, UC Berkeley, Spring 2014 Xu Chen Lecture 13: Internal Model Principle and Repetitive Control Big picture review of integral control in PID design example: 0 Es) C s) Ds) + + P s) Y s) where P s)

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

Lecture 2: Discrete-time Linear Quadratic Optimal Control

Lecture 2: Discrete-time Linear Quadratic Optimal Control ME 33, U Berkeley, Spring 04 Xu hen Lecture : Discrete-time Linear Quadratic Optimal ontrol Big picture Example onvergence of finite-time LQ solutions Big picture previously: dynamic programming and finite-horizon

More information

EE Control Systems LECTURE 9

EE Control Systems LECTURE 9 Updated: Sunday, February, 999 EE - Control Systems LECTURE 9 Copyright FL Lewis 998 All rights reserved STABILITY OF LINEAR SYSTEMS We discuss the stability of input/output systems and of state-space

More information

Computational Methods. Eigenvalues and Singular Values

Computational Methods. Eigenvalues and Singular Values Computational Methods Eigenvalues and Singular Values Manfred Huber 2010 1 Eigenvalues and Singular Values Eigenvalues and singular values describe important aspects of transformations and of data relations

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

INVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN. Leonid Lyubchyk

INVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN. Leonid Lyubchyk CINVESTAV Department of Automatic Control November 3, 20 INVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN Leonid Lyubchyk National Technical University of Ukraine Kharkov

More information

CDS Solutions to Final Exam

CDS Solutions to Final Exam CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space

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