Next Article in Journal
Isogeometric Schwarz Preconditioners with Generalized B-Splines for the Biharmonic Problem
Next Article in Special Issue
Families of Orbits Produced by Three-Dimensional Central and Polynomial Potentials: An Application to the 3D Harmonic Oscillator
Previous Article in Journal
The Difference of Zagreb Indices of Halin Graphs
Previous Article in Special Issue
Nonuniform Dichotomy with Growth Rates of Skew-Evolution Cocycles in Banach Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Dynamical Systems Solved by a Collocation Method Based on Refinable Spaces

by
Laura Pezza
1,*,† and
Simmaco Di Lillo
2,*,†
1
Department of Basic and Applied Sciences for Engineering (SBAI), Università di Roma “La Sapienza”, Via A. Scarpa 16, 00161 Rome, Italy
2
Department of Mathematics Guido Castelnuovo, Università di Roma “La Sapienza”, Piazzale Aldo Moro 5, 00185 Rome, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Axioms 2023, 12(5), 451; https://doi.org/10.3390/axioms12050451
Submission received: 24 March 2023 / Revised: 28 April 2023 / Accepted: 29 April 2023 / Published: 3 May 2023

Abstract

:
A dynamical system is a particle or set of particles whose state changes over time. The dynamics of the system is described by a set of differential equations. If the derivatives involved are of non-integer order, we obtain a fractional dynamical system. In this paper, we considered a fractional dynamical system with the Caputo fractional derivative. We collocated the fractional differential problem in dyadic nodes and used refinable functions as approximation functions to achieve a good degree of freedom in the choice of the regularity. The collocation method stands out as a particularly useful and attractive tool for solving fractional differential problems of various forms. A numerical result is presented to show that the numerical solution fits the analytical one very well. We collocated the fractional differential problem in dyadic nodes using refinable functions as approximation functions to achieve a good degree of freedom in the choice of regularity.

1. Introduction

For a long time, fractional calculus was considered to be a purely mathematical field with no real applications. Recently, however, this has changed. Think, for example, of all the problems where it is important to describe the behaviour and evolution of the function over a period of time and not just at a specific point in time. Think, for example, of materials with memory. Fractional calculus is, therefore, particularly important for studying and solving this type of problem.
Many researchers have found that fractional derivatives are suitable for mathematical modelling of physical problems such as diffusion in biological tissues [1,2,3], propagation in porous media [4], anomalous diffusion [5,6,7], viscoelasticity [8,9,10,11], and earthquakes [12]. Fractional calculus has also been considered for pandemic models [13,14,15] or in finance [16].
Since the fractional derivate involves an integral, many authors use quadrature formulae or discrete derivatives [17] to approximate the fractional derivative. Others use Galerkin methods [18,19] or spectral methods [20].
In our previous works [21,22], we proved that if we collocate the fractional differential problem in dyadic nodes and we use the B-splines as approximation functions, we obtain a collocation method that works very well. Moreover, the method takes advantage of the refinability and derivability properties of the B-splines in order to produce an accurate and efficient algorithm. Indeed, the fractional derivatives of the B-splines appearing in the collocation matrix satisfy a closed form involving fractional B-splines of lower orders. In this paper, we used this collocation method, generalised to a wider class of approximating refinable functions [23], and we used this method to numerically solve a fractional dynamical system. By using these fractional functions instead of fractional B-splines, we have more degrees of freedom in the choice of the regularity. For all these reasons, the collocation method stands out as a particularly useful and attractive tool for the solution of fractional differential problems of different forms [24].
The paper is structured as follows: In the first section, we describe the fractional dynamical systems and we define the fractional derivatives in several meaningful ways. In the second section, fractional refinable functions are introduced, as well as their main properties. The most-important property in the context of this paper, i.e., the fractional differential rule, is treated in the third section, as well as the way to differentiate the left-edge B-splines in the fractional sense. In the following section, a Multiresolution Analysis (MRA) on the semi-bounded interval [ 0 , + ) is introduced. The fractional collocation method is then introduced, specifying the entries of the collocation matrix, and a numerical test is provided in the last section. We compare one of the other results obtained using several B-splines and fractional B-splines in [18,21,25,26].

