A Subspace Framework for ${\mathcal L}_\infty$ Model Reduction

TL;DR

This paper introduces a novel subspace framework for ${ m L}_ fty$ model reduction that employs smooth optimization on smaller interpolated systems, achieving fast convergence and near-optimal solutions for large-scale descriptor systems.

Contribution

It proposes a new subspace-based approach using interpolation and smooth optimization techniques to efficiently solve large-scale ${ m L}_ fty$ model reduction problems, improving convergence and solution quality.

Findings

The method converges quadratically under certain conditions.

It produces locally optimal solutions efficiently for systems with thousands of states.

Numerical experiments demonstrate rapid convergence and high-quality reductions.

Abstract

We consider the problem of locating a nearest descriptor system of prescribed reduced order to a descriptor system with large order with respect to the ${\mathcal L}_\infty$ norm. Widely employed approaches such as the balanced truncation and best Hankel norm approximation for this ${\mathcal L}_\infty$ model reduction problem are usually expensive and yield solutions that are not optimal, not even locally. We propose approaches based on the minimization of the ${\mathcal L}_\infty$ objective by means of smooth optimization techniques. As we illustrate, direct applications of smooth optimization techniques are not feasible, since the optimization techniques converge at best at a linear rate requiring too many evaluations of the costly ${\mathcal L}_\infty$-norm objective to be practical. We replace the original large-scale system with a system of smaller order that interpolates the original system at points on the imaginary axis, and minimize the ${\mathcal L}_\infty$ objective after this replacement. The smaller system is refined by interpolating at additional imaginary points determined based on the local minimizer of the ${\mathcal L}_\infty$ objective, and the optimization is repeated. We argue the framework converges at a quadratic rate under smoothness and nondegeneracy assumptions, and describe how asymptotic stability constraints on the reduced system sought can be incorporated into our approach. The numerical experiments on benchmark examples illustrate that the approach leads to locally optimal solutions to the ${\mathcal L}_\infty$ model reduction problem, and the convergence occurs quickly for descriptors systems of order a few ten thousands.