Juncheng Li,Laizhong Song,and Chengzhi Liu

Abstract—A class of cubic trigonometric interpolation spline curves with two parameters is presented in this paper.The spline curves can automatically interpolate the given data points and become C2interpolation curves without solving equations system even if the interpolation conditions are fixed.Moreover,shape of the interpolation spline curves can be globally adjusted by the two parameters.By selecting proper values of the two parameters,the optimal interpolation spline curves can be obtained.

I.INTRODUCTION

IN geometric modeling,curves are usually constructed on the basis of polynomials.Recently,trigonometric polynomials have received much attention within geometric modeling.Some examples are the trigonometric polynomial B′ezier curves[1]-[4],the trigonometric polynomial B-spline curves[5]-[8],and the trigonometric Hermite spline curves[9],[10].Those trigonometric polynomial curves not only inherit similarities of the corresponding classical polynomial curves,but also have better performance ability in different aspects.

It is well known that the spline interpolation has a wide range of applications in geometric modeling.However,when the traditional cubic splines are used to construct C2interpolation curves,a linear equation system usually needs to be solved.In order to solve this problem,Su and Tan[11]and Yan and Liang[12]constructed two automatic interpolation splines on the basis of trigonometric polynomials.The main characteristic of these two splines is that they can automatically interpolate the given data points without solving equation systems.There are some deficiencies in these two splines,however,the shape of the spline presented in[11]cannot be adjusted when the data points and auxiliary points arefixed.Although the spline presented in[12]can achieve shape adjustment when the data points and auxiliary points are fixed,it can only become G2interpolation curve.The main purpose of this paper is to present a class of cubic trigonometric automatic interpolation spline curves with two parameters.The proposed spline curves can achieve shape adjustment by altering values of the two parameters.Specifically,the proposed spline curves can automatically interpolate the given data points and be C2without solving equations system even if the interpolation conditions are fixed.Furthermore,the optimal automatic interpolation spline curves can be obtained by selecting proper values of the two parameters.

The rest of this paper is organized as follows.In Section II,the cubic trigonometric automatic interpolation spline curves with two parameters are constructed,and properties of interpolation spline curves are discussed.In Section III,a method for determining the optimal automatic interpolation spline curves is presented,and some examples are given.A short conclusion is given in Section IV.

II.THECUBICTRIGONOMETRICAUTOMATIC INTERPOLATIONSPLINECURVES

A.The Basis Functions

Firstly,the cubic trigonometric spline basis functions with two parameters are defined as follows.

Definition 1:For 0≤ t≤ 1,α,β ∈R,the following four functions about t are defined as cubic trigonometric spline basis functions,

By simple deduction,the cubic trigonometric spline basis functions defined as(1)have the following properties at the endpoints,

From(1),the cubic trigonometric spline basis functions have two free parameters α and β.Hence,various shapes of curves of the cubic trigonometric spline basis functions can be obtained by altering values of the parameters α and β.Fig.1 shows curves of the cubic trigonometric spline basis functions for different values of α and β,where α and β are taken as(α,β)=(-0.5,0.5)(marked with dotted lines),(α,β)=(0,0)(marked with solid lines),(α,β)=(0.5,-0.5)(marked with dashed lines).

Fig.1.Curves of the proposed basis functions for different α and β.

B.Definition and Properties of the CTAI-spline Curves

Depending on the cubic trigonometric spline basis functions,the cubic trigonometric automatic interpolation spline curves with two parameters can be defined as follows.

Definition 2:Given a series of data points(xi,yi)(i=0,1,...,n),set hi= xi+1-xi(i= 0,1,...,n-1),the following piecewise curves are called cubic trigonometric automatic interpolation spline curves(CTAI-spline curves for short),

where x∈[xi,xi+1](i=1,2,...,n-2),0,1,2,3)are the cubic trigonometric spline basis functions defined according to(1).

Theorem 1:The CTAI-spline curves defined as(5)have the following properties,

a)Automatic interpolation property:For given data points(xi,yi)(i=0,1,...,n),the CTAI-spline curves Si(x)(i=1,2,...,n-2)automatically interpolate all the given data points except(x0,y0)and(xn,yn),i.e.,

b)C2continuity:When hi(i=0,1,...,n-1)are equal to a constant h,the CTAI-spline curves Si(x)(i=1,2,...,n-2)are C2,i.e.,

