Fractal Interpolation Using Harmonic Functions on the Koch Curve

: The Koch curve was ﬁrst described by the Swedish mathematician Helge von Koch in 1904 as an example of a continuous but nowhere differentiable curve. Such functions are now characterised as fractal since their graphs are in general fractal sets. Furthermore, it can be obtained as the graph of an appropriately chosen iterated function system. On the other hand, a fractal interpolation function can be seen as a special case of an iterated function system thus maintaining all of its characteristics. Fractal interpolation functions are continuous functions that can be used to model continuous signals. An in-depth discussion on the theory of afﬁne fractal interpolation functions generating the Koch Curve by using fractal analysis as well as its recent development including some of the research made by the authors is provided. We ensure that the graph of fractal interpolation functions on the Koch Curve are attractors of an iterated function system constructed by non-constant harmonic functions.


Introduction
The Koch curve appeared in a 1904 paper entitled "On a Continuous Curve Without Tangents, Constructible from Elementary Geometry" by the Swedish mathematician Helge von Koch. The Sierpiński gasket was introduced in 1915 by the Polish mathematician W. Sierpiński, about forty years after the discovery of the Cantor set. An interval, the Sierpinski Gasket, or SG for short, and the Koch Curve, or KC for short, are connected typical self-similar sets (see [1]). The functions or the curves that are continuous but nowhere differentiable are now called fractal functions or fractal curves; see [2], p. 46. An interval is not a fractal curve, but it is self-similar (see [3]). Although a fractal interpolation function, or FIF for short, is usually defined on line segments, there is insufficient discussion of a FIF defined on fractal curves like the SG or the KC.
In [4] the authors showed how one can construct space-filling curves by using hidden variable linear fractal interpolation functions. These curves resulted from the projection of the attractor of an iterated function system, or IFS for short. In [5] the authors showed how the theory of linear fractal interpolation functions together with the Deterministic Iteration Algorithm can be used to construct space-filling curves. Interpreting the polynomials of degree 1 as classical harmonic functions on an interval and replacing them on the KC by harmonic functions, the authors of [6] obtained an analogue of Theorem 2.2. Chapter VI of [7] for the KC. The authors of [8] showed how it is possible to generalise a fractal interpolation problem to certain post critically finite, or PCF for short, compact sets in R n by using harmonic functions to solve this fractal interpolation problem. Since non-constant harmonic functions on KC are not Lipschitz continuous (see [8]), to prove the uniqueness of invariant sets of IFSs on KC is not possible by using the classic fractal interpolation method. It is well known that fractal analysis is assessing fractal characteristics of data. The results of [8] and [6] enable us to study Hölder continuity of non-constant harmonic functions on KC and inspire us to ensure that graphs of FIFs generated on KC by non-constant harmonic functions of fractal analysis are attractors of some IFSs.
Fractal interpolation functions generated on some special affine fractal interpolation curve by harmonic functions of fractal analysis are given in [9]. The authors did not give all the details about the KC such as new and suitable IFS constructing KC, something very suitable in the proof of Hölder continuity of non-constant harmonic functions on KC. Although an affine fractal interpolation curve can be seen as the graph of some continuous function, the KC is not the graph of a continuous function but only a continuous image (more precisely, a continuous map) from a segment of a line to R 3 . However, it is possible to ensure that graphs of fractal interpolation functions on the KC are attractors of iterated function systems. So, fractal interpolation curves of R 3 generated on KC are much more important than the one on affine fractal interpolation curves. We should use another IFS for the generation of the KC because the proofs of existence and uniqueness of attractor of an IFS constructed on KC are dependent on the first point (or end point) of each curve. For this reason, there is an important difference between the proof of uniqueness of invariant set of IFS on KC and the one on a special affine fractal interpolation curve.
In this article, by using the important fact that the existence and uniqueness of invariant sets of IFSs on KC is closely related to the same suitable metric associated with Hölder exponent of non-constant harmonic functions on KC, that is, through strict accurate mathematical description for an IFS constructing KC, we improve and clarify some results for fractal interpolation on KC. The rest of this article is organised as follows and can be generally seen as an organised study and especially as an in-depth discussion of Corollary 4.8, Theorem 4.10 and Corollary 4.11 presented in [9]. In Section 2 we recall some already known results, and we give certain IFSs on KC. In Section 3 we give fractal interpolation as attractors of IFSs constructed on KC by non-constant harmonic functions of fractal analysis. In Section 4 we discuss and cite research areas presented in international journals dealing with similar topic. Finally, in Section 5 we draw our conclusions.

