Peter C. Müller. Introduction. - and the x 2 - subsystems are called the slow and the fast subsystem, correspondingly, cf. (Dai, 1989).
|
|
- Candice Strickland
- 6 years ago
- Views:
Transcription
1 Peter C. Müller Safety Control Engineering University of Wuppertal D-497 Wupperta. Germany Modified Lyapunov Equations for LI Descriptor Systems or linear time-invariant (LI) state space systems it is well-known that its asymptotic stability can be related to solution properties of the Lyapunov matrix equation according to so-called inertia theorems. he question now arises how analogous results can be obtained for LI descriptor systems (singular systems, differential-algebraic equations). he stability behaviour of a LI descriptor system is characterized by the eigenvalues of the related matrix pencil. Additionally, by a quadratic Lyapunov function the stability problem can be discussed by solution properties of a generalized Lyapunov matrix equation including a singular coefficient matrix. o overcome this difficult problem of singularity, the Lyapunov matrix equation will be modified such that a regular Lyapunov matrix equation appears and asymptotic stability is preserved. his aim can be reached by shifting the system matrices in a well defined manner. or that the a priori knowledge of an upper bound of the eigenvalues is assumed. It will be discussed how to get such bound. he paper ends with an inertia theorem where the solution properties of a regular modified Lyapunov matrix equation are uniquely related to the asymptotic stability of the LI descriptor system. Keywords: Descriptor systems, asymptotic stability, Lyapunov matrix equations, inertia theorem Introduction or linear time-invariant state space systems x ɺ (t) = Ax(t) it is well known that its asymptotic stability can be discussed by inertia theorems for the Lyapunov matrix equation P + AP = Q, cf. (Müller, 977). he question now arises how analogous results can be obtained for descriptor systems. In the last two decades the modelling of dynamical systems by descriptor systems (singular systems, differential-algebraic equations) became more and more familiar leading to. A E x(t) ɺ = Ax(t), rke = < n = dim x () for linear time-invariant systems where x is an n-dimensional generalized state vector (descriptor vector) and E is a quadratic singular matrix. It is assumed that the related matrix pencil ( λ E A) is regular, i.e. det( λ E A). hen system () can be decomposed in a slow and in a fast subsystem: there exist two regular matrices R, S such that I A RES =, RAS = () N I where N is a nilpotent matrix of degree k which is called the index of the descriptor system (). hen the decomposition is represented by x x = S x xɺ (t) = A x (t), Nxɺ (t) = x(t) dim( xi ) = ni, n + n = n, (3) where the x - and the x - subsystems are called the slow and the fast subsystem, correspondingly, cf. (Dai, 989). Mechanical Descriptor Systems Lagrange s equations of the first and second kinds are well established in analytical mechanics. hey describe the dynamic behaviour of discrete systems, particularly of multibody systems. he difference between the two kinds consists in the manipulation of the kinematic constraints. If a kinematic description of the system has been given in generalized coordinates which are consistent with the constraints, the Lagrange s equations of the second kind can be derived, leading to a set of differential equations only. However, if a redundant set of coordinates is used to describe kinematically the system containing still some constraints explicitly, then Lagrange s equations of the first kind follow. or LI multibody systems these equations result in M ɺ q(t) + Dq(t) ɺ + Kq(t) = λ + G λ (4) q (t) =, holonomic constraints, (5) G q(t) ɺ + Hq(t) =, nonholonomic constraints. (6) Here, q represents the q-dimensional displacement vector, M, D, K stand for the mass matrix, the dissipative and gyroscopic matrix, and the displacement and circulatory matrix, correspondingly. he f- and g-dimensional vectors λ and λ characterize the Lagrange s multipliers (constraint forces) due to the constraints (5) and (6). he first order representation of section. is obtained by defining the descriptor vector and the matrices x = q qɺ λ λ (7) Presented at XI DINAME International Symposium on Dynamic Problems of Mechanics, ebruary 8th - March 4th, 5, Ouro Preto. MG. Brazil. Paper accepted: June, 5. echnical Editors: J.R.. Arruda and D.A. Rade. 448 / Vol. XXVIII, No. 4, October-December 6 ABCM
2 Modified Lyapunov Equations for LI Descriptor Systems I q Iq M K D G E =, A = (8) H G he equations (4-6) are represented equivalently by Eq. () using the expressions (7,8). It should be mentioned that some properties of the system (4-6) can be characterized by certain system conditions. () he matrix pencil ( λe A) is assumed to be regular: det( λ E A). or that it is necessary that holds. rank = f + g G () he classification of independent holonomic and nonholonomic constraints requires the necessary and sufficient Conditions rank = f + g, rank = f + g G H (9) () (3) he existence of a unique static equilibrium q =, λ =, λ = is guaranteed by K G rank = q + f + g H or equivalently by Eq. () and () of Eq. (3), i.e. how inertia theorems can be developed for linear time-invariant descriptor systems. Because the matrix E is singular, the matrix Q can not be chosen as a regular matrix. he singularity of E results in a difficult problem for the solution of Eq. (3). Until now, only special inertia theorems are known. he results of Ishihara and erra () are restricted to systems of index k = (which excludes mechanical descriptor systems with holonomic constraints being of index k = 3), where (Ishihara and erra, ) gives some corrections of (Lewis, 986). he results of Müller (993), Owens and Debeljkovic (985), Stykel ( a, b) require the calculation of the transformation matrices R, S of Eq. (). But then A of the slow subsystem is also known and the stability behaviour can be discussed by A directly. Recently in (Wang et al, ) a new aspect for the stability discussion of system () has been introduced. If the fast subsystem, cf. Eq. (3), is replaced by an asymptotically stable subsystem with finite eigenvalues, then the stability properties of system () remain unchanged. his replacement can be achieved by changing E into a suitable regular matrix E. Unfortunately, this procedure for determining E again requires the matrices R, S which is disadvantageous as mentioned above. But nevertheless, we shall follow this idea. Modified System he aim is to modify the system () into a system with the matrix pencil ( λê Â) such that - Ê is regular leading to a "regular" Lyapunov matrix equation (9), - asymptotic stability is preserved, - the transformation matrices R, S are not required. Assuming a linear modification Ê = E α A, Â = A + βe (4) + + rank I = q ( f + g) q K I q () G G H H the infinite eigenvalues of system () correspond to ˆ i = < for >, i =,..., n α λ α (5) where (.)+ means the Moore-Penrose inverse of a matrix (.). he remarks (), (), (3) represent a series of requirements in ascending order. If Eq. () is satisfied, then the equations () and (9) are satisfied also. Mechanical descriptor systems (4-6) are a special but important application of the general descriptor systems () of which the asymptotic stability will be considered in detail. Asymptotic Stability he stability behaviour of Eq. () is defined by the eigenvalues of the matrix pencil ( λe A) which coincide with the eigenvalues of the system matrix A of the slow subsystem. By the quadratic Lyapunov function υ = x PEx the stability problem can be also discussed by solution properties of the generalized Lyapunov matrix equation A PE + E PA = Q. (3) Now the problem arises how the asymptotic stability of (), i.e. Re( λ i < ) for all eigenvalues of A, relates to solution properties of the modified system. he finite eigenvalues of system (), i.e. the eigenvalues of A, relate to eigenvalues ˆλ j of the matrix pencil λ Eˆ Aˆ by ( ) λ ˆ λ j β j =, αλˆ j αλ ˆ j +, j =,..., n (6) he half plane Re ˆ λ > is mapped by Eq. (6) into the λ -plane in the interior of the circle with the centre λ ( = β ) α α and the radius r = ( β ) α + α assuming α > according to Eq. (5). or < β < the circle crosses the imaginary axis in the points α ± j β α, cf. ig.. α β β Obviously the stability behaviours of the original system () and the modified system represented by the matrices (4) do not agree. Only in the limit case of,, /α an agreement is obtained. or practical applications J. of the Braz. Soc. of Mech. Sci. & Eng. Copyright 6 by ABCM October-December 6, Vol. XXVIII, No. 4 / 449
3 α = β small, β < (7) b are chosen where b is an upper bound of the eigenvalues λ j b, j =,..., n. (8) Combining Eq. (6) and Eq. (8), the asymptotic stability of the modified system leads to the asymptotic stability of the original system () where the finite eigenvalues λ j, j =,..., n, are located in a crescent-shaped region in the open left λ -plane. An illustration of this case is represented in ig..  PˆÊ + Ê Pˆ  = Qˆ. (9) heorem: If E,  are defined by Eq. (4) with α, β according to Eq. (7,8), then the descriptor system () is asymptotically stable if and only if for at least one symmetric, positive definite matrix Qˆ = Qˆ > (and then for all Qˆ = Qˆ > ) there is a symmetric, positive definite solution matrix Pˆ = Pˆ > of Eq. (9). he proof follows immediately from the standard Lyapunov matrix equation for the system matrix ÂÊ or equivalently for Ê Â. Remarks on the Eigenvalue Bound (8) General Remarks rom state space discussions a theoretical eigenvalue bound b is known according to the slow subsystem (3): b = A () igure. Illustration of mapping (6). where each matrix norm can be applied. However, this bound again has the disadvantage that the transformationmatrices R, S () have to be known. he problem arises to find a bound b without knowing R, S. his problem has been discussed by Müller (4) recently. Still there is no general result with respect to the matrices E, A, but for the special case of semi-explicit systems with uniform index nice results have been found. Semi-Explicit Systems with Uniform Index Very often the descriptor system () is represented in semiexplicit form xɺ (t) = Ax(t) + Ax (t) () = A x(t) + A x (t) () (Remark: Here the notation x, x is different to that of Eq. (3); in the following it will be referred to the new notation of Eq. (, ) only.) igure. Stability region of system (). Modified Lyapunov Matrix Equation Summarizing the results of section 3 the stability behaviour of () can be discussed by an inertia theorem of the modified Lyapunov matrix equation he index of a LI descriptor system, which was mentioned with Eq. (), indicates how often the algebraic equations () have to be differentiated to obtain ordinary differential equations for xɺ (t)... additionally. or one number less, k, the ( ) x = algebraic equations are solvable with respect to x. Each singular algebraic equation of the vector algebraic equation () may lead to an individual index. he utmost individual index represents the system index k. In the following it is assumed for simplicity of notation that the individual indices agree such that a uniform index is assumed. In the following the three cases of uniform indices k =, k =, k = 3 are considered. Uniform Index If system (, ) has index k = then A is regular and is valid. hen Eq. () runs as x (t) = A A x(t) (3) 45 / Vol. XXVIII, No. 4, October-December 6 ABCM
4 Modified Lyapunov Equations for LI Descriptor Systems ( A A A A ) x (t) xɺ (t) = (4) representing a state space system. his results immediately to the upper bound Uniform Index λ j b = A AAA (5) Here A and regularity of AA is assumed. Differentiating Eq. () leads to A xɺ (t) = A(Ax(t) + Ax (t)) and therefore x (t) = (AA ) AAx(t) (6) is obtained. Inserting Eq. (6) into Eq. () results in xɺ (t) = [A A (AA ) AA]x(t) (7) rom here the upper bound is obtained. λ j b = A A(AA ) AA (8) Uniform Index 3 In this case the assumptions on system (, ) are A =, AA =, AAA regular. hen x (t) = (A A A ) A A x (t) (9) holds resulting in the differential equation x he bound becomes ɺ (t) = [A A (A A A ) A A ]x (t) (3) j b = A A(AAA ) AA λ (3) In this special cases eigenvalue bounds are available. But with increasing index the calculation becomes more and more expensive. Example: Mechanical descriptor system with holonomic constraints or example the mechanical descriptor system (4) with holonomic constraints (5) is considered. hen a semi-explicit representation (, ) is given by x = [q qɺ ], x = λ; (3) Assuming a regular mass matrix M and independent holonomic constraints, rank = f, then it is easily shown that the mechanical descriptor system has uniform index k = 3: A, AA =, AAA = M = (34) According to Eq. (9) the constraint forces are determined by q(t) [ M K M D] q(t) λ (t) = (M ) (35) ɺ System (3) runs as with Iq xɺ (t) = x(t) (36) K D K = PM K, D = PM D, P = Iq M (M ) (37) he upper eigenvalue bound (3) is given by I b = q K D (38) In the case of symmetric matrices trices M = M, D = D, K = K this bound may be replaced by the bound ( D, K ) b = max (39) 8 8 which follows by a discussion of the Rayleigh quotients according to a second order representation of Eq. (36). he index s characterizes the spectral norm. Conclusions It has been shown that for a matrix pencil ( λe A) a modified matrix pencil ( Ê Â) λ can be assigned such that the asymptotic stability of the new matrix pencil includes the asymptotic stability of the original pencil. he advantage of this modification consists in the regularity of the matrix ˆE and thus in the regularity of the modified Lyapunov matrix equation (9). he asymptotic stability can be guaranteed by the inertia theorem of section 4. his result does not require any knowledge on the transformation matrices R, S for the Weierstrass-Kronecker representation () of the system, but it uses an upper bound (8) of the eigenvalues. his problem has been discussed in section 5 showing some first results but still simpler solutions are desired. A = M K A =, M Iq, M D A =, A = (33) References Dai, L., 989, Singular Control System, Lecture Notes in Control and Information Sciences, Vol. 8, Berlin- Heidelberg, 33p. Ishihara, J.Y. and erra, M.H.,, On the Lyapunov heorem for Singular Systems, IEEE rans. Autom. Control, Vol. 47, pp Lewis,.L., 986, A Survey of Linear Singular Systems, Circuits, Syst., Sig. Proc., Vol. 5, pp Müller, P.C., 977, Stability and Matrices (in German), Berlin- Heidelberg, p. J. of the Braz. Soc. of Mech. Sci. & Eng. Copyright 6 by ABCM October-December 6, Vol. XXVIII, No. 4 / 45
5 Müller, P.C., 993, Stability of Linear Mechanical Systems with Holonomic Constraints, Appl. Mech. Rev., Vol. 46, No., Part, pp. S6-S64. Müller, P.C., 4, Eigenvalue Bounds of Matrix Pencils, PAAM- Proc. Appl. Math. Mech., Vol. 4, to appear. Owens, D.H. and Debeljkovic, D.L., 985, Consistency and Lyapunov Stability of Linear Descriptor Systems: A General Analysis, IMA J. Math. Control & Inform., Vol., pp Stykel,., a, Stability and Inertia heorems for Generalized Lyapunov Equations, Linear Algebra and its Applications, Vol. 355, pp Stykel,., b, On Criteria for Asymptotic Stability of Differential- Algebraic Equations, Z. Angew. Math. Mech., Vol. 8, pp Wang, Q., Zhang, Q.L., Zhang, G.S. and Liu, W.Q.,, Lyapunov Equations with Positive Definite Solution for Descriptor Systems, Proc. 4th Asian Control Conf. (ASCC), WA9-6, Singapore, pp / Vol. XXVIII, No. 4, October-December 6 ABCM
The 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 informationDetermining the order of minimal realization of descriptor systems without use of the Weierstrass canonical form
Journal of Mathematical Modeling Vol. 3, No. 1, 2015, pp. 91-101 JMM Determining the order of minimal realization of descriptor systems without use of the Weierstrass canonical form Kamele Nassiri Pirbazari
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. III Estimation and Compensation of Nonlinear Perturbations by Disturbance Observers - Peter C.
ESTIMATION AND COMPENSATION OF NONLINEAR PERTURBATIONS BY DISTURBANCE OBSERVERS Peter C. Müller University of Wuppertal, Germany Keywords: Closed-loop control system, Compensation of nonlinearities, Disturbance
More informationStability and Inertia Theorems for Generalized Lyapunov Equations
Published in Linear Algebra and its Applications, 355(1-3, 2002, pp. 297-314. Stability and Inertia Theorems for Generalized Lyapunov Equations Tatjana Stykel Abstract We study generalized Lyapunov equations
More informationonly nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr
The discrete algebraic Riccati equation and linear matrix inequality nton. Stoorvogel y Department of Mathematics and Computing Science Eindhoven Univ. of Technology P.O. ox 53, 56 M Eindhoven The Netherlands
More informationAnalytical formulas for calculating the extremal ranks and inertias of A + BXB when X is a fixed-rank Hermitian matrix
Analytical formulas for calculating the extremal ranks and inertias of A + BXB when X is a fixed-rank Hermitian matrix Yongge Tian CEMA, Central University of Finance and Economics, Beijing 100081, China
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 informationDESCRIPTOR FRACTIONAL LINEAR SYSTEMS WITH REGULAR PENCILS
Int. J. Appl. Math. Comput. Sci., 3, Vol. 3, No., 39 35 DOI:.478/amcs-3-3 DESCRIPTOR FRACTIONAL LINEAR SYSTEMS WITH REGULAR PENCILS TADEUSZ KACZOREK Faculty of Electrical Engineering Białystok Technical
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 informationStructured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control
Electronic Journal of Linear Algebra Volume 34 Volume 34 08) Article 39 08 Structured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control Christian Mehl Technische Universitaet
More informationModel order reduction of electrical circuits with nonlinear elements
Model order reduction of electrical circuits with nonlinear elements Andreas Steinbrecher and Tatjana Stykel 1 Introduction The efficient and robust numerical simulation of electrical circuits plays a
More informationContour integral solutions of Sylvester-type matrix equations
Contour integral solutions of Sylvester-type matrix equations Harald K. Wimmer Mathematisches Institut, Universität Würzburg, 97074 Würzburg, Germany Abstract The linear matrix equations AXB CXD = E, AX
More informationELA THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES
Volume 22, pp. 480-489, May 20 THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES XUZHOU CHEN AND JUN JI Abstract. In this paper, we study the Moore-Penrose inverse
More informationMATH 425-Spring 2010 HOMEWORK ASSIGNMENTS
MATH 425-Spring 2010 HOMEWORK ASSIGNMENTS Instructor: Shmuel Friedland Department of Mathematics, Statistics and Computer Science email: friedlan@uic.edu Last update April 18, 2010 1 HOMEWORK ASSIGNMENT
More informationAnalytical formulas for calculating extremal ranks and inertias of quadratic matrix-valued functions and their applications
Analytical formulas for calculating extremal ranks and inertias of quadratic matrix-valued functions and their applications Yongge Tian CEMA, Central University of Finance and Economics, Beijing 100081,
More informationStability of Equilibrium Positions of Mechanical Systems with Switched Force Fields
SCIETIFIC PUBLICATIOS OF THE STATE UIVERSITY OF OVI PAZAR SER. A: APPL. MATH. IFORM. AD MECH. vol. 4, 2 2012, 35-39 Stability of Equilibrium Positions of Mechanical Systems with Switched Force Fields A.
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 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 informationLösung des linear-quadratischen Optimierungs- problems für lineare zeitinvariante Deskriptorsysteme. mittels Riccati-Matrizengleichungen
Lösung des linear-quadratischen Optimierungs- problems für lineare zeitinvariante Deskriptorsysteme mittels Riccati-Matrizengleichungen Peter C. Müller Sicherheitstechnische Regelungs- und Messtechnik
More informationA q x k+q + A q 1 x k+q A 0 x k = 0 (1.1) where k = 0, 1, 2,..., N q, or equivalently. A(σ)x k = 0, k = 0, 1, 2,..., N q (1.
A SPECTRAL CHARACTERIZATION OF THE BEHAVIOR OF DISCRETE TIME AR-REPRESENTATIONS OVER A FINITE TIME INTERVAL E.N.Antoniou, A.I.G.Vardulakis, N.P.Karampetakis Aristotle University of Thessaloniki Faculty
More informationFirst, we review some important facts on the location of eigenvalues of matrices.
BLOCK NORMAL MATRICES AND GERSHGORIN-TYPE DISCS JAKUB KIERZKOWSKI AND ALICJA SMOKTUNOWICZ Abstract The block analogues of the theorems on inclusion regions for the eigenvalues of normal matrices are given
More informationLyapunov Stability of Linear Predictor Feedback for Distributed Input Delays
IEEE TRANSACTIONS ON AUTOMATIC CONTROL VOL. 56 NO. 3 MARCH 2011 655 Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays Nikolaos Bekiaris-Liberis Miroslav Krstic In this case system
More informationStability of interval positive continuous-time linear systems
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 66, No. 1, 2018 DOI: 10.24425/119056 Stability of interval positive continuous-time linear systems T. KACZOREK Białystok University of
More informationThe symmetric minimal rank solution of the matrix equation AX=B and the optimal approximation
Electronic Journal of Linear Algebra Volume 18 Volume 18 (2009 Article 23 2009 The symmetric minimal rank solution of the matrix equation AX=B and the optimal approximation Qing-feng Xiao qfxiao@hnu.cn
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 informationAnalysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits
Computer Applications in Electrical Engineering Vol. 4 Analysis o the regularity pointwise completeness pointwise generacy o descriptor linear electrical circuits Tadeusz Kaczorek Białystok University
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 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 informationPerturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop
Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael
More informationRobust Input-Output Energy Decoupling for Uncertain Singular Systems
International Journal of Automation and Computing 1 (25) 37-42 Robust Input-Output Energy Decoupling for Uncertain Singular Systems Xin-Zhuang Dong, Qing-Ling Zhang Institute of System Science, Northeastern
More informationBlock companion matrices, discrete-time block diagonal stability and polynomial matrices
Block companion matrices, discrete-time block diagonal stability and polynomial matrices Harald K. Wimmer Mathematisches Institut Universität Würzburg D-97074 Würzburg Germany October 25, 2008 Abstract
More informationOPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY
published in IMA Journal of Numerical Analysis (IMAJNA), Vol. 23, 1-9, 23. OPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY SIEGFRIED M. RUMP Abstract. In this note we give lower
More informationOptimization problems on the rank and inertia of the Hermitian matrix expression A BX (BX) with applications
Optimization problems on the rank and inertia of the Hermitian matrix expression A BX (BX) with applications Yongge Tian China Economics and Management Academy, Central University of Finance and Economics,
More informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationBOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION
K Y BERNETIKA VOLUM E 46 ( 2010), NUMBER 4, P AGES 655 664 BOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION Guang-Da Hu and Qiao Zhu This paper is concerned with bounds of eigenvalues of a complex
More informationThe Skew-Symmetric Ortho-Symmetric Solutions of the Matrix Equations A XA = D
International Journal of Algebra, Vol. 5, 2011, no. 30, 1489-1504 The Skew-Symmetric Ortho-Symmetric Solutions of the Matrix Equations A XA = D D. Krishnaswamy Department of Mathematics Annamalai University
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 informationMatrix Inequalities by Means of Block Matrices 1
Mathematical Inequalities & Applications, Vol. 4, No. 4, 200, pp. 48-490. Matrix Inequalities by Means of Block Matrices Fuzhen Zhang 2 Department of Math, Science and Technology Nova Southeastern University,
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 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 informationResearch Article Stabilization Analysis and Synthesis of Discrete-Time Descriptor Markov Jump Systems with Partially Unknown Transition Probabilities
Research Journal of Applied Sciences, Engineering and Technology 7(4): 728-734, 214 DOI:1.1926/rjaset.7.39 ISSN: 24-7459; e-issn: 24-7467 214 Maxwell Scientific Publication Corp. Submitted: February 25,
More informationarxiv: v1 [cs.na] 7 Jul 2017
AE regularity of interval matrices Milan Hladík October 4, 2018 arxiv:1707.02102v1 [cs.na] 7 Jul 2017 Abstract Consider a linear system of equations with interval coefficients, and each interval coefficient
More informationTHE SPECTRAL DIAMETER IN BANACH ALGEBRAS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume»!. Number 1, Mav 1984 THE SPECTRAL DIAMETER IN BANACH ALGEBRAS SANDY GRABINER1 Abstract. The element a is in the center of the Banach algebra A modulo
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 informationA Delay-dependent Condition for the Exponential Stability of Switched Linear Systems with Time-varying Delay
A Delay-dependent Condition for the Exponential Stability of Switched Linear Systems with Time-varying Delay Kreangkri Ratchagit Department of Mathematics Faculty of Science Maejo University Chiang Mai
More informationResearch Article On Exponential Stability Conditions of Descriptor Systems with Time-Varying Delay
Applied Mathematics Volume 2012, Article ID 532912, 12 pages doi:10.1155/2012/532912 Research Article On Exponential Stability Conditions of Descriptor Systems with Time-Varying Delay S. Cong and Z.-B.
More informationDiagonal matrix solutions of a discrete-time Lyapunov inequality
Diagonal matrix solutions of a discrete-time Lyapunov inequality Harald K. Wimmer Mathematisches Institut Universität Würzburg D-97074 Würzburg, Germany February 3, 1997 Abstract Diagonal solutions of
More informationAdaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties
Australian Journal of Basic and Applied Sciences, 3(1): 308-322, 2009 ISSN 1991-8178 Adaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties M.R.Soltanpour, M.M.Fateh
More informationAppendix A Equations of Motion in the Configuration and State Spaces
Appendix A Equations of Motion in the Configuration and State Spaces A.1 Discrete Linear Systems A.1.1 Configuration Space Consider a system with a single degree of freedom and assume that the equation
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 informationNewton Method with Adaptive Step-Size for Under-Determined Systems of Equations
Newton Method with Adaptive Step-Size for Under-Determined Systems of Equations Boris T. Polyak Andrey A. Tremba V.A. Trapeznikov Institute of Control Sciences RAS, Moscow, Russia Profsoyuznaya, 65, 117997
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 informationSome inequalities for sum and product of positive semide nite matrices
Linear Algebra and its Applications 293 (1999) 39±49 www.elsevier.com/locate/laa Some inequalities for sum and product of positive semide nite matrices Bo-Ying Wang a,1,2, Bo-Yan Xi a, Fuzhen Zhang b,
More informationA BOUNDARY VALUE PROBLEM FOR LINEAR PDAEs
Int. J. Appl. Math. Comput. Sci. 2002 Vol.12 No.4 487 491 A BOUNDARY VALUE PROBLEM FOR LINEAR PDAEs WIESŁAW MARSZAŁEK ZDZISŁAW TRZASKA DeVry College of Technology 630 US Highway One North Brunswick N8902
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationTrajectory Planning from Multibody System Dynamics
Trajectory Planning from Multibody System Dynamics Pierangelo Masarati Politecnico di Milano Dipartimento di Ingegneria Aerospaziale Manipulators 2 Manipulator: chain of
More informationMultiplicative Perturbation Bounds of the Group Inverse and Oblique Projection
Filomat 30: 06, 37 375 DOI 0.98/FIL67M Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Multiplicative Perturbation Bounds of the Group
More informationAN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES
AN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES ZHONGYUN LIU AND HEIKE FAßBENDER Abstract: A partially described inverse eigenvalue problem
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationV.B. LARIN 1. Keywords: unilateral quadratic matrix equation, matrix sign function, parameters updating.
TWMS Jour. Pure Appl. Math., V.3, N.2, 2012, pp.202-209 THE UNILATERAL QUADRATIC MATRIX EQUATION AND PROBLEM OF UPDATING OF PARAMETERS OF MODEL V.B. LARIN 1 Abstract. The algorithm of construction of solutions
More informationSpectral Algorithms I. Slides based on Spectral Mesh Processing Siggraph 2010 course
Spectral Algorithms I Slides based on Spectral Mesh Processing Siggraph 2010 course Why Spectral? A different way to look at functions on a domain Why Spectral? Better representations lead to simpler solutions
More informationThe DMP Inverse for Rectangular Matrices
Filomat 31:19 (2017, 6015 6019 https://doi.org/10.2298/fil1719015m Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://.pmf.ni.ac.rs/filomat The DMP Inverse for
More informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
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 informationSUCCESSIVE POLE SHIFTING USING SAMPLED-DATA LQ REGULATORS. Sigeru Omatu
SUCCESSIVE POLE SHIFING USING SAMPLED-DAA LQ REGULAORS oru Fujinaka Sigeru Omatu Graduate School of Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Sakai, 599-8531 Japan Abstract: Design of sampled-data
More informationAn LMI Approach to Robust Controller Designs of Takagi-Sugeno fuzzy Systems with Parametric Uncertainties
An LMI Approach to Robust Controller Designs of akagi-sugeno fuzzy Systems with Parametric Uncertainties Li Qi and Jun-You Yang School of Electrical Engineering Shenyang University of echnolog Shenyang,
More informationLyapunov Stability Analysis: Open Loop
Copyright F.L. Lewis 008 All rights reserved Updated: hursday, August 8, 008 Lyapunov Stability Analysis: Open Loop We know that the stability of linear time-invariant (LI) dynamical systems can be determined
More informationLinear Quadratic Regulator (LQR) II
Optimal Control, Guidance and Estimation Lecture 11 Linear Quadratic Regulator (LQR) II Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Outline Summary o LQR design
More informationStructural Properties of LTI Singular Systems under Output Feedback
Structural Properties of LTI Singular Systems under Output Feedback Runyi Yu Dept. of Electrical and Electronic Eng. Eastern Mediterranean University Gazimagusa, Mersin 10, Turkey runyi.yu@emu.edu.tr Dianhui
More informationApproximation of the Linearized Boussinesq Equations
Approximation of the Linearized Boussinesq Equations Alan Lattimer Advisors Jeffrey Borggaard Serkan Gugercin Department of Mathematics Virginia Tech SIAM Talk Series, Virginia Tech, April 22, 2014 Alan
More informationb i (µ, x, s) ei ϕ(x) µ s (dx) ds (2) i=1
NONLINEAR EVOLTION EQATIONS FOR MEASRES ON INFINITE DIMENSIONAL SPACES V.I. Bogachev 1, G. Da Prato 2, M. Röckner 3, S.V. Shaposhnikov 1 The goal of this work is to prove the existence of a solution to
More informationIterative Solution of a Matrix Riccati Equation Arising in Stochastic Control
Iterative Solution of a Matrix Riccati Equation Arising in Stochastic Control Chun-Hua Guo Dedicated to Peter Lancaster on the occasion of his 70th birthday We consider iterative methods for finding the
More informationReview of Controllability Results of Dynamical System
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 4 Ver. II (Jul. Aug. 2017), PP 01-05 www.iosrjournals.org Review of Controllability Results of Dynamical System
More informationContents. 0.1 Notation... 3
Contents 0.1 Notation........................................ 3 1 A Short Course on Frame Theory 4 1.1 Examples of Signal Expansions............................ 4 1.2 Signal Expansions in Finite-Dimensional
More informationSolutions of a constrained Hermitian matrix-valued function optimization problem with applications
Solutions of a constrained Hermitian matrix-valued function optimization problem with applications Yongge Tian CEMA, Central University of Finance and Economics, Beijing 181, China Abstract. Let f(x) =
More informationIN THIS PAPER, we consider a class of continuous-time recurrent
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 51, NO. 4, APRIL 2004 161 Global Output Convergence of a Class of Continuous-Time Recurrent Neural Networks With Time-Varying Thresholds
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 information/$ IEEE
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 55, NO. 9, SEPTEMBER 2008 937 Analytical Stability Condition of the Latency Insertion Method for Nonuniform GLC Circuits Subramanian N.
More informationRank and inertia optimizations of two Hermitian quadratic matrix functions subject to restrictions with applications
Rank and inertia optimizations of two Hermitian quadratic matrix functions subject to restrictions with applications Yongge Tian a, Ying Li b,c a China Economics and Management Academy, Central University
More informationThere are six more problems on the next two pages
Math 435 bg & bu: Topics in linear algebra Summer 25 Final exam Wed., 8/3/5. Justify all your work to receive full credit. Name:. Let A 3 2 5 Find a permutation matrix P, a lower triangular matrix L with
More informationSingular Value Inequalities for Real and Imaginary Parts of Matrices
Filomat 3:1 16, 63 69 DOI 1.98/FIL16163C Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Singular Value Inequalities for Real Imaginary
More informationON COLISSIONS IN NONHOLONOMIC SYSTEMS
ON COLISSIONS IN NONHOLONOMIC SYSTEMS DMITRY TRESCHEV AND OLEG ZUBELEVICH DEPT. OF THEORETICAL MECHANICS, MECHANICS AND MATHEMATICS FACULTY, M. V. LOMONOSOV MOSCOW STATE UNIVERSITY RUSSIA, 119899, MOSCOW,
More informationPublications: Charles Fulton. Papers On Sturm-Liouville Theory
Publications: Charles Fulton Papers On Sturm-Liouville Theory References [1] C. Fulton, Parametrizations of Titchmarsh s m(λ)- Functions in the limit circle case, Trans. Amer. Math. Soc. 229, (1977), 51-63.
More informationMultiple Degree of Freedom Systems. The Millennium bridge required many degrees of freedom to model and design with.
Multiple Degree of Freedom Systems The Millennium bridge required many degrees of freedom to model and design with. The first step in analyzing multiple degrees of freedom (DOF) is to look at DOF DOF:
More informationOn Monoids over which All Strongly Flat Right S-Acts Are Regular
Æ26 Æ4 ² Vol.26, No.4 2006µ11Â JOURNAL OF MATHEMATICAL RESEARCH AND EXPOSITION Nov., 2006 Article ID: 1000-341X(2006)04-0720-05 Document code: A On Monoids over which All Strongly Flat Right S-Acts Are
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationA Cross-Associative Neural Network for SVD of Nonsquared Data Matrix in Signal Processing
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 12, NO. 5, SEPTEMBER 2001 1215 A Cross-Associative Neural Network for SVD of Nonsquared Data Matrix in Signal Processing Da-Zheng Feng, Zheng Bao, Xian-Da Zhang
More informationA note on the unique solution of linear complementarity problem
COMPUTATIONAL SCIENCE SHORT COMMUNICATION A note on the unique solution of linear complementarity problem Cui-Xia Li 1 and Shi-Liang Wu 1 * Received: 13 June 2016 Accepted: 14 November 2016 First Published:
More informationSTABILITY ANALYSIS FOR DISCRETE T-S FUZZY SYSTEMS
INERNAIONAL JOURNAL OF INFORMAION AND SYSEMS SCIENCES Volume, Number 3-4, Pages 339 346 c 005 Institute for Scientific Computing and Information SABILIY ANALYSIS FOR DISCREE -S FUZZY SYSEMS IAOGUANG YANG,
More informationANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3
ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3 ISSUED 24 FEBRUARY 2018 1 Gaussian elimination Let A be an (m n)-matrix Consider the following row operations on A (1) Swap the positions any
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationOn the existence of unramified p-extensions with prescribed Galois group. Osaka Journal of Mathematics. 47(4) P P.1165
Title Author(s) Citation On the existence of unramified p-extensions with prescribed Galois group Nomura, Akito Osaka Journal of Mathematics. 47(4) P.1159- P.1165 Issue Date 2010-12 Text Version publisher
More informationMATH 5640: Functions of Diagonalizable Matrices
MATH 5640: Functions of Diagonalizable Matrices Hung Phan, UMass Lowell November 27, 208 Spectral theorem for diagonalizable matrices Definition Let V = X Y Every v V is uniquely decomposed as u = x +
More informationA Note on Simple Nonzero Finite Generalized Singular Values
A Note on Simple Nonzero Finite Generalized Singular Values Wei Ma Zheng-Jian Bai December 21 212 Abstract In this paper we study the sensitivity and second order perturbation expansions of simple nonzero
More informationSpectral inequalities and equalities involving products of matrices
Spectral inequalities and equalities involving products of matrices Chi-Kwong Li 1 Department of Mathematics, College of William & Mary, Williamsburg, Virginia 23187 (ckli@math.wm.edu) Yiu-Tung Poon Department
More informationDesign of State Observer for a Class of Non linear Singular Systems Described by Takagi-Sugeno Model
Contemporary Engineering Sciences, Vol. 6, 213, no. 3, 99-19 HIKARI Ltd, www.m-hikari.com Design of State Observer for a Class of Non linear Singular Systems Described by Takagi-Sugeno Model Mohamed Essabre
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationSchur decomposition in the scaled boundary finite element method in elastostatics
IOP Conference Series: Materials Science and Engineering Schur decomposition in the scaled boundary finite element method in elastostatics o cite this article: M Li et al 010 IOP Conf. Ser.: Mater. Sci.
More informationOn the simultaneous diagonal stability of a pair of positive linear systems
On the simultaneous diagonal stability of a pair of positive linear systems Oliver Mason Hamilton Institute NUI Maynooth Ireland Robert Shorten Hamilton Institute NUI Maynooth Ireland Abstract In this
More informationOn the closures of orbits of fourth order matrix pencils
On the closures of orbits of fourth order matrix pencils Dmitri D. Pervouchine Abstract In this work we state a simple criterion for nilpotentness of a square n n matrix pencil with respect to the action
More information