A New Quasi Cubic Rational System with Two Parameters

: The purpose of this article is to develop a new system for the construction of curves and surfaces. Making the new system not only has excellent properties of the orthodox B é zier and the B-spline method but also has practical properties such as variation diminishing and local shape adjustability. First, a new set of the quasi-cubic rational (QCR) system with two parameters is given, which is veriﬁed on an optimal normalized totally positive system (B-system). The related QCR B é zier curve is deﬁned, and the de Casteljau-type algorithm are given. Next, a group of non-uniform QCR B-spline system is shown based on the linear combination of the proposed QCR system, the relative properties of the B-spline system are analyzed. Then, the deﬁnition and properties of non-uniform QCR B-spline curves are discussed in detail. Finally, the proposed QCR system is extended to the triangular domain, which is called the quasi-cubic rational Bernstein-B é zier (QCR-BB) system, and its related deﬁnition and properties of patches are given at length. The experimental image obtained by using MATLAB shows that the newly constructed system has excellent properties such as symmetry, totally positive, and C 2 continuity, and its corresponding curve has the properties of local shape adjustability and C 2 continuity. These extended systems in the extended triangular domain have symmetry, linear independence, etc. Hence, the methods in this article are suitable for the modeling design of complex curves and surfaces.


Introduction
Bézier and B-spline methods are mainstream methods for computer-aided geometric design (CAGD). The Bézier and B-Spline method have good properties such as symmetry, continuity, convexity, and geometric invariance. Both methods are generally used in geometric design because they have some good properties. The rapid development of the modern geometric industry has led to the difficulty that traditional Bézier or B-spline curves and surfaces could not meet the requirements for geometric designers. Thus, several rationalizing Bézier curves have been proposed, but they have the problem of gradualness. Moreover, improper use of the weight factor would destroy geometric models which have been designed [1][2][3][4][5]. To overcome the above problems, scholars have proposed several Bernstein and B-spline systems with parameters [6][7][8][9][10][11][12][13][14].
The curves and surfaces constructed in the algebraic polynomial space have the advantages of simple and intuitive structure, few calculations, and low computational complexity. To increase the ability of geometric modeling of the traditional Bézier and B-spline methods, researchers have studied the key techniques of Bézier and B-spline curves in the algebraic polynomial space with parameters, such as shape adjustment, joining, and rotary surfaces. In reference [15], a set of Bernstein systems is based on the algebraic polynomial space Span 1, 3t 2 − 2t 3 , (1 − λt)(1 − t) 3 , (1 − λ + λt)t 3 .
When α = β, the QCR-B system has a similar form as the literature [15]. In particular, when α = β = 0, the QCR-B system would degenerate to the classical cubic Bernstein system. The definition and properties (such as nonnegativity, symmetry, and end-point properties) of the QCR-B system had already been discussed in [15] at length. However, the totally positive property is missed. Thus, we could use the theoretical knowledge of the ECC and QEC space to discuss the totally positive property of the QCR-B system.
Next, we discuss that the QCR-B system will form a set of B-systems in space T α,β . For any α, β ∈ [0, 1], t ∈ [0, 1], the related mother-functions of space From Theorem 3.1 in reference [22], what we should do is to prove that the differential space

Proof. Let
By directly differentiating on u(t) and v(t), it follows that: Therefore, for any t ∈ [a, b], the related Wronskian could be expressed as follows: Regarding any t ∈ [a, b] ⊂ (0, 1), define the weight functions as follows: where λ i > 0(i = 1, 2, 3). Therefore, for t ∈ [a, b] ⊂ (0, 1), we have: Then, we consider the following ECC spaces from Definition 1: The formulas u 0 (t), u 1 (t), u 2 (t) that could be easily verified are linear combinations of the differential space Thus, the differentiating space DT α,β is an ECC space on the closed interval [a, b]. Moreover [a, b] is any subset of (0, 1), so it is an ECC space on interval (0, 1). Therefore, any linear combinations of DT α,β spaces have at most two zeros on interval (0, 1).
It is further proved that the differential space is the upper QEC space. Next, let us prove that the differentiating space DT α,β is a QEC space on [0, 1]. In this case, any linear combinations of DT α,β have at most two zeros on [0, 1]. The linear function can be considered as: Given that DT α,β is an ECC space on (0, 1). Hence, F(t) has at the utmost two zeros on (0, 1). If F(t) has one zero at the point t = 0, then µ 1 = 0. Under this condition: If µ 2 = 0, then F(t) has two zeros t = 0, t = 1.
Firstly, we prove that T i (t) is a linear independence. For ξ i ∈ R(i = 0, 1, 2, 3), the following linear combination: By directly differentiating on both sides of the above expression, we have: When t = 0, from Equations (2) and (3), we can easily get: Thus, we can get ξ 0 = ξ 1 = 0. When t = 1, we can get ξ 2 = ξ 3 = 0 in the same way. Thus, the system T i (t) is linear independence.
For convenience, we define the corresponding QCR-B system as T i (α, β; t). The images of the QCR-B system with different parameters are shown in Figure 1.

