The Convergence of the Geometric Interpolation Algorithm.
Hongwei Lin
DOI: https://doi.org/10.1016/j.cad.2010.01.006
IF: 3.652
2010-01-01
Computer-Aided Design
Abstract:The geometric interpolation algorithm is proposed by Maekawa et al. in [Maekawa T, Matsumoto Y, Namiki K. Interpolation by geometric algorithm. Computer-Aided Design 2007:39:313-23]. Without solving a system of equations, the algorithm generates a curve (surface) sequence, of which the limit curve (surface) interpolates the given data points. However, the convergence of the algorithm is a conjecture in the reference above, and demonstrated by lots of empirical examples. In this paper, we prove the conjecture given in the reference in theory, that is, the geometric interpolation algorithm is convergent for a blending curve (surface) with normalized totally positive basis, under the condition that the minimal eigenvalue lambda(min)(D(k)) of the collocation matrix D(k) of the totally positive basis in each iteration satisfies lambda(min)(D(k)) >= alpha > 0. As a consequence, the geometric interpolation algorithm is convergent for Bezier, B-spline, rational Bezier, and NURBS curve (surface) if they satisfy the condition aforementioned, since Bernstein basis and B-spline basis are both normalized totally positive. (C) 2010 Elsevier Ltd. All rights reserved.