Non-Stationary Fractal Interpolation

Abstract

We introduce the novel concept of a non-stationary iterated function system by considering a countable sequence of distinct set-valued maps ${\left\{{\mathcal{F}}_k\right\}}_{k\in \mathbb{N}}$ where each ${\mathcal{F}}_k$ maps $\mathcal{H}\left(X\right)\to \mathcal{H}\left(X\right)$ and arises from an iterated function system. Employing the recently-developed theory of non-stationary versions of fixed points and the concept of forward and backward trajectories, we present new classes of fractal functions exhibiting different local and global behavior and extend fractal interpolation to this new, more flexible setting.

1. Introduction

Contractive operators on complete function spaces play a crucial role in the theory of differential and integral equations and are essential for the development of iterative solvers. One class of such contractive operators is defined on the graphs of functions using a particular type of iterated function system (IFS). The fixed points of such IFSs are graphs of functions that exhibit fractal characteristics. There exists an extensive literature on IFSs and fractal functions including, for instance, [1,2,3].

Up to now, the construction of contractive operators on sets or functions uses primarily sequences of iterates of one operator. Recently, motivated by non-stationary subdivision algorithms, a more general class of sequences consisting of different contractive operators was introduced in [4] and their limit properties studied. These ideas were then extended in [5] to sequences of different contractive operators mapping between different spaces. Using different contractive operators provide one with the ability to construct limit attractors that have different shapes or features at different scales.

This article uses the aforementioned new ideas to introduce the novel concept of non-stationary IFS and non-stationary fractal interpolation. These new ideas widen the applicability of fractal functions and fractal interpolation as they now include scale- and location-dependent features.

The outline of this paper is as follows. After providing some necessary preliminaries in Section 2, some results from [4] are presented in Section 3. In Section 4, (stationary) fractal interpolation and the associated (stationary) IFSs are reviewed. Non-stationary fractal functions are constructed in Section 5, and non-stationary fractal interpolation is introduced in Section 6. The final Section 7 defines non-stationary fractal functions on the Bochner–Lebesgue ${\mathcal{L}}^p$-spaces with $0<p\le \infty$.

2. Preliminaries

Let $(X,d)$ be a complete metric space. Given a map $f:X\to X$, the Lipschitz constant associated with f is defined by:

$$\mathrm{Lip}\left(f\right)=\underset{x,y\in X,x\ne y}{sup}\frac{d\left(f\left(x\right),f\left(y\right)\right)}{d(x,y)}.$$

The map f is called Lipschitz if $\mathrm{Lip}\left(f\right)<+\infty$ and a contraction on X if $\mathrm{Lip}\left(f\right)<1$.

Definition 1.

Let $(X,d)$ be a complete metric space and $\mathcal{F}:=\{f_1,\dots ,f_n\}$ a finite family of contractions on X. Then, the pair $(X,\mathcal{F})$ is called a contractive iterated function system (IFS) on X.

Remark 1.

As we deal exclusively with contractive IFSs in this article, we drop the adjective “contractive” in the following.

Remark 2.

In order to avoid trivialities, we henceforth assume that the number of maps in an IFS is an integer greater than one.

With an IFS $(X,\mathcal{F})$ and its point maps $f\in \mathcal{F}$, we can associate a set-valued mapping, also denoted by $\mathcal{F}$, as follows. Let $\left(\mathcal{H}\right(X),h)$ be the hyperspace of all nonempty compact subsets of X endowed with the Hausdorff–Pompeiu metric:

$$h(S_1,S_2):=\max \{d(S_1,S_2),d(S_2,S_1)\},$$

where $d(S_1,S_2):=\underset{x\in S_1}{sup}d(x,S_2):=\underset{x\in S_1}{sup}\underset{y\in S_2}{inf}d(x,y)$.

Define the mapping $\mathcal{F}:\mathcal{H}\left(X\right)\to \mathcal{H}\left(X\right)$ by [1,6]:

$$\mathcal{F}\left(S\right):=\bigcup _{i=1}^n{f}_i\left(S\right).$$

(1)

It is known that for contractive mappings $f\in \mathcal{F}$, the set-valued map $\mathcal{F}$ defined by (1) is a contractive Lipschitz map on $\mathcal{H}\left(X\right)$ with Lipschitz constant $\mathrm{Lip}\left(\mathcal{F}\right)=\max \{\mathrm{Lip}\left(f_i\right):i\in {\mathbb{N}}_n\}$. Here, we set ${\mathbb{N}}_n:=\{1,\dots ,n\}$. Moreover, the completeness of $(X,d)$ implies the completeness of $\left(\mathcal{H}\right(X),h)$.

The next definition is motivated by the validity of the Banach fixed point theorem in the above setting.

Definition 2.

The unique fixed point $A\in \mathcal{H}\left(X\right)$ of the contractive set-valued map $\mathcal{F}$ is called the attractor of the IFS $(X,\mathcal{F})$.

Note that since A satisfies the self-referential equation:

$$A=\mathcal{F}\left(A\right)=\bigcup _{i=1}^n{f}_i\left(A\right),$$

(2)

the attractor is in general a fractal set.

It follows directly from the proof of the Banach fixed point theorem that the attractor A is obtained as the limit (in the Hausdorff–Pompeiu metric) of the iterative process $A_k:=\mathcal{F}\left(A_{k-1}\right)$, $k\in \mathbb{N}$:

$$A=\underset{k\to \infty }{lim}{A}_k=\underset{k\to \infty }{lim}{\mathcal{F}}^k\left(A_0\right),$$

(3)

for an arbitrary $A_0\in \mathcal{H}\left(X\right)$. Here, ${\mathcal{F}}^k$ denotes the k-fold composition of $\mathcal{F}$ with itself.

We refer to the element $A_k\in \mathcal{H}\left(X\right)$ as the $k^{\mathrm{th}}$ level approximant of A or as a pre-fractal of rank k [2].

3. Systems of Function Systems

In [4], a generalization of IFSs was presented. The idea for this generalization comes from the theory of subdivision schemes. Instead of using only one set-valued map $\mathcal{F}$ to obtain an iterative process ${\left\{A_n\right\}}_{n\in \mathbb{N}}$ with initial $A_0\in \mathcal{H}\left(X\right)$, a sequence of function systems consisting of different families $\mathcal{F}$ is considered.

To this end, let $(X,d)$ be a complete metric space and ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ a sequence of transformations $T_k:X\to X$.

Definition 3.

