Next Article in Journal
A Mixed Finite Volume Element Method for Time-Fractional Damping Beam Vibration Problem
Previous Article in Journal
Exact Solutions for the KMM System in (2+1)-Dimensions and Its Fractional Form with Beta-Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Some Results on a New Refinable Class Suitable for Fractional Differential Problems

1
Department of Basic and Applied Sciences for Engineering (SBAI), University La Sapienza of Rome, 00161 Rome, Italy
2
Facoltá di Scienze della Comunicazione, University of Teramo, 64100 Teramo, Italy
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(9), 521; https://doi.org/10.3390/fractalfract6090521
Submission received: 22 June 2022 / Revised: 28 July 2022 / Accepted: 2 August 2022 / Published: 15 September 2022
(This article belongs to the Topic Fractional Calculus: Theory and Applications)

Abstract

:
In recent years, we found that some multiscale methods applied to fractional differential problems, are easy and efficient to implement, when we use some fractional refinable functions introduced in the literature. In fact, these functions not only generate a multiresolution on R , but also have fractional (non-integer) derivative satisfying a very convenient recursive relation. For this reason, in this paper, we describe this class of refinable functions and focus our attention on their approximating properties.

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 R , 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 φ α , h , for α > 2 . The novelty of this paper, is that here we prove that these properties are also valid for α > 1 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 α > 2 to α > 1 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 α γ > 1 .
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 R 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:
c D t γ y ( t ) : = J ( k γ ) y ( k ) ( t ) ,         k 1 < γ < k ,         k N ,         t > 0 ,
where J ( β ) is the Riemann–Liouville integral operator
J ( β ) y ( t ) : = 1 Γ ( β )   0 t   y ( τ )   ( t τ ) β 1   d τ       β R ,
and Γ denotes Euler’s gamma function
Γ ( β ) : = 0   τ β 1   e τ   d τ   .
Hence,
c D t γ   y ( t ) : = 1 Γ ( k γ )   0 t   y ( k ) ( τ )   ( t τ ) k γ 1   d τ ,         k = γ .
For example, if γ = 0.5 then k = 1 and:
c D t 0.5 y ( t ) : = 1 Γ ( 0.5 )   0 t   y ( τ ) ( t τ )   d τ   .
If, for example, y ( t ) = t n then
c D t γ y ( t ) : = Γ ( n + 1 γ ) Γ ( n + 1 ) t n γ .
Riemann–Liouville definition is instead
R L D γ   y ( t ) : = d k d t k   ( γ ) y ( t )   ,         t > 0   .
They both reduce to the usual differential operator when γ N . In the general case, we have the following relation between the Caputo and the Riemann derivatives
c D γ   y ( t ) = R L D γ ( y ( t ) l = 0 k t l l ! y ( l ) ( 0 + ) ) .
The definitions coincide for homogenous boundary initial conditions.
In the Fourier domain one has
F D γ   , y ( t ) = ( i ω ) γ F ( y ( t ) )   ,       γ R +   , ω C
where F ( y ) is the Fourier transform of the function y.

3. MRA and Refinable Spaces

A sequence of functional spaces { V j , j Z } L 2 ( R ) , forms a multiresolution analysis (MRA) of L 2 ( R ) if
1.
V j V j + 1 , j Z ,
2.
j Z V j ¯ = L 2 ( R ) ;
3.
j Z V j = { 0 } ;
4.
f ( t ) V j f ( 2 t ) V j + 1 , j Z ;
5.
there exists a L 2 ( R ) -stable basis in V 0 .

MRA Based on Refinable Functions

An MRA can be generated by a refinable function ϕ, i.e., a function that satisfies a refinement functional equation
ϕ ( t ) = k Z   a k   ϕ ( 2   t k )   ,         t R   .
It is known that if the mask coefficients { a k ,   k Z } form a finite sequence and have some particular properties, then the existence of a unique solution to (10) in L 2 ( R ) , can be proved [12]. Moreover, the shifted refinable functions { ϕ ( t k ) , k Z } give rise to a stable basis in V 0 , i.e., the space they span.
It is important to associate (10) with its symbol
b n ( z ) = k a k   z k
When the mask is an infinity sequence, under suitable conditions the solution exists as proved in [8].
Now, we can define the spaces V j of the multiresolution:
V j : = span   { ϕ j k ( t ) : = ϕ ( 2 j   t k )   , k Z } ¯   ,         j Z   ,         t R   .
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 [ 0 , T ] , belonging to Ł 2 ( [ 0 , T ] ) .
Let us suppose that the support of ϕ is compact, i.e., s u p p   ϕ = [ 0 , σ ] . Then, we can define an MRA on the interval.
V j 0 [ 0 , T ] = span   { ϕ j k 0 ( t )   , k N j }   ,       j j 0   ,       t [ 0 , T ]   ¯ ,
where
ϕ j k 0 ( t ) : = { ϕ j k | [ 0 , T ] :   ϕ j k ( 0 ) = ϕ j k ( 0 ) = = ϕ j k ( n 1 ) ( 0 ) = 0 } ;
N j Z , with # N j = N j = 2 j + σ 1 , is the set of admissible index k and j 0 is the initial multiresolution scale, i.e., the minimal index such that supp   ϕ j 0 0 0 [ 0 , σ ] .

