Next Article in Journal
Diffusion in Heterogenous Media and Sorption—Desorption Processes
Next Article in Special Issue
Special Issue: Fractal Functions and Applications
Previous Article in Journal
Optimizing the First-Passage Process on a Class of Fractal Scale-Free Trees
Previous Article in Special Issue
A Mathematical Study of a Coronavirus Model with the Caputo Fractional-Order Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multivariate Fractal Functions in Some Complete Function Spaces and Fractional Integral of Continuous Fractal Functions

by
Kshitij Kumar Pandey
and
Puthan Veedu Viswanathan
*
Department of Mathematics, IIT Delhi, New Delhi 110016, India
*
Author to whom correspondence should be addressed.
Fractal Fract. 2021, 5(4), 185; https://doi.org/10.3390/fractalfract5040185
Submission received: 24 July 2021 / Revised: 18 October 2021 / Accepted: 20 October 2021 / Published: 25 October 2021
(This article belongs to the Special Issue Fractal Functions and Applications)

Abstract

:
There has been a considerable evolution of the theory of fractal interpolation function (FIF) over the last three decades. Recently, we introduced a multivariate analogue of a special class of FIFs, which is referred to as α -fractal functions, from the viewpoint of approximation theory. In the current note, we continue our study on multivariate α -fractal functions, but in the context of a few complete function spaces. For a class of fractal functions defined on a hyperrectangle Ω in the Euclidean space R n , we derive conditions on the defining parameters so that the fractal functions are elements of some standard function spaces such as the Lebesgue spaces L p ( Ω ) , Sobolev spaces W m , p ( Ω ) , and Hölder spaces C m , σ ( Ω ) , which are Banach spaces. As a simple consequence, for some special choices of the parameters, we provide bounds for the Hausdorff dimension of the graph of the corresponding multivariate α -fractal function. We shall also hint at an associated notion of fractal operator that maps each multivariate function in one of these function spaces to its fractal counterpart. The latter part of this note establishes that the Riemann–Liouville fractional integral of a continuous multivariate α -fractal function is a fractal function of similar kind.

1. Preamble

This note aims to offer a modest contribution to the field of fractal interpolation. In particular, we consider a special class of fractal interpolation functions referred to as the α -fractal function, which has played a considerable role in the theory of univariate fractal approximation. Our work in the current note seeks to show that a few results on the construction of univariate α -fractal functions in various function spaces and associated fractal operator (see, for instance, [1]) carry over to higher dimensions.
For a prescribed data set D = { ( x 0 , y 0 ) , ( x 1 , y 1 ) , , ( x k , y k ) } in R 2 with increasing abscissae, there are multitude of methods to construct a continuous function that maps each x i to y i —generally known as interpolation methods—available in the field of classical numerical analysis and approximation theory. Roughly speaking, the fractal interpolation function (FIF for short), as introduced by Barnsley in the original version [2], is a continuous function g : [ x 0 , x k ] R that interpolates D such that the graph of g, denoted by Gr ( g ) , is a self-referential set (fractal set). Here the word fractal or self-referential is used to indicate that Gr ( g ) is the attractor of an iterated function system [3]. That is, roughly, Gr ( g ) is a finite union of tranformed copies of itself. For a compendium of the theory of FIF and its applications in interpolation and approximation, the reader is referred to the book and monograph [4,5]; the recent articles [6,7,8] may also be of interest.
In her research works on fractal interpolation, Navascués emphasized a special class of univariate FIFs, named α -fractal functions, (see, for instance, [9,10]) which garnered a significant amount of research attention in fractal approximation theory. It is our opinion that the notion of α -fractal functions assisted the field of fractal interpolation to find connections and consequences in other branches of mathematics such as approximation theory, harmonic analysis, functional analysis and the theory of bases and frames; see, for instance, [11,12]. In the research works reported in [13,14], authors utilized α -fractal functions to demonstrate that FIFs can be applied in various constrained approximation problems.
Several extensions of FIF to higher dimensions, in particular, bivariate FIFs or fractal surfaces, have been studied in the literature; see, for example, [4,15,16,17,18,19]. Despite that the α -fractal function facilitated the theory of univariate FIF to merge seamlessly with various fields in mathematics, a similar approach to multivariate FIFs was not attempted except for a few research works on bivariate α -fractal functions reported lately in [20,21,22]. The aforementioned works on bivariate α -fractal functions find their origin, perhaps implicitly, in the general framework for the construction of fractal surfaces introduced in [23].
While an increasing amount of literature is being published in the field of univariate FIFs and fractal surfaces, the research in multivariate FIFs are still inadequate, especially in the framework of α -fractal functions. In the context of multivariate FIFs, the ingenious constructions appeared in [24,25], though worth mentioning, do not seem to be suitable for the implementation of the α -fractal function formalism. On the other hand, our acquaintance with the univariate and bivariate α -fractal functions revealed that the development of multivariate analogue of α -fractal function could be highly beneficial for the expansion of multivariate fractal approximation theory. Stimulated by the construction of fractal surface in [23], recently we put forward a satisfactory extension of the Barnsley’s theory of univariate FIF to the multivariate case [26].
In this note, we continue to explore the notion of multivariate α -fractal functions. In the first part, we define multivariate α -fractal functions in various function spaces such as the Lebesgue spaces L p , Sobolev spaces W m , p , and Hölder spaces C m , σ . We also hint at some elementary properties of the fractal operator associated with the notion of multivariate α -fractal functions.
Fractal dimension is an important parameter of fractal geometry providing information about the geometric structure of the objects that it deals with. There are different notions of fractal dimension, the two most commonly used being the Hausdorff dimension and box dimension [27]. In particular, the Hausdorff dimension and box dimension of the graphs of fractal interpolation functions have been investigated; see, for instance, [2,4,15,28]. Since the aforementioned fractal dimensions are scale-independent, they may not be useful for describing scale-dependent laws and more complicated phenomena in nature. To this end, a new definition of fractal dimension, referred to as the two-scale dimension, is broached in [29], and it is perhaps more akin to physics than mathematics. However, we are forced to settle for less in the framework of multivariate α -fractal functions considered herein, because the analysis for the fractal dimension of the general nonaffine case is subtle. We shall just mention bounds for the Hausdorff dimension of the graph of the multivariate α -fractal function as an immediate consequence of its Hölder continuity for suitable choice of parameters.
On the other hand, fractional calculus, which broadly deals with derivatives and integrals of fractional order, is rather an old subject. During the last decades, fractional calculus has opened its wings wider to cover several real world applications in science and engineering. Despite being an old subject, fractional calculus continues to be a hot topic of research, resulting in a substantial body of literature; we refer the reader to the informative surveys [30,31]. Some recent developments made in the direction of fractional PDEs and their applications deserve a special mention; see, for instance, [32,33,34,35]. Studies on the interconnection between fractional calculus and fractal geometry have gained significant attention in recent years. For some links between the two-scale problem mentioned previously and fractional calculus, the reader may consult [36]. In the second part of this note, our modest aim is to show that the fractional integral of the multivariate fractal function considered herein is again a fractal function of a similar kind.
Overall, this note discusses how some results in univariate fractal interpolation, to be specific α -fractal functions, fractal operator and fractional calculus of fractal functions, carry over to higher dimensions. We strongly believe that these research findings may assist efforts to find interesting interconnections between multivariate FIFs and the theory of PDEs.

2. Preparatory Facts

