City Research Online

An alternating projection algorithm for the “approximate” GCD calculation

Limantseva, O. ORCID: 0000-0003-2764-4415, Halikias, G. ORCID: 0000-0003-1260-1383 & Karcanias, N. (2020). An alternating projection algorithm for the “approximate” GCD calculation. In: IFAC-PapersOnLine. 21st IFAC World Congress, 11-17 Jul 2020, Berlin, Germany. doi: 10.1016/j.ifacol.2020.12.1630

Abstract

In the paper an approach is proposed for calculating the “best” approximate GCD of a set of coprime polynomials. The algorithm is motivated by the factorisation of the Sylvester resultant matrix of polynomial sets with nontrivial GCD. In the (generic) case of coprime polynomial sets considered here the aim is to minimise the norm of the residual error matrix of the inexact factorisation in order to compute the “optimal” approximate GCD. A least-squares alternating projection algorithm is proposed as an alternative to the solution of the corresponding optimisation problem via nonlinear programming techniques. The special structure of the problem in this case, however, means that the algorithm can be reduced to a sequence of standard subspace projections and hence no need arises to compute gradient vectors, Hessian matrices or optimal step-lengths. An estimate of the asymptotic convergence rate of the algorithm is finally established via the inclination of two subspaces.

Publication Type: Conference or Workshop Item (Paper)
Publisher Keywords: nonlinear least squares, Sylvester resultant matrix, approximate GCD, alternating projection algorithm, convergence analysis
Subjects: Q Science > QA Mathematics
Departments: School of Science & Technology > Engineering
[thumbnail of 1-s2.0-S240589632032228X-main.pdf]
Preview
Text - Published Version
Available under License Creative Commons Attribution Non-commercial No Derivatives.

Download (412kB) | Preview

Export

Add to AnyAdd to TwitterAdd to FacebookAdd to LinkedinAdd to PinterestAdd to Email

Downloads

Downloads per month over past year

View more statistics

Actions (login required)

Admin Login Admin Login