Abstract

Newton’s interpolation is a classical polynomial interpolation approach and plays a significant role in numerical analysis and image processing. The interpolation function of most classical approaches is unique to the given data. In this paper, univariate and bivariate parameterized Newton-type polynomial interpolation methods are introduced. In order to express the divided differences tables neatly, the multiplicity of the points can be adjusted by introducing new parameters. Our new polynomial interpolation can be constructed only based on divided differences with one or multiple parameters which satisfy the interpolation conditions. We discuss the interpolation algorithm, theorem, dual interpolation, and information matrix algorithm. Since the proposed novel interpolation functions are parametric, they are not unique to the interpolation data. Therefore, its value in the interpolant region can be adjusted under unaltered interpolant data through the parameter values. Our parameterized Newton-type polynomial interpolating functions have a simple and explicit mathematical representation, and the proposed algorithms are simple and easy to calculate. Various numerical examples are given to demonstrate the efficiency of our method.

1. Introduction

The interpolation problem has been the subject of many classic studies in approximation theory [15]. In the last several years, many researchers have been focusing on this subject and have obtained interesting results. The existing approaches can be divided into polynomial and rational interpolation methods, both of which are applicable to numerical approximation [1, 2], image interpolation [36], and arc structuring and surface modeling [710]. Newton’s polynomial interpolation and Thiele’s interpolating continued fractions can be incorporated to generate various interpolation schemes based on rectangular grids. As one can see, the (partial) inverse differences and the (partial) divided differences play a great role on the study of the polynomial and rational interpolation. Many scholars have constructed a variety of blending rational interpolation. Vonza [11] developed rational interpolation-based associated continued fractions. Varsamis and Karampetakis [12] presented recursive algorithms for the Newton polynomial interpolation of a given bivariate function. Tan et al. constructed Newton–Thiele [13] and Thiele–Newton blending rational interpolation [1]. However, blending rational interpolants strongly depend on the existence of so-called blending differences. By using the method of divide and conquer, Tan and Tang [14] proposed composite interpolation over triangular subgrids. By dividing the original set of support points into some blocks, Tan and Zhao proposed the block-based Thiele-like blending rational interpolation [15] and the block-based Newton-type blending rational interpolation [16]. Tang and Liang [17] developed the bivariate blending Thiele–Werner osculatory rational interpolation. Based on Thiele interpolating continued fraction, polynomial interpolation, and barycentric rational interpolation, many types of blending rational interpolation were studied. In recent years, Dyn and Floater [18] studied multivariate polynomial interpolation based on lower subsets. Cuyt and Salazar Celis [19] developed a generalized multivariate data fitting model to solve a variety of scientific computing problems, such as graphics, filtering, metamodeling, computational finance, and more, which captured linear as well as nonlinear phenomena. Interpolation was also applied to graphic image morphing and image processing [1, 36, 2023]. He et al. [3, 20] studied Thiele–Newton rational interpolation in the polar coordinates, applied it to image super-resolution, and obtained better performance. Karim and Saaban [21] proposed a novel rational bicubic ball function with one parameter by using tensor product approach and used it for image upscaling. Zhan et al. [22] developed a nonlocal and local image interpolation model based on nonlocal bounded variation regularization and local total variation and obtained better performance. He et al. [24] proposed an image inpainting algorithm by using continued fractions rational interpolation. In order to obtain better repaired results, Thiele’s rational interpolation was combined with Newton–Thiele rational interpolation to repair damaged images. In [25], the Newton interpolation method was adopted in capacitance gradient and the calibration of cantilever stiffness. Based on polynomial approximation of local surface, Rahman et al. [5] presented a novel local patch descriptor which is invariant to the changes in viewpoint, scale, orientation, and illumination.

The mathematical description and the construction method of the curve and surfaces are important problems in computer-aided geometric design [79, 26, 27]. There are many ways to handle this problem [1, 79, 28, 29], including Bézier and NURBS-based approaches as well as the polynomial spline method. These methods are used, for instance, for shaping the outer hull of aircraft and a ship [28]. However, the disadvantage of polynomial interpolation lies in its global property, that is, it is difficult to control the local constraint of a given interpolation point under the condition that the values of the interpolant points are fixed [29, 30]. To construct the polynomial spline, the derivative values of the interpolating data are required additionally to the function values. Unfortunately, in many practical fields, it is difficult to obtain these derivative values.