To begin with, we list pertinent definitions and notation for use throughout the remainder of this note.
The set of first n natural numbers shall be denoted by Σ n . For l = ( l 1 , , l n ) ( N { 0 } ) n , called a multi-index, let | l | = k = 1 n l k . Given two multi-indices l = ( l 1 , , l n ) and l = ( l 1 , , l n ) , we say that l l if l k l k for all k Σ n . For l l , we define:
l l = k = 1 n l k l k .
Let n 2 and Ω = k = 1 n I k be an n-dimensional hyperrectangle, where each I k = [ a k , b k ] is a closed and bounded interval in R . For a function g : Ω R and X 0 = ( x 1 0 , , x n 0 ) Ω , denoted by
D l ( g ) ( X 0 ) = l 1 + + l n g x 1 l 1 x n l n ( X 0 ) ,
provided the right-hand side exists.

2.1. Function Spaces

The purpose here is to provide a short presentation of various function spaces that are used in this note. We refer to Triebel [37] for more information.
Let C ( Ω ) denote the Banach space of all real-valued continuous functions defined on Ω , endowed with the sup-norm . . For a positive integer m, we consider the linear space C m ( Ω ) defined by
C m ( Ω ) = { g C ( Ω ) : D l ( g ) exists and it is continuous for each multi - index l with | l | m } .
For any g C m ( Ω ) , we define
g C m ( Ω ) = | l | m D l ( g ) .
It is well-known that C m ( Ω ) equipped with . C m ( Ω ) is a Banach space. Next, we recall the Lebesgue spaces. For 0 < p , let
L p ( Ω ) = g : Ω R such that g is measurable and g p < ,
where g p is defined as
g p = Ω | g ( X ) | p d X 1 p , for 0 < p < . ess sup X Ω | g ( X ) | , for p = .
It is a standard result in functional analysis that L p ( Ω ) , . p is a Banach space for 1 p . For 0 < p < 1 , . p is a quasi-norm, that is, in place of the triangle inequality one has
g + h p 2 1 p g p + h p ,
and L p ( Ω ) is a quasi-Banach space.
Let g L p ( Ω ) . For a multi-index l , a function h : Ω R is called the l t h -weak derivative of g if it satisfies
Ω g ( X ) D l ϕ ( X ) d X = ( 1 ) | l | Ω h ( X ) ϕ ( X ) d X ,
for all infinitely differentiable functions ϕ with compact support contained in Ω . By a slight abuse of notation, we write l t h -weak derivative of g as D l ( g ) = h .
For 1 p and a non-negative integer m , W m , p ( Ω ) denotes the Sobolev space with smoothness m and integrability p defined by
W m , p ( Ω ) = g : Ω R such that D l ( g ) L p ( Ω ) for all multi - index l with | l | m .
The linear space W m , p ( Ω ) endowed with the norm
g W m , p ( Ω ) = | l | m D l ( g ) p p 1 p , for 1 p < . max | l | m D l ( g ) , for p = .
is a Banach space. For p = 2 , it is a Hilbert space, which shall be denoted by H m ( Ω ) .
A function g : Ω R is Hölder continuous with exponent σ ( 0 , 1 ] (or σ - Hölder continuous) if
| g ( X ) g ( Y ) | C g X Y σ ,
for all X , Y Ω and some C g > 0 , called a Hölder constant of g. Given a Hölder continuous g : Ω R with exponent σ , the σ -Hölder semi-norm of g is defined by
[ g ] σ = sup X , Y Ω , X Y | g ( X ) g ( Y ) | X Y σ .
If m is a positive integer, then the Hölder space C m , σ ( Ω ) is defined as
C m , σ ( Ω ) = { g C m ( Ω ) : D l ( g ) is σ - Hölder continuous for all multi - index l with | l | = m } .
The space C m , σ ( Ω ) equipped with the norm
g C m , σ ( Ω ) = g C m ( Ω ) + | l | = m D l ( g ) σ
is a Banach space. Note that C 0 , σ ( Ω ) coincides with the space of all Hölder continuous functions with exponent σ .

2.2. Towards Multivariate FIF

Here we shall equip ourselves with a few rudiments needed for multivariate fractal functions that concern us. As mentioned previously, let I k = [ a k , b k ] , k = 1 , 2 , , n be compact intervals in R and Ω = k = 1 n I k be an n-dimensional hyperrectangle.
Let n 2 be an integer and Δ = { ( x 1 , i 1 , x 2 , i 2 , , x n , i n ) R n : i k = 0 , 1 , , N k ; k = 1 , , n } be such that a k = x k , 0 < x k , 1 < < x k , N k = b k for each k = 1 , 2 , , n . Note that a k = x k , 0 < x k , 1 < < x k , N k = b k determines a partition of I k into subintervals I k , i k = x k , i k 1 , x k , i k for i k = 1 , 2 , , N k 1 and I k , N k = x k , N k 1 , x k , N k . It is worth to note that I k = i k = 1 N k I k , i k and each knot point in the partition of I k is exactly in one of the subintervals I k , i k , i k = 1 , 2 , N k mentioned above. We call such a set Δ as a partition of Ω for an obvious reason.
For convenience, let us introduce the following notation. For a positive integer N ,
Σ N = { 1 , 2 , , N } , Σ N , 0 = { 0 , 1 , N } , Σ N , 0 = { 0 , N } , int Σ N , 0 = { 1 , 2 , , N 1 } .
For each i k Σ N k , let u k , i k : I k I k , i k be an affine map of the form
u k , i k ( x ) = a k , i k x + b k , i k ,
satisfying
u k , i k ( x k , 0 ) = x k , i k 1 and u k , i k ( x k , N k ) = x k , i k if i k is odd , u k , i k ( x k , 0 ) = x k , i k and u k , i k ( x k , N k ) = x k , i k 1 if i k is even .
When the interval I k , i k involved in the definition of the affine map is half-open, the above equation needs to be interpreted in terms of the one-sided limit. For instance, when i k { 1 , 2 , , N k 1 } is odd, u k , i k ( x k , N k ) = x k , i k in (1) actually means lim x x k , N k u k , i k ( x ) = x k , i k .
Note that
| u k , i k ( x ) u k , i k ( x ) | γ k , i k | x x | , x , x I k ,
for 0 γ k , i k = | a k , i k | < 1 . Using the definition of the map u k , i k , one can verify that
u k , i k 1 ( x k , i k ) = u k , i k + 1 1 ( x k , i k ) ,
for all i k int Σ N k , 0 .
Let τ : Z × { 0 , N 1 , N 2 , , N n } Z be defined by
τ ( i , 0 ) = i 1 , τ ( i , N k ) = i , if i is odd , τ ( i , 0 ) = i , τ ( i , N k ) = i 1 , if i is even .
Using the above notation, we see that u k , i k ( x k , j k ) = x k , τ ( i k , j k ) for all i k Σ N k , j k Σ N k , 0 and k Σ n .
It is easy to observe that the boundary of Ω in the usual metric of R n is
Ω = X = ( x 1 , , x k , j k , x n ) Ω : j k Σ N k , 0 , k Σ n .

3. Multivariate α -Fractal Functions in Some Complete Function Spaces