2. Fractional Dynamical Systems

Consider the following fractional dynamical system:
D t γ X ( t ) = A ( t ) X ( t ) , γ ( 0 , 1 ) , t > 0 X ( 0 ) = X 0
where X ( t ) = ( x 1 ( t ) , x 2 ( t ) , , x n ( t ) ) R n , and A ( t ) R n × n .
In this context, D t γ y denotes the Caputo fractional derivative with respect to the time t defined as
D t γ y ( t ) : = J ( k γ ) y ( k ) ( t ) , γ > 0
where k is an integer such that k 1 < γ < k and J ( γ ) is the Riemann–Liouville integral operator given by
J ( γ ) y ( t ) : = 1 Γ ( γ ) 0 t y ( τ ) ( t τ ) γ 1 d τ , t > 0 .
In the above equation, Γ denotes Euler’s gamma function:
Γ ( γ ) = 0 + x γ 1 e x d x for   γ > 0 , Γ ( γ ) = γ 1 Γ ( γ + 1 ) for   γ < 0 .
Through the Gamma function, it is possible to define the generalised binomial coefficients:
γ l = Γ ( γ + 1 ) Γ ( l + 1 ) Γ ( γ l + 1 ) , l Z .
If γ is a positive integer, the above definition and the usual definition of the binomial coefficient coincide, while they are infinitely supported if γ R + \ N or zero if l Z \ N .
However, they decay to infinity as [27]
γ l = O ( l γ 1 ) for l .
We note that, when the homogeneous initial conditions are imposed on the function y ( t ) , the Caputo definition (2) coincides with the Riemann–Liouville definition:
D t γ y ( t ) : = d k d t k J ( γ ) y ( t ) , t > 0 ,
which requires less regularity on y ( t ) , and both reduce to the usual differential operator in the case of γ N [28].
In the following, we will make use of D t γ f ( 2 j t ) , and so, we will need the following theorem, already proven in [22].
Theorem 1.
Let γ be a real number such that 0 < γ < 1 then
D t γ f ( 2 j t ) = 2 j γ D ( 2 j t ) γ f ( 2 j t ) .
Proof. 
Denote F ( t ) = f ( 2 j t ) , then F ( m ) ( t ) = 2 j m f ( m ) ( 2 j t ) , m N . Therefore,
D t γ F ( t ) = 1 Γ ( k γ ) 0 t F ( k ) ( τ ) ( t τ ) ( γ k + 1 ) d τ = = 1 Γ ( k γ ) 2 j k 0 t f ( k ) ( 2 j τ ) ( t τ ) ( γ k + 1 ) d τ .
After changing the variable in the integral, we obtain
D t γ F ( t ) = 1 Γ ( k γ ) 2 j k 0 2 j t f ( k ) ( τ ) ( t 2 j τ ) ( γ k + 1 ) 2 j d τ = = 1 Γ ( k γ ) 2 j k 2 j 2 j ( γ k + 1 ) 0 2 j t f ( k ) ( τ ) ( 2 j t τ ) ( γ k + 1 ) d τ = 2 j γ D ( 2 j t ) γ f ( 2 j t ) .
In the next section, we introduce a new collocation method, the multi-scale collocation method, which takes advantage of multi-scale techniques to produce an efficient and accurate algorithm.

3. Fractional Cardinal B-Splines and Fractional GP Functions

