Parameterization of Stabilizing Controllers for Interconnected Systems
|
|
- Hillary Stafford
- 5 years ago
- Views:
Transcription
1 Parameterization of Stabilizing Controllers for Interconnected Systems Michael Rotkowitz G D p D p D p D 1 G 2 G 3 G 4 p G 5 K 1 K 2 K 3 K 4 Contributors: Sanjay Lall, Randy Cogill K 5 Department of Signals, Sensors, and Systems Royal Institute of Technology (KTH) Control, Estimation, and Optimization of Interconnected Systems: From Theory to Industrial Applications CDC-ECC 05 Workshop, 11 December 2005
2 2 M. Rotkowitz, KTH Overview Motivation: Optimal Constrained Control Linear Time-Invariant Quadratic invariance Optimal Control over Networks Summation Nonlinear Time-Varying Iteration New condition
3 3 M. Rotkowitz, KTH Block Diagrams Two-Input Two-Output z y P 11 P 12 P 21 G u w z P 12 u K y P 11 P 21 w K G Classical K G
4 4 M. Rotkowitz, KTH Standard Formulation z y P 11 P 12 P 21 G K u w minimize subject to P 11 + P 12 K(I GK) 1 P 21 K stabilizes P
5 5 M. Rotkowitz, KTH Communicating Controllers G 1 G 2 G 3 G 4 G 5 D D D D D D D D K 1 K 2 K 3 K 4 K 5 Control design problem is to find K which is block tri-diagonal. u 1 K 11 DK y 1 u 2 u 3 u 4 = DK 21 K 22 DK y 2 0 DK 32 K 33 DK 34 0 y DK 43 K 44 DK 45 y 4 u DK 54 K 55 y 5
6 6 M. Rotkowitz, KTH General Formulation The set of K with a given decentralization constraint is a subspace S, called the information constraint. z P 11 P 12 P 21 G w y u We would like to solve K minimize subject to P 11 + P 12 K(I GK) 1 P 21 K stabilizes P K S For general P and S, there is no known tractable solution.
7 7 M. Rotkowitz, KTH Change of Variables - Stable Plant minimize subject to P 11 + P 12 K(I GK) 1 P 21 K stabilizes P Using the change of variables R = K(I GK) 1 we obtain the following equivalent problem minimize subject to P 11 + P 12 RP 21 R stable This is a convex optimization problem.
8 8 M. Rotkowitz, KTH Change of Variables - Stable Plant z P 11 w P 12 u K y R P 21 G Using the change of variables R = K(I GK) 1 we obtain the following equivalent problem z P 11 w P 12 R P 21 This is a convex optimization problem.
9 9 M. Rotkowitz, KTH Breakdown of Convexity - Stable Plant minimize subject to P 11 + P 12 K(I GK) 1 P 21 K stabilizes P K S Using the change of variables R = K(I GK) 1 we obtain the following equivalent problem minimize subject to P 11 + P 12 RP 21 R stable R(I + GR) 1 S
10 10 M. Rotkowitz, KTH Breakdown of Convexity - Stable Plant z P 11 w P 12 u K y R P 21 G Using the change of variables R = K(I GK) 1 we obtain the following equivalent problem minimize subject to P 11 + P 12 RP 21 R stable R(I + GR) 1 S
11 11 M. Rotkowitz, KTH Quadratic Invariance The set S is called quadratically invariant with respect to G if KGK S for all K S
12 12 M. Rotkowitz, KTH Quadratic Invariance The set S is called quadratically invariant with respect to G if KGK S for all K S Main Result S is quadratically invariant with respect to G if and only if K S K(I GK) 1 S
13 13 M. Rotkowitz, KTH Quadratic Invariance The set S is called quadratically invariant with respect to G if KGK S for all K S Main Result S is quadratically invariant with respect to G if and only if K S K(I GK) 1 S Parameterization {K K stabilizes P, K S} = { } R(I + GR) 1 R stable, R S
14 14 M. Rotkowitz, KTH Optimal Stabilizing Controller Suppose G R sp and S R p. We would like to solve minimize subject to P 11 + P 12 K(I GK) 1 P 21 K stabilizes P K S
15 15 M. Rotkowitz, KTH Optimal Stabilizing Controller Suppose G R sp and S R p. We would like to solve minimize subject to P 11 + P 12 K(I GK) 1 P 21 K stabilizes P K S If S is quadratically invariant with respect to G, we may solve which is convex. minimize subject to P 11 + P 12 RP 21 R stable R S
16 16 M. Rotkowitz, KTH Arbitrary Networks i t ij p ij j Assuming transmission delays satisfy the triangle inequality (i.e. transmissions take the quickest path) S is quadratically invariant with respect to G if t ij p ij for all i, j
17 17 M. Rotkowitz, KTH Distributed Control with Delays G D p D p D p D 1 G 2 G 3 G 4 p G 5 K 1 K 3 K 2 K 4 K 5 K 3 sees information from G 3 immediately G 2, G 4 after a delay of t G 1, G 5 after a delay of 2t G 3 is affected by inputs from K 3 immediately K 2, K 4 after a delay of p K 1, K 5 after a delay of 2p
18 18 M. Rotkowitz, KTH Distributed Control with Delays G D p D p D p D 1 G 2 G 3 G 4 p G 5 K 1 K 3 K 2 K 4 K 5 S is quadratically invariant with respect to G if t p
19 19 M. Rotkowitz, KTH Distributed Control with Delays G D p D p D p D 1 G 2 G 3 G 4 p G 5 D c D c D c D c D c D c D c K 1 K 3 K 2 K 4 K 5 K 3 sees information from G 3 after a delay of c G 2, G 4 after a delay of c + t G 1, G 5 after a delay of c + 2t G 3 is affected by inputs from K 3 immediately K 2, K 4 after a delay of p K 1, K 5 after a delay of 2p
20 20 M. Rotkowitz, KTH Distributed Control with Delays G D p D p D p D 1 G 2 G 3 G 4 p G 5 D c D c D c D c D c D c D c K 1 K 3 K 2 K 4 K 5 Without computational delay, S is quadratically invariant with respect to G if t p
21 21 M. Rotkowitz, KTH Distributed Control with Delays G D p D p D p D 1 G 2 G 3 G 4 p G 5 D c D c D c D c D c D c D c K 1 K 3 K 2 K 4 K 5 Without computational delay, S is quadratically invariant with respect to G if t p If computational delay is also present, then S is quadratically invariant with respect to G if t p + c n 1
22 22 M. Rotkowitz, KTH Two-Dimensional Lattice K 11 K 12 K 13 K 21 K 22 K 23 K 31 K 32 K 33 Assuming controllers communicate along edges
23 23 M. Rotkowitz, KTH Two-Dimensional Lattice K 11 K 12 K 13 G 11 D p G 12 D p G 13 D p D p D p K 21 G 21 D p K 22 G 22 D p K G D p D p D p K 31 K 32 K 33 G 31 D p G D 32 p G 33 Assuming controllers communicate along edges Assuming dynamics propagate along edges
24 24 M. Rotkowitz, KTH Two-Dimensional Lattice K 11 K 12 K 13 G 11 D p G 12 D p G 13 D p D p D p K 21 G 21 D p K 22 G 22 D p K G D p D p D p K 31 K 32 K 33 G 31 D p G D 32 p G 33 S is quadratically invariant with respect to G if t p
25 25 M. Rotkowitz, KTH Two-Dimensional Lattice K 11 K 12 K 13 K 21 K 22 K 23 K 31 K 32 K 33 Assuming controllers communicate along edges
26 26 M. Rotkowitz, KTH Two-Dimensional Lattice K 11 K 12 K 13 G 11 G 12 G 13 K 21 G 21 K 22 G 22 K G K 31 K 32 K 33 G 31 G 32 G 33 Assuming controllers communicate along edges Assuming dynamics propagate outward so that delay is proportional to geometric distance
27 27 M. Rotkowitz, KTH Two-Dimensional Lattice K 11 K 12 K 13 G 11 G 12 G 13 K 21 G 21 K 22 G 22 K G K 31 K 32 K 33 G 31 G 32 G 33 S is quadratically invariant with respect to G if t p 2
28 28 M. Rotkowitz, KTH Sparsity Example Suppose G S = K K For arbitrary K S KGK Hence S is quadratically invariant with respect to G.
29 29 M. Rotkowitz, KTH Same Structure Synthesis Suppose G S = K K Then KGK so S is not quadratically invariant with respect to G. In fact, K(I GK) 1
30 30 M. Rotkowitz, KTH Small Gain If a < 1 then (1 a) 1 = 1 + a + a 2 + a for example ( 1 1 ) 1 = more generally, if A < 1 then (1 A) 1 = I + A + A 2 + A
31 31 M. Rotkowitz, KTH Not So Small Gain Let Q = I + A + A 2 + A then AQ = A + A 2 + A subtracting we get (I A)Q = I and so Q = (I A) 1
32 32 M. Rotkowitz, KTH Convergence to Unstable Operator 25 Impulse Response 20 Amplitude Let W (s) = 2 s Time (sec) Plot shows impulse response of N k=0 W k for N = 1,..., 7 Converges to that of I I W = s 1 s+1 In what topology do the associated operators converge?
33 33 M. Rotkowitz, KTH Conditions for Convergence: Inert Very broad, reasonable class of plants and controllers Basically, impulse response must be finite at any time T Arbitrarily large Includes the case G R sp, K R p
34 34 M. Rotkowitz, KTH Proof Sketch K r y u G KGK S for all K S Suppose K S. Then K(I GK) 1 = K + KGK + K(GK) S
35 35 M. Rotkowitz, KTH Consider Nonlinear Let Q = I + A + A 2 + A then AQ = A + A 2 + A subtracting we get (I A)Q = I and so Q = (I A) 1
36 36 M. Rotkowitz, KTH Block Diagram Algebra K r y u G For a given r, we seek y, u such that y = r + Gu u = Ky and then define Y, R such that y = Y r = (I GK) 1 r u = Rr = K(I GK) 1 r
37 37 M. Rotkowitz, KTH Iteration of Signals K r y u G We can define the following iteration, commensurate with the diagram y (0) = r u (n) = Ky (n) y (n+1) = r + Gu (n)
38 38 M. Rotkowitz, KTH Iteration of Signals and Operators K r y u G We can define the following iterations, commensurate with the diagram y (0) = r u (n) = Ky (n) y (n+1) = r + Gu (n) Y (0) = I R (n) = KY (n) Y (n+1) = I + GR (n)
39 39 M. Rotkowitz, KTH Iteration of Signals K r y u G We then get the following recursions y (0) = r u (0) = Kr y (n+1) = r + GKy (n) u (n+1) = K(r + Gu (n) )
40 40 M. Rotkowitz, KTH Iteration of Signals K r y u G We then get the following recursions y (0) = r u (0) = Kr y (n+1) = r + GKy (n) u (n+1) = K(r + Gu (n) ).. y = lim n y(n) u = lim n u(n)
41 41 M. Rotkowitz, KTH Iteration of Operators K r y u G We then get the following recursions Y (0) = I R (0) = K Y (n+1) = I + GKY (n) R (n+1) = K(I + GR (n) ).. Y = lim n Y (n) R = lim n R(n)
42 42 M. Rotkowitz, KTH Conditions for Convergence Very broad, reasonable class of plants and controllers Arbitrarily large Includes the case G strictly causal, K causal Includes the case G R sp, K R p
43 43 M. Rotkowitz, KTH Parameterization - Nonlinear {K K stabilizes G} = { } R(I + GR) 1 R stable C.A. Desoer and R.W. Liu (1982) V. Anantharam and C.A. Desoer (1984)
44 44 M. Rotkowitz, KTH New Invariance Condition K 1 (I ± GK 2 ) S for all K 1, K 2 S
45 45 M. Rotkowitz, KTH New Invariance Condition K 1 (I ± GK 2 ) S for all K 1, K 2 S Invariance under Feedback If this condition is satisfied, then K S K(I GK) 1 S
46 46 M. Rotkowitz, KTH New Invariance Condition K 1 (I ± GK 2 ) S for all K 1, K 2 S Invariance under Feedback If this condition is satisfied, then K S K(I GK) 1 S Parameterization {K K stabilizes G, K S} = { } R(I + GR) 1 R stable, R S
47 47 M. Rotkowitz, KTH Proof Sketch K 1 (I ± GK 2 ) S K 1, K 2 S Suppose that K S. Then R (0) = K S. If we assume that R (n) S, then R (n+1) = K(I + GR (n) ) S Thus R (n) S for all n Z +. R = lim n R(n) S
48 48 M. Rotkowitz, KTH Conclusions Quadratic invariance allows parameterization of all (LTI) stablizing decentralized controllers. Similar condition allows parameterization of all (NLTV) stabilizing decentralized controllers These condition are satisfied when communications are faster than the propagation of dynamics.
49 49 M. Rotkowitz, KTH Quadratic Invariance The set S is called quadratically invariant with respect to G if KGK S for all K S????????? K 1 (I ± GK 2 ) S for all K 1, K 2 S
50 50 M. Rotkowitz, KTH Open Questions / Future Work Unstable plant Is there a weaker condition which achieves the same results? Perhaps K(I ± GK) S K S? When is the optimal (possibly nonlinear) controller linear? When (else) does it hold? What should we call it?
Information Structures Preserved Under Nonlinear Time-Varying Feedback
Information Structures Preserved Under Nonlinear Time-Varying Feedback Michael Rotkowitz Electrical Engineering Royal Institute of Technology (KTH) SE-100 44 Stockholm, Sweden Email: michael.rotkowitz@ee.kth.se
More informationMichael Rotkowitz 1,2
November 23, 2006 edited Parameterization of Causal Stabilizing Controllers over Networks with Delays Michael Rotkowitz 1,2 2006 IEEE Global Telecommunications Conference Abstract We consider the problem
More informationDecentralized Control Subject to Communication and Propagation Delays
Decentralized Control Subject to Communication and Propagation Delays Michael Rotkowitz 1,3 Sanjay Lall 2,3 IEEE Conference on Decision and Control, 2004 Abstract In this paper, we prove that a wide class
More informationA Characterization of Convex Problems in Decentralized Control
1984 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 51, NO. 2, FEBRUARY 2006 A Characterization of Convex Problems in Decentralized Control Michael Rotkowitz and Sanjay Lall Abstract We consider the problem
More informationarxiv: v1 [math.oc] 28 Sep 2011
1 On the Nearest Quadratically Invariant Information Constraint Michael C. Rotkowitz and Nuno C. Martins arxiv:1109.6259v1 [math.oc] 28 Sep 2011 Abstract Quadratic invariance is a condition which has been
More information274 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 51, NO. 2, FEBRUARY A Characterization of Convex Problems in Decentralized Control
274 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 51, NO 2, FEBRUARY 2006 A Characterization of Convex Problems in Decentralized Control Michael Rotkowitz Sanjay Lall Abstract We consider the problem of
More informationConvexity of optimal control over networks with delays and arbitrary topology
30 Int. J. Systems, Control and Communications, Vol. 2, Nos. 1/2/3, 2010 Convexity of optimal control over networks with delays and arbitrary topology Michael Rotkowitz* Department of Electrical and Electronic
More informationarxiv: v1 [cs.sy] 24 May 2013
Convexity of Decentralized Controller Synthesis Laurent Lessard Sanjay Lall arxiv:35.5859v [cs.sy] 4 May 3 Abstract In decentralized control problems, a standard approach is to specify the set of allowable
More informationTRACTABLE PROBLEMS IN OPTIMAL DECENTRALIZED CONTROL
TRACTABLE PROBLEMS IN OPTIMAL DECENTRALIZED CONTROL A DISSERTATION SUBMITTED TO THE DEPARTMENT OF AERONAUTICS AND ASTRONAUTICS AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT
More informationStructured 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 informationA Partial Order Approach to Decentralized Control
Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-11, 2008 A Partial Order Approach to Decentralized Control Parikshit Shah and Pablo A. Parrilo Department of Electrical
More informationDecentralized Control of Stochastic Systems
Decentralized Control of Stochastic Systems Sanjay Lall Stanford University CDC-ECC Workshop, December 11, 2005 2 S. Lall, Stanford 2005.12.11.02 Decentralized Control G 1 G 2 G 3 G 4 G 5 y 1 u 1 y 2 u
More informationA Convex Characterization of Multidimensional Linear Systems Subject to SQI Constraints
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 62, NO 6, JUNE 2017 2981 A Convex Characterization of Multidimensional Linear Systems Subject to SQI Constraints Şerban Sabău, Nuno C Martins, and Michael C
More informationAn Optimization-based Approach to Decentralized Assignability
2016 American Control Conference (ACC) Boston Marriott Copley Place July 6-8, 2016 Boston, MA, USA An Optimization-based Approach to Decentralized Assignability Alborz Alavian and Michael Rotkowitz Abstract
More informationA State-Space Approach to Control of Interconnected Systems
A State-Space Approach to Control of Interconnected Systems Part II: General Interconnections Cédric Langbort Center for the Mathematics of Information CALIFORNIA INSTITUTE OF TECHNOLOGY clangbort@ist.caltech.edu
More informationLecture 2 Discrete-Time LTI Systems: Introduction
Lecture 2 Discrete-Time LTI Systems: Introduction Outline 2.1 Classification of Systems.............................. 1 2.1.1 Memoryless................................. 1 2.1.2 Causal....................................
More informationA Convex Characterization of Distributed Control Problems in Spatially Invariant Systems with Communication Constraints
A Convex Characterization of Distributed Control Problems in Spatially Invariant Systems with Communication Constraints Bassam Bamieh Petros G. Voulgaris Revised Dec, 23 Abstract In this paper we consider
More informationOptimal H Control Design under Model Information Limitations and State Measurement Constraints
Optimal H Control Design under Model Information Limitations and State Measurement Constraints F. Farokhi, H. Sandberg, and K. H. Johansson ACCESS Linnaeus Center, School of Electrical Engineering, KTH-Royal
More informationZ-Transform. x (n) Sampler
Chapter Two A- Discrete Time Signals: The discrete time signal x(n) is obtained by taking samples of the analog signal xa (t) every Ts seconds as shown in Figure below. Analog signal Discrete time signal
More informationUNIT-II Z-TRANSFORM. This expression is also called a one sided z-transform. This non causal sequence produces positive powers of z in X (z).
Page no: 1 UNIT-II Z-TRANSFORM The Z-Transform The direct -transform, properties of the -transform, rational -transforms, inversion of the transform, analysis of linear time-invariant systems in the -
More informationDigital Signal Processing, Homework 1, Spring 2013, Prof. C.D. Chung
Digital Signal Processing, Homework, Spring 203, Prof. C.D. Chung. (0.5%) Page 99, Problem 2.2 (a) The impulse response h [n] of an LTI system is known to be zero, except in the interval N 0 n N. The input
More informationAM 205: lecture 19. Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods Optimality Conditions: Equality Constrained Case As another example of equality
More informationAM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods Quasi-Newton Methods General form of quasi-newton methods: x k+1 = x k α
More informationLARGE-SCALE networked systems have emerged in extremely
SUBMITTED TO IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. XX, JANUARY 2017 1 Separable and Localized System Level Synthesis for Large-Scale Systems Yuh-Shyang Wang, Member, IEEE, Nikolai Matni,
More informationControllers design for two interconnected systems via unbiased observers
Preprints of the 19th World Congress The nternational Federation of Automatic Control Cape Town, South Africa. August 24-29, 214 Controllers design for two interconnected systems via unbiased observers
More informationAustralian National University WORKSHOP ON SYSTEMS AND CONTROL
Australian National University WORKSHOP ON SYSTEMS AND CONTROL Canberra, AU December 7, 2017 Australian National University WORKSHOP ON SYSTEMS AND CONTROL A Distributed Algorithm for Finding a Common
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 information(Refer Slide Time: )
Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi FIR Lattice Synthesis Lecture - 32 This is the 32nd lecture and our topic for
More informationCosc 3451 Signals and Systems. What is a system? Systems Terminology and Properties of Systems
Cosc 3451 Signals and Systems Systems Terminology and Properties of Systems What is a system? an entity that manipulates one or more signals to yield new signals (often to accomplish a function) can be
More informationLecture 19 IIR Filters
Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class
More informationContents. 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 informationProblems in VLSI design
Problems in VLSI design wire and transistor sizing signal delay in RC circuits transistor and wire sizing Elmore delay minimization via GP dominant time constant minimization via SDP placement problems
More informationMIMO analysis: loop-at-a-time
MIMO robustness MIMO analysis: loop-at-a-time y 1 y 2 P (s) + + K 2 (s) r 1 r 2 K 1 (s) Plant: P (s) = 1 s 2 + α 2 s α 2 α(s + 1) α(s + 1) s α 2. (take α = 10 in the following numerical analysis) Controller:
More informationDecentralized LQG Control of Systems with a Broadcast Architecture
Decentralized LQG Control of Systems with a Broadcast Architecture Laurent Lessard 1,2 IEEE Conference on Decision and Control, pp. 6241 6246, 212 Abstract In this paper, we consider dynamical subsystems
More informationK 1 K 2. System Level Synthesis: A Tutorial. John C. Doyle, Nikolai Matni, Yuh-Shyang Wang, James Anderson, and Steven Low
System Level Synthesis: A Tutorial John C. Doyle, Nikolai Matni, Yuh-Shyang Wang, James Anderson, and Steven Low Abstract This tutorial paper provides an overview of the System Level Approach to control
More informationTHE PROBLEM of designing control systems where multiple
308 IEEE TRANSACTIONS ON CONTROL OF NETWORK SYSTEMS, VOL. 1, NO. 4, DECEMBER 2014 An Algebraic Approach to the Control of Decentralized Systems Laurent Lessard, Member, IEEE, and Sanjay Lall, Senior Member,
More information3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE
3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3.0 INTRODUCTION The purpose of this chapter is to introduce estimators shortly. More elaborated courses on System Identification, which are given
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 informationIII. Time Domain Analysis of systems
1 III. Time Domain Analysis of systems Here, we adapt properties of continuous time systems to discrete time systems Section 2.2-2.5, pp 17-39 System Notation y(n) = T[ x(n) ] A. Types of Systems Memoryless
More informationParametrization of All Strictly Causal Stabilizing Controllers of Multidimensional Systems single-input single-output case
Parametrization of All Strictly Causal Stabilizing Controllers of Multidimensional Systems single-input single-output case K. Mori Abstract We give a parametrization of all strictly causal stabilizing
More informationHere 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 informationDifferential and Difference LTI systems
Signals and Systems Lecture: 6 Differential and Difference LTI systems Differential and difference linear time-invariant (LTI) systems constitute an extremely important class of systems in engineering.
More informationLecture 11 FIR Filters
Lecture 11 FIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/4/12 1 The Unit Impulse Sequence Any sequence can be represented in this way. The equation is true if k ranges
More informationModule 09 Decentralized Networked Control Systems: Battling Time-Delays and Perturbations
Module 09 Decentralized Networked Control Systems: Battling Time-Delays and Perturbations Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/
More informationDiscrete-Time Systems
FIR Filters With this chapter we turn to systems as opposed to signals. The systems discussed in this chapter are finite impulse response (FIR) digital filters. The term digital filter arises because these
More informationDigital Signal Processing Lecture 4
Remote Sensing Laboratory Dept. of Information Engineering and Computer Science University of Trento Via Sommarive, 14, I-38123 Povo, Trento, Italy Digital Signal Processing Lecture 4 Begüm Demir E-mail:
More information1.7 Inequalities. Copyright Cengage Learning. All rights reserved.
1.7 Inequalities Copyright Cengage Learning. All rights reserved. Objectives Solving Linear Inequalities Solving Nonlinear Inequalities Absolute Value Inequalities Modeling with Inequalities 2 Inequalities
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 informationStability of Switched Linear Hyperbolic Systems by Lyapunov Techniques
2196 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 59, NO. 8, AUGUST 2014 Stability of Switched Linear Hyperbolic Systems by Lyapunov Techniques Christophe Prieur, Antoine Girard, Emmanuel Witrant Abstract
More informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More informationNonlinear Optimization for Optimal Control
Nonlinear Optimization for Optimal Control Pieter Abbeel UC Berkeley EECS Many slides and figures adapted from Stephen Boyd [optional] Boyd and Vandenberghe, Convex Optimization, Chapters 9 11 [optional]
More informationRobust Control. 2nd class. Spring, 2018 Instructor: Prof. Masayuki Fujita (S5-303B) Tue., 17th April, 2018, 10:45~12:15, S423 Lecture Room
Robust Control Spring, 2018 Instructor: Prof. Masayuki Fujita (S5-303B) 2nd class Tue., 17th April, 2018, 10:45~12:15, S423 Lecture Room 2. Nominal Performance 2.1 Weighted Sensitivity [SP05, Sec. 2.8,
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationNAME: 13 February 2013 EE301 Signals and Systems Exam 1 Cover Sheet
NAME: February EE Signals and Systems Exam Cover Sheet Test Duration: 75 minutes. Coverage: Chaps., Open Book but Closed Notes. One 8.5 in. x in. crib sheet Calculators NOT allowed. This test contains
More informationZ Transform (Part - II)
Z Transform (Part - II). The Z Transform of the following real exponential sequence x(nt) = a n, nt 0 = 0, nt < 0, a > 0 (a) ; z > (c) for all z z (b) ; z (d) ; z < a > a az az Soln. The given sequence
More informationSolutions to homework assignment #7 Math 119B UC Davis, Spring for 1 r 4. Furthermore, the derivative of the logistic map is. L r(x) = r(1 2x).
Solutions to homework assignment #7 Math 9B UC Davis, Spring 0. A fixed point x of an interval map T is called superstable if T (x ) = 0. Find the value of 0 < r 4 for which the logistic map L r has a
More informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
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 informationTeam Decision Theory and Information Structures in Optimal Control Problems
Team Decision Theory and Information Structures in Optimal Control Problems YU-CHI HO and K AI-CHING CHU Presented by: Ali Zaidi, Comm. Theory, KTH October 4, 2011 YU-CHI HO and K AI-CHING CHU Presented
More informationA Guaranteed Cost LMI-Based Approach for Event-Triggered Average Consensus in Multi-Agent Networks
A Guaranteed Cost LMI-Based Approach for Event-Triggered Average Consensus in Multi-Agent Networks Amir Amini, Arash Mohammadi, Amir Asif Electrical and Computer Engineering,, Montreal, Canada. Concordia
More informationSolution 10 July 2015 ECE301 Signals and Systems: Midterm. Cover Sheet
Solution 10 July 2015 ECE301 Signals and Systems: Midterm Cover Sheet Test Duration: 60 minutes Coverage: Chap. 1,2,3,4 One 8.5" x 11" crib sheet is allowed. Calculators, textbooks, notes are not allowed.
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 informationLTI Systems (Continuous & Discrete) - Basics
LTI Systems (Continuous & Discrete) - Basics 1. A system with an input x(t) and output y(t) is described by the relation: y(t) = t. x(t). This system is (a) linear and time-invariant (b) linear and time-varying
More informationUniqueness of Generalized Equilibrium for Box Constrained Problems and Applications
Uniqueness of Generalized Equilibrium for Box Constrained Problems and Applications Alp Simsek Department of Electrical Engineering and Computer Science Massachusetts Institute of Technology Asuman E.
More informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture 1: Course Overview; Discrete-Time Signals & Systems Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 2008 K. E.
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 informationEEL3135: Homework #4
EEL335: Homework #4 Problem : For each of the systems below, determine whether or not the system is () linear, () time-invariant, and (3) causal: (a) (b) (c) xn [ ] cos( 04πn) (d) xn [ ] xn [ ] xn [ 5]
More informationThe z-transform Part 2
http://faculty.kfupm.edu.sa/ee/muqaibel/ The z-transform Part 2 Dr. Ali Hussein Muqaibel The material to be covered in this lecture is as follows: Properties of the z-transform Linearity Initial and final
More informationCh 2: Linear Time-Invariant System
Ch 2: Linear Time-Invariant System A system is said to be Linear Time-Invariant (LTI) if it possesses the basic system properties of linearity and time-invariance. Consider a system with an output signal
More informationLifted approach to ILC/Repetitive Control
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
More informationQ-Parametrization and an SDP for H -optimal Decentralized Control
4th IFAC Workshop on Distributed Estimation and Control in Networked Systems September 25-26, 2013. Rhine-Moselle-Hall, Koblenz, Germany Q-Parametrization and an SDP for H -optimal Decentralized Control
More informationChapter 2 Time-Domain Representations of LTI Systems
Chapter 2 Time-Domain Representations of LTI Systems 1 Introduction Impulse responses of LTI systems Linear constant-coefficients differential or difference equations of LTI systems Block diagram representations
More informationOptimal Control of Spatially Distributed Systems
University of Pennsylvania ScholarlyCommons Lab Papers (GRASP) General Robotics, Automation, Sensing and Perception Laboratory 8-2008 Optimal Control of Spatially Distributed Systems Nader Motee University
More informationBasic concepts in DT systems. Alexandra Branzan Albu ELEC 310-Spring 2009-Lecture 4 1
Basic concepts in DT systems Alexandra Branzan Albu ELEC 310-Spring 2009-Lecture 4 1 Readings and homework For DT systems: Textbook: sections 1.5, 1.6 Suggested homework: pp. 57-58: 1.15 1.16 1.18 1.19
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationStrong stability of neutral equations with dependent delays
Strong stability of neutral equations with dependent delays W Michiels, T Vyhlidal, P Zitek, H Nijmeijer D Henrion 1 Introduction We discuss stability properties of the linear neutral equation p 2 ẋt +
More informationExamples. 2-input, 1-output discrete-time systems: 1-input, 1-output discrete-time systems:
Discrete-Time s - I Time-Domain Representation CHAPTER 4 These lecture slides are based on "Digital Signal Processing: A Computer-Based Approach, 4th ed." textbook by S.K. Mitra and its instructor materials.
More informationStructured Coprime Factorizations Description of Linear and Time Invariant Networks
52nd IEEE Conference on Decision and Control December 10-13, 2013. Florence, Italy Structured Coprime Factorizations Description of Linear and Time Invariant Networks Şerban Sabău, Cristian ară, Sean Warnick
More informationAdaptive Filtering Part II
Adaptive Filtering Part II In previous Lecture we saw that: Setting the gradient of cost function equal to zero, we obtain the optimum values of filter coefficients: (Wiener-Hopf equation) Adaptive Filtering,
More informationEE 3054: Signals, Systems, and Transforms Spring A causal discrete-time LTI system is described by the equation. y(n) = 1 4.
EE : Signals, Systems, and Transforms Spring 7. A causal discrete-time LTI system is described by the equation Test y(n) = X x(n k) k= No notes, closed book. Show your work. Simplify your answers.. A discrete-time
More informationGeometric Interpolation by Planar Cubic Polynomials
1 / 20 Geometric Interpolation by Planar Cubic Polynomials Jernej Kozak, Marjeta Krajnc Faculty of Mathematics and Physics University of Ljubljana Institute of Mathematics, Physics and Mechanics Avignon,
More informationECE 301 Division 1 Exam 1 Solutions, 10/6/2011, 8-9:45pm in ME 1061.
ECE 301 Division 1 Exam 1 Solutions, 10/6/011, 8-9:45pm in ME 1061. Your ID will be checked during the exam. Please bring a No. pencil to fill out the answer sheet. This is a closed-book exam. No calculators
More informationLecture 6. Chapter 8: Robust Stability and Performance Analysis for MIMO Systems. Eugenio Schuster.
Lecture 6 Chapter 8: Robust Stability and Performance Analysis for MIMO Systems Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 6 p. 1/73 6.1 General
More informationLyapunov Stability of Linear Predictor Feedback for Distributed Input Delays
IEEE TRANSACTIONS ON AUTOMATIC CONTROL VOL. 56 NO. 3 MARCH 2011 655 Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays Nikolaos Bekiaris-Liberis Miroslav Krstic In this case system
More informationSYMMETRIC SUBGROUP ACTIONS ON ISOTROPIC GRASSMANNIANS
1 SYMMETRIC SUBGROUP ACTIONS ON ISOTROPIC GRASSMANNIANS HUAJUN HUANG AND HONGYU HE Abstract. Let G be the group preserving a nondegenerate sesquilinear form B on a vector space V, and H a symmetric subgroup
More informationSignals and Systems Chapter 2
Signals and Systems Chapter 2 Continuous-Time Systems Prof. Yasser Mostafa Kadah Overview of Chapter 2 Systems and their classification Linear time-invariant systems System Concept Mathematical transformation
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 informationMathematics for Control Theory
Mathematics for Control Theory Outline of Dissipativity and Passivity Hanz Richter Mechanical Engineering Department Cleveland State University Reading materials Only as a reference: Charles A. Desoer
More informationDynamic measurement: application of system identification in metrology
1 / 25 Dynamic measurement: application of system identification in metrology Ivan Markovsky Dynamic measurement takes into account the dynamical properties of the sensor 2 / 25 model of sensor as dynamical
More informationTheory in Model Predictive Control :" Constraint Satisfaction and Stability!
Theory in Model Predictive Control :" Constraint Satisfaction and Stability Colin Jones, Melanie Zeilinger Automatic Control Laboratory, EPFL Example: Cessna Citation Aircraft Linearized continuous-time
More informationIntroduction to DSP Time Domain Representation of Signals and Systems
Introduction to DSP Time Domain Representation of Signals and Systems Dr. Waleed Al-Hanafy waleed alhanafy@yahoo.com Faculty of Electronic Engineering, Menoufia Univ., Egypt Digital Signal Processing (ECE407)
More informationNumerical Sequences and Series
Numerical Sequences and Series Written by Men-Gen Tsai email: b89902089@ntu.edu.tw. Prove that the convergence of {s n } implies convergence of { s n }. Is the converse true? Solution: Since {s n } is
More informationDistributed and Real-time Predictive Control
Distributed and Real-time Predictive Control Melanie Zeilinger Christian Conte (ETH) Alexander Domahidi (ETH) Ye Pu (EPFL) Colin Jones (EPFL) Challenges in modern control systems Power system: - Frequency
More informationEigenvalue problems and optimization
Notes for 2016-04-27 Seeking structure For the past three weeks, we have discussed rather general-purpose optimization methods for nonlinear equation solving and optimization. In practice, of course, we
More informationVariable-gain output feedback control
7. Variable-gain output feedback control 7.1. Introduction PUC-Rio - Certificação Digital Nº 611865/CA In designing control laws, the usual first step is to describe the plant at a given operating point
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 18
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 18 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory May 31, 2012 Andre Tkacenko
More informationLet 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 informationConsensus Problems on Small World Graphs: A Structural Study
Consensus Problems on Small World Graphs: A Structural Study Pedram Hovareshti and John S. Baras 1 Department of Electrical and Computer Engineering and the Institute for Systems Research, University of
More informationFurther Results on Model Structure Validation for Closed Loop System Identification
Advances in Wireless Communications and etworks 7; 3(5: 57-66 http://www.sciencepublishinggroup.com/j/awcn doi:.648/j.awcn.735. Further esults on Model Structure Validation for Closed Loop System Identification
More informationVU Signal and Image Processing
052600 VU Signal and Image Processing Torsten Möller + Hrvoje Bogunović + Raphael Sahann torsten.moeller@univie.ac.at hrvoje.bogunovic@meduniwien.ac.at raphael.sahann@univie.ac.at vda.cs.univie.ac.at/teaching/sip/18s/
More information