This section targets to construct fractal functions (self-referential functions) in the complete function spaces C m ( Ω ) , L p ( Ω ) , W m , p ( Ω ) and C m , σ ( Ω ) , which we recalled in the previous section. To this end, let X be any of the function space from the list { C m ( Ω ) , L p ( Ω ) , W m , p ( Ω ) , C m , σ ( Ω ) } and f X , be a fixed function, which we shall refer to as the germ function. Let b X be a fixed function, called the base function.
For each g X , and X = ( x 1 , , x n ) k = 1 n I k , i k , ( i 1 , , i n ) k = 1 n Σ N k , we define T f ( g ) as
T f ( g ) ( X ) = f ( X ) + α i 1 i n ( g b ) u i 1 i n 1 ( X ) ,
where u i 1 i n 1 ( X ) = u 1 , i 1 1 ( x 1 ) , , u 1 , i 1 1 ( x n ) and α i 1 i n are real numbers such that
max { | α i 1 i n | : ( i 1 , , i n ) k = 1 n Σ N k } < 1 .
The k = 1 n N k -tuple comprised of the real numbers α i 1 i n is called the scaling vector and it is denoted by α . We define
α = max { | α i 1 i n | : ( i 1 , , i n ) k = 1 n Σ N k } .
The main objective in this section is to choose the scale vector α and base function b in (5) so that the Read-Bajraktarević (RB) operator T f is a well-defined map, and, in fact, T f is a contraction map on the function X or a suitable subspace of X . It is worth to emphasize that throughout the current note, a partition Δ of the hyperrectangle Ω is chosen as mentioned in the previous section.
Theorem 1.
Let f C m ( Ω ) and define
C f m ( Ω ) : = g C m ( Ω ) : g ( X ) = f ( X ) X Ω .
Suppose that the scaling vector α is so chosen that
max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | m < 1
and b C f m ( Ω ) .
Then the following hold.
  • The map T f given in (5) is well-defined on C f m ( Ω ) .
  • In fact, T f : C f m ( Ω ) C f m ( Ω ) C m ( Ω ) is a contraction map.
  • As a consequence, by the Banach fixed point theorem, there exists a unique function f Δ , b α C f m ( Ω ) such that
    D l ( f Δ , b α ) ( x 1 , i 1 , x 2 , i 2 , , x n , i n ) = D l ( f ) ( x 1 , i 1 , x 2 , i 2 , , x n , i n ) ,
    for all ( i 1 , , i n ) k = 1 n Σ N k , 0 and multi-index l with | l | m . Moreover, the function f Δ , b α and its derivatives satisfy the self-referential equations given by
    D l ( f Δ , b α ) ( X ) = D l ( f ) ( X ) + α i 1 i n k = 1 n a k , i k l k D l ( f Δ , b α b ) u i 1 i n 1 ( X ) ,
    for all X k = 1 n I k , i k , ( i 1 , , i n ) k = 1 n Σ N k and multi-index l with | l | m .
Proof. 
We shall first show that T f is well-defined on C f m ( Ω ) , that is, we show that T f ( g ) C f m ( Ω ) for all g C f m ( Ω ) .
Let g C f m ( Ω ) and X = ( x 1 , , x r , , x n ) k = 1 n I k , i k be such that x r I r , i r I r , i r + 1 for some r Σ n and i r int Σ N r , 0 .
Note that this is possible only when x r = x r , i r and in that case, by (3), we have
u 1 , i 1 1 ( x 1 ) , , u r , i r 1 ( x r , i r ) , , u n , i n 1 ( x n ) = u 1 , i 1 1 ( x 1 ) , , u r , i r + 1 1 ( x r , i r ) , , u n , i n 1 ( x n ) ,
and u i 1 i n 1 ( X ) Ω . So, by the specified choice of b, we have
D l ( g b ) u 1 , i 1 1 ( x 1 ) , , u r , i r 1 ( x r , i r ) , , u 1 , i 1 1 ( x n ) = 0 ,
for all multi-index l with | l | m . Thus,
D l ( T f ( g ) ) ( x 1 , , x r , , x n ) = D l ( f ) ( x 1 , , x r , , x n ) + α i 1 i n k = 1 n a k , i k l k [ D l ( g b ) ] u 1 , i 1 1 ( x 1 ) , , u r , i r 1 ( x r , i r ) , , u 1 , i 1 1 ( x n ) = D l ( f ) ( x 1 , , x r , , x n ) + α i 1 i n k = 1 n a k , i k l k [ D l ( g b ) ] u 1 , i 1 1 ( x 1 ) , , u r , i r 1 ( x r , i r + 1 ) , , u 1 , i 1 1 ( x n ) = D l ( f ) ( x 1 , , x r , , x n ) .
That is, D l T f ( g ) ( X ) = D l ( f ) ( X ) irrespective of whether x r is considered as a point in I r , i r or as a point in I r , i r + 1 . The above observation also yields the following:
  • D l T f ( g ) ( X ) = D l ( f ) ( X ) for all X Ω and | l | m .
  • D l T f ( g ) ( x 1 , i 1 , , x n , i n ) = D l ( f ) ( x 1 , i 1 , , x n , i n ) for all ( i 1 , , i n ) k = 1 n Σ N k and | l | m .
