A Generalized Quasi Cubic Trigonometric Bernstein Basis Functions and Its B-Spline Form

: In this paper, under the framework of Extended Chebyshev space, four new generalized quasi cubic trigonometric Bernstein basis functions with two shape functions α ( t ) and β ( t ) are constructed in a generalized quasi cubic trigonometric space span { 1,sin 2 t , ( 1 − sin t ) 2 α ( t ) , ( 1 − cos t ) 2 β ( t ) } , which includes lots of previous work as special cases. Sufficient conditions concerning the two shape functions to guarantee the new construction of Bernstein basis functions are given, and three specific examples of the shape functions and the related applications are shown. The corresponding generalized quasi cubic trigonometric Bézier curves and the corner cutting algorithm are also given. Based on the new constructed generalized quasi cubic trigonometric Bernstein basis functions, a kind of new generalized quasi cubic trigonometric B-spline basis functions with two local shape functions α i ( t ) and β i ( t ) is also constructed in detail. Some important properties of the new generalized quasi cubic trigonometric B-spline basis functions are proven, including partition of unity, nonnegativity, linear independence, total positivity and C 2 continuity. The shape of the parametric curves generated by the new proposed B-spline basis functions can be adjusted flexibly. For equidistant knots, the R i ( u ) and the knot vector U are called a uniform generalized quasi cubic trigonometric B-spline basis function and a uniform vector. For non-equidistant knots, R i ( u ) and U are called that of non-uniform ones, respectively. Figure 7 gives some plots of generalized quasi cubic trigonometric B-spline basis functions with different shape parameters with α ( t ) and β ( t ) given in Equation (13) and Equation (14).


Introduction
The construction of basis functions with good properties, including partion of unity, nonnegativity and total positivity, is a basic subject within computer aided geometric design (GAGD) and Computer Graphics (CG). Basis functions with desirable properties have an important role in curves and surfaces construction.
In order to control the shape of the curves flexibly, various spline curves possessing shape parameters have been proposed. In [1], a kind of G 2 continuous rational cubic spline curves with tension shape parameters was developed. In [2], a kind of G 2 continuous Beta-spline curves with local bias and tension parameters was proposed. In [3], in the polynomial space {1, t, (1 − t) p , t q }, a kind of variable degree spline curves was constructed. In [4], a kind of changeable degree spline curves was constructed. In [5], new cubic rational B-spline curves with two shape parameters were constructed. In the space span 1, 3t 2 − 2t 3 , (1 − t) α , t β , a class of αβ-Bernstein basis functions and the related B-spline form with two local shape parameters was proposed in [6]. In [7], a kind of λµ-Bernstein basis functions with two exponential functions was developed in the space span{1, 3t 2 − 2t 3 , (1 − t) 3 e −λt , t 3 e −µ(1−t) }.
Blossom theory plays a significant role in the construction of B-splines, which is first stated in [8]. In [9], Mazure and Laurent proposed a general notion of blossom which can be extended to Chebyshev splines. In [10], the necessary and sufficient conditions for the existence of blossoms for splines with connection matrices and with sections in different four-dimensional quasi-Chebyshev spaces were explored. In [11], Mazure showed how the existence of blossoms in some spaces automatically leads to optimal bases. The approximation power and the conditions of existence of a normalized B-basis were further 2 of 25 studied in [12]. In [13], the dimension elevation process from Quasi Extended Chebyshev spaces was analyzed and a new insight into some special cases was given. In [14], a new space of variable degree polynomial was proposed, where a Bernstein basis was shown and some generalized tension properties were given.
The purpose of this paper is to propose a new kind of generalized quasi cubic trigonometric Bernstein basis functions with two shape functions α(t) and β(t) constructed in the space T α,β := span 1, sin 2 t, (1 − sin t) 2 α(t), (1 − cos t) 2 β(t) , which form a normalized optimal totally positive basis and include the functions in [15][16][17][18][19][20] as special cases. In [23][24][25], the concept of normalized B-basis related to the normalized optimal totally positive basis is mentioned, which corresponds to the optimal shape preserving properties. This kind of new generalized quasi cubic trigonometric Bernstein basis functions has good generality and wide applications, as different shape functions α(t) and β(t) lead to different kinds of spaces mentioned in [15][16][17][18][19][20]. In order to construct more new basis functions conveniently, we propose some simplified constructive conditions. Based on those simplified constructive conditions, we construct three specific examples of the shape functions α(t) and β(t), and give the corresponding applications of those three examples. What is more, after verification, the space mentioned in [19] is a special case that satisfies the simplified constructive conditions. Moreover, the corresponding quasi cubic trigonometric B-spline basis functions with two local shape functions α i (t) and β i (t) are also constructed, and the shape functions α i (t) and β i (t) satisfy some additional conditions. The rest of this paper is organized as follows. In Section 2, the meterial and methods are introduced. In Section 3, the construction and properties of the generalized quasi trigonometric Bernstein basis and the corresponding quasi trigonometric Bézeir curves are given, including the simplified constructive conditions that the two shape functions satisfy.
Three specific examples are given. Section 4 gives the definitions and properties of the generalized quasi cubic trigonometric B-spline basis functions and corresponding curves, including the additional conditions that the two shape functions satisfy. In Section 5, our conclusions are given.