In recent years, many scholars have studied parametric spline interpolation. Zhang et al. [29, 30] proposed new types of weighted blending spline interpolation. By selecting appropriate parameters and different coefficients, the value of the spline interpolation function can be modified at any point in the interpolant interval, under the condition that the values of the interpolant points are fixed, so that the geometric surfaces can be adjusted. The drawback is that the computation is complicated. The bivariate rational interpolation with parameters, based only on the function values, has been studied in [29]. Based on function values and partial derivatives of the function being interpolated, Duan et al. [31] proposed a new bivariate rational interpolation, which had a simple and explicit rational mathematical representation with parameters. Based on function values, Duan et al. [32] constructed a rational cubic spline with parameters. This spline is C1-continuous in the interpolating interval. By selecting suitable parameters, it can also become C2-continuous. Reddy et al. [33] presented univariate and bivariate rational cubic fractal interpolation functions with shape parameters, which can be used for surface modeling when the partial derivatives of the given function are irregular. In [34], the authors presented a weighted bivariate rational bicubic spline interpolation based on function values, which is C1-continuous for any positive parameters. The shape of the interpolating surface can be modified by changing the parameters under the condition that the values of the interpolating nodes are fixed, which is convenient in engineering and useful in computer-aided geometric design [2834]. Tang and Zou [35] studied and presented general interpolant frames with many classical interpolation formulas. Univariate and bivariate Thiele-like rational interpolation continued fractions with parameters were studied in [36]. There, the interpolation function can deal with interpolation problems where inverse differences do not exist. Through the selection of proper parameter values, the value of the proposed interpolation function can be changed at any point in the interpolant region under the condition that the values of the interpolating nodes are fixed. However, how to select appropriate parameters is a hard problem. It is also difficult to define such an interpolant function. Polynomials have many advantages, such as a simple structure, easy differentiation, integration, and having derivatives of any order. The polynomial function is a good tool for approximating smooth function. In computer-aided geometric design, there is still a great demand for complicated models and the integration of design and fabrication, but there are still two unresolved problems [29, 30]:(1)How can one construct a suitable polynomial interpolation function with an explicit mathematical representation, simple calculation, which is also convenient to use and to study theoretically?(2)For given data, how can one remodel the curves or surfaces shape to meet the actual requirements in the computer-aided geometric design?

Given points and values , in the traditional interpolation method, the classical interpolating polynomial of degree less or equal to is unique to the given interpolation points, so it cannot be an answer to the above questions. At present, Newton’s polynomial interpolation is in the center of research on polynomial interpolation methods. While it has shown good interpolation performance, its disadvantage is that it is unique to the given interpolation points [28]. To overcome the above shortcomings, we propose a novel Newton interpolation polynomial only based on (partial) divided differences by introducing one or multiple parameters, which can be seen as a new approach to interpolation method of Newton polynomials. Similar to the interpolation functions in [2832], the proposed method allows for modifying the curve or surface shape under the condition that the values of the interpolant points are fixed. At the same time, the proposed method has many advantages, such as a simple and explicit mathematical representation and better performance. The interpolant function only depends on the values of the function being interpolated and the (partial) divided differences, so the computation is simple.

The rest of the paper is organized as follows. First, we develop one univariate parameterized Newton-type polynomial interpolation and give the interpolation algorithms and theorems in Section 2. Then, we present one kind of bivariate parameterized Newton-type polynomial interpolation and give the interpolation algorithms, theorems, and the dual interpolation polynomial in Section 3. In Section 4, to illustrate the effectiveness of the proposed algorithms, numerical examples are given.

2. Parameterized Univariate Newton-Type Polynomial Interpolation

We now develop the univariate Newton-type polynomial interpolation problems in one variable. Suppose function and we are given the support points set on interval . We obtain the following classic Newton polynomial representation [1]:where represent the divided differences of .

If the interpolation polynomial function exists for the given interpolation data points, then classical interpolating polynomial by of degree less or equal to is unique. Zhao and Tan [16] developed the block-based Newton-type blending rational interpolation. For a given set of interpolant data, they divided the interpolation into many subsets (blocks); based on different block point data, many types of polynomial or rational interpolation can be constructed. The block-based divided difference is calculated based on the block interpolation points first. Then, the block-based Newton-type blending rational interpolation over all blocks is constructed. This requires a lot of calculation, which is inconvenient for interpolation application. Hence, we construct a novel Newton interpolation polynomial with a single parameter based on one virtual point. The point is adjusted strategically, and the univariate Newton-type polynomial interpolation with one parameter is constructed.