In particular, T f ( g ) C f m ( Ω ) .
Next, let g , h C f m ( Ω ) and l be a multi-index with | l | m . Then
| D l [ T f ( g ) T f ( h ) ] ( X ) | = | α i 1 i n k = 1 n a k , i k l k D l ( g h ) ( u i 1 i n 1 ( X ) ) | max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | m D l ( g h ) .
Taking sum over all | l | m , we get
T f ( g ) T f ( h ) C m ( Ω ) max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | m g h C m ( Ω ) .
Since, max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | m < 1 , the map T f : C f m ( Ω ) C f m ( Ω ) is a contraction. Rest of the claim follows by a simple application of the Banach fixed point theorem. □
Example 1.
Let us consider the surface in R 3 defined by the bivariate function f ( x , y ) = 1 x 2 + y 2 + 1 for all ( x , y ) [ 1 , 1 ] × [ 1 , 1 ] and a mesh partition Δ = { 1 , 0.5 , 0 , 0.5 , 1 } × { 1 , 0.5 , 0 , 0.5 , 1 } of the square [ 1 , 1 ] × [ 1 , 1 ] . Fractal functions f Δ , b α corresponding tofassociated with different choices of scale vector α and base function b are shown below.
Let us consider two base functions as follows:
b 1 ( x , y ) = 1 ( x 1 ) ( x + 1 ) ( y 1 ) ( y + 1 ) f ( x , y ) ,
and
b 2 ( x , y ) = 1 2 ( x + 1 ) f ( 1 , y ) ( x 1 ) f ( 1 , y ) + ( y + 1 ) f ( x , 1 ) ( y 1 ) f ( x , 1 ) 1 4 ( x + 1 ) ( y + 1 ) f ( 1 , 1 ) ( x 1 ) ( y + 1 ) f ( 1 , 1 ) 1 4 ( x 1 ) ( y 1 ) f ( 1 , 1 ) ( x + 1 ) ( y 1 ) f ( 1 , 1 ) .
Figure 1a is the graph of the germ function f ( x , y ) = 1 x 2 + y 2 + 1 . Figure 1b is the graph of fractal perturbation f Δ , b α with base function b = b 1 and uniform scale vector α , where α i 1 i 2 = 0.7 for all i k Σ 4 for k = 1 , 2 . Figure 1c depicts the graph of f Δ , b α with base function b = b 2 and uniform scale vector α as taken previously. Finally, Figure 1d displays the graph of f Δ , b α with base function b = b 2 and uniform scale vector α , where α i 1 i 2 = 0.01 . In this case, the parameters satisfy the conditions prescribed in Theorem 1, for m = 1 . Thus, Figure 1a,b corroborate the technique demonstrated for the construction of smoothness preserving fractal functions in Theorem 1.
Theorem 2.
Let f C m , σ ( Ω ) and define
C f m , σ ( Ω ) : = g C m , σ ( Ω ) : g ( X ) = f ( X ) X Ω .
Choose the scale vector satisfying
max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | m + σ < 1
and the base function b C f m , σ ( Ω ) . Then the RB operator T f defined in (5) is a contraction on C f m , σ ( Ω ) , and its unique fixed point f Δ , b α satisfies the self-referential Equation (6).
Proof. 
Using Theorem 1 we see that T f ( g ) C f m ( Ω ) for all g C f m , σ ( Ω ) . We shall show that for all multi-index l with | l | = m , D l ( T f ( g ) ) is Hölder continuous with exponent σ . Towards this, let X , Y k = 1 n I k , i k be two points in the same rectangular mesh. We have
| D l ( T f ( g ) ) ( X ) D l ( T f ( g ) ) ( Y ) | = | α i 1 i n k = 1 n a k , i k l k D l ( g b ) ( u i 1 i n 1 ( X ) ) D l ( g b ) ( u i 1 i n 1 ( Y ) ) | | α i 1 i n | k = 1 n | a k , i k | l k | D l ( g ) ( u i 1 i n 1 ( X ) ) D l ( g ) ( u i 1 i n 1 ( Y ) ) | + | α i 1 i n | k = 1 n | a k , i k | l k | D l ( b ) ( u i 1 i n 1 ( X ) ) D l ( b ) ( u i 1 i n 1 ( Y ) ) | C D l ( g ) | α i 1 i n | k = 1 n | a k , i k | l k u i 1 i n 1 ( X ) u i 1 i n 1 ( Y ) σ + C D l ( b ) | α i 1 i n | k = 1 n | a k , i k | l k u i 1 i n 1 ( X ) u i 1 i n 1 ( Y ) σ = C D l ( g ) | α i 1 i n | k = 1 n | a k , i k | l k + C D l ( b ) | α i 1 i n | k = 1 n | a k , i k | l k u i 1 i n 1 ( X ) u i 1 i n 1 ( Y ) σ max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | m + σ C D l ( g ) + C D l ( b ) X Y σ ,
where C D l ( g ) and C D l ( b ) denote the Hölder constants of D l ( g ) and D l ( b ) , respectively. If X and Y lie in two distinct but in adjacent meshes, then by taking point on their common boundary and repeating the above steps we get
| D l ( T f ( g ) ) ( X ) D l ( T f ( g ) ) ( Y ) | 2 max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | m + σ [ C D l ( g ) + C D l ( b ) ] X Y σ .
Since the total number of rectangular meshes is k = 1 n N k , for any X , Y Ω , we have
| D l ( T f ( g ) ) ( X ) D l ( T f ( g ) ) ( Y ) | ( k = 1 n N k ) max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | m + σ [ C D l ( g ) + C D l ( b ) ] X Y σ ,
which shows that T f C f m , σ ( Ω ) C f m , σ ( Ω ) .
A similar computation reveals also that the map T f : C f m , σ ( Ω ) C f m , σ ( Ω ) is a contraction map, completing the proof. □
Corollary 1.
Let f : Ω R be a Hölder continuous function with exponent σ . Assume that a scaling vector α is so chosen that
max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | k = 1 n | a k , i k | σ < 1
and the parameter mapbis a Hölder continuous function with exponent σ and b ( X ) = f ( X ) for all X Ω . Then the Hausdorff dimension of the graph of the corresponding self-referential function f Δ , b α satisfies
n dim H Gr ( f Δ , b α ) n + 1 σ .
Proof. 
With the stated hypotheses on α and b, it follows from the previous theorem (with m = 0 ) that the self-referential counterpart f Δ , b α of f is a Hölder continuous function with exponent σ . Define a map A : Gr ( f Δ , b α ) Ω by
A x , f Δ , b α ( x ) = x ,
where we endow Gr ( f Δ , b α ) R n + 1 and Ω R n with the usual Euclidean norm. It is plain to see that A is a surjective Lipschitz map. From fundamental properties of the Hausdorff dimension given in ([38], Theorem 2, Items (5), (8)) we have
n = dim H ( Ω ) = dim H A Gr ( f Δ , b α ) dim H Gr ( f Δ , b α ) .
For the desired upper bound, let us recall that the Hausdorff dimension of the graph of a Hölder continuous function with Hölder exponent s ( 0 , 1 ] whose domain is a compact subset of R n with the Hausdorff dimension equal to d is less than or equal to min { d + 1 s , d s } ([39], Chapter 10). Therefore,
dim H Gr ( f Δ , b α ) n + 1 σ ,
completing the proof. □
Theorem 3.
Let f L p ( Ω ) for 0 < p . Suppose the scaling vector α is so chosen that
i n = 1 N n i 1 = 1 N 1 ( k = 1 n | a k , i k | ) | α i 1 i n | p 1 / p < 1 , for 1 p < . α < 1 , for p = . i n = 1 N n i 1 = 1 N 1 ( k = 1 n | a k , i k | ) | α i 1 i n | p < 1 , for 0 < p < 1 .
Then T f defined in (5) maps L p ( Ω ) to L p ( Ω ) . Further, T f is a contraction map and hence by the Banach fixed point theorem, there exists a unique f Δ , b α L p ( Ω ) such that
f Δ , b α ( X ) = f ( X ) + α i 1 i n ( f Δ , b α b ) u i 1 i n 1 ( X ) ,
for X k = 1 n I k , i k , and ( i 1 , , i n ) k = 1 n Σ N k .
Proof. 
Using the stated hypotheses, it is easy to verify that the operator T f : L p ( Ω ) L p ( Ω ) is well-defined. What remains is to show that T f is a contraction map. To this end, let g , h L p ( Ω ) , 1 p < . We have
T f ( g ) T f ( h ) p p = Ω | ( T f ( g ) T f ( h ) ) ( X ) | p d X = i n = 1 N n i 1 = 1 N 1 X k = 1 n I k , i k | α i 1 i n ( g h ) u i 1 i n 1 ( X ) | p d X = i n = 1 N n i 1 = 1 N 1 ( k = 1 n | a k , i k | ) | α i 1 i n | p Ω | ( g h ) ( X ˜ ) | p d X ˜ = i n = 1 N n i 1 = 1 N 1 ( k = 1 n | a k , i k | ) | α i 1 i n | p g h p p
Thus,
T f ( g ) T f ( h ) p = i n = 1 N n i 1 = 1 N 1 ( k = 1 n | a k , i k | ) | α i 1 i n | p 1 / p g h p ,
proving the claim for the case 1 p < . The other cases can be dealt similarly. □
Next, let us construct self-referential functions associated with a function f W m , p ( Ω ) . First, let us recall the following result, popularly known as the Leibniz theorem.
If f W m , p ( Ω ) and ϕ is infinitely differentiable on Ω , then ϕ f W m , p ( Ω ) and
D l ( ϕ f ) = q = 0 l l p D p ( ϕ ) D l q ( f ) , | l | m .
Theorem 4.
Let f W m , p ( Ω ) for 1 p . Suppose that the base function b W m , p ( Ω ) and the scaling vector is chosen so that
i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | m p 1 1 / p < 1 , for 1 p < . max { | α i 1 i n | k = 1 n | a k , i k | m : ( i 1 , , i n ) k = 1 n Σ N k } < 1 , for p = .
Then the RB operator T f given in (5) is a contraction map on W m , p ( Ω ) . Consequently, T f has a unique fixed point f Δ , b α .
Proof. 
A routine computation yields that the RB operator is well-defined and it maps the space W m , p ( Ω ) into itself. We shall just show that it is a contraction on W m , p ( Ω ) . To this end, let 1 p < and l be a multi-index with | l | m . We note that
D l T f ( g ) T f ( h ) p p = Ω | D l [ T f ( g ) T f ( h ) ] ( X ) | p d X = Ω | D l [ T f ( g ) T f ( h ) ] ( X ) | p d X = i n = 1 N n i 1 = 1 N 1 k = 1 n I k , i k | D l [ T f ( g ) T f ( h ) ] ( X ) | p d X = i n = 1 N n i 1 = 1 N 1 k = 1 n I k , i k | D l α i 1 i n ( g h ) ( u i 1 i n 1 ( X ) ) | p d X = i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n I k , i k | 1 k = 1 n a k , i k l k D l ( g h ) ( u i 1 i n 1 ( X ) ) | p d X = i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | p l k 1 Ω | D l ( g h ) ( X ) | p d X = i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | p l k 1 D l ( g h ) p p .
Thus, for a multi-index l with | l | m we have
D l [ T f ( g ) T f ( h ) ] p = i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | p l k 1 1 p D l ( g h ) p .
Hence,
T f ( g ) T f ( h ) W m , p ( Ω ) = | l | m i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | p l k 1 1 p D l ( g h ) p i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | p m 1 1 p | l | m D l ( g h ) p i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | p m 1 1 p g h W m , p ( Ω )
The rest of the theorem follows from the Banach fixed point theorem and the assumption on the scale vector. The case p = can be worked out similarly. □