Bernstein Basis Functions and Bézeir Curves
Definition 1. Given control points P i ∈ R 2 or R 3 (i = 0,1,..., n), the Bézier curve can be defined as where the Bernstein basis functions B i,n (t), i = 0, 1, 2, ..., n are defined as follows Some important properties of the Bernstein basis functions B i,n (t), i = 0, 1, 2, ..., n are listed as follows. ( (F) Total positivity. For any sequence of points 0 ≤ t 0 < t 1 < ... < t n ≤ 1, the collocation matrix (B j,n (t i )) 0≤i,j≤n is totally positive, which means that all its minors are nonnegative.
From the properties of Bernstein basis functions above, the following properties of Bézier curves can be obtained.
(A) Affine invariance. Since the Bernstein basis functions have the property of partition of unity, the corresponding Bézier curve has the property of affine invariance.
(B) End-point property. The starting point and end point of Bézier curve coincide with P 0 and P n .
(C) Convex hull property. Since the Bernstein basis functions have nonnegativity and partition of unity, the Bézier curve must lie inside the control polygon spanned by the control points.
(D) Symmetry. A new Bézier curve is constructed from the control vertex P i = P n−i , i = 0, 1, . . . , n, which has the same shape as the original Bézier curve and the opposite direction.
(E) Variation diminishing property. Since the Bernstein basis functions form a group of optimal normalized totally positive bases, the corresponding Bézier curve possesses variation diminishing property, that is no plane intersects a Bézier curve more often than it intersects the corresponding control polygon.

Extended Chebyshev Space and Blossom Theory
In this subsection, we shall give a brief presentation of Chebyshev spaces and blossom theory, which will be used further in the following construction. Definition 2. Given a real interval I, we call a k-dimensional space U as an Extended Chebyshev space (EC space) on I if any nonzero element of U has at most k − 1 zeros (counted with multiplicities) in I. It is said to be a complete extended Chebyshev space (ECC space) on I if there exists a nested sequence where U i is an i-dimensional EC space on I, for i = 1, 2, ..., k. The related proofs are given in [26].

Remark 2.
A given k-dimensional subspace U is an ECC space on I iff there exists k positive functions w i (i = 1, 2, ...k) such that U = KerD • L k , where D stands for the ordinary differentiation and w i (i = 1, 2, ...k) are called weight functions associated with U. For L 1 , ..., L k , the differential operators are defined by Blossom is a useful tool for curves and surfaces construction, and it is especially powerful when used to analyze splines. More details are shown in [8]. Using means of intersections of convenient osculating flats, it is possible to define a notion of blossom.

Definition 3.
A geometrically regular function of order k, Φ : I → A, is said to be a piecewise smooth Chebyshev function of order k on I, if for all distinct points τ 1 , τ 2 , ..., τ r ∈ I and all positive integers µ 1 , µ 2 , ..., µ r whose sum is equal to k, the affine flat r i=1 Osc k−µ i Φ(τ i ) consists of a single point. The function ϕ : I k → A so defined is called the blossom of Φ. It is a symmetric function and it satisfies Let ε be the associated space of Φ. For any affine map h: affIm(Φ) → C, the blossom of the ε-function F = h • Φ will be defined by The detailed proofs and properties are given in [27].
In order to construct more new basis functions conveniently, here we give some simplified constructive conditions that α(t), The expression of the associated mother-function is Firstly, we are about to show that the following space is a 3-dimensional EC-space on [0, π/2]. Then by Theorem 2.1 of [28], we can see the existence of blossom in T α,β , which implies that the new space T α,β is suited for curves design.
Next, we are about to prove that the space DT α,β is also an EC-space on [0, π/2]. We only need to prove that any nonzero element of the space DT α,β has at most two zeroes (counting multiplicities as far as possible up to 2) on [0, π/2]. For t ∈ [0, π/2], consider the nonzero function Since the space DT α,β is an ECC-space in (0, π/2), F(t) has at most 2 zeroes in (0, π/2). If we assume that the function F(t) vanishes at 0, then we can obtain that C 1 = 0. In this case, if C 2 = 0, then F(t) has a singular zero at 0 and a singular zero at 1. If C 0 = 0, it is obvious that 0 is a double zero of F(t). If C 0 C 2 > 0, F(t) vanishes only one time at 0 and it does not vanish anywhere else on (0, π/2]. If C 0 C 2 < 0, F(t) vanishes only one time at 0 and it does not vanish at π/2. For t ∈ (0, π/2), the computation of the second derivative of F(t) gives that is a monotonic function in (0,π/2). By Rolle's theorem, we can conclude that F(t) has at most one zero in (0,π/2). Similarly, by analyzing the case that F(t) vanishes at π/2, we can also obtain the conclusion that F(t) has at most two zeroes on [0, π/2]. Summarizing, the space DT α,β is an EC-space Since the space DT α,β is an EC-space on [0, π/2], by Theorem 4.1 of [28], we can obtain that the space T α,β has a kind of normalized Bernstein basis on [0, π/2]. In the following theorem, we shall give the expression of Chebyshev-Bézier points of the mother-function Φ(t) defined in Equation (2) and the associated generalized quasi cubic trigonometric Bernstein basis Theorem 2. The four Chebyshev-Bézier points P i := P i (0, π/2) of the mother-function Φ(t) defined in Equation (2) are given by Moreover, the associated four generalized quasi cubic trigonometric Bernstein basis functions Proof of Theorem 2. From the definition of the mother-function Φ(t) given in Equation (2), we have Thus, by simply computing, we have we can easily deduce the expressions of the generalized quasi cubic trigonometric Bernstein basis B i (t), i = 0, 1, 2, 3.

Remark 3.
By analyzing the expressions of the generalized quasi cubic trigonometric Bernstein basis functions B i (t) (i = 0, 1, 2, 3) given in Equation (5), it is easy to conclude that these basis functions possess the properties of linear independence, nonnegativity and partition of unity. What is more, for i = 0, 1, 2, 3, the end-point properties are given as follows Moreover, B 0 (t) vanishes two times at π/2 and B 3 (t) vanishes two times at 0 (counting multiplicities as far as possible up to 2). Therefore, by Definition 2.10 of [28], we can see that the system (B 0 (t), B 1 (t), B 2 (t), B 3 (t)) is precisely the normalized Bernstein basis of the space span T α,β , which indicates that the generalized quasi cubic trigonometric Bernstein basis is totally positive and possesses optimal shape preserving properties. In [23][24][25], the concept of normalized B-basis conincides with the normalized Bernstein basis, which corresponds to the optimal shape preserving properties.
These shape preserving properties mainly include monotonicity preserving and variation diminishing, with more details stated in [23][24][25]29,30]. On the one hand, monotonicity preserving indicates that the curve has the same monotonicity as the monotone control points, and the length of the generated parametric curves is bounded above by the length of their control polygon. If the control polygon is planar and convex, then the curve is also planar and convex. On the other hand, the basis functions have variation diminishing property, so the curve never crosses a hyperplane more often than its control polygon does.
There are many shape functions that satisfy Equation (1). Here we analyze three specific examples and show some plots of these examples. In order to state conveniently, in the following discussion, we will also denote the four generalized quasi cubic trigonometric Bernstein basis functions as Example 1. We consider the simple example without parameters. The expressions of the shape functions are Next, we analyze that these two shape functions of Equation (6) satisfy Equation (1) mentioned above. By substituting the endpoint value 0 and π/2 and calculating the corresponding derivation, Then, we will prove α(t) = 2 π−2 (π/2 − t − cos t) ≥ 0. By simple analysis, we have For this, we only need to prove that So f 1 (t) reaches the maximum value when t = arc sin( 1 4 ) and reaches the minimum value when t=π/2, since  and β(t) given in Equation (6).