In the paper [27], the extension of the class of the cardinal B-splines to the fractional ones is provided by introducing the concept of the fractional finite difference operator. The starting point is the definition of the classical cardinal B-spline, i.e., the B-spline on integer knots. Consider the truncated power function:
T n ( x ) : = ( max ( 0 , x ) ) n , x R , n N
and the finite difference operator:
Δ n f ( x ) : = k = 0 n n k ( 1 ) k f ( x k ) , x R .
Then, the classical B-spline of integer-order n is defined by the formula:
B n ( x ) : = 1 ( n + 1 ) ! Δ n + 1 T n ( x ) .
From this definition, the Fourier transform of B n becomes
B ^ n ( ω ) = 1 e i ω i ω n + 1 .
Due to the definition of generalised binomial coefficients, the definition (4) can be extended to a real index α . Hence, using the fractional finite difference operator:
Δ γ f ( x ) : = k 0 γ k ( 1 ) k f ( x k ) ,
we can define the fractional B-spline of order α :
B α ( x ) : = 1 Γ ( α + 1 ) Δ α + 1 T α ( x ) .
It is also easy to prove that
Δ γ B n ( x ) = Δ n + 1 T n γ ( x ) Γ ( n γ + 1 ) .
Thus, due to the finite fractional differential rule [22]:
D γ B n ( x ) = Δ γ B n γ ( x )
we can say that the fractional derivative of a B-spline is a fractional B-spline.
It easy to prove that the Fourier transform of B α becomes
B ^ α ( ω ) = 1 e i ω i ω α + 1 .
A similar theory can be developed for the fractional derivatives of a class of refinable functions compactly supported in [ 0 , n + 1 ] , called Gori–Pitolli (GP) refinable functions, depending on a real parameter h n . These functions are denoted by φ n , h . Their fractional derivatives are fractional GP refinable functions, i.e., GP refinable functions with non-integer-order α : φ α , h (see Figure 1). See [23,25,26,29] for more details.
In Figure 1a, the fractional GP functions are plotted for α = 2:0.25:6 and h = α + 3 . It is possible to see the continuous dependence of the family on the index α and the interpolation with the fractional B-splines of [27] when α = h .
Instead, when h > > α , each B-spline looks like the corresponding one of degrees less than 2.
The same observations can be made for Figure 1b, where h = α + 0.5 ; here, the fractional functions are more similar to the fractional B-splines.
In both cases, the support of the fractional GP function looks like [ 0 , σ ] , where σ = γ + 1 .

4. Multiresolution Analysis on R and on [ 0 , )

Let { V j , j Z } be a sequence of closed subspaces of L 2 ( R ) . The sequence forms a multiresolution analysis (MRA) of L 2 ( R ) if:
(i)
V j V j + 1 , j Z ;
(ii)
j Z V j ¯ = L 2 ( R ) ;
(iii)
j Z V j = { 0 } ;
(iv)
f ( t ) V j f ( 2 t ) V j + 1 , j Z ;
(v)
There exists an L 2 ( R ) -stable basis in V 0 .
An MRA can be generated by a refinable function, i.e., a function defined through a refinement mask a = { a k R , k Z } and a refinement equation, as follows:
φ ( t ) = k Z a k φ ( 2 t k ) , t R .
Suitable conditions on the mask coefficients { a k } ensure the existence of a unique function φ solution to (6), which belongs to L 2 ( R ) and such that their integer translates { φ ( t k ) , k Z } form an L 2 ( R ) -stable basis in V 0 (see [30] for details). As a consequence, the refinable function φ generates all the spaces V j through dilation and translation, i.e.,
V j = span { φ j k ( t ) : = φ ( 2 j t k ) , k Z , t [ 0 , ) } , j 0 .
Let V j 0 [ 0 , + ) denote the restriction of V j to the semi-infinite interval [ 0 , + ) . Thus, V j 0 [ 0 , + ) is generated by suitable functions φ j k , i.e., with φ j k ( 0 ) = 0 :
V j 0 [ 0 , + ) = span { φ j k ( t ) , k Z , t [ 0 , ] } , j 0 .
It is possible to prove that V j 0 [ 0 , + ) still generate an MRA on [ 0 , + ) [22].