([4], Definition 3.6) Let ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ be a sequence of transformations $T_k:X\to X$. A subset $\mathcal{I}$ of X is called an invariant set of the sequence ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ if:

$$\forall k\in \mathbb{N}\forall x\in \mathcal{I}:T_k\left(x\right)\in \mathcal{I}.$$

A criterion for obtaining an invariant domain for a sequence ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ of transformations on X is given below.

Proposition 1.

([4], Lemma 3.7) Let ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ be a sequence of transformations on $(X,d)$. Suppose there exists a $q\in X$ such that for all $x\in X$:

$$d(T_k\left(x\right),q)\le \mu d(x,q)+M,$$

for some $\mu \in [0,1)$ and $M>0$. Then, the ball $B_r\left(q\right)$ of radius $r=M/(1-\mu )$ centered at q is an invariant set for ${\left\{T_k\right\}}_{k\in \mathbb{N}}$.

Proof. 

For the proof, we refer the interested reader to [4]. □

Now, suppose that ${\left\{{\mathcal{F}}_k\right\}}_{k\in \mathbb{N}}$ is a sequence of set-valued maps ${\mathcal{F}}_k:\mathcal{H}\left(X\right)\to \mathcal{H}\left(X\right)$ defined by:

$${\mathcal{F}}_k\left(A_0\right):=\bigcup _{i=1}^{n_k}{f}_{i,k}\left(A_0\right),A_0\in \mathcal{H}\left(X\right),$$

(4)

where ${\mathcal{F}}_k=\{f_{i,k}:i\in {\mathbb{N}}_{n_k}\}$ is a family of contractions constituting an IFS on a complete metric space $(X,d)$. Setting $s_{i,k}:=\mathrm{Lip}\left(f_{i,k}\right)$, we obtain that $\mathrm{Lip}\left({\mathcal{F}}_k\right)=\max \{s_{i,k}:i\in {\mathbb{N}}_{n_k}\}<1$.

The following definitions are taken from [4] (Section 4).

Definition 4.

Let $A_0\in \mathcal{H}\left(X\right)$. The sequences:

$${\mathsf{\Phi }}_k\left(A_0\right):={\mathcal{F}}_k\circ {\mathcal{F}}_{k-1}\circ \cdots \circ {\mathcal{F}}_1\left(A_0\right)$$

(5)

and:

$${\mathsf{\Psi }}_k\left(A_0\right):={\mathcal{F}}_1\circ {\mathcal{F}}_2\circ \cdots \circ {\mathcal{F}}_k\left(A_0\right)$$

(6)

are called the forward and backward trajectories of $A_0$, respectively.

For our current setting, it was shown in ([4], Corollary 4.2) that if:

(i)

$n:=n_k$, for all $k\in \mathbb{N}$;

(ii)

there exists a common nonempty compact invariant set $\mathcal{I}\subseteq X$ for the maps $\left\{f_{i,k}\right\}$, $i\in {\mathbb{N}}_n$, $k\in \mathbb{N}$, such that ${\left\{f_{i,k}\right\}}_{k\in \mathbb{N}}$ converges uniformly on $\mathcal{I}$ to $f_i$ as $k\to \infty$;

(iii)

the IFS $(X,\mathcal{F})$ with $\mathcal{F}=\{f_i:i\in {\mathbb{N}}_n\}$ is contractive on $(X,d)$;

then the forward trajectory $\left\{{\mathsf{\Phi }}_k\left(A_0\right)\right\}$ converges for an arbitrary $A_0\subseteq \mathcal{I}$ to the unique attractor of $(X,\mathcal{F})$.

It was observed in [4] that the limits of forward trajectories do not lead to new classes of fractals. On the other hand, backward trajectories converge under rather mild conditions, even when forward trajectories do not converge to a (contractive) IFS, and generate new types of fractal sets.

As the convergence of backward trajectories is important for this article, we summarize the result in the next theorem, whose proof the reader can find in [4].

Theorem 1.

([4], Corollary 4.4) Let ${\left\{{\mathcal{F}}_k\right\}}_{k\in \mathbb{N}}$ be a family of set-valued maps of the form (4) whose elements are collections ${\mathcal{F}}_k=\{f_{i,k}:i\in {\mathbb{N}}_{n_k}\}$ of contractions constituting IFSs on a complete metric space $(X,d)$. Suppose that:

(i)

there exists a nonempty closed invariant set $\mathcal{I}\subseteq X$ for $\left\{f_{i,k}\right\}$, $i\in {\mathbb{N}}_{n_k}$, $k\in \mathbb{N}$;

(ii)

and:

$$\sum _{k=1}^{\infty }\prod _{j=1}^k\mathrm{Lip}\left({\mathcal{F}}_j\right)<\infty .$$

(7)

Then, the backward trajectories $\left\{{\mathsf{\Psi }}_k\left(A_0\right)\right\}$ converge for any initial $A_0\subseteq \mathcal{I}$ to a unique attractor $A\subseteq \mathcal{I}$.

Remark 3.

In [4] (Proposition 3.11), it was required that the invariant set $\mathcal{I}$ be compact. However, it suffices to only require that $\mathcal{I}$ be closed as $(X,d)$ is complete (see the proof of Proposition 3.11 in [4]).

Remark 4.

The conditions for convergence of the forward and backward trajectories are more general in [4]. For our purposes and setting, the above criteria are however sufficient.

Remark 5.

Fractals generated by backwards trajectories allow for more flexibility in their shapes. By a proper choice of IFSs, one can construct fractals exhibiting different local behavior (cf. [4]). This is due to the fact that in the sequence:

$${\mathcal{F}}_1\circ {\mathcal{F}}_2\circ \cdots {\mathcal{F}}_{k-1}\circ {\mathcal{F}}_k\left(A_0\right),A_0\in H\left(X\right),$$

the global shape of the attractor is determined by the initial maps ${\mathcal{F}}_1\circ {\mathcal{F}}_2\dots$, whereas the local shape is given by the final maps ${\mathcal{F}}_{k-1}\circ {\mathcal{F}}_k\dots$. Thus, scaling the attractor by $\mathrm{Lip}\left({\mathsf{\Psi }}_k\right)$, ${\mathsf{\Psi }}_k={\mathcal{F}}_1\circ {\mathcal{F}}_2\circ \cdots {\mathcal{F}}_{k-1}\circ {\mathcal{F}}_k$, reveals the behavior of the attractor of ${\left\{{\mathcal{F}}_m\right\}}_{m>k}$. See also [4] (Example 5.1).

Remark 6.

