Next Article in Journal
A Robust High-Dimensional Test for Two-Sample Comparisons
Next Article in Special Issue
Revised and Generalized Results of Averaging Principles for the Fractional Case
Previous Article in Journal
Generalized Partially Functional Linear Model with Interaction between Functional Predictors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Positivity-Preserving Rational Cubic Fractal Interpolation Function Together with Its Zipper Form

by
Shamli Sharma
1,
Kuldip Katiyar
1,
Gadug Sudhamsu
2,
Manjinder Kaur Wratch
3 and
Rohit Salgotra
4,5,6,*
1
Department of Mathematics, Chandigarh University, Gharuan, Mohali 140413, Punjab, India
2
Department of Computer Science and Engineering, School of Engineering and Technology, JAIN (Deemed to be University), Bangalore 560069, Karnataka, India
3
Department of Computer Science Engineering, Chandigarh Engineering College, Chandigarh Group of Colleges, Jhanjeri, Mohali 140307, Punjab, India
4
Faculty of Physics and Applied Computer Science, AGH University of Krakow, 30-059 Krakow, Poland
5
MEU Research Unit, Middle East University, Amman 11813, Jordan
6
Data Science Institute, University of Technology Sydney, Sydney, NSW 2007, Australia
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(9), 584; https://doi.org/10.3390/axioms13090584
Submission received: 17 June 2024 / Revised: 8 August 2024 / Accepted: 9 August 2024 / Published: 28 August 2024

Abstract

:
In this paper, a novel class of rational cubic fractal interpolation function (RCFIF) has been proposed, which is characterized by one shape parameter and a linear denominator. In interpolation for shape preservation, the proposed rational cubic fractal interpolation function provides a simple but effective approach. The nature of shape preservation of the proposed rational cubic fractal interpolation function makes them valuable in the field of data visualization, as it is crucial to maintain the original data shape in data visualization. Furthermore, we discussed the upper bound of error and explored the mathematical framework to ensure the convergence of RCFIF. Shape parameters and scaling factors are constraints to obtain the desired shape-preserving properties. We further generalized the proposed RCFIF by introducing the concept of signature, giving its construction in the form of a zipper-rational cubic fractal interpolation function (ZRCFIF). The positivity conditions for the rational cubic fractal interpolation function and zipper-rational cubic fractal interpolation function are found, which required a detailed analysis of the conditions where constraints on shape parameters and scaling factor lead to the desired shape-preserving properties. In the field of shape preservation, the proposed rational cubic fractal interpolation function and zipper fractal interpolation function both represent significant advancement by offering a strong tool for data visualization.

1. Introduction

Interpolation is the process of finding the value of a function at a point from its neighboring points, or it is a method for generating new data points within the range of a discrete set of known data points. In numerous areas of scientific computation, especially computer graphics, the science of biology and physical science, and the field of geology and numerical computation, interpolation techniques are necessary, specifically when it comes to repeatedly visualizing discrete data. In the field of science and engineering, visualizing discrete scientific data in continuous form is crucial. Based on their graphical representation, the collected data from any experiment are broadly categorized as positive, monotonic, concave, complex, or constrained by their combination, curve, and surface. For example, in chemical experiments, the quantity of obtained product is always positive, and in the case of metals, the resistivity with temperature monotonically increases in nature. Also, with a certain given initial angle and velocity, the trajectory of a projectile is always concave, and in the case of alternating current, the amplitude during some specific time interval shows both concave and convex properties. For representing discrete data as a continuous curve, splines have proven to be extremely valuable. For the initial development of spline theory for various types of data polynomial splines studied, see [1,2,3,4].
Classical interpolants are suitable for approximating functions with irregular data sets because some of the derivatives of classical interpolants are globally smooth, and some are piecewise; however, it is not possible to preserve the hidden shape properties of the data. This is why it is not an ideal method for designing irregular curves and surfaces. When non-desirable oscillations occur, like a pendulum’s motion on a cart placed in an electrochemical system [5] or a spherical ball falling in a hot micellar solution [6], the observed results are extremely misguided, violating the inherited nature of the data; meanwhile, fractal interpolations are ideal methods for irregularity. Apart from the theoretical significance of fractal interpolation, fractal interpolation has a self-similarity property and the ability to possess fractality in the functions or their respective derivatives, which helps them to approximate nonlinear phenomena accurately. Fractal interpolation is an advanced and modern technique that can be used to analyze numerous scientifically obtained data from some scientific phenomena and unknown complicated functions. Fractal interpolation not only helps construct rough structures but also ensures that smooth structures can be constructed. Depending on the values of parameters of iterated function systems, such structures can have a non-integer dimension [7,8,9].
The word “fractal” was coined by the mathematician Mandelbrot in 1975 to describe sophisticated and irregular natural objects and to better understand some scientific experiments that have the presence of self-similarity. The roots of fractal theory began to form in the 17th century when Leibniz mistakenly considered the case that self-similarity only occurs in straight lines. Karl Weierstrass, on 18 July 1872, presented a paper in which he gave the first definition for a function that was continuous everywhere but differentiable nowhere. Soon afterward, Georg Cantor, in 1883, presented examples of a subset of the real line, known as the Cantor set, which had some irregular properties that made it a fractal. At the end of that century, Henri Poincare and Felix Klein introduced another type of fractal, known as a self-inverse fractal. In 1904, Helge Von Koch extended Poincare’s idea and was also not satisfied with the abstract and analytic definition given by Weierstrass. He presented a definition that was more geometric, along with a hand drawn image of such a function; these functions are now known as Koch snowflakes. Wacław Sierpinski made his contribution in 1915 by presenting the triangle and carpet. The fractal interpolation function was presented by Barnsley [8,10] in 1986. There, he introduced a real-valued interpolation function, which was defined on a compact interval of R. This function works well with objects that are present in nature, with some geometrical self-similarity. Ever since, fractal interpolation has been attracting more attention, as it can be used in almost every field of life science, including medicine, visualization, engineering science, finance, and other areas. Fractal interpolation has a wide area of use in the real world; it is used in movies, medical science, image compression, etc. Knowledge of fractals and ecosystems together is used to determine the acidrain spread of smoke. The best method to construct the fractal interpolation function is using the theory of IFS [11]. Harrington [12] introduced the k-times differentiable fractal interpolation function (FIF). The smooth FIF [13] was shown by Navascués in 2006. Many researchers have worked on shape-preserving smooth FIFs since they are useful in many scientific and economic areas such as medicine, stock markets, etc. Rational splines are one of the most used shape-preserving interpolants as they are more flexible [14,15,16,17,18]. M. A. Navascués introduced the concept of α fractal function [19], and M. Nasim Akhtar, together with M. A. Navascués [20], introduced the local non-affine α fractal interpolation function. In recent years, many authors have worked on different areas of fractal theory, such as fractal dimensions and fractal surfaces [21,22,23,24,25,26,27,28]. Recently, Bilel Selmi discussed the dimension of fractal functions on the product of the Sierpinski gasket [29]. M. A. Navascués [30] discussed a new type of contraction and the existence and uniqueness theorem of this map. In Ref. [31], the fractal interpolation function using the Kannan contraction is discussed, and the generation of fractals using the IFS of the Kannan contraction in a controlled metric space is presented in [32].
Shape preservation is one of the most important properties of fractal interpolation. The geometrical characteristics of the interpolating data are maintained using shape preservation. Shape preservation involves three key properties of functions: positivity, monotonicity, and convexity. All these three properties are essential in numerous real-world applications as they can define or maintain the curvature and nature of functions or data within their domain. Positivity is important in various real-world applications, including physical quantities such as mass, length, probability distributions, time series analysis, and more. The positivity of various classical and fractal interpolation functions has been discussed in the literature. For detailed studies, see [33,34,35,36,37] and the references therein. Similarly, monotonicity is another crucial property. This property ensures that the function follows a consistent trend (non-decreasing or non-increasing). Monotonicity is useful in preserving the trend and order of functions or data, which is crucial in maintaining the integrity of time series data and other applications. The monotonicity of various classical and fractal interpolation functions has been discussed in the literature; for detailed studies, see [38,39,40] and the references therein. The convexity is also equally crucial for many real-world applications, such as in the stock market. It helps in predicting market trends and behavior. The convexity property ensures that the function maintains a specific smoothness in the shape of the function or data. The convexity of various classical and fractal interpolation functions has been discussed in the literature; for detailed studies, see [41,42,43,44] and the references therein.
Basic constructions of the zipper have been reported [45,46]. The fractal interpolation functions using the zipper provide greater flexibility due to the presence of a binary vector, also known as a signature. Sangita Jha discussed the zipper quadratic fractal interpolation function and its positivity [47]. They also ensured that the discussed function is non-negative in the defined interval. Vijay and Chand [48] proposed the zipper fractal interpolation function considering variable scaling. By considering variable scaling, one can enhance the adaptability of the function. The convexity of the zipper cubic fractal interpolation function with a quadratic denominator has been discussed by Vijay and Chand [49]. In this paper, we present the construction of a zipper-rational cubic fractal interpolation function. Interpolation with constrained data has a wide application in real-world applications, such as removing undesired wiggles in prominent lines of automobile roofs, in engineering, preventing oscillations between points, and more.
The aim of this paper is to construct a family of shape-preserving rational cubic fractal interpolation functions, including their zipper form, and to demonstrate their advantage over the classical interpolation functions. The qualitative features of an interpolant curve depend mainly on the scaling factor and shape parameters. The constrained RCFIF is designed to ensure that it lies above the straight line for all real values that share this characteristic. This paper is organized as follows: In the Section 2, we provide a short review of the fractal interpolation function (FIF) and the zipper fractal interpolation function (ZFIF). First, we discussed the fundamental concepts and mathematical construction of the fractal interpolation function, along with its uniqueness and existence theorems. This section also highlights the difference between classical interpolation functions and fractal interpolation functions. Then, we discuss the zipper fractal interpolation function, including its basic construction and the presence of a signature, which enhances the flexibility of the interpolating function or data. In the Section 3, we present the construction of the rational cubic fractal interpolation function, and the zipper-rational cubic fractal interpolation function. For the rational cubic fractal interpolation function, we consider a rational function of the form u j ( θ ) v j ( θ ) , while for zipper-rational cubic fractal interpolation, the rational function is u j , ε ( θ ) v j , ε ( θ ) , where the numerator is a cubic polynomial, and the denominator is a linear polynomial with one shape parameter. In the Section 4, we discuss the convergence of the rational cubic fractal interpolation function and the zipper-rational cubic fractal interpolation function regarding the data-generating function. Their respective theorems are given with proofs. The convergence analysis includes the error bound and conditions that guarantee the precision and reliability of the considered interpolant. In the Section 5, we constrain the RCFIF and ZRCFIF so that they lie above the straight line. In the Section 6, sufficient conditions are met to obtain the positivity conditions for the rational cubic fractal interpolation function and zipper-rational cubic fractal interpolation function. We derive the mathematical framework to provide the necessary and sufficient conditions for the function to be non-negative. To support the theoretical part, we provide an example with graphical representations. The respective theorems for providing sufficient conditions for the rational cubic fractal interpolation function and zipper-rational cubic fractal interpolation function to be positive are discussed. In the Section 7, we conclude the work done in previous sections and discuss their results. We also suggest potential future work in this area, which can be done by extending the given function to higher-order dimensions and other possible extensions.