4. The Collocation Method and the Galerkin Method

We use the MRA to approximate each fractional differential problem by the collocation method and Galerkin method. For both methods, we pose the solution in the following form:
y j ( t ) = k N j   c j k   ζ j k ( t )
where ζ j k is a refinable function generating an MRA.
  • In the collocation method, we substitute (13) in the differential problem and we collocate it in the dyadic nodes { t p = p / 2 s , p = 0 , , N s , s j } . So, we obtain a linear algebraic system in N s equations and N j unknowns { c j k } k = 1 N j . Usually, we solve the system by least-squares method.
  • In the Galerkin method, we rewrite the differential problem in a weak form, and we substitute (13), using ζ j k as trial and test functions. In this way, the resulting linear algebraic system, will contain as the coefficient matrix, the integrals between u and the test (trial) functions ζ j k (Stiffness matrix).
In this way, the differential problem is converted into a system of algebraic equations that is suitable for computer programming.
Note.
If u also depends on x, i.e., u ( t , x ) , then the coefficients are c j k = c j k ( x ) [5,6].

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 [ 0 , n + 1 ] , in each interval of the partition are polynomials of degree at most n and in R have regularity C n 1 ( R ) . The Fourier transform of the classical B-Splines is:
B ^ n ( ω ) = ( 1 e i ω i ω ) n + 1 ,         n = 0 , 1 ,
We can define a fractional B-Spline starting with its Fourier transform obtained introducing a fractional (non-integer) exponent in (14):
B ^ α ( ω ) = ( 1 e i ω i ω ) α + 1 ,     α > 1
It is proven that for α > 1 , the antitransform B α is in L 1 ( R ) , while B α is in L 2 ( R ) for α > 1 / 2 [9].
In the time domain, the cardinal B-Splines B n , are defined in the following way. Let t + : = max ( 0 , t ) be the usual truncated power function and the finite difference operator
n   v ( t ) : = k N 0   ( 1 ) k   ( n k )   v ( t   k )   ,         n N
Then, B n ( t ) can be defined as:
B n ( t ) : = n + 1   t + n ( n + 1 ) ! ,
whose symbol is
b n ( z ) = 1 2 n ( 1 + z ) n + 1
In the non-integer case, we define the generalized finite difference operator
γ   v ( t ) : = k N 0   ( 1 ) k   ( γ k )   v ( t   k )   ,       γ R + .
When γ N , { ( γ k ) } k is a compactly supported sequence and we get the usual finite difference operator.
On the other hand, when γ R + \ N , then
( γ k ) : = Γ ( γ + 1 ) k !   Γ ( γ k + 1 ) = O ( k γ 1 )   ,       k N 0 ,     γ R +
and thus the sequence { ( γ k ) } k 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:
B α ( t ) : = α + 1   t + α Γ ( α + 1 )   ,       α > 1 2   ,
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,
D γ B n ( x ) = n + 1 t + n γ ( x ) Γ ( n + 1 γ ) = γ B n γ ( x )
In fact, one has
D γ B n ( x ) = D γ n + 1 t + n ( x ) ( n + 1 ) ! = n + 1 D γ t + n ( x ) ( n + 1 ) ! =
Γ ( n + 1 ) Γ ( n + 1 γ ) n + 1 t + n γ ( x ) ( n + 1 ) ! = n + 1 t + n γ ( x ) Γ ( n + 1 γ )
Proof. 
Now, for the rule of the difference finite operator composition
γ n γ + 1 = n + 1 ,
it is easy to verify that
n + 1 t + n γ ( x ) Γ ( n + 1 γ ) = γ n γ + 1 t + n γ ( x ) Γ ( n + 1 γ ) = γ B n γ ( x ) .
The theorem is proved.  ☐
It is also worthwhile to define the symbol b α of B α , i.e.,
b α ( z ) = 1 2 α ( 1 + z ) α + 1

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 B α belong to L 1 ( R ) , for α > 1 and to the Sobolev space W 2 r ( R ) , 0 r < α + 1 2 , for α > 1 2 ; where W 2 r ( R ) , represents the Banach subspace of L 2 ( R ) , equipped with the norm
| | f | | r = | | f | | 2 + | | D r f | | 2
Proposition 2.
When α > 1 / 2 , 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 L 2 ( R )
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 α + 1 , instead of α + 1 .
For the CAGD and isogeometric context, it is important to know that they form a partition of unity for α > 1 .
Proposition 4.
It is also important to consider the following fractional derivation rule that is a generalization of the Formula (20)
D γ   ( B α ) = γ   B α γ
where D γ 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 L 2 ( R ) , of order of approximation α. Then, ϕ α can be factorized as
ϕ α = B α ϕ 0 ,
α 0 and ϕ 0 is a distribution such that ϕ 0 = 1 [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 ϕ 0 C 0 ( R ) , then ϕ α C α ( R ) and
D γ   ϕ α = D γ   ( B α ϕ 0 ) = γ   ( B α γ ϕ 0 ) = γ ϕ α γ   ,         0 < γ α   .
The claim follows from some results in [10].
For shifted functions ϕ α , k ( t ) , we obtain a similar result.
Proposition 6.
Let ϕ α , k ( t ) : = ϕ α ( t k ) . Then,
D t γ   ϕ α , k = γ   ϕ α γ , k   ,         0 < γ α   .
Let us note that since ϕ L 2 ( R ) and the generalized binomial coefficients decay similar to k γ 1 as k + , thus the series in (24) converges. Thus, in practical computation, γ ϕ α γ , k 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 [ 0 , 2 ] ) introduced in [11]. We consider
ϕ α , h = 1 2 B α 2 ϕ 1 , h ^   ,       0 α h   ,
where ϕ h ^ L 2 ( R ) is the elementary refinable function, solution of the refinement equation
ϕ 1 , h ^ ( t ) = k = 0 2   a h ^ , k   ϕ h ( 2 t k )   ,       t R   ,
with mask coefficients in [11] and h ^ = h α + 1 . h is a real shape parameter that controls the shape of ϕ α , h . The symbol of ϕ n , h in general is
b h n ( z ) = 1 2 h [ ( 1 + z ) n + 1 + 4 ( 2 h n 1 ) z ( 1 + z ) n 1 ] .
that, for n = 1 reduces to
b h 1 ( z ) = 1 2 h [ ( 1 + z ) 2 + 4 ( 2 h 1 1 ) z ] .
In the Fourier domain, the definition of ϕ α , h becomes:
F ( ϕ α , h ) ( ω ) = ( 1 e i ω i ω ) α 1 ϕ ^ h ^ ( ω )
We observe that when α N , α 0 , then ϕ α , h is compactly supported, belongs to C α 1 ( R ) and is a GP function as in [11]; in particular for h = α it reduces to a cardinal B-Spline. Instead, when α is not an integer but h = α , then ϕ α , α is a fractional B-spline in [9].
It is easy to show that ϕ α , h can be also obtained by placing a fractional index in the mask of ϕ n , h , i.e., a α , h , k = 1 h ( α + 1 k ) + 4 ( 2 h α 1 ) ( α 1 k 1 ) and, in this case ϕ α , h becomes:
b h α ( z ) = 1 2 h [ ( 1 + z ) α + 1 + 4 ( 2 h α 1 ) z ( 1 + z ) α 1 ]
Therefore, it is not difficult to prove that:
b h α ( z ) = 1 2 b α 2 ( z )   b h ^ 1 ( z )
where
b α 2 ( z ) = 1 2 α 2 ( 1 + z ) α 1     a n d   h ^ = h α + 1
In fact,
b h α ( z ) = 1 2 h ( 1 + z ) α 1 [ ( 1 + z ) 2 + ( 2 h α + 2 ) 2 2 ) z ] =
= 1 2 α 2 1 2 h α + 2 ( 1 + z ) α 1 [ ( z 2 + ( 2 h α + 2 ) 2 ) z + 1 ] .
Observe that from (29), we deduce that B α 2 carries all the approximation properties of φ α , h . In fact, since ϕ h ^ is summable, the convolution preserves all the properties of B α 2 . So, we have the following theorem,
Theorem 2.
For any admissible α and h, ϕ α , h belongs to C α 2 ( R ) (and decays to the infinity rather rapidly so that in practice it can be assumed compactly supported).
Moreover, it has derivative α 1 , but it is only bounded, not necessary continuous; one says that it is α- continous. As for the order of approximation,   ϕ α , h has order of approximation α 1 and order of polynomial reproducibility α 1 ; so it does not verify the Strang and Fix theory.
Finally, the differentiation rule is specified in
D t γ   ϕ α , h ( t ) = γ   ϕ α γ , h α + 2 ( t ) = k N 0   ( 1 ) k   ( α k )   ϕ α γ , h α + 2 ( t k ) ,     0 < γ α   .
Proof. 
The properties of ϕ α , h are the same properties of B α 2 [9] that are preserved through the convolution Formula (27) since ϕ 1 , h ^ 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 ϕ h α with respect to the GP refinable function, is that, in practice, due to the rapid decay of ϕ h α , their supports appear strictly contained in the supports [ 0 , n + 1 ] of ϕ n , h , but the order of exactness is the same, i.e., n 1 . 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 ϕ h α , with α = 1.5 , then the order of polynomial reproducibility is α 1 = 1 , that is the straight line can be reproduced, in the same manner as classical GP refinable, when n = 1 , and support [ 0 , 2 ] .

