Source author record

Vanni Noferini

Vanni Noferini appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

Perturbation theory of transfer function matrices

Zeros of rational transfer function matrices $R(λ)$ are the eigenvalues of associated polynomial system matrices $P(λ)$, under minimality conditions. In this paper we define a structured condition number for a simple eigenvalue $λ_0$ of a (locally) minimal polynomial system matrix $P(λ)$, which in turn is a simple zero $λ_0$ of its transfer function matrix $R(λ)$. Since any rational matrix can be written as the transfer function of a polynomial system matrix, our analysis yield a structured perturbation theory for simple zeros of rational matrices $R(λ)$. To capture all the zeros of $R(λ)$, regardless of whether they are poles or not, we consider the notion of root vectors. As corollaries of the main results, we pay particular attention to the special case of $λ_0$ being not a pole of $R(λ)$ since in this case the results get simpler and can be useful in practice. We also compare our structured condition number with Tisseur's unstructured condition number for eigenvalues of matrix polynomials, and show that the latter can be unboundedly larger. Finally, we corroborate our analysis by numerical experiments.

preprint2022arXiv

The limit empirical spectral distribution of Gaussian monic complex matrix polynomials

We define the empirical spectral distribution (ESD) of a random matrix polynomial with invertible leading coefficient, and we study it for complex $n \times n$ Gaussian monic matrix polynomials of degree $k$. We obtain exact formulae for the almost sure limit of the ESD in two distinct scenarios: (1) $n \rightarrow \infty$ with $k$ constant and (2) $k \rightarrow \infty$ with $n$ constant. The main tool for our approach is the replacement principle by Tao, Vu and Krishnapur. Along the way, we also develop some auxiliary results of potential independent interest: we slightly extend a result by Bürgisser and Cucker on the tail bound for the norm of the pseudoinverse of a non-zero mean matrix, and we obtain several estimates on the singular values of certain structured random matrices.

preprint2021arXiv

Nearest $Ω$-stable matrix via Riemannian optimization

We study the problem of finding the nearest $Ω$-stable matrix to a certain matrix $A$, i.e., the nearest matrix with all its eigenvalues in a prescribed closed set $Ω$. Distances are measured in the Frobenius norm. An important special case is finding the nearest Hurwitz or Schur stable matrix, which has applications in systems theory. We describe a reformulation of the task as an optimization problem on the Riemannian manifold of orthogonal (or unitary) matrices. The problem can then be solved using standard methods from the theory of Riemannian optimization. The resulting algorithm is remarkably fast on small-scale and medium-scale matrices, and returns directly a Schur factorization of the minimizer, sidestepping the numerical difficulties associated with eigenvalues with high multiplicity.

preprint2021arXiv

The limit empirical spectral distribution of complex matrix polynomials

We study the empirical spectral distribution (ESD) for complex n x n matrix polynomials of degree k. We obtain exact formulae for the almost sure limit of the ESD in two distinct scenarios: (1) n -> \infty with k constant and (2) k -> \infty with n bounded by O(k^P) for some P>0. The main tools used are the logarithmic potential of some measure related to the matrix polynomial, and some classical estimates on the singular values of full random matrices with i.i.d. entries.

preprint2020arXiv

Wilkinson's bus: Weak condition numbers, with an application to singular polynomial eigenproblems

We propose a new approach to the theory of conditioning for numerical analysis problems for which both classical and stochastic perturbation theory fail to predict the observed accuracy of computed solutions. To motivate our ideas, we present examples of problems that are discontinuous at a given input and have infinite classical and stochastic condition number, but where the solution is still computed to machine precision without relying on structured algorithms. Stimulated by the failure of classical and stochastic perturbation theory in capturing such phenomena, we define and analyse a weak worst-case and a weak stochastic condition number. This new theory is a more powerful predictor of the accuracy of computations than existing tools, especially when the worst-case and the expected sensitivity of a problem to perturbations of the input is not finite. We apply our analysis to the computation of simple eigenvalues of matrix polynomials, including the more difficult case of singular matrix polynomials. In addition, we show how the weak condition numbers can be estimated in practice.