A comparison to V-variable fractals [7] was also undertaken in [4] (Section 4.1), showing that systems of function systems (SFSs) have weaker prerequisites than V-variable fractals.

4. Fractal Interpolation

Before introducing the new concept of non-stationary fractal interpolation, we need to briefly recall the rudimentaries of (stationary) fractal interpolation and (stationary) fractal functions. This is the purpose of the current section.

4.1. Stationary Fractal Interpolation

Suppose we are given a finite family ${\left\{l_i\right\}}_{i=1}^n$ of injective contractions $X\to X$ generating a partition of X in the sense that:

$$\begin{array}{c}X=\bigcup _{i=1}^n{l}_i\left(X\right);\hfill \end{array}$$

(8)

$$\begin{array}{c}{l}_i\left(X\right)\cap l_j\left(X\right)=\mathsf{\varnothing },\forall i,j\in {\mathbb{N}}_n,i\ne j.\hfill \end{array}$$

(9)

Let $(Y,d_Y)$ be a complete metric space with metric $d_Y$. A map $g:X\to Y$ is called bounded (with respect to the metric $d_Y$) if there exists a constant $M>0$ such that for all $x_1,x_2\in X$, $d_Y(g\left(x_1\right),g\left(x_2\right))<M$.

Recall that the set $\mathcal{B}(X,Y):=\{g:X\to Y:g\mathrm{is}\mathrm{bounded}\}$ becomes a complete metric space when endowed with the metric:

$${\displaystyle d(g,h):=\underset{x\in X}{sup}{d}_Y(g\left(x\right),h\left(x\right)).}$$

(10)

Remark 7.

Under the usual addition and scalar multiplication of functions, the space $\mathcal{B}(X,Y)$ is actually a metric linear space, i.e., a vector space for which the operations of vector addition and scalar multiplication are continuous (see, for instance, [8]).

For $i\in {\mathbb{N}}_n$, let $F_i:X\times Y\to Y$ be a mapping that is uniformly contractive in the second variable, i.e., there exists a $c\in [0,1)$ so that for all $y_1,y_2\in Y$:

$$d_Y(F_i(x,y_1),F_i(x,y_2))\le c{d}_Y(y_1,y_2),\forall x\in X,\forall i\in {\mathbb{N}}_n.$$

(11)

Define an operator $T:\mathcal{B}(X,Y)\to \mathcal{B}(X,Y)$ by:

$$Tg\left(x\right):=\sum _{i=1}^n{F}_i(l_i^{-1}\left(x\right),g\circ l_i^{-1}\left(x\right)){\chi }_{l_i\left(X\right)}\left(x\right),$$

(12)

where ${\chi }_M$ denotes the characteristic function of a set M. Such operators are referred to as Read–Bajractarević (RB) operators. The operator T is well defined, and since g is bounded and each $F_i$ contractive in the second variable, $Tg\in \mathcal{B}(X,Y)$.

Equivalently, (12) can also be expressed in the form:

$$(Tg\circ l_i)\left(x\right):=F_i(x,g\left(x\right)),x\in X,i\in {\mathbb{N}}_n.$$

(13)

Moreover, (11) implies that T is contractive on $\mathcal{B}(X,Y)$:

$$\begin{array}{cc}\hfill d(Tg,Th)& =\underset{x\in X}{sup}{d}_Y(Tg\left(x\right),Th\left(x\right))\hfill \\ & =\underset{x\in X}{sup}{d}_Y(F(l_i^{-1}\left(x\right),g\left(l_i^{-1}\left(x\right)\right)),F(l_i^{-1}\left(x\right),h\left(l_i^{-1}\left(x\right)\right)))\hfill \\ & \le c\underset{x\in X}{sup}{d}_Y(g\circ l_i^{-1}\left(x\right),h\circ l_i^{-1}\left(x\right))\le c{d}_Y(g,h).\hfill \end{array}$$

(14)

To achieve notational simplicity, we set $F(x,y):={\displaystyle \sum _{i=1}^n}{F}_i(x,y){\chi }_X\left(x\right)$ in the above equation.

Therefore, by the Banach fixed point theorem, T has a unique fixed point $f^{*}$ in $\mathcal{B}(X,Y)$. This unique fixed point is called the bounded fractal function (generated by T), and it satisfies the self-referential equation:

$$f^{*}\left(x\right)=\sum _{i=1}^n{F}_i(l_i^{-1}\left(x\right),f^{*}\circ l_i^{-1}\left(x\right)){\chi }_{l_i\left(X\right)}\left(x\right),$$

(15)

or, equivalently,

$$f^{*}\circ l_i\left(x\right)=F_i(x,f^{*}\left(x\right)),x\in X,i\in {\mathbb{N}}_n.$$

(16)

The fixed point $f^{*}\in \mathcal{B}(X,Y)$ is obtained as the limit of the sequence of mappings:

$$T^k\left(f_0\right)\to f^{*},\mathrm{as}k\to \infty ,$$

(17)

where $f_0\in \mathcal{B}(X,Y)$ is arbitrary.

Next, we would like to consider a special choice for the mappings $F_i$. To this end, we require the concept of an F-space. We recall that a metric $d:Y\times Y\to \mathbb{R}$ is called complete if every Cauchy sequence in Y converges (with respect to d) to a point of Y, and translation-invariant if:

$$d(x+a,y+a)=d(x,y),\mathrm{for}\mathrm{all}x,y,a\in Y.$$

Now, assume that Y is an F-space, i.e., a topological vector space whose topology is induced by a complete translation-invariant metric d and, in addition, that this metric is homogeneous. This setting allows us to consider mappings $F_i$ of the form:

$$F_i(x,y):=q_i\left(x\right)+S_i\left(x\right)y,i\in {\mathbb{N}}_n,$$

(18)

where $q_i\in \mathcal{B}(X,Y)$ and $S_i:X\to \mathbb{R}$ is a function.

The metric $d_Y$ is homogeneous, and therefore, the mappings (18) satisfy condition (11) if the functions $S_i$ are assumed to be bounded on X with bounds in $[0,1)$. For then:

$$\begin{array}{cc}\hfill d_Y(q_i\left(x\right)+S_i\left(x\right)y_1,q_i\left(x\right)+S_i\left(x\right)y_2)& =d_Y(S_i\left(x\right)y_1,S_i\left(x\right)y_2)\hfill \\ & =|S_i\left(x\right)|d_Y(y_1,y_2)\le {\parallel S_i\parallel }_{\infty }{d}_Y(y_1,y_2)\le s{d}_Y(y_1,y_2).\hfill \end{array}$$