Example 2.
We consider the shape functions with denominator parameters. The expression of the shape functions are where λ, µ ≥ 1/2, p, q = 0, 1. Next, we analyze that these two shape functions satisfy Equation (1) mentioned above. When p = q = 0, it has been stated in [6] in detail. Then we prove the case of p = q = 1. The other cases are the same. It is obvious that Then, we further prove that α (t) − 1 3 α(t) ≥ 0. By simple calculations, we only need to prove The expression of the partial derivative is The axis of symmetry of we know that has only one zero on [0, 1], and the function f 2 (x, λ) increases monotonically first and then decreases monotonically. The minimum value of the function f 2 (x, λ) is By simple analysis, the axis of symmetry of f 2 (0, λ) is λ = − 1 2 , and the axis of symmetry of So we have proved that f 2 (x, λ) > 0. Then we will further prove α (t) ≤ 0.
We only need to prove Then, we have By analyzing the structure of the function   (7) and p = q = 1. Example 3. We consider the shape functions with exponential parameters. The expression of the shape functions are where λ, µ ≥ 1, p, q = 0, 1, 2. Next, we are about to prove that these two shape functions of Equation (8) also satisfy Equation (1). We take p = q = 2 as an example, and other cases are the same. It's obvious that Then we will further prove α (t) ≤ 0. By simple calculations, combining with the range of parameters, we have So α(t) satisfies Equation (1). The analysis of β(t) is the same. The figure of the corresponding generalized quasi cubic trigonometric Bernstein basis functions with exponential parameters defined in Equation (5) are shown in Figure 3.