2. Basics of Fractal Interpolation Functions and Zipper Fractal Interpolation Functions

In this section, we provide a short review of fractal interpolation functions (FIFs) and zipper fractal interpolation functions (ZFIFs) together with their fundamental equation and uniqueness and existence theorems. First, we discuss the fundamental concepts and mathematical construction of the fractal interpolation function. This section also highlights the difference between classical interpolation functions and fractal interpolation functions. Then, we discuss the zipper fractal interpolation function, including its basic construction and the presence of a signature, which enhances the flexibility of the interpolating function or data.

2.1. Fractal Interpolation Function

Let ( y m , z m ) J × K , m M M = 1 , 2 , 3 M be the dataset, where K is a compact set of R and J = [ y 1 , y M ] satisfying y 1 < y 2 < y 3 < y M and J j = [ y j , y j + 1 ] for j M M 1 , M M 1 = 1 , 2 , , M 1 . Consider L j : J J j the contractive homomorphism defined as L j ( y ) = a j y + b j such that y , y * J and for some 0 l j < 1 .
L j ( y ) L j ( y * ) l j y y * ; L j ( y 1 ) = y j , L j ( y M ) = y j + 1 , j M M 1
Define F j : K : = J × c , d R such that, for and let δ j 1 , 1 ,
F j ( y , z ) F j y , z δ j z z * ; F j y 1 , z 1 = z j , F j y M , z M = z j + 1 , j M M 1
Now construct a function
μ j : J × R J j × R ,
μ j ( y , z ) = ( L j ( y ) , F j ( y , z ) ) , j M M 1 .
Theorem 1
(see [8,10]). The iterated function system  J × K ; μ j , j M M 1 possesses a unique attractor H , which is also the graph of the continuous function g : J R that satisfies g ( y j ) = z j j .
The above function g is known as the fractal interpolation function corresponding to the IFS J × K ; μ j , j M M 1 . The most widely studied fractal interpolations and their applications are taken from the following iterated function system:
J × K ; μ j y , z = L j y , F j y , z ; j M M 1
L j y = a j y + b j , F j y , z = δ j z + Q j y , where Q j : J R is a continuous function. The multiplier δ j is known as the vertical scaling factor of a continuous map μ j .
Let C r J be the collection of all r -times continuously differentiable real-valued functions on J . The differentiable fractal interpolation functions are given by the following theorem:
Theorem 2
[12]. Let ( y m , z m ) J × K : m = 1 , 2 , , M be the dataset, where y 1 < y 2 < < y M . Let L j y = a j y + b j , F j y , z = δ j z + Q j y , where Q j is a continuous function for j M M 1 . Suppose that for some integer r q 0 , δ j < a j q , j M M 1 . Let F j , r y , z = δ j z + Q j r y a j r , Q j r y be the rth derivative of Q j y
z 1 , r = Q 1 r y 1 a 1 r δ 1 , z M , k = Q M 1 r y M a M 1 r δ M 1 , r = 1 , 2 , , q .
If F j 1 , r y M , z M , r = F j , r y 1 , z 1 , r for j M M 2   &   r = 1 , 2 , , q then a fractal interpolation function g is determined by an iterated function system J × K ; μ j y , z = L j y , F j y , z , j M M 1 and the fractal interpolation function g r C q y 1 , y M is determined by an iterated function system J × K , μ j y , z = L j y , F j , r y , z , j M M 1 .

2.2. Zipper Fractal Interpolation Function

In this subsection, we review the construction of the zipper IFS for a given data set. First, we recall the definition of a zipper.
Definition 1.
Consider a non-surjective map   μ j , j M M 1  on a complete metric space  X , d . Then, the system  = X ; μ j : j M M 1  is known as the zipper-iterated function system with the vertices  v 1 , v 2 , , v M  and a signature  ε = ε 1 , ε 2 , , ε m 0 , 1 M 1 ,  if  μ j  satisfies  μ j v 1 = v j + ε j , μ j v M = v j + 1 ε j   j M M 1 . Let  B X  be any compact set satisfying the self-referential equation 
B = j = 1 M 1 μ j B
This is known as the attractor corresponding to the zipper IFS  . The definition of IFS with vertices  v 1 , v 2 , , v M  and the definition of zipper IFS coincide if we take the signature vector  ε j = 0 j M M 1 . The following review of the zipper fractal interpolation function is based on [10,46,50].
Let  y m , z m J × R , m M M  be the dataset with  y 1 < y 2 < < y M . Consider  J = [ y 1 , y M ]   &   J j = [ y j , y j + 1 ] , j M M 1 . Let  L j : J J j  the affine map given by  L j ( y ) = a j y + b j , j M M 1  such that
L j ( y 1 ) = y j + ε j , L j ( y M ) = y j + 1 ε j
a j = y j + 1 ε j y j + ε j y M y 1 , b j = y M y j + ε j y 1 y j + 1 ε j y M y 1 .
Let  j M M 1 ,   F j : K : = J × c , d R  be a function defined as 
F j ( y , z ) = δ j z + Q j ( y )
where  Q j : J R  is the continuous function and  δ j  be the vertical scaling factor such that  δ = max δ j : j M M 1 < 1 , and the following condition holds: 
F j ( y 1 , z 1 ) = z j + ε j , F j ( y M , z M ) = z j + 1 ε j , j M M 1 .
Now, define a map μ j : K J j × R , j M M 1 as
μ j y , z = L j y , F j y , z y , z K .
The system = K ; μ j , j M M 1 is a zipper-iterated function system with vertices y m , z m , m M M and a signature ε = ε 1 , ε 2 , , ε M 1 . The construction of Zipper-FIF is presented in the following theorem from [51].
Theorem 3.
The following conditions hold for the IFS  K ; μ j , j M M 1
  • There is a unique compact set  H K  such that  H = j = 1 M 1 μ j H .
  • The continuous function  g ε δ : J R  that interpolates the dataset  { ( y m , z m ) J × R : m M M } , i.e.,  H = y , g ε δ y : y J   &   g ε δ y j = z j j .

3. Construction of Rational Cubic Fractal Interpolation Function and Zipper-Rational Cubic Fractal Interpolation Function

This section is divided into two parts. In the first subsection, we discuss the construction of the rational cubic fractal interpolation function, and in the second, we discuss the construction of the zipper-rational cubic fractal interpolation function:

3.1. C1—Rational Cubic Fractal Interpolation Function