4. Fractal Operator on Function Spaces

Let f X , where X is a fixed function space from the list
C m ( Ω ) , L p ( Ω ) , W m , p ( Ω ) , C m , σ ( Ω ) .
The results established in the previous section provide a self-referential counterpart to each f X , and consequently provide an operator. That is, for a prescribed set of parameters such as the partition, scale vector and the base function, there exists a fractal operator F Δ , b α : X X defined by F Δ , b α ( f ) = f Δ , b α . This section intends to record a few elementary properties of the operator F Δ , b α : X X , what we call a multivariate self-referential operator (fractal operator); see also [1]. We shall provide the details only for X = W m , p ( Ω ) , as the other spaces can be similarly dealt with. For future reference, we introduce the notation
K α , m , p = i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | m p 1 1 / p , for 1 p < . max | α i 1 i n | k = 1 n | a k , i k | m : i k Σ N k , k Σ n , for p = .
Proposition 1.
(Perturbation Error) Let f W m , p ( Ω ) . Suppose that a partition Δ of the hyperrectangle Ω , base function b W m , p ( Ω ) , and scale vector α be chosen as in Theorem 4. Then
f Δ , b α f W m , p ( Ω ) K α , m , p f Δ , b α b W m , p ( Ω ) .
Proof. 
Let us recall the self-referential equations satisfied by the fractal counterpart f Δ , b α W m , p ( Ω ) and its derivatives
D l ( f Δ , b α ) ( X ) = D l ( f ) ( X ) + α i 1 i n k = 1 n a k , i k l k D l ( f Δ , b α b ) ( u i 1 i n 1 ( X ) ) ,
for all X k = 1 n I k , i k , ( i 1 , , i n ) k = 1 n Σ N k and multi-index l with | l | m . Assume that 1 p < . By simple calculations
D l f Δ , b α f p = i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | p l k 1 1 p D l ( f Δ , b α b ) p .
Therefore,
f Δ , b α f W m , p ( Ω ) = | l | m D l ( f Δ , b α f ) p p 1 p i n = 1 N n i 1 = 1 N 1 | α i 1 i n | p k = 1 n | a k , i k | m p 1 1 / p f Δ , b α b W m , p ( Ω ) .
Similar analysis for p = . □
Now, let us take the multivariate base function b : Ω R used in the construction of the self-referential function f Δ , b α through a suitable operator L : X X . That is, we take b = L ( f ) so that the conditions required for b are satisfied. In this case, the multivariate fractal operator will be denoted by F Δ , L α : X X . In what follows, we intend to record some elementary properties of the multivariate fractal operator F Δ , L α .
The following proposition provides a counterpart to the linearity property of the fractal operator well explored in the setting of univariate α -fractal functions on various function spaces; see, for instance, [11]. The proof follows almost verbatim, and hence omitted.
Proposition 2.
Let X = W m . p ( Ω ) and L : X X be a linear operator. Choose the base functionbin the construction of fractal function f Δ , b α via this operatorLso that b = L ( f ) . Then the corresponding fractal operator, which shall be denoted by F Δ , L α : X X , defined by F Δ , L α ( f ) = f Δ , L ( f ) α is linear.
Let X be a Banach space and A : X X be a bounded linear operator such that I A < 1 , where I is the identity operator on X. Then, it is well-known that A is bijective and A 1 is bounded; see, for instance, [40]. The following result available in [41] is a generalization of the aforementioned Neumann’s lemma.
Lemma 1.
([41], Lemma 1) Let A : X X be a linear operator on a Banach space X such that
A ( x ) x λ 1 x + λ 2 A ( x ) x X ,
for some λ 1 and λ 2 [ 0 , 1 ) . Then A is a topological automorphism (a bounded, invertible map that possesses a bounded inverse). Furthermore,
1 λ 1 1 + λ 2 x A ( x ) 1 + λ 1 1 λ 2 x , 1 λ 2 1 + λ 1 x A 1 ( x ) 1 + λ 2 1 λ 1 x , x X .
Proposition 3.
Let L : W m , p ( Ω ) W m , p ( Ω ) be a bounded linear operator and the scale vector α be chosen such that max K α , m , p , L K α , m , p < 1 . Then the linear operator F Δ , L α : W m , p ( Ω ) W m , p ( Ω ) is a topological automorphism.
Proof. 
Recall that here the base function b = L ( f ) so that by Proposition 1 we have
F Δ , L α ( f ) f W m , p ( Ω ) = f Δ , L ( f ) α f W m , p ( Ω ) K α , m , p f Δ , L ( f ) α L ( f ) W m , p ( Ω ) K α , m , p f Δ , L ( f ) α W m , p ( Ω ) + K α , m , p L f W m , p ( Ω ) = K α , m , p F Δ , L α ( f ) W m , p ( Ω ) + K α , m , p L f W m , p ( Ω ) .
The assertion is now immediate from the previous proposition. □
The existence of Schauder bases consisting of appropriate functions for the Sobolev spaces is quite desirable in analysis of PDEs, for instance, for demonstrating the existence of solutions of various non-linear boundary value problems. We have the following result giving a Schauder basis consisting of self-referential functions for the Sobolev space W m , p ( Ω ) . The heart of the matter is an elementary result in the theory of bases, which states that a topological isomorphism preserves Schauder bases; see, for instance, [37].
Corollary 2.
The Banach space W m , p ( Ω ) has a Schauder basis consisting of multivariate self-referential functions.
Proof. 
Let { f m } m N be a Schauder basis of W m , p ( Ω ) whose existence is established and reported, for instance, in [42,43,44]. Choose the scale function α and operator L as in the previous proposition so that the fractal operator F Δ , L α is a topological automorphism. As an isomorphism, in particular, an automorphism, preserves Schauder bases, we conclude that { ( f m ) Δ , L α } m N , where ( f m ) Δ , L α = F Δ , L α ( f m ) is a Schauder basis consisting of self-referential functions for the Banach space W m , p ( Ω ) . □