Definition 4. For control points
as a QCR-Bézier curve with two parameters β α, .
, the QCR-Bézier curve in reference [12] has proved to have convex hull, convexity, and end-point properties. Given that the system has the totally positive property, the related QCR-Bézier curve was variation diminishing as well.
Next, we give the de Casteljau-type algorithm of the QCR-Bézier curve, which is stable and efficient for computing each point of the QCR-Bézier curve. For any Then, Equation (4) can be rewritten as follows: Moreover, let
When α, β ∈ [0, 1], the QCR-Bézier curve in reference [12] has proved to have convex hull, convexity, and end-point properties. Given that the system has the totally positive property, the related QCR-Bézier curve was variation diminishing as well.
Next, we give the de Casteljau-type algorithm of the QCR-Bézier curve, which is stable and efficient for computing each point of the QCR-Bézier curve. For any t ∈ [0, 1], let Then, Equation (4) can be rewritten as follows: Moreover, let (5) can be rewritten as follows: Equations (5) and (6) are stable and efficient de Casteljau-type algorithms for generating the QCR-Bézier curve. For convenience, we also can rewrite Equation (4) as Equation (7). Figure 2a shows the detailed process of the algorithm, and Figure 2b-d shows the curve generated by the de Casteljau-type algorithm when taking different parametric values.
type algorithm when taking different parametric values.
( ) (a) The detailed process of the algorithm (b) QCR-Bézier curve when

Structure of the Non-Uniform QCR-B Spline System
Given any knot sequence , the non-uniform QCR-B spline system has the following form:

Structure of the Non-Uniform QCR-B Spline System
Given any knot sequence u 0 < u 1 < · · · < u n+4 , let us write it down as U = (u 0 , u 1 , · · · , u n+4 ). For α i , β i , the non-uniform QCR-B spline system has the following form: where By easily calculating, the related coefficient values are shown as follows: (8) as a set of the QCR-B spline system concerning parameters α and β.
The case of vector U in equidistant node vectors is different from the case of vector U in non-equidistant node vectors. In the former case, the QCR-B spline system is referred to as a uniform QCR-B spline system. Under the latter condition, the QCR-B spline system is referred to as a non-uniform QCR-B spline system. Figure 3 shows the images of the uniform QCR-B-spline system images under different parameters.  Direct calculations yield the following lemma, which will be very useful for subsequent discussion.

Lemma 1. For any
Z i ∈ , the following equations are true:  Direct calculations yield the following lemma, which will be very useful for subsequent discussion.
could form a set of normalized totally positive system of T α i ,β i .
Proof. For u ∈ [u i , u i+1 ], i = 3, 4, · · · , n, we have: where, and t i (u) = (u − u i )/h i , i = 0, 1, 2, · · · , n + 3. From Theorem 2, (T 0 (t i; α i ), T 1 (t i; α i ), T 2 (t i; β i ), T 3 (t i; β i )) is a B-system of space T α i ,β i . By Theorem 4.2 in reference [27], H i is a non-singular random totally positive matrix. Thus, the system (B i−3 (u), B i−2 (u), B i−1 (u), B i (u)) is from a totally positive system in space For any α i , β i , a i , b i,j , c i,j , d i > 0 is evident. From Lemma 1, we can get that the sum of each row of the matric H i is 1. That is, H i is a stochastic matrix. Moreover, by directly calculating, it follows that: Thus, the matrix H i is a non-singular random totally positive matrix.
Proof. For any α i , β i ∈ [0, 1], it could conclude: From Lemma 1 and the above discussion, we have:

Definition and Properties
Definition 6. When the knot vector U and control points P i (i = 0, 1, . . . , n) of non-uniform knots vector are given. A linear combination of the non-uniform QCR-B spline system and control front can generate the non-uniform QCR-B spline curve. For α i , β i ∈ [0, 1], i ≥ 3, u ∈ [u 3 , u n+1 ], we call: is the non-uniform QCR-B spline curve with parameters α i and β i .
For u ∈ [u i , u i+1 ], the non-uniform QCR-B spline segment of curve Q i (u) can be expressed as follows: From the properties of the non-uniform QCR-B spline system (1) and (2), for u ∈ [u i , u i+1 ], the non-uniform QCR-B spline curve Q i (u) has affine invariance. Furthermore, from the properties of the non-uniform QCR-B spline system (6), the non-uniform QCR-B spline curve Q i (u) has variation diminishing. Therefore, the non-uniform QCR-B spline curve keeps all the essential properties of the traditional B-spline. Theorem 3. Given a non-uniform knot vector U, for any α i , β i the non-uniform QCR-B spline curve has C 2 continuity at each control point.