5. The Fractional Derivative of B-Splines

In this section, we want to analyse the application of the fractional relation (5) to the B-splines in I = [ 0 , ) , in particular for the left-edge functions. Let [ 0 , n + 1 ] be the support of the cardinal B-spline B n . Let us consider
B k n ( x ) : = B n ( x k ) .
Note that, for the B k n whose support is entirely contained in I, the index k is greater than or equal to zero, while the left-edge B-splines have index k = n , , 1 . Thus, the following two theorems state the adapted differentiation relation for both types. First, we state a theorem for the exact calculation of some integrals involved in the integration of the B-splines derivative. By (6), we have
B 0 n ( x ) = 1 n ! r = 0 n + 1 ( 1 ) r n + 1 r ( x r ) + n
and for the derivative,
d d x B 0 n ( x ) = 1 ( n 1 ) ! r = 0 n + 1 ( 1 ) r n + 1 r ( x r ) + n 1 .
If the index k = n , , 1 , we have the left-edge B-splines. To determine the fractional derivative of these functions, for x > 0 and n k 1 , we use the difference between the whole integral k D x α B k n ( x ) , which is known, and the part before zero, 0 D x α B k n ( x ) :
0 D x α B 0 n ( x ) : = 1 Γ ( 1 α ) 0 x B k ( t ) ( x t ) α d t = k D x α B k n ( x ) k D 0 α B k n ( x )
where
k D x α B k n ( x ) : = 1 Γ ( 1 α ) k x B k ( s + k ) ( x k s ) α d s .
By summing and subtracting the index k and using the change of the variable t k = s , we obtain the following equation:
k D 0 α B k n ( x ) = 1 Γ ( 1 α ) k 0 B k ( s + k ) ( x k s ) α d s = = 1 Γ ( 1 α ) 0 k B 0 ( s ) ( x k s ) α d s .
Therefore, (7) can be rewritten as
0 D x α B k n ( x ) = Δ n + 1 x + n α Γ ( n + 1 α ) x k 1 Γ ( 1 α ) 0 k B 0 ( s ) ( x k s ) α d s k = n , , 1 .
The calculus of the fractional derivatives can be reduced to evaluating k D 0 α B k n ( x ) .
Theorem 2.
For 0 α 1 , for the B-splines with all the support in I = [ 0 , T ] , one has
k D x α B k n ( x ) = 0 D x k α B 0 n ( x k ) , k = 0 , , b a ,
while for the left B-splines,
0 D x α B k n ( x ) = Δ n + 1 ( x ) + n α Γ ( n + 1 α ) x k 1 Γ ( 1 α ) 0 k B 0 ( s ) ( x k s ) α d s , k = n , , 1
whereby the second term has the following recurrence relation:
1 Γ ( 1 α ) 0 k B 0 ( s ) ( x k s ) α d s = = 1 Γ ( n + 1 α ) r = 0 k 1 ( 1 ) r n + 1 r ( x k r ) ( n α ) + x 1 α p = 0 n 1 ( 1 ) n p ( k r ) n 1 p ( x k r ) p ( n 1 p ) ! s = 1 n 1 p ( α s ) .
Proof. 
For the B-splines contained in [ 0 , T ] , we have
k D x α B k n ( x ) = 1 Γ ( 1 α ) k x B k ( t ) ( x t ) α d t .
By summing and subtracting k and by setting t k = s , we have
k D x α B k n ( x ) = 1 Γ ( 1 α ) 0 x k B k ( s + k ) ( x k s ) α d s .
Since B k ( x ) = B 0 ( x k ) , we have
k D x α B k n ( x ) = 1 Γ ( 1 α ) 0 x k B 0 ( s ) ( x k s ) α d s = 0 D x k α B 0 n ( x k ) ;
This is the first thesis.
As for the second thesis, we remember that B 0 ( x ) is written as
B 0 ( x ) = d d x Δ n + 1 x + n n ! = 1 ( n 1 ) ! r = 0 n + 1 ( 1 ) r n + 1 r ( x r ) + n 1 .
Putting it in (7) and setting s r = t , we obtain
k D 0 α B k n ( x ) = 1 Γ ( 1 α ) ( n 1 ) ! r = 0 k 1 ( 1 ) r n + 1 r r k r t + n 1 ( x k r t ) α d t = 1 Γ ( 1 α ) ( n 1 ) ! r = 0 k 1 ( 1 ) r n + 1 r 0 k r t + n 1 ( x k r t ) α d t
which is the sum of only the first k terms.
The integration rules for rational functions are
0 k t n 1 ( x k t ) α d t = ( n 1 ) ! s = 1 n ( α s ) ( x k t ) ( 1 α ) r = 0 n 1 ( k x ) r ( t ) n 1 r ( n 1 r ) ! s = 1 n 1 r ( α s ) t = 0 k = ( n 1 ) ! s = 1 n ( α s ) ( x k ) ( n α ) + x 1 α r = 0 n 1 ( 1 ) n r ( k ) n 1 r ( x k ) r ( n 1 r ) ! · · s = 1 n 1 r ( α s ) .

