A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case
|
|
- Arthur Kennedy
- 5 years ago
- Views:
Transcription
1 A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case ATARU KASE Osaka Institute of Technology Department of Electrical and Electronic Systems Engineering miya, Asahi-ku, Osaka JAPAN Abstract: An interactor matrix plays several important roles in the control systems theory In this paper, we present a simple method to derive the right interactor for tall transfer function matrices using Moore-Penrose pseudoinverse By the presented method, all zeros of the interactor lie at the origin The method will be applied to the inner-outer factorization It will be shown that the stability of the interactor is necessary to calculate the factorization Key ords: interactor matrix, inner-outer factorization, strictly proper plant, continuous-time systems 1 Introduction The inner-outer factorization for stable transfer function matrices is a useful tool for robust controller analysis and synthesis It is known as a reasonable method for computer calculation based on the state space representation 1 However, the method is only useful for proper transfer function In this paper, we will propose to use the interactor 2 to compensate a strictly proper plant In the several control problems, an interactor matrix plays some important roles Although some papers were shown along with the above idea, those method were very complex since the derivations of the interactor were also complex 3-5 Some new calculation methods were presented (6-9) However, they were not easy On the other hand, it is shown that the coefficient matrices of the interactor can be obtained by solving a linear matrix equation of certain type 1 The authors presented a method to solve the equation by using Moore-Penrose pseudoinverse 11 But the method can be only used for square transfer function matrices In this paper, the method presented in 11 will be extended to the case nonsquare transfer function matrices For our application to the inner-outer factorization, we will consider the derivation of the right interactor, unlike our previous work 11 It will be shown that the proposed interactor has all-pass property in the discrete-time, so all of its zeros lie at the origin However, it can not be used since the origin is not stable for the continuous-time systems Thus, we will also show the zeros assignment of the interactor The paper is organized as follows In the next section, the derivation of an interactor will be presented For this purpose, Moore-Penrose pseudoinverse will be employed In Section 3, an application to the innerouter factorization will be shown The result is very simple but some properties should be proved, which will be shown in Appendix Concluding remarks will be given in Section 4 2 Simple Derivation of Interactor Matrix For a given m p (m p) strictly proper transfer function matrix G(s), there exists a p p polynomial matrix L(s), which satisfies the following equation: lim G(s)L(s) = K (full column rank) (1) z Such an L(s) is called a right interactor matrix of G(s) 1 Let (A, B, C) denote a minimal realization of G(s) Setting the coefficient of L(s) by L(s) = sl s w L w = ss w 1 (s)l, (2) 1 Although the definition in 2 is restricted the structure of L(s) (lower triangular), we do not consider such a restriction since it is not essential in this paper ISBN:
2 then the the following equation holds as in 1: T w 1 L = J w 1 (3) L 1 L 2 K L =, J w 1 =, m(w 1) p L i R p p (i = 1, 2,, w) L w CB CAB CA w 1 B CB CA w 2 B T w 1 = CB SI w p (s) = s s w The integer w will be determined later In 11, the special solution ξ of eqn(3) is given by ξ = T w 1 J w 1 using the pseudoinverse T w 1 of T w 1, p = m and K = However, it can not be assumed the special structure of K since p m At this time, considering the structure of J w 1, set L = T w 1 J w 1 = T w 1 (:, 1 : m)k, (4) T w 1 (:, 1 : m) denote the submatrix constituted of the first m-th columns of T w 1 Substituting the above equation to eqn(3), T w 1 T w 1 (:, 1 : m)k = J w 1 (5) Define Λ by Λ = CB CAB CA w 1 B T w 1 (:, 1 : m), (6) the first m-th rows of eqn(5) can be written by ΛK = K (7) That is, if eqn(3) is solvable, its special solution is given by eqn(4) and K must satisfy eqn(7) Eqn(7) means that K is the set of eigenvectors of Λ which correspond to the eigenvalues at λ = 1 Since the definition of T w 1, T w 1T w 1 is a real symmetric matrix, and thus Λ is also a real symmetric Therefore, the geometric multiplicity of the eigenvalue one in Λ equals to the algebraic multiplicity and thus K can be found Therefore, w is the least integer when Λ has the eigenvalue at λ = 1 with multiplicity p K is constituted of corresponding eigenvectors Note that from the singular value decomposition of T w 1 and T w 1, every eigenvalue of T w 1T w 1 is one or zero In the case the eigenvalues of Λ is determined, it should be used the eigenvalues and corresponding eigenvectors of T w 1 T w 1 Lemma 1 For the interactor L(s) derived in the above, the following equation holds: L (s)l(s) = L (8) L (s) = (s 1 ) (9) = s s s w w (Proof) Since T w 1 is the pseudoinverse of T w 1, the following relation must hold Partitioning T w 1 by T w 1 T w 1 = (T w 1 T w 1) T (1) T w 1 = T w 1 (:, 1 : m) P, P R pw m(w 1) (11) and substituting the above into eqn(1), we have T w 1 CB (:, 1 : m) CAB CA w 1 B +P m(w 1) p T w 2 = Tw 1(T T w 1 )T (12) Left multiplying to the both side of the above equation, we have T w 1 CB (:, 1 : m) CAB CA w 1 B + P m(w 1) p T w 2 = Jw 1(T T w 1 )T = K T (T w 1 (:, 1 : m))t = (13) From the definition of an interactor matrix, it follows that L w L w 1 L 2 L w L 3 T w 2 = (14) L w Thus from eqn(13), we have L w L w 1 L 1 L w L 2 L w ISBN:
3 = ( T w 1 (:, 1 : m) CB CAB CA w 1 B L w L w 1 L 1 ) + L w L 2 P m(w 1) p T w 2 L w = T w 1 (:, 1 : m) m p(w 1) K = p p(w 1) L (15) Since the above equation can be rewritten by L w L w 1 L w = LT 1 L w = = 1 L w L w 1 L = = 1 LT 2 L =, the following equation also holds: L 1 L 2 L 1 L w 1 L w 2 L 1 = (16) Therefore, using eqns(15) and (16), L (s)l(s) can be described by L (s)l(s) s 1 L(s) = s 2 L(s) s w L(s) s w 1 s w 2 L w L w 1 L 1 = L w L 2 L 1 I p L w L w 1 L s 1 1 s 1 w s w 1 = p p(w 1) L p p(w 1) s 1 w = L (17) If eqn(8) holds, L(s) is having all-pass property in the discrete-time Therefore, all zeros of the interactor lie at the origin But the origin is not stable for the continuous-time systems Therefore, the zeros of interactor should be assigned to the stable region Theorem 2 Define H := AB A 2 B A w B L Let (Ã, B, C) denote a minimal representation of (A HK C, B, K C) and denote an observer gain which makes the eigenvalues of à := à C be desired ones Then, there exists a matrix L such that T,w L = n p J w 1 B à B à w B C B T,w := CÃw 1 B C B Moreover, defining (18) ˆL(s) := S w (s)l, (19) the zeros of detˆl(s) are equivalent to the eigenvalues of à (Proof) Since H is the maximum uncontrollable observer gain for the system (A, B, K C), there are no finite zeros in (Ã, B, C) Therefore, the inverse system of (Ã, B, C) can be presented by polynomial matrix, say ˆL(s) That is, the zeros of det ˆL(s) are equivalent to the eigenvalues of à Since a right interactor for (A, B, C) is also the one of (Ã, B, K C) 12, ˆL(s) can be represented the polynomial matrix with w-th degrees as shown in eqn(19) Thus, the following relation holds: à B C S w (s)l = ISBN:
4 Expanding the transfer function matrix using Markov parameters, it can be obtained the following equation: C B Cà B CÃw 1 B C B Cà w 2 B L = J w 1 (2) C B C Cà CÃ2 B à B à w 1 B L = Since ( C, à ) is observable, B à B à w 1 B L = and combining the above equations, it can be obtained eqn(19) Example 1 Consider the following transfer function matrix: s + 1 s + 2 G(s) = s + 3 s + 4 s + 5 s + 6 Then, a state space realization (A, B, C) is given by A =, B = C = 1 1 In this example, w = 3 is the first integer when Λ = has the multiple eigenvalue at λ = 1 Then K is given by 1 K = Thus, the right interactor L(s) can be calculated by L(s) = si 2 s 2 I 2 s 3 I 2 L 1 L 2 L = ssi 2 2 (s), det L(s) = 5s 4 Next, design another interactor ˆL(s) which has the desired zeros at s = 12, 14, 16, 18 For the realization, let denote an observer gain which makes the closed-loop s poles be 12, 14, 16 and 18 Setting à := à C, ˆL(s) is given by B à B Ã2 B Ã3 B ˆL(s) = SI 3 2 (s) C B Cà B CÃ2 B C B Cà B C B = I 2 si 2 s 2 I 2 s 3 I I Application to Inner-Outer Factorization For a given m p (m p) proper and stable transfer function matrix G(s) which has no invariant zeros on the imaginary axis, it is known the inner-outer factorization using state space representation if eqn(1) holds with L(s) = In this section, it will be proposed the method applicable for the case L(s), which can be summarized by the following procedures: step 1 Calculate a right interactor L(s) for a given G(s) step 2 Calculate a spectral factorization for G(s)L(s), ie, calculate a stable rational ISBN:
5 function matrix H(s) satisfying L (s)g (s)g(s)l(s) = H (s)h(s)(21) H 1 (s) is stable, G (s) := G T ( s) step 3 Set inner matrix H i (s) and outer matrix H o (s) as follows: H i (s) := G(s)L(s)H 1 (s), H o (s) := H(s)L 1 (s) (22) It is clear that H i (s)h o (s) = G(s) and since H i (s)h i (s) = {G(s)L(s)H 1 (s)} G(s)L(s)H 1 (s) = {H 1 (s)} H (s)h(s)h 1 (s) = (23) hold, H i (z) is an inner Since it is hard to carry out the above calculations based on transfer function matrix, we will present the method based on the state space representation Although we only consider the case G(s) is strictly proper, same idea can be used for the proper case Let (A, B, C) denote a minimal state space realization of G(s) Then, the state space representation of G(s)L(s) is given by (A, ˆB, C, K), Γ := B AB A w 1 B L ˆB := AΓ + BL, K = CΓ, (24) if the coefficient matrix L(s) is given by eqn(2) Consider the following Riccati equation: A T X + XA + C T C (25) (XB + C T K)R 1 (B T X + K T C) = R := K T K (26) Let λ and x T denote an eigenvalue and corresponding eigenvector of A Then, x T B, since (A, B) is controllable On the other hand, x T ˆB = x T B(L + λl λ w L w ) = x T BL(λ) Since all zeros of detl(s) lie in the left half plane, detł(λ) = only if Reλ < Thus, the pair (A, ˆB) is stabilizable and there exists a positive semi definite solution X for Riccati equation (25) Moreover, setting F := R 1 ( ˆB T X + K T C), (27) it is known that X makes A ˆBF be stable Define M(s) and N(s) by M(s) := (si A) 1 ˆB, N(s) := F M(s) (28) Then from Riccati equation, the following relation holds: { + N (s)}r{ + N(s)} = M C (s) I T C C T K M(s) p K T C K T K = {G(s)L(s)} G(s)L(s) (29) Now, define H(s) by H(s) = R 1/2 { + N(s)} = R 1/2 A ˆB (3) F H 1 A ˆBF ˆB (s) = R 1/2, A F := A F ˆBF Then from eqn(22), we can obtain H i (z) = A ˆBF ˆB R 1/2 (31) C KF K As shown in 14, L(s) is also a right interactor for F (si A) 1 B if detl(s) is Hurwitz So, ˆB is the observer gain of the inverted interactorizing for (A, B, F ) (see 13 for detail), ie, L 1 (s) = AF B (32) F Substituting the above equation and eqn(3) to eqn(22), we can obtain H o (s) = R 1/2 A B (33) F Example 2 Same plant is considered as shown in Example 1 Using the solution X of Riccati equation (25), the outer matrix is given by H o (s) = ISBN:
6 as the inner matrix is given by H i (s) = Note that all zeros of detˆl(s) are included as poles of H i (s) However, zeros of the interactor are canceled and H i (s) yields to H i (s) = ie, the parameters of the interactor does not appear 4 Conclusions A simple derivation of a right interactor was given It was shown that the proposed interactor has all-pass property in the discrete-time It was also shown that the zeros assignment method of the interactor Applying the result, it was presented a method of inner-outer factorization based on the state space representation The method was independent on the choice of the interactor This fact is quite natural since the interactor is not an essential but a technical tool for the factorization The result for the discrete-time systems can be found in 14 References: 1 B A Francis, A course in H control theory, Lecture Notes in Control and Information Science, vol88, Springer-Verlag, NY, A olovich and P L Falb, Invariants and canonical form under dynamic compensations, SIAM Journal of Control and Optimization, vol14, pp996-18, X Xin and T Mita, Inner-outer factorization for non-square proper functions with infinite and finite jω-axis zeros, International Journal of Control, vol71, pp , A Varga, Computations of inner-outer factorizations of rational matrices, IEEE Transactions on Automatic Control, vol43, pp , C Oara and A Varga, Computations of general inner-outer and spectral factorizations, IEEE Transactions Automatic Control, vol45, pp , 2 6 R Rogozinski, A P Paplinski and M J Gibbard, An algorithm for the calculation of a nilpotent interactor matrix for linear multivariable systems, IEEE Transactions on Automatic Control, vol32, pp , S Bittani, P Colaneri and M F Mongiovi, Singular filtering via spectral interactor matrix, IEEE Transactions on Automatic Control, vol4, pp , B Huang, S L Shah and H Fujii, The unitary interactor matrix and its estimation using closed-loop data, Journal of Process Control, vol bf 7, pp195-27, T Mita, T K Nam and X Xin, Sliding mode control for invertible systems based on a direct design of interactors, Asian Journal of Control, 5, pp242-25, 22 1 Y Mutoh and R Ortega, Interactor structure estimation for adaptive control of discrete-time multivariable nondecouplable systems, Automatica, vol29, pp , Kase and Y Mutoh, A simple derivation of interactor matrix and its applications, International Journal of Systems Science, vol4, pp , Y Mutoh, A note on the invariant properties of interactors, International Journal of Control, vol62, pp , Y Mutoh and P N Nikiforuk, Inversed interactorizing and triangularizing with an arbitrary pole assignment using the state feedback, IEEE Transactions on Automatic Control, vol37, pp63-633, Kase and Y Mutoh, A simple derivation of right interactor for tall plant and its application to inner-outer factorization, Proceedings of the 19th Mediterranean Conference on Control and Automation, pp , 211 ISBN:
Performance assessment of MIMO systems under partial information
Performance assessment of MIMO systems under partial information H Xia P Majecki A Ordys M Grimble Abstract Minimum variance (MV) can characterize the most fundamental performance limitation of a system,
More informationMultivariable ARMA Systems Making a Polynomial Matrix Proper
Technical Report TR2009/240 Multivariable ARMA Systems Making a Polynomial Matrix Proper Andrew P. Papliński Clayton School of Information Technology Monash University, Clayton 3800, Australia Andrew.Paplinski@infotech.monash.edu.au
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 informationComputing generalized inverse systems using matrix pencil methods
Computing generalized inverse systems using matrix pencil methods A. Varga German Aerospace Center DLR - Oberpfaffenhofen Institute of Robotics and Mechatronics D-82234 Wessling, Germany. Andras.Varga@dlr.de
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 informationOn Positive Real Lemma for Non-minimal Realization Systems
Proceedings of the 17th World Congress The International Federation of Automatic Control On Positive Real Lemma for Non-minimal Realization Systems Sadaaki Kunimatsu Kim Sang-Hoon Takao Fujii Mitsuaki
More information6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski. Solutions to Problem Set 1 1. Massachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Solutions to Problem Set 1 1 Problem 1.1T Consider the
More informationClosed-Loop Structure of Discrete Time H Controller
Closed-Loop Structure of Discrete Time H Controller Waree Kongprawechnon 1,Shun Ushida 2, Hidenori Kimura 2 Abstract This paper is concerned with the investigation of the closed-loop structure of a discrete
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 informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More informationInfinite elementary divisor structure-preserving transformations for polynomial matrices
Infinite elementary divisor structure-preserving transformations for polynomial matrices N P Karampetakis and S Vologiannidis Aristotle University of Thessaloniki, Department of Mathematics, Thessaloniki
More informationDescriptor system techniques in solving H 2 -optimal fault detection problems
Descriptor system techniques in solving H 2 -optimal fault detection problems Andras Varga German Aerospace Center (DLR) DAE 10 Workshop Banff, Canada, October 25-29, 2010 Outline approximate fault detection
More informationMatrix Mathematics. Theory, Facts, and Formulas with Application to Linear Systems Theory. Dennis S. Bernstein
Matrix Mathematics Theory, Facts, and Formulas with Application to Linear Systems Theory Dennis S. Bernstein PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Special Symbols xv Conventions, Notation,
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationAn LQ R weight selection approach to the discrete generalized H 2 control problem
INT. J. CONTROL, 1998, VOL. 71, NO. 1, 93± 11 An LQ R weight selection approach to the discrete generalized H 2 control problem D. A. WILSON², M. A. NEKOUI² and G. D. HALIKIAS² It is known that a generalized
More informationThe Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII
CONTENTS VOLUME VII Control of Linear Multivariable Systems 1 Katsuhisa Furuta,Tokyo Denki University, School of Science and Engineering, Ishizaka, Hatoyama, Saitama, Japan 1. Linear Multivariable Systems
More informationOUTPUT REGULATION OF RÖSSLER PROTOTYPE-4 CHAOTIC SYSTEM BY STATE FEEDBACK CONTROL
International Journal in Foundations of Computer Science & Technology (IJFCST),Vol., No., March 01 OUTPUT REGULATION OF RÖSSLER PROTOTYPE-4 CHAOTIC SYSTEM BY STATE FEEDBACK CONTROL Sundarapandian Vaidyanathan
More informationA Necessary and Sufficient Condition for High-Frequency Robustness of Non-Strictly-Proper Feedback Systems
A Necessary and Sufficient Condition for High-Frequency Robustness of Non-Strictly-Proper Feedback Systems Daniel Cobb Department of Electrical Engineering University of Wisconsin Madison WI 53706-1691
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 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 informationThe Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals
The Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals Timo Hämäläinen Seppo Pohjolainen Tampere University of Technology Department of Mathematics
More informationAn Internal Stability Example
An Internal Stability Example Roy Smith 26 April 2015 To illustrate the concept of internal stability we will look at an example where there are several pole-zero cancellations between the controller and
More informationDiscussion on: Measurable signal decoupling with dynamic feedforward compensation and unknown-input observation for systems with direct feedthrough
Discussion on: Measurable signal decoupling with dynamic feedforward compensation and unknown-input observation for systems with direct feedthrough H.L. Trentelman 1 The geometric approach In the last
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 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 informationOutput Regulation of the Tigan System
Output Regulation of the Tigan System Dr. V. Sundarapandian Professor (Systems & Control Eng.), Research and Development Centre Vel Tech Dr. RR & Dr. SR Technical University Avadi, Chennai-6 6, Tamil Nadu,
More informationlinearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice
3. Eigenvalues and Eigenvectors, Spectral Representation 3.. Eigenvalues and Eigenvectors A vector ' is eigenvector of a matrix K, if K' is parallel to ' and ' 6, i.e., K' k' k is the eigenvalue. If is
More informationA DESIGN METHOD FOR SIMPLE REPETITIVE CONTROLLERS WITH SPECIFIED INPUT-OUTPUT CHARACTERISTIC
International Journal of Innovative Computing, Information Control ICIC International c 202 ISSN 349-498 Volume 8, Number 7(A), July 202 pp. 4883 4899 A DESIGN METHOD FOR SIMPLE REPETITIVE CONTROLLERS
More informationBUMPLESS SWITCHING CONTROLLERS. William A. Wolovich and Alan B. Arehart 1. December 27, Abstract
BUMPLESS SWITCHING CONTROLLERS William A. Wolovich and Alan B. Arehart 1 December 7, 1995 Abstract This paper outlines the design of bumpless switching controllers that can be used to stabilize MIMO plants
More informationMultivariable MRAC with State Feedback for Output Tracking
29 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 1-12, 29 WeA18.5 Multivariable MRAC with State Feedback for Output Tracking Jiaxing Guo, Yu Liu and Gang Tao Department
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 informationChapter 1. Matrix Algebra
ST4233, Linear Models, Semester 1 2008-2009 Chapter 1. Matrix Algebra 1 Matrix and vector notation Definition 1.1 A matrix is a rectangular or square array of numbers of variables. We use uppercase boldface
More informationIndex. book 2009/5/27 page 121. (Page numbers set in bold type indicate the definition of an entry.)
page 121 Index (Page numbers set in bold type indicate the definition of an entry.) A absolute error...26 componentwise...31 in subtraction...27 normwise...31 angle in least squares problem...98,99 approximation
More informationECE 388 Automatic Control
Controllability and State Feedback Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage:
More informationPARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT
PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT Hans Norlander Systems and Control, Department of Information Technology Uppsala University P O Box 337 SE 75105 UPPSALA, Sweden HansNorlander@ituuse
More informationFinite-Time Behavior of Inner Systems
1134 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 7, JULY 2003 Finite-Time Behavior of Inner Systems Jobert H. A. Ludlage, Member, IEEE, Siep Weiland, Anton A. Stoorvogel, Senior Member, IEEE,
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 informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2017 LECTURE 5
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 17 LECTURE 5 1 existence of svd Theorem 1 (Existence of SVD) Every matrix has a singular value decomposition (condensed version) Proof Let A C m n and for simplicity
More informationLecture 7: Positive Semidefinite Matrices
Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.
More informationA design method for two-degree-of-freedom multi-period repetitive controllers for multiple-input/multiple-output systems
A design method for two-degree-of-freedom multi-period repetitive controllers for multiple-input/multiple-output systems Zhongxiang Chen Kou Yamada Tatsuya Sakanushi Iwanori Murakami Yoshinori Ando Nhan
More information3 Stabilization of MIMO Feedback Systems
3 Stabilization of MIMO Feedback Systems 3.1 Notation The sets R and S are as before. We will use the notation M (R) to denote the set of matrices with elements in R. The dimensions are not explicitly
More informationZero controllability in discrete-time structured systems
1 Zero controllability in discrete-time structured systems Jacob van der Woude arxiv:173.8394v1 [math.oc] 24 Mar 217 Abstract In this paper we consider complex dynamical networks modeled by means of state
More informationProperties of Zero-Free Spectral Matrices Brian D. O. Anderson, Life Fellow, IEEE, and Manfred Deistler, Fellow, IEEE
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 54, NO 10, OCTOBER 2009 2365 Properties of Zero-Free Spectral Matrices Brian D O Anderson, Life Fellow, IEEE, and Manfred Deistler, Fellow, IEEE Abstract In
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 informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationSpectral factorization and H 2 -model following
Spectral factorization and H 2 -model following Fabio Morbidi Department of Mechanical Engineering, Northwestern University, Evanston, IL, USA MTNS - July 5-9, 2010 F. Morbidi (Northwestern Univ.) MTNS
More information3.1 Overview 3.2 Process and control-loop interactions
3. Multivariable 3.1 Overview 3.2 and control-loop interactions 3.2.1 Interaction analysis 3.2.2 Closed-loop stability 3.3 Decoupling control 3.3.1 Basic design principle 3.3.2 Complete decoupling 3.3.3
More informationOn linear quadratic optimal control of linear time-varying singular systems
On linear quadratic optimal control of linear time-varying singular systems Chi-Jo Wang Department of Electrical Engineering Southern Taiwan University of Technology 1 Nan-Tai Street, Yungkung, Tainan
More informationSMITH MCMILLAN FORMS
Appendix B SMITH MCMILLAN FORMS B. Introduction Smith McMillan forms correspond to the underlying structures of natural MIMO transfer-function matrices. The key ideas are summarized below. B.2 Polynomial
More informationTHE DESIGN OF ACTIVE CONTROLLER FOR THE OUTPUT REGULATION OF LIU-LIU-LIU-SU CHAOTIC SYSTEM
THE DESIGN OF ACTIVE CONTROLLER FOR THE OUTPUT REGULATION OF LIU-LIU-LIU-SU CHAOTIC SYSTEM Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University
More informationFall 線性系統 Linear Systems. Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian. NTU-EE Sep07 Jan08
Fall 2007 線性系統 Linear Systems Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian NTU-EE Sep07 Jan08 Materials used in these lecture notes are adopted from Linear System Theory & Design, 3rd.
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More informationMaths for Signals and Systems Linear Algebra in Engineering
Maths for Signals and Systems Linear Algebra in Engineering Lectures 13 15, Tuesday 8 th and Friday 11 th November 016 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE
More informationALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA
ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA Kent State University Department of Mathematical Sciences Compiled and Maintained by Donald L. White Version: August 29, 2017 CONTENTS LINEAR ALGEBRA AND
More informationJ-SPECTRAL FACTORIZATION
J-SPECTRAL FACTORIZATION of Regular Para-Hermitian Transfer Matrices Qing-Chang Zhong zhongqc@ieee.org School of Electronics University of Glamorgan United Kingdom Outline Notations and definitions Regular
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 informationOn the Stabilization of Neutrally Stable Linear Discrete Time Systems
TWCCC Texas Wisconsin California Control Consortium Technical report number 2017 01 On the Stabilization of Neutrally Stable Linear Discrete Time Systems Travis J. Arnold and James B. Rawlings Department
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 informationMultivariable Control. Lecture 05. Multivariable Poles and Zeros. John T. Wen. September 14, 2006
Multivariable Control Lecture 05 Multivariable Poles and Zeros John T. Wen September 4, 2006 SISO poles/zeros SISO transfer function: G(s) = n(s) d(s) (no common factors between n(s) and d(s)). Poles:
More information10 Transfer Matrix Models
MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important
More informationReview problems for MA 54, Fall 2004.
Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on
More informationMIT Final Exam Solutions, Spring 2017
MIT 8.6 Final Exam Solutions, Spring 7 Problem : For some real matrix A, the following vectors form a basis for its column space and null space: C(A) = span,, N(A) = span,,. (a) What is the size m n of
More informationw T 1 w T 2. w T n 0 if i j 1 if i = j
Lyapunov Operator Let A F n n be given, and define a linear operator L A : C n n C n n as L A (X) := A X + XA Suppose A is diagonalizable (what follows can be generalized even if this is not possible -
More informationEigenvalues and Eigenvectors
5 Eigenvalues and Eigenvectors 5.2 THE CHARACTERISTIC EQUATION DETERMINANATS n n Let A be an matrix, let U be any echelon form obtained from A by row replacements and row interchanges (without scaling),
More informationAchievable performance of multivariable systems with unstable zeros and poles
Achievable performance of multivariable systems with unstable zeros and poles K. Havre Λ and S. Skogestad y Chemical Engineering, Norwegian University of Science and Technology, N-7491 Trondheim, Norway.
More informationStabilization, Pole Placement, and Regular Implementability
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 47, NO. 5, MAY 2002 735 Stabilization, Pole Placement, and Regular Implementability Madhu N. Belur and H. L. Trentelman, Senior Member, IEEE Abstract In this
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationStability of Hybrid Control Systems Based on Time-State Control Forms
Stability of Hybrid Control Systems Based on Time-State Control Forms Yoshikatsu HOSHI, Mitsuji SAMPEI, Shigeki NAKAURA Department of Mechanical and Control Engineering Tokyo Institute of Technology 2
More informationLecture Notes: Eigenvalues and Eigenvectors. 1 Definitions. 2 Finding All Eigenvalues
Lecture Notes: Eigenvalues and Eigenvectors Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk 1 Definitions Let A be an n n matrix. If there
More informationLearn2Control Laboratory
Learn2Control Laboratory Version 3.2 Summer Term 2014 1 This Script is for use in the scope of the Process Control lab. It is in no way claimed to be in any scientific way complete or unique. Errors should
More informationContents. Preface for the Instructor. Preface for the Student. xvii. Acknowledgments. 1 Vector Spaces 1 1.A R n and C n 2
Contents Preface for the Instructor xi Preface for the Student xv Acknowledgments xvii 1 Vector Spaces 1 1.A R n and C n 2 Complex Numbers 2 Lists 5 F n 6 Digression on Fields 10 Exercises 1.A 11 1.B Definition
More informationWEAK STRUCTURE AT INFINITY AND ROW-BY-ROW DECOUPLING FOR LINEAR DELAY SYSTEMS
KYBERNETIKA VOLUME 40 (2004), NUMBER 2, PAGES 181 195 WEAK STRUCTURE AT INFINITY AND ROW-BY-ROW DECOUPLING FOR LINEAR DELAY SYSTEMS Rabah Rabah and Michel Malabre We consider the row-by-row decoupling
More informationThe Cayley-Hamilton Theorem and the Jordan Decomposition
LECTURE 19 The Cayley-Hamilton Theorem and the Jordan Decomposition Let me begin by summarizing the main results of the last lecture Suppose T is a endomorphism of a vector space V Then T has a minimal
More informationDefinition (T -invariant subspace) Example. Example
Eigenvalues, Eigenvectors, Similarity, and Diagonalization We now turn our attention to linear transformations of the form T : V V. To better understand the effect of T on the vector space V, we begin
More informationQUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY
QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY D.D. Olesky 1 Department of Computer Science University of Victoria Victoria, B.C. V8W 3P6 Michael Tsatsomeros Department of Mathematics
More informationFull-State Feedback Design for a Multi-Input System
Full-State Feedback Design for a Multi-Input System A. Introduction The open-loop system is described by the following state space model. x(t) = Ax(t)+Bu(t), y(t) =Cx(t)+Du(t) () 4 8.5 A =, B =.5.5, C
More informationMin-Max Output Integral Sliding Mode Control for Multiplant Linear Uncertain Systems
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July -3, 27 FrC.4 Min-Max Output Integral Sliding Mode Control for Multiplant Linear Uncertain
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationQuestion: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI?
Section 5. The Characteristic Polynomial Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Property The eigenvalues
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Basic Feedback Analysis & Design
AERO 422: Active Controls for Aerospace Vehicles Basic Feedback Analysis & Design Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University Routh s Stability
More informationINVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN. Leonid Lyubchyk
CINVESTAV Department of Automatic Control November 3, 20 INVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN Leonid Lyubchyk National Technical University of Ukraine Kharkov
More informationInternational Journal of Engineering Research & Science (IJOER) ISSN: [ ] [Vol-2, Issue-11, November- 2016]
Synthesis of Discrete Steady-State Error Free Modal State Controller Based on Predefined Pole Placement Area and Measurable State Variables Mariela Alexandrova Department of Automation, Technical University
More informationLinGloss. A glossary of linear algebra
LinGloss A glossary of linear algebra Contents: Decompositions Types of Matrices Theorems Other objects? Quasi-triangular A matrix A is quasi-triangular iff it is a triangular matrix except its diagonal
More informationLinear Algebra Practice Problems
Math 7, Professor Ramras Linear Algebra Practice Problems () Consider the following system of linear equations in the variables x, y, and z, in which the constants a and b are real numbers. x y + z = a
More informationOutput Regulation of the Arneodo Chaotic System
Vol. 0, No. 05, 00, 60-608 Output Regulation of the Arneodo Chaotic System Sundarapandian Vaidyanathan R & D Centre, Vel Tech Dr. RR & Dr. SR Technical University Avadi-Alamathi Road, Avadi, Chennai-600
More informationA Method to Teach the Parameterization of All Stabilizing Controllers
Preprints of the 8th FAC World Congress Milano (taly) August 8 - September, A Method to Teach the Parameterization of All Stabilizing Controllers Vladimír Kučera* *Czech Technical University in Prague,
More informationThe parameterization of all. of all two-degree-of-freedom strongly stabilizing controllers
The parameterization stabilizing controllers 89 The parameterization of all two-degree-of-freedom strongly stabilizing controllers Tatsuya Hoshikawa, Kou Yamada 2, Yuko Tatsumi 3, Non-members ABSTRACT
More informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationOUTPUT REGULATION OF THE SIMPLIFIED LORENZ CHAOTIC SYSTEM
OUTPUT REGULATION OF THE SIMPLIFIED LORENZ CHAOTIC SYSTEM Sundarapandian Vaidyanathan Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University Avadi, Chennai-600 06, Tamil Nadu, INDIA
More information1030 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 56, NO. 5, MAY 2011
1030 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 56, NO 5, MAY 2011 L L 2 Low-Gain Feedback: Their Properties, Characterizations Applications in Constrained Control Bin Zhou, Member, IEEE, Zongli Lin,
More informationMANY adaptive control methods rely on parameter estimation
610 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 52, NO 4, APRIL 2007 Direct Adaptive Dynamic Compensation for Minimum Phase Systems With Unknown Relative Degree Jesse B Hoagg and Dennis S Bernstein Abstract
More informationCONSIDER the linear discrete-time system
786 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 42, NO. 6, JUNE 1997 Optimal Damping of Forced Oscillations in Discrete-Time Systems Anders Lindquist, Fellow, IEEE, Vladimir A. Yakubovich Abstract In
More information[k,g,gfin] = hinfsyn(p,nmeas,ncon,gmin,gmax,tol)
8 H Controller Synthesis: hinfsyn The command syntax is k,g,gfin = hinfsyn(p,nmeas,ncon,gmin,gmax,tol) associated with the general feedback diagram below. e y P d u e G = F L (P, K) d K hinfsyn calculates
More informationReview of some mathematical tools
MATHEMATICAL FOUNDATIONS OF SIGNAL PROCESSING Fall 2016 Benjamín Béjar Haro, Mihailo Kolundžija, Reza Parhizkar, Adam Scholefield Teaching assistants: Golnoosh Elhami, Hanjie Pan Review of some mathematical
More informationISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION
ISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION Harald K. Wimmer 1 The set of all negative-semidefinite solutions of the CARE A X + XA + XBB X C C = 0 is homeomorphic
More informationProf. Dr.-Ing. Armin Dekorsy Department of Communications Engineering. Stochastic Processes and Linear Algebra Recap Slides
Prof. Dr.-Ing. Armin Dekorsy Department of Communications Engineering Stochastic Processes and Linear Algebra Recap Slides Stochastic processes and variables XX tt 0 = XX xx nn (tt) xx 2 (tt) XX tt XX
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 informationStructured State Space Realizations for SLS Distributed Controllers
Structured State Space Realizations for SLS Distributed Controllers James Anderson and Nikolai Matni Abstract In recent work the system level synthesis (SLS) paradigm has been shown to provide a truly
More information