Subspace method for the estimation of large-scale structured real stability radius

We consider the autonomous dynamical system x '= This linear dynamical system is asymptotically stable if all of the eigenvalues of A lie in the open left-half of the complex plane. In this case, the matrix A is said to be Hurwitz stable or shortly a stable matrix. In practice, the stability of a system can be violated because of perturbations such as modeling errors. In such cases, one deals with the robust stability of the system rather than its stability. The system above is said to be robustly stable if the system, as well as all of its perturbations from a certain perturbation class, are stable. To measure the robustness of the system subject to perturbations, a quantity of interest is the stability radius or in other words the distance to instability. In this paper, we focus on the estimation of the structured real stability radius for large-scale systems. We propose a subspace framework to estimate the structured real stability radius and prove that our new method converges at a quadratic rate in theory. Our method benefits from a one-sided interpolatory model order reduction technique, in the sense that the left and the right subspaces are the same. The quadratic convergence of the method is due to the certain Hermite interpolation properties between the full and reduced problems. The proposed framework estimates the structured real stability radius for large-scale systems efficiently. The efficiency of the method is demonstrated on several numerical experiments.