Robust Multivariable Control
|
|
- Ann Cynthia Stewart
- 5 years ago
- Views:
Transcription
1 Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet
2 Today s topics
3 Today s topics Norms
4 Today s topics Norms Representation of dynamic systems
5 Today s topics Norms Representation of dynamic systems Lyapunov equations
6 Today s topics Norms Representation of dynamic systems Lyapunov equations Gramians
7 Today s topics Norms Representation of dynamic systems Lyapunov equations Gramians Balancing
8 Today s topics Norms Representation of dynamic systems Lyapunov equations Gramians Balancing Model reduction
9 Today s topics Norms Representation of dynamic systems Lyapunov equations Gramians Balancing Model reduction H 2 and H norms
10 Small gain theorem G If G is stable and G < 1 then the closed loop system is stable if 1. Which norm should we use?
11 Norms p-norms for vectors: x p = ( ) 1/p x i p i Most common: p = 1: sum of the absolute values of the elements; p = 2: Euclidian norm of the vector x; p = : the maximum of absolute values. Induced norms maximum singular value. A = sup x =0 Ax 2 x 2 = σ(a)
12 Singular values Singular value decomposition where A = UΣV, σ 1 σ 2 Σ = U U = I, V V = I... Rn n σn assuming σ 1 σ 2... σ n 0. (A A)V = (UΣV ) UΣV V = VΣU UΣ = VΣ 2, which means that Σ 2 are the eigenvalues of A A. Numerically, this is not a good method to compute Σ.
13 Singular values Also, [ 0 ] [ A U U A 0 V V [ UΣ UΣ = VΣ VΣ ] [ 0 UΣV = ] [ U U VΣU 0 V V ] [ ] [ ] U U Σ 0 = V V 0 Σ ] and, consequently, ±Σ are the eigenvalues of [ 0 ] A A 0
14 Singular values and LMIs We can also use a Linear Matrix Inequality (LMI) to find the maximum singluar value. Find the minimum γ for which [ ] γi A A 0. γi Here 0 means positive semidefinite (all eigenvalues 0). Also, min{γ : γi A A 0} = σ 2 (A) γ
15 Representation of dynamic systems State-space representation ẋ = Ax + Bu y = Cx + Du Transfer function G(s) = C(sI A) 1 B + D (1) Sometimes we write this as [ ] [ ẋ A B = y C D or [ ] A B G = C D ] [ x u ] (2) (3)
16 Serial and parallel connection Serial connection: Parallel connection: [ ] A1 B G 1 = 1 C 1 D 1 G 1 G 2 = G 1 + G 2 = [ ] A2 B G 2 = 2 C 2 D 2 A 1 B 1 C 2 B 1 D 2 0 A 2 B 2 C 1 D 1 C 2 D 1 D 2 A 1 0 B 1 0 A 2 B 2 C 1 C 2 D 1 + D 2 (4) (5) (6)
17 Lyapunov equaions The linear equation A T X + XA + Q = 0 (7) is a Lyapunov equation. It has a unique solution if A och A do not have any common eigenvalues; for instance if A has all its eigenvalues in the left hand plane. In Matlab we can solve the equation using X = lyap (A, Q);
18 Observability gramian L o The observability gramian is defined by Regard the following system (u = 0): {ẋ = Ax A T L o + L o A + C T C = 0. (8) y = Cx Introduce V o = x T L o x and take its time derivative V o = ẋ T L o x + x T L o ẋ = x T A T L o x + x T L o Ax = x T (A T L o + L o A) x = x T C T Cx = y T y = y 2 2 }{{} C T C
19 Observability gramian L o Integrate V o T 0 T V o dt = V o (x(t)) V o (x(0)) = y(t) 2 2 dt. 0 If we assume that the system is stable and T then y(t) approaches zero. Thus V o (x(0)) = x T (0)L o x(0) = T 0 y(t) 2 2 dt The gramian, L o, includes information about how much a certain state becomes visible (as energy) in the output signal.
20 Controllabilty gramian L c The controllability gramian is defined by AL c + L c A T + BB T = 0. (9) Introduce V c (x) = x T Lc 1 x. This is a measure of the smallest energy needed in the input to obtain a certain state x. As before, take the derivative of V c with ẋ = Ax + Bu: V c (x) = ẋ T Lc 1 x + x T Lc 1 ẋ = (Ax + Bu) T Lc 1 x + x T Lc 1 (Ax + Bu) = u T B T Lc 1 x + x T Lc 1 Bu + x T (A T Lc 1 + Lc 1 A) x }{{} = (B T L 1 c Lc 1 BB T Lc 1 x + u) T (B T L 1 x + u) + u T u u T u c
21 Controllability gramian L c The control law u opt = B T Lc 1 x gives the smallest control energy. The energy needed for reaching a certain state (starting in x( ) = 0) is given by u 2 = V c (x) + u + B T L 1 c x 2 (10) Thus, u 2 V c (x) = x T L 1 c x, with equality when u = u opt.
22 Balancing The representation of a system is not unique. [ ] A B G = C D Using a similarity transformation, it can also be written as [ ] [ Â ˆB TAT 1 ] TB G = = Ĉ ˆD CT 1 D (11) (12) The transformation also affects the gramians ˆL o = T T L o T 1 (13) and ˆL c = TL c T T (14)
23 Balancing Choose T in a clever way such that σ 1 σ 2 ˆL o = ˆL c = Σ =... σn (15) Then σ 2 i = λ i (ˆL o ˆL c ) = λ i (L o L c ). (16) To obtain this factor L o och L c, for instance using Cholesky decomposition: L o = Uo T U o and L c = Uc T U c. Then compute UΣV T = U o Uc T using singular value decomposition, where UU T = VV T = I. The transformation matrix is given by T = Σ 2 1 U T U o.
24 An example of balancing G(s) = (s + 1)(s + 5) = = with Σ = diag [σ 1, σ 2 ] = diag [0.1124, ] In Matlab this is done by g = tf(1,[1 1])*tf(1,[1 5]); [gb,sig] = balreal(g); Here σ 2 is small compared to σ 1. If we remove the corresponding state, x 2, the H error is limited by 2σ 2 : G Ĝ 1 2σ 2, Ĝ 1 (s) = [ ] (17)
25 An example of balancing Singular Values G G 1 G G 1 2σ 2 25 Singular Values (db) Frequency (rad/sec)
26 Model reduction Model reduction can be done by truncating the states that have the smallest Hankel singular values: G n Ĝ r 2 n i=r+1 σ i (18) This upper bound can be halved by using a different scheme that also modifes the D matrix: σ r+1 G n Ĝ r n σ i (19) i=r+1
27 With modified D optimal Hankel Singular Values G G 1 G G 1 =σ Singular Values (db) Frequency (rad/sec)
28 H 2 norm e G z Let e = δ (dirac pulse) or white noise Energy in the output: ẋ = Ax + Bδ(t 0) }{{} x=b V(x) = x T (0)L o x(0) = B T L o B = G 2 2 = CL c C T
29 H 2 norm e G z In MIMO case G 2 2 = tr B T L o B = tr CL c C T The H 2 norm can be computed by solving a linear matrix equation.
30 H norm H norm is a measure of the maximum gain of a stable system, G, over all frequencies, ω. That is to say L norm also applies to unstable systems. G = max G(jω) (20) ω H norm cannot be computed direct (as the H 2 norm), but we can test if G < γ.
31 H norm [ A B G(s) = C D ] = D + C(sI A) 1 B Assume that G < γ. If G is SISO then Φ(jω) = γ 2 σ 2 (G(jω)) = γ 2 G(jω) G(jω) > 0, for all ω. If γ G then Φ(jω) has at least one zero for ω R.
32 H norm where This gives Φ(s) = γ 2 G (s)g(s) = γ 2 G( s) T G(s) G (s) = [ A T C T B T D T A 0 B Φ(s) = C T C A T C T D D T C B T γ 2 I D T D Requirement 1: γ > σ(d) = λ max (D T D). ] If Φ(s) has a zero on the imaginary axis then Φ 1 (s) must have a pole on the imaginary axis.
33 Inverse of systems We need to compute the inverse of a state space system. The inverse of a system can be written as [ ] 1 [ G 1 A B A BD 1 C BD 1 (s) = = C D D 1 C D 1 ] since...
34 Inverse of systems G 1 (s)g(s) = I I 0 = 0 I I A BD 1 C BD 1 C BD 1 D 0 A B D 1 C D 1 C D 1 D A BD 1 C BD 1 C B 0 A B D 1 C D 1 C I A BD 1 C A BD 1 C 0 0 A B = D 1 C D 1 C I A BD 1 C 0 0 = 0 A B I D 1 C 0 I I I 0 0 I I I I 0 0 I I
35 H norm The inverse of the system can be written as ] 1 [ A BD 1 C BD 1 = C [ A B ] D D 1 C D 1 Let H be the A matrix of the inverse of Φ(s): [ A 0 H = C T C A T [ = ] [ A + BR 1 D T C C T C + C T DR 1 D T C B C T D ] ( γ 2 I D D) T 1 [ D T C B T ] } {{ } R BR 1 B T A T C T DR 1 B T If H has eigenvalues on the imaginary axis then γ G. ]
36 Exemple of the H norm Let G(s) = s+s 2 : γ eig H 10 ±1j, ±0.9950j 10.1 ± ± j ±ɛ ± j 1/(s 2 +s/10+1) Singular Values (db) Frequency (rad/sec)
37 Induced norms The H norm can also be described as an induced L 2 or l 2 norm. For a discrete-time system, if γ > G then max u T t=0 ( ) y T (t)y(t) γ 2 u T (t)u(t) < 0 for the system [ x(t + 1) y(t) ] [ A B = C D ] [ x(t) u(t) ].
38 Induced norms Introduce a Lyapunov function, V(x) = x T Px, where P = P T 0 is positive definite: = = = T ( ) y T (t)y(t) γ 2 u T (t)u(t) t=0 T t=0 T t=0 T t=0 ( ) y T (t)y(t) γ 2 u T (t)u(t) + V(x(T + 1)) V(x(0)) + V(x(0)) V(x(T + 1)) ( ) y T (t)y(t) γ 2 u T (t)u(t) + V(x(t + 1)) V(x(t)) + V(x(0)) V(x(T + 1)) [ x(t) u(t) ] T ([ A T PA P A T PB B T PA B T PB γ 2 I + V(x(0)) V(x(T + 1)) ] [ C T + D T ] [ C D ] ) [ x(t) u(t) ] T
39 LMIs Assure that [ A T PA P A T PB B T PA B T PB γ 2 I ] [ C T + D T ] [ C D ] 0 by choosing an appropriate matrix P = P T 0. The smallest γ that we can find such a P for gives the upper bound of G. This is a linear matrix inequality (LMI), which is a convex problem.
40 LMIs For continuous-time systems [ PA + A T P PB B T P γ 2 I ] [ C T + D T ] [ C D ] 0 The smallest γ for which we can find such a P = P T 0 gives the upper bound of G.
Model reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More 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 informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
More 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 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 informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationCDS 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 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 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 informationProblem set 5 solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION Problem set 5 solutions Problem 5. For each of the stetements below, state
More 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 informationTopic # Feedback Control Systems
Topic #20 16.31 Feedback Control Systems Closed-loop system analysis Bounded Gain Theorem Robust Stability Fall 2007 16.31 20 1 SISO Performance Objectives Basic setup: d i d o r u y G c (s) G(s) n control
More 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 informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
More informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 2 Homework 7: State Estimators (a) For the same system as discussed in previous slides, design another closed-loop state estimator,
More informationModel Reduction for Linear Dynamical Systems
Summer School on Numerical Linear Algebra for Dynamical and High-Dimensional Problems Trogir, October 10 15, 2011 Model Reduction for Linear Dynamical Systems Peter Benner Max Planck Institute for Dynamics
More informationIterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations
Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Klajdi Sinani Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University
More 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 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 informationH 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
H 2 Optimal State Feedback Control Synthesis Raktim Bhattacharya Aerospace Engineering, Texas A&M University Motivation Motivation w(t) u(t) G K y(t) z(t) w(t) are exogenous signals reference, process
More 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 informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More 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 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 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 informationModern 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 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 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 informationECEEN 5448 Fall 2011 Homework #4 Solutions
ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with
More informationControllability, Observability, Full State Feedback, Observer Based Control
Multivariable Control Lecture 4 Controllability, Observability, Full State Feedback, Observer Based Control John T. Wen September 13, 24 Ref: 3.2-3.4 of Text Controllability ẋ = Ax + Bu; x() = x. At time
More 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 informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of
More 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 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 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 informationBalanced Model Reduction
1 Balanced Realization Balanced Model Reduction CEE 629. System Identification Duke University, Fall 17 A balanced realization is a realization for which the controllability gramian Q and the observability
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 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 informationLecture 8. Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control. Eugenio Schuster.
Lecture 8 Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture
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 informationControl Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli
Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More informationLecture 4: Analysis of MIMO Systems
Lecture 4: Analysis of MIMO Systems Norms The concept of norm will be extremely useful for evaluating signals and systems quantitatively during this course In the following, we will present vector norms
More informationOrder Reduction for Large Scale Finite Element Models: a Systems Perspective
Order Reduction for Large Scale Finite Element Models: a Systems Perspective William Gressick, John T. Wen, Jacob Fish ABSTRACT Large scale finite element models are routinely used in design and optimization
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 informationOrder Reduction of a Distributed Parameter PEM Fuel Cell Model
Order Reduction of a Distributed Parameter PEM Fuel Cell Model María Sarmiento Carnevali 1 Carles Batlle Arnau 2 Maria Serra Prat 2,3 Immaculada Massana Hugas 2 (1) Intelligent Energy Limited, (2) Universitat
More informationTopic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback
Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall
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 informationMAE 143B - Homework 8 Solutions
MAE 43B - Homework 8 Solutions P6.4 b) With this system, the root locus simply starts at the pole and ends at the zero. Sketches by hand and matlab are in Figure. In matlab, use zpk to build the system
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 informationH 2 optimal model reduction - Wilson s conditions for the cross-gramian
H 2 optimal model reduction - Wilson s conditions for the cross-gramian Ha Binh Minh a, Carles Batlle b a School of Applied Mathematics and Informatics, Hanoi University of Science and Technology, Dai
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 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 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 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 informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationGramians of structured systems and an error bound for structure-preserving model reduction
Gramians of structured systems and an error bound for structure-preserving model reduction D.C. Sorensen and A.C. Antoulas e-mail:{sorensen@rice.edu, aca@rice.edu} September 3, 24 Abstract In this paper
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 informationModel reduction of coupled systems
Model reduction of coupled systems Tatjana Stykel Technische Universität Berlin ( joint work with Timo Reis, TU Kaiserslautern ) Model order reduction, coupled problems and optimization Lorentz Center,
More informationA Singular Value Decomposition Based Closed Loop Stability Preserving Controller Reduction Method
A Singular Value Decomposition Based Closed Loop Stability Preserving Controller Reduction Method Kin Cheong Sou and Anders Rantzer Abstract In this paper a controller reduction method which preserves
More information6 Optimal Model Order Reduction in the Hankel Norm
6 Optimal Model Order Reduction in the Hankel Norm In this lecture, we discuss optimal model reduction in the Hankel norm. There is no known solution to the optimal model reduction problem in the H -norm.
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 informationMoving Frontiers in Model Reduction Using Numerical Linear Algebra
Using Numerical Linear Algebra Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory Group Magdeburg, Germany Technische Universität Chemnitz
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 informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 5. Input-Output Stability DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Input-Output Stability y = Hu H denotes
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 informationHuman Friendly Control : an application to Drive by Wire
Human Friendly Control : an application to Drive by Wire By Emmanuel Witrant Master Thesis, 2001-2002 Laboratoire d Automatique de Grenoble with Carlos Canudas-de-Wit 2 Juillet 2002 Typeset by FoilTEX
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 16.323 Principles of Optimal Control Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 16.323 Lecture
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 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 informationTopic # Feedback Control
Topic #5 6.3 Feedback Control State-Space Systems Full-state Feedback Control How do we change the poles of the state-space system? Or,evenifwecanchangethepolelocations. Where do we put the poles? Linear
More informationMODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS
MODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS Ulrike Baur joint work with Peter Benner Mathematics in Industry and Technology Faculty of Mathematics Chemnitz University of Technology
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 informationEL2520 Control Theory and Practice
So far EL2520 Control Theory and Practice r Fr wu u G w z n Lecture 5: Multivariable systems -Fy Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden SISO control revisited: Signal
More informationSiMpLIfy: A Toolbox for Structured Model Reduction
SiMpLIfy: A Toolbox for Structured Model Reduction Martin Biel, Farhad Farokhi, and Henrik Sandberg ACCESS Linnaeus Center, KTH Royal Institute of Technology Department of Electrical and Electronic Engineering,
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 informationLecture 8. Applications
Lecture 8. Applications Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering May 8, 205 / 3 Logistics hw7 due this Wed, May 20 do an easy problem or CYOA hw8 (design problem) will
More informationNonlinear Control. Nonlinear Control Lecture # 18 Stability of Feedback Systems
Nonlinear Control Lecture # 18 Stability of Feedback Systems Absolute Stability + r = 0 u y G(s) ψ( ) Definition 7.1 The system is absolutely stable if the origin is globally uniformly asymptotically stable
More informationProfessor Fearing EE C128 / ME C134 Problem Set 7 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley
Professor Fearing EE C8 / ME C34 Problem Set 7 Solution Fall Jansen Sheng and Wenjie Chen, UC Berkeley. 35 pts Lag compensation. For open loop plant Gs ss+5s+8 a Find compensator gain Ds k such that the
More informationSingular Value Decomposition (SVD)
School of Computing National University of Singapore CS CS524 Theoretical Foundations of Multimedia More Linear Algebra Singular Value Decomposition (SVD) The highpoint of linear algebra Gilbert Strang
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 informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationarxiv: v1 [math.oc] 17 Oct 2014
SiMpLIfy: A Toolbox for Structured Model Reduction Martin Biel, Farhad Farokhi, and Henrik Sandberg arxiv:1414613v1 [mathoc] 17 Oct 214 Abstract In this paper, we present a toolbox for structured model
More informationDECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION
Journal of ELECTRICAL ENGINEERING, VOL. 58, NO. 6, 2007, 307 312 DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Szabolcs Dorák Danica Rosinová Decentralized control design approach based on partial
More informationQuadratic Stability of Dynamical Systems. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
.. Quadratic Stability of Dynamical Systems Raktim Bhattacharya Aerospace Engineering, Texas A&M University Quadratic Lyapunov Functions Quadratic Stability Dynamical system is quadratically stable if
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 informationPOLE PLACEMENT. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 19
POLE PLACEMENT Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 19 Outline 1 State Feedback 2 Observer 3 Observer Feedback 4 Reduced Order
More informationWeighted balanced realization and model reduction for nonlinear systems
Weighted balanced realization and model reduction for nonlinear systems Daisuke Tsubakino and Kenji Fujimoto Abstract In this paper a weighted balanced realization and model reduction for nonlinear systems
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 informationControl Systems. Frequency domain analysis. L. Lanari
Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic
More informationEECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 5-8 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3
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 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 informationModel Reduction of Inhomogeneous Initial Conditions
Model Reduction of Inhomogeneous Initial Conditions Caleb Magruder Abstract Our goal is to develop model reduction processes for linear dynamical systems with non-zero initial conditions. Standard model
More informationModel Reduction using a Frequency-Limited H 2 -Cost
Technical report from Automatic Control at Linköpings universitet Model Reduction using a Frequency-Limited H 2 -Cost Daniel Petersson, Johan Löfberg Division of Automatic Control E-mail: petersson@isy.liu.se,
More informationFluid flow dynamical model approximation and control
Fluid flow dynamical model approximation and control... a case-study on an open cavity flow C. Poussot-Vassal & D. Sipp Journée conjointe GT Contrôle de Décollement & GT MOSAR Frequency response of an
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
More informationRobust Control 5 Nominal Controller Design Continued
Robust Control 5 Nominal Controller Design Continued Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 4/14/2003 Outline he LQR Problem A Generalization to LQR Min-Max
More information5. Observer-based Controller Design
EE635 - Control System Theory 5. Observer-based Controller Design Jitkomut Songsiri state feedback pole-placement design regulation and tracking state observer feedback observer design LQR and LQG 5-1
More 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 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 information