Definition 4. Given control points Q i
is called a generalized quasi cubic trigonometric Bézier curve with two shape functions α(t) and β(t).
In [30], Bernstein-Bezier methods are emphasized. As the new generalized quasi cubic trigonometric Bernstein basis given in Equation (5) possesses the properties of partition of unity, nonnegativity and total positivity, we can conclude that the corresponding generalized quasi cubic trigonometric Bézier curve given in Equation (9) possesses the properties of affine invariance, convex hull and variation diminishing, which are crucial properties in the field of curve design and surface construction. With exact computation, the end-point properties of the generalized quasi cubic trigonometric Bézier curve are as follows The above end-point properties imply that for any α(t) and β(t) satisfying Equation (1), the generalized quasi cubic trigonometric Bézier curve has end-point interpolation property and Q 0 Q 1 , Q 2 Q 3 are the tangent lines of the curve at the points Q 0 and Q 3 , respectively. From these equations, we can conclude that the generalized quasi cubic trigonometric Bézier curve has some properties analogous to the similar properties of the cubic Bézier curve.
Corner cutting algorithm is an efficient and stable method of computing the generalized quasi cubic trigonometric Bézier curve and it is formed by convex combinations. In order to develop such an algorithm, we rewrite the generalized quasi cubic trigonometric Bézier curve in Equation (9) in the following matrix form From these, we can easily obtain a corner cutting algorithm for computing the generalized quasi cubic trigonometric Bernstein-Bézier curve, and more details for the Bernstein-Bézeier curve is stated in detail in [29]. While dealing with the approximation of control polygon, shape preserving approximations are considered in [31]. Analogous to these approximation methods mentioned in [31], Figure 4 gives the structure of the corner cutting algorithm and the corresponding applications of three specific examples, respectively. Moreover, for t ∈ [0, π/2], we rewrite Equation (9) as the following form From Equation (10), we can check that the shape functions α(t) and β(t) only affect the curve on the control edge Q 0 − Q 1 and Q 3 − Q 2 , respectively. Moreover, for any fixed t ∈ (0, π/2), B 0 (t; α(t)) decreases with the increase of α(t), which implies that the curve moves in the same direction of the edge Q 0 − Q 1 as α(t) increases. Inversely, when α(t) decreases, the curve moves in the opposite direction to the edge Q 0 − Q 1 . The effect of the shape function β(t) on the edge Q 3 − Q 2 is similar. When the shape functions satisfy α(t) = β(π/2 − t), with the increase or the decrease of α(t), the curve moves in the same direction or the opposite direction to the edge Q 2 − Q 1 , respectively. Furthermore, when α(t) and β(t) increase at the same time, the curve tends to the edge Q 2 − Q 1 . Thus, we can see that the shape functions α(t) and β(t) can be seen as local tension parameters. Figure 5 shows the effect of the shape parameters on the generalized quasi cubic trigonometric Bézier curves. Figure 6 shows the comparison among the new proposed generalized quasi cubic trigonometric Bézier curves (blue dashed lines, black solid lines and black dashed lines), the classical cubic rational Bézier curves (red dashed lines) and λµ-Bézier curves (red solid lines) defined by the λµ-Bernstein basis functions given in [7], constructed in the space span{1, 3t 2 − 2t 3 , (1 − t) 3 e −λt , t 3 e −µ(1−t) }. Clearly, for the same weights, the classical cubic rational Bézier curve is the farmost away from the control polygon among the three curves, and the proposed generalized quasi trigonometric Bézier curves are closer to the control polygon than the λµ-Bézier curves given in [7] under the same control points.
(a) The structure of the corner cutting algorithm.

