English

On skew-Hamiltonian Matrices and their Krylov-Lagrangian Subspaces

Numerical Analysis 2020-12-01 v2 Numerical Analysis

Abstract

It is a well-known fact that the Krylov space Kj(H,x)\mathcal{K}_j(H,x) generated by a skew-Hamiltonian matrix HR2n×2nH \in \mathbb{R}^{2n \times 2n} and some xR2nx \in \mathbb{R}^{2n} is isotropic for any jNj \in \mathbb{N}. For any given isotropic subspace LR2n\mathcal{L} \subset \mathbb{R}^{2n} of dimension nn - which is called a Lagrangian subspace - the question whether L\mathcal{L} can be generated as the Krylov space of some skew-Hamiltonian matrix is considered. The affine variety HK\mathbb{HK} of all skew-Hamiltonian matrices HR2n×2nH \in \mathbb{R}^{2n \times 2n} that generate L\mathcal{L} as a Krylov space is analyzed. Existence and uniqueness results are proven, the dimension of HK\mathbb{HK} is found and skew-Hamiltonian matrices with minimal 22-norm and Frobenius norm in HK\mathbb{HK} are identified. In addition, a simple algorithm is presented to find a basis of HK\mathbb{HK}.

Keywords

Cite

@article{arxiv.1910.12904,
  title  = {On skew-Hamiltonian Matrices and their Krylov-Lagrangian Subspaces},
  author = {Philip Saltenberger and Michel-Niklas Senn},
  journal= {arXiv preprint arXiv:1910.12904},
  year   = {2020}
}