6. The Fractional Collocation Method

Let j be fixed, and for simplicity, we ignore the dependence of j.
In the theory, we are looking for a vector of approximating functions in V j [ 0 , ) , i.e.,
x ˜ i ( t ) = k Z c i k φ j k ( t ) , i = 1 , , n
which solves the differential problem (1) on a set of collocation points.
In practice, since φ j k have a fast decay, the series reduces to a finite number of terms:
x i ( t ) = k = K 1 K 2 c i k φ j k ( t ) , i = 1 , , n
where K 1 and K 2 are given by
K 1 = min { k Z | supp ( φ j k ) [ 0 , 1 ] }
and
K 2 = max { k Z | supp ( φ j k ) [ 0 , 1 ] } .
Let
C = c 1 K 1 c 1 K 2 c n K 1 c n K 2
be the coefficient matrix and
Φ ( t ) = φ j K 1 ( t ) φ j K 2 ( t ) .
Therefore, we can rewrite (8) in matrix form:
X = C Φ ( t ) .
Using the linearity of D γ and (8), we obtain
D γ X ( t ) = C D ( t )
where
D ( t ) = D γ φ j K 1 ( t ) D γ φ j K 2 ( t ) .
If we choose as collocation points the dyadic nodes that are defined by { t p = p / 2 s , p = 0 , , 2 s } , then
D γ X ( t p ) = A ( t p ) X ( t p ) , p = 0 , , 2 s X ( 0 ) = X 0 .
By (9) and (10), (11) becomes
C D ( t p ) = A ( t p ) C Φ ( t p ) , p = 0 , , 2 s C Φ ( 0 ) = X 0 .
Using the Kronecker product, we obtain the following linear algebraic system:
( D ( t p ) T I n Φ ( t p ) T A ( t p ) ) vec ( C ) = 0 , p = 0 , , 2 s ( Φ ( 0 ) T I n ) vec ( C ) = X 0 .
where vec ( X ) denotes the vectorisation of the matrix X, formed by stacking the columns of X into a single column vector.
In the next section, we will solve a specific test problem numerically.

7. Numerical Results

