CDS Solutions to Final Exam
|
|
- Aubrie Tyler
- 5 years ago
- Views:
Transcription
1 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 representation of the system with transfer function G(s) Assuming u( ) = y( ) = as usual, we have: s 3 y(s) + 2s 2 y(s) + sy(s) + y(s) = u(s) y (t) + 2ÿ(t) + ẏ(t) + y(t) = u(t) Defining x = y, x 2 = ẏ, x 3 = ÿ, the corresponding state-space representation is given by: ẋ x ẋ 2 = x 2 + u ẋ 3 2 x 3 y = [ Since A is Hurwitz (one way of verifying this is using the Routh-Hurwitz criterion), the Lyapunov equation x x 2 x 3 AL + LA + BB = has a unique solution L = L, and G 2 = CLC What is left is to solve the Lyapunov equation Letting: a a a 2 L = a a 3 a 4 a 2 a 4 a 5 and substituting, we get the following system of equations: 2a = a 2 + a 3 = 2a 4 = a 4 a a 2a 2 = a 5 a a 3 2a 4 = 2[a 2 + a 4 + 2a 5 =
2 The solution to this system of equations is: a =, a = a 4 =, a 2 = 2, a 3 = a 5 = 2 By direct substitution, the H 2 norm of G can now be computed to be equal to (b) We begin by noting that the relation between inputs u and u 2 and outputs y and y 2 given by: [ [ [ y (s) /s /s u (s) = y 2 (s) /s u 2 (s) can be represented in terms of the block diagram in Figure u (s) /s x (s) + + y (s) u 2 (s) /s x 2 (s) y 2 (s) Figure : Block diagram for problem (b) Let x (s) and x 2 (s) be the signals at the output of the two integrators in Figure, we have: sx (s) = u (s) sx 2 (s) = u 2 (s) y (s) = x (s) + x 2 (s) y 2 (s) = x 2 (s) The corresponding state-space representation is given by: [ [ [ [ [ ẋ x u = + ẋ 2 x 2 u 2 [ [ [ y x = y 2 Since both B and C have full rank, (A,B) and (C,A) are reachable and observable, respectively, and the realization is indeed minimal What is left is to verify (as a sanity check) that: x 2 C(sI A) B = C(sI) I = H(s) which is readily done by direct substitution 2
3 (c) By definition, the poles of the transfer matrix are those of the individual entries Hence, we have a pole at s = 2 and two repeated poles at s = To verify that, note that the number of poles of the system is equal to the order of the minimal realization of the system, which in this case is 3 (One way of seeing that is by drawing a block diagram as in part (b) and noting that we need exactly 3 first order, proper transfer blocks of order to obtain the desired input-output map) The zeroes of the transfer matrix are those values of s (finite or infinite) at which H(s) loses rank By inspection, note that: [ H() = 2 2 has rank, hence s = is a zero of the system Also note that: [ lim H(s) = 2 s also has rank, hence s = is an (infinite) zero of the system Finally, note that for [ 2(s 2) u(s) = (s + )(s ) we have: and: [ lim u(s) = s 2 3 [ lim H(s)u(s) = s 2 Hence s = 2 is also a zero of the system Note that there cannot be any more zeroes since the number of zeroes cannot exceed the number of poles Remark: Alternatively, the poles and finite zeroes of H(s) could have been computed by finding a minimal state-space realization of the system, computing the eigenvalues of A to determine the poles of the system, and finding the values of s C at which the system matrix [ A si B C D loses rank to compute the finite zeroes of the system (d) In order to verify stability of the equilibrium point at (,), we begin by finding the set of equilibrium points of the system The solutions to the system of equations: { 2x x 2 2 = 2x 2 x 2 4x 3 2 = 3
4 are the set of points (α,), α R, ie the x axis Since every neighborhood of the origin contains an equilibrium point distinct from (, ), the origin cannot be asymptotically stable To verify stability isl, consider the candidate Lyapunov function V (x,x 2 ) = x 2 + x 2 2 Note that V (x), x R 2, and that V (x) = x = x 2 = Also note that V (x) = [ V x [ V ẋ x 2 ẋ 2 = 2x 2 2(3x 2 + 2x 2 2), x R 2 V(x) is thus indeed a Lyapunov function for the system, and the origin is stable isl Problem 2 (a) Omitted (b) FALSE The system has a pole at s = 2, and is hence unstable To see that, consider the input { t [, u(t) = otherwise Note that u 2 2 = u 2 (t)dt = Also note that y() = and that y(t) = e 2(t ) y(), for t Hence y 2 2 = e 2( T) u(t)dt = e 2 e 2T dt = 2 (e2 ) y 2 (t)dt y 2 (t)dt = ( e 2 y() ) 2 ( e 2t ) 2 dt = Thus, the system can infinitely amplify an input signal with bounded L 2 norm (c) TRUE Let M(s) = s s + 5 and let (s) = e st We begin by s rewriting H(s), by simple algebraic manipulation, as: H(s) = s + 5 M(s) (s) 4
5 Since transfer function M(s) (s) is itself stable is stable, H(s) is stable iff the transfer function s + 5 Note that M = sup M(jw) = sup w w w = By the small gain w theorem, the feedback interconnection of M and is guaranteed to be internally stable if < In particular, the transfer function from u to y in Figure is then guaranteed to be stable u + + (s) M(s) y Figure 2: Block diagram for problem 2(c) Noting that the transfer function from u to y as shown in Figure is T u y = M(s) (s) it is clear that < is a sufficient condition for H(s) to be stable What is thus left to show is that < for all t [,) We have: and = sup w e jwt w e jwt = e j w w = tw = t w where the inequality uses the hint Hence for all t [,), we have: e jwt w t < and which concludes the proof sup w e jwt w < 5
6 (d) TRUE Problem 3 It is possible to derive a weak lower bound for T from the given bound on S(jw) in the interval w [ 4,4 Indeed, the algebraic constraint on S and T: S + T = and the triangle inequality allow us to conclude: = S + T S + T T S S(jw) w 4 = 9 A slightly stronger lower bound on T can be obtained by considering the effect of right half plane poles and zeroes of P (which are also right half plane poles and zeroes of the open loop transfer function L = PC, since there can be no unstable pole/zero cancellations between P and C) In particular, note that the pole at s = 3 and the zero at s = 2 impose the following interpolation constraints on S and T: S(2) = T(2) = S(3) = T(3) = We can immediately conclude, noting that T is stable hence analytic on the closed right half plane, and using the Maximum Modulus Principle that: T = sup T(jw) T(3) = w The lower bound can be improved further Since S is a stable non-minimum phase (S(3) = ) transfer function, it can be written as the product of a stable minimum phase system S MP and a stable all-pass filter S AP of the form: S AP = s 3 s + 3 S AP where S AP is a stable all-pass filter to account for other (potential) non-minimum phase zeroes of S: S = s 3 s + 3 S APS MP We have: S(jw) = S AP (jw)s MP (jw) = S AP (jw) S MP (jw) = S MP (jw), w We also have: S(2) = 5 S AP(2)S MP (2) 6
7 Hence: = S(2) = 5 S AP(2) S MP (2) S MP (2) = 5 S AP (2) 5 Finally, noting that for any stable, rational transfer function S, log S is analytic and bounded in the closed right half plane, we can apply Poisson s Integral Formula We have: log S MP (2) = 2 log S MP (jw) π 2 2 dw log(5) + w2 2 2 log S(jw) π 2 2 dw log(5) + w2 4 π 4 log S(jw) w 2 dw + 4 π 4 log S(jw) 2 2 dw log(5) + w2 Letting α = sup S(jw), and recalling that S(jw) for w < 4, we have: w >4 Noting that: and: We get: 4 4 π log() w 2 dw + 4 π log(α) dw log(5) + w w 2 dw = 2 arctan w 4 = arctan (2) w 2 dw = 2 arctan w 2 4 = 2 arctan (2) + π 4 log(α) log(5) 2 π log()arctan(2) 2 π arctan(2) 39 Hence S α Once again using the constraint S + T = and the triangle inequality, we have: S T S + T = T S 229 which is a much stronger lower bound Problem 4 Since the open loop system is reachable, for any choice of terminal state x o, there exists an input sequence, say u, u,, u n taking the system from the origin to x o in n steps Let, x, x 2,,x n be the state trajectory of the open 7
8 loop system corresponding to this input sequence, and note that for the closed loop system under the given non-linear feedback, we have: By setting the exogenous input to be: x(k + ) = Ax(k) + B (w(k) + f(x(k))) w (k) = u k f(x k) the state evolution of the closed loop system (starting from the origin) under w is identical to that of the open loop system under u ; in particular, x(n) = x o Since the choice of x o was arbitrary, the system remains reachable under arbitrary non-linear state feedback Problem 5 We begin by noting that the relation between the initial state and the inputs and outputs of the system up until time step n is given by: Y = y() y() y(n ) = C CA CA n = θx() + MU x() + D CB D CA n 2 B CB D The initial state can be uniquely determined from U and Y iff the observability matrix θ has full column rank (ie is left invertible), in which case we have: where θ L is the left inverse of θ x() = θ L (Y MU) (a) When the system is observable, the initial state can be uniquely determined from the output measurements alone iff: θ L M = () The given conditions are sufficient to ensure that () holds since M = under these conditions (b) When the output y is a scalar, the observability matrix is a square matrix of size n In this case, for an observable system, the observability matrix has a left and a right inverse, θ L = θ R = θ, and condition () holds iff M =, from which it follows that the conditions in part (a) are both necessary and sufficient u() u(n ) 8
9 for the initial state to be uniquely determined from output measurements alone (c) The transfer function of the system is the z-transform of its impulse response, assuming zero initial conditions In this case, we have: y() y() D CB D = CAB CB D Hence G(z) = z n y(n) n= = D + z CB + CAB + z2 = D + z C[I + A z + A2 z 2 + B = D + z C(I A z ) B = D + C(zI A) B We wish to prove that in the SISO case, the conditions of part (a) are both necessary and sufficient conditions for G(z) to have no finite zeroes To prove necessity, note that if G(z) has no finite zeroes, then G(z) exists and we can write: G(z) y(z) = u(z) Thus the input corresponding to a given ouput signal is unique (for zero initial conditions), and the initial state can be uniquely determined from output measurements alone It follows from parts (a) and (b) that the conditions on the Markov parameters of the system stated in part (a) hold The proof of sufficiency is more involved We begin by making some basic observations First, consider the characteristic polynomial of matrix zi A: a (z) + a (z)λ + + λ n = Recall that by the Cayley-Hamilton Theorem, we have: Recalling the expansion: a (z)i + a (z)(zi A) + + (zi A) n = (2) (a b) i = i j= ( i j ) a i j ( b) j 9
10 we note that (zi A) i can be written as a sum of terms involving powers of A (and z) up to order i In particular, when i n 2, under the conditions in part (a), we have: C(zI A) i B = while: C(zI A) n B = ( ) n CA n B Now suppose that z o is a finite zero of G(z) and that the conditions of part (a) hold; that is G(z o ) = C(z o I A) B = (i) If a (z o ), by rearranging (2) and multiplying by C and B, then using the above observations, we have: G(z o ) = C(z o I A) B [ n 2 = a (z o ) C (z o I A) n a i (z o )(z o I A) i B = ( )n a (z o ) CAn B = Equivalently, CA n B =, θb =, and B = Hence, barring the trivial case where B = (in which case the system is unforced and its transfer function is identically ) the transfer function cannot have finite zeroes (ii) When a (z o ) =, it can be shown by a similar argument that z o is both a pole and zero of the transfer function and hence cancels out (hence it is not a zero of the transfer function, although it is one of the given non-minimal state-space realization) i= Problem 6 The plant [ can be represented as shown in Figure δ Let = The transfer function M(s) with inputs w δ, w 2 and outputs z, z 2 is given by: M(s) = [ CNδ N δ CD δ D δ Note that M is stable by design Moreover, note that M is a rank one matrix, and can be decomposed as follows: [ M(s) = Nδ D δ We are asked to compute: = ab [ C γ o = inf w γ(w)
11 z δ w N δ D w δ 2 δ z N o D o y K Figure 3: Block Diagram of problem 6 where We have: γ(w) = inf {max γ i : det(i M(jw) (jw)) = } i det(i M ) = det(i ab ) = det( b a) = b a = [ δ δ n Let c i = D δi CN δi We can rewrite γ(w) as: D δ CN δ D δn CN δ γ(w) = inf δ R n{max γ i : i n γ i c i = } i= This problem is substantially more difficult in this case than in the case where δ is complex To gain some intuition, consider the three dimensional case (n = 3, the axis is (δ,δ 2,δ 3 )); the generalization to higher n follows At each frequency w, define: P = {δ : P 2 = {δ : 3 δ i I(c i (jw)) = ) i= 3 δ i R(c i (jw)) = ) i=
12 L = P P 2 In general, L R 3 is a line This makes our problem equivalent to finding the smallest β such that the cube B β = {δ R n : max i δ i β} touches the line L That can be done by looking at the following projections: For every j such that I(c j ), do the following: Look for the smallest β j such that max i δ i = β j and δ i (R(c i ) α j I(c i )) = i where α j = R(cj) I(c For each such j, it can be shown that the minimum β j) j is Èi R(ci) αji(ci) and the optimal δ is δ i = β j sgn(r(c i ) α j I(c i )) By doing this, we have found the smallest β j such that the projection of B βj in δ j = touches the projection of L Among È all of the above candidate solutions, the only admissible is j i j such that δ j = I(ci)δi c j satisfies δ j β j Thus, at each frequency w, γ(w) is given by: γ(w) = i R(c i(jw)) α j I(c i (jw)) 2
CDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
More informationCDS Final Exam
CDS 22 - Final Exam Instructor: Danielle C. Tarraf December 4, 2007 INSTRUCTIONS : Please read carefully! () Description & duration of the exam: The exam consists of 6 problems. You have a total of 24
More informationECE504: 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 informationMASSACHUSETTS 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.24: Dynamic Systems Spring 20 Homework 9 Solutions Exercise 2. We can use additive perturbation model with
More informationDissipativity. Outline. Motivation. Dissipative Systems. M. Sami Fadali EBME Dept., UNR
Dissipativity M. Sami Fadali EBME Dept., UNR 1 Outline Differential storage functions. QSR Dissipativity. Algebraic conditions for dissipativity. Stability of dissipative systems. Feedback Interconnections
More informationCONTROL 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 informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationFinal: Signal, Systems and Control (BME )
Final: Signal, Systems and Control (BME 580.) Instructor: René Vidal May 0th 007 HONOR SYSTEM: This examination is strictly individual. You are not allowed to talk, discuss, exchange solutions, etc., with
More informationEEE582 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 informationL2 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 informationSome 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 informationControl 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 informationChapter 30 Minimality and Stability of Interconnected Systems 30.1 Introduction: Relating I/O and State-Space Properties We have already seen in Chapt
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter
More informationIntroduction. 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 informationProblem 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 informationBalanced 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 informationChapter 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 informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear
More informationKalman 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 informationMultivariable 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 informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More informationZeros 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 informationState will have dimension 5. One possible choice is given by y and its derivatives up to y (4)
A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationDesign 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 informationIntro. 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 informationTheorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).
Linear Systems Notes Lecture Proposition. A M n (R) is positive definite iff all nested minors are greater than or equal to zero. n Proof. ( ): Positive definite iff λ i >. Let det(a) = λj and H = {x D
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationLaplace Transform Analysis of Signals and Systems
Laplace Transform Analysis of Signals and Systems Transfer Functions Transfer functions of CT systems can be found from analysis of Differential Equations Block Diagrams Circuit Diagrams 5/10/04 M. J.
More informationModule 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 informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016 1. Please write your 1- or 2-digit exam number on
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More information4F3 - Predictive Control
4F3 Predictive Control - Discrete-time systems p. 1/30 4F3 - Predictive Control Discrete-time State Space Control Theory For reference only Jan Maciejowski jmm@eng.cam.ac.uk 4F3 Predictive Control - Discrete-time
More information6.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 informationLecture 4 Stabilization
Lecture 4 Stabilization This lecture follows Chapter 5 of Doyle-Francis-Tannenbaum, with proofs and Section 5.3 omitted 17013 IOC-UPC, Lecture 4, November 2nd 2005 p. 1/23 Stable plants (I) We assume that
More informationECE 602 Exam 2 Solutions, 3/23/2011.
NAME: ECE 62 Exam 2 Solutions, 3/23/211. You can use any books or paper notes you wish to bring. No electronic devices of any kind are allowed. You can only use materials that you bring yourself. You are
More informationECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008
ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28-Oct-2008 Worcester Polytechnic Institute D. Richard Brown III 28-Oct-2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are
More informationME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms
I/O Stability Preliminaries: Vector and function norms 1. Sup norms are used for vectors for simplicity: x = max i x i. Other norms are also okay 2. Induced matrix norms: let A R n n, (i stands for induced)
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationAutomatic 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 informationProperties 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 informationModule 4 : Laplace and Z Transform Problem Set 4
Module 4 : Laplace and Z Transform Problem Set 4 Problem 1 The input x(t) and output y(t) of a causal LTI system are related to the block diagram representation shown in the figure. (a) Determine a differential
More informationProjection Theorem 1
Projection Theorem 1 Cauchy-Schwarz Inequality Lemma. (Cauchy-Schwarz Inequality) For all x, y in an inner product space, [ xy, ] x y. Equality holds if and only if x y or y θ. Proof. If y θ, the inequality
More informationChapter Stability Robustness Introduction Last chapter showed how the Nyquist stability criterion provides conditions for the stability robustness of
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 Stability
More informationLinear Systems. Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output
Linear Systems Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output Our interest is in dynamic systems Dynamic system means a system with memory of course including
More informationMASSACHUSETTS 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 informationAnalysis 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 informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2014
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2014 1. Please write your 1- or 2-digit exam number on
More informationAnalysis 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 informationLecture 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 informationProblem 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 informationControl 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 informationLet H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 )
Review: Poles and Zeros of Fractional Form Let H() = P()/Q() be the system function of a rational form. Let us represent both P() and Q() as polynomials of (not - ) Then Poles: the roots of Q()=0 Zeros:
More informationLTI system response. Daniele Carnevale. Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata
LTI system response Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Fondamenti di Automatica e Controlli Automatici A.A. 2014-2015 1 / 15 Laplace
More information6.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 informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More informationModelling 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 informationTopic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis
Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability of
More informationTopic # Feedback Control
Topic #11 16.31 Feedback Control State-Space Systems State-space model features Observability Controllability Minimal Realizations Copyright 21 by Jonathan How. 1 Fall 21 16.31 11 1 State-Space Model Features
More informationAdvanced Control Theory
State Space Solution and Realization chibum@seoultech.ac.kr Outline State space solution 2 Solution of state-space equations x t = Ax t + Bu t First, recall results for scalar equation: x t = a x t + b
More informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
More informationTime-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationL 1 Adaptive Output Feedback Controller to Systems of Unknown
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 27 WeB1.1 L 1 Adaptive Output Feedback Controller to Systems of Unknown Dimension
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationDenis 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 informationNonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability
Nonlinear Control Lecture # 6 Passivity and Input-Output Stability Passivity: Memoryless Functions y y y u u u (a) (b) (c) Passive Passive Not passive y = h(t,u), h [0, ] Vector case: y = h(t,u), h T =
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More information2.3. VECTOR SPACES 25
2.3. VECTOR SPACES 25 2.3 Vector Spaces MATH 294 FALL 982 PRELIM # 3a 2.3. Let C[, ] denote the space of continuous functions defined on the interval [,] (i.e. f(x) is a member of C[, ] if f(x) is continuous
More informationBIBO STABILITY AND ASYMPTOTIC STABILITY
BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More information1. 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 informationMassachusetts 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 informationThe Important State Coordinates of a Nonlinear System
The Important State Coordinates of a Nonlinear System Arthur J. Krener 1 University of California, Davis, CA and Naval Postgraduate School, Monterey, CA ajkrener@ucdavis.edu Summary. We offer an alternative
More informationA 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 informationTheory of Linear Systems Exercises. Luigi Palopoli and Daniele Fontanelli
Theory of Linear Systems Exercises Luigi Palopoli and Daniele Fontanelli Dipartimento di Ingegneria e Scienza dell Informazione Università di Trento Contents Chapter. Exercises on the Laplace Transform
More informationLecture 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 informationAutomatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 18
More informationSolution 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 informationExam in Systems Engineering/Process Control
Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 27-6-2 Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total
More information10 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 informationMULTIVARIABLE ZEROS OF STATE-SPACE SYSTEMS
Copyright F.L. Lewis All rights reserved Updated: Monday, September 9, 8 MULIVARIABLE ZEROS OF SAE-SPACE SYSEMS If a system has more than one input or output, it is called multi-input/multi-output (MIMO)
More informationStability of Parameter Adaptation Algorithms. Big picture
ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about
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)
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 informationModule 9: State Feedback Control Design Lecture Note 1
Module 9: State Feedback Control Design Lecture Note 1 The design techniques described in the preceding lectures are based on the transfer function of a system. In this lecture we would discuss the state
More informationHankel 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 informationLyapunov Stability Analysis: Open Loop
Copyright F.L. Lewis 008 All rights reserved Updated: hursday, August 8, 008 Lyapunov Stability Analysis: Open Loop We know that the stability of linear time-invariant (LI) dynamical systems can be determined
More informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
More informationESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Discrete-Time Systems Hanz Richter, Professor Mechanical Engineering Department Cleveland State University Discrete-Time vs. Sampled-Data Systems A continuous-time
More informationEE 380. Linear Control Systems. Lecture 10
EE 380 Linear Control Systems Lecture 10 Professor Jeffrey Schiano Department of Electrical Engineering Lecture 10. 1 Lecture 10 Topics Stability Definitions Methods for Determining Stability Lecture 10.
More informationNonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control
Nonlinear Control Lecture # 14 Input-Output Stability L Stability Input-Output Models: y = Hu u(t) is a piecewise continuous function of t and belongs to a linear space of signals The space of bounded
More informationEngineering Tripos Part IIB Nonlinear Systems and Control. Handout 4: Circle and Popov Criteria
Engineering Tripos Part IIB Module 4F2 Nonlinear Systems and Control Handout 4: Circle and Popov Criteria 1 Introduction The stability criteria discussed in these notes are reminiscent of the Nyquist criterion
More information1 (30 pts) Dominant Pole
EECS C8/ME C34 Fall Problem Set 9 Solutions (3 pts) Dominant Pole For the following transfer function: Y (s) U(s) = (s + )(s + ) a) Give state space description of the system in parallel form (ẋ = Ax +
More informationLinear 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 information7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system
7 Stability 7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system ẋ(t) = A x(t), x(0) = x 0, A R n n, x 0 R n. (14) The origin x = 0 is a globally asymptotically
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More information6.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 informationMEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction
MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/2/214 Outline Solving State Equations Variation
More informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationControl 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 informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More information