By introducing a new parameter , an arbitrary point of the original points is regarded as a virtual double point, while the multiplicity of the other points remains the same.

Let

When , for ,

For ,

When , for ,

Based on equations (2)–(5), the novel divided differences are given as shown in Table 1.

Compared with the classical Newton polynomial interpolation, it is easy to see that the method in this paper is simple and convenient and is the same as that of the classical Newton polynomial interpolation.

The Newton-type polynomial interpolation with a parameter can be constructed:where

Theorem 1. Given the interpolation data ,

Proof. Suppose .
Equation (6) is the classical Newton polynomial interpolation, so For ,So, we haveTherefore, the results follow.
Similarly, we can also consider the Newton-type polynomial interpolation with two parameters with the proposed algorithm. It also can be divided into two cases. One case is that an arbitrary point from the original points is treated as a virtual treble point; another case is that any two points in the original points are regarded as two virtual double points. The Newton interpolation polynomial with more parameters can be constructed similarly.

3. Bivariate Parameterized Newton-Type Polynomial Interpolation

Now, we generalize some applications of the previous results to the multivariate Newton-type polynomial interpolation.

Suppose is a set of two dimensional points on the rectangular region D, is a bivariate real function defined on the rectangular region D, and let

The bivariate Newton polynomial interpolation is as follows:and represent the bivariate partial divided differences.

Theorem 2 (see [2]). Suppose is a set of two dimensional points on rectangular region D and is a real function defined on rectangular region D; then,

3.1. Bivariate Newton-Type Polynomial Interpolation with a Single Parameter

By introducing a new parameter , we use any point of the original points as a virtual double point, while the multiplicity of the other points remains 1. The Newton-type polynomial interpolation based on parameter can be constructed using Algorithm 1.

Algorithm 1. Step 1: initialization.Step 2: for ,Step 3: ,Step 4: by introducing parameter into the formula and using formulas (2)–(5), the final results areStep 5: using the data in formulas (16) and (17), the Newton-type polynomial interpolation with a single parameter with regard to can be constructed:Step 6: letThen, is a bivariate Newton-type polynomial interpolation with a single parameter.

Theorem 3. Given the interpolation points , the bivariate Newton-type polynomial interpolation with a single parameter satisfies

Proof. For an arbitrary point , if , it holds thatIf ,Regardless of , from Theorem 1, we can deriveFrom Theorem 4, ,Then, we have proved Theorem 3.
Similarity, we can also consider the Newton-type polynomial interpolation with two parameters. There are three cases: one case is that a point was viewed as a virtual treble point; another case is that two points were viewed as virtual double points in the same column; the third case is taking two points as virtual double points in different columns. The bivariate Newton-type polynomial interpolation with three or more parameters can be derived in a similar fashion.

3.2. Dual Bivariate Newton-Type Polynomial Interpolation with a Single Parameter

From Algorithm 1, according to the partial divided difference of x first and then concerning the partial divided difference of y, the novel bivariate Newton-type polynomial interpolation was constructed. In fact, we can also construct it by first regarding the partial divided difference of and then the partial divided difference of . By introducing a new parameter , any point in the original interpolant points is regarded as a virtual double point. Then, the Newton-type polynomial interpolation with parameter can be obtained using Algorithm 2.

Algorithm 2. Step 1: initialization.Step 2: if , letStep 3: for ,Step 4: by introducing a parameter into the and calculating them using formulas (2)–(5), the final results areStep 5: by using the data in formulas (27) and (28), the univariate Newton-type polynomial interpolation with a single parameter in regard to x can be constructed as follows:Step 6: letThen, is a dual bivariate Newton-type polynomial interpolation with a single parameter.
We call as the dual interpolation of . As a matter of fact, the bivariate Newton-type polynomial interpolation with a single parameter can also be calculated by the information matrix method. This algorithm borrows elementary transformation thought of linear algebra, so it is easy to implement. We can calculate the partial divided differences in Algorithms 1 and 2 in a manner of the information matrix method. Here, we only give the information matrix method of the dual bivariate Newton-type polynomial interpolation with a single parameter. Algorithm 2 can also be described by Algorithm 3 based on information matrix as follows.

Algorithm 3. Step 1: initialization. LetOne can define the following initial information matrix:Step 2: for , letThe above recursive process aims to change the initial information matrix into the matrix by using row transformation:Step 3: for , letThe above recursive process aims to change the columns in the initial information matrix into the matrix by using column transformation:By introducing a parameter into the in the information matrix and calculating them using formulas (2)–(5), the final results areStep 4: by using the data in formulas (36) and (37), the univariate Newton-type polynomial interpolation with a single parameter in regard to x can be constructed as follows:Step 5: letThen, is the dual bivariate Newton-type polynomial interpolation with a single parameter.

