LMI Methods in Optimal and Robust Control
|
|
- Anissa Griffith
- 5 years ago
- Views:
Transcription
1 LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 14: LMIs for Robust Control in the LF Framework
2 ypes of Uncertainty In this Lecture, we will cover Unstructured, Dynamic, norm-bounded: Structured, Static, norm-bounded: := { L(L 2 ) : H < 1} := {diag(δ 1,, δ K, 1, N ) : δ i < 1, σ( i ) < 1} Structured, Dynamic, norm-bounded: := { 1, 2, L(L 2 ) : i H < 1} Unstructured, Static, norm-bounded: Parametric, Polytopic: := { R n n : 1} := { R n n : = i Parametric, Interval: := { i α i H i, α i 0, i i δ i : δ i [δ i, δ+ i ]} Each of these can be ime-varying or ime-invariant! α i = 1} M. Peet Lecture 14: 2 / 28
3 Back to the Linear Fractional ransformation he interval and polytopic cases rely on Linearity of the uncertain parameters. ẋ(t) = (A 0 + (t))x(t) he Linear-Fractional ransformation, however [ ] [ ] [ ] ẋ1 = z(t) S(P, x1 (t) ) = (P F (t) 22 +P 21 (I P 11 ) 1 x(t) P 12 ) F (t) is an arbitrary rational function. We focus on two results: he S-Procedure for Unstructured Uncertainty Sets he Structured Singular Value for Structured Uncertainty Sets. p M q M. Peet Lecture 14: 3 / 28
4 Robust Stability p M q Questions: Is S(, M) stable for all? Is I M 11 invertible for all? M. Peet Lecture 14: 4 / 28
5 Redefine Robust and Quadratic Stability Suppose we have the system Definition 1. [ ] M11 M M = 12 M 21 M 22 he pair (M, ) is Robustly Stable if (I M 22 ) is invertible for all. [ẋ(t) ] [ ] Alternatively, if = z(t) S(M, x(t) ) w(t) Definition 2. he pair (M, ) is Robustly Stable if for some β > 0, ρ(m 22 + M 21 (I M 11 ) 1 M 12 ) + βi is Hurwitz for all. [ ] [ ] Alternatively, if xk+1 = z S(M, xk ) k w k Definition 3. he pair (M, ) is Robustly Stable if ρ(m 22 + M 21 (I M 11 ) 1 M 12 ) = β < 1 for all. M. Peet Lecture 14: 5 / 28
6 Quadratic Stability - Parametric Uncertainty Focus on the 1,1 block of S(M, ): If ẋ(t) = S(M, )x(t), Definition 4. he pair (M, ) is Quadratically Stable if there exists a P > 0 such that S(M, ) P + P S(M, ) < βi for all Alternatively, if x k+1 = S(M, )x k, Definition 5. he pair (M, ) is Quadratically Stable if there exists a P > 0 such that S(M, ) P S(M, ) P < βi for all for all. M. Peet Lecture 14: 6 / 28
7 Parametric, Norm-Bounded ime-varying Uncertainty Consider the state-space representation: ẋ(t) = Ax(t) + Mp(t), q(t) = Nx(t) + Qp(t), p(t) = (t)q(t), (t) Parametric, Norm-Bounded Uncertainty: := { R n n : 1} M. Peet Lecture 14: 7 / 28
8 Parametric, Norm-Bounded Uncertainty Quadratic Stability: here exists a P > 0 such that P (Ax(t) + Mp) + (Ax(t) + Mp) P < 0 for all p {p : p = q, q = Nx + Qp} heorem 6. he system ẋ(t) = Ax(t) + Mp(t), p(t) = (t)q(t), q(t) = Nx(t) + Qp(t), := { R n n : 1} is quadratically stable if and only if there exists some P > 0 such that [ ] [ ] [ ] x A P + P A P M x p M < 0 P 0 p [ ] {[ ] [ ] [ ] [ ] x x x N for all : N N Q x p p p Q N I Q Q p } 0 M. Peet Lecture 14: 8 / 28
9 Parametric, Norm-Bounded Uncertainty If. If [ ] [ x A P + P A P M p M P 0 then for all x, p such that ] [ x < 0 p] [ ] x for all p {[ ] [ ] [ [ ] } x x N : N N Q x p p Q N I Q Q] 0 p x P (Ax + Mp) + (Ax + Mp) P x < 0 p 2 Nx + Qp 2 herefore, since p = q implies p q, we have quadratic stability. he only if direction is similar. M. Peet Lecture 14: 9 / 28
10 he S-Procedure A Significant LMI for your oolbox Quadratic stability here requires positivity of a matrix on a subset. his is Generally a very hard problem NP-hard to determine if x F x 0 for all x 0. (Matrix Copositivity) S-procedure to the rescue! he S-procedure asks the question: Is z F z 0 for all z {x : x Gx 0}? Corollary 7 (S-Procedure). z F z 0 for all z {x : x Gx 0} if there exists a scalar τ 0 such that F τg 0. he S-procedure is Necessary if {x : x Gx > 0}. M. Peet Lecture 14: 10 / 28
11 Parametric, Norm-Bounded Uncertainty heorem 8 (Dual Version). he system ẋ(t) = Ax(t) + Mp(t), p(t) = (t)q(t), q(t) = Nx(t) + Qp(t), := { R n n : 1} is quadratically stable if and only if there exists some µ 0 and P > 0 such that [ ] [ ] AP + P A P N MM MQ + µ NP 0 QM QQ < 0} I Noting that the LMI can be written as [ ] [ ] [ ] AP + P A P N M M + µ < 0 NP µi Q Q or AP + P A P N M NP µi Q < 0 M Q 1 µ I we see that this condition is simply an H gain condition on the nominal system H < 1. M. Peet Lecture 14: 11 / 28
12 Necessity of the Small-Gain Condition his leads to the interesting result: p M q If := { L(L 2 ) : 1}, then S(P, ) H if and only if P 11 H < 1 he small gain condition is necessary and sufficient for stability. Quadratic Stability is equivalent to stability. Holds for Dynamic and Parametric Uncertainty M. Peet Lecture 14: 12 / 28
13 Quadratic Stability and Equivalence to Robust Stability Consider Quadratic Stability in Discrete-ime: x k+1 = S l (M, )x k. Definition 9. (S l, ) is QS if S l (M, ) P S l (M, ) P < 0 for all heorem 10 (Packard and Doyle). Let M R (n+m) (n+m) be given with ρ(m 11 ) 1 and σ(m 22 ) < 1. hen the following are equivalent. 1. he pair (M, = R m m ) is quadratically stable. 2. he pair (M, = C m m ) is quadratically stable. 3. he pair (M, = C m m ) is robustly stable. M. Peet Lecture 14: 13 / 28
14 Quadratically Stabilizing Controllers with Parametric Norm-Bounded Uncertainty However, we can add controllers: heorem 11. he system with u(t) = Kx(t) and ẋ(t) = A 0 x(t) + Bu(t) + Mp(t), p(t) = (t)q(t), q(t) = Nx(t) + Qp(t) + D 12 u(t), := { R n n : 1} is quadratically stable if and only if there exists some µ 0 and P > 0 such that [ ] [ ] (A + BK)P + P (A + BK) P (N + D 12 K) MM MQ +µ (N + D 12 K)P 0 QM QQ < 0} I Of course, this is bilinear in P and K, so we make the change of variables Z = KP. M. Peet Lecture 14: 14 / 28
15 An LMI for Quadratically Stabilizing Controllers with Parametric Norm-Bounded Uncertainty heorem 12. here exists a K such that the system with u(t) = Kx(t) ẋ(t) = Ax(t) + Bu(t) + Mp(t), p(t) = (t)q(t), q(t) = Nx(t) + Qp(t) + D 12 u(t), := { R n n : 1} is quadratically stable if and only if there exists some µ 0, Z and P > 0 such that [ ] [ ] AP + BZ + P A + Z B P N + Z D12 MM MQ +µ NP + D 12 Z 0 QM QQ < 0}. I hen K = ZP 1 is a quadratically stabilizing controller. We can also extend this result to optimal control in the H norm. M. Peet Lecture 14: 15 / 28
16 An LMI for H -Optimal Quadratically Stabilizing Controllers with Parametric Norm-Bounded Uncertainty In this case, we set Q = 0. heorem 13. here exists a K such that the system with u(t) = Kx(t) ẋ(t) = Ax(t) + Bu(t) + Mp(t) + B 2 w(t), p(t) = (t)q(t), q(t) = Nx(t) + D 12 u(t), := { R n n : 1} y(t) = Cx(t) + D 22 u(t) satisfies y L2 γ u L2 if there exists some µ 0, Z and P > 0 such that AP + BZ + P A + Z B + µmm (CP + D 22 Z) P N + Z D12 CP + D 22 Z γ 2 I 0 < 0. NP + D 12 Z 0 µi hen K = ZP 1 is the corresponding controller. M. Peet Lecture 14: 16 / 28
17 Structure, Norm-Bounded Uncertainty For the case of structured parametric uncertainty, we define the structured set = { = diag(δ 1 I n1,, δ s I ns, s+1,, s+f ) : δ i R, R n k n k } δ and represent unknown parameters. s is the number of scalar parameters. f is the number of matrix parameters. M. Peet Lecture 14: 17 / 28
18 he Structured Singular Value For the case of structured parametric uncertainty, we define the structured singular value. Definition 14. Given system M L(L 2 ) and set as above, we define the Structured Singular Value of (M, ) as µ(m, ) = inf 1 I M is singular Of course, S(M, ) is stable if and only if µ(m 11, ) < 1. Obviously, µ(m, ) < M For := { L(L 2 ) : 1}, µ(m, ) = M µ(αm, ) = α µ(m, ) 1 Can increase M by a factor µ(m, ) before losing stability. In general, computing µ is NP-hard unless uncertainty is unstructured or block-diagonal. M. Peet Lecture 14: 18 / 28
19 Scalings and he Structured Singular Value Suppose Θ = {Θ : Θ = Θ for all } hen µ(m, ) = inf Θ Θ ΘMΘ 1. Θ is the set of scalings. M. Peet Lecture 14: 19 / 28
20 Scalings and he Structured Singular Value = { = diag(δ 1 I n1,, δ s I ns, s+1,, s+f ) : δ i R, R n k n k } Define the set of scalings PΘ := {diag(θ 1,, Θ s, θ s+1 I,, θ s+f I : Θ i > 0, θ j > 0} heorem 15. Suppose system M has transfer function ˆM(s) = C(sI A) 1 B + D with ˆM H. he following are equivalent here exists Θ Θ such that ΘMΘ 1 2 < γ. here exists Θ PΘ and X > 0 such that [ ] A X + XA XB B + 1 [ ] C X Θ γ 2 D Θ [ C D ] < 0 Note: o minimize γ, you must use bisection. M. Peet Lecture 14: 20 / 28
21 An LMI for Stability of Structured, Norm-Bounded Uncertainty his allows us to generalize the S-procedure to structured uncertainty heorem 16. he system ẋ(t) = Ax(t) + Mp(t), p(t) = (t)q(t), q(t) = Nx(t) + Qp(t),, 1 is quadratically stable if and only if there exists some Θ PΘ and P > 0 such that [ ] [ ] AP + P A P N MΘM MΘQ + NP 0 QΘM QΘQ < 0} Θ his is an LMI in Θ and P. he constraint Θ PΘ is linear PΘ := {diag(θ 1,, Θ s, θ s+1 I,, θ s+f I) : Θ i > 0, θ j > 0} M. Peet Lecture 14: 21 / 28
22 An LMI for Stability with Structured, Norm-Bounded Uncertainty o prove the theorem, we can take a closer look at the scalings: Since = for PΘ, the system can equivalently be written as ẋ(t) = Ax(t) + M 1 p(t), p(t) = (t)q(t), q(t) = Nx(t) + Q 1 p(t),, 1 for any PΘ. hen [ ] [ ] AP + P A P N MM MQ + NP 0 QM QQ < 0 I becomes [ AP + P A P N NP 0 Pre- and Post-multiplying by the LMI condition. ] [ M + 2 M M 2 Q Q 2 M Q 2 Q I [ I ] < 0} ], and using Θ = 2 PΘ, we recover M. Peet Lecture 14: 22 / 28
23 An LMI for Stabilizing State-Feedback Controllers with Structured Norm-Bounded Uncertainty heorem 17. here exists a K such that the system with u(t) = Kx(t) ẋ(t) = Ax(t) + Bu(t) + Mp(t), p(t) = (t)q(t), q(t) = Nx(t) + Qp(t) + D 12 u(t),, 1 is quadratically stable if and only if there exists some Θ PΘ, P > 0 and Z such that [ ] [ ] AP + BZ + P A + Z B P N + Z D12 MΘM MΘQ + NP + D 12 Z 0 QΘM QΘQ < 0. Θ hen K = ZP 1 is a quadratically stabilizing controller. We can also extend this result to optimal control in the H norm. M. Peet Lecture 14: 23 / 28
24 An LMI for Optimal State-Feedback Controllers with Structured Norm-Bounded Uncertainty In this case, we set Q = 0. heorem 18. here exists a K such that the system with u(t) = Kx(t) ẋ(t) = Ax(t) + Bu(t) + Mp(t) + B 2 w(t), q(t) = Nx(t) + D 12 u(t),, 1 y(t) = Cx(t) + D 22 u(t) p(t) = (t)q(t), satisfies y L2 γ u L2 if there exists some Θ PΘ, Z and P > 0 such that AP + BZ + P A + Z B + MΘM (CP + D 22 Z) P N + Z D12 CP + D 22 Z γ 2 I 0 < 0. NP + D 12 Z 0 Θ hen K = ZP 1 is the corresponding controller. M. Peet Lecture 14: 24 / 28
25 An LMI for Optimal State-Feedback Controllers with Structured Norm-Bounded Uncertainty Using the scaled system we get ẋ(t) = Ax(t) + Bu(t) + M 1 p(t) + B 2 w(t), q(t) = Nx(t) + D 12 u(t),, 1 y(t) = Cx(t) + D 22 u(t) p(t) = (t)q(t), AP + BZ + P A + Z B + M 2 M (CP + D 22 Z) P N + Z D 12 CP + D 22 Z γ 2 I 0 < 0. NP + D 12 Z 0 I I 0 0 Pre- and Post-multiplying by 0 I 0, and using Θ = 2 PΘ, we recover the LMI condition. M. Peet Lecture 14: 25 / 28
26 Output-Feedback Robust Controller Synthesis How to Solve the Output Feedback Case??? inf K sup S( S(G, ), K) H M. Peet Lecture 14: 26 / 28
27 D-K Iteration A Heuristic for Dynamic Output Feedback Synthesis Finally, we mention a Heuristic for Output-Feedback Controller synthesis. Initialize: Θ = I. Define: A B 1 Θ 1 2 B 2 Ĝ Θ (s) = Θ 1 2 C 1 Θ 1 2 D 11 Θ 1 2 Θ 1 2 D 12 C 2 D 21 Θ Step 1: Fix Θ and solve inf K S(G Θ, K) H Step 2: Fix K and minimize γ such that there exists Θ PΘ ( or Θ PΘ I if you include the regulated output channel.) and X > 0 such that [ ] A cl X + XA cl XB cl Bcl X Θ + 1 [ ] C cl γ 2 Θ [ ] C cl D cl < 0 D cl where A cl, B cl, C cl, D cl define S(G I, K). (Requires Bisection). Step 3: GOO Step 1 M. Peet Lecture 14: 27 / 28
28 A Word on D-K Iteration with Static Uncertainty A Heuristic for Dynamic Output Feedback Synthesis he D-K iteration outlined in this lecture is only valid for Dynamic Uncertainty: (t). Our Scalings Θ are time-invariant. For Static uncertainties, we should search for Dynamic Scaling Factors Θ(s) is a ransfer Function his is much harder to represent as an LMI (Or by any other method!). Matlab has built-in functionality, but it is hard to use. We will return to µ analysis for static uncertainties when we consider more advanced forms of optimization. M. Peet Lecture 14: 28 / 28
Modern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
More informationLinear Matrix Inequalities in Robust Control. Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002
Linear Matrix Inequalities in Robust Control Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002 Objective A brief introduction to LMI techniques for Robust Control Emphasis on
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 22: H 2, LQG and LGR Conclusion To solve the H -optimal state-feedback problem, we solve min γ such that γ,x 1,Y 1,A n,b n,c n,d
More informationLMIs 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 informationModern 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 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 informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 21: Optimal Output Feedback Control connection is called the (lower) star-product of P and Optimal Output Feedback ansformation (LFT).
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 4: LMIs for State-Space Internal Stability Solving the Equations Find the output given the input State-Space:
More informationFEL3210 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 informationarzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.1 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS STABILITY ANALYSIS Didier HENRION www.laas.fr/ henrion henrion@laas.fr Denis ARZELIER www.laas.fr/
More informationStructured singular value and µ-synthesis
Structured singular value and µ-synthesis Robust Control Course Department of Automatic Control, LTH Autumn 2011 LFT and General Framework z P w z M w K z = F u (F l (P,K), )w = F u (M, )w. - Last week
More informationOptimization of Polynomials
Optimization of Polynomials Matthew M. Peet Arizona State University Thanks to S. Lall and P. Parrilo for guidance and supporting material Lecture 03: Optimization of Polynomials Overview In this lecture,
More informationLecture 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 informationRobust Stability. Robust stability against time-invariant and time-varying uncertainties. Parameter dependent Lyapunov functions
Robust Stability Robust stability against time-invariant and time-varying uncertainties Parameter dependent Lyapunov functions Semi-infinite LMI problems From nominal to robust performance 1/24 Time-Invariant
More informationGrammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology
Grammians Matthew M. Peet Illinois Institute of Technology Lecture 2: Grammians Lyapunov Equations Proposition 1. Suppose A is Hurwitz and Q is a square matrix. Then X = e AT s Qe As ds is the unique solution
More informationRobust Multi-Objective Control for Linear Systems
Robust Multi-Objective Control for Linear Systems Elements of theory and ROMULOC toolbox Dimitri PEAUCELLE & Denis ARZELIER LAAS-CNRS, Toulouse, FRANCE Part of the OLOCEP project (includes GloptiPoly)
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 informationTwo Challenging Problems in Control Theory
Two Challenging Problems in Control Theory Minyue Fu School of Electrical Engineering and Computer Science University of Newcastle, NSW 2308 Australia Email: minyue.fu@newcastle.edu.au Abstract This chapter
More informationOutline. Linear Matrix Inequalities in Control. Outline. System Interconnection. j _jst. ]Bt Bjj. Generalized plant framework
Outline Linear Matrix Inequalities in Control Carsten Scherer and Siep Weiland 7th Elgersburg School on Mathematical Systems heory Class 3 1 Single-Objective Synthesis Setup State-Feedback Output-Feedback
More informationEE363 homework 8 solutions
EE363 Prof. S. Boyd EE363 homework 8 solutions 1. Lyapunov condition for passivity. The system described by ẋ = f(x, u), y = g(x), x() =, with u(t), y(t) R m, is said to be passive if t u(τ) T y(τ) dτ
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 20: LMI/SOS Tools for the Study of Hybrid Systems Stability Concepts There are several classes of problems for
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 informationMulti-Objective Robust Control of Rotor/Active Magnetic Bearing Systems
Multi-Objective Robust Control of Rotor/Active Magnetic Bearing Systems İbrahim Sina Kuseyri Ph.D. Dissertation June 13, 211 İ. Sina Kuseyri (B.U. Mech.E.) Robust Control of Rotor/AMB Systems June 13,
More informationLOW ORDER H CONTROLLER DESIGN: AN LMI APPROACH
LOW ORDER H CONROLLER DESIGN: AN LMI APPROACH Guisheng Zhai, Shinichi Murao, Naoki Koyama, Masaharu Yoshida Faculty of Systems Engineering, Wakayama University, Wakayama 640-8510, Japan Email: zhai@sys.wakayama-u.ac.jp
More informationModern Control Systems
Modern Control Systems Matthew M. Peet Arizona State University Lecture 09: Observability Observability For Static Full-State Feedback, we assume knowledge of the Full-State. In reality, we only have measurements
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 informationSemidefinite 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 informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:
More informationStability and Robustness 1
Lecture 2 Stability and Robustness This lecture discusses the role of stability in feedback design. The emphasis is notonyes/notestsforstability,butratheronhowtomeasurethedistanceto instability. The small
More informationME 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 informationChapter Robust Performance and Introduction to the Structured Singular Value Function Introduction As discussed in Lecture 0, a process is better desc
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 Robust
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 informationAnalysis of systems containing nonlinear and uncertain components by using Integral Quadratic Constraints
Analysis of systems containing nonlinear and uncertain components by using Integral Quadratic Constraints Tamas Peni MTA-SZTAKI October 2, 2012 Tamas Peni (MTA-SZTAKI) Seminar October 2, 2012 1 / 42 Contents
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 informationIntegral action in state feedback control
Automatic Control 1 in state feedback control Prof. Alberto Bemporad University of Trento Academic year 21-211 Prof. Alberto Bemporad (University of Trento) Automatic Control 1 Academic year 21-211 1 /
More informationMultiobjective Optimization Applied to Robust H 2 /H State-feedback Control Synthesis
Multiobjective Optimization Applied to Robust H 2 /H State-feedback Control Synthesis Eduardo N. Gonçalves, Reinaldo M. Palhares, and Ricardo H. C. Takahashi Abstract This paper presents an algorithm for
More informationACM/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 informationFEL3210 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 informationH 2 Suboptimal Estimation and Control for Nonnegative
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 FrC20.3 H 2 Suboptimal Estimation and Control for Nonnegative Dynamical Systems
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 informationControl 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 informationAnalysis of Systems with State-Dependent Delay
Analysis of Systems with State-Dependent Delay Matthew M. Peet Arizona State University Tempe, AZ USA American Institute of Aeronautics and Astronautics Guidance, Navigation and Control Conference Boston,
More informationLinear Matrix Inequality (LMI)
Linear Matrix Inequality (LMI) A linear matrix inequality is an expression of the form where F (x) F 0 + x 1 F 1 + + x m F m > 0 (1) x = (x 1,, x m ) R m, F 0,, F m are real symmetric matrices, and the
More informationRobust Input-Output Energy Decoupling for Uncertain Singular Systems
International Journal of Automation and Computing 1 (25) 37-42 Robust Input-Output Energy Decoupling for Uncertain Singular Systems Xin-Zhuang Dong, Qing-Ling Zhang Institute of System Science, Northeastern
More informationStatic Output Feedback Stabilisation with H Performance for a Class of Plants
Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,
More informationROBUST STABILITY AND PERFORMANCE ANALYSIS* [8 # ]
ROBUST STABILITY AND PERFORMANCE ANALYSIS* [8 # ] General control configuration with uncertainty [8.1] For our robustness analysis we use a system representation in which the uncertain perturbations are
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 information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
More informationResearch Article An Equivalent LMI Representation of Bounded Real Lemma for Continuous-Time Systems
Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 28, Article ID 67295, 8 pages doi:1.1155/28/67295 Research Article An Equivalent LMI Representation of Bounded Real Lemma
More informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
More informationTopic # 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 informationController Coefficient Truncation Using Lyapunov Performance Certificate
Controller Coefficient Truncation Using Lyapunov Performance Certificate Joëlle Skaf Stephen Boyd Information Systems Laboratory Electrical Engineering Department Stanford University European Control Conference,
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 02: Optimization (Convex and Otherwise) What is Optimization? An Optimization Problem has 3 parts. x F f(x) :
More informationRobust Anti-Windup Controller Synthesis: A Mixed H 2 /H Setting
Robust Anti-Windup Controller Synthesis: A Mixed H /H Setting ADDISON RIOS-BOLIVAR Departamento de Sistemas de Control Universidad de Los Andes Av. ulio Febres, Mérida 511 VENEZUELA SOLBEN GODOY Postgrado
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 informationLecture: Convex Optimization Problems
1/36 Lecture: Convex Optimization Problems http://bicmr.pku.edu.cn/~wenzw/opt-2015-fall.html Acknowledgement: this slides is based on Prof. Lieven Vandenberghe s lecture notes Introduction 2/36 optimization
More informationEvent-triggered control subject to actuator saturation
Event-triggered control subject to actuator saturation GEORG A. KIENER Degree project in Automatic Control Master's thesis Stockholm, Sweden 212 XR-EE-RT 212:9 Diploma Thesis Event-triggered control subject
More informationLecture 4 Continuous time linear quadratic regulator
EE363 Winter 2008-09 Lecture 4 Continuous time linear quadratic regulator continuous-time LQR problem dynamic programming solution Hamiltonian system and two point boundary value problem infinite horizon
More informationLinear 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 informationECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77
1/77 ECEN 605 LINEAR SYSTEMS Lecture 7 Solution of State Equations Solution of State Space Equations Recall from the previous Lecture note, for a system: ẋ(t) = A x(t) + B u(t) y(t) = C x(t) + D u(t),
More informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
More information9 The LQR Problem Revisited
9 he LQR Problem Revisited Problem: Compute a state feedback controller u(t) = Kx(t) that stabilizes the closed loop system and minimizes J := lim E [ z(t) z(t) ] t for the LI system Assumptions: a) D
More informationIN many practical systems, there is such a kind of systems
L 1 -induced Performance Analysis and Sparse Controller Synthesis for Interval Positive Systems Xiaoming Chen, James Lam, Ping Li, and Zhan Shu Abstract This paper is concerned with the design of L 1 -
More information8 A First Glimpse on Design with LMIs
8 A First Glimpse on Design with LMIs 8.1 Conceptual Design Problem Given a linear time invariant system design a linear time invariant controller or filter so as to guarantee some closed loop indices
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 informationLecture Note 5: Semidefinite Programming for Stability Analysis
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 5: Semidefinite Programming for Stability Analysis Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio State
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 6: Robust stability and performance in MIMO systems [Ch.8] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 6: Robust Stability and Performance () FEL3210
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More information4. Convex optimization problems
Convex Optimization Boyd & Vandenberghe 4. Convex optimization problems optimization problem in standard form convex optimization problems quasiconvex optimization linear optimization quadratic optimization
More informationLinear Algebra. P R E R E Q U I S I T E S A S S E S S M E N T Ahmad F. Taha August 24, 2015
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS P R E R E Q U I S I T E S A S S E S S M E N T Ahmad F. Taha August 24, 2015 The objective of this exercise is to assess
More informationSYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev
SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1 P. V. Pakshin, S. G. Soloviev Nizhny Novgorod State Technical University at Arzamas, 19, Kalinina ul., Arzamas, 607227,
More informationRiccati Equations and Inequalities in Robust Control
Riccati Equations and Inequalities in Robust Control Lianhao Yin Gabriel Ingesson Martin Karlsson Optimal Control LP4 2014 June 10, 2014 Lianhao Yin Gabriel Ingesson Martin Karlsson (LTH) H control problem
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 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 information1 Some Facts on Symmetric Matrices
1 Some Facts on Symmetric Matrices Definition: Matrix A is symmetric if A = A T. Theorem: Any symmetric matrix 1) has only real eigenvalues; 2) is always iagonalizable; 3) has orthogonal eigenvectors.
More informationModel 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 informationSystem Identification by Nuclear Norm Minimization
Dept. of Information Engineering University of Pisa (Italy) System Identification by Nuclear Norm Minimization eng. Sergio Grammatico grammatico.sergio@gmail.com Class of Identification of Uncertain Systems
More informationRobust Multivariable Control
Lecture 1 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Addresses email: anders.helmerson@liu.se mobile: 0734278419 http://users.isy.liu.se/rt/andersh/teaching/robkurs.html
More informationCDS 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 informationSYSTEMTEORI - 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 information3 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 informationGeometric Control of Patterned Linear Systems
Geometric Control of Patterned Linear Systems Sarah C Hamilton and Mireille E Broucke Abstract We introduce and study a new class of linear control systems called patterned systems Mathematically, this
More informationExam. 135 minutes, 15 minutes reading time
Exam August 6, 208 Control Systems II (5-0590-00) Dr. Jacopo Tani Exam Exam Duration: 35 minutes, 5 minutes reading time Number of Problems: 35 Number of Points: 47 Permitted aids: 0 pages (5 sheets) A4.
More informationLecture 10: Linear Matrix Inequalities Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Matrix Inequalities A linear matrix inequality (LMI)
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 informationI. D. Landau, A. Karimi: A Course on Adaptive Control Adaptive Control. Part 9: Adaptive Control with Multiple Models and Switching
I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 1 Adaptive Control Part 9: Adaptive Control with Multiple Models and Switching I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 2 Outline
More informationModel Uncertainty and Robust Stability for Multivariable Systems
Model Uncertainty and Robust Stability for Multivariable Systems ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Devron Profile Control Solutions Outline Representing model uncertainty.
More informationMarcus Pantoja da Silva 1 and Celso Pascoli Bottura 2. Abstract: Nonlinear systems with time-varying uncertainties
A NEW PROPOSAL FOR H NORM CHARACTERIZATION AND THE OPTIMAL H CONTROL OF NONLINEAR SSTEMS WITH TIME-VARING UNCERTAINTIES WITH KNOWN NORM BOUND AND EXOGENOUS DISTURBANCES Marcus Pantoja da Silva 1 and Celso
More informationRobust 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 informationRecent robust analysis and design results. for simple adaptive control
Recent robust analysis and design results for simple adaptive control Dimitri PEAUCELLE LAAS-CNRS - Université de Toulouse - FRANCE Cooperation program between CNRS, RAS and RFBR A. Fradkov, B. Andrievsky,
More informationBalanced realization and model order reduction for nonlinear systems based on singular value analysis
Balanced realization and model order reduction for nonlinear systems based on singular value analysis Kenji Fujimoto a, and Jacquelien M. A. Scherpen b a Department of Mechanical Science and Engineering
More informationControl 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 informationRelation between Eigenvalues and Singular Values in the Problem of Stability Maintenance of Ellipsoidal Estimates
Relation between Eigenvalues and Singular Values in the Problem of Stability Maintenance of Ellipsoidal Estimates Taalaibek A. Akunov and Anatoly V. Ushakov Department of Automatic and Remote Control St.-Petersburg
More informationProblem 1 Cost of an Infinite Horizon LQR
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K # 5 Ahmad F. Taha October 12, 215 Homework Instructions: 1. Type your solutions in the LATEX homework
More informationOutput Stabilization of Time-Varying Input Delay System using Interval Observer Technique
Output Stabilization of Time-Varying Input Delay System using Interval Observer Technique Andrey Polyakov a, Denis Efimov a, Wilfrid Perruquetti a,b and Jean-Pierre Richard a,b a - NON-A, INRIA Lille Nord-Europe
More informationRobust Control of Time-delay Systems
Robust Control of Time-delay Systems Qing-Chang Zhong Distinguished Lecturer, IEEE Power Electronics Society Max McGraw Endowed Chair Professor in Energy and Power Engineering Dept. of Electrical and Computer
More informationInput-output finite-time stabilization for a class of hybrid systems
Input-output finite-time stabilization for a class of hybrid systems Francesco Amato 1 Gianmaria De Tommasi 2 1 Università degli Studi Magna Græcia di Catanzaro, Catanzaro, Italy, 2 Università degli Studi
More informationSolving non-linear systems of equations
Solving non-linear systems of equations Felix Kubler 1 1 DBF, University of Zurich and Swiss Finance Institute October 7, 2017 Felix Kubler Comp.Econ. Gerzensee, Ch2 October 7, 2017 1 / 38 The problem
More information