New Generalized Quasi Cubic Trigonometric B-Spline Basis with Two Shape Functions
In this section, similar with the construction of cubic B-spline basis starting from cubic Bernstein basis, we will propose a new kind of generalized quasi cubic trigonometric B-spline basis starting from the generalized quasi cubic trigonometric Bernstein basis given in Equation (5).

New Generalized Quasi Cubic Trigonometric B-Spline Basis Functions
Given knots u 0 < u 1 < . . . < u n+4 , we denote U = (u 0 , u 1 , . . . , u n+4 ) as a knot vector. Let h j = u j+1 − u j , and t j (u) = π u − u j 2h j , j = 0, 1, . . . , n + 3. For any functions α i (t), β i (t) (i = 0, 1, . . . , n) satisfying Equation (1) and some additional conditions we want to construct the associated generalized quasi cubic trigonometric B-spline basis functions as follows where B j (t i ; α i (t), β i (t)), j = 0, 1, 2, 3 are the generalized quasi cubic trigonometric Bernstein basis functions given in Equation (5). C 2 continuity at each of knots and the property of partition of unity on the interval [u 3 , u n+1 ] are considered. From these conditions, we will calculate the expression of the coefficients a i , b i,j , c i,j , d i . Function having C 2 continuity at each of knots needs to satisfy that the function, the first derivative and the second derivative are continuous at this node, respectively. Obviously, R i (u) has C 2 continuity at u i . Then when R i (u) has C 2 continuity at u i+1 , we have When R i (u) has C 2 continuity at u i+3 , we have Further simplification When R i (u) has C 2 continuity at u i+2 , we have Further simplification The above calculation results are summarized as follows To form a partition of unity on the interval [u 3 , u n+1 ], these generalized quasi cubic trigonometric B-spline basis functions need to satisfy the equation as follows For any t ∈ [0, π/2], from ∑ 3 i=0 B i (t) = 1, we can deduce the following equations about the relation among the coefficients a i , b i,j , c i,j , d i (i = 0, 1, . . . , n) for other coefficients and bring them into the equations, we get To solve the system of equations, we get Combine the above calculations, when these generalized quasi cubic trigonometric B-spline basis functions have C 2 continuity at each of knots and form a partition of unity on the interval [u 3 , u n+1 ], we will get the expression of the coefficients a i , b i,j , c i,j , d i as follows Definition 5. Given a knot vector U, for any functions α i (t), β i (t) satisfying Equation (1) and additional Equation (11) and the coefficients a i , b i,j , c i,j , d i given in the above expressions, the expressions of Equation (12) are defined to be the associated generalized quasi cubic trigonometric B-spline basis functions.

Remark 4.
By changing the range of some parameters, Equation (7) and Equation (8) also satisfy the additional Equation (11).
For equidistant knots, the R i (u) and the knot vector U are called a uniform generalized quasi cubic trigonometric B-spline basis function and a uniform vector. For non-equidistant knots, R i (u) and U are called that of non-uniform ones, respectively. Figure 7 gives some plots of generalized quasi cubic trigonometric B-spline basis functions with different shape parameters with α(t) and β(t) given in Equation (13) and Equation (14). satisfy the conditions of Equation (1) and Equation (11) simultaneously. Moreover, with λ ≥ 1, µ ≥ 3, it is obvious that the following two shape functions α and β(t) satisfy the conditions of Equation (1) and Equation (11) simultaneously. For equidistant knots, the R i (u) and the knot vector U are called a uniform generaliz quasi cubic trigonometric B-spline basis function and a uniform vector. For non-equidista knots, R i (u) and U are called that of non-uniform ones, respectively. Figure 7 gives som plots of generalized quasi cubic trigonometric B-spline basis functions with different sha parameters with α(t) and β(t) given in Equation (13) and Equation (14).