In this subsection, we provide the construction of the rational cubic fractal interpolation function with one shape parameter. Let y m , z m J × K ; m M M be the given set of data for an original function φ such that y 1 < y 2 < < y M . Let J × K ; μ j y , z = L j y , F j y , z j M M 1 be the IFS, where L j y = a j y + b j ,   F j y , z = δ j z + Q j y , Q j y = u j y v j y , u j y is a cubic polynomial and v j y is a linear polynomial, v j y 0 y J , δ j < a j , j M M 1 . Let
F j 1 y , d = δ j d + Q j 1 y a j ,
where Q j 1 y is the first-order derivative of Q j y , y y 1 , y M . The following join-up conditions are satisfied by F j y , S δ
F j y 1 , S 1 δ = z j , F j y M , S M δ = z j + 1 , F j 1 y 1 , d 1 = d j , F j 1 y M , d M = d j + 1 ,
where d j is the first-order derivative of S δ regarding y at knot y j . The attractor obtained by the above-iterated function system will be the graph of the C 1 -rational cubic fractal interpolation function. Now, from (3), it can be observed that the FIF can be written as:
S δ L j y = δ j S δ y + Q j y = δ j S δ y + u j θ v j θ ,
u j θ = A 1 θ 3 + B 1 θ 2 θ + C 1 θ θ 2 + D θ 3 ,
v j θ = 1 + s j 1 θ , θ = y y 1 l , l = y M y 1 y J ,
and s j is a positive shape parameter. The following conditions are imposed to ensure that the rational cubic fractal interpolation function is C 1 -continuous:
S δ L j y 1 = z j , S δ L j y M = z j + 1 , S δ L j y 1 = d j , S δ L j y M = d j + 1 , j M M 1 .
From Equations (7) and (8), it is evident that if we take y = y 1 in Equation (7), then we obtain
S δ L j y 1 = z j z j = δ j z 1 + A 1 + s j A = 1 + s j z j δ j z 1 .
Similarly, to calculate the value of D , if we take y = y M in Equation (7), then we obtain
S δ L j y M = z j + 1 z j + 1 = δ j z M + D 1 D = z j + 1 δ j z M .
Differentiate Equation (7) with respect to the y,
S δ 1 L j y L j   1 y = δ j S δ 1 y + Q j   1 y
Now, to calculate the value of B , if we take y = y 1 in S δ L j y and use Equation (8), we obtain
B = 3 + 2 s j z j δ j z 1 + l j 1 + s j a j d j δ j d 1
Similarly, to calculate the value of C , take y = y M in S δ L j y and use Equation (8). Then, we have
S δ L j y M = d j + 1 a j d j + 1 = δ j d M + C + 3 D D s j l j
C = 3 + s j z j + 1 δ j z M l j a j d j + 1 δ j d M
Now, by substituting the values of A , B , C   &   D into Equation (7), we obtain the required well-defined C 1 -rational cubic fractal interpolation function, whose numerator is:
u j θ = 1 + s j z j δ j z 1 1 θ 3 + 3 + 2 s j z j δ j z 1 + l j 1 + s j a j d j δ j d 1 1 θ 2 θ + 3 + s j z j + 1 δ j z M l j a j d j + 1 δ j d M 1 θ θ 2 + z j + 1 δ j z M θ 3 .
The derivatives d m , m M M can be calculated either using numerical methods or from the given set of data if it is not given. We calculated it using the arithmetic mean method on a given set of data.
Remark 1.
If  δ j = 0 , then the rational cubic fractal interpolation function will be reduced to the rational cubic interpolant, which is of the form
S y = u j * t v j * t , y y j , y j + 1 ,
with
u j * t = 1 + s j z j 1 t 3 + 3 + 2 s j z j + j 1 + s j d j 1 t 2 t + 3 + s j z j + 1 j d j + 1 1 t t 2 + z j + 1 t 3 , v j * t = 1 + s j 1 t , t = y y 1 y M y 1 , y J
Example 1.
To verify the diversity and flexibility of the rational cubic fractal interpolation function over the classical cubic fractal interpolation function, we provide an example considering the interpolating points  0 , 0 , 2 , 10 , 4 , 15 , 7 , 38 , 10 , 50 . Throughout Figure 1, the shape parameter  s j  is considered  50 , 50 , 50 , 50 . The rational cubic fractal interpolation function is shown in Figure 1. Now, if we compare the figures with each other, the change is visible and is due to different values of the scaling factor. We have taken different values for each figure. In Figure 1a, we have taken the scaling factor  δ = 0 , 0 , 0 , 0  and  δ = 0.2 , 0.25 , 0.45 , 0.5  in Figure 1b, comparing Figure 1a,b, Figure 1a is less detailed, whereas Figure 1b is more detailed and dense due to the increased scaling factor. Similarly, in Figure 1c,d, the scaling factors are  δ = 0.3 , 0.24 , 0.6 , 0.72  and  δ = 0.04 , 0.15 , 0.35 , 0.75 , respectively. Figure 1d is more basic compared to Figure 1c, where the shape is more filled in and intricate. Figure 1e with scaling factor  δ = ( 0.14 , 0.42 , 0.56 , 0.92 )  and Figure 1f with scaling factor  δ = 0.04 , 0.15 , 0.35 , 0.75  can be compared similarly. In this way, we can compare all these figures based on the scaling factor.

3.2. C1-Zipper-Rational Cubic Fractal Interpolation Function

In this subsection, we construct the zipper-rational cubic fractal interpolation function:
Consider the zipper-rational iterated function system, where the rational function is
Q j , ε y = u j , ε θ v j , ε θ
where
u j , ε θ = A j , ε 1 θ 3 + B j , ε 1 θ 2 θ + C j , ε 1 θ θ 2 + D j , ε θ 3 ,
v j , ε θ = 1 + s j 1 θ
For 0 θ = y y 1 y M y 1 1 , s j is a non-negative shape parameter. Suppose C ˜ J = g : J R , g y 1 = z 1 , g y M = z M . Then, C ˜ J is the subspace of space C J and C ˜ J is complete with respect to the uniform norm. Let ε j 0 , 1 M 1 be fixed and define the Read–Bajraktarevic (R–B) operator T ε δ : C ˜ J C ˜ J as:
T ε δ L j y = F j y , g y
= δ j g y + A j , ε 1 θ 3 + B j , ε 1 θ 2 θ + C j , ε 1 θ θ 2 + D j , ε θ 3 1 + s j 1 θ
T ε δ is a contraction on complete metric space C ˜ J , . . By the Banach contraction, theorem, T ε δ possesses a fixed point S ε δ that satisfies:
S ε δ L j y = δ j S ε δ y + Q j , ε y
= δ j S ε δ y + A j , ε 1 θ 3 + B j , ε 1 θ 2 θ + C j , ε 1 θ θ 2 + D j , ε θ 3 1 + s j 1 θ
The coefficients in this expression Q j , ε y can be calculated using conditions given in Equation (2), which is an operation similar to calculating them from the relation S ε δ y j = y j + ε j   &   S ε δ y j + 1 = y j + 1 ε j . According to [12], for a differentiable fractal interpolation function, the scaling factor is selected as δ j < a j j M M 1 . For the construction of the C1 zipper-rational cubic fractal interpolation function, we need the following interpolating conditions:
S ε δ L j y 1 = z j + ε j , S ε δ L j y M = z j + 1 ε j & S ε δ 1 L j y 0 = d j + ε j , S ε δ 1 L j y M = d j + 1 ε j ,
where d j + ε j is the first-order derivative with respect to “y” at knots y j , respectively. To calculate the value of A j , ε , take y = y 1 in Equation (10) and we obtain
A j , ε = 1 + s j z j + ε j δ j z 1
Similarly, to calculate the value of D j , ε , take y = y M in Equation (10) and we obtain
D j , ε = z j + 1 ε j δ j z M
By taking the first-order derivative of the functional equation, we obtain
S ε δ 1 L j y L j   1 y = δ j S ε δ 1 y + Q j , ε   1 y
to calculate the value of B j , ε , take y = y 1 in the above equation, and we obtain
B j , ε = 3 + 2 s j z j + ε j + 1 + s j h j d j + ε j δ j 3 + 2 s j z 1 + 1 + s j y M y 1 d 1
where h j = y j + 1 ε j y j + ε j . Again, to calculate the value of C j , ε , taking y = y M in the above equation, we obtain
C j , ε = 3 + s j z j + 1 ε j h j d j + 1 ε j δ j 3 + s j z M y M y 1 d M .
Thus, we obtain the desired zipper-rational cubic fractal interpolation function, which is
S ε δ L j y = δ j S ε δ y + u j , ε θ v j , ε θ
The values of u j , ε θ and v j , ε θ are provided below:
u j , ε θ u j , ε θ = 1 + s j z j + ε j δ j z 1 1 θ 3 + z j + 1 ε j δ j z M θ 3 + 3 + 2 s j z j + ε j + 1 + s j h j d j + ε j δ j 3 + 2 s j z 1 + 1 + s j y M y 1 d 1 1 θ 2 θ + 3 + s j z j + 1 ε j h j d j + 1 ε j δ j 3 + s j z M y M y 1 d M 1 θ θ 2 ,
v j , ε θ = 1 + s j 1 θ
Remark 2.
If the signature is  ε m = 0 , then the zipper-rational cubic fractal interpolation function reduces to the rational cubic fractal interpolation function
S δ L j y = δ j S δ y + Q j y , Q j y = u j θ v j θ
where the values of   u j θ  and  v j θ  are given in Section 3.1.
Remark 3.
If  ε j = 0 = δ j , then the zipper-rational cubic fractal interpolation function reduces to the rational cubic interpolant, discussed in Remark 1 of Section 3.1.
Example 2.
To verify the diversity and flexibility of the zipper-rational cubic fractal interpolation function compared to the rational cubic fractal interpolation function, we provide an example considering the interpolating points  0 , 0 , 2 , 10 , 4 , 15 , 7 , 38 , 10 , 50 . Throughout Figure 2, the shape parameter  s j  is considered  50 , 50 , 50 , 50 . In Figure 2a, we have taken the scaling factor  δ = 0 , 0 , 0 , 0  and signature as  ε = 0 , 0 , 0 , 0 , and in Figure 2b, the signature as  ε = 0 , 0 , 0 , 0  and scaling factor as  δ = 0.2 , 0.25 , 0.45 , 0.5 . Comparing Figure 2a,b, Figure 2a is less detailed, whereas Figure 2b is more is more detailed and dense due to the increased scaling factor. Since we have taken the signature as  ε = 0 , 0 , 0 , 0 , these figures are the same as Figure 1a and Figure 1b. Comparing Figure 2b,c, we have taken the same scaling factor as  δ = 0.2 , 0.25 , 0.45 , 0.5  for both figures and the signature as  ε = 0 , 0 , 0 , 0  and  ε = 1 , 1 , 1 , 0 , respectively, Figure 2c is more detailed or enhanced as it is producing a negative graph due to the change in signature value. Figure 2d,e can be compared as we have taken the same scaling factor  δ = 0.04 , 0.15 , 0.35 , 0.75  in both figures and signature as  ε = 0 , 1 , 1 , 0  and  ε = 1 , 1 , 1 , 1 , respectively. Comparing Figure 2e with scaling factor  δ = 0.04 , 0.15 , 0.35 , 0.75  and signature as  ε = 1 , 1 , 1 , 1  to Figure 2f with scaling factor  δ = ( 0.14 , 0.42 , 0.56 , 0.92 )  and signature as  ε = 1 , 1 , 1 , 1 , Figure 2f is more detailed compared to Figure 2e, which is due to the increase in the scaling factor.
The changes in the graphs of Figure 1 and Figure 2 are clearly visible. We have taken the same scaling factor in some of the graphs, but in Figure 2, due to the presence of signature, the curve of function is changed.

4. Convergence of Rational Cubic Fractal Interpolation Function and Zipper-Rational Cubic Fractal Interpolation Function