In the test, we used the cubic B-spline B 3 as an approximation function.
We set
φ ˜ j k ( t ) : = 2 j 2 B 3 ( 2 j t k ) .
Therefore, we define the spaces:
V j [ 0 , 1 ] = span { φ ˜ j k ( t ) , k N j , t [ 0 , 1 ] } ¯ j j 0
where N j Z is the set of admissible indices k and j 0 is the initial multiresolution scale.
It is not difficult to prove that these spaces produce an MRA on [ 0 , 1 ] . See [25] for more details.
Now, we consider the following example of a fractional dynamical system in two dimensions:
D γ x ( t ) = 2 x ( t ) y ( t ) D γ y ( t ) = 4 x ( t ) 3 y ( t ) , t [ 0 , 1 ] x ( 0 ) = 1.2 y ( 0 ) = 4.2
where 0 < γ < 1 and whose exact solution is [31]
x ( t ) = 1 5 E γ ( t γ ) + E γ ( 2 t γ ) y ( t ) = 1 5 E γ ( t γ ) + 4 E γ ( 2 t γ ) .
E γ is the Mittag–Leffler function, i.e.,
E γ ( z ) = k = 0 z k Γ ( γ k + 1 ) .
We note that the matrix associated with the dynamical system has eigenvalues λ I such that min | arg λ i | > γ π / 2 , so that the stability of the system is guaranteed [32].
Using the collocation method with s = 8 and j = 8 , we numerically solved the dynamic system in the interval [ 0 , 1 ] , and we compared it with the exact solution (Figure 2). As expected, the error is of the order of the machine precision.

8. Conclusions and Future Work

We constructed a collocation method that uses fractional refinable functions as approximate functions, to solve a fractional non-stationary dynamical system. We provided an explicit formula that allowed us to evaluate the fractional derivatives of the approximate functions in an accurate and simple way. This formula uses a linear combination of fractional refinable functions of minor order to obtain the prescribed derivative. The method can be efficiently implemented by using standard multiscale techniques to evaluate the fractional B-splines and least mean squares for the final rectangular system. The numerical results showed that this collocation method approximates the solution of a test fractional system with good accuracy. In the future, we can use, as approximation functions, the fractional refinable functions described in [29]. The central point of the work is that this efficient technique can be extended to any fractional differential problem.

Author Contributions