Here, ${\parallel ·\parallel }_{\infty }$ denotes the supremum norm and $s:=\max \{\parallel S_i{\parallel }_{\infty }:i\in {\mathbb{N}}_n\}$. Henceforth, we will assume that all functions $S_i$ are bounded above by $s\in [0,1)$.

With the choice (18), the RB operator T becomes an affine operator on $\mathcal{B}(X,Y)$ of the form:

$$\begin{array}{cc}\hfill Tg& =\sum _{i=1}^n{q}_i\circ l_i^{-1}{\chi }_{l_i\left(X\right)}+\sum _{i=1}^n{S}_i\circ l_i^{-1}·g\circ l_i^{-1}{\chi }_{l_i\left(X\right)}\hfill \end{array}$$

(19)

$$\begin{array}{c}=T\left(0\right)+\sum _{i=1}^n{S}_i\circ l_i^{-1}·g\circ l_i^{-1}{\chi }_{l_i\left(X\right)}.\hfill \end{array}$$

(20)

Next, we exhibit the relationship between the graph $G\left(f^{*}\right)$ of the fixed point $f^{*}$ of the operator T given by (12) and the attractor of an associated contractive IFS.

To this end, consider the complete metric space $X\times Y$, and define mappings $w_i:X\times Y\to X\times Y$ by:

$$w_i(x,y):=(l_i\left(x\right),F_i(x,y)),i\in {\mathbb{N}}_n.$$

(21)

Assume that the mappings $F_i$ in addition to being uniformly contractive in the second variable are also uniformly Lipschitz continuous in the first variable, i.e., there exists a constant $L>0$ such that for all $y\in Y$,

$$d_Y(F_i(x_1,y),F_i(x_2,y))\le L{d}_X(x_1,x_2),\forall x_1,x_2\in X,\forall i=1,\dots ,n.$$

Denote by $a:=\max \{a_i:i\in {\mathbb{N}}_n\}$ the largest of the contractivity constants of the $l_i$, and let $\theta :=\frac{1-a}{2L}$. Then, the mapping $d_{\theta }:(X\times Y)\times (X\times Y)\to \mathbb{R}$ given by:

$$d_{\theta }:=d_X+\theta d_Y$$

is a metric on $X\times Y$ compatible with the product topology on $X\times Y$.

The next theorem is a special case of a result presented in [9].

Theorem 2.

The family ${\mathcal{F}}_T:=(X\times Y,w_1,w_2,\dots ,w_n)$ is a contractive IFS in the metric $d_{\theta }$, and the graph $G\left(f^{*}\right)$ of the fractal function $f^{*}$ generated by the RB operator T given by (12) is the unique attractor of ${\mathcal{F}}_T$. Moreover,

$$G\left(Tg\right)={\mathcal{F}}_T\left(G\left(g\right)\right),\forall g\in B(X,Y),$$

(22)

where ${\mathcal{F}}_T$ denotes the set-valued operator (1).

Equation (22) can be represented by the following commutative diagram:

$$\begin{array}{ccc}X\times Y& \stackrel{{\mathcal{F}}_T}{\to }& X\times X\\ \phantom{G}\uparrow G& & \phantom{G}\uparrow G\\ \mathcal{B}(X,Y)& \stackrel{T}{\to }& \mathcal{B}(X,Y)\end{array}$$

(23)

where G is the mapping $\mathcal{B}(X,Y)\ni g\mapsto G\left(g\right)=\left\{\right(x,g\left(x\right)):x\in X\}\in X\times Y$.

Conversely, assume that $\mathcal{F}=(X\times Y,w_1,w_2,\dots ,w_n)$ is an IFS whose maps $w_i$ are of the form (21) where the functions $l_i$ are contractive injections satisfying (8) and (9), and the mappings $F_i$ are uniformly Lipschitz continuous in the first variable and uniformly contractive in the second variable. Thus, an RB operator $T_{\mathcal{F}}$ of the form (12) can be associated with the IFS $\mathcal{F}$. The attractor $A_{\mathcal{F}}$ of $\mathcal{F}$ is then the graph $G\left(f\right)$ of the fixed point f of $T_{\mathcal{F}}$ (this was the original approach in [10] to define a fractal interpolation function on a compact interval in $\mathbb{R}$). The diagram (23) is then commutative with ${\mathcal{F}}_T$ replaced by $\mathcal{F}$ and T replaced by $T_{\mathcal{F}}$.

We now specialize even further and choose arbitrary $f,b\in \mathcal{B}(X,Y)$ and set:

$$q_i:=f\circ l_i-S_i·b.$$

(24)

Then, the RB operator T becomes:

$$Tg=f+(S_i\circ l_i^{-1})·(g-b)\circ l_i^{-1},\mathrm{on}{l}_i\left(X\right),i\in {\mathbb{N}}_n,$$

(25)

and under the assumption that $s<1$, its unique fixed point $f^{*}\in \mathcal{B}(X,Y)$ satisfies the self-referential equation:

$$f^{*}=f+(S_i\circ l_i^{-1})·(f^{*}-b)\circ l_i^{-1},\mathrm{on}{l}_i\left(X\right),i\in {\mathbb{N}}_n.$$

(26)

Remark 8.

The functions f and b are referred to as the seed and base function, respectively.

Remark 9.

The fixed point $f^{*}$ in (26) clearly depends on the seed function f, the base function b, and the scaling functions $S_i$. Fixing f and b, but varying $S_i$ generate an uncountable family of fractal functions $f^{*}=f^{*}(S_1,\dots ,S_n)$ originating from $f=f^{*}(0,\dots ,0)$.

In the case of univariate fractal interpolation on the real line with $X:=[a,b]$, $-\infty <a<b<+\infty$, the base function b can be chosen to be the affine function whose graph connects the points $(a,f(a\left)\right)$ and $(b,f(b\left)\right)$.

If we consider the complete metric space of continuous functions $\left(\mathcal{C}\right(X,\mathbb{R}),d)$ instead of $\left(\mathcal{B}\right(X,\mathbb{R}),d)$, define:

$$x_0:=a,x_n:=b,\mathrm{and}{x}_i:=l_i\left(b\right),i\in {\mathbb{N}}_n,$$

and impose the join-up conditions:

$$Tf(x_j-)=Tf(x_j+),j\in {\mathbb{N}}_{n-1},$$

(27)

