Bivariate Free Knot Splines

We consider the least squares approximation of gridded 2D data by tensor product splines with free knots. The smoothing functional to be minimized—a generalization of the univariate Schoenberg functional—is chosen in such a way that the solution of the bivariate problem separates into the solution of a sequence of univariate problems in case of fixed knots. The resulting optimization problem is a constrained separable least squares problem with tensor product structure. Based on some ideas developed by the authors for the univariate case, an efficient method for solving the specially structured 2D problem is proposed, analyzed and tested on hand of some examples from the literature.

[1]  R. E. Carlson,et al.  An algorithm for monotone piecewise bicubic interpolation , 1989 .

[2]  R. E. Carlson,et al.  Monotone Piecewise Cubic Interpolation , 1980 .

[3]  P. Dierckx.  An algorithm for surface-fitting with spline functions , 1981 .

[4]  Douglas M. Bates,et al.  Loosely coupled nonlinear least squares , 1992 .

[5]  P. Dierckx.  A Fast Algorithm for Smoothing Data on a Rectangular Grid while Using Spline Functions , 1982 .

[6]  D. Fausett,et al.  Large Least Squares Problems Involving Kronecker Products , 1994, SIAM J. Matrix Anal. Appl..

[7]  Günther Nürnberger,et al.  Bivariate segment approximation and free knot splines: Research problems 96-4 , 1996 .

[8]  Thomas F. Coleman,et al.  Optimization Toolbox User's Guide , 1998 .

[9]  Larry L. Schumaker,et al.  Bivariate Natural Spline Smoothing , 1985 .

[10]  Larry L. Schumaker,et al.  Complete spline smoothing , 1986 .

[11]  Hubert Schwetlick,et al.  Constrained approximation by splines with free knots , 1997 .

[12]  Hubert Schwetlick,et al.  Least squares approximation by splines with free knots , 1995 .

[13]  Michael A. Saunders,et al.  User''s guide for NPSOL (Ver-sion 4.0): A FORTRAN package for nonlinear programming , 1984 .

[14]  J. Navarro-Pedreño.  Numerical Methods for Least Squares Problems , 1996 .

[15]  Linda Kaufman.  Tradeoffs in Algorithms for Separable Nonlinear Least Squares , 1990 .

[16]  Linda Kaufman,et al.  Structured Linear Least-Squares Problems in System Identification and Separable Nonlinear Data Fitting , 1994, SIAM J. Optim..

[17]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[18]  Hubert Schwetlick,et al.  Spline smoothing under constraints on derivatives , 1993 .

[19]  Linda Kaufman,et al.  Separable Nonlinear Least Squares with Multiple Right-Hand Sides , 1992, SIAM J. Matrix Anal. Appl..

[20]  Marion Walther,et al.  Gridded Data Interpolation with Restrictions on the First Order Derivatives , 1997 .

[21]  Günther Nürnberger,et al.  Bivariate segment approximation and splines , 1996, Adv. Comput. Math..