In this section, we determine the error bounds for uniform distance between the data-generating function φ , which belongs to C 2 y 1 , y M , and the rational cubic fractal interpolation function. It is difficult to calculate the uniform error bound φ S δ using the standard techniques since the rational cubic fractal interpolation function S δ has an implicit expression. Hence, we use the classical rational cubic function to drive the upper bound of the uniform error:
φ S δ φ S + S S δ
To show the convergence of the C1-zipper-rational cubic fractal interpolation function S ε δ toward the φ , the data-generating function belongs to C 2 y 1 , y M . We need to find a uniform distance between them. It will be difficult to calculate φ S ε δ using the standard techniques due to the implicit nature of S ε δ . Let S ε δ be the C1 zipper-rational cubic fractal interpolation function and φ data-generating function, respectively. We will derive an upper bound of the error using S ε δ S δ .
By the triangular inequality, we have
φ S ε δ φ S δ + S ε δ S δ
Now, we use the convergence result of the classical rational cubic interpolation function S .
Theorem 4.
The error between the data-generating function  φ C 2  and the classical rational cubic function is 
φ y S y 1 2 φ 2 j 2 e j , y [ y j , y j + 1 ] ,
where  e j = max 0 t j 1 σ s j , t s j 0  and for any positive value of  s j , the error optimal constant  e j  is bounded by  0 e j 0.2685 .
Theorem 5.
Let  φ  be the data-generating function for the given set of data  y m , z m J × K , m M M  and  S δ  be the  C 1 -rational cubic fractal interpolation function. At knot  y j ,  d j , j M M 1  is the bounded first-order derivative. Let
= max 1 j M 1 j , Ξ = φ + j t 1 t Ξ 1 , Ξ * = Z + j t 1 t Ξ 2
Ξ 1 = max d m 1 m M , Z = max z 1 , z M , Ξ 2 = max d 1 , d M
then, the condition holds
φ S δ e j 2 φ 2 + δ 1 δ Ξ + Ξ
where  e = max 1 j M 1 e j , e j = max 0 t 1 σ s j , t .
Proof. 
Consider  C ˜ J  as a space
C ˜ J = g : J R , | g y 1 = z 1 , g y M = z M , g y 1 = d 1 , g y M = d M
from Equations (2) and (7), the R–B operator
T ε δ : C ˜ J C ˜ J
for rational cubic fractal interpolation function, this can be written as
T δ * g y = δ j g L j 1 y + u j L j 1 y , δ j v j L j 1 y , y J j , j M M 1 .
Clearly, the interpolant S δ is a fixed point of T δ * with δ 0 and g being a fixed point of T δ * with δ = 0 . Since the R–B operator T δ * is a contractive operator with δ as its contractivity factor, we have
T δ * S δ T δ * S δ S δ S .
Using the mean value theorem of a function with multiple variables and from Equation (15), we have
T δ * S y T 0 * S y δ S + u j L j 1 y , τ j v j L j 1 y δ j ; τ j 0 , δ j .
Now, we calculate the right-hand side of Equation (17). From the classical rational cubic interpolation function in Equation (9), we have
S y = ω 1 s j , t z j + ω 2 s j , t z j + 1 + ω 3 s j , t d j ω 4 s j , t d j + 1 ,
where
ω 1 s j , t = 1 + s j 1 t 3 + 3 + 2 s j 1 t 2 t 1 + s j 1 t 0
ω 2 s j , t = t 3 + 3 + s j 1 t t 2 1 + s j 1 t 0
ω 3 s j , t = j 1 + s j 1 t 2 t 1 + s j 1 t 0   and ω 4 s j , t = j 1 t t 2 1 + s j 1 t 0
we can easily verify that ω 1 s j , t + ω 2 s j , t = 1 and ω 3 s j , t ω 4 s j , t = h j t 1 t . From Equation (18), we obtain
S y max m = j , j + 1 z m + h j t 1 t max m = j , j + 1 d m .
The above inequality is true j M M 1 ; we obtain the following expression:
S Ξ : = φ + j t 1 t Ξ 1 .
As v j t is independent of δ j , from the first term of the right-hand side of Equation (17),
u j L j 1 y , τ j v j L j 1 y δ j = ω 1 s j , t z 1 + ω 2 s j , t z M + ω 3 s j , t d 1 ω 4 s j , t d M .
Now, using the same arguments we used for the first part, we obtain the following bound:
u j L j 1 y , τ j v j L j 1 y δ j Ξ * ; = z + j t 1 t Ξ 2
substituting Equations (19) and (20) in Equation (17), we have
T δ * S y T 0 * S y δ Ξ + Ξ * , y y j , y j + 1
Since the above result is valid in each subinterval, then we obtain
T δ * S T 0 * S δ Ξ + Ξ * .
Using Equations (16) and (21) in
S δ S = T δ * S δ T 0 * S T δ * S δ T δ * S + T δ * S T 0 * S ,
from the above inequality, we have the following estimate:
S δ S δ Ξ + Ξ 1 δ .
As we have φ S δ φ S + S S δ using the Theorem 4 and Equation (22) in the above inequality, we obtain the desired upper bound
φ S δ e j 2 φ 2 + δ 1 δ Ξ + Ξ
Theorem 6.
Let  S ε δ  and  S δ  be the zipper-rational cubic fractal interpolation function and rational cubic fractal interpolation function, respectively, for the dataset  y m , z m J × K ; m M M  with vertical scaling  δ j < a j , and then for fixed signature  ε 0 , 1 M 1 ,
S ε δ S δ ε λ 1 δ 2 + s 4 ω φ , + ω φ ,
where
= max h j ; j M M 1 , ω φ , = sup y 1 y 2 h φ y 1 φ y 2 ; y 1 , y 2 J , λ = min λ j : j M M
λ j = min v j , ε θ : 0 θ 1   a n d   ε = 0 , i f ε j = 0 1 , i f ε j 0 .
Proof. 
Using the self-referential equation for the zipper-rational fractal interpolation function S ε δ and rational fractal interpolation function S δ , y J j , we have
S ε δ L j y S δ L j y = δ j S ε δ y S δ y + 1 + s j z j + ε j z j v j , ε θ 1 θ 3 + z j + 1 ε j z j + 1 v j , ε θ θ 3 + 3 + 2 s j z j + ε j + 1 + s j h j d j + ε j 3 + 2 s j z j 1 + s j h j d j v j , ε θ 1 θ 2 θ + 3 + s j z j + 1 ε j h j d j + 1 ε j 3 + s j z j + 1 + h j d j + 1 v j , ε θ 1 θ θ 2
using the properties of modulus, the above inequality can be written as
δ S ε δ y S δ y + 1 + s j z j + ε j z j v j , ε θ + 3 + 2 s j z j + ε j z j + 1 + s j h j d j + ε j d j v j , ε θ + 3 + s j z j + 1 ε j z j + 1 h j d j + 1 ε j d j + 1 v j , ε θ + z j + 1 ε j z j + 1 v j , ε θ
δ S ε δ S δ + ε j z j + 1 z j + s z j + 1 z j λ + ε j 3 z j + 1 z j + 2 s z j + 1 z j + d j + 1 d j + s d j + 1 d j λ + ε j 3 z j + 1 z j + s z j + 1 z j + d j + 1 d j λ + ε j z j + 1 z j λ
δ S ε δ S δ + ε 4 s z j + 1 z j + 8 z j + 1 z j + s d j + 1 d j + 2 d j + 1 d j λ
δ S ε δ S δ + ε 4 s ω φ , + 8 ω φ , + s ω φ , + 2 ω φ , λ
δ S ε δ S δ + ε 2 + s 4 ω φ , + ω φ , λ .
which is the required result. □
Theorem 7.
Let  S ε δ  be the zipper-rational cubic fractal interpolation function for the considered dataset  y m , z m J × K ; m M M  with fixed signature  ε 0 , 1 M  and the vertical scaling  δ j < a j , then
φ S δ e j 2 φ 2 + δ 1 δ Ξ + Ξ + ε λ 1 δ 2 + s 4 ω φ , + ω φ , .
Proof. 
Now, since in Theorems 5 and 6, we have calculated the values of φ S δ and S ε δ S δ , by substituting these values in the following inequality:
φ S ε δ φ S δ + S ε δ S δ
we obtain our desired result. □
Convergence Result: Since we have δ j a j , this implies δ h . From the above theorem, we can say that the zipper-rational cubic fractal interpolation function converges to the data-generating function φ as 0 , and in the case of the rational cubic fractal interpolation function δ l , then from Theorem 6, we can say that the rational cubic fractal interpolation function converges to the data-generating original function φ as 0 , if we take δ j a j 2 , then φ S ε δ = O 2 and φ S δ = O 2 as 0 .

5. Constrained Rational Cubic Fractal Interpolation Function and Zipper-Rational Cubic Fractal Interpolation Function

