Generalized Pellet's theorem for analytic matrix-valued functions

Abstract

Generalized Pellet's theorem is extended to analytic matrix-valued functions using Descartes' rule of signs for infinite series. Cauchy's estimate is incorporated to allow numerical computation of the bounds. To address the issue that Pellet's theorem and its generalizations may not always guarantee the existence of bounds, for the case of matrix size one, we employ the Dandelin- Lobachevsky-Graeffe method, which ensures the existence of the bounds for the roots and further improves them as the order of the transform increases. Without the need for any a priori information on the function, these theorems allow us to compute mathematically verified bounds for the roots in a certain region using an interval arithmetic implementation of the function. Numerical examples and an implementation are provided to illustrate the effectiveness of the proposed method.