On computing with the Hilbert spline transform

  • Authors:
  • Charles A. Micchelli;Yuesheng Xu;Bo Yu

  • Affiliations:
  • Department of Mathematics, City University of Hong Kong, Hong Kong, China and Department of Mathematics and Statistics, State University of New York, The University at Albany, Albany, USA 12222;Department of Mathematics, Syracuse University, Syracuse, USA 13244 and Guangdong Province Key Lab of Computational Science, Sun Yat-sen University, Guangzhou, People's Republic of China 510275;Department of Mathematics, College of Science, China Three Gorges University, Yichang, China 443002

  • Venue:
  • Advances in Computational Mathematics
  • Year:
  • 2013

Quantified Score

Hi-index 0.00

Visualization

Abstract

We develop a fast algorithm for computing the Hilbert transform of a function from a data set consisting of n function values and prove that the complexity of the proposed algorithm is O(n log n). Our point of view is fundamentally based on a B-spline series approximation constructed from available data. In this regard, we obtain new formulas for the Hilbert transform of a B-spline as a divided difference. For theoretical simplicity and computational efficiency we give a detailed description of our algorithm, as well as provided optimal approximation order, only in the case of quadratic splines. However, if higher accuracy is required, extensions of our method to spline approximation of any prescribed degree readily follows the pattern of the quadratic case. Numerical experiments have confirmed that our algorithm has superior performance than previously available methods which we briefly survey in Section 2.