In this section, we discuss the constrained RCFIF and ZRCFIF, whose graph lies above the straight line. Due to the arbitrary choice of IFS parameters, the RCFIF and ZRCFIF may not lie above the straight line. Therefore, to avoid such circumstances, we have deduced the sufficient conditions on the shape parameter s j and scaling factor δ j .
Theorem 8.
Let  { ( y m , z m ) J × R : m M M }  be the given dataset and  S δ  be the rational cubic fractal interpolation function. Suppose that the given dataset lies above the straight line  L ¯ = Υ = n y + k , the parametric form of  L ¯  is  Υ L j y = 1 θ β + θ γ  on  y 1 , y M , then the C1 rational cubic fractal interpolation function defined here will lie above the straight line if the conditions given below hold:
δ j 0 , ς 0 , a , ς < a , a < ς
where
ς = min a j , z j β j z 1 , z j β j γ j z 1 , z j γ j z M
and conditions over shape parameters are
c o n d 1 < s j < c o n d 2
c o n d 1 = 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1 2 z j δ j z 1 β j γ j
c o n d 2 = max 0 , 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1 2 z j δ j z 1 β j γ j , z j + 1 δ j z M l j a j d j + 1 δ j d M z j + 1 δ j z M γ j .
Proof. 
Let { ( y m , z m ) J × R : m M M } be the given positive data which lie above a straight line. The parametric form of the straight line is Υ L j y = 1 θ β + θ γ and if we take y = y 1 then β j = Υ j = n y j + k . Similarly, if we take y = y M then β j + 1 = Υ j + 1 = n y j + 1 + k .
Now the C1 rational cubic fractal interpolation function lies above a straight line if
S δ L j y > Υ L j y δ j S δ y + u j θ v j θ > 1 θ β + θ γ
this implies
δ j S δ y + u j θ v j θ 1 θ β + θ γ > 0
δ j S δ y + u ^ j θ v j θ > 0
u ^ j θ = A ^ 1 θ 3 + B ^ 1 θ 2 θ + C ^ 1 θ θ 2 + D ^ θ 3
here,
A ^ = 1 + s j z j δ j z 1 β j & D ^ = z j + 1 δ j z M γ j
B ^ = s j 2 z j δ j z 1 β j γ j + 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1
C ^ = 2 + s j z j + 1 δ j z M γ j + z j δ j z M l j a j d j + 1 δ j d M
Now, to show that Equation (24) holds, we need to prove that Equation (25) is true for y = y j . As we know, S δ y j > 0 because of the scaling factor δ j 0 . The shape parameter s j is non-negative, provided v j θ > 0 . Now, all we are left with is u ^ j θ and u ^ j θ 0 if all the free parameters are non-negative. Therefore, let
A ^ 0 δ j z j β j z 1 and D ^ 0 δ j z j + 1 γ j z M
Let B ^ 0 . Then, we have
B ^ = s j 2 z j δ j z 1 β j γ j + 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1 if we take a j d j δ j d 1 0 , then s j 0 and 2 z j δ j z 1 β j γ j 0 δ j z j β j γ j z 1 , provided B ^ 0 . Otherwise, to validate that B ^ 0 , we choose
s j 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1 2 z j δ j z 1 β j γ j
If 2 z j δ j z 1 β j γ j > 0 , then
s j > 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1 2 z j δ j z 1 β j γ j
and if 2 z j δ j z 1 β j γ j < 0 , then
s j < 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1 2 z j δ j z 1 β j γ j
Similarly, let C ^ 0 . Then, we have
C ^ = 2 + s j z j + 1 δ j z M γ j + z j + 1 δ j z M l j a j d j + 1 δ j d M
If we take a j d j + 1 δ j d M 0 , then s j 0 and z j + 1 δ j z M γ j 0 δ j z j + 1 γ j z M , provided C ^ 0 . Otherwise, to validate that C ^ 0 , we choose
s j z j + 1 δ j z M l j a j d j + 1 δ j d M z j + 1 δ j z M γ j .
The above calculation can be written as
δ j 0 , ς 0 , a , ς < a , a < ς
where
ς = min a j , z j β j z 1 , z j β j γ j z 1 , z j γ j z M
and conditions over shape parameters are
c o n d 1 < s j < c o n d 2
c o n d 1 = 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1 2 z j δ j z 1 β j γ j
c o n d 2 = max 0 , 3 z j δ j z 1 2 β j γ j + 1 + s j l j a j d j δ j d 1 2 z j δ j z 1 β j γ j , z j + 1 δ j z M l j a j d j + 1 δ j d M z j + 1 δ j z M γ j .
Hence, we obtain the desired conditions for the rational fractal interpolation function to lie above the straight line.
Similarly, we discuss the sufficient condition for the C1-zipper-rational cubic fractal interpolation function to lie above the straight line.
Theorem 9.
Let  { ( y m , z m ) J × R : m M M }  be the given dataset and   S ε δ  be the zipper-rational cubic fractal interpolation function. Suppose that the given dataset lies above the straight line  L ¯ = Υ = n y + k , the parametric form of  L ¯  is  Υ L j y = 1 θ β + θ γ  on  y 1 , y M , then the C1 zipper-rational cubic fractal interpolation function defined here lies above the straight line if the conditions given below hold:
δ j 0 , ς 0 , a , ς < a , a < ς
where 
ς = min a j , z j + ε j β j z 1 , z j + ε j β j γ j z 1 , z j + 1 ε j γ j z M
and conditions over shape parameters are  
c o n d 1 < s j < c o n d 2
c o n d 1 = 3 z j + ε j δ j z 1 2 β j γ j + 1 + s j h j d j + ε j δ j y M y 1 d 1 2 z j + ε j δ j z 1 β j γ j
c o n d 2 = max 0 , 3 z j + ε j δ j z j 2 β j γ j + 1 + s j h j d j + ε j δ j y M y 1 d 1 2 z j + ε j δ j z 1 β j γ j , z j + 1 ε j δ j z M h j d j + 1 ε j + δ j y M y 1 d M z j + 1 ε j δ j z M γ j .
Proof. 
Let { ( y m , z m ) J × R : m M M } be the given positive data, which lies above the straight line. The parametric form of the straight line is Υ L j y = 1 θ β + θ γ and if we take y = y 1 then β j = Υ j = n y j + k . Similarly, if we take y = y M , then β j + 1 = Υ j + 1 = n y j + 1 + k . Now, the C1-zipper-rational cubic fractal interpolation function lies above a straight line if
S ε δ L j y > Υ L j y δ j S ε δ y + u j , ε θ v j , ε θ > 1 θ β + θ γ
δ j S ε δ y + u j , ε θ v j , ε θ 1 θ β + θ γ > 0
δ j S ε δ y + u j , ε θ v j , ε θ > 0
u j , ε θ = A ^ j , ε 1 θ 3 + B ^ j , ε 1 θ 2 θ + C ^ j , ε 1 θ θ 2 + D ^ j , ε θ 3
here,
A ^ j , ε = 1 + s j z j + ε j δ j z 1 β j , D ^ j , ε = z j + 1 ε j δ j z M γ j
B ^ j , ε = s j 2 z j + ε j δ j z 1 β j γ j + 3 z j + ε j δ j z 1 2 β j γ j
+ 1 + s j h j d j + ε j δ j y M y 1 d 1
C ^ j , ε = 2 + s j z j + 1 ε j δ j z M γ j + z j + 1 ε j δ j z M h j d j + 1 ε j + δ j y M y 1 d M
Now, to show that Equation (26) holds, we need to prove that Equation (27) is true for y = y j . As we know, S ε δ y j > 0 because of the scaling factor δ j 0 . In addition, the shape parameter s j is non-negative provided v j , ε θ > 0 . Now, all we are left with is u j , ε θ and u j , ε θ 0 if all the free parameters are non-negative. Therefore, let
A ^ j , ε 0 δ j z j + ε j β j z 1   and   D ^ j , ε 0 δ j z j + 1 ε j γ j z M
Let B ^ j , ε 0 , we have
B ^ j , ε = s j 2 z j + ε j δ j z 1 β j γ j + 3 z j + ε j δ j z 1 2 β j γ j + 1 + s j h j d j + ε j δ j y M y 1 d 1 if we take h j d j + ε j δ j y M y 1 d 1 0 , then s j 0 and 2 z j + ε j δ j z 1 β j γ j 0 δ j z j + ε j β j γ j z 1 , provided B ^ j , ε 0 . Otherwise, to validate that B ^ j , ε 0 , we choose
s j 3 z j + ε j δ j z 1 2 β j γ j + 1 + s j h j d j + ε j δ j y M y 1 d 1 2 z j + ε j δ j z 1 β j γ j
If 2 z j + ε j δ j z 1 β j γ j > 0 then
s j > 3 z j + ε j δ j z 1 2 β j γ j + 1 + s j h j d j + ε j δ j y M y 1 d 1 2 z j + ε j δ j z 1 β j γ j
and if 2 z j + ε j δ j z 1 β j γ j < 0 then
s j < 3 z j + ε j δ j z 1 2 β j γ j + 1 + s j h j d j + ε j δ j y M y 1 d 1 2 z j + ε j δ j z 1 β j γ j
Similarly, let C ^ j , ε 0 , we have
C ^ j , ε = 2 + s j z j + 1 ε j δ j z M γ j + z j + 1 ε j δ j z M h j d j + 1 ε j + δ j y M y 1 d M
if we take h j d j + 1 ε j + δ j y M y 1 d M 0 , then s j 0 and z j + 1 ε j δ j z M γ j 0 δ j z j + 1 ε j γ j z M , provided C ^ j , ε 0 . Otherwise, to validate that C ^ j , ε 0 , we choose
s j z j + 1 ε j δ j z M h j d j + 1 ε j + δ j y M y 1 d M z j + 1 ε j δ j z M γ j
Hence, we obtain the desired conditions for the rational fractal interpolation function to lie above the straight line. □

6. Positivity of Rational Cubic Fractal Interpolation Function and Zipper-Rational Cubic Fractal Interpolation Function