Proof:a)By(2)and(5),we have

Equation(6)shows that the CTAI-spline curves Si(x)(i=1,2,...,n-2)automatically interpolate all the given data points except(x0,y0)and(xn,yn).

b)When hi(i=0,1,...,n-1)are equal to a constant h,by(3)and(5),we have

By(4)and(5),we have

From(6)-(8),we have

Equation(9)shows that the CTAI-spline curves Si(x)(i=1,2,...,n-2)are C2.

From Theorem 1,if two auxiliary points(x-1,y-1)and(xn+1,yn+1)are added to the given data points,the CTAI-spline curves Si(x)(i=0,1,...,n-1)interpolating all the data points(xi,yi)(i=0,1,...,n)would be naturally obtained.The two auxiliary points can be added arbitrarily,and they only have effects on the starting and the end segments of the interpolation curves.For convenience,(x-1,y-1)and(xn+1,yn+1)can be taken as

It is clear that there exist two degree of freedom in the CTAI-spline curves,even if the data points and auxiliary points are fixed.Different shapes of the interpolation curves can be obtained by altering the parameters α and β.

Example 1:Given data points

and two auxiliary points(x-1,y-1),(x5,y5)are added to the data points according to(10).Then,the C2CTAI-spline curves with different values of parameters α and β are shown in Fig.2,where the parameters are taken as(α,β)=(1,-1)(marked with dotted lines),(α,β)=(0,0)(marked with solid lines),(α,β)=(-1,1)(marked with dashed lines).

Fig.2.Different shapes of the C2CTAI-spline curves for data points.

Example 2:Given a function y=r(x)=1/(1+x2)(-5≤x≤5),set xi=-5+i,yi=r(xi)(i=0,1,...,10),and two auxiliary points(x-1,y-1)and(x11,y11)are added to the data points according to(10).Then,the curve of the given function(marked with solid lines)and the C2CTAI-spline curves with different values of parameters α and β are shown in Fig.3,where the parameters are taken as(α,β)=(-1,1)(marked with dotted lines),(α,β)=(1,-1)(marked with dashed lines).

Fig.3.Different shapes of the C2CTAI-spline curves for a function.

Example 3:Given proper data points and added two auxiliary points,some special curves generated by the G2CTAI-spline curves can be obtained.The shape of the curves also can be adjusted by altering values of the parameters α and β,see in Figs.4 and 5.

Remark 1:When the traditional cubic polynomial splines are used to construct C2interpolation curves,the general way is to reverse the linear equations system.However,due to the automatic interpolation property and continuity of the CTAI-spline,the C2interpolation curves can be naturally obtained without solving a system of equations.On the other hand,when the data points and auxiliary points are fixed,the traditional cubic polynomial spline interpolation curves are unique,while shapes of the CTAI-spline curves can be adjusted by altering values of the parameters α and β.

Fig.4.The note-like curves generated by the CTAI-spline curves.

Fig.5.The vane-like curves generated by the CTAI-spline curves.

Remark 2:On one hand,compared with some trigonometric spline curves(such as[6]-[8]),the main characteristic of the CTAI-spline curves is that they can automatically interpolate the given data points without solving systems of equations.On the other hand,compared with the two trigonometric automatic interpolation splines presented in[11]and[12],the CTAI-spline has the following characteristics,

a)The shape of the trigonometric automatic interpolation spline presented in[11]is unique when the data points and auxiliary points are fixed,while the CTAI-spline curves can achieve shape adjustment by parameters α and β,even if the data points and auxiliary points are unchanged.

b)Although the shape of the trigonometric automatic interpolation spline presented in[12]can be adjusted when the data points and auxiliary points are fixed,the interpolation curves are only G2.However,the CTAI-spline curves are not only adjustable in shape,but also can become C2interpolation curves.

III.THEOPTIMALCTAI-SPLINECURVES

As mentioned above,suppose a function y=r(x)(x0≤x≤xn),set xi=x0+hi(h＞0),yi=r(xi)(i=0,1,...,n),and two auxiliary points(x-1,y-1),(xn+1,yn+1)are added to the data points according to(10),the corresponding C2CTAI-spline curves would be completely determined by the parameters α and β.Hence,bad interpolation curves would be generated if the parameters α and β are chosen improperly.