Theorem 4. Given the interpolation data , the dual bivariate Newton-type polynomial interpolation with a single parameter satisfies

Proof. For an arbitrary point , if , we haveIf ,Regardless of , from Theorem 1, we haveSo, we haveThen, we can obtain the result.
In addition, it can also construct other dual bivariate parameterized Newton-type polynomial interpolations, similar to the proposed algorithm in Section 3.2. Furthermore, according to the previous discussion, we can also construct parameterized Newton-type polynomial interpolation in both and directions. The formulas of every new parameterized Newton-type polynomial interpolation depend on the selected parameters. We can construct different parameterized Newton-type polynomial interpolation functions according to the application needs. At the same time, the value of the interpolant function can be changed at any point in the interpolant region under the condition that the values of the interpolant points are fixed by selecting appropriate parameters.
As a matter of fact, one can also construct a one-parameter family of interpolating polynomials:with being the classical interpolating polynomial and .
The proposed parameterized Newton-like interpolation polynomial is another approach to the one-parameter family of interpolating polynomials. In this paper, the proposed algorithms are only based on (partial) divided differences, which seems a little complicated. In fact, the calculation is the same as that of the classical Newton interpolation polynomial. We can also discuss the univariate case with more parameters and bivariate case with one or more parameters similarly. The algorithm given in this paper describes a new approach method for Newton interpolation polynomial.
The parameterized Newton-type polynomial interpolation defined in this paper is a local interpolation in each subinterval or subrectangle, which only depends on the function values. There are some parameters in each interpolating region, and when the parameters vary, the interpolant function will be changed under the condition that the values of the interpolant points are fixed. Therefore, the interpolation curves and surfaces will change with different parameters. Thus, the shape of the interpolant curves or surface can be adjusted for the given interpolation points. However, the question of how to select the optimal parameters of the proposed interpolant algorithms is still a challenging problem, which we will consider in our future work.

4. Numerical Examples

In this section, we give some examples to demonstrate how the proposed algorithms are implemented and how the parameters can be adjusted to change the curve or surface shape.

Example 1. Let the interpolation data be as given in Table 2.
From the interpolation data in the Table 2, we get the Newton interpolation polynomial asTake as a virtual double point, and the multiplicity of the other points remains 1. The Newton-type polynomial interpolation with one parameter derived via Algorithm 1 isIt is easy to verify thatIf the interpolant points are given, the shape of interpolating curves is fixed. This is because of the uniqueness of the interpolation function. The parameterized Newton-type polynomial interpolation methods defined in this paper are parametric. By changing the parameters, the interpolation function can be adjusted while observing the condition that the values of the interpolant points are fixed. The locality constraint in the interpolation curves is to make the function value at a specific point equal to a given real number in actual design [29, 30]. The central point in the interpolating curve is more important than other points for shape control [29, 30], so we consider the central point value control problem.
In fact, if we take M = 1.5, based on the parameterized Newton-type polynomial interpolation, we can get , and the function value of the central point is . Figure 1(b) shows the plot of the interpolation . If we take M = 0, it is easy to prove , and Figure 1(c) shows the plot of the interpolation . In this case, the function value of central point is . If we take M = −1, it is easy to prove , and Figure 1(d) shows the plot of the interpolation ; the function value of central point is .
By choosing different values of parameter , the function invariably satisfies the given interpolation data. The value of the interpolation function can be changed in the interpolant region by choosing appropriate parameter .

