Technische Universit at Chemnitz

Total Page:16

File Type:pdf, Size:1020Kb

Technische Universit at Chemnitz Technische Universitat Chemnitz Sonderforschungsbereich Numerische Simulation auf massiv parallelen Rechnern Peter Benner Ralph Byers Volker Mehrmann Hongguo Xu Numerical Computation of Deating Subspaces of Embedded Hamiltonian Pencils Preprint SFB PreprintReihe des Chemnitzer SFB SFB June Contents Introduction and Preliminaries Embedding of Extended Matrix Pencils Hamiltonian Triangular Forms SkewHamiltonianHamiltonian Matrix Pencils Algorithms Error and Perturbation Analysis Conclusion A Algorithmic Details A Schur Forms for SkewHamiltonianHamiltonian Matrix Pencils and SkewHa miltonianskewHamiltonian Matrix Pencils A Subroutines Required by Algorithm Authors Peter Benner Ralph Byers Zentrum f ur Technomathematik Department of Mathematics Fachbereich Mathematik und Informatik University of Kansas Universitat Bremen Lawrence KS D Bremen FRG USA bennermathunibremende byersmathukansedu Volker Mehrmann Hongguo Xu Deparment of Mathematics Fakultat f ur Mathematik Case Western University Technische Universitat Chemnitz Cleveland OH D Chemnitz FRG USA mehrmannmathematiktuchemnitzde hxxpocwruedu Numerical Computation of Deating Subspaces of Embedded Hamiltonian Pencils 2 3 Peter Benner Ralph Byers Volker Mehrmann Hongguo Xu Abstract We discuss the numerical solution of structured generalized eigenvalue problems that arise from linearquadratic optimal control problems H optimization multibo dy sys 1 tems and many other areas of applied mathematics physics and chemistry The classical approach for these problems requires computing invariant and deating subspaces of ma trices and matrix p encils with Hamiltonian andor skewHamiltonian structure We ex tend the recently developed metho ds for Hamiltonian matrices and matrix p encils to the general case of embedded matrix p encils The rounding error and p erturbation analysis of the resulting algorithms is favorable Keywords eigenvalue problem deating subspace algebraic Riccati equation Hamiltonian matrix skewHamiltonian matrix skewHamiltonianHamiltonian matrix p encil embedded matrix p encils AMS sub ject classication N F B B C Introduction and Preliminaries In this pap er we study eigenvalue and invariant subspace computations involving matrices and matrix p encils with the following algebraic structures i h I 0 n where I is the n n identity matrix Denition Let J n I 0 n 2n2n H a A matrix H C is Hamiltonian if HJ HJ The Lie Algebra of Hamiltonian 2 n2n matrices in C is denoted by H 2n 2n2n H b A matrix H C is skewHamiltonian if HJ H J The Jordan algebra of 2 n2n skewHamiltonian matrices in C is denoted by SH 2n 2 n2n c A matrix pencil S H C is skewHamiltonianHamiltonian if S SH and 2n H H 2n 2 n2n H d A matrix S C is symplectic if SJ S J The Lie group of symplectic matrices 2 n2n in C is denoted by S 2n H 2n2n H I The J and U U e A matrix U C is unitary symplectic if U J U 2n d d d d d 2 n2n compact Lie group of unitary symplectic matrices in C is denoted by US 2n 2n H f A subspace L of C is called Lagrangian if it has dimension n and x J y for al l x y L Working title was Numerical Computation of Deating Subspaces of Embedded Hamiltonian and Sym plectic Pencils All authors were partially supp orted by Deutsche Forschungsgemeinschaft Research Grant Me and Sonderforschungsbereich Numerische Simulation auf massiv parallelen Rechnern 2 This author was partially supp orted by National Science Foundation awards CCR MRI and by the NSF EPSCoRKSTAR program through the Center for Advanced Scientic Computing 3 This work was completed while this author was with the TU Chemnitz There is little dierence b etween the structure of complex skewHamiltonian matrices and complex Hamiltonian matrices If S SH and H H then iS H and iH 2n 2 n 2n SH Similarly there is little dierence b etween the structure of complex skewHamiltoni 2n anHamiltonian matrix p encils complex skewHamiltonianskewHamiltonian and complex HamiltonianHamiltonian matrix p encils For S SH and H H the matrix p encil 2n 2 n S H is skewHamiltonianHamiltonian the matrix p encil S iH is skewHamil tonianskewHamiltonian and the matrix p encil iS H is HamiltonianHamiltonian Here i and i However real skewHamiltonian matrices are not real scalar multiples of Hamiltonian matrices There is a greater dierence in the structure of real skewHamiltonian and real Hamiltonian matrices than in the complex case The structures in Denition arise in linearquadratic optimal control H optimization and several other areas of applied mathematics computational physics and chemistry eg gyroscopic systems numerical simulation of elastic deformation and linear resp onse theory Here we fo cus on applications in linearquadratic optimal control and H optimization First we consider the continuous time innite horizon linearquadratic optimal control problem Choose a control function ut to minimize the cost functional Z H xt Q S xt S dt c H ut S R ut t 0 subject to the linear dierentialalgebraic system descriptor control system 0 E x Ax B u xt x 0 m n nn nm H nn H mm Here ut C xt C A E C B C Q Q C R R C and nm S C If the m n m n weighting matrix Q S R H S R is Hermitian and p ositive semidenite then the problem is wellposed Note that R Q andor R may b e singular Typically in addition to minimizing the control ut must make xt asymptotically stable Of course if R is p ositive denite then asymptotic stability of xt is automatically achieved by any control for which the cost functional is nite Application of the maximum principle yields as a necessary optimality condition that the control u satises the twopoint b oundary value problem of EulerLagrange equations x x 0 H E A xt x lim E t c c 0 t u u with the matrix p encil A B E H H Q A S E E A c c H H S B R The matrix p encil do es not have matrices with one of the structures of Denition Nevertheless many of the prop erties of Hamiltonian matrices carry over As reviewed b e low if E andor R are nonsingular then and reduce to an ordinary dierential equation with Hamiltonian matrix co ecients or a dierentialalgebraic equation with a skewHamilto nianHamiltonian matrix p encil co ecients In Section we circumvent the nonsingularity assumptions by embedding the matrix p encil into a skewHamiltonianHamiltonian ma trix p encil of larger dimension H If b oth E and R are nonsingular then with E reduces to the twopoint b oundary value problem x x 0 H xt x lim t 0 t with the Hamiltonian matrix 1 1 H 1 1 H H F G E A BR S E BR B E H 1 H 1 1 H H H Q SR S E A BR S H F This is the classical formulation found in many textb o oks on linearquadratic optimal control like The solution of this b oundary value problem can b e obtained in many dierent ways For example let Y b e a symmetric solution if it exists of the asso ciated continuoustime algebraic Riccati equation H H Y F F Y Y GY h i h i h i h i I 0 x I 0 x Multiplying from the left by the matrix and changing variables to Y I Y I one obtains the decoupled system x F GY G x xt x lim t 0 0 H t F Y G In the desired solution is identically zero In that case x is the solution to x F GY x 0 H with initial condition x x Y x E and the control that minimizes 1 H H H is u R S B E Y x If E is singular and R is nonsingular then do es not simplify quite so much b ecause it is a dierentialalgebraic equation with nontrivial linear constraints Substituting ut 1 H H R S xt B t system do es simplify to x x 0 H S H xt x lim E t 0 t with the reduced skewHamiltonianHamiltonian matrix p encil 1 H 1 H E A BR S BR B S H H 1 H 1 H E SR S Q A BR S In this case the corresp onding generalized algebraic Riccati equation is 1 H H 1 Q SR S E Y A BR S 1 H H 1 H A BR S Y E E Y BR B Y E but in general the relationship with solutions of the optimal control problem is lost or hidden See for details But even for E nonsingular the formulation in is preferable from a numerical p oint of view if E is illconditioned with resp ect to inversion If R is singular and E is nonsingular then the situation b ecomes still more complicated Although the b oundary value problem remains well dened the Riccati equation do es not This case has b een recently studied At this writing the case in which b oth E and R are singular has not b een analyzed in full generality The EulerLagrange and Riccati equations their solvability and their numerical solution have b een the sub ject of numerous publications in recent years see eg and the references therein In most numerical metho ds the Riccati solutions are obtained through the computation of deating or invariant subspaces of asso ciated matrix p encils see A key observation in this context is that the Riccati solutions need not b e formed explicitly It suces to work with bases of the deating subspaces of the matrix p encil For example supp ose E A in has an ndimensional deating subspace asso ciated c c with eigenvalues in the left half plane This can only b e the case if E is invertible Let this subspace b e spanned by the columns of a matrix U partitioned conformally with as U 1 U U 2 U 3 1 If U is invertible then the optimal control is a linear feedback of the form ut U U xt 1 3 1 1 1 and the solution of the asso ciated Riccati equation is Y U U E see Observe 2 1 that the optimal control may
Recommended publications
  • Orthogonal Reduction 1 the Row Echelon Form -.: Mathematical
    MATH 5330: Computational Methods of Linear Algebra Lecture 9: Orthogonal Reduction Xianyi Zeng Department of Mathematical Sciences, UTEP 1 The Row Echelon Form Our target is to solve the normal equation: AtAx = Atb ; (1.1) m×n where A 2 R is arbitrary; we have shown previously that this is equivalent to the least squares problem: min jjAx−bjj : (1.2) x2Rn t n×n As A A2R is symmetric positive semi-definite, we can try to compute the Cholesky decom- t t n×n position such that A A = L L for some lower-triangular matrix L 2 R . One problem with this approach is that we're not fully exploring our information, particularly in Cholesky decomposition we treat AtA as a single entity in ignorance of the information about A itself. t m×m In particular, the structure A A motivates us to study a factorization A=QE, where Q2R m×n is orthogonal and E 2 R is to be determined. Then we may transform the normal equation to: EtEx = EtQtb ; (1.3) t m×m where the identity Q Q = Im (the identity matrix in R ) is used. This normal equation is equivalent to the least squares problem with E: t min Ex−Q b : (1.4) x2Rn Because orthogonal transformation preserves the L2-norm, (1.2) and (1.4) are equivalent to each n other. Indeed, for any x 2 R : jjAx−bjj2 = (b−Ax)t(b−Ax) = (b−QEx)t(b−QEx) = [Q(Qtb−Ex)]t[Q(Qtb−Ex)] t t t t t t t t 2 = (Q b−Ex) Q Q(Q b−Ex) = (Q b−Ex) (Q b−Ex) = Ex−Q b : Hence the target is to find an E such that (1.3) is easier to solve.
    [Show full text]
  • 1 Introduction
    Real and Complex Hamiltonian Square Ro ots of SkewHamiltonian Matrices y z HeikeFab ender D Steven Mackey Niloufer Mackey x Hongguo Xu January DedicatedtoProfessor Ludwig Elsner on the occasion of his th birthday Abstract We present a constructive existence pro of that every real skewHamiltonian matrix W has a real Hamiltonian square ro ot The key step in this construction shows how one may bring anysuch W into a real quasiJordan canonical form via symplectic similarityWe show further that every W has innitely many real Hamiltonian square ro ots and givealower b ound on the dimension of the set of all such square ro ots Extensions to complex matrices are also presented This report is an updated version of the paper Hamiltonian squareroots of skew Hamiltonian matrices that appearedinLinear Algebra its Applicationsv pp AMS sub ject classication A A A F B Intro duction Any matrix X suchthatX A is said to b e a square ro ot of the matrix AF or general nn complex matrices A C there exists a welldevelop ed although somewhat complicated theory of matrix square ro ots and a numb er of algorithms for their eective computation Similarly for the theory and computation of real square ro ots for real matrices By contrast structured squareroot problems where b oth the matrix A Fachb ereich Mathematik und Informatik Zentrum fur Technomathematik Universitat Bremen D Bremen FRG email heikemathunibremende y Department of Mathematics Computer Science Kalamazo o College Academy Street Kalamazo o MI USA z Department of Mathematics
    [Show full text]
  • Eigenvalue Perturbation Theory of Structured Real Matrices and Their Sign Characteristics Under Generic Structured Rank-One Perturbations∗
    Eigenvalue perturbation theory of structured real matrices and their sign characteristics under generic structured rank-one perturbations∗ Christian Mehlx Volker Mehrmannx Andr´eC. M. Ran{ Leiba Rodmank July 9, 2014 Abstract An eigenvalue perturbation theory under rank-one perturbations is developed for classes of real matrices that are symmetric with respect to a non-degenerate bilinear form, or Hamiltonian with respect to a non-degenerate skew-symmetric form. In contrast to the case of complex matrices, the sign characteristic is a cru- cial feature of matrices in these classes. The behavior of the sign characteristic under generic rank-one perturbations is analyzed in each of these two classes of matrices. Partial results are presented, but some questions remain open. Appli- cations include boundedness and robust boundedness for solutions of structured systems of linear differential equations with respect to general perturbations as well as with respect to structured rank perturbations of the coefficients. Key Words: real J-Hamiltonian matrices, real H-symmetric matrices, indefinite inner product, perturbation analysis, generic perturbation, rank-one perturbation, T - even matrix polynomial, symmetric matrix polynomial, bounded solution of differential equations, robustly bounded solution of differential equations, invariant Lagrangian subspaces. Mathematics Subject Classification: 15A63, 15A21, 15A54, 15B57. xInstitut f¨urMathematik, TU Berlin, Straße des 17. Juni 136, D-10623 Berlin, Germany. Email: fmehl,[email protected]. {Afdeling Wiskunde, Faculteit der Exacte Wetenschappen, Vrije Universiteit Amsterdam, De Boele- laan 1081a, 1081 HV Amsterdam, The Netherlands and Unit for BMI, North West University, Potchef- stroom, South Africa. E-mail: [email protected] kCollege of William and Mary, Department of Mathematics, P.O.Box 8795, Williamsburg, VA 23187-8795, USA.
    [Show full text]
  • On the Eigenvalues of Euclidean Distance Matrices
    “main” — 2008/10/13 — 23:12 — page 237 — #1 Volume 27, N. 3, pp. 237–250, 2008 Copyright © 2008 SBMAC ISSN 0101-8205 www.scielo.br/cam On the eigenvalues of Euclidean distance matrices A.Y. ALFAKIH∗ Department of Mathematics and Statistics University of Windsor, Windsor, Ontario N9B 3P4, Canada E-mail: [email protected] Abstract. In this paper, the notion of equitable partitions (EP) is used to study the eigenvalues of Euclidean distance matrices (EDMs). In particular, EP is used to obtain the characteristic poly- nomials of regular EDMs and non-spherical centrally symmetric EDMs. The paper also presents methods for constructing cospectral EDMs and EDMs with exactly three distinct eigenvalues. Mathematical subject classification: 51K05, 15A18, 05C50. Key words: Euclidean distance matrices, eigenvalues, equitable partitions, characteristic poly- nomial. 1 Introduction ( ) An n ×n nonzero matrix D = di j is called a Euclidean distance matrix (EDM) 1, 2,..., n r if there exist points p p p in some Euclidean space < such that i j 2 , ,..., , di j = ||p − p || for all i j = 1 n where || || denotes the Euclidean norm. i , ,..., Let p , i ∈ N = {1 2 n}, be the set of points that generate an EDM π π ( , ,..., ) D. An m-partition of D is an ordered sequence = N1 N2 Nm of ,..., nonempty disjoint subsets of N whose union is N. The subsets N1 Nm are called the cells of the partition. The n-partition of D where each cell consists #760/08. Received: 07/IV/08. Accepted: 17/VI/08. ∗Research supported by the Natural Sciences and Engineering Research Council of Canada and MITACS.
    [Show full text]
  • 4.1 RANK of a MATRIX Rank List Given Matrix M, the Following Are Equal
    page 1 of Section 4.1 CHAPTER 4 MATRICES CONTINUED SECTION 4.1 RANK OF A MATRIX rank list Given matrix M, the following are equal: (1) maximal number of ind cols (i.e., dim of the col space of M) (2) maximal number of ind rows (i.e., dim of the row space of M) (3) number of cols with pivots in the echelon form of M (4) number of nonzero rows in the echelon form of M You know that (1) = (3) and (2) = (4) from Section 3.1. To see that (3) = (4) just stare at some echelon forms. This one number (that all four things equal) is called the rank of M. As a special case, a zero matrix is said to have rank 0. how row ops affect rank Row ops don't change the rank because they don't change the max number of ind cols or rows. example 1 12-10 24-20 has rank 1 (maximally one ind col, by inspection) 48-40 000 []000 has rank 0 example 2 2540 0001 LetM= -2 1 -1 0 21170 To find the rank of M, use row ops R3=R1+R3 R4=-R1+R4 R2ØR3 R4=-R2+R4 2540 0630 to get the unreduced echelon form 0001 0000 Cols 1,2,4 have pivots. So the rank of M is 3. how the rank is limited by the size of the matrix IfAis7≈4then its rank is either 0 (if it's the zero matrix), 1, 2, 3 or 4. The rank can't be 5 or larger because there can't be 5 ind cols when there are only 4 cols to begin with.
    [Show full text]
  • Workshop on Applied Matrix Positivity International Centre for Mathematical Sciences 19Th – 23Rd July 2021
    Workshop on Applied Matrix Positivity International Centre for Mathematical Sciences 19th { 23rd July 2021 Timetable All times are in BST (GMT +1). 0900 1000 1100 1600 1700 1800 Monday Berg Bhatia Khare Charina Schilling Skalski Tuesday K¨ostler Franz Fagnola Apanasovich Emery Cuevas Wednesday Jain Choudhury Vishwakarma Pascoe Pereira - Thursday Sharma Vashisht Mishra - St¨ockler Knese Titles and abstracts 1. New classes of multivariate covariance functions Tatiyana Apanasovich (George Washington University, USA) The class which is refereed to as the Cauchy family allows for the simultaneous modeling of the long memory dependence and correlation at short and intermediate lags. We introduce a valid parametric family of cross-covariance functions for multivariate spatial random fields where each component has a covariance function from a Cauchy family. We present the conditions on the parameter space that result in valid models with varying degrees of complexity. Practical implementations, including reparameterizations to reflect the conditions on the parameter space will be discussed. We show results of various Monte Carlo simulation experiments to explore the performances of our approach in terms of estimation and cokriging. The application of the proposed multivariate Cauchy model is illustrated on a dataset from the field of Satellite Oceanography. 1 2. A unified view of covariance functions through Gelfand pairs Christian Berg (University of Copenhagen) In Geostatistics one examines measurements depending on the location on the earth and on time. This leads to Random Fields of stochastic variables Z(ξ; u) indexed by (ξ; u) 2 2 belonging to S × R, where S { the 2-dimensional unit sphere { is a model for the earth, and R is a model for time.
    [Show full text]
  • Doubling Algorithms with Permuted Lagrangian Graph Bases
    DOUBLING ALGORITHMS WITH PERMUTED LAGRANGIAN GRAPH BASES VOLKER MEHRMANN∗ AND FEDERICO POLONI† Abstract. We derive a new representation of Lagrangian subspaces in the form h I i Im ΠT , X where Π is a symplectic matrix which is the product of a permutation matrix and a real orthogonal diagonal matrix, and X satisfies 1 if i = j, |Xij | ≤ √ 2 if i 6= j. This representation allows to limit element growth in the context of doubling algorithms for the computation of Lagrangian subspaces and the solution of Riccati equations. It is shown that a simple doubling algorithm using this representation can reach full machine accuracy on a wide range of problems, obtaining invariant subspaces of the same quality as those computed by the state-of-the-art algorithms based on orthogonal transformations. The same idea carries over to representations of arbitrary subspaces and can be used for other types of structured pencils. Key words. Lagrangian subspace, optimal control, structure-preserving doubling algorithm, symplectic matrix, Hamiltonian matrix, matrix pencil, graph subspace AMS subject classifications. 65F30, 49N10 1. Introduction. A Lagrangian subspace U is an n-dimensional subspace of C2n such that u∗Jv = 0 for each u, v ∈ U. Here u∗ denotes the conjugate transpose of u, and the transpose of u in the real case, and we set 0 I J = . −I 0 The computation of Lagrangian invariant subspaces of Hamiltonian matrices of the form FG H = H −F ∗ with H = H∗,G = G∗, satisfying (HJ)∗ = HJ, (as well as symplectic matrices S, satisfying S∗JS = J), is an important task in many optimal control problems [16, 25, 31, 36].
    [Show full text]
  • Rotation Matrix - Wikipedia, the Free Encyclopedia Page 1 of 22
    Rotation matrix - Wikipedia, the free encyclopedia Page 1 of 22 Rotation matrix From Wikipedia, the free encyclopedia In linear algebra, a rotation matrix is a matrix that is used to perform a rotation in Euclidean space. For example the matrix rotates points in the xy -Cartesian plane counterclockwise through an angle θ about the origin of the Cartesian coordinate system. To perform the rotation, the position of each point must be represented by a column vector v, containing the coordinates of the point. A rotated vector is obtained by using the matrix multiplication Rv (see below for details). In two and three dimensions, rotation matrices are among the simplest algebraic descriptions of rotations, and are used extensively for computations in geometry, physics, and computer graphics. Though most applications involve rotations in two or three dimensions, rotation matrices can be defined for n-dimensional space. Rotation matrices are always square, with real entries. Algebraically, a rotation matrix in n-dimensions is a n × n special orthogonal matrix, i.e. an orthogonal matrix whose determinant is 1: . The set of all rotation matrices forms a group, known as the rotation group or the special orthogonal group. It is a subset of the orthogonal group, which includes reflections and consists of all orthogonal matrices with determinant 1 or -1, and of the special linear group, which includes all volume-preserving transformations and consists of matrices with determinant 1. Contents 1 Rotations in two dimensions 1.1 Non-standard orientation
    [Show full text]
  • §9.2 Orthogonal Matrices and Similarity Transformations
    n×n Thm: Suppose matrix Q 2 R is orthogonal. Then −1 T I Q is invertible with Q = Q . n T T I For any x; y 2 R ,(Q x) (Q y) = x y. n I For any x 2 R , kQ xk2 = kxk2. Ex 0 1 1 1 1 1 2 2 2 2 B C B 1 1 1 1 C B − 2 2 − 2 2 C B C T H = B C ; H H = I : B C B − 1 − 1 1 1 C B 2 2 2 2 C @ A 1 1 1 1 2 − 2 − 2 2 x9.2 Orthogonal Matrices and Similarity Transformations n×n Def: A matrix Q 2 R is said to be orthogonal if its columns n (1) (2) (n)o n q ; q ; ··· ; q form an orthonormal set in R . Ex 0 1 1 1 1 1 2 2 2 2 B C B 1 1 1 1 C B − 2 2 − 2 2 C B C T H = B C ; H H = I : B C B − 1 − 1 1 1 C B 2 2 2 2 C @ A 1 1 1 1 2 − 2 − 2 2 x9.2 Orthogonal Matrices and Similarity Transformations n×n Def: A matrix Q 2 R is said to be orthogonal if its columns n (1) (2) (n)o n q ; q ; ··· ; q form an orthonormal set in R . n×n Thm: Suppose matrix Q 2 R is orthogonal. Then −1 T I Q is invertible with Q = Q . n T T I For any x; y 2 R ,(Q x) (Q y) = x y.
    [Show full text]
  • The Jordan Canonical Forms of Complex Orthogonal and Skew-Symmetric Matrices Roger A
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Elsevier - Publisher Connector Linear Algebra and its Applications 302–303 (1999) 411–421 www.elsevier.com/locate/laa The Jordan Canonical Forms of complex orthogonal and skew-symmetric matrices Roger A. Horn a,∗, Dennis I. Merino b aDepartment of Mathematics, University of Utah, 155 South 1400 East, Salt Lake City, UT 84112–0090, USA bDepartment of Mathematics, Southeastern Louisiana University, Hammond, LA 70402-0687, USA Received 25 August 1999; accepted 3 September 1999 Submitted by B. Cain Dedicated to Hans Schneider Abstract We study the Jordan Canonical Forms of complex orthogonal and skew-symmetric matrices, and consider some related results. © 1999 Elsevier Science Inc. All rights reserved. Keywords: Canonical form; Complex orthogonal matrix; Complex skew-symmetric matrix 1. Introduction and notation Every square complex matrix A is similar to its transpose, AT ([2, Section 3.2.3] or [1, Chapter XI, Theorem 5]), and the similarity class of the n-by-n complex symmetric matrices is all of Mn [2, Theorem 4.4.9], the set of n-by-n complex matrices. However, other natural similarity classes of matrices are non-trivial and can be characterized by simple conditions involving the Jordan Canonical Form. For example, A is similar to its complex conjugate, A (and hence also to its T adjoint, A∗ = A ), if and only if A is similar to a real matrix [2, Theorem 4.1.7]; the Jordan Canonical Form of such a matrix can contain only Jordan blocks with real eigenvalues and pairs of Jordan blocks of the form Jk(λ) ⊕ Jk(λ) for non-real λ.We denote by Jk(λ) the standard upper triangular k-by-k Jordan block with eigenvalue ∗ Corresponding author.
    [Show full text]
  • Structure-Preserving Flows of Symplectic Matrix Pairs
    STRUCTURE-PRESERVING FLOWS OF SYMPLECTIC MATRIX PAIRS YUEH-CHENG KUO∗, WEN-WEI LINy , AND SHIH-FENG SHIEHz Abstract. We construct a nonlinear differential equation of matrix pairs (M(t); L(t)) that are invariant (Structure-Preserving Property) in the class of symplectic matrix pairs X12 0 IX11 (M; L) = S2; S1 X = [Xij ]1≤i;j≤2 is Hermitian ; X22 I 0 X21 where S1 and S2 are two fixed symplectic matrices. Furthermore, its solution also preserves de- flating subspaces on the whole orbit (Eigenvector-Preserving Property). Such a flow is called a structure-preserving flow and is governed by a Riccati differential equation (RDE) of the form > > W_ (t) = [−W (t);I]H [I;W (t) ] , W (0) = W0, for some suitable Hamiltonian matrix H . We then utilize the Grassmann manifolds to extend the domain of the structure-preserving flow to the whole R except some isolated points. On the other hand, the Structure-Preserving Doubling Algorithm (SDA) is an efficient numerical method for solving algebraic Riccati equations and nonlinear matrix equations. In conjunction with the structure-preserving flow, we consider two special classes of sym- plectic pairs: S1 = S2 = I2n and S1 = J , S2 = −I2n as well as the associated algorithms SDA-1 k−1 and SDA-2. It is shown that at t = 2 ; k 2 Z this flow passes through the iterates generated by SDA-1 and SDA-2, respectively. Therefore, the SDA and its corresponding structure-preserving flow have identical asymptotic behaviors. Taking advantage of the special structure and properties of the Hamiltonian matrix, we apply a symplectically similarity transformation to reduce H to a Hamiltonian Jordan canonical form J.
    [Show full text]
  • Hamiltonian Matrix Representations for the Determination of Approximate Wave Functions for Molecular Resonances
    Journal of Applied and Fundamental Sciences HAMILTONIAN MATRIX REPRESENTATIONS FOR THE DETERMINATION OF APPROXIMATE WAVE FUNCTIONS FOR MOLECULAR RESONANCES Robert J. Buenker Fachbereich C-Mathematik und Naturwissenschaften, Bergische Universität Wuppertal, Gaussstr. 20, D-42097 Wuppertal, Germany *For correspondence. ([email protected]) Abstract: Wave functions obtained using a standard complex Hamiltonian matrix diagonalization procedure are square integrable and therefore constitute only approximations to the corresponding resonance solutions of the Schrödinger equation. The nature of this approximation is investigated by means of explicit calculations using the above method which employ accurate diabatic potentials of the B 1Σ+ - D’ 1Σ+ vibronic resonance states of the CO molecule. It is shown that expanding the basis of complex harmonic oscillator functions gradually improves the description of the exact resonance wave functions out to ever larger internuclear distances before they take on their unwanted bound-state characteristics. The justification of the above matrix method has been based on a theorem that states that the eigenvalues of a complex-scaled Hamiltonian H (ReiΘ) are associated with the energy position and linewidth of resonance states (R is an internuclear coordinate and Θ is a real number). It is well known, however, that the results of the approximate method can be obtained directly using the unscaled Hamiltonian H (R) in real coordinates provided a particular rule is followed for the evaluation of the corresponding matrix elements. It is shown that the latter rule can itself be justified by carrying out the complex diagonalization of the Hamiltonian in real space via a product of two transformation matrices, one of which is unitary and the other is complex orthogonal, in which case only the symmetric scalar product is actually used in the evaluation of all matrix elements.
    [Show full text]