Remark 6.
The coefficients of the B-spline stated in the Equation (12) are positive, so the gene ized quasi cubic trigonometric B-spline is also positive. Moreover, with the nonnegativity of generalized quasi cubic trigonometric Bernstein basis functions, the corresponding corner cutt algorithm is stable, then the convex combination is formed, so the overall stability is guarante Instability may only occur when caculating the trigonometric functions.

Remark 6.
The coefficients of the B-spline stated in Equation (12) are positive, so the generalized quasi cubic trigonometric B-spline is also positive. Moreover, with the nonnegativity of the generalized quasi cubic trigonometric Bernstein basis functions, the corresponding corner cutting algorithm is stable, then the convex combination is formed, so the overall stability is guaranteed. Instability may only occur when caculating the trigonometric functions.
With straightforward caculation, the following lemma is given, which is effecient when working with the partition of unity and the continuity of the generalized quasi cubic trigonometric B-spline basis functions.

Lemma 1.
For all possible i ∈ Z + , the following equations hold

Properties of the New Generalized Quasi Cubic Trigonometric B-Spline Basis Functions
In this subsection, we are about to show that the new generalized quasi cubic trigonometric B-spline basis given in Equation (12) possesses some significant properties for curves and surfaces construction, such as partition of unity (Theorem 3), nonnegativity (Theorem 4), linear independence (Theorem 5), total positivity (Theorem 6) and C 2 continuity (Theorem 7).

Proof of Theorem 3.
For u ∈ [u i , u i+1 ), i = 3, 4, . . . , n, it is easy to check that R j (u) = 0 for This indicates that the new proposed generalized quasi cubic trigonometric B-spline basis has the property of partition of unity. Theorem 4. If functions α i (t), β i (t) satisfy Equation (1) and Equation (11), then R i (u) > 0 for u i < u < u i+4 .
Proof of Theorem 5. For ξ i ∈ R (i = 0, 1, . . . , n), u ∈ [u 3 , u n+1 ], let For functions α i (t), β i (t) satisfying Equation (1) and Equation (11), direct computation gives that where i = 3, 4, . . . , n + 1. Thus, we can conclude the following linear systems of equations Notice that a i + b i0 + c i0 = 1, in order to compute the determinant of the coefficient matrix M i given by the above linear systems of equations, we add the first and the third column to the second column, respectively. For the specific expression of the determinant, we have Therefore, we can conclude that ξ i−3 = ξ i−2 = ξ i−1 = 0 for i = 3, 4, . . . , n + 1.
Proof of Theorem 6. For u ∈ [u i , u i+1 ], i = 3, 4, . . . , n, it can be easy to check that and functions α i (t), β i (t) satisfy Equation (1) and Equation (11), ) is the normalized B-basis of the space T α i ,β i . By Theorem 4.2 of [25], it is sufficient to conclude that H i is a nonsingular stochastic and totally positive matrix. For any functions α i (t), β i (t) satisfying Equations (1) and (11), for all possible i ∈ Z + , j = 0, 1, 2, 3, it is easy to check that a i , b i,j c i,j , d i > 0. In addition, from Lemma 1, we can conclude that H i is stochastic. We want to prove that H i is a totally positive matrix, we need to check the nonnegativity of all its minors. With direct computation, we have From these, we can deduce that H i is nonsingular and the nonnegativity of all its remaining minors is proved, which implies the theorem.

Theorem 7.
With a non-uniform knot vector, the generalized quasi cubic trigonometric B-spline function R i (u) with shape functions α i (t), β i (t) satisfying Equations (1) and (11) at each of the knots has C 2 continuity.
Proof of Theorem 7. Firstly, consider the continuity at the knot u i+1 . For any functions α i (t), β i (t) satisfying Equation (1) and Equation (11) From these and Lemma 1, the theorem follows at the knot u i+1 . Other knots can be analyzed in the same way.