Author Contributions

Conceptualization, L.P. and L.T.; methodology, L.P.; software, L.T.; validation, L.P. and L.T.; formal analysis, L.P.; investigation, L.P.; resources, L.P.; data curation, L.P. and Luca Tallini; writing—original draft preparation, L.P. and L.T.; writing—review and editing, L.P. and L.T.; visualization, L.P. and L.T.; supervision, L.P. and L.T.; project administration, L.P. and Luca Tallini; funding acquisition, L.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We would like to thank Francesca Pitolli and our daughter Irene Tallini for all the precious suggestions they gave us.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kimeu, J.M. Fractional Calculus: Definitions and Applications. Ph.D. Thesis, Western Kentucky University, Bowling Green, KY, USA, 2009. [Google Scholar]
  2. Lino, P.; Maione, G.; Garrappa, R.; Holm, S. An approach to optimal integer and fractional-order modeling of electro-injectors in compression-ignition engines. Control. Eng. Pract. 2021, 115, 104890. [Google Scholar] [CrossRef]
  3. Garrappa, R.; Popolizio, M. Evaluation of generalized Mittag–Leffler functions on the real line. Adv. Comput. Math. 2013, 39, 205–225. [Google Scholar] [CrossRef]
  4. Pedas, A.; Tamme, E. On the convergence of spline collocation methods for solving fractional differential equations. J. Comput. Appl. Math. 2011, 235, 3502–3514. [Google Scholar] [CrossRef]
  5. Pezza, L.; Pitolli, F. A fractional spline collocation-Galerkin method, for the fractional diffusion equation. CAIM 2018, 9, 104–120. [Google Scholar] [CrossRef]
  6. Pezza, L.; Pitolli, F. A multiscale collocation method for fractional differential problems. Math. Comput. Simul. 2018, 147, 210–219. [Google Scholar] [CrossRef]
  7. Pezza, L.; Pitolli, F. A Fractional Spline Collocation Method for the Fractional Order Logistic Equation. In Approximation Theory XV: San Antonio 2016, Proceedings in Mathematics and Statistics; Fasshauer, G., Schumaker, L., Eds.; Springer: Berlin/Heidelberg, Germany, 2017; Volume 201, pp. 307–318. [Google Scholar]
  8. Pezza, L. Fractional GP Refinable Functions. Rend. Mat. Ser. VII 2007, 27, 73–87. [Google Scholar]
  9. Unser, M.; Blu, T. Fractional splines and wavelets. SIAM Rev. 2000, 42, 43–67. [Google Scholar] [CrossRef]
  10. Unser, M.; Blu, T. Wavelet theory demystified. IEEE Trans. Sig. Proc. 2003, 51, 470–483. [Google Scholar] [CrossRef]
  11. Gori, L.; Pitolli, F. A class of totally positive refinable functions. Rend. Mat. Ser. VII 2000, 20, 305–322. [Google Scholar]
  12. Dahmen, W.; Micchelli, C.A. Using the refinement equation for evaluating integrals of wavelets. SIAMJ Numer. Anal. 1993, 30, 507–537. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pezza, L.; Tallini, L. Some Results on a New Refinable Class Suitable for Fractional Differential Problems. Fractal Fract. 2022, 6, 521. https://doi.org/10.3390/fractalfract6090521

AMA Style

Pezza L, Tallini L. Some Results on a New Refinable Class Suitable for Fractional Differential Problems. Fractal and Fractional. 2022; 6(9):521. https://doi.org/10.3390/fractalfract6090521

Chicago/Turabian Style

Pezza, Laura, and Luca Tallini. 2022. "Some Results on a New Refinable Class Suitable for Fractional Differential Problems" Fractal and Fractional 6, no. 9: 521. https://doi.org/10.3390/fractalfract6090521

APA Style

Pezza, L., & Tallini, L. (2022). Some Results on a New Refinable Class Suitable for Fractional Differential Problems. Fractal and Fractional, 6(9), 521. https://doi.org/10.3390/fractalfract6090521

Article Metrics

Back to TopTop