Example 4:Given a function y=r(x)=(sinx)/x(1≤x≤5),set xi=1+i/2,yi=r(xi)(i=0,1,...,8),and two auxiliary points(x-1,y-1)and(x9,y9)are added to the data points according to(10).Fig.6 shows the curve of the given function(marked with dotted lines)and the C2CTAI-spline curves with different values of parameters α and β.

It is obvious that the interpolation curves in Fig.3(a)are more satisfactory than the curves in Fig.3(b).Hence,how to determine proper values of the parameters α and β is of great importance when constructing the C2CTAI-spline curves.A method for determining values of the parameters α and β is presented as follows.

Fig.6.Effects of the parameters on the C2CTAI-spline curves.

Let

where α, β ∈ R,Si(x)is the corresponding CTAI-spline curves defined in the interval[xi,xi+1](i=0,1,2,...,n-1).

In order to obtain the minimum interpolation error,we have

Set

Then,(5)can be rewritten as follows:

where

By(13),(11)can be rewritten as follows:

where

From(14),(12)can be rewritten as follows:

If C1C2-/=0,we have

Remark 3:If C1C2-=0,then(15)has no unique solution.In this case,the data points could be changed to guarantee that C1C2-/=0.

After the optimal parameters α and β are calculated from(16),the optimal C2CTAI-spline curves Si(x)(i=0,1,...,n-1)interpolating the given function y=r(x)can be obtained.

Since the trigonometric automatic interpolation spline curves presented by[11]are C2,and the trigonometric automatic interpolation spline curves presented by[12]are G2,so,the spline curves presented by[11]are chosen to compare with the optimal C2CTAI-spline curves in the following two examples.

Example 5:Given a function y=r(x)=1/(1+x2)(-5≤x≤5),set(i=0,1,...,n),and two auxiliary points(x-1,y-1)and(xn+1,yn+1)are added to the data points according to(10).For different n(where n+1 is the number of the data points),the optimal parameters of the C2CTAI-spline curves are calculated from(16).Interpolation errors of the optimal CTAI-spline curves and the spline curves presented by[11]are shown in Table I.

TABLE I INTERPOLATIONRESULTS FORDIFFERENTnOFEXAMPLE5

The curve of the given function(marked with solid lines),the optimal C2CTAI-spline curves(n=10,marked with dotted lines)and the spline curves presented by[11](n=10,marked with dashed lines)are shown in Fig.7.To facilitate observation,the curve of the given function is increased by 0.1 units and the spline curves presented by[11]are increased 0.2 units in Fig.7.

Fig.7.Interpolation results of Example 5(n=10).

Example 6:Given a function y=r(x)=(sinx)/x(1≤x≤5),set(i=0,1,...,n),and two auxiliary points(x-1,y-1)and(xn+1,yn+1)are added to the data points according to(10).For different n(where n+1 is the number of the data points),the optimal parameters of theC2CTAI-spline curves are calculated from(16).Interpolation errors of the optimal CTAI-spline curves and the spline curves presented by[11]are shown in Table II.

TABLE II INTERPOLATIONRESULTS FORDIFFERENTnOFEXAMPLE6

The curve of the given function(marked with solid lines),the optimal C2CTAI-spline curves(n=8,marked with dotted lines),and the spline curves presented by[11](n=8,marked with dashed lines)are shown in Fig.8.To facilitate observation,the curve of the given function is increased 0.1 units and the spline curves presented by[11]are increased 0.2 units in Fig.8.

Fig.8.Interpolation results of Example 6(n=8).

Example 5 and Example 6 show that the interpolation result of the optimal C2CTAI-spline curves is better than the trigonometric automatic interpolation spline curves presented by[11],which is due to the parameters α and β of the CTAI-spline curves observing their optimal values.

IV.CONCLUSIONS

A class of cubic trigonometric automatic interpolation spline curves with two parameters is presented in this paper.The spline curves can not only automatically interpolate the given data points and be C2without solving a system of equations,but can also achieve shape adjustment by altering values of the two parameters,even when the interpolation conditions are fixed.The examples illustrate the feasibility of this method.Furthermore,cubic trigonometric automatic interpolation spline surfaces with shape parameters can be constructed,and they enjoy the same characteristics as the interpolation spline curves.The corresponding parametric spline curves and surfaces are discussed as well.