5. Fractional Integral of Continuous Multivariate α -Fractal Function

As mentioned in the introductory section, exploration of interconnection between fractional calculus and fractal geometry has always been of interest. Our purpose in this section is limited; we shall observe that the Riemann–Liouville fractional integral of the continuous multivariate α -fractal function is also a fractal function. A similar result regarding univariate FIF can be found in [28].
Definition 1.
[45] Let f be a continuous function on the closed and bounded hyperrectangle Ω in R n . The left-hand-sided mixed Riemann–Liouville fractional integral of f of order γ is defined as
I a γ f ( X ) = 1 k = 1 n Γ ( γ k ) a 1 x 1 a n x n ( x 1 s 1 ) γ 1 1 ( x n s n ) γ n 1 . f ( s 1 , s 2 , , s n ) d s 1 d s 2 d s n ,
where a = ( a 1 , , a n ) is a fixed point, X = ( x 1 , x 2 , x n ) and γ = ( γ 1 , , γ n ) with a k x k , γ k > 0 for each k Σ n .
Let f C ( Ω ) . We write f ( X ) = f ( x 1 , x 2 , , x n ) . From Theorem 1 it follows that by choosing b C f ( Ω ) and scaling vector α such that
α : = max ( i 1 , , i n ) k = 1 n Σ N k | α i 1 i n | < 1 ,
the fractal counterpart f Δ , b α of f belongs to C ( Ω ) . Furthermore, since f Δ , b α is the fixed point of the RB operator T f : C ( Ω ) C ( Ω ) defined by
T f ( g ) ( X ) = f ( X ) + α i 1 i n ( g b ) ( u i 1 i n 1 ( X ) ) ,
for all X k = 1 n I k , i k , ( i 1 , , i n ) k = 1 n Σ N k . Consequently, f Δ , b α satisfies the functional equation
f Δ , b α ( u i 1 i n ( X ) ) = f ( u i 1 i n ( X ) ) + α i 1 i n ( f Δ , b α b ) ( X ) = α i 1 i n f Δ , b α ( X ) + f ( u i 1 i n ( X ) ) α i 1 i n b ( X ) .
Let us define
q i 1 i n ( X ) = f u i 1 i n ( X ) α i 1 i n b ( X )
so that the self-referential equation for f Δ , b α becomes
f Δ , b α u i 1 i n ( X ) = α i 1 i n f Δ , b α ( X ) + q i 1 i n ( X ) .
Since the multivariate fractal function f Δ , b α is continuous, we can talk about its Riemann–Liouville fractional integral. In what follows, we establish that the Riemann–Liouville fractional integral of f Δ , b α is again a fractal function.
For the sake of convenience, we shall deal with the uniform scaling factor, that is, α i 1 i n = α for all ( i 1 , , i n ) k = 1 n Σ N k . Then, with a slight abuse of notation, the above equation reduces to
f Δ , b α u i 1 i n ( X ) = α f Δ , b α ( X ) + q i 1 i n ( X ) .
Theorem 5.
Let Δ be a partition of the hyperrectangle Ω in R n and f C ( Ω ) . Assume that b : Ω R is continuous and b ( X ) = f ( X ) for all X Ω , the boundary of Ω. Choose a scaling vector α such that α < 1 . Then I a γ f Δ , b α , the left-hand-sided mixed Riemann–Liouville fractional integral of order γ of the self-referential function f Δ , b α , satisfies the following equation:
I a γ f Δ , b α u i 1 i n ( X ) = α k = 1 n a k , i k γ k I a γ f Δ , b α ( X ) + q ^ i 1 i n ( X ) ,
where
q ^ i 1 , i n ( X ) = 1 k = 1 n Γ ( γ k ) k = 0 n 1 j = 0 k a n j , i n j γ n j a 1 u 1 , i 1 ( x 1 ) a n k u n k , i n k ( x n k ) a n k + 1 x n k + 1 a n x n u 1 , i 1 ( x 1 ) s 1 γ 1 1 u n k , i n k ( x n k ) s n k γ n k 1 ( x n k + 1 s n k + 1 ) γ n k + 1 1 ( x n s n ) γ n 1 f Δ , b α s 1 , , s n k , u n k + 1 , i n k + 1 ( s n k + 1 ) , , u n , i n ( s n ) d s 1 d s n + k = 1 n a k , i k γ k I a γ q i 1 i n ( X ) .
Proof. 
According to Theorem 1, it follows that f Δ , b α is continuous on Ω and satisfies the equation
f Δ , b α u i 1 i n ( X ) = α f Δ , b α ( X ) + q i 1 i n ( X ) , X Ω .
Hence,
I a γ f Δ , b α u i 1 i n ( X ) = 1 k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) a n u n , i n ( x n ) u 1 , i 1 ( x 1 ) s 1 γ 1 1 u n , i n ( x n ) s n γ n 1 f Δ , b α ( s 1 , , s n ) d s 1 d s n = 1 k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) a n u n , i n ( a n ) ( u 1 , i 1 ( x 1 ) s 1 ) γ 1 1 ( u n , i n ( x n ) s n ) γ n 1 f Δ , b α ( s 1 , , s n ) d s 1 d s n + 1 k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) u n , i n ( a n ) u n , i n ( x n ) u 1 , i 1 ( x 1 ) s 1 γ 1 1 u n , i n ( x n ) s n γ n 1 f Δ , b α ( s 1 , , s n ) d s 1 d s n .
Let us write
E 0 = 1 k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) a n u n , i n ( a n ) u 1 , i 1 ( x 1 ) s 1 γ 1 1 u n , i n ( x n ) s n γ n 1 f Δ , b α ( s 1 , , s n ) d s 1 d s n , J 1 = 1 k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) u n , i n ( a n ) u n , i n ( x n ) u 1 , i 1 ( x 1 ) s 1 γ 1 1 u n , i n ( x n ) s n γ n 1 f Δ , b α ( s 1 , , s n ) d s 1 d s n .
so that
I a γ f Δ , b α u i 1 i n ( X ) = E 0 + J 1 .
Turning our attention to J 1 , let us change the variable s n using the transformation s n = u n , i n ( t n ) . We have
J 1 = a n , i n γ n k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) a n x n u 1 , i 1 ( x 1 ) s 1 γ 1 1 ( x n t n ) γ n 1 f Δ , b α s 1 , , u n , i n ( t n ) d s 1 d t n = a n , i n γ n k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) a n 1 u n 1 , i n 1 ( x n 1 ) a n x n u 1 , i 1 ( x 1 ) s 1 γ 1 1 ( x n s n ) γ n 1 f Δ , b α s 1 , , u n , i n ( s n ) d s 1 d s n .
Applying similar process to the variable s n 1 we get
J 1 = a n , i n γ n k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) a n 1 u n 1 , i n 1 ( a n 1 ) a n x n ( u 1 , i 1 ( x 1 ) s 1 ) γ 1 1 ( x n s n ) γ n 1 f Δ , b α s 1 , , u n , i n ( s n ) d s 1 d s n + a n , i n γ n k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) u n 1 , i n 1 ( a n 1 ) u n 1 , i n 1 ( x n 1 ) a n x n ( u 1 , i 1 ( x 1 ) s 1 ) γ 1 1 ( x n s n ) γ n 1 f Δ , b α s 1 , , u n , i n ( s n ) d s 1 d s n = : E 1 + J 2 .
In J 2 , let us perform a change of variable using s n 1 = u n 1 , i n 1 ( t n 1 ) so that
J 2 = a n 1 , i n 1 γ n 1 a n , i n γ n k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) a n 1 x n 1 a n x n ( u 1 , i 1 ( x 1 ) s 1 ) γ 1 1 ( x n 1 s n 1 ) γ n 1 1 ( x n s n ) γ n 1 f Δ , b α s 1 , , u n 1 , i n 1 ( s n 1 ) , u n , i n ( s n ) d s 1 d s n .
Consequently,
I a γ f Δ , b α u i 1 i n ( X ) = E 0 + E 1 + J 2 .
Proceeding in the same fashion, at the n t h step we get
I a γ f Δ , b α u i 1 i n ( X ) = k = 0 n 1 E k + k = 1 n a k , i k γ k k = 1 n Γ ( γ k ) a 1 x 1 a n 1 x n 1 a n x n ( x 1 s 1 ) γ 1 1 ( x n 1 s n 1 ) γ n 1 1 ( x n s n ) γ n 1 f Δ , b α u i 1 i n ( S ) d s 1 d s n ,
where
E k = j = 0 k a n j , i n j γ n j k = 1 n Γ ( γ k ) a 1 u 1 , i 1 ( x 1 ) a n k u n k , i n k ( x n k ) a n k + 1 x n k + 1 a n x n u 1 , i 1 ( x 1 ) s 1 γ 1 1 u n k , i n k ( x n k ) s n k γ n k 1 ( x n k + 1 s n k + 1 ) γ n k + 1 1 ( x n s n ) γ n 1 f Δ , b α s 1 , , s n k , u n k + 1 , i n k + 1 ( s n k + 1 ) , , u n , i n ( s n ) d s 1 d s n
for k = 0 , 1 , , n 1 , with the assumption that j = 0 0 a n j , i n j γ n j = 1 .
Finally, using the functional equation
f Δ , b α u i 1 i n ( S ) = α f Δ , b α ( S ) + q i 1 i n ( S ) ,
for all S Ω , we get
I a γ f Δ , b α u i 1 i n ( X ) = k = 0 n 1 E k + k = 1 n a k , i k γ k k = 1 n Γ ( γ k ) a 1 x 1 a n 1 x n 1 a n x n ( x 1 s 1 ) γ 1 1 ( x n 1 s n 1 ) γ n 1 1 ( x n s n ) γ n 1 α f Δ , b α ( S ) + q i 1 i n ( S ) d s 1 d s n = α k = 1 n a k , i k γ k k = 1 n Γ ( γ k ) a 1 x 1 a n 1 x n 1 a n x n ( x 1 s 1 ) γ 1 1 ( x n 1 s n 1 ) γ n 1 1 ( x n s n ) γ n 1 f Δ , b α ( S ) d s 1 d s n + ( k = 0 n 1 E k + k = 1 n a k , i k γ k k = 1 n Γ ( γ k ) a 1 x 1 a n 1 x n 1 a n x n ( x 1 s 1 ) γ 1 1 ( x n 1 s n 1 ) γ n 1 1 ( x n s n ) γ n 1 q i 1 i n ( S ) d s 1 d s n ) = α k = 1 n a k , i k γ k I a γ f Δ , b α ( X ) + q ^ i 1 i n ( X ) ,
as desired. □

