Unstable dynamical systems can be viewed from a variety of perspectives. We discuss the potential of an input-output map associated with an unstable system to represent a bounded map from ℒ 2 (ℝ) to itself and then develop criteria for optimal reduced order approximations to the original (unstable) system with respect to an ℒ 2 -induced Hilbert-Schmidt norm. Our optimality criteria extend the Meier-Luenberger interpolation conditions for optimal ℋ 2  approximation of stable dynamical systems. Based on this interpolation framework, we describe an iteratively corrected rational Krylov algorithm for ℒ 2  model reduction. A numerical example involving a hard-to-approximate full-order model illustrates the effectiveness of the proposed approach.
['Caleb Magruder', 'Christopher A. Beattie', 'Serkan Gugercin']
Rational Krylov methods for optimal ℒ 2 model reduction