1. Introduction
In the last decades, fractional calculus has increased in popularity, owing to the awareness that many physical problems, such as viscoelasticity, Brownian motion, medical issues and so forth, require fractional derivatives to be modeled appropriately. For a better understanding, please see [
1,
2].
Analytical solutions for certain problems have been found. They are expressed through the Mittag–Leffler function [
3], which is a series expansion, and thus require numerical tools to be computed. For this issue and for the other unsolved problems, the literature provides many ways to numerically solve fractional differential problems. Most of the methods employ the quadrature rule to compute the fractional derivatives [
4]; others use spectral or Galerkin methods [
5].
In recent papers [
5,
6,
7], the authors proved that the multiscale collocation methods are easy and efficient to implement, when using certain fractional refinable functions introduced in [
6,
8]. In fact, these functions not only generate a multiresolution on
, but also satisfy a fractional derivative convenient formula. Moreover, the collocation technique allows one to obtain an algebraic system from a differential problem.
The coefficient matrix is given by the collocation of basis functions into the collocation nodes. In this way, the result is given by the solution (often in a least-squares sense) of a linear algebraic system. The goal of this paper is to prove further approximating properties of this class of fractional refinable functions with respect to [
6,
8], suitable to the solution of fractional differential problems.
More precisely, in [
6,
8], we proved the basis properties of the class
, for
. The novelty of this paper, is that here we prove that these properties are also valid for
and that other important approximating and smoothing properties can be proved, e.g., the order of polynomial reproducibility. In this way, we enlarge the class of fractional refinable functions from
to
and thus, also its applicability to
a wider class of fractional differential problems. Furthermore, we prove that all the properties derive from a suitable convolution formula. Note that when we apply these functions to a differential problem with fractional derivative
, we have to choose refinable functions of approximation order
such that
.
The paper is organized as folllows.
Section 3 introduces some fractional derivative definitions, that can be computed by numerical quadrature rules. We choose the Caputo derivative for several reasons: computational efficiency, minor regularity required, stability [
2].
Section 4 explains Multiresolution Analysis (MRA) properties on
and on the interval.
Section 5 describes the collocation and the Galerkin methods constructed with MRA.
Section 6 lists the main properties of the fractional B-splines, introduced in [
9,
10], emphasizing the fractional derivative properties.
Section 7 describes the new class of fractional refinable functions constructed introduced by [
6,
8], through a convolution formula involving the functions in [
11] and with a continuos dependence from a parameter
h. We prove that these functions satisfy new properties that are similar to those of the fractional B-splines, such as, for example, the polynomial reproducibility. Furthermore, we prove a differentiation formula that makes them particularly interesting in the fractional derivative context. In the conclusions, we explain all the advantages of this new class of fractional refinable functions, including an example on polynomial reproducibility.
2. Fractional Derivatives
The fractional derivative can be defined in many ways: for example, in the Caputo sense or in the Riemann Liouville way.
The Caputo definition of the
fractional derivative is:
where
is the
Riemann–Liouville integral operator
and
denotes Euler’s gamma function
For example, if
then
and:
If, for example,
then
Riemann–Liouville definition is instead
They both reduce to the usual differential operator when
. In the general case, we have the following relation between the Caputo and the Riemann derivatives
The definitions coincide for homogenous boundary initial conditions.
In the Fourier domain one has
where
is the Fourier transform of the function
y.
3. MRA and Refinable Spaces
A sequence of functional spaces , forms a multiresolution analysis (MRA) of if
- 1.
, ,
- 2.
;
- 3.
;
- 4.
, ;
- 5.
there exists a -stable basis in .
MRA Based on Refinable Functions
An MRA can be generated by a
refinable function ϕ, i.e., a function that satisfies a
refinement functional equationIt is known that if the mask coefficients
form a finite sequence and have some particular properties, then the existence of a unique solution to (
10) in
, can be proved [
12]. Moreover, the shifted refinable functions
give rise to a stable basis in
, i.e., the space they span.
It is important to associate (
10) with its symbol
When the mask is an infinity sequence, under suitable conditions the solution exists as proved in [
8].
Now, we can define the spaces
of the multiresolution:
Since we are taking into account differential problems of order n with initial conditions, it is also important to define an MRA on an interval , belonging to .
Let us suppose that the support of
is compact, i.e.,
. Then, we can define an MRA on the interval.
where
, with , is the set of admissible index k and is the initial multiresolution scale, i.e., the minimal index such that supp .
5. Fractional B-Splines
A particular class of refinable functions is provided by the cardinal B-Splines of degree n, i.e., functions that are positive and compactly supported in
, in each interval of the partition are polynomials of degree at most n and in
have regularity
. The Fourier transform of the classical B-Splines is:
We can define a
fractional B-Spline starting with its Fourier transform obtained introducing a fractional (non-integer) exponent in (
14):
It is proven that for
, the antitransform
is in
, while
is in
for
[
9].
In the time domain, the cardinal B-Splines
, are defined in the following way. Let
be the usual truncated power function and the finite difference operator
Then,
can be defined as:
whose
symbol is
In the
non-integer case, we define the
generalized finite difference operatorWhen , is a compactly supported sequence and we get the usual finite difference operator.
On the other hand, when
, then
and thus the sequence
is absolutely summable and the limit of the series (
18) exists under suitable hypothesis on
v [
9].
The
fractional B-spline, i.e., the B-spline of non-integer order, in the time domain is defined as:
The following theorem writes, with a different proof with respect to [
9].
Theorem 1. The fractional derivative of a B-Spline is a fractional B-Spline. More precisely, Proof. Now, for the rule of the difference finite operator composition
it is easy to verify that
The theorem is proved. ☐
It is also worthwhile to define the symbol
of
, i.e.,
5.1. Main Properties of Fractional B-Splines
In the study by [
9], fractional B-splines are introduced for the first time and their main properties are proved. We summarize these properties in the following propositions avoiding the proof.
Proposition 1. The fractional B-Splines belong to , for and to the Sobolev space , for ; where , represents the Banach subspace of , equipped with the norm Proposition 2. When , the fractional B-Splines are α-order continuous, i.e., they can be derived up to the order α but is in general only bounded.
Moreover, they generate an MRA of
Proposition 3. The fractional B-splines reproduce polynomials up to degree , but they do not satisfy Strang and Fix theory. In fact, they have fractional approximation order , instead of .
For the CAGD and isogeometric context, it is important to know that they form a partition of unity for
Proposition 4. It is also important to consider the following fractional derivation rule that is a generalization of the Formula (20)where is the usual derivative of order γ. There is also a formula that allows us to assume that a fractional B-spline preserves the order of approximation of any refinable function of order .
Proposition 5. Let be a refinable function generating an MRA in , of order of approximation α. Then, can be factorized as and is a distribution such that [10]. Let us observe that all the previous propositions can be proved by starting from Proposition 5.
5.2. Fractional Derivative of Refinable Functions
If we consider a generic function
of order
, it is possible to generalize the differentiation rule (
21).
In fact, let it be that
, then
and
The claim follows from some results in [
10].
For shifted functions , we obtain a similar result.
Proposition 6. Let . Then, Let us note that since
and the generalized binomial coefficients decay similar to
as
, thus the series in (
24) converges. Thus, in practical computation,
is a finite sum.
6. Fractional GP Refinable Functions
We present here the main results regarding a new class of refinable functions of fractional order
, obtained starting by a suitable refinable function (of support
) introduced in [
11]. We consider
where
is the
elementary refinable function, solution of the refinement equation
with mask coefficients in [
11] and
.
h is a real
shape parameter that controls the shape of
. The symbol of
in general is
that, for
reduces to
In the Fourier domain, the definition of
becomes:
We observe that when
,
, then
is compactly supported, belongs to
and is a GP function as in [
11]; in particular for
it reduces to a cardinal B-Spline. Instead, when
is not an integer but
, then
is a fractional B-spline in [
9].
It is easy to show that
can be also obtained by placing a fractional index in the mask of
, i.e.,
and, in this case
becomes:
Therefore, it is not difficult to prove that:
where
Observe that from (29), we deduce that carries all the approximation properties of . In fact, since is summable, the convolution preserves all the properties of . So, we have the following theorem,
Theorem 2. For any admissible α and h, belongs to (and decays to the infinity rather rapidly so that in practice it can be assumed compactly supported).
Moreover, it has derivative , but it is only bounded, not necessary continuous; one says that it is α- continous. As for the order of approximation, has order of approximation and order of polynomial reproducibility ; so it does not verify the Strang and Fix theory.
Finally, the differentiation rule is specified in Proof. The properties of
are the same properties of
[
9] that are preserved through the convolution Formula (
27) since
is summable. ☐
7. Conclusions
Since, as in the classical B-spline case, the fractional derivative of a GP refinable function is a GP fractional refinable function, we deal in this paper with fractional GP functions stemming from the fractional derivative of GP refinable functions. In this way, we obtain a class of refinable functions, closed with respect to the fractional derivative.
Another advantage of these fractional GP refinable functions
with respect to the GP refinable function, is that, in practice, due to the rapid decay of
, their supports appear strictly contained in the supports
of
, but the order of exactness is the same, i.e.,
. This property, in addition to derivative Formulas (
23) and (
24), renders them highly suitable for solving fractional differential problems, as shown in [
5,
6].
More precisely, if, for example, we consider , with , then the order of polynomial reproducibility is , that is the straight line can be reproduced, in the same manner as classical GP refinable, when , and support .