An interpolation approach is pursued to solve a bilateral matrix polynomial equation frequently arising in control and signal processing. It results in eecient and numerically reliable resolution methods. Abstract New numerical procedures are proposed to solve the symmetric matrix polynomial equation A T (?s)X(s) + X T (?s)A(s) = 2B(s) that is frequently encountered in control and signal proces...

We present our public-domain software for the following tasks in sparse (or toric) elimination theory, given a well-constrained polynomial system. First, C code for computing the mixed volume of the system. Second, Maple code for defining an overconstrained system and constructing a Sylvester-type matrix of its sparse resultant. Third, C code for a Sylvester-type matrix of the sparse resultant ...

This paper describes the application of the Sylvester resultant matrix to image deblurring. In particular, an image is represented as a bivariate polynomial and it is shown that operations on polynomials, specifically greatest common divisor (GCD) computations and polynomial divisions, enable the point spread function to be calculated and an image to be deblurred. The GCD computations are perfo...

The task of determining the approximate greatest common divisor (GCD) of univariate polynomials with inexact coefficients can be formulated as computing for a given Sylvester matrix a new Sylvester matrix of lower rank whose entries are near the corresponding entries of that input matrix. We solve the approximate GCD problem by a new method based on structured total least norm (STLN) algorithms...

We present an eecient design and implementation for the solution of the Sylvester Observer matrix equation on a distributed memory computer. It consists of two main computational intensive segments, and one of a lesser complexity: the reduction of the system matrix A to a Hessenberg form, the solution of n shifted Hessenberg systems and a number of matrix-matrix multiplies. The Hessenberg reduc...

In this talk we suggest a new formula for the residual of Galerkin projection onto rational Krylov spaces applied to a Sylvester equation, and establish a relation to three different underlying extremal problems for rational functions. These extremal problems enable us to compare the size of the residual for the above method with that obtained by ADI. In addition, we may deduce several new a pr...

Sylvester equations play a central role in many areas of applied mathematics and in particular in systems and control theory. Here we will show how low-rank solutions for stable Sylvester equations can be computed based on the matrix sign function method. We discuss applications in model reduction as well as in observer design. Numerical experiments demonstrate the efficiency of the proposed me...