Fractal Interpolation on a Line Segment
In this section, we review some already known results about fractal interpolation on a closed interval in order to derive interpolation functions as attractors of IFSs constructed on KC by non-constant harmonic functions of fractal analysis; see [10], pp. 44-45 or [2], Definition 2.2, p. 44.
Let N be a positive integer greater than 1, where |s n | < 1 (see [11], p. 344, see [7], p. 214, see [12], p. 308). Then each h n is a unique classic nonconstant harmonic function on an interval [x 0 , x N ] such that, see [13], Definition 1. Let (X, ρ) be a metric space and α ∈ (0, 1]. A function f : X → X is said to be Hölder continuous on X with respect to ρ and α, if there exists a constant k ≥ 0 such that for all x, y ∈ X, The above constant k is called the Hölder constant of a function f . Let Then we call α f ∈ (0, 1] the Hölder exponent of a function f (cf. [14], p. 36). If α f = 1, then a function f is said to be Lipschitz continuous on X. A function f is said to be Banach contraction, if α f = 1 and 0 ≤ k < 1.
Proof. The proof is similar to the one of [3] (or [1]) and is omitted here.
The graphs of harmonic functions on K 2 , K 3 and Koch curve (see  are given below. Compare and contrast them with Figure 4 of page 3233 in [6], where a graph of 1-nonconstant harmonic function on KC can be seen. Without loss of generality, we may assume that n = 1 (see [15], p. 314).
The following motivation will be the key in the proof of Theorem 6.
(1) Piecewise linear interpolation functions on an interval are special FIFs that are Lipschitz continuous (see [7], p. 212, p. 214), and so Lipschitz continuity with respect to the first variable x of functions F n (x, y) and metrically equivalent metric are used in the proof of the uniqueness of invariant of IFS (see [7], p. 217, Theorem 1). (2) Harmonic functions on KC are special FIFs on KC (cf. [6], cf. [8]).
(3) The graph of a constant harmonic function f (x, y) ≡ c, (x, y) ∈ KC is KC itself, and Hausdorff dimension of KC is log 4 log 3 , (see [2], p. 135, see [3]), and so box dimensions of graphs of non-constant harmonic functions on KC can be non-integers, and nonconstant harmonic functions on KC can be Hölder continuous. (4) Since non-constant harmonic functions on KC are not Lipschitz continuous (see [8], p. 36), to obtain FIFs on KC as attractors of some IFSs is not possible from Barnsley's fractal interpolation method. (5) The assertion that two metrics d 1 and d 2 are metrically equivalent is much stronger than the statement that they are topologically equivalent: to be metrically equivalent there must exist constants c 1 and c 2 such that for all x, y ∈ R n , and to be topologically equivalent must ensure that a sequence {(x n , y n )} ∞ n=1 ⊂ R n which is d 1 -convergent to (x 0 , y 0 ) ∈ R n is also d 2 -convergent to (x 0 , y 0 ) ∈ R n , and a sequence {(x n , y n )} ∞ n=1 ⊂ R n which is d 2 -convergent to (x 0 , y 0 ) ∈ R n is also d 1 -convergent to (x 0 , y 0 ) ∈ R n . (6) There can exist some metric d on R n , topologically equivalent (not necessarily metrically equivalent) to the Euclidean metric, such that IFS has a unique invariant set because the existence of metrically equivalent metric is a sufficient condition so that functions w i (x, y) are Banach contractions with respect to some metric. (7) The results of [6,8] inspire us to ensure that graphs of FIFs generated on KC by non-constant harmonic functions of fractal analysis are attractors of some IFSs.
Our idea, based on the above motivation, is to use Hölder continuity of non-constant harmonic functions on KC and use some suitable metric which is topologically equivalent to the Euclidean metric but is not metrically equivalent to that.
(1) In the case of affine FIFs, functions F n (x, y) are Lipschitz continuous with respect to the first variable x and are Banach contractions with respect to the second variable y because and |F n (x, y ) − F n (x, y )| = |s n y + h n (x) − s n y + h n (x)| ≤ |s n | |y − y |.
(2) In the formulation and proof of Theorem 2, it is enough to assume that functions F n (x, y) are continuous with respect to the first variable x and are Banach contractions with respect to the second variable y (see Theorem 2 in p. 218 of [7]), and in the formulation and proof of Theorem 5, it is enough to assume that functions F n (x, y, z) are continuous with respect to the first variable x and second variable y and are Banach contractions with respect to the third variable z. {KC × R; w n , n = 1, 2, 3, 4} have unique invariant sets. The uniqueness of invariant sets is determined explicitly in Theorems 3 and 6 (cf. the Theorem 1 in p. 217 of [7]). (4) In the formulation and proof of Theorem 6, it is needed to consider functions F n (x, y, z) that are Hölder continuous with respect to the first variable x and second variable y and are Banach contractions with respect to the third variable z. (5) Theorem 3 is dedicated to show a sufficient condition for functions F n (x, y) such that the IFS{[x 0 , x N ] × R; w n , n = 1, 2, . . . , N} has a unique invariant set, and Theorem 6 is dedicated to show an essential sufficient condition for functions F n (x, y, z) such that the IFS {KC × R; w n , n = 1, 2, 3, 4} has a unique invariant set.

Discussion
One of the objectives would be to investigate how to optimise such a process in order to enhance mechanical properties of materials produced using the afore-mentioned method. The article [17] reports on the results of three studies investigating visual interest, visual preference, and mood responses elicited by varying complexities of fractal light patterns projected on walls and floors of an interior space. The article [18] presents a study on the design, and microstructural and mechanical characterization of additively manufactured reinforcing elements for composite materials exhibiting fractal geometry, with a focus on the flexural reinforcement of cement-matrix composites.

Conclusions
In this article, we used Hölder continuity of nonconstant harmonic functions on KC and a suitable metric, equivalent to the Euclidean metric, to prove Theorem 6. Moreover, we showed that it is possible to ensure that the graph of an FIF on KC by nonconstant harmonic functions of fractal analysis is the attractor of some IFS. An important fact is that the proof of the existence of fractal interpolation functions on the Koch curve is closely related to the same suitable iterated function system that generates the Koch curve, whereas a key result is that the proof of the uniqueness of FIFs depends not only on Hölder continuity of harmonic functions on the Koch curve but also on the topologically equivalent metric.