We analyze robust stability and performance of dynamical systems with real uncertain parameters. We compute criteria like the distance to instability, the worst-case spectral abscissa, or the worst-case H∞-norm, which quantify the degree of robustness of such a system when parameters vary in a given set ∆. As computing these indices is NP-hard, we present a heuristic which finds good lower bounds fast and reliably. Posterior certification is then obtained by an intelligent global strategy. A test bench of 87 systems with up to 70 states 39 uncertain parameters with up to 11 repetitions demonstrates the potential of our approach.1 Problem Specification.Robustness specifications limit the loss of performance and stability in a system where differences between the mathematical model and reality crop up. Robustness against real parametric uncertainties is among the most challenging calls in this regard. Already deciding whether a given system with uncertain real parameters δ is robustly stable over a given parameter range δ ∈ ∆ is NP-hard, and this is aggravated when it comes to deciding whether a given level of H 2 -or H ∞ -performance is guaranteed over that range. In this work we address this type of uncertainty by computing three key criteria, which quantify the degree of parametric robustness of a system. These are (a) the worst-case H ∞ -norm, and (b) the worst-case spectral abscissa over a given parameter range, and (c) the distance to instability, or stability margin, of a system with uncertain parameters.Consider a Linear Fractional Transform [23] with real parametric uncertainties as in Figure 1, whereand x ∈ R n is the state, w ∈ R m1 a vector of exogenous inputs, and z ∈ R p1 a vector of regulated outputs. The uncertainty channel is defined as p = ∆q, where the time-invariant uncertain matrix ∆ has the blockdiagonal form2) * ONERA, Control System Department, Toulouse, France † Université de Toulouse, Institut de Mathématiques with δ 1 , . . . , δ m representing real uncertain parameters, and r i giving the number of repetitions of δ i . Here we assume without loss that δ = 0 ∈ ∆ represents the nominal parameter value, and we consider δ ∈ ∆ in one-to-one correspondence with the matrix ∆ in (1.2). For practical applications it is generally sufficient to consider the caseTo analyze the performance of (1.1) in the presence of the uncertain δ ∈ R m we compute the worst-casewhere · ∞ is the H ∞ -norm, and where T wz (s, δ) is the transfer function z(s) = F u (P (s), ∆)w(s), obtained by closing the loop between (1.1) and p = ∆q with (1.2) in Figure 1. The solution δ * ∈ ∆ of (1.3) represents a worst possible choice of the parameters δ ∈ ∆, which may be an important element in analyzing performance and robustness of the system, see e.g. [1].Our second criterion is similar in nature, as it allows to verify whether the uncertain system (1.1) is robustly stable over a given parameter range ∆. This can be tested by maximizing the spectral abscissa of the system A-matrix over the parameter range−1 C q , and w...