Fractional splines and wavelets

We extend Schoenberg's family of polynomial splines with uniform knots to all fractional degrees alpha > -1. These splines, which involve linear combinations of the one-sided power functions x(+)(alpha) = max(0,x)(alpha), are alpha-Holder continuous for alpha > 0. We construct the corresponding B-splines by taking fractional finite differences and provide an explicit characterization in both time and frequency domains. We show that these functions satisfy most of the properties of the traditional B-splines. including the convolution property, and a generalized fractional differentiation rule that involves finite differences only. We characterize the decay of the B-splines that are not compactly supported for nonintegral alpha's. Their most astonishing feature tin reference to the Strang-Fix theory) is that they have a fractional order of approximation alpha + 1 while they reproduce the polynomials of degree inverted right perpendicular alpha inverted left perpendicular. For alpha > -1/2, they satisfy all the requirements for a multiresolution analysis of L-2 (Riesz bounds, two-scale relation) and may therefore be used to build new families of wavelet bases with a continuously varying order parameter. Our construction also yields symmetrized fractional B-splines which provide the connection with Duchon's general theory of radial (m, s)-splines (including thin-plate splines). In particular, we show that the symmetric version of our splines can be obtained as the solution of a variational problem involving the norm of a fractional derivative.