New Generalized Quasi Cubic Trigonometric B-Spline Curves
Definition 6. Given a knot vector U and control points Q i (i = 0, 1, . . . , n) in R 2 or R 3 , then, for any functions α i (t), β i (t) satisfying Equations (1) and (11) and arbitrary real numbers n ≥ 3, u ∈ [u 3 , u n+1 ], is called a new generalized quasi cubic trigonometric B-spline curve with two local shape functions α i (t), β i (t).
Obviously, for u ∈ [u i , u i+1 ], i = 3, 4, . . . , n, the curve R(u) can be represented by the following curve segment From Theorems 3 and 4, we know that the generalized quasi cubic trigonometric Bspline basis has the properties of partition of unity and nonnegativity, so the corresponding new generalized quasi cubic trigonometric B-spline curve R(u) has affine invariance. Furthermore, for u ∈ [u i , u i+1 ], R(u) lies in the convex hull of the points Q i−3 , Q i−2 , Q i−1 , Q i . Similarly, from Theorem 6, as the generalized quasi cubic trigonometric B-spline basis has the property of total positivity, the new generalized quasi cubic trigonometric B-spline curve R(u) has variation diminishing property, which implies that R(u) is suited for the shape control. As a direct result of Theorem 7, we have the following theorem.

Theorem 8.
With a non-uniform knot vector, at each of the knots, the new generalized quasi cubic trigonometric B-spline curve R(u) has C 2 continuity. Moreover, if u i is a simple knot, then for i = 3, 4, . . . , n + 1, we have Example 4. Figure 8 shows the open spline curves with α i (t) and β i (t) given in Equation (13) with different shape parameters, which implies the shape of a rabbit. The first one is constructed by setting all λ i = 0.5, µ i = 2.4. The second one is constructed by changing one µ i from 2.4 to 6 based on the first one. The third one is constructed by changing other two µ i from 2.4 to 4 based on the second one.  (14) also satisfies the additional conditions. Figure 9 shows the open spline curves with α i (t) and β i (t) given in Equation (14) with different shape parameters, which describes the shape of a note. The first one is constructed by setting all λ i = 0.5, µ i = 2.5. The second one is constructed by changing one µ i from 2.5 to 6 based on the first one. The third one is constructed by changing another one µ i from 2.5 to 6 based on the second one. . Close curves with uniform knot vector of the shape functions with α i (t) and β i (t) given in Equation (14). Figure 10 shows the comparison among the new proposed generalized quasi cubic trigonometric B-spline curves (blue dashed lines with α i (t) and β i (t) given in Equation (13) and red dashed lines with α i (t) and β i (t) given in Equation (14)) and the trigonometric B-spline basis functions defined in [18] (blue solid lines). It is easy to conclude that for the same weights, the proposed generalized quasi trigonometric B-spline curves are closer to the control polygon than trigonometric B-spline basis functions defined in [18] under the same control points and the same parameters. What is more, the shape of the new proposed curves are more smooth. Figure 10. Comparisons among the new generalized quasi cubic trigonometric B-spline curves (blue dashed lines with α i (t) and β i (t) given in Equation (13), red dashed lines with α i (t) and β i (t) given in Equation (14)), trigonometric B-spline basis functions defined in [18] (blue solid lines) where λ i = µ i = 4.

Conclusions
We have proposed a new class of generalized quasi cubic trigonometric Bernstein basis with two shape functions, which has totally positive property and optimal shape preserving property. These two shape functions satisfy some conditions with wide range of use, and three specific examples are given. The corresponding generalized quasi cubic trigonometric Bézeir curve and the corresponding corner cutting algorithm are shown. Based on the basis, with additional conditions, a new class of generalized quasi cubic trigonometric B-spline basis functions is developed, which has partition of unity, nonnegativity, linear independence, total positivity and C 2 continuity. The corresponding generalized quasi cubic trigonometric B-spline curve and its applications are given. There are also some worthwhile work to study further, such as exploring more practical applications of those generalized basis with two shape functions. These will be our future work. Acknowledgments: The authors acknowledge the associate editor and the referees for their insightful comments and invaluable suggestions.

Conflicts of Interest:
The authors declare that they have no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:

EC
Extended Chebychev ECC Extended Complete Chebychev