preprint2016arXiv

Numerical instability of resultant methods for multidimensional rootfinding

Hidden-variable resultant methods are a class of algorithms for solving multidimensional polynomial rootfinding problems. In two dimensions, when significant care is taken, they are competitive practical rootfinders. However, in higher dimensions they are known to miss zeros, calculate roots to low precision, and introduce spurious solutions. We show that the hidden variable resultant method based on the Cayley (Dixon or Bézout) matrix is inherently and spectacularly numerically unstable by a factor that grows exponentially with the dimension. We also show that the Sylvester matrix for solving bivariate polynomial systems can square the condition number of the problem. In other words, two popular hidden variable resultant methods are numerically unstable, and this mathematically explains the difficulties that are frequently reported by practitioners. Regardless of how the constructed polynomial eigenvalue problem is solved, severe numerical difficulties will be present. Along the way, we prove that the Cayley resultant is a generalization of Cramer's rule for solving linear systems and generalize Clenshaw's algorithm to an evaluation scheme for polynomials expressed in a degree-graded polynomial basis.

preprint2016arXiv

Vector spaces of linearizations for matrix polynomials: a bivariate polynomial approach

We revisit the landmark paper [D. S. Mackey, N. Mackey, C. Mehl, and V. Mehrmann, SIAM J. Matrix Anal. Appl., 28 (2006), pp.~971--1004] and, by viewing matrices as coefficients for bivariate polynomials, we provide concise proofs for key properties of linearizations for matrix polynomials. We also show that every pencil in the double ansatz space is intrinsically connected to a Bézout matrix, which we use to prove the eigenvalue exclusion theorem. In addition our exposition allows for any polynomial basis and for any field. The new viewpoint also leads to new results. We generalize the double ansatz space by exploiting its algebraic interpretation as a space of Bézout pencils to derive new linearizations with potential applications in the theory of structured matrix polynomials. Moreover, we analyze the conditioning of double ansatz space linearizations in the important practical case of a Chebyshev basis.

preprint2012arXiv

Locating the eigenvalues of matrix polynomials

Some known results for locating the roots of polynomials are extended to the case of matrix polynomials. In particular, a theorem by A.E. Pellet [Bulletin des Sciences Mathématiques, (2), vol 5 (1881), pp.393-395], some results of D.A. Bini [Numer. Algorithms 13:179-200, 1996] based on the Newton polygon technique, and recent results of M. Akian, S. Gaubert and M. Sharify (see in particular [LNCIS, 389, Springer p.p.291-303] and [M. Sharify, Ph.D. thesis, École Polytechnique, ParisTech, 2011]). These extensions are applied for determining effective initial approximations for the numerical computation of the eigenvalues of matrix polynomials by means of simultaneous iterations, like the Ehrlich-Aberth method. Numerical experiments that show the computational advantage of these results are presented.

preprint2012arXiv

The behaviour of the complete eigenstructure of a polynomial matrix under a generic rational transformation

Given a polynomial matrix P(x) of grade g and a rational function $x(y) = n(y)/d(y)$, where $n(y)$ and $d(y)$ are coprime nonzero scalar polynomials, the polynomial matrix $Q(y) :=[d(y)]^gP(x(y))$ is defined. The complete eigenstructures of $P(x)$ and $Q(y)$ are related, including characteristic values, elementary divisors and minimal indices. A Theorem on the matter, valid in the most general hypotheses, is stated and proved.

preprint2011arXiv

The Ehrlich-Aberth method for palindromic matrix polynomials represented in the Dickson basis

An algorithm based on the Ehrlich-Aberth root-finding method is presented for the computation of the eigenvalues of a T-palindromic matrix polynomial. A structured linearization of the polynomial represented in the Dickson basis is introduced in order to exploit the symmetry of the roots by halving the total number of the required approximations. The rank structure properties of the linearization allow the design of a fast and numerically robust implementation of the root-finding iteration. Numerical experiments that confirm the effectiveness and the robustness of the approach are provided.