Conceptualization, L.P. and S.D.L.; methodology, L.P. and S.D.L.; software, S.D.L.; writing—original draft preparation, L.P. and S.D.L.; writing—review and editing, L.P. and S.D.L.; funding acquisition, L.P. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was payed using the fund “00300_22_RS_PEZZA_SEED_PNR_2022-PEZZA-ATENEO SEED PNR 2022”.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Magin, R.L. Fractional calculus models of complex dynamics in biological tissues. Comput. Math. Appl. 2010, 59, 1586–1593. [Google Scholar] [CrossRef]
  2. Magin, R.L.; Ingo, C.; Colon-Perez, L.; Triplett, W.; Mareci, T.H. Characterization of anomalous diffusion in porous biological tissues using fractional order derivatives and entropy. Microporous Mesoporous Mater. 2013, 178, 39–43. [Google Scholar] [CrossRef] [PubMed]
  3. Hodaei, M.; Rabbani, V.; Maghoul, P. Transient acoustic wave propagation in bone-like porous materials using the theory of poroelasticity and fractional derivative: A sensitivity analysis. Acta Mech. 2020, 231, 179–203. [Google Scholar] [CrossRef]
  4. Fellah, M.; Fellah, Z.; Depollier, C. Transient wave propagation in inhomogeneous porous materials: Application of fractional derivatives. Signal Process. 2006, 86, 2658–2667. [Google Scholar] [CrossRef]
  5. Magin, R.L.; Abdullah, O.; Baleanu, D.; Zhou, X.J. Anomalous diffusion expressed through fractional order differential operators in the Bloch–Torrey equation. J. Magn. Reson. 2008, 190, 255–270. [Google Scholar] [CrossRef]
  6. Chen, W.; Sun, H.; Zhang, X.; Korošak, D. Anomalous diffusion modeling by fractal and fractional derivatives. Comput. Math. Appl. 2010, 59, 1754–1758. [Google Scholar] [CrossRef]
  7. Xiao-Jun, Y.; Tenreiro Machado, J.; Baleanu, D. Anomalous diffusion models with general fractional derivatives within the kernels of the extended Mittag-Leffler type functions. Rom. Rep. Phys. 2017, 69, 115. [Google Scholar]
  8. Bagley, R.L.; Torvik, P. A theoretical basis for the application of fractional calculus to viscoelasticity. J. Rheol. 1983, 27, 201–210. [Google Scholar] [CrossRef]
  9. Torvik, P.; Bagley, R. On the Appearance of the Fractional Derivative in the Behavior of Real Materials. J. Appl. Mech. 1984, 51, 294–298. [Google Scholar] [CrossRef]
  10. Heymans, N.; Bauwens, J.C. Fractal rheological models and fractional differential equations for viscoelastic behavior. Rheol. Acta 1994, 33, 210–219. [Google Scholar] [CrossRef]
  11. Chatterjee, A. Statistical origins of fractional derivatives in viscoelasticity. J. Sound Vib. 2005, 284, 1239–1245. [Google Scholar] [CrossRef]
  12. Pelap, F.; Tanekou, G.; Fogang, C.; Kengne, R. Fractional-order stability analysis of earthquake dynamics. J. Geophys. Eng. 2018, 15, 1673–1687. [Google Scholar] [CrossRef]
  13. Zeraick Monteiro, N.; Rodrigues Mazorche, S. Fractional derivatives applied to epidemiology. Trends Comput. Appl. Math. 2021, 22, 157–177. [Google Scholar] [CrossRef]
  14. Butt, A.I.K.; Imran, M.; Batool, S.; Nuwairan, M.A. Theoretical Analysis of a COVID-19 CF-Fractional Model to Optimally Control the Spread of Pandemic. Symmetry 2023, 15, 380. [Google Scholar] [CrossRef]
  15. Hanif, A.; Kashif Butt, A.I.; Ahmad, W. Numerical approach to solve Caputo-Fabrizio-fractional model of corona pandemic with optimal control design and analysis. Math. Methods Appl. Sci. 2023. [Google Scholar] [CrossRef]
  16. Ali, I.; Khan, S.U. A Dynamic Competition Analysis of Stochastic Fractional Differential Equation Arising in Finance via Pseudospectral Method. Mathematics 2023, 11, 1328. [Google Scholar] [CrossRef]
  17. Zeng, F.; Li, C.; Liu, F.; Turner, I. The use of finite difference/element approaches for solving the time-fractional subdiffusion equation. SIAM J. Sci. Comput. 2013, 35, A2976–A3000. [Google Scholar] [CrossRef]
  18. Pezza, L.; Pitolli, F. A fractional wavelet Galerkin method for the fractional diffusion problem. IMACS Ser. Comput. Appl. Math 2017, 20, 1–10. [Google Scholar]
  19. Li, X.; Xu, C. Existence and uniqueness of the weak solution of the space-time fractional diffusion equation and a spectral method approximation. Commun. Comput. Phys. 2010, 8, 1016. [Google Scholar]
  20. 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]
  21. Pitolli, F.; Pezza, L. A fractional spline collocation method for the fractional-order logistic equation. In Approximation Theory XV: San Antonio 2016; Springer: Cham, Switzerland, 2017; pp. 307–318. [Google Scholar]
  22. Pellegrino, E.; Pezza, L.; Pitolli, F. A collocation method in spline spaces for the solution of linear fractional dynamical systems. Math. Comput. Simul. 2020, 176, 266–278. [Google Scholar] [CrossRef]
  23. Pezza, L.; Tallini, L. Some Results on a New Refinable Class Suitable for Fractional Differential Problems. Fractal Fract. 2022, 6, 521. [Google Scholar] [CrossRef]
  24. Mazza, M.; Donatelli, M.; Manni, C.; Speleers, H. On the matrices in B-spline collocation methods for Riesz fractional equations and their spectral properties. Numer. Linear Algebra Appl. 2023, 30, e2462. [Google Scholar] [CrossRef]
  25. Pezza, L.; Pitolli, F. A multiscale collocation method for fractional differential problems. Math. Comput. Simul. 2017, 147, 210–219. [Google Scholar] [CrossRef]
  26. Pezza, L.; Pitolli, F. A fractional spline collocation-Galerkin method for the time-fractional diffusion equation. Commun. Appl. Ind. Math. 2018, 9, 104–120. [Google Scholar] [CrossRef]
  27. Unser, M.; Blu, T. Fractional splines and wavelets. SIAM Rev. 2000, 42, 43–67. [Google Scholar] [CrossRef]
  28. Podlubny, I. Fractional Differential Equations, Volume 198 of Mathematics in Science and Engineering; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  29. Gori, L.; Pitolli, F. A class of totally positive refinable functions. Rend. Di Mat. Ser. VII 2000, 20, 305–322. [Google Scholar]
  30. Mallat, S. A Wavelet Tour of Signal Processing; Academic Press, Inc.: San Diego, CA, USA, 1998; p. xxiv+577. [Google Scholar]
  31. Garrappa, R.; Popolizio, M. Evaluation of generalized Mittag–Leffler functions on the real line. Adv. Comput. Math. 2013, 39, 205–225. [Google Scholar] [CrossRef]
  32. Kaczorek, T.; Rogowski, K. Fractional Linear Systems and Electrical Circuits; Springer: Cham, Switzerland, 2015. [Google Scholar]