Local Adjustment
Since the QCR-B spline curve Q(u) has parameters α i , β i . The shape of the Q(u) could be locally modified by changing the value of α i , β i when keeping the control polygon fixed. From Equation (13), when u ∈ [u i , u i+1 ], the shape of the Q(u) could be affected by α i , β i . Moreover, α i can control the shape of the four segments on [u i−2 , u i+2 ], and β i can control the shape of the four segments on [u i−1 , u i+3 ].
Trends of the QCR-B spline curve Q(u) can be predicted by altering the values of the α i , β i . When the value of α i , β i increases, the coefficient of P i−3 , P i decreases and the coefficient of the P i−2 , P i−1 increases. Therefore, when α i , β i increases synchronously, the curve Q(u) approaches to the edge P i−2 P i−1 . As α i increases, the curve Q(u) approaches to P i−2 . As β i increases, the curve Q(u) approaches to P i−1 . Figure 5 shows the QCR-B spline curve. Figure 5a displays the case where the parameters α i , β i are different. Figure 5b displays that the curve has a locally adjustable property.  Locally shape adjustable QCR-B spline curve.

QCR-BB System over Triangular Domain
, the ten polynomials are called a QCR-BB system concerning three parameters β α , , and γ over the triangular domain D: 3 3
Proof. We will proof properties (2) and (5) as follows, and the remaining properties are easy to be verified.
combination of the non-uniform QCR-BB system and the control mesh can generate a QCR-BB patch. We call: Figure 6. The images of QCR-BB system.
(4) Boundary property: When w = 0, R(u, v, w) degenerated into a QCR-Bézier curve with two parameters α, β. When u = 0, R(u, v, w) degenerated into a QCR-Bézier curve with β, γ. When v = 0, R(u, v, w) degenerated into a QCR-Bézier curve with α, γ. With the values of α, β, and γ increasing, the QCR-BB patch will be approached to the control mesh. Hence, the parameters α, β, and γ have a tension effect. (5) Shape adjustability: The shape of the R(u, v, w) can be turned-up by modifying the value of the α, β, and γ when the control mesh is stabled. With the values of α, β, and γ increasing, the R(u, v, w) would approach to the control mesh. Hence, it is easy for us to get the parameters α, β, and γ to have a tension effect. According to the boundary property of the R(u, v, w), each boundary curve R(u, 0, w), R(0, v, w), and R(u, v, 0) only have two related parameters. Thus, changing a parameter can only affect the shape of two boundary curves. Figure 7 shows the effects of different parameter values on the QCR-BB patch when the control mesh is fixed.

De Casteljau-Type Algorithm
Next, given is a high-efficiency de Casteljau-type algorithm for generating QCR-BB patches. For , it can: From the above discussion, the following de Casteljau-type algorithm can be obtained:
From the above discussions, the following theorem can be obtained.

Theorem 4.
For any α i , β, γ ∈ [0, 1], i = 1, 2, if the control points satisfy Equation (15) and Equation (19), simultaneously, then the two QCR-BB patches are G 1 continuity. Figure 8 shows the G 1 continuity of two QCR-BB patches with different parameters. , if the control points satisfy Equation (15) and Equation (19), simultaneously, then the two QCR-BB patches are 1 G continuity. Figure 8 shows the 1 G continuity of two QCR-BB patches with different parameters.

Conclusions
Within the theoretical framework of the QEC space, a set of B-systems with two parameters was constructed. The novel system overcomes the shortcomings of the traditional literature that only discusses curve flexibility but neglects important properties (variation diminishing and totally positive). Based on the new system, the related QCR Bézier curves are constructed. A large experimental show of the non-uniform B-spline curves based on the B-system not only retains the excellent properties of the traditional B-spline method but also practical properties such as local shape adjustability and 2 C continuity. Moreover, the new system is also extended to the triangular domain, and the QCR-BB patches are defined, which is proved to have flexible shape adjustability. Under certain conditions, we give the conditions of the 1 G continuity joining patches of two QCR-BB patches. Thus, the curve and patch constructed in this work have a strong ability to geometric modeling. Our future work will pay attention to the precise quantitative analysis of the influence of parameters on a B-spline curve and the analysis of the cusps and inflection points of the proposed curve is necessary to meet the needs of complex modeling design.

Conclusions
Within the theoretical framework of the QEC space, a set of B-systems with two parameters was constructed. The novel system overcomes the shortcomings of the traditional literature that only discusses curve flexibility but neglects important properties (variation diminishing and totally positive). Based on the new system, the related QCR Bézier curves are constructed. A large experimental show of the non-uniform B-spline curves based on the B-system not only retains the excellent properties of the traditional B-spline method but also practical properties such as local shape adjustability and C 2 continuity. Moreover, the new system is also extended to the triangular domain, and the QCR-BB patches are defined, which is proved to have flexible shape adjustability. Under certain conditions, we give the conditions of the G 1 continuity joining patches of two QCR-BB patches. Thus, the curve and patch constructed in this work have a strong ability to geometric modeling. Our future work will pay attention to the precise quantitative analysis of the influence of parameters on a B-spline curve and the analysis of the cusps and inflection points of the proposed curve is necessary to meet the needs of complex modeling design.