A. K. B. Chand,P. Viswanathan
Abstract:Fractal interpolation functions (FIFs) developed through iterated function systems (IFSs) prove more versatile than classical interpolants. However, the applications of FIFs in the domain of `shape preserving interpolation' are not fully addressed so far. Among various techniques available in the classical numerical analysis, rational interpolation schemes are well suited for the shape preservation problems and shape modification analysis. In this paper, the capability of FIFs to generalize smooth classical interpolants, and the effectiveness of rational function models in shape preservation are intertwingly exploited to provide a new solution to the shape preserving interpolation problem in fractal perspective. As a common platform for these two techniques to work together, we introduce rational cubic spline FIFs involving tension parameters for the first time in literature. Suitable conditions on parameters of the associated IFS are developed so that the rational fractal interpolant inherits fundamental shape properties such as monotonicity, convexity, and positivity present in the given data. With some suitable hypotheses on the original function, the convergence analysis of the $\mathcal{C}^1$-rational cubic spline FIF is carried out. Due to the presence of the scaling factors in the rational cubic spline fractal interpolant, our approach generalizes the classical results on the shape preserving rational interpolation by Delbourgo and Gregory [SIAM J. Sci. Stat. Comput., 6 (1985), pp. 967-976]. The effectiveness of the shape preserving interpolation schemes are illustrated with suitably chosen numerical examples and graphs, which support the practical utility of our methods.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to combine fractal interpolation functions (FIFs) with shape - preserving interpolation to deal with data sets with complex and irregular structures. Specifically:
1. **Limitations of traditional interpolation methods**: For complex and irregular data sets (such as financial sequences, seismic data, bio - electrical signals, etc.), traditional interpolation methods may not provide satisfactory results.
2. **Requirement for shape - preserving interpolation**: In order to obtain an effective interpretation of physical processes, it is necessary to develop interpolation schemes that can inherit certain properties (such as monotonicity, convexity, positivity, etc.) in a given data set.
3. **Advantages of fractal interpolation**: Fractal interpolation functions (FIFs) are defined by iterative function systems (IFSs) and have self - similarity and flexibility, and can better fit data with different - scale details or some degree of self - similarity.
4. **Combining shape - preserving and fractal interpolation**: This paper introduces a new class of rational cubic spline fractal interpolation functions (rational cubic spline FIFs). These functions contain tension parameters and can utilize the flexibility of fractal interpolation while maintaining shape characteristics.
5. **Theoretical contributions**:
- Developed parameters suitable for the conditions, making rational fractal interpolation remain valid.
- Through convergence analysis, it is proved that the proposed rational cubic spline FIF can converge to the original data generation function under certain assumptions.
- Proposed a method for automatically selecting parameters to ensure that the interpolation results preserve the monotonicity and convexity of the data.
6. **Practical applications**: The effectiveness of the proposed method is demonstrated through numerical examples and graphics, especially when dealing with derivatives with varying irregularities, it is superior to traditional interpolation methods.
In summary, this paper aims to provide a new interpolation method for complex data sets, which can both maintain the shape characteristics of the data and utilize the flexibility of fractal interpolation, thereby improving the quality and applicability of interpolation results.
### Formula display
The formulas involved in the paper include, but are not limited to, the following:
- Definition of fractal interpolation function:
\[
f(x) = F_i\left(L_i^{-1}(x), f \circ L_i^{-1}(x)\right), \quad x \in I_i, \quad i \in J.
\]
- Form of rational cubic spline fractal interpolation function:
\[
S(L_i(x)) = \alpha_i S(x) + \frac{P_i(x)}{Q_i(x)},
\]
where,
\[
P_i(x) \equiv P_i^*(\theta) = (y_i - \alpha_i y_1)(1 - \theta)^3+\left[r_i y_i + h_i d_i-\alpha_i (r_i y_1 + d_1 (x_N - x_1))\right] \theta(1 - \theta)^2+\left[r_i y_{i + 1}-h_i d_{i + 1}-\alpha_i (r_i y_N - d_N (x_N - x_1))\right] \theta^2(1 - \theta)+(y_{i + 1}-\alpha_i y_N) \theta^3,
\]
\[
Q_i(x) \equiv Q_i^*(\theta) = 1+(r_i - 3) \theta(1 - \theta),
\]
\[
\theta=\frac{x - x_1}{x_N - x_1}.
\]
- Error estimate in convergence analysis:
\[
\|f - S\|_\infty \leq \frac{h^4 c \max_{i \in J}\left\{\vert y_i^{(1)} - d_i\vert,\vert y_{i + 1}^{(1)} - d_{i + 1}\vert\right\}}{384}+\frac{\vert h\vert^4 \|f^{(4)}\|_\infty