Polynomial least squares fitting in the Bernstein basis

From MaRDI portal
(Redirected from Publication:5962279)



Abstract: The problem of polynomial regression in which the usual monomial basis is replaced by the Bernstein basis is considered. The coefficient matrix A of the overdetermined system to be solved in the least squares sense is then a rectangular Bernstein-Vandermonde matrix. In order to use the method based on the QR decomposition of A, the first stage consists of computing the bidiagonal decomposition of the coefficient matrix A. Starting from that bidiagonal decomposition, an algorithm for obtaining the QR decomposition of A is the applied. Finally, a triangular system is solved by using the bidiagonal decomposition of the R-factor of A. Some numerical experiments showing the behavior of this approach are included.


Based on the fact that in the least squares method used in the polynomial fitting, the rectangular Vandermonde matrix corresponding to the nodes, associated with the overdetermined system of the normal equations is usually an ill-conditioned matrix and the use of QR factorization doesn't solve too much of this problem, the paper is focused on the idea to take, instead of the classical monomial basis, a polynomial basis which leads to a matrix of the nodes with smaller condition number than the Vandermonde matrix.NEWLINENEWLINEThus, the Bernstein basis of polynomials, widely used in computer aided geometric design is considered. Moreover, the nodes are taken strictly increasing between 0 and 1. To solve this type of polynomial fitting, the authors use some results related to total positivity of the rectangular Bernstein-Vandermonde matrix, defined for stricly increasing nodes, between 0 and 1.NEWLINENEWLINEThe bidiagonal factorization is realized as starting point for the computation of QR decomposition for the Bernstein-Vandermonde matrix of the nodes, then a triangular system is solved by using the bidiagonal decomposition of the R factor. NEWLINENEWLINEIn this order, the authors propose an efficient algorithm with high accuracy dedicated to solve the problem of polynomial regression in the Bernstein basis. Finally, two numerical experiments illustrate the efficiency of the proposed algorithm showing that the accurcay is keeped when the condition number of the associated Bernstein-Vandermonde matrix increases.




Cited in
(22)


Describes a project that uses

Uses Software






This page was built for publication: Polynomial least squares fitting in the Bernstein basis

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5962279)