In this section, we discuss the conditions over the scaling factor and shape parameter to preserve the positivity conditions for the rational cubic fractal interpolation function and zipper-rational cubic fractal interpolation function. In shape preservation, we deal with three properties of functions, which are positivity, monotonicity, and convexity. Positivity is one of the critical properties of shape preservation. This is important in various real-world applications such as physical quantities like mass, length etc., probability distributions, time series analysis, and so on. We constrained the scaling factor and shape parameter to obtain sufficient conditions for positivity. The shape parameter plays a crucial role in the curvature of a function and in defining the behavior of a function in its domain. In Theorem 10, we discuss the positivity of the rational cubic fractal interpolation function, and in Theorem 11, we discuss the positivity of the zipper-rational cubic fractal interpolation function.
Theorem 10.
Let  y m , z m J × K ; m M M  be the given positive dataset, then the rational cubic fractal interpolation function  S δ  is positive if the following conditions hold for the shape parameter and scaling factor:
0 δ j < min a j , z j z 1 , z j + 1 z M
s j > max 0 , l j a j d j δ j d 1 z j δ j z 1 , l j a j d j + 1 δ j d M z j + 1 δ j z M , j M M 1
Proof. 
To prove that the rational cubic fractal interpolation function S δ is positive, it will be sufficient to show that S δ L j y 0 . By the hypothesis of the iterated function system, we can say that S δ L j y will be positive if u j θ 0   θ 0 , 1 . To prove that u j θ 0 , we need to show that A , B , C and D 0 , so we have A 0 1 + s j z j δ j z 1 0 , as s j is a non-negative shape parameter, the scaling factor is constrained as
δ j z j z 1
Similarly, D 0 z j + 1 δ j z M 0 . Then, we obtain
δ j z j + 1 z M
Now, B 0 if 3 + 2 s j z j δ j z 1 + l j 1 + s j a j d j δ j d 1 0
s j l j a j d j δ j d 1 z j δ j z 1 ; j M M 1 ,   provided   z j δ j z 1 > 0
Similarly, C 0 if 3 + s j z j + 1 δ j z M l j a j d j + 1 δ j d M 0
s j l j a j d j + 1 δ j d M z j + 1 δ j z M ; j M M 1 ,   provided   z j + 1 δ j z M > 0 .
Example 3.
Let  0 , 0 , 1 , 5 , 2 , 7 , 3 , 9 , 4 , 10  be the positive interpolating points, and with these data points, we generate the positive rational cubic fractal interpolation function with different values of shape parameters and scaling factors. Now, for the positive rational cubic fractal interpolation function, we have calculated the shape parameters according to the above theorem. In Figure 3a, we have taken the shape parameter, i.e.,  s j = 2 , 4 , 4.2 , 3  and the scaling factor  δ = 0.15 , 0.5 , 0.35 , 0.47 . The corresponding graph is positive. In Figure 3b and Figure 3c, we used the same shape parameter, i.e.,   s j = 7 , 9.8 , 10.2 , 11  but different scaling factors   δ = 0.24 , 0.12 , 0.32 , 0.55  and  δ = 0.15 , 0.25 , 0.32 , 0.46 , respectively. We can easily witness the change in shape due to the different vertical scaling factors. In Figure 3d and Figure 3e, we used the same scaling factor, i.e.,   δ = 0.1 , 0.12 , 0.23 , 0.5  but different shape parameters, which are  s j = 5 , 5.8 , 6.2 , 7  and   s j = 7 , 9.8 , 10.2 , 11 , respectively. The change in the graph of the positive rational cubic fractal interpolation function is due to the different values for signature and shape parameters. Similarly, by taking the shape parameter   s j = 9 , 11.2 , 12.8 , 13  and scaling factor,   δ = 0.15 , 0.25 , 0.35 , 0.45 , Figure 3f is generated, which is another positive graph for the rational cubic fractal interpolation function.
Theorem 11.
Let  y m , z m J × K ; m M M  be the considered positive dataset, then the zipper-rational cubic fractal interpolation function  S ε δ  is positive if the following conditions hold for the shape parameter and scaling factor:
0 δ j < min a j , z j + ε j z 1 , z j + 1 ε j z M
s j > max 0 , h j d j + ε j δ j y M y 1 d 1 z j + ε j δ j z 1 , h j d j + 1 ε j δ j y M y 1 d M z j + 1 ε j δ j z M , j M M 1
Proof. 
To prove that the zipper-rational cubic fractal interpolation function S ε δ is positive, it will be sufficient to show that S ε δ L j y 0 , and by the hypothesis of the iterated function system, we can say that S ε δ L j y will be positive if u j , ε θ 0 θ 0 , 1 . To prove that u j , ε θ 0 we need to show that A j , ε , B j , ε , C j , ε and D j , ε 0 , so we have A j , ε 0 1 + s j z j + ε j δ j z 1 0 , as s j is a non-negative shape parameter, so
δ j z j + ε j z 1
Similarly, D j , ε 0 z j + 1 ε j δ j z M 0 , we obtain
δ j z j + 1 ε j z M
Now, B j , ε 0 if 3 + 2 s j z j + ε j + 1 + s j h j d j + ε j δ j 3 + 2 s j z 1 + 1 + s j y M y 1 d 1 0
s j h j d j + ε j δ j y M y 1 d 1 z j + ε j δ j z 1 ; j M M 1 ,   provided   z j + ε j δ j z 1 > 0
Similarly, C j , ε 0 if 3 + s j z j + 1 ε j h j d j + 1 ε j δ j 3 + s j z M y M y 1 d M 0
s j h j d j + 1 ε j δ j y M y 1 d M z j + 1 ε j δ j z M ; j M M 1 ,   provided   z j + 1 ε j δ j z M > 0 .
Example 4.
Let  0 , 0 , 1 , 5 , 2 , 7 , 3 , 9 , 4 , 10  be the positive interpolating points, and with these data points, we will be generating the positive zipper-rational cubic fractal interpolation function with different shape parameters, signatures, and scaling factors. In Figure 4a, we used the arbitrary shape parameters, i.e.,   s j = 2 , 4 , 4.2 , 3   and the scaling factor and signature are  δ = 0.15 , 0.5 , 0.35 , 0.47  , and   ε = 1 , 1 , 1 , 1 , respectively. The corresponding graph is non-positive. For the positive zipper-rational cubic fractal interpolation function, we calculated the shape parameters according to the above theorem. In Figure 4b and Figure 4c, we used the same shape parameter, i.e.,   s j = 7 , 9.8 , 10.2 , 11  but different scaling factor and signature, i.e.,  δ = 0.24 , 0.12 , 0.32 , 0.55 ,   δ = 0.15 , 0.25 , 0.32 , 0.46 , and  ε = 0 , 0 , 0 , 0 ,   ε = 1 , 1 , 1 , 1 , respectively. We can clearly witness the change in shape due to the signature and scaling factor. In Figure 4d and Figure 4e, we used the same scaling factor and signature, i.e.,   δ = 0.1 , 0.12 , 0.23 , 0.5 , ε = ( 0 , 1 , 0 , 1 )   but different shape parameters, which are   s j = 5 , 5.8 , 6.2 , 7  and  s j = 7 , 9.8 , 10.2 , 11 , respectively. We can witness the change in the graph of the positive zipper-rational cubic fractal interpolation function due to the different values for signature and shape parameters.
If we see Figure 3a and Figure 4a, we have taken the same values for the scaling factor, but for the C1 rational cubic fractal interpolation function, it is a positive graph. For the C1 zipper-rational cubic fractal interpolation function, it is a negative graph.

7. Conclusions and Future Work

In this paper, the C1 rational cubic fractal interpolation function (RCFIF) and C1 zipper cubic fractal interpolation function (ZRCFIF) with a linear denominator and one shape parameter are constructed. The proposed ZRCFIF becomes the classical rational cubic interpolation function, if we take zero values of both signature and scaling factor. Also, if only the values of the scaling factor are taken to be zero, then ZRCFIF reduces to the classical C1 zipper-rational cubic interpolation function, and if only the signature is taken to be zero, then it reduces to RCFIF. The uniform convergence of the C1 rational cubic fractal interpolation function and the C1 zipper-rational cubic fractal interpolation function to the original function is also discussed. The results are also supported by numerical examples for both the rational cubic fractal interpolation function and the zipper-rational cubic fractal interpolation function. The graphs of the interpolating function are demonstrated for various selections of signature, scaling factor, and shape parameters. If one suitably chooses (as stated in Theorems 10 and 11) the range of scaling factors and shape parameters, then the proposed RQFIF and ZRQFIF preserve positivity.
Figure 1, Figure 2, Figure 3 and Figure 4 demonstrate the sensitivity of the proposed class of fractal interpolation functions with respect to the signature, scaling factors, and shape parameters. We can also see that the zipper fractal interpolation function is more flexible than the fractal interpolation function due to the presence of a signature. Duan et al. [52] demonstrated that cubic spline gives a better approximation with a linear denominator than with quadratic or cubic denominators. In that way, the zipper-rational cubic fractal interpolation function with a linear denominator will give a better approximation as compared to the zipper fractal interpolation function with quadratic or cubic denominators. The proposed C1 RCFIFs and ZRCFIFs can be used for data visualization and in physics.
For future work, the vertical scaling factor can be replaced by a variable scaling for generating a C1 rational cubic fractal interpolation function and a C1 zipper-rational cubic fractal interpolation function with variable scaling.

Author Contributions

S.S.: Conceptualization, Writing—Original Draft, K.K.: Writing—Review and Editing, Supervision, G.S.: Methodology, Investigation, M.K.W.: Writing—Review and Editing, and R.S.: Supervision, Visualization, Writing—Review and Editing, Project Administration, Investigation, Formal Analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used in the current study is available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this research.

