Computing Extremal Points of Symplectic Pseudospectra and Solving Symplectic Matrix Nearness Problems

We study differential equations that lead to extremal points in symplectic pseudospectra. In a two-level approach, where on the inner level we compute extremizers of the symplectic $\varepsilon$-pseudospectrum for a given $\varepsilon$ and on the outer level we optimize over $\varepsilon$, this is used to solve symplectic matrix nearness problems such as the following: For a symplectic matrix with eigenvalues of unit modulus, we aim to determine the nearest complex symplectic matrix such that some or all eigenvalues leave the complex unit circle. Conversely, for a symplectic matrix with all eigenvalues lying off the unit circle, we consider the problem of computing the nearest symplectic matrix that has an eigenvalue on the unit circle.

[1]  Francesca Bonizzoni,et al.  Perturbation Analysis for the Darcy Problem with Log-Normal Permeability , 2014, SIAM/ASA J. Uncertain. Quantification.

[2]  Daniel Kressner,et al.  On the computation of structured singular values and pseudospectra , 2010, Syst. Control. Lett..

[3]  Brian D. O. Anderson,et al.  Riccati Equations, Network Theory and Brune Synthesis: Old Solutions for Contemporary Problems , 1999 .

[4]  A. Bunse-Gerstner.  Matrix factorizations for symplectic QR-like methods , 1986 .

[5]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[6]  Marcus J. Grote,et al.  Finite Element Heterogeneous Multiscale Method for the Wave Equation: Long-Time Effects , 2014, Multiscale Model. Simul..

[7]  E. Celledoni.  Lie group methods , 2009 .

[8]  D. Hinrichsen,et al.  Mathematical Systems Theory I: Modelling, State Space Analysis, Stability and Robustness , 2010 .

[9]  L. Trefethen,et al.  Spectra and Pseudospectra , 2020 .

[10]  N. Guglielmi,et al.  Differential equations for real-structured (and unstructured) defectivity measures , 2014 .

[11]  Assyr Abdulle,et al.  High Order Numerical Approximation of the Invariant Measure of Ergodic SDEs , 2014, SIAM J. Numer. Anal..

[12]  Diederich Hinrichsen,et al.  Modelling, state space analysis, stability and robustness , 2005 .

[13]  Meiyue Shao.  On the finite section method for computing exponentials of doubly-infinite skew-Hermitian matrices , 2014 .

[14]  Anders Rantzer,et al.  Real Perturbation Values and Real Quadratic Forms in a Complex Vector Space , 1998 .

[15]  A. Laub.  A schur method for solving algebraic Riccati equations , 1978, 1978 IEEE Conference on Decision and Control including the 17th Symposium on Adaptive Processes.

[16]  G. Stewart,et al.  Matrix Perturbation Theory , 1990 .

[17]  Nicola Guglielmi,et al.  Differential Equations for Roaming Pseudospectra: Paths to Extremal Points and Boundary Tracking , 2011, SIAM J. Numer. Anal..

[18]  Albert Cohen,et al.  Discrete least squares polynomial approximation with random evaluations − application to parametric and stochastic elliptic PDEs , 2015 .

[19]  Daniel Kressner,et al.  Low rank differential equations for Hamiltonian matrix nearness problems , 2015, Numerische Mathematik.

[20]  A. Quarteroni,et al.  Thermodynamically consistent orthotropic activation model capturing ventricular systolic wall thickening in cardiac electromechanics , 2014 .

[21]  Robert E. Mahony,et al.  Optimization Algorithms on Matrix Manifolds , 2007 .

[22]  Othmar Koch,et al.  Dynamical Low-Rank Approximation , 2007, SIAM J. Matrix Anal. Appl..

[23]  S. Rump.  EIGENVALUES, PSEUDOSPECTRUM AND STRUCTURED PERTURBATIONS , 2006 .

[24]  Nicola Guglielmi,et al.  Low-Rank Dynamics for Computing Extremal Points of Real Pseudospectra , 2013, SIAM J. Matrix Anal. Appl..

[25]  A. Iserles,et al.  Lie-group methods , 2000, Acta Numerica.

[26]  H. Faßbender.  Symplectic Methods for the Symplectic Eigenproblem , 2002, Springer US.

[27]  Volker Mehrmann,et al.  Perturbation Theory for Hamiltonian Matrices and the Distance to Bounded-Realness , 2011, SIAM J. Matrix Anal. Appl..

[28]  B. Datta.  Numerical methods for linear control systems : design and analysis , 2004 .

[29]  Stef Graillat,et al.  STRUCTURED CONDITION NUMBERS AND BACKWARD ERRORS IN SCALAR PRODUCT SPACES , 2006 .

[30]  Federico Poloni,et al.  Enforcing solvability of a nonlinear matrix equation and estimation of multivariate ARMA time series , 2013 .

[31]  Michael Karow,et al.  Structured Pseudospectra for Small Perturbations , 2011, SIAM J. Matrix Anal. Appl..

[32]  Daniel Kressner,et al.  Algorithm 953 , 2015 .

[33]  Gianluigi Rozza,et al.  Multilevel and weighted reduced basis method for stochastic optimal control problems constrained by Stokes equations , 2016, Numerische Mathematik.

[34]  Peter Benner,et al.  DAREX - A Collection of Benchmark Examples for Discrete-Time Algebraic Riccati Equations (Version , 1999 .

[35]  André Uschmajew,et al.  On Convergence of the Maximum Block Improvement Method , 2015, SIAM J. Optim..

[36]  G. Rozza,et al.  A combination between the reduced basis method and the ANOVA expansion: On the computation of sensitivity indices , 2013 .