the fixed point $f^{*}$ will be a continuous function whose graph interpolates the set $\{(x_j,f\left(x_j\right)):j=0,1,\dots ,n\}$. Such functions are usually referred to as fractal interpolation functions [6,10]. As the RB operator is the same at each level of recursion (17), we refer to this as stationary fractal interpolation.

5. Non-Stationary Fractal Functions

Here, we introduce non-stationary versions of the concepts of fractal functions as presented in the previous section.

To this end, consider a doubly-indexed family of injective contractions $\{l_{i_k,k}:i_k\in {\mathbb{N}}_{n_k},k\in \mathbb{N}\}$ from $X\to X$ generating a partition of X for each $k\in \mathbb{N}$ in the sense of (8) and (9).

Suppose that Y is an F-space, $\{q_{i_k,k}:i_k\in {\mathbb{N}}_{n_k},k\in \mathbb{N}\}\subset \mathcal{B}(X,Y)$ and $\{S_{i_k,k}:i_k\in {\mathbb{N}}_{n_k},k\in \mathbb{N}\}\subset \mathcal{B}(X,\mathbb{R})$ is such that:

$$s:=\underset{k\in \mathbb{N}}{sup}\underset{i_k\in {\mathbb{N}}_k}{max}{\parallel S_{i_k,k}\parallel }_{\infty }<1.$$

For each $k\in \mathbb{N}$, define an RB operator $T_k:\mathcal{B}(X,Y)\to \mathcal{B}(X,Y)$ by:

$$\begin{array}{cc}\hfill T_{k}f& :=\sum _{i_k=1}^{n_k}{q}_{i_k,k}\circ l_{i_k,k}^{-1}{\chi }_{l_{i_k,k}\left(X\right)}+\sum _{i_k=1}^{n_k}{S}_{i_k,k}\circ l_{i_k,k}^{-1}·f\circ l_{i_k,k}^{-1}{\chi }_{l_{i_k,k}\left(X\right)}\hfill \end{array}$$

(28)

$$\begin{array}{c}=T_k\left(0\right)+\sum _{i_k=1}^{n_k}{S}_{i_k,k}\circ l_{i_k,k}^{-1}·f\circ l_{i_k,k}^{-1}{\chi }_{l_{i_k,k}\left(X\right)}.\hfill \end{array}$$

(29)

It is straightforward to verify that each RB operator $T_k$ is a contraction on $\mathcal{B}(X,Y)$ with Lipschitz constant:

$$\mathrm{Lip}\left(T_k\right)=\underset{i_k\in {\mathbb{N}}_k}{max}{\parallel S_{i_k,k}\parallel }_{\infty }\le s<1.$$

(30)

Proposition 2.

Let ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ be a sequence of RB operators of the form (28) on $\left(\mathcal{B}\right(X,Y),d)$. Suppose that the elements of $\{q_{i_k,k}:i_k\in {\mathbb{N}}_{n_k},k\in \mathbb{N}\}$ satisfy:

$$\underset{k\in \mathbb{N}}{sup}\underset{i_k\in {\mathbb{N}}_k}{max}d(q_{i_k,k},0)\le M,$$

(31)

for some $M>0$. Then, the ball $B_r\left(0\right)$ of radius $r=M/(1-s)$ centered at $0\in \mathcal{B}(X,Y)$ is an invariant set for ${\left\{T_k\right\}}_{k\in \mathbb{N}}$.

Proof. 

Note that since Y is an F-space, we have for all $a,b\in Y$,

$$d_Y(a+b,0)\le d_Y(a+b,b)+d_Y(b,0)=d_Y(a,0)+d_Y(b,0).$$

Now, let $x\in X$. Then, there exists an $i_k\in {\mathbb{N}}_{n_k}$ with $x\in l_{i_k,k}\left(X\right)$. Thus, for any $f\in \mathcal{B}(X,Y)$,

$$\begin{array}{cc}\hfill d_Y(T_{k}f\left(x\right),0)& \le d_Y(S_{i_k,k}\circ l_{i_k,k}^{-1}\left(x\right)·f\circ l_{i_k,k}^{-1}\left(x\right),0)+d_Y(T_k\left(0\right),0).\hfill \end{array}$$

By (31), $T_k\left(0\right)$ is uniformly bounded in $\mathcal{B}(X,Y)$ by $M>0$. As the metric $d_Y$ is homogeneous,

$$d_Y(S_{i_k,k}\circ l_{i_k,k}^{-1}\left(x\right)·f\circ l_{i_k,k}^{-1}\left(x\right),0)\le s{d}_Y(f\circ l_{i_k,k}^{-1}\left(x\right),0),$$

which shows, after taking the sup over $x\in X$, that $d(T_{k}f,0)\le sd(f,0)+M$. Proposition 1 now yields the statement. □

Considering the backward trajectories ${\left\{{\mathsf{\Psi }}_k\right\}}_{k\in \mathbb{N}}$ of the sequence ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ of RB operators defined above and using Theorem 1, we obtain the next result.

Theorem 3.

The backwards trajectories ${\left\{{\mathsf{\Psi }}_k\right\}}_{k\in \mathbb{N}}$ converge for any initial $f_0\in \mathcal{I}$ to a unique attractor $f^{*}\in \mathcal{I}$, where $\mathcal{I}$ is the closed ball in $\mathcal{B}(X,Y)$ of radius $M/(1-s)$ centered at zero.

Proof. 

By Theorem 1, it remains to show that ${\displaystyle \sum _{k=1}^{\infty }}{\displaystyle \prod _{j=1}^k}\mathrm{Lip}\left(T_j\right)$ converges. This, however, follows directly from (30):

$$\prod _{j=1}^k\mathrm{Lip}\left(T_j\right)\le s^k\mathrm{and}\sum _{k=1}^{\infty }{s}^k=\frac{s}{1-s}.$$

 □

A fixed point $f^{*}$ generated by a sequence $\left\{T_k\right\}$ of different RB operators will be called a non-stationary fractal function (of class $\mathcal{B}(X,Y)$).

Remark 10.

Item (b) in Remark 3, of course, also applies to a sequence of RB operators $\left\{T_k\right\}$, thus allowing the construction of more general fractal functions exhibiting different local behavior at different scales.

Example 1.

Let $X:=[0,1]$ and $Y:=\mathbb{R}$. Consider the two RB operators:

