165 research outputs found
Structured total least norm and approximate GCDs of inexact polynomials
The determination of an approximate greatest common divisor (GCD) of two inexact polynomials f=f(y) and g=g(y) arises in several applications, including signal processing and control. This approximate GCD can be obtained by computing a structured low rank approximation S*(f,g) of the Sylvester resultant matrix S(f,g). In this paper, the method of structured total least norm (STLN) is used to compute a low rank approximation of S(f,g), and it is shown that important issues that have a considerable effect on the approximate GCD have not been considered. For example, the established works only yield one matrix S*(f,g), and therefore one approximate GCD, but it is shown in this paper that a family of structured low rank approximations can be computed, each member of which yields a different approximate GCD. Examples that illustrate the importance of these and other issues are presented
A non-linear structure preserving matrix method for the low rank approximation of the Sylvester resultant matrix
A non-linear structure preserving matrix method for the computation of a structured low rank approximation S((f) over bar , (g) over bar) of the Sylvester resultant matrix S(f , g) of two inexact polynomials f = f(y) and g = g(y) is considered in this paper. It is shown that considerably improved results are obtained when f (y) and g(y) are processed prior to the computation of S((f) over bar , (g) over bar), and that these preprocessing operations introduce two parameters. These parameters can either be held constant during the computation of S((f) over bar , (g) over bar), which leads to a linear structure preserving matrix method, or they can be incremented during the computation of S((f) over bar, (g) over bar), which leads to a non-linear structure preserving matrix method. It is shown that the non-linear method yields a better structured low rank approximation of S((f) over bar , (g) over bar) and that the assignment of f (y) and g(y) is important because S((f) over bar , (g) over bar) may be a good structured low rank approximation of S(f, g), but S((f) over bar , (g) over bar) may be a poor structured low rank approximation of S (g f) because its numerical rank is not defined. Examples that illustrate the differences between the linear and non-linear structure preserving matrix methods, and the importance of the assignment off (y) and g(y), are shown. (C) 2010 Elsevier B.V. All rights reserved
Over-constrained Weierstrass iteration and the nearest consistent system
We propose a generalization of the Weierstrass iteration for over-constrained
systems of equations and we prove that the proposed method is the Gauss-Newton
iteration to find the nearest system which has at least common roots and
which is obtained via a perturbation of prescribed structure. In the univariate
case we show the connection of our method to the optimization problem
formulated by Karmarkar and Lakshman for the nearest GCD. In the multivariate
case we generalize the expressions of Karmarkar and Lakshman, and give
explicitly several iteration functions to compute the optimum.
The arithmetic complexity of the iterations is detailed
Generalized companion matrix for approximate GCD
We study a variant of the univariate approximate GCD problem, where the
coefficients of one polynomial f(x)are known exactly, whereas the coefficients
of the second polynomial g(x)may be perturbed. Our approach relies on the
properties of the matrix which describes the operator of multiplication by gin
the quotient ring C[x]=(f). In particular, the structure of the null space of
the multiplication matrix contains all the essential information about GCD(f;
g). Moreover, the multiplication matrix exhibits a displacement structure that
allows us to design a fast algorithm for approximate GCD computation with
quadratic complexity w.r.t. polynomial degrees.Comment: Submitted to MEGA 201
Multivariate Subresultants in Roots
We give rational expressions for the subresultants of n+1 generic polynomials
f_1,..., f_{n+1} in n variables as a function of the coordinates of the common
roots of f_1,..., f_n and their evaluation in f_{n+1}. We present a simple
technique to prove our results, giving new proofs and generalizing the
classical Poisson product formula for the projective resultant, as well as the
expressions of Hong for univariate subresultants in roots.Comment: 22 pages, no figures, elsart style, revised version of the paper
presented in MEGA 2005, accepted for publication in Journal of Algebr
Comparison of algorithms for calculation of the greatest common divisor of several polynomials
summary:The computation of the greatest common divisor (GCD) has many applications in several disciplines including computer graphics, image deblurring problem or computing multiple roots of inexact polynomials. In this paper, Sylvester and Bézout matrices are considered for this purpose. The computation is divided into three stages. A rank revealing method is shortly mentioned in the first one and then the algorithms for calculation of an approximation of GCD are formulated. In the final stage the coefficients are improved using Gauss-Newton method. Numerical results show the efficiency of proposed last two stages
- …