Problem Set 5 Solutions 1

Size: px
Start display at page:

Download "Problem Set 5 Solutions 1"

Transcription

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 norm approximations and Q-parameterization, and their applications in feedback design. Problem 5.T Let G(s) = s + + s + a, where a > 0 is a real parameter. (a) Find an analytical expression for the Hankel norm of G. Answer: Reasoning: a state space model ( 4 + 4a ) 2 ( ) 4a 4a. (a + ) 2 ẋ = Ax + Bw, y = Cx + Dw for G is given by A = [ 0 0 a Version of December, 20. ] [, B = ], C = [ ], D = 0.

2 2 The corresponding controllability and observability Gramians are equal: ] W c = W o =. [ 2 a+ Hence Hankel norm of G, as the square root of the largest eigenvalue of W o W c, is the larger root of the characteristic polynomial ( δ(t) = t ) ( ) t + 4a 4a. (a + ) 2 This can be checked by MATLAB code like this: a+ 2a function ps5(a) s=tf( s ); sig=hsvd(/(s+)+/(s+a)); p=/4+/(4*a);q=/(4*a)-/(a+)^2;d=sqrt(p^2-q); siga=[p+d;p-d]; fprintf( ps5 error for a=%f: %e [%f %f]\n,... a,norm(sig-siga),siga(),siga(2)) (b) Find an analytical expression for minimal Hankel norm of the difference G Ĝ where Ĝ(s) = b/(s + c) for some b R and c > 0. What happens when a? Answer: the minimal equals ( 4 + 4a 4 + ) 2 4a ( ) 4a. (a + ) 2 Reasoning: since Ĝ(s) ranges over all stable transfer functions of order, the minimum of the Hankel norm of G Ĝ equals the second largest Hankel singular value of G, which is computed in (a). Since G has order when a = the minimum converges to zero as a. Problem 5.2T Consider the setup shown on Figure, where P is a given DT LTI state space model defined by equations x(t + ) = Ax(t) + Bu(t), q(t) = Cx(t),

3 3 and K is stabilizing DT LTI feedback, to be designed as a state space model x f (t + ) = A f x f (t) + B f y(t), u(t) = C f x f (t) + D f y(t). The signals u, w, q are assumed to be scalar. Let F, L be the real vectors (F is a row vector, L is a column vector) such that A+BF and A+LC are Schur matrices (i.e. all eigenvalues λ satisfy λ < ). Let G 0, G, G 2 be the corresponding transfer functions from the standard Q parameterization G = G 0 +G QG 2 of all achievable stabilized closed loop transfer functions from w to u. w + q K P u Figure : Feedback Setup for Problem 5.2T (a) Give state space models for G 0, G, G 2 in terms of A, B, C, F, L. Answer: G 0 : x + = (A + BF )x BF, + = (A + LC) + Lw, e = F x F, G : η + = (A + BF )η + Bv, ξ = F η + ξ, G 2 : + = (A + LC) + Lw, δ = C + w. Reasoning: apply the general Q-parametrization formulae with A = A, B = 0, B 2 = B, C = 0, D = 0, D 2 = I, C 2 = C, D 2 = I. (b) Let P (z) = P 0 (z z )(z z 2 )... (z z k ) (z p )(z p 2 )... (z p n ).

4 4 Assume for simplicity that z i z l and p i p l for i l. Express zeros of G and G 2 in terms of P 0, z i, p l. Answer: a number z C such that z is a zero of G (or G 2 ) if and only if z = p i for some i. Locations of stable zeros of G (and G 2 ) cannot be determined based on the assumptions. Reasoning: for z C such that the matrix zi A BF is invertible (in particular, for z ) we have G (z) = + F (zi A BF ) B = det( + F (zi A BF ) B) = det(i + (zi A BF ) BF ) = det((zi A BF ) ) det(zi A BF + BF ) = det(zi A) det(zi A BF ), G 2 (z) = + C(zI A LC) L =... det(zi A) = det(zi A LC). Since all uncontrollable or unobservable modes of the state space model for P must be stable, this establishes locations of the unstable zeros of G and G 2. Since uncontrollable or unobservable modes of a state space model do not show up in its transfer function, locations of stable zeros of G and G 2 cannot be determined based on the information provided. (c) Find lim z pi G 0 (z)/(z p i ) in terms of lim z pi (z p i )P (z). Answer: for p i we have lim G 0 (z)/(z p i ) = z p i lim z pi (z p i )P (z). For p i <, the limit cannot be determined based on the information provided. Reasoning: let K = K(s) be the transfer function of the stabilizing controller which yields closed loop transfer function G = G 0. Then z p K(z) K(z)P (z) = K(z)P (z) K(z)P (z) (z p)p (z).

5 5 When p = p i and p i, we have K(z)P (z) as z p, and therefore KP/( KP ), which yields the expressing in the answer. When p = p i and p i <, it is admissible for K to have a zero at z = p i, and the limit value of KP/( KP ) at p can be arbitrary. What does this all tell about the set of all possible closed loop transfer functions G from w to u? A transfer function G is achievable in a stabilized closed loop if and only if it is stable and, for every p i with p i, satisfies the interpolation constraints G(p i ) = 0, Ġ(p i ) = d dz P (z). z=pi Problem 5.3P In general, it is impossible to compute singular numbers of the Hankel operator associated with a given transfer function G H by just looking at the values of G(jω). However, the situation is different when G is all-pass, i.e. G(jω) = for all ω R. Let G be an all-pass CT stable transfer function of order n, or, equivalently, G(s) = s a s + a s a 2 s + a 2... s a n s + a n, where a i > 0 are some positive constants. Assume for simplicity (though the assumption is not really necessary) that a i a k for i k. Let H G denote the Hankel operator associated with G. Since G has order n, the rank of H G equals n as well. (a) Run some numerical experiments with randomly generated a k to formulate a hypothesis on what are the singular numbers of H G. Hypotheses: σ = σ 2 = = σ n =. Sample code (in MATLAB): function ps53(n) if nargin<, n=5; end q=poly(-rand(n,)); p=((-).^(0:n)).*q; h=hsvd(tf(p,q)); close(gcf);bar(h);grid

6 6 (b) Is it possible for a singular number of H G to be greater than? Explain. Answer: no. Reasoning: the linear transformation L G : L 2 (R) L 2 (R) is defined by L G w = w+g w, where denotes convolution, and g = g(t) be the inverse Fourier transform of G. Naturally, if y = L G w and Y, W are Fourier transforms of w, y respectively then Y = GW. Since G(jω) for all ω, we have Y = W, hence, by the Parceval identity, w(t) 2 dt = 2π W (jω) 2 dω = Y (jω) 2 dω = w(t) 2 dt. 2π Since y h = H G w is defined by zeroing the values of y = L G w at all t < 0, L2 norm of H G w is not larger than that of w. Hence the largest singular number of H G is not larger than. (c) Find some non-zero square integrable functions w(t) which are anticausal (i.e. w(t) = 0 for t 0), and Y (jω) = G(jω)W (jω), where y = H G w, and Y and W are the Fourier transforms of y and w respectively. (Hint: use numbers a i as a guidance.) Answer: w(t) = where c i are arbitrary real coefficients. { n i= c ie a it, t < 0, 0, t 0, Reasoning: Fourier transform W = W (s) of w = w(t) is given by W (s) = n i= c i a i s. Since GW is a strictly proper rational function with no unstable poles, it is an element of H 2, hence L G w = H G w. (d) Use the observation made in (c) to prove the hypothesis formulated in (a). Proof: since H G w w for every w, and H G w = w over a linear subspace U = {w} of dimension n, the first n singular numbers of H G must be equal to.

7 7 So, does it make much sense to apply model reduction to all-pass transfer functions? Answer: no, or, at least, not with the objective to minimize L2 gain of the error. Reasoning: the best approximation of G by a transfer function Ĝ of order less than n will have L2 gain error G Ĝ of at least σ n (G) =, which is readily achieved by zero approximation Ĝ 0. K(s) e τs P 0 (s) Figure 2: Feedback Setup for Problem 5.4P Problem 5.4P This task calls for the use of approximation and model reduction in designing linear feebdack. To keep things simple, the attention is limited to the objective of feedback stabilization. Consider the feedback setup shown on Figure 2, where P 0 (s) = s s 2 (s 2) is a given plant, τ [0, T ] models the unknown feedback loop delay, T > 0 is a known upper bound for τ, and K is the finite order controller to be designed, to guarantee robust stabilization, understood as finite L2 gain feedback loop stability for all τ [0, T ]. (a) Compute analytically (as a function of τ 0) L2 gain of the LTI system defined by transfer function M τ (s) = e τs. s It is advisable to calculate M τ numerically first (e.g. by sampling the values of M τ on the imaginary axis), to get an idea of what the answer should be.

8 8 Answer: M τ = τ. Reasoning: the inverse Fourier transform h τ of M τ is given by {, 0θ < τ, h τ (t) = 0, otherwise. Hence M τ (jω) = τ which proves that M τ = τ. (b) Use the decomposition 0 τ e jωt h τ (t)dt h τ (t) dt = τ = M τ (0), e τs P 0 (s) = 0 ( ) s s + M τ(s) s(s 2), together with the result from (a) and the small gain theorem, to claim that a feedback K which is stabilizing for T = 0 achieves robust stabilization for some T > 0. Use H-Infinity optimization with hinfsyn.m to design a third order K which achieves robust stabilization with T as large as possible 2. Conclusion: stabilization guaranteed for T 0.. Figure 3: Design setup for Problem 5.4 Approach: we use design setup from Figure 3, and test setup from Figure 3, both controlled by ps54.m, a slightly modified version of which is shown below: 2 This does not mean finding the absolutely largest T, which is quite difficiult. Just try to do the best possible job using the approach.

9 9 function K=ps54(T) d=e-3; s=tf( s ); P=(s-)/(s*(s-2)); L=/s; r=t; assignin( base, d,d) % export variables assignin( base, P,P) assignin( base, L,L) p=linmod( ps54des ); % extract open loop p=ss(p.a,p.b,p.c,p.d); K=hinfsyn(p,,); % optimize controller assignin( base, K,K) % export controller p=linmod( ps54test ); % extract closed loop S=ss(p.a,p.b,p.c,p.d); fprintf( [T=%f] stability: %f<0, small gain: %f<%f\n,... T,max(real(eig(p.a))),norm(S,Inf),/r) Figure 4: Test setup for Problem 5.4 (c) The result in (b) was achieved by using the approximation M τ (s) 0, which has L2 gain error M τ. When τ > 0 is known, better designs can be made by finding a higher quality approximation M τ ˆM τ (s), where ˆM τ is a stable rational transfer function, and the error is quantified by M τ ˆM τ. A rather straightforward way of doing this is to start with the rough (inaccurate at high frequencies) approximation ( ) n τs/2n e τs Êτ(s) =, + τs/2n

10 0 where n is large enough, and to have ˆM τ defined as the reduced model ˆM τ (s) Êτ(s), s of order significantly smaller than n. Apply this approach (with hankelmr.m used for model reduction) to find a good third order approximation ˆM τ of M τ for a given τ > 0. Implementation: the following code (a part of ps54.m) shows the use of hankelmr.m: s=tf( s ); [A,B,C,D]=ssdata(ss((-s*(0.5*T/n))/(+s*(0.5*T/n)))^n); L0=hankelmr(ss(A,B,C/A,0),m); Here L0 is a Hankel optimal reduced model of (Êτ )/s. To avoid using minreal.m on a high order model, a state space of Êτ is formed first, and then we use the fact that if A is an invertible matrix such that ẋ = Ax + Bw, is a state space model for G(s) then ẋ = Ax + Bw, y = Cx + Dw v = CA (Ax + Bw) is a state space model for (G(s) G(0))/s. Since, apparently, hankelmr.m does not have an option for optimizing the D matrix of the reduced model, an alternative code using the old hankmr.m can do better: sys=pck(a,b,c/a,0); [sys,sig]=sysbal(sys); [A,B,C,D]=unpck(hankmr(sys,sig,m, d )); L0=ss(A,B,C,D); Though in general optimization of D can lower model reduction error significantly, in this particular application it only yields a minor improvement. (d) Use the results from (c) to design K of order 6 to achieve feebdack stabilization with a fixed τ which is as large as possible. Conclusion: stabilization is guaranteed for τ 0.36 (with D optimization using hankmr.m) or τ 0.32 (using hankelmr.m). The corresponding code is contained in the full version of ps54.m.

Problem Set 4 Solutions 1

Problem Set 4 Solutions 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6245: MULTIVARIABLE CONTROL SYSTEMS by A Megretski Problem Set 4 Solutions Problem 4T For the following transfer

More information

Problem Set 8 Solutions 1

Problem Set 8 Solutions 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6245: MULTIVARIABLE CONTROL SYSTEMS by A Megretski Problem Set 8 Solutions 1 The problem set deals with H-Infinity

More information

Hankel Optimal Model Reduction 1

Hankel Optimal Model Reduction 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Hankel Optimal Model Reduction 1 This lecture covers both the theory and

More information

Problem Set 2 Solutions 1

Problem Set 2 Solutions 1 Massachusetts Institute of echnology Department of Electrical Engineering and Computer Science 6.45: MULIVARIABLE CONROL SYSEMS by A. Megretski Problem Set Solutions Problem. For each of the statements

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

Chapter 5. Standard LTI Feedback Optimization Setup. 5.1 The Canonical Setup

Chapter 5. Standard LTI Feedback Optimization Setup. 5.1 The Canonical Setup Chapter 5 Standard LTI Feedback Optimization Setup Efficient LTI feedback optimization algorithms comprise a major component of modern feedback design approach: application problems involving complex models

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

Problem set 5 solutions 1

Problem set 5 solutions 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION Problem set 5 solutions Problem 5. For each of the stetements below, state

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

10 Transfer Matrix Models

10 Transfer Matrix Models MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important

More information

6.241 Dynamic Systems and Control

6.241 Dynamic Systems and Control 6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May

More information

6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski. Solutions to Problem Set 1 1. Massachusetts Institute of Technology

6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski. Solutions to Problem Set 1 1. Massachusetts Institute of Technology Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Solutions to Problem Set 1 1 Problem 1.1T Consider the

More information

L2 gains and system approximation quality 1

L2 gains and system approximation quality 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION L2 gains and system approximation quality 1 This lecture discusses the utility

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

Introduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31

Introduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31 Introduction Classical Control Robust Control u(t) y(t) G u(t) G + y(t) G : nominal model G = G + : plant uncertainty Uncertainty sources : Structured : parametric uncertainty, multimodel uncertainty Unstructured

More information

Balanced Truncation 1

Balanced Truncation 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [

More information

Module 07 Controllability and Controller Design of Dynamical LTI Systems

Module 07 Controllability and Controller Design of Dynamical LTI Systems Module 07 Controllability and Controller Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October

More information

Lecture 15: H Control Synthesis

Lecture 15: H Control Synthesis c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control for Linear Systems: Lecture 15 p. 1/14 Lecture 15: H Control Synthesis Example c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control

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

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 3. Chapter 4: Elements of Linear System Theory. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Lecture 3. Chapter 4: Elements of Linear System Theory. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University. Lecture 3 Chapter 4: Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 3 p. 1/77 3.1 System Descriptions [4.1] Let f(u) be a liner operator, u 1 and u

More information

Multivariable Control. Lecture 03. Description of Linear Time Invariant Systems. John T. Wen. September 7, 2006

Multivariable Control. Lecture 03. Description of Linear Time Invariant Systems. John T. Wen. September 7, 2006 Multivariable Control Lecture 3 Description of Linear Time Invariant Systems John T. Wen September 7, 26 Outline Mathematical description of LTI Systems Ref: 3.1-3.4 of text September 7, 26Copyrighted

More information

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems

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

More information

to have roots with negative real parts, the necessary and sufficient conditions are that:

to have roots with negative real parts, the necessary and sufficient conditions are that: THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 543 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 7 Sebastian A. Nugroho November 6, 7 Due date of the homework is: Sunday, November 6th @ :59pm.. The following

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

Some solutions of the written exam of January 27th, 2014

Some solutions of the written exam of January 27th, 2014 TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following

More information

Robust Multivariable Control

Robust Multivariable Control Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms

More information

Linear System Theory

Linear System Theory Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability

More information

Control Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017

More information

An Internal Stability Example

An Internal Stability Example An Internal Stability Example Roy Smith 26 April 2015 To illustrate the concept of internal stability we will look at an example where there are several pole-zero cancellations between the controller and

More information

FEL3210 Multivariable Feedback Control

FEL3210 Multivariable Feedback Control FEL3210 Multivariable Feedback Control Lecture 5: Uncertainty and Robustness in SISO Systems [Ch.7-(8)] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 5:Uncertainty and Robustness () FEL3210 MIMO

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

1. Find the solution of the following uncontrolled linear system. 2 α 1 1 Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +

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

FEL3210 Multivariable Feedback Control

FEL3210 Multivariable Feedback Control FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction

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

sc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11

sc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11 sc46 - Control Systems Design Q Sem Ac Yr / Mock Exam originally given November 5 9 Notes: Please be reminded that only an A4 paper with formulas may be used during the exam no other material is to be

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

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

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

Modern Control Systems

Modern Control Systems Modern Control Systems Matthew M. Peet Illinois Institute of Technology Lecture 18: Linear Causal Time-Invariant Operators Operators L 2 and ˆL 2 space Because L 2 (, ) and ˆL 2 are isomorphic, so are

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

3 Gramians and Balanced Realizations

3 Gramians and Balanced Realizations 3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations

More information

Singular Value Decomposition Analysis

Singular Value Decomposition Analysis Singular Value Decomposition Analysis Singular Value Decomposition Analysis Introduction Introduce a linear algebra tool: singular values of a matrix Motivation Why do we need singular values in MIMO control

More information

Linear Quadratic Gausssian Control Design with Loop Transfer Recovery

Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Leonid Freidovich Department of Mathematics Michigan State University MI 48824, USA e-mail:leonid@math.msu.edu http://www.math.msu.edu/

More information

Topic # Feedback Control Systems

Topic # Feedback Control Systems Topic #20 16.31 Feedback Control Systems Closed-loop system analysis Bounded Gain Theorem Robust Stability Fall 2007 16.31 20 1 SISO Performance Objectives Basic setup: d i d o r u y G c (s) G(s) n control

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

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

The Generalized Laplace Transform: Applications to Adaptive Control*

The Generalized Laplace Transform: Applications to Adaptive Control* The Transform: Applications to Adaptive * J.M. Davis 1, I.A. Gravagne 2, B.J. Jackson 1, R.J. Marks II 2, A.A. Ramos 1 1 Department of Mathematics 2 Department of Electrical Engineering Baylor University

More information

Module 08 Observability and State Estimator Design of Dynamical LTI Systems

Module 08 Observability and State Estimator Design of Dynamical LTI Systems Module 08 Observability and State Estimator Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha November

More information

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition

More information

CME 345: MODEL REDUCTION

CME 345: MODEL REDUCTION CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of

More information

H 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University

H 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University H 2 Optimal State Feedback Control Synthesis Raktim Bhattacharya Aerospace Engineering, Texas A&M University Motivation Motivation w(t) u(t) G K y(t) z(t) w(t) are exogenous signals reference, process

More information

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016 ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe

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

Topic # Feedback Control Systems

Topic # Feedback Control Systems Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the

More information

Denis ARZELIER arzelier

Denis ARZELIER   arzelier COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15

More information

Definition of the Laplace transform. 0 x(t)e st dt

Definition of the Laplace transform. 0 x(t)e st dt Definition of the Laplace transform Bilateral Laplace Transform: X(s) = x(t)e st dt Unilateral (or one-sided) Laplace Transform: X(s) = 0 x(t)e st dt ECE352 1 Definition of the Laplace transform (cont.)

More information

Model Reduction for Unstable Systems

Model Reduction for Unstable Systems Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model

More information

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods

More information

Control Systems Design

Control Systems Design ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline

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

6.241 Dynamic Systems and Control

6.241 Dynamic Systems and Control 6.241 Dynamic Systems and Control Lecture 12: I/O Stability Readings: DDV, Chapters 15, 16 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology March 14, 2011 E. Frazzoli

More information

Controllability, Observability, Full State Feedback, Observer Based Control

Controllability, Observability, Full State Feedback, Observer Based Control Multivariable Control Lecture 4 Controllability, Observability, Full State Feedback, Observer Based Control John T. Wen September 13, 24 Ref: 3.2-3.4 of Text Controllability ẋ = Ax + Bu; x() = x. At time

More information

ECEEN 5448 Fall 2011 Homework #5 Solutions

ECEEN 5448 Fall 2011 Homework #5 Solutions ECEEN 5448 Fall 211 Homework #5 Solutions Professor David G. Meyer December 8, 211 1. Consider the 1-dimensional time-varying linear system ẋ t (u x) (a) Find the state-transition matrix, Φ(t, τ). Here

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

6.302 Feedback Systems Recitation 16: Compensation Prof. Joel L. Dawson

6.302 Feedback Systems Recitation 16: Compensation Prof. Joel L. Dawson Bode Obstacle Course is one technique for doing compensation, or designing a feedback system to make the closed-loop behavior what we want it to be. To review: - G c (s) G(s) H(s) you are here! plant For

More information

Feedback Stabilization over Signal-to-Noise Ratio Constrained Channels

Feedback Stabilization over Signal-to-Noise Ratio Constrained Channels Feedback Stabilization over Signal-to-Noise Ratio Constrained Channels Julio H. Braslavsky, Rick H. Middleton and Jim S. Freudenberg Abstract There has recently been significant interest in feedback stabilization

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

Full-State Feedback Design for a Multi-Input System

Full-State Feedback Design for a Multi-Input System Full-State Feedback Design for a Multi-Input System A. Introduction The open-loop system is described by the following state space model. x(t) = Ax(t)+Bu(t), y(t) =Cx(t)+Du(t) () 4 8.5 A =, B =.5.5, C

More information

Laplace Transforms and use in Automatic Control

Laplace Transforms and use in Automatic Control Laplace Transforms and use in Automatic Control P.S. Gandhi Mechanical Engineering IIT Bombay Acknowledgements: P.Santosh Krishna, SYSCON Recap Fourier series Fourier transform: aperiodic Convolution integral

More information

Richiami di Controlli Automatici

Richiami di Controlli Automatici Richiami di Controlli Automatici Gianmaria De Tommasi 1 1 Università degli Studi di Napoli Federico II detommas@unina.it Ottobre 2012 Corsi AnsaldoBreda G. De Tommasi (UNINA) Richiami di Controlli Automatici

More information

Semidefinite Programming Duality and Linear Time-invariant Systems

Semidefinite Programming Duality and Linear Time-invariant Systems Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,

More information

Control Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho Tel: Fax:

Control Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho Tel: Fax: Control Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 02-760-4253 Fax:02-760-4435 Introduction In this lesson, you will learn the following : How to determine

More information

Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection

Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute

More information

Introduction & Laplace Transforms Lectures 1 & 2

Introduction & Laplace Transforms Lectures 1 & 2 Introduction & Lectures 1 & 2, Professor Department of Electrical and Computer Engineering Colorado State University Fall 2016 Control System Definition of a Control System Group of components that collectively

More information

The norms can also be characterized in terms of Riccati inequalities.

The norms can also be characterized in terms of Riccati inequalities. 9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements

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

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

Module 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control

Module 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control Module 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control Ahmad F. Taha EE 3413: Analysis and Desgin of Control Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/

More information

Solution of Linear State-space Systems

Solution of Linear State-space Systems Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state

More information

Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations

Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Klajdi Sinani Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University

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

EE 3054: Signals, Systems, and Transforms Summer It is observed of some continuous-time LTI system that the input signal.

EE 3054: Signals, Systems, and Transforms Summer It is observed of some continuous-time LTI system that the input signal. EE 34: Signals, Systems, and Transforms Summer 7 Test No notes, closed book. Show your work. Simplify your answers. 3. It is observed of some continuous-time LTI system that the input signal = 3 u(t) produces

More information

LMI Methods in Optimal and Robust Control

LMI Methods in Optimal and Robust Control LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 14: LMIs for Robust Control in the LF Framework ypes of Uncertainty In this Lecture, we will cover Unstructured,

More information

A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case

A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case ATARU KASE Osaka Institute of Technology Department of

More information

1 Controller Optimization according to the Modulus Optimum

1 Controller Optimization according to the Modulus Optimum Controller Optimization according to the Modulus Optimum w G K (s) F 0 (s) x The goal of applying a control loop usually is to get the control value x equal to the reference value w. x(t) w(t) X(s) W (s)

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 2 Laplace Transform I 1/52

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52 1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n

More information

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

CDS 101: Lecture 5-1 Reachability and State Space Feedback CDS 11: Lecture 5-1 Reachability and State Space Feedback Richard M. Murray 23 October 26 Goals: Define reachability of a control system Give tests for reachability of linear systems and apply to examples

More information

Massachusetts Institute of Technology

Massachusetts Institute of Technology Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.11: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 21 QUESTION BOOKLET

More information

Model reduction for linear systems by balancing

Model reduction for linear systems by balancing Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,

More information

EE221A Linear System Theory Final Exam

EE221A Linear System Theory Final Exam EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,

More information

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

PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT

PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT Hans Norlander Systems and Control, Department of Information Technology Uppsala University P O Box 337 SE 75105 UPPSALA, Sweden HansNorlander@ituuse

More information

Lecture 7 (Weeks 13-14)

Lecture 7 (Weeks 13-14) Lecture 7 (Weeks 13-14) Introduction to Multivariable Control (SP - Chapters 3 & 4) Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 7 (Weeks 13-14) p.

More information

Problem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations:

Problem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations: EEE58 Exam, Fall 6 AA Rodriguez Rules: Closed notes/books, No calculators permitted, open minds GWC 35, 965-37 Problem (Dynamic Augmentation: State Space Representation) Consider a dynamical system consisting

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

GATE EE Topic wise Questions SIGNALS & SYSTEMS

GATE EE Topic wise Questions SIGNALS & SYSTEMS www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)

More information