References

  1. Schmidt, J.W.; Heß, W. Positivity of cubic polynomials on intervals and positive spline interpolation. BIT Numer. Math. 1988, 28, 340–352. [Google Scholar] [CrossRef]
  2. Butt, S.; Brodlie, K.W. Preserving positivity using piecewise cubic interpolation. Comput. Graph. 1993, 17, 55–64. [Google Scholar] [CrossRef]
  3. Sarfraz, M.; Hussain, M.Z. Data visualization using rational spline interpolation. J. Comput. Appl. Math. 2006, 189, 513–525. [Google Scholar] [CrossRef]
  4. Hussain, M.Z.; Sarfraz, M. Positivity-preserving interpolation of positive data by rational cubics. J. Comput. Appl. Math. 2008, 218, 446–458. [Google Scholar] [CrossRef]
  5. Chernous’ ko, F.L.; Ananievski, I.M.; Reshmin, S.A. Control of Nonlinear Dynamical Systems: Methods and Applications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  6. Jayaraman, A.; Belmonte, A. Oscillations of a solid sphere falling through a wormlike micellar fluid. Phys. Rev. E. 2003, 67, 065301. [Google Scholar] [CrossRef]
  7. Massopust, P. Interpolation and Approximation with Splines and Fractals; Oxford University Press, Inc.: Oxford, UK, 2010. [Google Scholar]
  8. Barnsley, M.F. Fractals Everywhere; Academic Press: Cambridge, MA, USA, 2014. [Google Scholar]
  9. Gibert, S.; Massopust, P.R. The exact Hausdorff dimension for a class of fractal functions. J. Math. Anal. Appl. 1992, 168, 171–183. [Google Scholar] [CrossRef]
  10. Barnsley, M.F. Fractal functions and interpolation. Constr. Approx. 1986, 2, 303–329. [Google Scholar] [CrossRef]
  11. Hutchinson, J.E. Fractals and self similarity. Indiana Univ. Math. J. 1981, 30, 713–747. [Google Scholar] [CrossRef]
  12. Barnsley, M.F.; Harrington, A.N. The calculus of fractal interpolation functions. J. Approx. Theory 1989, 57, 14–34. [Google Scholar] [CrossRef]
  13. Navascués, M.A.; Sebastián, M.V. Smooth fractal interpolation. J. Inequalities Appl. 2006, 2006, 78734. [Google Scholar] [CrossRef]
  14. Prasad, B.; Katiyar, K.; Singh, B. Trigonometric quadratic fractal interpolation functions. Int. J. Appl. Eng. Res. 2015, 37, 290–302. [Google Scholar]
  15. Garg, S.; Katiyar, K. A new type of zipper fractal interpolation surfaces and associated bivariate zipper fractal operator. J. Anal. 2023, 31, 3021–3043. [Google Scholar] [CrossRef]
  16. Gautam, S.S.; Katiyar, K. Alpha fractal rational quintic spline with shape preserving properties. Int. J. Comput. Appl. Math. Comput. Sci. 2023, 3, 113–121. [Google Scholar] [CrossRef]
  17. Chand, A.K.B.; Vijender, N.; Navascués, M.A. Shape preservation of scientific data through rational fractal splines. Calcolo 2014, 51, 329–362. [Google Scholar] [CrossRef]
  18. Reddy, K.M.; Chand, A.K.; Viswanathan, P. Data visualization by rational fractal function based on function values. J. Anal. 2020, 28, 261–277. [Google Scholar] [CrossRef]
  19. Navascués, M.A. Fractal polynomial interpolation. Z. Anal. Ihre Anwendungen 2005, 24, 401–418. [Google Scholar] [CrossRef]
  20. Banerjee, A.; Akhtar, M.N.; Navascués, M.A. Local α-fractal interpolation function. Eur. Phys. J. Spec. Top. 2023, 232, 1043–1050. [Google Scholar] [CrossRef]
  21. Feng, Z.; Sun, X. Box-counting dimensions of fractal interpolation surfaces derived from fractal interpolation functions. J. Math. Anal. Appl. 2014, 412, 416–425. [Google Scholar] [CrossRef]
  22. Akhtar, M.N.; Prasad, M.G.; Navascués, M.A. Box dimensions of α-fractal functions. Fractals 2016, 24, 1650037. [Google Scholar] [CrossRef]
  23. Nayak, S.R.; Mishra, J.; Palai, G. Analysing roughness of surface through fractal dimension: A review. Image Vis. Comput. 2019, 89, 21–34. [Google Scholar] [CrossRef]
  24. Ri, S. New types of fractal interpolation surfaces. Chaos Solitons Fractals 2019, 119, 291–307. [Google Scholar] [CrossRef]
  25. Navascués, M.A.; Akhtar, M.N.; Mohapatra, R. Fractal frames of functions on the rectangle. Fractal Fract. 2021, 5, 42. [Google Scholar] [CrossRef]
  26. Pandey, M.; Som, T.; Verma, S. Fractal dimension of Katugampola fractional integral of vector-valued functions. Eur. Phys. J. Spec. Top. 2021, 230, 3807–3814. [Google Scholar] [CrossRef]
  27. Navascués, M.A. Fractal curves on Banach algebras. Fractal Fract. 2022, 6, 722. [Google Scholar] [CrossRef]
  28. He, C.H.; Liu, C. Fractal dimensions of a porous concrete and its effect on the concrete’s strength. Facta Univ. Ser. Mech. Eng. 2023, 21, 137–150. [Google Scholar] [CrossRef]
  29. Lal, R.; Selmi, B.; Verma, S. On Dimension of Fractal Functions on Product of the Sierpiński Gaskets and Associated Measures. Results Math. 2024, 79, 73. [Google Scholar] [CrossRef]
  30. Navascués, M.A. Approximation of fixed points and fractal functions by means of different iterative algorithms. Chaos Solitons Fractals 2024, 180, 114535. [Google Scholar] [CrossRef]
  31. Chandra, S.; Verma, S.; Abbas, S. Construction of fractal functions using Kannan mappings and smoothness analysis. arXiv 2023, arXiv:2301.03075. [Google Scholar]
  32. Thangaraj, C.; Easwaramoorthy, D.; Selmi, B.; Chamola, B.P. Generation of fractals via iterated function system of Kannan contractions in controlled metric space. Math. Comput. Simul. 2024, 222, 188–198. [Google Scholar] [CrossRef]
  33. Garg, S.; Katiyar, K. Positivity and monotonicity shape preserving using rational quintic fractal interpolation functions. Adv. Math. Sci. J. 2020, 9, 5511–5520. [Google Scholar]
  34. Tyada, K.R.; Chand, A.K.B.; Sajid, M. Shape preserving rational cubic trigonometric fractal interpolation functions. Math. Comput. Simulation. 2021, 90, 866–891. [Google Scholar] [CrossRef]
  35. Ibraheem, F.; Hussain, M.; Hussain, M.Z.; Bhatti, A.A. Positive data visualization using trigonometric function. J. Appl. Math. 2012, 2012, 247120. [Google Scholar] [CrossRef]
  36. Dube, M.; Rana, P.S. Positivity preserving interpolation of positive data by rational quadratic trigonometric spline. IOSR J. Math. 2014, 10, 42–47. [Google Scholar] [CrossRef]
  37. Hussain, M.Z.; Ali, J. Positivity-preserving piecewise rational cubic interpolation. Matematika 2006, 22, 147–153. [Google Scholar]
  38. Piah, A.R.M.; Unsworth, K. Improved sufficient conditions for monotonic piecewise rational quartic interpolation. Sains Malays. 2011, 40, 1173–1178. [Google Scholar]
  39. Abbas, M.; Majid, A.A.; Ali, J.M. Monotonicity-preserving C2 rational cubic spline for monotone data. Appl. Math. Comput. 2012, 219, 2885–2895. [Google Scholar] [CrossRef]
  40. Chand, A.K.B.; Vijender, N. Monotonicity preserving rational quadratic fractal interpolation functions. Adv. Numer. Anal. 2014, 2014, 504825. [Google Scholar] [CrossRef]
  41. Dube, M.; Tiwari, P. Convexity preserving C2 rational quadratic trigonometric spline. Int. J. Sci. Res. Publ. 2012, 3, 401. [Google Scholar]
  42. Sharma, S.; Katiyar, K. Preserving convexity through C2 rational quintic fractal interpolation function. InAIP Conf. Proc. 2023, 2735, 040015. [Google Scholar]
  43. Brodlie, K.W.; Butt, S. Preserving convexity using piecewise cubic interpolation. Comput. Graph. 1991, 15, 15–23. [Google Scholar] [CrossRef]
  44. Viswanathan, P.; Chand, A.K.B.; Agarwal, R.P. Preserving convexity through rational cubic spline fractal interpolation function. J. Comput. Appl. Math. 2014, 263, 262–276. [Google Scholar] [CrossRef]
  45. Aseev, V.V.; Tetenov, A.V.; Kravchenko, A.S. On selfsimilar Jordan curves on the plane. Sib. Math. J. 2003, 44, 379–386. [Google Scholar] [CrossRef]
  46. Aseev, V.V.; Tetenov, A.V. On the self-similar Jordan arcs admitting structure parametrization. Sib. Math. J. 2005, 46, 581–592. [Google Scholar] [CrossRef]
  47. Jha, S.; Chand, A.K.B. Zipper rational quadratic fractal interpolation functions. In Proceedings of the Fifth International Conference on Mathematics and Computing: ICMC; Springer: Singapore, 2020; pp. 229–241. [Google Scholar]
  48. Chand, A.K.B. Zipper fractal functions with variable scalings. Adv. Theory Nonlinear Anal. Appl. 2022, 6, 481–501. [Google Scholar]
  49. Vijay; Chand, A.K.B. Convexity-preserving rational cubic zipper fractal interpolation curves and surfaces. Math. Comput. Appl. 2023, 28, 74. [Google Scholar]
  50. Chand, A.K.B.; Vijender, N.; Viswanathan, P.; Tetenov, A.V. Affine zipper fractal interpolation functions. BIT Numer. Math. 2020, 60, 319–344. [Google Scholar] [CrossRef]
  51. Aseev, V.V. On the regularity of self-similar zippers. In Proceedings of the 6-th Russian-Korean International Symposium on Science and Technology, KORUS-2002, Novosibirsk, Russia, 24–30 June 2002; p. 167. [Google Scholar]
  52. Duan, Q.; Djidjeli, K.; Price, W.G.; Twizell, E.H. The approximation properties of some rational cubic splines. Int. J. Comput. Math. 1999, 72, 155–166. [Google Scholar] [CrossRef]