Figure 1. (a) The graphs of the fractional GP for α = 2:0.25:6, h = α + 3 , classical GP for α = 2 , , 6 , and classical B-splines for h = α integer. (b) The graphs of the fractional GP for α = 2:0.25:6, h = α + 0.5 , classical GP for α = 2 , , 6 , and classical B-splines for h = α integer.
Figure 1. (a) The graphs of the fractional GP for α = 2:0.25:6, h = α + 3 , classical GP for α = 2 , , 6 , and classical B-splines for h = α integer. (b) The graphs of the fractional GP for α = 2:0.25:6, h = α + 0.5 , classical GP for α = 2 , , 6 , and classical B-splines for h = α integer.
Axioms 12 00451 g001
Figure 2. The approximate solutions x j , y j with j = 8 (red line) obtained with the cubic B-spline φ 3 , 3 B 3 and the exact solutions x ( t ) , y ( t ) (blue dashed line). We consider four different example obtained whit different value for γ . In (a) we use γ = 0.10 , in (b) γ = 0.25 , in (c) γ = 0.5 and in (d) γ = 0.75 .
Figure 2. The approximate solutions x j , y j with j = 8 (red line) obtained with the cubic B-spline φ 3 , 3 B 3 and the exact solutions x ( t ) , y ( t ) (blue dashed line). We consider four different example obtained whit different value for γ . In (a) we use γ = 0.10 , in (b) γ = 0.25 , in (c) γ = 0.5 and in (d) γ = 0.75 .
Axioms 12 00451 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pezza, L.; Di Lillo, S. Fractional Dynamical Systems Solved by a Collocation Method Based on Refinable Spaces. Axioms 2023, 12, 451. https://doi.org/10.3390/axioms12050451

AMA Style

Pezza L, Di Lillo S. Fractional Dynamical Systems Solved by a Collocation Method Based on Refinable Spaces. Axioms. 2023; 12(5):451. https://doi.org/10.3390/axioms12050451

Chicago/Turabian Style

Pezza, Laura, and Simmaco Di Lillo. 2023. "Fractional Dynamical Systems Solved by a Collocation Method Based on Refinable Spaces" Axioms 12, no. 5: 451. https://doi.org/10.3390/axioms12050451

APA Style

Pezza, L., & Di Lillo, S. (2023). Fractional Dynamical Systems Solved by a Collocation Method Based on Refinable Spaces. Axioms, 12(5), 451. https://doi.org/10.3390/axioms12050451

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop