DIFFERENTIAL EQUATIONS FOR REAL-STRUCTURED DEFECTIVITY MEASURES

作者：Butta P^*; Guglielmi N; Manetta M; Noschese S

来源：SIAM Journal on Matrix Analysis and Applications, 2015, 36(2): 523-548.

DOI：10.1137/140964631

• pseudospectrum
• real-structured pseudospectrum
• low-rank differential equations
• defective eigenvalue
• distance to defectivity
• Wilkinson problem

Let A be a real matrix with all distinct eigenvalues. We propose a new method for the computation of the distance w(R)(A) of the matrix A from the set of real defective matrices, i.e., the set of those real matrices with at least one multiple eigenvalue with algebraic multiplicity larger than its geometric multiplicity. For 0 < epsilon <= w(R)(A), this problem is closely related to the computation of the most ill-conditioned epsilon-pseudoeigenvalues of A, that is, points in the epsilon-pseudospectrum of A characterized by the highest condition number. The method we propose couples a system of differential equations on a low-rank manifold which determines the epsilon-pseudoeigenvalue closest to coalescence, with a fast Newton-like iteration aiming to determine the minimal value epsilon such that an epsilon-pseudoeigenvalue becomes defective. The method has a local behavior; this means that, in general, we find upper bounds for w(R)(A). However, these bounds usually provide good approximations, in those (simple) cases where we can check this. The methodology can be extended to a structured matrix, where it is required that the distance be computed within some manifold defining the structure of the matrix. In this paper we extensively examine the case of real matrices. As far as we know, there do not exist methods in the literature able to compute such distance.

• 出版日期2015