Introduction to Splines 2.0 2003/05/05 2003/05/05 Alena Ixmocane Scott oetting@stat.rice.edu Alena Ixmocane Scott oetting@stat.rice.edu Sujesh Sreedharan sujesh@rice.edu Rebecca Willett willett@rice.edu Adan Galvan jago@rice.edu Charlet Reedstrom charlet@rice.edu B-splines interpolation piecewise polynomials smoothing Splines wavelets Basic introduction to splines. Link to other modules. Suppose that you are given a data set D x 1 y 1 … x m y m in 2 and that you must find a smooth curve to these points. There are many solutions to such a problem. Most simply, you could "eyeball it" and draw in a curve that looks good; however, this solution is too subjective, not very precise. You could just connect the dots, but this will probably not give a smooth line. You could fit a polynomial, but this fitted polynomial may wiggle too much. Or you could use splines. Splines are piecewise polynomials with pieces smoothly connected together. The joining points of the polynomial pieces are called knots. Knots do not have to be evenly spaced. When each segment of a spline is a polynomial of degree n, we say that the spline is a spline of degree n. We need to add some constraints to ensure smoothness. For a spline of degree n, we require that the spline has continuous derivatives up to order n 1 at each of the knots, i.e. a spline of degree n is in C n 1 . Remember our problem of fitting a curve to the data? There are generally two ways of approaching this problem: interpolation and smoothing. We can fit a spline to interpolate the data; i.e. i i 1 … m f x i y i We can fit a smoothing spline: i.e. find f to minimize 1 n i 1 n y i f x i 2 λ u x 1 x m f x m u 2 The first term measures the closeness of the fitted function to the data, while the second penalizes the curvature in the function. λ established the trade off between the two. For 0 λ , this constraint is minimized by a natural spline of degree m 1 ( Schoenberg). If λ 0 , then f can be any function which interpolates the data. If λ and m 2 , then this is the simple least squares line fit, since no second derivative can be tolerated.