Singular Sylvester Equation and Its Applications

Total Page:16

File Type:pdf, Size:1020Kb

Singular Sylvester Equation and Its Applications University of Niˇs Faculty of Sciences and Mathematics Bogdan D. Dord¯evi´c- Singular Sylvester equation and its applications DOCTORAL DISSERTATION Niˇs,2020. Univerzitet u Niˇsu Prirodno-matematiˇckifakultet Bogdan D. Dord¯evi´c- Singularna Silvesterova jednaˇcinai njene primene DOKTORSKA DISERTACIJA Niˇs,2020. Data on Doctoral Dissertation Doctoral Nebojša Dinčić, full professor at Faculty of Sciences and mathematics, University of Niš supervisor: Title: Singular Sylvester equation and its applications Abstract: This thesis concerns singular Sylvester operator equations, that is, equations of the form AX- XB=C, under the premise that they are either unsolvable or have infinitely many solutions. The equation is studied in different cases, first in the matrix case, then in the case when A, B and C are bounded linear operators on Banach spaces, and finally in the case when A and B are closed linear operators defined on Banach or Hilbert spaces. In each of these cases, solvability conditions are derived and then, under those conditions, the initial equation is solved. Exact solutions are obtained in their closed forms, and their classification is conducted. It is shown that all solutions are obtained in the manner illustrated in this thesis. Special attention is dedicated to approximation schemes of the solutions. Obtained results are illustrated on some contemporary problems from operator theory, among which are spectral problems of bounded and unbounded linear operators, Sturm-Liouville inverse problems and some operator equations from quantum mechanics. Scientific Mathematics field: Scientific Functional analysis discipline: Sylvester equation; operator equations; operator algebras; spectral theory of operators; closed Key words: operators; Fredholm theory; Banach algebras. UDC: 517.983.23; 517.983.24; 517.983.5; 517.986.3; 517.984; 517.984.3; 517.984.4 CERIF P140: Series, Fourier analysis, functional analysis classification: Creative CC BY-NC-ND common licence: Подаци о докторској дисертацији Ментор: Небојша Динчић, редовни професор Природно-математичког факултета Универзитета у Нишу Наслов: Сингуларна Силвестерова једначина и њене примене Резиме: У овој дисертацији се изучава сингуларна Силвестерова операторска једначина, односно, операторска једначина облика AX-XB=C, под претпоставком да је она или нерешива или да има бесконачно много решења. Једначина се посматра у више разчличитих случаја, најпре у матричном случају, затим у случају када су у питању ограничени линеарни оператори на Банаховим просторима и коначно у случају када су у питању затворени линеарни оператори на Банаховим или Хилбертовим просторима. У сваком од поменутих сценарија се прво изводе довољни услови решивости полазне једначине, а онда се под тим претпоставкама прелази на њено решавање. Долази се до егзактних решења у затвореној форми, те се прелази на њихову класификацију и карактеризацију, односно, показује се да су изведеним посупцима обухваћена сва могућа решења сингуларне Силвестерове једначине. Посебна пажња је посвећена апроксимацијама решења. Добијени резултати су илустровани на неким савременим проблемима из теорије оператора, као што су спектрални проблеми ограничених и неограничених линеарних оператора, инверзни проблеми Штурм-Лиувилове теорије и операторске једначине које се јављају у квантној механици. Научна Математичке науке област: Научна Функционална анализа дисциплина: Силвестерова једначина; операторске једначине; алгебре оператора; спектрална теорија Кључне оператора; затворени оператори; Фредхолмова теорија; Банахове алгебре. речи: УДК: 517.983.23; 517.983.24; 517.983.5; 517.986.3; 517.984; 517.984.3; 517.984.4 CERIF P140: Класе, Фуријеова анализа, функционална анализа класифи- кација: Тип CC BY-NC-ND лиценце креативне заједнице: Podaci o mentoru i ˇclanovima komisije Mentor: dr NebojˇsaDinˇci´c Redovni profesor Prirodno-matematiˇckifakultet Univerzitet u Niˇsu Clanoviˇ komisije: 1. dr Vladimir Rakoˇcevi´c Dopisni ˇclanSANU Redovni profesor Prirodno-matematiˇckifakultet Univerzitet u Niˇsu 2. Akademik dr Stevan Pilipovi´c Redovni profesor Prirodno-matematiˇckifakultet Univerzitet u Novom Sadu 3. dr Peter Semrlˇ Redovni profesor Fakultet za matematiku i fiziku Univerzitet u Ljubljani 4. dr Sneˇzana Zivkovi´c-Zlatanovi´cˇ Redovni profesor Prirodno-matematiˇckifakultet Univerzitet u Niˇsu 5. dr Dijana Mosi´c Redovni profesor Prirodno-matematiˇckifakultet Univerzitet u Niˇsu i Abstract The main goal of this doctoral dissertation is to investigate behavior of sin- gular Sylvester equations, i. e. behavior of operator equations AX − XB = C; under the assumption that they are either unsolvable or have infinitely many solutions. Once solvability conditions are derived, characterization and clas- sification of the solutions is conducted and an explicit general formula for those solutions is provided, thus forming the general solution of the given Sylvester equation. Standard techniques, such as the generalized inverses, are omitted, because the assumption that A and B have closed ranges which are complemented in the corresponding Banach spaces is dropped. Instead, new and original methods are developed for solving this problem, and they are the original scientific contribution of the author, published in papers [24]{[29]. The dissertation is broken down into several chapters. Chapter 1 is the in- troductory chapter, where regular Sylvester equation (which has a unique solution) is introduced and solved. Some important applications of the equa- tion are mentioned. Chapter 2 concerns the singular case where A, B and C are matrices. The re- sults are obtained by the shared-eigenvalue discussion for matrices A and B, and by the analysis of the corresponding eigenspaces. Generalized commu- tators of matrices A and B are characterized, and the solutions are approxi- mated when possible. Perturbation analysis is conducted, using majorization theory for matrices. The main results in this chapter were obtained by the author and his PhD advisor in their joint works [28] and partially [29], and by the author in his individual paper [27]. Chapter 3 concerns the singular case when A, B and C are bounded linear operators on Banach spaces. Since the spectra of A and B do not necessarily consist of eigenvalues only, an alternative approach is required. First, a special operator algebra is introduced, which is not a Banach algebra per se, but still allows a functional calculus of its elements. This algebra gives a different form of the general solution to the given Sylvester equation, and solves every basic operator equation AX − XB = C; AXB = C; X − AXB = C; ii in the same manner, discarding regularity of the equations (only their solv- ability is required). The advantage of this method compared to the gener- alized inverses techniques (which are commonly used in singular equations) is that it does not require complementedness of the appropriate ranges and null-spaces, but rather solves each equation directly. This algebra is intro- duced and studied in detail by the author in [25]. Afterwards, this algebra is used to solve the initial Sylvester equation, with help from Fredholm theory. The author obtained these results in [26]. Applications to some contempo- rary problems in Banach spaces are illustrated as well. Chapter 4 concerns the singular case when A and B are densely defined closed operators on Banach spaces, and C is a densely defined linear operator. The initial premise is that the point spectra of A and B intersect, and in that case, weak solutions X are obtained, which are defined on appropriate eigenspaces of B. Techniques used involve decompositions of the given operators and spaces. Further, the results are extended to Schauder bases when possible, and are applied to Sturm-Liouville operators. These results were achieved by the author and his PhD advisor in the joint work [29]. Afterwards, a special case is analyzed, where A and B are densely defined self-adjoint operators on Hilbert spaces, while the point-spectrum-intersection assumption is dropped. In that case, the weak solutions X obtained in [29] are extended to the largest domains possible, which are constructed by the Spectral mapping theorem for self-adjoint operators and by the Berberian-Buoni-Harte-Wickstead con- struction. These results were obtained by the author in [24] and they are illustrated on an example which stems from quantum mechanics. iii Abstrakt Glavni cilj ove doktorske disertacije je ispitivanje prirode singularne Silves- terove jednaˇcine,odnosno, operatorske jednaˇcineoblika AX − XB = C; pod pretpostavkom da je ona ili nereˇsiva, ili da ima beskonaˇcnomnogo reˇsenja. Najpre bi se obezbedili dovoljni uslovi reˇsivosti jednaˇcine,a po- tom bi se sprovela karakterizacija i klasifikacija reˇsenja. Ta reˇsenjabi se zatim izvela analitˇckimi egzaktnim metodima, u zatvorenom obliku, ˇcimebi formirala opˇstereˇsenjepolazne jednaˇcine. Za razliku od standarndih metoda koriˇs´cenihza reˇsavanje singularnih op- eratorskih jednaˇcina,poput uopˇstenihinverza, u ovoj disertaciji se, izmed¯u ostalog, posmatraju i sluˇcajevi u kojima dati operatori nisu uopˇstenoin- vertibilni, odnonso, njihova jezgra i njihove slike ne moraju biti zatvoreni sa topoloˇskimkomplementima u odgovaraju´cimprostorima. Stoga se dati problem analizira na nov i originalan naˇcin,ˇstoje ujedno i nauˇcnidoprinos autora ovoj temi. Originalni rezultati autora, na kojima se i bazira ova dis- ertacija, publikovani su u radovima [24]{[29]. Sama disertacija je podeljena u nekoliko glava. Glava 1 je uvodnog karaktera,
Recommended publications
  • How and Why to Solve the Operator Equation Axâ‹™Xb
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Publications of the IAS Fellows HOW AND WHY TO SOLVE THE OPERATOR EQUATION AX®XB ¯ Y RAJENDRA BHATIA PETER ROSENTHAL 1. Introduction The entities A, B, X, Y in the title are operators, by which we mean either linear transformations on a finite-dimensional vector space (matrices) or bounded (¯ continuous) linear transformations on a Banach space. (All scalars will be complex numbers.) The definitions and statements below are valid in both the finite- dimensional and the infinite-dimensional cases, unless the contrary is stated. The simplest kind of equation involving operators is a linear equation of the form AX ¯ Y or of the form XB ¯ Y. The condition that is both necessary and sufficient that such an equation have a unique solution is that A (or B) should be inertible (bijective). Let σ(A) denote the spectrum of A; in the finite-dimensional case this is just the set of eigenvalues of A. In any case, invertibility of A is the statement 0 a σ(A). The equation ax ¯ y in numbers can always be solved when a 1 0. The analogous condition for the operator equation AX ¯ Y is 0 a σ(A). This analogy provides guidance in formulating potential theorems and guessing solutions of AX®XB ¯ Y. When does AX®XB ¯ Y have a unique solution X for each given Y? In the scalar case of ax®xb ¯ y, the answer is obvious: we must have a®b 1 0. The answer is almost as obvious in another case: if the matrices A and B are both diagonal, with diagonal entries ²λ",...,λn´ and ²µ",...,µn´, respectively, then the above equation has the solution xij ¯ yij}(λi®µj) provided λi®µj 1 0 for all i, j.
    [Show full text]
  • Computational Methods for Linear Matrix Equations∗
    COMPUTATIONAL METHODS FOR LINEAR MATRIX EQUATIONS∗ V. SIMONCINI† Abstract. Given the square matrices A,B,D,E and the matrix C of conforming dimensions, we consider the linear matrix equation AXE + DXB = C in the unknown matrix X. Our aim is to provide an overview of the major algorithmic developments that have taken place in the past few decades in the numerical solution of this and of related problems, which are becoming a reliable numerical tool in the formulation and solution of advanced mathematical models in engineering and scientific computing. Key words. Sylvester equation, Lyapunov equation, Stein equation, multiple right-hand side, generalized matrix equations. Schur decomposition. Large scale computation. AMS subject classifications. 65F10, 65F30, 15A06 Contents 1. Introduction 2 2. Notationandpreliminarydefinitions 4 3. Applications 5 4. Continuous-timeSylvesterequation 10 4.1. StabilityoftheSylvesterequation 12 4.2. Small scale computation 13 4.3. Computation with large A and small B 15 4.4. Computation with large A and large B 19 4.4.1. Projectionmethods 21 4.4.2. ADI iteration 25 4.4.3. Datasparseandothermethods 27 5. Continuous-time Lyapunov equation 28 5.1. Small scale computation 29 5.2. Large scale computation 31 5.2.1. Projectionmethods 31 5.2.2. ADI methods 35 5.2.3. Spectral, Sparse format and other methods 39 6. ThediscreteSylvesterandSteinequations 42 7. Generalizedlinearequations 44 7.1. The Generalized Sylvester and Lyapunov equations 44 7.2. Bilinear, constrained, and other linear equations 46 7.3. Sylvester-like and Lyapunov-like equations 50 8. Software and high performance computation 52 9. Concluding remarks and future outlook 53 10.
    [Show full text]
  • Arxiv:1903.05360V1 [Math.NA]
    Relation between the T-congruence Sylvester equation and the generalized Sylvester equation Yuki Satake,∗ Masaya Oozawa∗, Tomohiro Sogabe∗ Yuto Miyatake,† Tomoya Kemmochi∗, Shao-Liang Zhang∗ Abstract The T-congruenceSylvester equation is the matrix equation AX +XT B = C, where A ∈ Rm×n, B ∈ Rn×m, and C ∈ Rm×m are given, and X ∈ Rn×m is to be determined. Recently, Oozawa et al. discovered a transformation that the matrix equation is equivalent to one of the well-studied matrix equations (the Lyapunov equation); however, the condition of the transformation seems to be too limited because matrices A and B are assumed to be square matrices (m = n). In this paper, two transformations are provided for rectangular matrices A and B. One of them is an extension of the result of Oozawa et al. for the case m ≥ n, and the other is a novel transformation for the case m ≤ n. 1 Introduction We consider the T-congruence Sylvester equation of the form AX + XTB = C, (1) where A ∈ Rm×n, B ∈ Rn×m, and C ∈ Rm×m are given, and X ∈ Rn×m is to be determined. The matrix equation (1) appears in the industrial and applied mathematical field: palindromic eigenvalue problems for studying the vibration of rail tracks under the excitation arising from high-speed trains [1, 2] and also in an abstract mathematical field: the tangent spaces of the Stiefel manifold [3]. In the past decade, the numerical solvers for solving equation (1) in view of direct methods and iterative methods have been well developed, (see e.g., [4, 5, 6]); however, the theoretical features of the matrix equation itself such as the uniqueness of the solution have not been explored extensively.
    [Show full text]
  • SHIFTED KRONECKER PRODUCT SYSTEMS ∗ 1. Introduction. Matrix Problems with Replicated Block Structure Abound in Signal and Imag
    SHIFTED KRONECKER PRODUCT SYSTEMS ∗ CARLA D. MORAVITZ MARTIN † AND CHARLES F. VAN LOAN ‡ Abstract. A fast method for solving a linear system of the form (A(p) ⊗···⊗ A(1) − λI)x = b (i) is given where each A is an ni-by-ni matrix. The first step is to convert the problem to triangular form (T (p) ⊗···⊗ T (1) − λI)y = c by computing the (complex) Schur decompositions of the A(i). This is followed by a recursive back-substitution process that fully exploits the Kronecker structure and requires just O(N(n1 +···+np)) flops where N = n1 ··· np. A similar method is employed when the real Schur Decomposition is used to convert each A(i) to quasi-triangular form. The numerical properties of these new methods are the same as if we explicitly formed (T (p) ⊗···⊗ T (1) − λI) and used conventional back-substitution to solve for y. Key words. Linear Systems, Schur decomposition, Back-substitution, Kronecker products AMS subject classifications. 15A06, 65F05, 65G50 1. Introduction. Matrix problems with replicated block structure abound in signal and image processing, semidefinite programming, control theory, and many other application areas. In these venues fast algorithms have emerged that exploit the rich algebra of the Kronecker product. Perhaps the best example of this is the fast Fourier transform which can be described using the “language” of sparse, matrix factorizations and the Kronecker product. This operation is surfacing more and more as cheap memory prompts the assembly of huge, multidimensional datasets. When techniques for problems of low dimension are generalized or “tensored” together to address a high-dimensional, multilinear problem, then one typically finds a computa- tional challenge that involves the Kronecker product.
    [Show full text]
  • Graph Minor from Wikipedia, the Free Encyclopedia Contents
    Graph minor From Wikipedia, the free encyclopedia Contents 1 2 × 2 real matrices 1 1.1 Profile ................................................. 1 1.2 Equi-areal mapping .......................................... 2 1.3 Functions of 2 × 2 real matrices .................................... 2 1.4 2 × 2 real matrices as complex numbers ............................... 3 1.5 References ............................................... 4 2 Abelian group 5 2.1 Definition ............................................... 5 2.2 Facts ................................................. 5 2.2.1 Notation ........................................... 5 2.2.2 Multiplication table ...................................... 6 2.3 Examples ............................................... 6 2.4 Historical remarks .......................................... 6 2.5 Properties ............................................... 6 2.6 Finite abelian groups ......................................... 7 2.6.1 Classification ......................................... 7 2.6.2 Automorphisms ....................................... 7 2.7 Infinite abelian groups ........................................ 8 2.7.1 Torsion groups ........................................ 9 2.7.2 Torsion-free and mixed groups ................................ 9 2.7.3 Invariants and classification .................................. 9 2.7.4 Additive groups of rings ................................... 9 2.8 Relation to other mathematical topics ................................. 10 2.9 A note on the typography ......................................
    [Show full text]
  • Al-Kashi Background Statistics
    Al-Kashi Background Statistics http://www.ar-php.org/stats/al-kashi/ PDF generated using the open source mwlib toolkit. See http://code.pediapress.com/ for more information. PDF generated at: Sat, 01 Feb 2014 07:11:01 UTC Contents Articles Summary Statistics 1 Mean 1 Median 6 Mode (statistics) 15 Variance 20 Standard deviation 32 Coefficient of variation 46 Skewness 49 Kurtosis 55 Ranking 61 Graphics 65 Box plot 65 Histogram 68 Q–Q plot 74 Ternary plot 79 Distributions 84 Normal distribution 84 Student's t-distribution 113 F-distribution 128 Feature scaling 131 Correlation and Regression 134 Covariance 134 Correlation and dependence 138 Regression analysis 144 Path analysis (statistics) 154 Analysis 156 Moving average 156 Student's t-test 162 Contingency table 171 Analysis of variance 174 Principal component analysis 189 Diversity index 205 Diversity index 205 Clustering 210 Hierarchical clustering 210 K-means clustering 215 Matrix 225 Matrix (mathematics) 225 Matrix addition 247 Matrix multiplication 249 Transpose 263 Determinant 266 Minor (linear algebra) 284 Adjugate matrix 287 Invertible matrix 291 Eigenvalues and eigenvectors 298 System of linear equations 316 References Article Sources and Contributors 326 Image Sources, Licenses and Contributors 332 Article Licenses License 335 1 Summary Statistics Mean In mathematics, mean has several different definitions depending on the context. In probability and statistics, mean and expected value are used synonymously to refer to one measure of the central tendency either of a probability distribution or of the random variable characterized by that distribution. In the case of a discrete probability distribution of a random variable X, the mean is equal to the sum over every possible value weighted by the probability of that value; that is, it is computed by taking the product of each possible value x of X and its probability P(x), and then adding all these products together, giving .[1] An analogous formula applies to the case of a continuous probability distribution.
    [Show full text]
  • Sylvester Equations M×N P×Q Goal Represent Linear Functions R → R
    Sylvester equations m×n p×q Goal represent linear functions R → R . For instance, to deal with problems like the following one. Sylvester equation AX − XB = C m×m m×n n×n A ∈ C , C, X ∈ C , B ∈ C . This must be a mn × mn linear system, right? Vectorization gives an explicit way to map it to a vector. Vectorization: definition x11 x 21 . . xm1 x12 x x11 x12 ... x1n 22 . x21 x22 ... x2n . vec X = vec . := . . .. xm2 x x ... x . m1 m2 mn . x1n x2n . . xmn Vectorization: comments Column-major order: leftmost index ‘changes more often’. Matches Fortran, Matlab standard (C/C++ prefer row-major instead). Converting indices in the matrix into indices in the vector: (X)ij = (vec X)i+mj 0-based, (X)ij = (vec X)i+m(j−1) 1-based. vec(AXB) First, we will work out the representation of a simple linear map, X 7→ AXB (for fixed matrices A, B of compatible dimensions). m×n p×q If X ∈ R , AXB ∈ R , we need the pq × mn matrix that maps vec X to vec(AXB). X X X (AXB)hl = (AX)hj (B)jl = Ahi Xij Bjl j j i = Ah1B1l Ah2B1l ... AhmB1l Ah1B2l Ah2B2l ... AhmB2l ... Ah1Bnl Ah2Bnl AhmBnl ] vec X Kronecker product: definition b11A b21A ... bn1A b12A b22A ... bn2A vec(AXB) = . vec X . .. . b1qA b2qA ... bnqA Each block is a multiple of A, with coefficient given by the corresponding entry of B>. Definition x11Y x12Y ... x1nY . . x21Y x22Y ... X ⊗ Y := . . . .. .. . xm1Y xm2Y ... xmnY so the matrix above is B> ⊗ A. Properties of Kronecker products x11Y x12Y ..
    [Show full text]
  • On the Solution of the Nonsymmetric T-Riccati Equation∗
    Electronic Transactions on Numerical Analysis. Volume 54, pp. 68–88, 2021. ETNA Kent State University and Copyright © 2021, Kent State University. Johann Radon Institute (RICAM) ISSN 1068–9613. DOI: 10.1553/etna_vol54s68 ON THE SOLUTION OF THE NONSYMMETRIC T-RICCATI EQUATION∗ PETER BENNERy AND DAVIDE PALITTAy Abstract. The nonsymmetric T-Riccati equation is a quadratic matrix equation where the linear part corresponds to the so-called T-Sylvester or T-Lyapunov operator that has previously been studied in the literature. It has applications in macroeconomics and policy dynamics. So far, it presents an unexplored problem in numerical analysis, and both theoretical results and computational methods are lacking in the literature. In this paper we provide some sufficient conditions for the existence and uniqueness of a nonnegative minimal solution, namely the solution with component- wise minimal entries. Moreover, the efficient computation of such a solution is analyzed. Both the small-scale and large-scale settings are addressed, and Newton-Kleinman-like methods are derived. The convergence of these procedures to the minimal solution is proven, and several numerical results illustrate the computational efficiency of the proposed methods. Key words. T-Riccati equation, M-matrices, minimal nonnegative solution, Newton-Kleinman method AMS subject classifications. 65F30, 15A24, 49M15, 39B42, 40C05 1. Introduction. In this paper, we consider the nonsymmetric T-Riccati operator n×n n×n T T RT : R ! R ; RT (X) := DX + X A − X BX + C; where A; B; C; D 2 Rn×n, and we provide sufficient conditions for the existence and unique- n×n ness of a minimal solution Xmin 2 R to (1.1) RT (X) = 0: The solution of the nonsymmetric T-Riccati equation (1.1) plays a role in solving dynamics generalized equilibrium (DSGE) problems [9, 22, 25].
    [Show full text]
  • Solving the Sylvester Equation Ax − Xb = C When Σ(A) ∩ Σ(B) 6= ∅ ∗
    Electronic Journal of Linear Algebra, ISSN 1081-3810 A publication of the International Linear Algebra Society Volume 35, pp. 1-23, January 2019. SOLVING THE SYLVESTER EQUATION AX − XB = C WHEN σ(A) \ σ(B) 6= ; ∗ NEBOJSAˇ C.ˇ DINCIˇ C´ y Abstract. The method for solving the Sylvester equation AX−XB = C in the complex matrix case, when σ(A)\σ(B) 6= ;, by using Jordan normal form, is given. Also, the approach via the Schur decomposition is presented. Key words. Sylvester equation, Jordan normal form, Schur decomposition. AMS subject classifications. 15A24. 1. An introduction. The Sylvester equation (1.1) AX − XB = C; where A 2 Cm×m;B 2 Cn×n and C 2 Cm×n are given, and X 2 Cm×n is an unknown matrix, is a very popular topic in linear algebra, with numerous applications in control theory, signal processing, image restoration, engineering, solving ordinary and partial differential equations, etc. It is named for famous mathematician J. J. Sylvester, who was the first to prove that this equation in matrix case has unique solution if and only if σ(A) \ σ(B) = ; [17]. One important special case of the Sylvester equation is the continuous-time Lyapunov matrix equation (AX + XA∗ = C). AC In 1952, Roth [15] proved that for operators A; B on finite-dimensional spaces is similar to 0 B A 0 if and only if the equation AX − XB = C has a solution X. Rosenblum [14] proved that if A 0 B and B are operators such that σ(A) \ σ(B) = ;, then the equation AX − XB = Y has a unique solution X for every operator Y .
    [Show full text]
  • Solution Formulas for Differential Sylvester and Lyapunov Equations
    Solution Formulas for Differential Sylvester and Lyapunov Equations Maximilian Behr Peter Benner Jan Heiland November 21, 2018 Abstract The differential Sylvester equation and its symmetric version, the differential Lyapunov equation, appear in different fields of applied mathematics like control theory, system theory, and model order reduction. The few available straight-forward numerical approaches if applied to large- scale systems come with prohibitively large storage requirements. This shortage motivates us to summarize and explore existing solution formulas for these equations. We develop a unifying approach based on the spectral theorem for normal operators like the Sylvester operator (X) = AX + XB and derive a formula for its norm using an induced operator norm based on S the spectrum of A and B. In view of numerical approximations, we propose an algorithm that identifies a suitable Krylov subspace using Taylor series and use a projection to approximate the solution. Numerical results for large-scale differential Lyapunov equations are presented in the last sections. Contents 1. Introduction 2 2. Preliminaries 3 3. Spectral Decomposition of the Sylvester Operator4 4. Variation of Constants7 5. Solution as Taylor Series7 6. Feasible Numerical Solution Approaches9 arXiv:1811.08327v1 [math.NA] 20 Nov 2018 7. Numerical Results 12 7.1. Setup . 12 7.2. Projection Approach . 12 7.3. Computational Time . 15 7.4. Comparison with Backward Differentiation Formulas . 16 8. Conclusions 17 A. Numerical Results Projection Approach 18 A.1. Results for the Differential Lyapunov Equation . 18 A.2. Results for the Transposed Differential Lyapunov Equation . 22 1 B. Backward Differentiation Formulas 25 B.1. Results for the Differential Lyapunov Equation .
    [Show full text]
  • Krylov Methods for Low-Rank Commuting Generalized Sylvester Equations
    Krylov methods for low-rank commuting generalized Sylvester equations Elias Jarlebring ,∗ Giampaolo Mele ,∗ Davide Palitta ,y Emil Ringh ∗ April 14, 2017 Abstract We consider generalizations of the Sylvester matrix equation, consist- ing of the sum of a Sylvester operator and a linear operator Π with a particular structure. More precisely, the commutator of the matrix coef- ficients of the operator Π and the Sylvester operator coefficients are as- sumed to be matrices with low rank. We show (under certain additional conditions) low-rank approximability of this problem, i.e., the solution to this matrix equation can be approximated with a low-rank matrix. Projec- tion methods have successfully been used to solve other matrix equations with low-rank approximability. We propose a new projection method for this class of matrix equations. The choice of subspace is a crucial ingredi- ent for any projection method for matrix equations. Our method is based on an adaption and extension of the extended Krylov subspace method for Sylvester equations. A constructive choice of the starting vector/block is derived from the low-rank commutators. We illustrate the effectiveness of our method by solving large-scale matrix equations arising from appli- cations in control theory and the discretization of PDEs. The advantages of our approach in comparison to other methods are also illustrated. Keywords Generalized Sylvester equation · Low-rank commutation · Krylov subspace · projection methods · Iterative solvers · Matrix equation Mathematics Subject Classification
    [Show full text]
  • Affine Transformation
    Affine transformation From Wikipedia, the free encyclopedia Contents 1 2 × 2 real matrices 1 1.1 Profile ................................................. 1 1.2 Equi-areal mapping .......................................... 2 1.3 Functions of 2 × 2 real matrices .................................... 2 1.4 2 × 2 real matrices as complex numbers ............................... 3 1.5 References ............................................... 4 2 3D projection 5 2.1 Orthographic projection ........................................ 5 2.2 Weak perspective projection ..................................... 5 2.3 Perspective projection ......................................... 6 2.4 Diagram ................................................ 8 2.5 See also ................................................ 8 2.6 References ............................................... 9 2.7 External links ............................................. 9 2.8 Further reading ............................................ 9 3 Affine coordinate system 10 3.1 See also ................................................ 10 4 Affine geometry 11 4.1 History ................................................. 12 4.2 Systems of axioms ........................................... 12 4.2.1 Pappus’ law .......................................... 12 4.2.2 Ordered structure ....................................... 13 4.2.3 Ternary rings ......................................... 13 4.3 Affine transformations ......................................... 14 4.4 Affine space .............................................
    [Show full text]