On numerical methods for discrete least-squares approximation by trigonometric polynomials
Abstract: Fast, efficient and reliable algorithms for discrete least-squares approximation of a real-valued function given at arbitrary distinct nodes in by trigonometric polynomials are presented. The algorithms are based on schemes for the solution of inverse unitary eigenproblems and require only arithmetic operations as compared to operations needed for algorithms that ignore the structure of the problem. An algorithm which solves this problem with real-valued data and real-valued solution using only real arithmetic is given. Numerical examples are presented that show that the proposed algorithms produce consistently accurate results that are often better than those obtained by general QR decomposition methods for the least-squares problem.
Affiliation: Universität Bremen, Fachbereich 3 Mathematik und Informatik, 28334 Bremen, Germany
Keywords: Trigonometric approximation, unitary Hessenberg matrix, Schur parameter
Received by editor(s): March 29, 1995
Received by editor(s) in revised form: January 31, 1996
Article copyright: © Copyright 1997 American Mathematical Society