6. Conclusions

The α -fractal formalism of fractal interpolation function is proved to be beneficial in expanding the applications of univariate fractal approximation theory. Through the construction of multivariate α -fractal functions on a few complete function spaces which are ubiquitous in the theory of partial differential equations and harmonic analysis, the present work intends to be a step forward in the theory of multivariate fractal approximation. The construction of self-referential analogue for each germ function in a complete function space under consideration leads naturally to an operator, referred to as the multivariate fractal operator. We have studied a few elementary properties of the fractal operator. The multivariate fractal operator introduced and studied here enabled us, in particular, to construct Schauder bases consisting of self-referential functions for the function spaces. Further, taking a slight detour from the main theme, it is shown that the Riemann–Liouville fractional integral of a self-referential counterpart of the given multivariate germ function will also be a self-referential function under some suitable conditions.

Author Contributions

Conceptualization, P.V.V.; Methodology, K.K.P. and P.V.V.; Supervision, P.V.V.; Writing—original draft, K.K.P.; Writing—review & editing, P.V.V. All authors have read and agreed to the published version of the manuscript.

Funding

The second author is thankful to the project CRG/2020/002309 from the Science and Engineering Research Board (SERB), Government of India.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the referees for their divers remarks and suggestions. In particular, we thank one of the referees for bringing the two-scale fractal (dimension) theory and references [32,33,34,35] to our attention.

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. Viswanathan, P.; Navascués, M.A. A Fractal Operator on Some Standard Spaces of Functions. Proc. Edinb. Math. Soc. 2017, 60, 771–786. [Google Scholar] [CrossRef] [Green Version]
  2. Barnsley, M.F. Fractal functions and interpolation. Constr. Approx. 1986, 2, 303–329. [Google Scholar] [CrossRef]
  3. Hutchinson, J. Fractals and self-similarity. Indiana Univ. Math. J. 1981, 30, 713–747. [Google Scholar] [CrossRef]
  4. Massopust, P.R. Fractal Functions, Fractal Surfaces and Wavelets; Academic Press: Cambridge, MA, USA, 2014. [Google Scholar]
  5. Massopust, P.R. Interpolation and Approximation with Splines and Fractals; Oxford University Press: New York, NY, USA, 2010. [Google Scholar]
  6. Kok, C.W.; Tam, W.S. Fractal Image Interpolation: A Tutorial and New Result. Fractal Fract. 2019, 3, 7. [Google Scholar] [CrossRef] [Green Version]
  7. Navascués, M.A.; Chand, A.K.B.; Veedu, V.P.; Sebastián, M.V. Fractal Interpolation Functions: A Short Survey. Appl. Math. 2014, 5, 1834–1841. [Google Scholar] [CrossRef] [Green Version]
  8. Ri, S., II; Drakopoulos, V.; Nam, S.-M. Fractal Interpolation Using Harmonic Functions on the Koch Curve. Fractal Fract. 2021, 5, 28. [Google Scholar] [CrossRef]
  9. Navascués, M.A. Fractal polynomial interpolation. Z. Anal. Anwend. 2005, 25, 401–418. [Google Scholar] [CrossRef]
  10. Navascués, M.A. Fractal trigonometric approximation. Electron. Trans. Numer. Anal. 2005, 20, 64–74. [Google Scholar]
  11. Navascués, M.A. Fractal approximation. Complex Anal. Oper. Theory 2010, 4, 953–974. [Google Scholar] [CrossRef]
  12. Navascués, M.A. Fractal bases for Lp-spaces. Fractals 2012, 20, 141–148. [Google Scholar] [CrossRef]
  13. Viswanathan, P.; Chand, A.K.B.; Navascués, M. Fractal perturbation preserving fundamental shapes: Bounds on the scale factors. J. Math. Anal. Appl. 2014, 419, 804–817. [Google Scholar] [CrossRef]
  14. Viswanathan, P.; Navascués, M.A.; Chand, A.K.B. Fractal Polynomials and Maps in Approximation of Continuous Functions. Numer. Funct. Anal. Optim. 2015, 37, 106–127. [Google Scholar] [CrossRef]
  15. Bouboulis, P.; Dalla, L.; Drakopoulos, V. Construction of recurrent bivariate fractal interpolation surfaces and computation of their box-counting dimension. J. Approx. Theory 2006, 141, 99–117. [Google Scholar] [CrossRef] [Green Version]
  16. Chand, A.K.B.; Kapoor, G.P. Hidden Variable Bivariate Fractal Interpolation Surfaces. Fractals 2003, 11, 277–288. [Google Scholar] [CrossRef]
  17. Dalla, L. Bivariate fractal interpolation functions on grids. Fractals 2002, 10, 53–58. [Google Scholar] [CrossRef]
  18. Feng, Z. Variation and Minkowski dimension of fractal interpolation surfaces. J. Math. Anal. Appl. 1993, 176, 561–586. [Google Scholar] [CrossRef] [Green Version]
  19. Massopust, P.R. Fractal surfaces. J. Math. Anal. Appl. 1990, 151, 275–290. [Google Scholar] [CrossRef] [Green Version]
  20. Navascues, M.A.; Mohapatra, R.N.; Akhtar, M.N. Construction of Fractal Surfaces. Fractals 2020, 28, 13. [Google Scholar] [CrossRef]
  21. Verma, S.; Viswanathan, P. A Fractal Operator Associated with Bivariate Fractal Interpolation Functions on Rectangular Grids. Results Math. 2020, 75, 28. [Google Scholar] [CrossRef]
  22. Verma, S.; Viswanathan, P. Parameter Identification for a Class of Bivariate Fractal Interpolation Functions and Constrained Approximation. Numer. Funct. Anal. Optim. 2020, 41, 1109–1148. [Google Scholar] [CrossRef]
  23. Ruan, H.-J.; Xu, Q. Fractal interpolation surfaces on rectangular grids. Bull. Aust. Math. Soc. 2015, 91, 435–446. [Google Scholar] [CrossRef]
  24. Bouboulis, P.; Dalla, L. A general construction of fractal interpolation functions on grids of R n. Eur. J. Appl. Math. 2007, 18, 449–476. [Google Scholar] [CrossRef] [Green Version]
  25. Hardin, D.; Massopust, P.R. Fractal interpolation function from R n to R m and their projections. Z. Anal. Anw. 1993, 12, 561–586. [Google Scholar]
  26. Pandey, K.K.; Viswanathan, P. Multivariate fractal interpolation functions: Some approximation aspects and an associated fractal interpolation operator. arXiv 2021, arXiv:2104.02950V1. [Google Scholar]
  27. Falconer, K. Fractal Geometry—Mathematical Foundations and Applications, 3rd ed.; John Wiley: Hoboken, NJ, USA, 2014. [Google Scholar]
  28. Ruan, H.-J.; Su, W.-Y.; Yao, K. Box dimension and fractional integral of linear fractal interpolation functions. J. Approx. Theory 2009, 161, 187–197. [Google Scholar] [CrossRef] [Green Version]
  29. Ain, Q.T.; He, J.-H. On two-scale dimension and its applications. Therm. Sci. 2019, 23, 1707–1712. [Google Scholar] [CrossRef] [Green Version]
  30. Machado, J.A.T.; Kiryakova, V.; Mainardi, F. Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1140–1153. [Google Scholar] [CrossRef] [Green Version]
  31. Sun, H.G.; Zhang, Y.; Baleanu, D.; Chen, W.; Chen, Y. A new collection of rel world applications of fractional calculus in science and engineering. Commun. Nonlinear Sci. Numer. Simul. 2018, 64, 213–231. [Google Scholar] [CrossRef]
  32. Dai, D.-D.; Ban, T.-T.; Wang, Y.-L.; Zhang, W. The piecewise reproducing kernel method for the time variable fractional order advection-reaction-diffusion equations. Therm. Sci. 2021, 25, 1261–1268. [Google Scholar] [CrossRef]
  33. Orovio, A.B.; Kay, D.; Burrage, K. Fourier spectral methods for fractional-in-space reaction-diffusion equations. BIT Numer. Math. 2014, 54, 937–954. [Google Scholar] [CrossRef]
  34. Tian, Y.; Liu, J. A modified exp-function method for fractional partial differential equations. Therm. Sci. 2021, 25, 1237–1241. [Google Scholar] [CrossRef]
  35. Tian, Y.; Liu, J. Direct algebraic method for solving fractional Fokas equation. Therm. Sci. 2021, 25, 2235–2244. [Google Scholar] [CrossRef]
  36. He, J.-H.; Ji, F.-Y. Two-scale mathematics and fractional calculus for thermodynamics. Therm. Sci. 2019, 23, 2131–2133. [Google Scholar] [CrossRef]
  37. Triebel, H. Theory of Function Spaces; Birkhaüser: Basel, Switzerland, 1992. [Google Scholar]
  38. Schleicher, D. Hausdorff dimension, its properties and its surprises. Amer. Math. Monthly 2007, 114, 509–528. [Google Scholar] [CrossRef] [Green Version]
  39. Kahane, J.P. Some Random Series of Functions, 2nd ed.; Cambridge Univ. Press: Cambridge, UK, 1985. [Google Scholar]
  40. Bachman, G.; Narici, L. Functional Analysis; Dover: New York, NY, USA, 1998. [Google Scholar]
  41. Casazza, P.G.; Christensen, O. Perturbation of operators and application to frame theory. J. Fourier Anal. Appl. 1997, 3, 543–557. [Google Scholar] [CrossRef]
  42. Fučík, S.; John, O.; Nečas, J. On the existence of Schauder bases in Sobolev spaces. Comment. Math. Univ. Carolinae 1972, 13, 163–175. [Google Scholar]
  43. Garrigós, G.; Seeger, A.; Ullrich, T. The Haar system as a Schauder basis in the spaces of Hardy-Sobolev type. J. Fourier Anal. Appl. 2018, 24, 1319–1339. [Google Scholar] [CrossRef] [Green Version]
  44. Singer, I. Bases in Banach Spaces I; Springer: Berlin, Germany, 1970. [Google Scholar]
  45. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives, Theory and Applications; Gordon and Breach: Yverdon, Switzerland, 1993. [Google Scholar]
Figure 1. Fractal functions corresponding to the seed function f associated with different choices of base function and scale vector.
Figure 1. Fractal functions corresponding to the seed function f associated with different choices of base function and scale vector.
Fractalfract 05 00185 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pandey, K.K.; Viswanathan, P.V. Multivariate Fractal Functions in Some Complete Function Spaces and Fractional Integral of Continuous Fractal Functions. Fractal Fract. 2021, 5, 185. https://doi.org/10.3390/fractalfract5040185

AMA Style

Pandey KK, Viswanathan PV. Multivariate Fractal Functions in Some Complete Function Spaces and Fractional Integral of Continuous Fractal Functions. Fractal and Fractional. 2021; 5(4):185. https://doi.org/10.3390/fractalfract5040185

Chicago/Turabian Style

Pandey, Kshitij Kumar, and Puthan Veedu Viswanathan. 2021. "Multivariate Fractal Functions in Some Complete Function Spaces and Fractional Integral of Continuous Fractal Functions" Fractal and Fractional 5, no. 4: 185. https://doi.org/10.3390/fractalfract5040185

Article Metrics

Back to TopTop