Linear Matrix Inequality (LMI)
|
|
- Nathaniel McKinney
- 5 years ago
- Views:
Transcription
1 Linear Matrix Inequality (LMI) A linear matrix inequality is an expression of the form where F (x) F 0 + x 1 F x m F m > 0 (1) x = (x 1,, x m ) R m, F 0,, F m are real symmetric matrices, and the inequality > 0 in (1) means positive definite, i.e., u T F (x)u > 0 for all u R n, u 0. Equivalently, the smallest eigenvalue of F (x) is positive. 1
2 Definition[Linear matrix inequality(lmi)] A linear matrix inequality is F (x) > 0 (2) where F is an affine function mapping a finite dimensional vector space to the set S n {M : M = M T R n n }, n > 0, of real matrices. remark Recall, from definition, that an affine mapping F : V S n necessarily takes the form F (x) = F 0 +T (x) where F 0 S n and T : V S n is a linear transformation. Thus if V is of dimension m, and {e 1,, e m } constitutes a basis for V, then we can write T (x) = m j=1 x j F j where the elements {x 1,, x m } are such that x = m j=1 x j e j and F j = T (e j ) for j = 1,, m. Hence we obtain (1) as a special case. 2
3 Remark. The same remark applies to mappings F : R m 1 m 2 S n where m 1, m 2 Z +. A simple example where m 1 = m 2 is the Lyapunov inequality F (X) = A T X + XA + Q > 0. Here, A, Q R m m are assumed to be given, Q is symmetric, and X R m m is the unknown matrix. In this case, the domain V of F in definition is equal to S m. We can view this LMI as a special case of (1) by defining a basis E 1,, E m of S m and writing X = m j=1 x j E j : F (X) = F m j=1 x j E j = F 0 + which is of the form (1). = F 0 + m j=1 m j=1 x j F (E j ) x j F j 3
4 Remark. The LMI F (x) = F 0 + xf x m F m defines a convex constraint on x = (x 1,, x m ). i.e., the set F {x : F (x) > 0} is convex. Indeed, if x 1, x 2 F and α (0, 1) then F (αx 1 +(1 α)x 2 ) = αf (x 1 )+(1 α)f (x 2 ) > 0 Convexity has an important consequence: even though the LMI has no analytical solution in general, it can be solved numerically with guarantees of fining a solution when one exists. Although the LMI may seem special, it turns out that many convex sets can be represented in this way. 4
5 1. Note that a system of LMIs (i.e. a finite set of LMIs) can be written as a single LMI since F 1 (x) < 0. F K (x) < 0 is equivalent to F (x) diag[f 1 (x),, F K (x)] < 0 2. Combined constraints (in the unknown x) of the form { F (x) > 0 Ax = b or { F (x) > 0 x = Ay + b for some y where the affine function F : R m S n and matrices A R n m and b R n are given can be lumped into one LMI. More generally, the combined equations { F (x) > 0 (3) x M 5
6 where M is an affine subset of R n, i.e. M = x 0 + M 0 = {x 0 + m m M 0 } with x 0 R n and M 0 a linear subspace of R n, can be written in the form of one single LMI. In order to see this, let e 1,, e k R n be a basis of M 0 and let F (x) = F 0 + T (x) be decomposed as in remark. Then (3) can be rewritten as 0 < F (x) = F 0 + T (x 0 + k j=1 = F 0 + T (x 0 ) + }{{} constant part x j e j ) k j=1 = F 0 + x 1 F x k F k F ( x) x j T (e j ) }{{} linear part where F 0 = F 0 + T (x 0 ), F j = T (e j ) and x = (x 1,, x k ). This implies that x R n satisfies (3) if and only if F (x) > 0. Note
7 that the dimension of x is smaller than the dimension of x. 3. (Schur Complement) Let F : V S n be an affine function partitioned to [ ] F11 (x) F F (x) = 12 (x) F 21 (x) F 22 (x) where F 11 (x) is square. Then F (x) > 0 iff { F11 (x) > 0 F 22 (x) F 21 (x)f11 1 (x)f 12(x) > 0 (4) Note that the second inequality in (4) is a nonlinear matrix inequality in x. It follows that nonlinear matrix inequalities of the form (4) can be converted to LMIs, and nonlinear inequalities (4) define a convex constraint on x.
8 Types of LMI problems Suppose that F, G : V S n 1 and H : V S n 2 are affine functions. Feasibility: The test whether or not there exist solutions x of F (x) > 0 is called a feasibility problem. The LMI is called nonfeasible if no solutions exist. Optimization: Let f : S R and suppose that S = {x F (x) > 0}. The problem to determine V opt = inf x S f(x) is called an optimization problem with an LMI constraint. Generalized eigenvalue problem: Minimize a scalar λ R subject to λf (x) G(x) > 0 F (x) > 0 H(x) > 0 6
9 What are LMIs good for? Many optimization problems in control design, identification, and signal processing can be formulated using LMIs. Example. Asymptotic stability of the LTI system ẋ = Ax, A R n n (5) Lyapunov said, asymptotically stable iff there exists X S n such that X > 0, A T X + XA < 0 i.e. equivalent to feasibility of the LMI [ ] X 0 0 A T > 0 X XA 7
10 Example. Determine a diagonal matrix D such that DMD 1 < 1 where M is some given matrix. Since DMD 1 < 1 D T M T D T DMD 1 < I M T D T DM < D T D X M T XM > 0 where X := D T D > 0 we see that the existence of such a matrix means the feasibility of LMI. 8
11 Example. Let F be an affine function and consider the problem of minimizing f(x) σ max (F (x)) over x. λ max (F T (x)f (x)) < γ γi F T (x)f (x) > 0 [ if we define [ ] x x, F ( x) γ γi F (x) [ γi F (x) F T (x) I ] F T (x) I > 0 ], f( x) γ, then F is an affine function of x and the problem to minimize the maximum eigenvalue of F (x) is equivalent to determining inf f( x) subject to the LMI F ( x) > 0. Hence, this is an optimization problem with a linear objective function f and an LMI constraint. 9
12 Example(Simultaneous stabilization) Consider k LTI systems with n-dim state space and m-dim input space: ẋ = A i x + B i u where A i R n n and B i R n m, i 1,, k. We d like to find a state feedback law u = F x, F R m n such that the eigenvalues λ(a i + B i F ) lie on the LHP for i 1,, k. From the example above, this is solved when we find matrices F and X i, i 1,, k such that for i 1,, k, { Xi > 0 (A i + B i F ) T (6) X i + X i (A i + B i F ) < 0 Note that this is not a system of LMIs in X i and F. If we introduce Y i = X 1 i and K = F Y i, then (6) becomes { Yi > 0 A i Y i + Y i A T i + B ik + K T i B i < 0, which can be further simplified by assuming 10
13 the existence of a joint Lyapunov function, i.e. X i = = X k = X. The joint stabilization problem has a solution if this system of LMIs is feasible.
14 H nominal performance Consider x = Ax + Bu (7) y = Cx + Du (8) with state space X = R n, input space U = R m and output space Y = R p. proposition If the system (7) is asymptotically stable then G < γ whenever there exists a solution K = K T > 0 to the LMI [ A T K + KA + C T C KB + C T ] D B T K + D T C D T D γ 2 < 0. (9) I Can compute the H norm of the transfer function by minimizing γ > 0 over all variables γ and K > 0 that satisfy the LMI. 11
15 H 2 nominal performance We take impulsive inputs of the form u(t) = δ(t)e i with e i the i th basis vector in the standard basis of the input space R m. (i = 1 m). With zero initial conditions, the corresponding output y i L 2 and is given by y i (t) = C exp(at)be i for t > 0 De i δ(t) for t = 0 0 for t < 0.. Only if D = 0, the sum of the squared norms of all such impulse responses m i=1 y i 2 2 is well defined and given by m i=1 y i 2 2 = trace = trace 0 BT exp(a t )C T C exp(at)b dt 0 C exp(at)bbt exp(a T t)c T dt = trace G(jω)G (jω) dω where G is the transfer function of the system. 12
16 proposition Suppose that the system (7) is asymptotically stable (and D = 0), then the following statements are equivalent. (a) G 2 < γ (b) there exists K = K T > 0 and Z such that [ A T ] [ ] K + KA KB K C T B T < 0; > 0; K I C Z (10) trace(z) < γ 2 (11) (c) there exists K = K T > 0 and Z such that [ AK + KA T KC T ] [ ] K B < 0; CK I B T > 0; Z (12) trace(z) < γ 2 (13) 13
17 pf. note that G 2 < γ is equivalent to requiring that the controllability gramian W c := 0 exp(at)bb T exp(a T t) dt satisfies trace(cw C T ) < γ 2. Since the controllability gramian is the unique positive definite solution to the Lyapunov equation AW + W A T + BB T = 0 this is equivalent to saying that there exists X > 0 such that AX + XA T + BB T < 0; trace(cxc T ) < γ 2. With a change of variables K := X 1, this is equivalent to the existence of K > 0 and Z such that A T K + KA + KBB T K < 0; CK 1 C T < Z; and trace(z) < γ 2. 14
18 Now, using Schur complements for the first two inequalities yields that G 2 < γ is equivalent to the existence of K > 0 and Z such that [ A T ] [ ] K + KA KB K C T B T < 0; > 0; K I C Z and trace(z) < γ 2. The equivalence with (12) is obtained by the observation that G 2 = G T 2. Therefore, the smallest possible upper bound of the H2-norm of the transfer function can be calculated by minimizing the criterion trace(z) over the variables K > 0 and Z that satisfy the LMIs defined by the first two inequalities in (10) or (12).
19 Controller Synthesis Let and ẋ = Ax + B 1 w + B 2 u z = C x + D 1 w + D 2 u z 2 = C 2 x + D 21 w + D 22 u y = C y x + D y1 w ẋ K = A K x K + B K y u = C K x K + D K y be state-space realizations of the plant P (s) and the controller K(s) respectively. 15
20 Denoting by T (s) and T 2 (s) the CL TF from w to z and z 2, respectively, we consider the following multi-objective synthesis problem: Design an output feedback controller u = K(s)y such that H performance: maintains the H norm of T below γ 0. H 2 performance: of T 2 below ν 0. maintains the H 2 norm Multi-objective H 2 /H controller design: minimizes the trade-off criterion of the form α T 2 + β T with some α, β 0. Pole placement: places the CL poles in some prescribed LMI region D. 16
21 Let the following denote the corresponding CL state-space eqns, x cl = A cl x cl + B cl w z = C cl1 x cl + D cl1 w z 2 = C cl2 x cl + D cl2 w then our design objectives can be expressed as follows: H performance: the CL RMS gain from w to z does not exceed γ iff there exists a symmetric matrix X such that A cl X + X A T cl B cl X Ccl1 T Bcl T I Dcl1 T C cl1 X D cl1 γ 2 I < 0 X > 0 H 2 performance: the LQG cost from w to z 2 does not exceed ν iff D cl2 = 0 and there 17
22 exists a symmetric matrices X 2 and Q such that [ Acl X 2 + X 2 A T ] cl B cl < 0 I [ B T cl Q Ccl2 T X 2 X 2 Ccl2 T > 0 X 2 trace(q) < ν 2 ] Pole placement: the CL poles lie in the LMI region D := {z C : L + Mz + M T z < 0} with L = L T = [λ ij ] 1 i,j m and M = [µ ij ] 1 i,j m iff there exists a symmetric matrix X pol such that [λ ij X pol + µ ij A cl X pol + µ ji X pol A T cl ] 1 i,j m < 0 X pol > 0.
23 For tractability, we seek a single Lyapunov matrix X := X = X 2 = X pol that enforces all three sets of constraints. Factorizing X as [ ] [ ] 1 R I 0 S X = M T 0 I N T and introducing the transformed controller variables: B K := NB K + SB 2 D K C K := C K M T + D K C y R A K := NA K M T + NB K C y R + SB 2 C K M T +S(A + B 2 D K C y )R, the inequality constraints on X are turned into LMI constraints in the variables R, S, Q, A K, B K, C K and D K. And we have the following suboptimal LMI formulation of our multi-objective synthesis problem: 18
24 Minimize αγ 2 +βtrace(q) over R, S, Q, A K, B K, C K, D K and γ 2 satisfying: AR + RA T + B 2 C K + CK T BT 2 C R + D 2 C K [ [ R λ ij I A K + A + B 2 D K C y A T S + SA + B K C y + Cy T BT K C + D 2 D K C y I S ] Q [ AR + B2 C + µ K A + ij A K SA µ ji [ (AR + B2 C K (A + B 2 D K C 19
25 Given optimal solutions γ, Q of this LMI problem, the closed loop performances are bounded by T γ, T 2 trace(q ). This has been implemented by the matlab command hinfmix. 20
26 Reference Boyd S, El Ghaoui L, Feron E, Balakrishnan V. Linear matrix inequalities in system and control theory, vol. 15 ed. Scherer C, Weiland S. Linear matrix inequalities in control. Lecture notes of DISC Course LMI Control Toolbox, Gahinet, Nemirovski, Laub, Chilali, Mathworks 21
27 Affine combinations of linear systems Often models uncertainty about specific parameters is reflected as uncertainty in specific entries of the state space matrices A, B, C, D. Let p = (p 1,..., p n ) denote the parameter vector which expresses the uncertain quantities in the system and suppose that this parameter vector belongs to some subset P R n. Then the uncertain model can be thought of as being parameterized by p P through its state space representation ẋ = A(p)x + B(p)u (14) y = C(p)x + D(p)u. (15) One way to think of equations of this sort is to view them as a set of linear time-invariant systems as parameterized by p P. However, if p is time, then (14) defines a linear time varying dynamical system and it can therefore also be viewed as such. If components of p are 22
28 time varying and coincide with state components then (14) is better viewed as a nonlinear system. Of particular interest will be those systems in which the system matrices affinely depend on p. This means that A(p) = A 0 + p 1 A p n A n (16) B(p) = B 0 + p 1 B p n B n (17) C(p) = C 0 + p 1 C p n C n (18) D(p) = D 0 + p 1 D p n D n. (19) Or, written in a more compact form where S(p) = S 0 + p 1 S p n S n S(p) = [ A(p) B(p) C(p) D(p) is the system matrix associated with (14). We call these models affine parameter dependent ]
29 models. In MATLAB such a system is represented with the routines psys and pvec. For n = 2 and a parameter box P {(p 1, p 2 ) p 1 [p min 1, p max 1 ], p 2 [p min 2, p max 2 ]} the syntax is affsys = psys( p, [s0, s1, s2] ); p = pvec( box, [p1min p1max ; p2min p2max]) where p is the parameter vector whose i-th component ranges between pimin and pimax. Bounds on the rate of variations, ṗ i (t) can be specified by adding a third argument rate when calling pvec.
30 See also the following routines: pdsimul for time simulations of affine parameter models aff2pol to convert an affine model to an equivalent polytopic model pvinfo to inquire about the parameter vector 23
Static 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 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 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 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 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 informationConvex Optimization Approach to Dynamic Output Feedback Control for Delay Differential Systems of Neutral Type 1,2
journal of optimization theory and applications: Vol. 127 No. 2 pp. 411 423 November 2005 ( 2005) DOI: 10.1007/s10957-005-6552-7 Convex Optimization Approach to Dynamic Output Feedback Control for Delay
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 Output Feedback Controller Design via Genetic Algorithms and LMIs: The Mixed H 2 /H Problem
Robust Output Feedback Controller Design via Genetic Algorithms and LMIs: The Mixed H 2 /H Problem Gustavo J. Pereira and Humberto X. de Araújo Abstract This paper deals with the mixed H 2/H control problem
More informationIntroduction to linear matrix inequalities Wojciech Paszke
Introduction to linear matrix inequalities Wojciech Paszke Institute of Control and Computation Engineering, University of Zielona Góra, Poland e-mail: W.Paszke@issi.uz.zgora.pl Outline Introduction to
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 informationOn Bounded Real Matrix Inequality Dilation
On Bounded Real Matrix Inequality Dilation Solmaz Sajjadi-Kia and Faryar Jabbari Abstract We discuss a variation of dilated matrix inequalities for the conventional Bounded Real matrix inequality, and
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 informationAn LMI Approach to the Control of a Compact Disc Player. Marco Dettori SC Solutions Inc. Santa Clara, California
An LMI Approach to the Control of a Compact Disc Player Marco Dettori SC Solutions Inc. Santa Clara, California IEEE SCV Control Systems Society Santa Clara University March 15, 2001 Overview of my Ph.D.
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 informationOptimization based robust control
Optimization based robust control Didier Henrion 1,2 Draft of March 27, 2014 Prepared for possible inclusion into The Encyclopedia of Systems and Control edited by John Baillieul and Tariq Samad and published
More informationAn LMI Optimization Approach for Structured Linear Controllers
An LMI Optimization Approach for Structured Linear Controllers Jeongheon Han* and Robert E. Skelton Structural Systems and Control Laboratory Department of Mechanical & Aerospace Engineering University
More informationA new robust delay-dependent stability criterion for a class of uncertain systems with delay
A new robust delay-dependent stability criterion for a class of uncertain systems with delay Fei Hao Long Wang and Tianguang Chu Abstract A new robust delay-dependent stability criterion for a class of
More informationOn Computing the Worst-case Performance of Lur'e Systems with Uncertain Time-invariant Delays
Article On Computing the Worst-case Performance of Lur'e Systems with Uncertain Time-invariant Delays Thapana Nampradit and David Banjerdpongchai* Department of Electrical Engineering, Faculty of Engineering,
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 informationA New Strategy to the Multi-Objective Control of Linear Systems
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 25 Seville, Spain, December 12-15, 25 TuC8.6 A New Strategy to the Multi-Objective Control of Linear
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 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 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 informationAppendix A Solving Linear Matrix Inequality (LMI) Problems
Appendix A Solving Linear Matrix Inequality (LMI) Problems In this section, we present a brief introduction about linear matrix inequalities which have been used extensively to solve the FDI problems described
More informationLinear Matrix Inequalities in Control
Linear Matrix Inequalities in Control Delft Center for Systems and Control (DCSC) Delft University of Technology The Netherlands Department of Electrical Engineering Eindhoven University of Technology
More informationStability of linear time-varying systems through quadratically parameter-dependent Lyapunov functions
Stability of linear time-varying systems through quadratically parameter-dependent Lyapunov functions Vinícius F. Montagner Department of Telematics Pedro L. D. Peres School of Electrical and Computer
More informationGraph and Controller Design for Disturbance Attenuation in Consensus Networks
203 3th International Conference on Control, Automation and Systems (ICCAS 203) Oct. 20-23, 203 in Kimdaejung Convention Center, Gwangju, Korea Graph and Controller Design for Disturbance Attenuation in
More informationParameterized Linear Matrix Inequality Techniques in Fuzzy Control System Design
324 IEEE TRANSACTIONS ON FUZZY SYSTEMS, VOL. 9, NO. 2, APRIL 2001 Parameterized Linear Matrix Inequality Techniques in Fuzzy Control System Design H. D. Tuan, P. Apkarian, T. Narikiyo, and Y. Yamamoto
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 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 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 informationRank-one LMIs and Lyapunov's Inequality. Gjerrit Meinsma 4. Abstract. We describe a new proof of the well-known Lyapunov's matrix inequality about
Rank-one LMIs and Lyapunov's Inequality Didier Henrion 1;; Gjerrit Meinsma Abstract We describe a new proof of the well-known Lyapunov's matrix inequality about the location of the eigenvalues of a matrix
More informationThe model reduction algorithm proposed is based on an iterative two-step LMI scheme. The convergence of the algorithm is not analyzed but examples sho
Model Reduction from an H 1 /LMI perspective A. Helmersson Department of Electrical Engineering Linkoping University S-581 8 Linkoping, Sweden tel: +6 1 816 fax: +6 1 86 email: andersh@isy.liu.se September
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 informationRobust multi objective H2/H Control of nonlinear uncertain systems using multiple linear model and ANFIS
Robust multi objective H2/H Control of nonlinear uncertain systems using multiple linear model and ANFIS Vahid Azimi, Member, IEEE, Peyman Akhlaghi, and Mohammad Hossein Kazemi Abstract This paper considers
More informationRobust Anti-Windup Compensation for PID Controllers
Robust Anti-Windup Compensation for PID Controllers ADDISON RIOS-BOLIVAR Universidad de Los Andes Av. Tulio Febres, Mérida 511 VENEZUELA FRANCKLIN RIVAS-ECHEVERRIA Universidad de Los Andes Av. Tulio Febres,
More information7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system
7 Stability 7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system ẋ(t) = A x(t), x(0) = x 0, A R n n, x 0 R n. (14) The origin x = 0 is a globally asymptotically
More informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
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 informationLecture Notes DISC Course on. Linear Matrix Inequalities in Control. Carsten Scherer and Siep Weiland
Lecture Notes DISC Course on Linear Matrix Inequalities in Control Carsten Scherer and Siep Weiland Version 2.: April, 1999 2 Preface Objectives In recent years linear matrix inequalities LMI s have emerged
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 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 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 informationApplication of LMI for design of digital control systems
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 56, No. 4, 28 Application of LMI for design of digital control systems A. KRÓLIKOWSKI and D. HORLA Institute of Control and Information
More informationOn the solving of matrix equation of Sylvester type
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 7, No. 1, 2019, pp. 96-104 On the solving of matrix equation of Sylvester type Fikret Ahmadali Aliev Institute of Applied
More informationLINEAR ALGEBRA W W L CHEN
LINEAR ALGEBRA W W L CHEN c W W L Chen, 994, 28. This chapter is available free to all individuals, on the understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,
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 informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More 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 informationIntroduction to Linear Matrix Inequalities (LMIs)
ECE 680 Fall 2017 Introduction to Linear Matrix Inequalities (LMIs) by Stanislaw H. Żak and Guisheng Zhai October 18, 2017 Henry Ford offered customers the Model T Ford in any color the customer wants,
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 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 informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Input-Output and Input-State Linearization Zero Dynamics of Nonlinear Systems Hanz Richter Mechanical Engineering Department Cleveland State University
More informationLMI based output-feedback controllers: γ-optimal versus linear quadratic.
Proceedings of the 17th World Congress he International Federation of Automatic Control Seoul Korea July 6-11 28 LMI based output-feedback controllers: γ-optimal versus linear quadratic. Dmitry V. Balandin
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 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 informationRobust and Optimal Control, Spring 2015
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) D. Linear Matrix Inequality D.1 Convex Optimization D.2 Linear Matrix Inequality(LMI) D.3 Control Design and LMI Formulation
More informationControl for stability and Positivity of 2-D linear discrete-time systems
Manuscript received Nov. 2, 27; revised Dec. 2, 27 Control for stability and Positivity of 2-D linear discrete-time systems MOHAMMED ALFIDI and ABDELAZIZ HMAMED LESSI, Département de Physique Faculté des
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 informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
More informationState feedback gain scheduling for linear systems with time-varying parameters
State feedback gain scheduling for linear systems with time-varying parameters Vinícius F. Montagner and Pedro L. D. Peres Abstract This paper addresses the problem of parameter dependent state feedback
More informationON POLE PLACEMENT IN LMI REGION FOR DESCRIPTOR LINEAR SYSTEMS. Received January 2011; revised May 2011
International Journal of Innovative Computing, Information and Control ICIC International c 2012 ISSN 1349-4198 Volume 8, Number 4, April 2012 pp. 2613 2624 ON POLE PLACEMENT IN LMI REGION FOR DESCRIPTOR
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 informationOptimization Based Output Feedback Control Design in Descriptor Systems
Trabalho apresentado no XXXVII CNMAC, S.J. dos Campos - SP, 017. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics Optimization Based Output Feedback Control Design in
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 informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationFrom Convex Optimization to Linear Matrix Inequalities
Dep. of Information Engineering University of Pisa (Italy) From Convex Optimization to Linear Matrix Inequalities eng. Sergio Grammatico grammatico.sergio@gmail.com Class of Identification of Uncertain
More information1 Continuous-time Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationEfficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph
Efficient robust optimization for robust control with constraints p. 1 Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph Efficient robust optimization
More informationCourse Outline. FRTN10 Multivariable Control, Lecture 13. General idea for Lectures Lecture 13 Outline. Example 1 (Doyle Stein, 1979)
Course Outline FRTN Multivariable Control, Lecture Automatic Control LTH, 6 L-L Specifications, models and loop-shaping by hand L6-L8 Limitations on achievable performance L9-L Controller optimization:
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 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 informationDesign of hybrid control systems for continuous-time plants: from the Clegg integrator to the hybrid H controller
Design of hybrid control systems for continuous-time plants: from the Clegg integrator to the hybrid H controller Luca Zaccarian LAAS-CNRS, Toulouse and University of Trento University of Oxford November
More informationJune Engineering Department, Stanford University. System Analysis and Synthesis. Linear Matrix Inequalities. Stephen Boyd (E.
Stephen Boyd (E. Feron :::) System Analysis and Synthesis Control Linear Matrix Inequalities via Engineering Department, Stanford University Electrical June 1993 ACC, 1 linear matrix inequalities (LMIs)
More informationJosé C. Geromel. Australian National University Canberra, December 7-8, 2017
5 1 15 2 25 3 35 4 45 5 1 15 2 25 3 35 4 45 5 55 Differential LMI in Optimal Sampled-data Control José C. Geromel School of Electrical and Computer Engineering University of Campinas - Brazil Australian
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 informationConvex Optimization 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Convex Optimization 1 Many optimization objectives generated
More informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
More informationEconomic sensor/actuator selection and its application to flexible structure control
To appear in Proc. SPIE Int. Soc. Opt. Eng. 2004 Economic sensor/actuator selection and its application to flexible structure control Robert E. Skelton a and Faming Li a a Department of Mechanical and
More informationLecture 7 : Generalized Plant and LFT form Dr.-Ing. Sudchai Boonto Assistant Professor
Dr.-Ing. Sudchai Boonto Assistant Professor Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Quadratic Gaussian The state space
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 informationEXERCISE SET 5.1. = (kx + kx + k, ky + ky + k ) = (kx + kx + 1, ky + ky + 1) = ((k + )x + 1, (k + )y + 1)
EXERCISE SET 5. 6. The pair (, 2) is in the set but the pair ( )(, 2) = (, 2) is not because the first component is negative; hence Axiom 6 fails. Axiom 5 also fails. 8. Axioms, 2, 3, 6, 9, and are easily
More informationInterval solutions for interval algebraic equations
Mathematics and Computers in Simulation 66 (2004) 207 217 Interval solutions for interval algebraic equations B.T. Polyak, S.A. Nazin Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya
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 informationMulti-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures
Preprints of the 19th World Congress The International Federation of Automatic Control Multi-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures Eric Peterson Harry G.
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationCharacterizing Robust Solution Sets of Convex Programs under Data Uncertainty
Characterizing Robust Solution Sets of Convex Programs under Data Uncertainty V. Jeyakumar, G. M. Lee and G. Li Communicated by Sándor Zoltán Németh Abstract This paper deals with convex optimization problems
More informationStability, Pole Placement, Observers and Stabilization
Stability, Pole Placement, Observers and Stabilization 1 1, The Netherlands DISC Course Mathematical Models of Systems Outline 1 Stability of autonomous systems 2 The pole placement problem 3 Stabilization
More informationLecture 10 - Eigenvalues problem
Lecture 10 - Eigenvalues problem Department of Computer Science University of Houston February 28, 2008 1 Lecture 10 - Eigenvalues problem Introduction Eigenvalue problems form an important class of problems
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 informationMath Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88
Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More informationProblem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations:
EEE58 Exam, Fall 6 AA Rodriguez Rules: Closed notes/books, No calculators permitted, open minds GWC 35, 965-37 Problem (Dynamic Augmentation: State Space Representation) Consider a dynamical system consisting
More informationFixed-Order Robust H Controller Design with Regional Pole Assignment
SUBMITTED 1 Fixed-Order Robust H Controller Design with Regional Pole Assignment Fuwen Yang, Mahbub Gani, and Didier Henrion Abstract In this paper, the problem of designing fixed-order robust H controllers
More informationTheorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).
Linear Systems Notes Lecture Proposition. A M n (R) is positive definite iff all nested minors are greater than or equal to zero. n Proof. ( ): Positive definite iff λ i >. Let det(a) = λj and H = {x D
More informationRobust PID Controller Design for Nonlinear Systems
Robust PID Controller Design for Nonlinear Systems Part II Amin Salar 8700884 Final Project Nonlinear Control Course Dr H.D. Taghirad 1 About the Project In part one we discussed about auto tuning techniques
More informationLinear Systems with Saturating Controls: An LMI Approach. subject to control saturation. No assumption is made concerning open-loop stability and no
Output Feedback Robust Stabilization of Uncertain Linear Systems with Saturating Controls: An LMI Approach Didier Henrion 1 Sophie Tarbouriech 1; Germain Garcia 1; Abstract : The problem of robust controller
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More information2nd Symposium on System, Structure and Control, Oaxaca, 2004
263 2nd Symposium on System, Structure and Control, Oaxaca, 2004 A PROJECTIVE ALGORITHM FOR STATIC OUTPUT FEEDBACK STABILIZATION Kaiyang Yang, Robert Orsi and John B. Moore Department of Systems Engineering,
More informationCSC Linear Programming and Combinatorial Optimization Lecture 10: Semidefinite Programming
CSC2411 - Linear Programming and Combinatorial Optimization Lecture 10: Semidefinite Programming Notes taken by Mike Jamieson March 28, 2005 Summary: In this lecture, we introduce semidefinite programming
More information