Abstract:Abstract. The ε-pseudospectral abscissa αε and radius ρε of an n × n matrix are, respectively, the maximal real part and the maximal modulus of points in its ε-pseudospectrum, defined using the spectral norm. It was proved in [A.S. Lewis and C.H.J. Pang, SIAM J. Optim., 19 (2008), pp. 1048-1072] that for fixed ε > 0, αε and ρε are Lipschitz continuous at a matrix A except when αε and ρε are attained at a critical point of the norm of the resolvent (in the nonsmooth sense), and it was conjectured that the point… Show more
“…Remark 4. Numerical algorithms for computing the stability radii for ODEs are proposed in a number of works, e.g., see [4,6,12,34,38,39,40,42,45,71,72,73,81]. Some extensions to DAEs are discussed in [5].…”
This paper presents a survey of recent results on the robust stability analysis and the distance to instability for linear time-invariant and time-varying differential-algebraic equations (DAEs). Different stability concepts such as exponential and asymptotic stability are studied and their robustness is analyzed under general as well as restricted sets of real or complex perturbations. Formulas for the distances are presented whenever these are available and the continuity of the distances in terms of the data is discussed. Some open problems and challenges are indicated.
“…Remark 4. Numerical algorithms for computing the stability radii for ODEs are proposed in a number of works, e.g., see [4,6,12,34,38,39,40,42,45,71,72,73,81]. Some extensions to DAEs are discussed in [5].…”
This paper presents a survey of recent results on the robust stability analysis and the distance to instability for linear time-invariant and time-varying differential-algebraic equations (DAEs). Different stability concepts such as exponential and asymptotic stability are studied and their robustness is analyzed under general as well as restricted sets of real or complex perturbations. Formulas for the distances are presented whenever these are available and the continuity of the distances in terms of the data is discussed. Some open problems and challenges are indicated.
“…Like the spectral radius, the pseudospectral radius is nonconvex, nonsmooth, and continuously differentiable almost everywhere. However, in contrast to the spectral radius, the pseudospectral radius is locally Lipschitz [GO12] and is thus potentially an easier function to optimize. For example, the known convergence rates for gradient sampling hold for minimizing the pseudospectral radius but not the spectral radius.…”
We propose an algorithm for solving nonsmooth, nonconvex, constrained optimization problems as well as a new set of visualization tools for comparing the performance of optimization algorithms. Our algorithm is a sequential quadratic optimization method that employs Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton Hessian approximations and an exact penalty function whose parameter is controlled using a steering strategy. While our method has no convergence guarantees, we have found it to perform very well in practice on challenging test problems in controller design involving both locally Lipschitz and non-locally-Lipschitz objective and constraint functions with constraints that are typically active at local minimizers. In order to empirically validate and compare our method with available alternatives-on a new test set of 200 problems of varying sizes-we employ new visualization tools which we call relative minimization profiles. Such profiles are designed to simultaneously assess the relative performance of several algorithms with respect to objective quality, feasibility, and speed of progress, highlighting the trade-offs between these measures when comparing algorithm performance.
The stability radius of an n × n matrix A (or distance to instability) is a well-known measure of robustness of stability of the linear stable dynamical systemẋ = Ax. Such a distance is commonly measured either in the 2-norm or in the Frobenius norm. Even if the matrix A is real, the distance to instability is most often considered with respect to complex valued matrices (in such case the two norms turn out to be equivalent) and restricting the distance to real matrices makes the problem more complicated, and in the case of Frobenius norm-to our knowledge-unresolved. Here we present a novel approach to approximate real stability radii, particularly well-suited for large sparse matrices. The method consists of a two level iteration, the inner one aiming to compute the εpseudospectral abscissa of a low-rank (1 or 2) dynamical system, and the outer one consisting of an exact Newton iteration. Due to its local convergence property it generally provides upper bounds for the stability radii but in practice usually computes the correct values. The method requires the computation of the rightmost eigenvalue of a sequence of matrices, each of them given by the sum of the original matrix A and a low-rank one. This makes it particularly suitable for large sparse problems, for which several existing methods become inefficient, due to the fact that they require to solve full Hamiltonian eigenvalue problems and/or compute multiple SVDs.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.