Figure 1. Rational cubic fractal interpolation functions S δ with the change in scaling factor. (a) Rational cubic fractal interpolation with scaling factor δ = ( 0 , 0 , 0 , 0 ) . (b) Rational cubic fractal interpolation with scaling factor δ = 0.2 , 0.25 , 0.45 , 0.5 . (c) Rational cubic fractal interpolation with scaling factor δ = 0.3 , 0.24 , 0.6 , 0.72 . (d) Rational cubic fractal interpolation with scaling factor δ = 0.04 , 0.15 , 0.35 , 0.75 . (e) Rational cubic fractal interpolation with scaling factor δ = 0.04 , 0.15 , 0.35 , 0.75 . (f) Rational cubic fractal interpolation with scaling factor δ = 0.14 , 0.42 , 0.56 , 0.92 .
Figure 1. Rational cubic fractal interpolation functions S δ with the change in scaling factor. (a) Rational cubic fractal interpolation with scaling factor δ = ( 0 , 0 , 0 , 0 ) . (b) Rational cubic fractal interpolation with scaling factor δ = 0.2 , 0.25 , 0.45 , 0.5 . (c) Rational cubic fractal interpolation with scaling factor δ = 0.3 , 0.24 , 0.6 , 0.72 . (d) Rational cubic fractal interpolation with scaling factor δ = 0.04 , 0.15 , 0.35 , 0.75 . (e) Rational cubic fractal interpolation with scaling factor δ = 0.04 , 0.15 , 0.35 , 0.75 . (f) Rational cubic fractal interpolation with scaling factor δ = 0.14 , 0.42 , 0.56 , 0.92 .
Axioms 13 00584 g001
Figure 2. Zipper-rational cubic fractal interpolation functions S ε δ with the different values of signature ε and scaling factor δ . (a) zipper-rational cubic fractal interpolation with signature ε = ( 0 , 0 , 0 , 0 ) and scaling factor δ = ( 0 , 0 , 0 , 0 ) . (b) zipper-rational cubic fractal interpolation with signature ε = ( 0 , 0 , 0 , 0 ) and scaling factor δ = 0.2 , 0.25 , 0.45 , 0.5 . (c) zipper-rational cubic fractal interpolation with signature ε = ( 1 , 1 , 1 , 0 ) and scaling factor δ = 0.2 , 0.25 , 0.45 , 0.5 . (d) zipper-rational cubic fractal interpolation with signature ε = ( 0 , 1 , 1 , 0 ) and scaling factor δ = 0.04 , 0.15 , 0.35 , 0.75 . (e) zipper-rational cubic fractal interpolation with signature ε = ( 1 , 1 , 1 , 1 ) and scaling factor δ = 0.04 , 0.15 , 0.35 , 0.75 . (f) zipper-rational cubic fractal interpolation with signature ε = ( 1 , 1 , 1 , 1 ) and scaling factor δ = 0.14 , 0.42 , 0.56 , 0.92 .
Figure 2. Zipper-rational cubic fractal interpolation functions S ε δ with the different values of signature ε and scaling factor δ . (a) zipper-rational cubic fractal interpolation with signature ε = ( 0 , 0 , 0 , 0 ) and scaling factor δ = ( 0 , 0 , 0 , 0 ) . (b) zipper-rational cubic fractal interpolation with signature ε = ( 0 , 0 , 0 , 0 ) and scaling factor δ = 0.2 , 0.25 , 0.45 , 0.5 . (c) zipper-rational cubic fractal interpolation with signature ε = ( 1 , 1 , 1 , 0 ) and scaling factor δ = 0.2 , 0.25 , 0.45 , 0.5 . (d) zipper-rational cubic fractal interpolation with signature ε = ( 0 , 1 , 1 , 0 ) and scaling factor δ = 0.04 , 0.15 , 0.35 , 0.75 . (e) zipper-rational cubic fractal interpolation with signature ε = ( 1 , 1 , 1 , 1 ) and scaling factor δ = 0.04 , 0.15 , 0.35 , 0.75 . (f) zipper-rational cubic fractal interpolation with signature ε = ( 1 , 1 , 1 , 1 ) and scaling factor δ = 0.14 , 0.42 , 0.56 , 0.92 .
Axioms 13 00584 g002
Figure 3. Positive rational cubic fractal interpolation with different values of shape parameter and scaling factor S δ . (a) Positive rational cubic fractal interpolation with shape parameter s j = ( 2 , 4 , 4.2 , 3 ) and scaling factor δ = ( 0.15 , 0.5 , 0.35 , 0.47 ) . (b) Positive rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 and scaling factor δ = ( 0.24 , 0.12 , 0.32 , 0.55 ) . (c) Positive rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 and scaling factor δ = 0.15 , 0.25 , 0.32 , 0.46 . (d) Positive rational cubic fractal interpolation with shape parameter s j = 5 , 5.8 , 6.2 , 7 and scaling factor δ = 0.1 , 0.12 , 0.23 , 0.5 . (e) Positive rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 and scaling factor δ = 0.1 , 0.12 , 0.23 , 0.5 . (f) Positive rational cubic fractal interpolation with shape parameter s j = 9 , 11.2 , 12.8 , 13 and scaling factor δ = 0.15 , 0.25 , 0.35 , 0.45 .
Figure 3. Positive rational cubic fractal interpolation with different values of shape parameter and scaling factor S δ . (a) Positive rational cubic fractal interpolation with shape parameter s j = ( 2 , 4 , 4.2 , 3 ) and scaling factor δ = ( 0.15 , 0.5 , 0.35 , 0.47 ) . (b) Positive rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 and scaling factor δ = ( 0.24 , 0.12 , 0.32 , 0.55 ) . (c) Positive rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 and scaling factor δ = 0.15 , 0.25 , 0.32 , 0.46 . (d) Positive rational cubic fractal interpolation with shape parameter s j = 5 , 5.8 , 6.2 , 7 and scaling factor δ = 0.1 , 0.12 , 0.23 , 0.5 . (e) Positive rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 and scaling factor δ = 0.1 , 0.12 , 0.23 , 0.5 . (f) Positive rational cubic fractal interpolation with shape parameter s j = 9 , 11.2 , 12.8 , 13 and scaling factor δ = 0.15 , 0.25 , 0.35 , 0.45 .
Axioms 13 00584 g003aAxioms 13 00584 g003b
Figure 4. Non-positive and positive zipper-rational cubic fractal interpolation S ε δ with different values of the shape parameter s j , signature ε , and scaling factor δ . (a) Non-positive zipper-rational cubic fractal interpolation with shape parameter s j = ( 2 , 4 , 4.2 , 3 ) , signature ε = ( 1 , 1 , 1 , 1 ) , and scaling factor δ = ( 0.15 , 0.5 , 0.35 , 0.47 ) . (b) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 signature ε = ( 0 , 0 , 0 , 0 ) , and scaling factor δ = ( 0.24 , 0.12 , 0.32 , 0.55 ) . (c) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 , signature ε = ( 1 , 1 , 1 , 1 ) , and scaling factor δ = 0.15 , 0.25 , 0.32 , 0.46 . (d) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 5 , 5.8 , 6.2 , 7 , signature ε = ( 0 , 1 , 0 , 1 ) and scaling factor δ = 0.1 , 0.12 , 0.23 , 0.5 . (e) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 , signature ε = ( 0 , 1 , 0 , 1 ) and scaling factor δ = 0.1 , 0.12 , 0.23 , 0.5 . (f) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 9 , 11.2 , 12.8 , 13 , signature ε = ( 1 , 1 , 1 , 1 ) and scaling factor δ = 0.15 , 0.25 , 0.35 , 0.45 .
Figure 4. Non-positive and positive zipper-rational cubic fractal interpolation S ε δ with different values of the shape parameter s j , signature ε , and scaling factor δ . (a) Non-positive zipper-rational cubic fractal interpolation with shape parameter s j = ( 2 , 4 , 4.2 , 3 ) , signature ε = ( 1 , 1 , 1 , 1 ) , and scaling factor δ = ( 0.15 , 0.5 , 0.35 , 0.47 ) . (b) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 signature ε = ( 0 , 0 , 0 , 0 ) , and scaling factor δ = ( 0.24 , 0.12 , 0.32 , 0.55 ) . (c) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 , signature ε = ( 1 , 1 , 1 , 1 ) , and scaling factor δ = 0.15 , 0.25 , 0.32 , 0.46 . (d) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 5 , 5.8 , 6.2 , 7 , signature ε = ( 0 , 1 , 0 , 1 ) and scaling factor δ = 0.1 , 0.12 , 0.23 , 0.5 . (e) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 7 , 9.8 , 10.2 , 11 , signature ε = ( 0 , 1 , 0 , 1 ) and scaling factor δ = 0.1 , 0.12 , 0.23 , 0.5 . (f) Positive zipper-rational cubic fractal interpolation with shape parameter s j = 9 , 11.2 , 12.8 , 13 , signature ε = ( 1 , 1 , 1 , 1 ) and scaling factor δ = 0.15 , 0.25 , 0.35 , 0.45 .
Axioms 13 00584 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sharma, S.; Katiyar, K.; Sudhamsu, G.; Wratch, M.K.; Salgotra, R. Positivity-Preserving Rational Cubic Fractal Interpolation Function Together with Its Zipper Form. Axioms 2024, 13, 584. https://doi.org/10.3390/axioms13090584

AMA Style

Sharma S, Katiyar K, Sudhamsu G, Wratch MK, Salgotra R. Positivity-Preserving Rational Cubic Fractal Interpolation Function Together with Its Zipper Form. Axioms. 2024; 13(9):584. https://doi.org/10.3390/axioms13090584

Chicago/Turabian Style

Sharma, Shamli, Kuldip Katiyar, Gadug Sudhamsu, Manjinder Kaur Wratch, and Rohit Salgotra. 2024. "Positivity-Preserving Rational Cubic Fractal Interpolation Function Together with Its Zipper Form" Axioms 13, no. 9: 584. https://doi.org/10.3390/axioms13090584

APA Style

Sharma, S., Katiyar, K., Sudhamsu, G., Wratch, M. K., & Salgotra, R. (2024). Positivity-Preserving Rational Cubic Fractal Interpolation Function Together with Its Zipper Form. Axioms, 13(9), 584. https://doi.org/10.3390/axioms13090584

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

Article Metrics

Back to TopTop