$$T_{1}f\left(x\right):=\left\{\begin{array}{cc}2x+\frac{1}{2}f\left(2x\right),\hfill & x\in [0,\frac{1}{2}),\hfill \\ 2-2x+\frac{1}{2}f(2x-1),\hfill & x\in [\frac{1}{2},1],\hfill \end{array}\right.$$

and:

$$T_{2}f\left(x\right):=\left\{\begin{array}{cc}2x+\frac{1}{4}f\left(2x\right),\hfill & x\in [0,\frac{1}{2}),\hfill \\ 2-2x+\frac{1}{4}f(2x-1),\hfill & x\in [\frac{1}{2},1].\hfill \end{array}\right.$$

For both operators, $l_i\left(x\right):=\frac{1}{2}(x+i-1)$, $i=1,2$.

It is known that $T_1^{k}f\to \tau$, where τ denotes the Takagi function [11], and that $T_2^k\to q$, where $q\left(x\right)=4x(1-x)$.

Consider the alternating sequence ${\left\{T_i\right\}}_{i\in \mathbb{N}}$ of RB operators given by:

$$T_i:=\left\{\begin{array}{cc}{T}_1,\hfill & 10(j-1)<i\le 10j-5,\hfill \\ T_2,\hfill & 10j-5<i\le 10j,\hfill \end{array}\right.j\in \mathbb{N}.$$

Two images of this hybrid attractor of the backward trajectory ${\mathsf{\Psi }}_k$ starting with $f_0\equiv 0$ are shown in Figure 1.

6. Non-Stationary Fractal Interpolation

Let us now consider the case $X:=[0,1]$ and $Y:=\mathbb{R}$. Both spaces are metrizable under the usual Euclidean distance. In the following, we consider a sequence $\left\{T_k\right\}$ of RB operators of the form (25) acting on an appropriate metric subspace of $\mathcal{B}[0,1]:=\mathcal{B}\left(\right[0,1],\mathbb{R})$. Our emphasis here lies in the construction of attractors that are continuous functions on $[0,1]$. For this purpose, we need to impose conditions on the RB operators that guarantee global continuity of the iterates across $[0,1]$.

For $k\in \mathbb{N}$, let $\{l_{i_k,k}:i_k\in {\mathbb{N}}_{n_k},k\in \mathbb{N}\}$ be a family of injections from $[0,1]\to [0,1]$ generating a partition of $[0,1]$ in the sense of (8) and (9). Assume w.l.o.g. that $l_{1,k}\left(0\right)=0$ and $l_{n_k,k}\left(1\right)=1$, and define:

$$\begin{array}{c}\hfill x_{i_k-1,k}:=l_{i_k,k}\left(0\right),x_{i_k,k}:=l_{i_k,k}\left(1\right),i_k\in {\mathbb{N}}_{n_k},\end{array}$$

where $x_{0,k}:=0$ and $x_{n_k,k}:=1$. By relabeling, if necessary, we may assume that $0=x_{0,k}<\cdots <x_{i_k-1,k}<x_{i_k,k}<\cdots x_{n_k,k}=1$.

Let $f\in \mathcal{C}[0,1]$ be arbitrary. Define a metric subspace of $\mathcal{C}[0,1]$ by:

$${\mathcal{C}}_{*}[0,1]:=\{g\in \mathcal{C}[0,1]:g\left(0\right)=f\left(0\right)\wedge g\left(1\right)=f\left(1\right)\}$$

and note that ${\mathcal{C}}_{*}[0,1]$ becomes a complete linear metric space when endowed with the metric induced by the sup-norm on continuous functions. Additionally, let $b\in {\mathcal{C}}_{*}[0,1]$ be the unique affine function whose graph connects the points $(0,f(0\left)\right)$ and $(1,f(1\left)\right)$:

$$b\left(x\right)=\left(f\right(1)-f(0\left)\right)x+f\left(0\right).$$

(32)

Further, let ${\left\{{\mathcal{P}}_k\right\}}_{k\in \mathbb{N}}$, where ${\mathcal{P}}_k:=\{(x_{j_k},f\left(x_{j,k}\right))\in [0,1]\times \mathbb{R}:j=0,1,\dots ,n\}$, be a family of sets of points in $[0,1]\times \mathbb{R}$. For $k\in \mathbb{N}$, define an RB operator $T_k:{\mathcal{C}}_{*}[0,1]\to {\mathcal{C}}_{*}[0,1]$ by:

$$T_{k}g=f+\sum _{i_k=1}^{n_k}{S}_{i_k,k}\circ l_{i_k,k}^{-1}·(g-b)\circ l_{i_k,k}^{-1}{\chi }_{l_{i_k,k}[0,1]},$$

(33)

where ${\left\{S_{i_k,k}\right\}}_{i_k=1}^{n_k}\subset \mathcal{C}[0,1]$ such that

$$\underset{k\in \mathbb{N}}{sup}\underset{i_k\in {\mathbb{N}}_{i_k}}{\max }{\parallel S_{i_k,k}\parallel }_{\infty }<1.$$

Note that we have the continuity of $T_{k}g$ at the points $x_{i_k,k}\in [0,1]$:

$$T_{k}g(x_{i_k,k}-)=T_{k}g(x_{i_k,k}+),\forall i_k\in \{1,\dots ,n-1\}.$$

For,

$$\begin{array}{cc}\hfill T_{k}g(x_{i_k,k}-)& =f(x_{i_k,k}-)+S_{i_k,k}\circ l_{i_k,k}^{-1}(x_{i_k,k}-)·(g-b)\circ l_{i_k,k}^{-1}(x_{i_k,k}-)\hfill \\ & =f\left(x_{i_k,k}\right)+S_{i_k,k}\left(1\right)·(f-b)\left(1\right)=f\left(x_{i_k,k}\right)\hfill \end{array}$$

and

$$\begin{array}{cc}\hfill T_{k}g(x_{i_k,k}+)& =f(x_{i_k,k}+)+S_{i_k+1,k}\circ l_{i_k+1,k}^{-1}(x_{i_k,k}+)·(g-b)\circ l_{i_k+1,k}^{-1}(x_{i_k,k}+)\hfill \\ & =f\left(x_{i_k,k}\right)+S_{i_k+1,k}\left(0\right)·(f-b)\left(0\right)=f\left(x_{i_k,k}\right).\hfill \end{array}$$

Therefore, $T_{k}g\in {\mathcal{C}}_{*}[0,1]$, and $T_{k}g$ interpolates ${\mathcal{P}}_k$ in the sense that:

$$T_{k}g\left(x_{i_k,k}\right)=f\left(x_{i_k,k}\right),\forall i_k\in {\mathbb{N}}_{n_k}.$$

Remark 11.

Denote by $([0,1],{\mathcal{L}}_k)$ the IFS given by the maps ${\mathcal{L}}_k:=\{l_{i_k,k}:i_k\in {\mathbb{N}}_{n_k}\}$, and observe that, for each $k\in \mathbb{N}$, the attractor of $([0,1],{\mathcal{L}}_k)$ is the interval $[0,1]$. The invariant set, in $\mathcal{H}\left(\right[0,1\left]\right)$, for ${\mathcal{L}}_k$ is given by $[0,1]$. Hence, all backward trajectories ${\mathcal{L}}_1\circ \cdots \circ {\mathcal{L}}_k$ converge to $[0,1]$ as $k\to \infty$ (as do all forward trajectories).

Proposition 3.

A nonempty closed invariant set for ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ is given by the closed ball in ${\mathcal{C}}_{*}[0,1]$,

$$\mathcal{I}=\left\{g\in {\mathcal{C}}_{*}[0,1]:{\parallel g\parallel }_{\infty }\le \frac{{\parallel f\parallel }_{\infty }+s{\parallel b\parallel }_{\infty }}{1-s}\right\},$$

(34)

where s is given by (30).

Proof. 

Using the form (24) for the functions $q_{i_k,k}$, we obtain from (31) the estimate $\parallel q_{i_k,k}{\parallel }_{\infty }\le {\parallel f\parallel }_{\infty }+s{\parallel b\parallel }_{\infty }$, which by Proposition 2 yields the result. □

In connection with Theorem 3, the above arguments prove the next result.

Theorem 4.

Let ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ be a sequence of RB operators of the form (33), each of whose elements acts on the complete metric space $({\mathcal{C}}_{*}[0,1],d)$ where $f\in C_{*}[0,1]$ is arbitrary and b is given by (32). Further, let the family of functions $\left\{S_{i_k,k}\right\}\subset \mathcal{C}[0,1]$ satisfy (30). Then, the backward trajectories ${\mathsf{\Psi }}_k\left(f_0\right)$ converge to a function $f^{*}\in \mathcal{I}$, for any $f_0\in \mathcal{I}$. As $f_0$ one may choose f or b.

We refer to the fixed point $f^{*}\in C_{*}[0,1]$ as a continuous non-stationary fractal interpolation function.

To illustrate the above results, we refer to Remark 3(c) and present the following example.

Example 2.

Here, we consider the two RB operators $T_i:C[0,1]\to C[0,1]$, $i=1,2$, given by:

$$\left(T_{1}f\right)\left(x\right)=\left\{\begin{array}{cc}-\frac{1}{2}f\left(4x\right),\hfill & x\in [0,\frac{1}{4}),\hfill \\ -\frac{1}{2}+\frac{1}{2}f(4x-1),\hfill & x\in [\frac{1}{4},\frac{1}{2}),\hfill \\ \frac{1}{2}f(4x-2),\hfill & x\in [\frac{1}{2},\frac{3}{4}),\hfill \\ \frac{1}{2}+\frac{1}{2}f(4x-3),\hfill & x\in [\frac{3}{4},1],\hfill \end{array}\right.$$

and

$$\left(T_{2}f\right)\left(x\right):=\left\{\begin{array}{cc}\frac{3}{4}f\left(2x\right),\hfill & x\in [0,\frac{1}{2}),\hfill \\ \frac{3}{4}+\frac{1}{4}f(2x-1),\hfill & x\in [\frac{1}{2},1].\hfill \end{array}\right.$$

The RB operators $T_1$ and $T_2$ generate Kiesswetter’s fractal function [12] and a Casino function [13], respectively.

Consider again the alternating sequence ${\left\{T_i\right\}}_{i\in \mathbb{N}}$ of RB operators given by:

$$T_i:=\left\{\begin{array}{cc}{T}_1,\hfill & 10(j-1)<i\le 10j-5,\hfill \\ T_2,\hfill & 10j-5<i\le 10j,\hfill \end{array}\right.j\in \mathbb{N}.$$

Two images of the hybrid attractor of the backward trajectory ${\mathsf{\Psi }}_k$ starting with the function $f_0\left(x\right)=x$, $x\in [0,1]$, are shown below in Figure 2.

Remark 12.

Theorem 2 holds in the case of non-stationary fractal functions as well. For $k\in \mathbb{N}$, a non-stationary IFS is associated with $T_k$ by setting:

$$w_{i_k,k}(x,y):=(l_{i_k,k}\left(x\right),f\circ l_{i_k,k}\left(x\right)+S_{i_k,k}\left(x\right)·(y-b)).$$

The conditions imposed on $S_{i_k,k}$ and the form of the second component allow the immediate transfer of the proof of Theorem 2. Hence, even in the non-stationary case, one may choose the geometry (IFS) or the analytic (RB operator) approach when defining non-stationary fractal functions.

7. Non-Stationary Fractal Functions in Bochner–Lebesgue Spaces

In this section, we construct non-stationary fractal functions in the Bochner–Lebesgue spaces ${\mathcal{L}}^p$ with $0<p\le \infty$. To this end, suppose that X is a closed subspace of a Banach space $\mathsf{X}$ and $\mathbb{X}:=(X,\Sigma ,\mu )$ a measure space. Further, suppose that $(\mathsf{Y},\parallel ·{\parallel }_{\mathsf{Y}})$ is a Banach space.

Recall that the Bochner–Lebesgue space ${\mathcal{L}}^p(\mathbb{X},\mathsf{Y})$, $1\le p\le \infty$ is comprised of all Bochner measurable functions $f:X\to \mathsf{Y}$ such that:

$${\parallel f\parallel }_{{\mathcal{L}}^p(\mathbb{X},\mathsf{Y})}:={\left({\int }_X{\parallel f\left(x\right)\parallel }_{\mathsf{Y}}^{p}d\mu \left(x\right)\right)}^{1/p}<\infty ,1\le p<\infty ,$$

and

$${\parallel f\parallel }_{{\mathcal{L}}^{\infty }(\mathbb{X},\mathsf{Y})}:={ess\; sup}_{x\in X}{\parallel f\left(x\right)\parallel }_{\mathsf{Y}}<\infty ,p=\infty .$$

For $0<p<1$, the spaces ${\mathcal{L}}^p(\mathbb{X},\mathsf{Y})$ are defined using a metric instead of a norm to obtain completeness. More precisely, for $0<p<1$, define $d_p:{\mathcal{L}}^p(\mathbb{X},\mathsf{Y})\times {\mathcal{L}}^p(\mathbb{X},\mathsf{Y})\to \mathbb{R}$ by:

$$d_p(f,g):={\parallel f-g\parallel }_{\mathsf{Y}}^p.$$

Then, $({\mathcal{L}}^p(\mathbb{X},\mathsf{Y}),d_p)$ becomes an F-space (note that the inequality ${(a+b)}^p\le a^p+b^p$ holds for all $a,b\ge 0$). For more details, we refer to [14,15].

In order to work in both cases simultaneously, we define ${\rho }_p:{\mathcal{L}}^p(\mathbb{X},\mathsf{Y})\times {\mathcal{L}}^p(\mathbb{X},\mathsf{Y})\to \mathbb{R}$ by:

$${\rho }_p(g,h):=\left\{\begin{array}{cc}{\parallel g-h\parallel }_{{\mathcal{L}}^p(\mathbb{X},\mathsf{Y})},\hfill & 1\le p\le \infty ,\hfill \\ {\parallel g-h\parallel }_Y^p,\hfill & 0<p<1,\hfill \end{array}\right.$$

with the usual modification for $p=\infty$.

We use the notation and terminology of Section 5 and assume that:

(A1)

$\{q_{i_k,k}:i_k\in {\mathbb{N}}_{n_k},k\in \mathbb{N}\}\subset {\mathcal{L}}^p(\mathbb{X},\mathsf{Y})$;

(A2)

$\{S_{i_k,k}:i_k\in {\mathbb{N}}_{n_k},k\in \mathbb{N}\}\subset {\mathcal{L}}^p(\mathbb{X},\mathbb{R})$;

(A3)

$\{l_{i_k,k}:i_k\in {\mathbb{N}}_{n_k},k\in \mathbb{N}\}$ is a family of $\mu$-measurable diffeomorphisms $X\to X$ generating for each $k\in \mathbb{N}$ a partition of X in the sense of (8) and (9).

If we define for each $k\in \mathbb{N}$ an RB operator $T_k$ on ${\mathcal{L}}^p(\mathbb{X},\mathsf{Y})$ of the form (28), whose maps satisfy assumptions (A1), (A2), and (A3), then a straightforward computation shows that $T_k$ has the following Lipschitz constants on ${\mathcal{L}}^p(\mathbb{X},\mathsf{Y})$:

$$\begin{array}{c}\hfill {\rho }_p(T_{k}g,T_{k}h)\le \left\{\begin{array}{c}{\left({\displaystyle \sum _{i_k=1}^{n_k}}{\parallel S_{i_k,k}\parallel }_{{\mathcal{L}}^p(\mathbb{X},\mathsf{Y})}^p·L_{i_k,k}\right)}^{1/p},(1\le p<\infty )\\ \underset{i_k\in {\mathbb{N}}_{n_k}}{max}{\parallel S_{i_k,k}\parallel }_{{\mathcal{L}}^{\infty }(\mathbb{X},\mathsf{Y})},(p=\infty )\\ {\displaystyle \sum _{i_k=1}^{n_k}}{\parallel S_{i_k,k}\parallel }_{{\mathcal{L}}^p(\mathbb{X},\mathsf{Y})}^p·L_{i_k,k},(0<p<1)\end{array}\right\}{\rho }_p(g,h),\end{array}$$

where $L_{i_k,k}$ denotes the Lipschitz constant of $D{l}_{i_k,k}^{-1}$ and D the Fréchet derivative on X.

Now, set:

$${\gamma }_p:=\left\{\begin{array}{cc}\underset{k\in \mathbb{N}}{sup}{\left({\displaystyle \sum _{i_k=1}^{n_k}}{\parallel S_{i_k,k}\parallel }_{{\mathcal{L}}^p(\mathbb{X},\mathsf{Y})}^p·L_{i_k,k}\right)}^{1/p},\hfill & 1\le p<\infty \hfill \\ \underset{k\in \mathbb{N}}{sup}\underset{i_k\in {\mathbb{N}}_{n_k}}{max}{\parallel S_{i_k,k}\parallel }_{{\mathcal{L}}^{\infty }(\mathbb{X},\mathsf{Y})},\hfill & p=\infty \hfill \\ \underset{k\in \mathbb{N}}{sup}\left({\displaystyle \sum _{i_k=1}^{n_k}}{\parallel S_{i_k,k}\parallel }_{{\mathcal{L}}^p(\mathbb{X},\mathsf{Y})}^p·L_{i_k,k}\right),\hfill & 0<p<1.\hfill \end{array}\right.$$

(35)

Imposing the condition:

$$\underset{k\in \mathbb{N}}{sup}\underset{i_k\in {\mathbb{N}}_{n_k}}{max}{\rho }_p(q_{i_k,k},0)<M,$$

(36)

for some $M>0$ and further requiring that:

$$\mathrm{Lip}{T}_k\le {\gamma }_p<1,\forall k\in \mathbb{N},$$

(37)

yield by Proposition 2 an invariant set for ${\left\{T_k\right\}}_{k\in \mathbb{N}}$, namely the closed ${\mathcal{L}}^p$-ball:

$$\mathcal{I}=B_r\left(0\right)\mathrm{with}r=M/(1-{\gamma }_p).$$

The above elaborations now prove the following theorem.

Theorem 5.

Let ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ be a sequence of RB operators of the from (28) mapping ${\mathcal{L}}^p(\mathbb{X},\mathsf{Y})$ into itself. Further, suppose that the Lipschitz constant of $T_k$ satisfies (37) and that the maps $\left\{q_{i_k,k}\right\}$ fulfill (36). Then, the backward trajectories ${\left\{{\mathsf{\Psi }}_k\right\}}_{k\in \mathbb{N}}$ of ${\left\{T_k\right\}}_{k\in \mathbb{N}}$ converge for any initial $f_0\in \mathcal{I}$ to a unique attractor $f^{*}\in \mathcal{I}$, where $\mathcal{I}$ is the ball in $\mathcal{L}(\mathbb{X},\mathsf{Y})$ of radius $M/(1-{\gamma }_p)$ centered at zero.

Proof. 

Only (7) needs to be established. This, however, carries over directly from the proof of Theorem 3 with ${\gamma }_p$ instead of s. □

The attractor $f^{*}:X\to \mathsf{Y}$ whose existence is guaranteed by Theorem 5 is called a non-stationary fractal function of class ${\mathcal{L}}^p(\mathbb{X},\mathsf{Y})$.