Example 2. Let the interpolation data be as given in Table 3.
We obtain the bivariate Newton polynomial interpolation by the method in [2] asFollowing Algorithm 1 in this paper, we increase the multiplicity of node to 2, the multiplicity of the other points remains 1, and then we construct the Hermite interpolation with first-order derivative in point by introducing parameter . We obtain the Newton-type polynomial interpolation with parameter with regard to :Following Algorithm 2, we derive the Newton-type polynomial interpolation with parameter with regard to :Following Algorithms 1 and 2, we derive the Newton-type polynomial interpolation with parameters with regard to and :We can also obtain the bivariate Newton polynomial interpolation by adding times with parameter :It is easy to verifyIt is easy to verify that the function invariably satisfies the given interpolation data with the different values of parameter . Figure 2 shows the plot of the classical bivariate Newton polynomial interpolation .
Similar to interpolant curves, we consider the local point control problem of the interpolation surface, that is, how to make the function value of the interpolation at a point equal to a given real number S in the actual design. We cannot adjust the shape of the classical Newton interpolation polynomial . We can use the algorithm in this article to solve this problem. Compared with other points, the center point is the key area of the interpolation patch. The function value of the central point is . If we set S = −1, based on the parameterized Newton-type polynomial interpolation, we can get , and the shape of the surface will “drop.” In this case, the function value of the central point is . If we set S = 4, it is easy to prove , and the shape of the surface “rises.” In this case, the function value of central point is . Figures 3(a) and 3(b) show the plot of the parameterized Newton-type polynomial interpolation with and , respectively. Similarly, we can set S = −1 and S = 4 for and ; we can discuss polynomial interpolation , , and . Figures 4(a) and 4(b) show the plot of the parameterized Newton-type polynomial interpolation with and , respectively. Figures 5(a)5(d) show the plot of the parameterized Newton-type polynomial interpolation with and and and , respectively. In fact, we can choose other combinations of parameters that satisfy certain conditions, so that and . Figures 6(a) and 6(b) show the plot of the parameterized Newton-type polynomial interpolation with and , respectively.
For the given interpolation points, we get four Newton polynomial interpolations with one or two parameters. The polynomial has a higher degree than the proposed method in this paper. In terms of the central point value control problem, the proposed method in this paper has a large adjustment space as shown in Figures 35. Meanwhile, we can adjust the shape of surface in one or in two directions or add parameters to adjust in both directions according to actual needs. The interpolation function can be adjusted under the condition that the values of the interpolant points are fixed. This is achieved by adjusting the parameters of our Newton-type polynomial interpolation methods. We can also modify the new Newton-type polynomial interpolation by parameter , so our methods give a new choice for the application and a new method for studying the polynomial interpolation theory. According to the other algorithms in this paper, we can also construct many other novel Newton-type polynomial interpolations with more parameters similarly.

5. Conclusions and Future Work

We presented a new approach to interpolation method of Newton polynomials. We constructed several kinds of univariate and bivariate parameterized Newton-type polynomial interpolation. We discussed the interpolation theorem, the algorithms, the dual bivariate parameterized Newton-type polynomial interpolation, and algorithm based on information matrix. Our methods are easy to be used and theoretically sound. The value of the Newton-type polynomial interpolant function can be adjusted in the interpolant region by choosing appropriate parameter values, while observing the condition that the values of the interpolant points are fixed. According to the actual geometric design needs, the shape of the interpolation curves or surfaces can be adjusted, which is convenient in engineering and useful in computer-aided geometric design [28, 30]. The proposed method has advantages, such as a simple and explicit mathematical representation and ease of computation. But how to select the optimal parameter for application is a difficult problem. Our future work will concentrate on the following aspects:(1)How to select the optimal parameter values and suitably adjust the shape of the interpolant curves or surface according to the actual design requirements.(2)Investigating how the presented parameterized Newton-type polynomial interpolation method can be applied in other pixel-level image processing steps, such as super-resolution reconstruction, image upscaling, image rotation, removal of salt and pepper noise, image metamorphosis, and image inpainting.(3)How to construct multivariate composite interpolation and blending rational interpolation over lacunary, triangular, pyramid-type, and other irregularity grids based on the parameterized Newton-type polynomial interpolation, Thiele continued fraction interpolation, and other classical interpolation methods. The proposed method can also be generalized to Hermite interpolation and osculatory rational interpolation.

In summary, the parameterized Newton-type polynomial interpolation algorithm has many advantages, such as a simple and explicit mathematical representation and ease of computation; the value of the Newton-type polynomial interpolant function can be adjusted in the interpolant region by choosing appropriate parameter values; according to the actual geometric design needs, the shape of the interpolation curves or surfaces can be adjusted.

In addition, with the help of the Samlson generalized inverse, it is easy to generalize the parameterized Newton-type polynomial interpolation method to vector-valued cases or matrix-valued cases [1, 15, 16].

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (nos. 61672204, 61673359, and 61806068), the Natural Science Foundation of Anhui Province (nos. 1908085MF184 and 1908085QF285), in part by the Key Scientific Research Foundation of Education Department of Anhui Province (nos. KJ2018A0555 and KJ2019A0833), and in part by the Key Technologies R&D Program of Anhui Province (no. 1804a09020058).