(ProQuest: ... denotes non-US-ASCII text omitted.)
Academic Editor:Nicola Mastronardi
1, School of Mathematics and Statistics, Central South University, Changsha 410083, China
Received 12 September 2013; Revised 9 February 2014; Accepted 12 February 2014; 17 March 2014
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Trigonometric splines and polynomials have attracted widespread interest within computer aided geometric design (CAGD), particularly within curve design. In [1], the recurrence relation for the trigonometric B-splines of arbitrary order was established. Later, in [2], it was further shown that the trigonometric B-splines of odd order form a partition of a constant in the case of equidistant knots, and thus the associated trigonometric B-spline curve possesses the convex hull property. In [3], a family of trigonometric polynomials was introduced, which contains the trigonometric Lagrange and Bernstein polynomials. The study for the C-Bézier and C-B-spline curves can be found in [4-6]. In [7], the totally positive property of the C-B-spline was proved. In [8-11], some quadratic trigonometric polynomial splines with shape parameters were shown. In [12], a class of cubic trigonometric Bézier (T-Bézier, for short) basis with a shape parameter was proposed. In [13], the cubic T-Bézier basis was further extended to possess two shape parameters. Recently, based on the theory of envelop and topological mapping, shape analysis of the cubic T-Bézier curve with shape parameters was given in [14, 15]. For the problems of shape preserving interpolation, the cubic and quadratic T-Bézier bases show great potential applications; see [16-19]. Blossom is a powerful tool for studying Bézier-like bases and B-spline-like bases. Quasi Extended Chebyshev (QEC) space is the largest class of spaces of sufficient regularity suited for curve design in the sense that they do possess blossom. See the recently developed theory concerning blossom and QEC-space given in [20-24].
Many shape preserving properties are obtained when the normalized bases are totally positive; see [25]. For instance, the length, number of inflections, and angular variation of the curve defined by normalized totally positive bases are bounded above by those of the control polygon. The optimal normalized totally positive basis (i.e., the normalized B-basis) is the unique normalized basis of a space with optimal shape preserving properties. It is well known that the Bernstein basis is the optimal normalized totally positive basis of the space of polynomials of degree less than or equal to n ; see [26]. Corner cutting algorithm plays a vital role in CADG since it provides simple geometric constructions of curves; see [27]. A relevant example is the de Casteljau algorithm for the evaluation of a polynomial curve using the Bernstein basis. The purpose of this paper is to prove that the cubic trigonometric Bézier basis with two shape parameters λ and μ given in [13] forms an optimal normalized totally positive basis for λ , μ ∈ ( - 2,1 ] and to show a new developed corner cutting algorithm for computing the corresponding cubic trigonometric Bézier curves.
The rest of this paper is organized as follows. Section 2 gives the proof of the optimal normalized totally positive property of the cubic trigonometric Bézier basis within the general framework of QEC-space. In Section 3, a new developed corner cutting algorithm for computing the corresponding cubic trigonometric Bézier curves is shown. Conclusions are given in Section 4.
2. Cubic Trigonometric Bézier Basis Functions
2.1. Preliminaries
Here, we give the necessary background on Extended Completed Chebyshev (ECC) space and Quasi Extended Chebyshev (QEC) space. Few results are described briefly for a good understanding to this paper, and more details can be found in [20-24, 28].
Let I denote a given closed bounded interval [ a , b ] , with a < b . We call function space ( u 0 , ... , u n ) a n + 1 -dimension ECC-space in canonical form generated by positive weight functions w j ∈ C n - j ( I ) provided that [figure omitted; refer to PDF]
It is well known that ( n + 1 ) -dimension function space ( u 0 , ... , u n ) ⊂ C n ( I ) is an ECC-space on I if and only if for any k , 0 ...4; k ...4; n , any nonzero element of the subspace ( u 0 , ... , u k ) has at most k zeros (counting multiplicities). And ( n + 1 ) -dimension function space ( u 0 , ... , u n ) ⊂ C n - 1 ( I ) is a QEC-space on I if any nonzero element of the space vanishes at most n times in I , counting multiplicities as far as possible for functions in C n - 1 ( I ) , that is, up to n ; see [21-24].
A basis ( u 0 , ... , u n ) is said to be normalized on [ a , b ] if ∑ i = 0 n u i ( t ) = 1 for any t ∈ [ a , b ] . And a basis ( u 0 , ... , u n ) is said to be totally positive on [ a , b ] if, for any sequence of points a ...4; t 0 < t 1 < ... < t n ...4; b , the collocation matrix ( u j ( t i ) ) 0 ...4; i , j ...4; n is totally positive; that is, all its minors are nonnegative. For a given function space possessing totally positive basis, the optimal normalized totally positive basis (i.e., the normalized B-basis) is a normalized totally positive basis from which all other totally positive bases can be deduced by multiplying with a (regular) totally positive matrix. Such optimal normalized totally positive basis is unique and has optimal shape preserving properties (see [25]) in the sense that the curves defined by such basis best imitate the corresponding control polygons (e.g., the monotonicity and convexity of the control polygon are inherited by the curves defined by such basis).
2.2. Definition of the Cubic Trigonometric Bézier Basis
Here, we recall the definition of the cubic trigonometric Bézier basis with two shape parameters given in [13] as follows.
Definition 1.
Let λ , μ ∈ [ - 2,1 ] , for t ∈ [ 0 , π / 2 ] ; the following four functions are defined to be the cubic trigonometric Bézier (T-Bézier) basis functions, with two shape parameters λ and μ : [figure omitted; refer to PDF]
Specially, for λ = μ , this kind of T-Bézier basis with a shape parameter was also proposed in [12].
By rewriting the expressions of T 1 ( t ) and T 2 ( t ) as the following forms: [figure omitted; refer to PDF] we can easily see that the cubic T-Bézier basis (2) is normalized and is actually a basis of the following trigonometric function space: [figure omitted; refer to PDF]
2.3. Optimal Totally Positive Property of the Basis
Firstly, we want to show that for any λ , μ ∈ ( - 2,1 ] , the following space [figure omitted; refer to PDF] is a 3-dimensional QEC-space on [ 0 , π / 2 ] .
Theorem 2.
For any real numbers λ , μ ∈ ( - 2,1 ] , the space D T λ , μ is a 3-dimensional QEC-space on [ 0 , π / 2 ] .
Proof.
For any ξ i ∈ ... , λ , μ ∈ ( - 2,1 ] , and t ∈ [ 0 , π / 2 ] , we consider a linear combination: [figure omitted; refer to PDF]
For t = 0 , from (6), we can immediately obtain ξ 1 = 0 . Similarly, for t = π / 2 , from (6), we have ξ 2 = 0 . And finally, we have ξ 0 = 0 . Thus the space D T λ , μ is a 3-dimensional space.
Now, we want to prove that the space D T λ , μ is a 3-dimensional ECC-space in ( 0 , π / 2 ) . For any t ∈ [ a , b ] ⊂ ( 0 , π / 2 ) , let [figure omitted; refer to PDF] For any λ , μ ∈ ( - 2 , 1 ] , it is obvious that u ( t ) > 0 and v ( t ) > 0 . By directly computing, we have [figure omitted; refer to PDF] Therefore, for the Wronskian of u and v , we have [figure omitted; refer to PDF] For t ∈ [ a , b ] , we define the following weight functions: [figure omitted; refer to PDF] where A , B , and C are three arbitrary positive real numbers. Obviously, these weight functions w i ( t ) ( i = 0,1 , 2 ) are bounded, positive, and C ∞ on [ a , b ] . Consider the following ECC-space defined by the weight functions w i ( t ) ( i = 0,1 , 2 ) : [figure omitted; refer to PDF] After some simple computations, we can see that these functions u i , i = 0,1 , 2 are in fact some linear combinations of the three functions 2 sin t cos ... t , - cos ... t ( 1 - sin t ) ( 2 + λ - 3 λ sin t ) , and sin t ( 1 - cos ... t ) ( 2 + μ - 3 μ cos ... t ) , which indicates that the space D T λ , μ is a ECC-space on [ a , b ] . Since [ a , b ] is arbitrary subinterval of ( 0 , π / 2 ) , we can conclude that the space D T λ , μ is an ECC-space in ( 0 , π / 2 ) . Now, we want to further prove that the space D T λ , μ is also a QEC-space on [ 0 , π / 2 ] . For this purpose, we need to prove that any nonzero element of the space D T λ , μ has at most 2 zeroes on [ 0 , π / 2 ] (counting multiplicities as far as possible up to 2). Consider any nonzero function: [figure omitted; refer to PDF] where t ∈ [ 0 , π / 2 ] . Since the space D T λ , μ is an ECC-space in ( 0 , π / 2 ) , the function F ( t ) has at most two zeroes in ( 0 , π / 2 ) . Let us consider a case that F ( t ) vanishes at 0 ; then, we have C 1 = 0 . For this case, if C 2 = 0 , then F ( t ) has a singular zero at 0 and a singular zero at π / 2 . If C 0 = 0 , it can be easily checked that 0 is a double zero of F ( t ) (counting multiplicities as far as possible up to 2). If C 0 C 2 > 0 , F ( t ) vanishes exactly one time at 0 and it does not vanish on ( 0 , π / 2 ] . If C 0 C 2 < 0 , F ( t ) vanishes exactly one time at 0 and it does not vanish at π / 2 . Moreover, consider the following function: [figure omitted; refer to PDF] Direct computation gives that G [variant prime] ( t ) = sin t K ( t ) , where K ( t ) = - 2 C 0 + C 2 ( 2 + 4 μ - 6 μ cos ... t ) . Obviously, K ( t ) is a monotonic function on [ 0 , π / 2 ] , and thus K ( t ) has at most one zero in ( 0 , π / 2 ) , which implies that G [variant prime] ( t ) also has at most one zero in ( 0 , π / 2 ) . From these together with G ( 0 ) G ( π / 2 ) = 2 ( 2 + μ ) C 0 C 2 < 0 , we can see that G ( t ) has exactly one zero in ( 0 , π / 2 ) , and thus we can immediately conclude that F ( t ) = sin t G ( t ) (notice that C 1 = 0 for the current case) has exactly one zero in ( 0 , π / 2 ) too. Similarly, for the case that F ( t ) vanishes at π / 2 , we can also deduce that the function F ( t ) has at most 2 zeroes on [ 0 , π / 2 ] (counting multiplicities as far as possible up to 2). To summarize, the space D T λ , μ is a QEC-space on [ 0 , π / 2 ] .
Remark 3.
By Theorem 3.1 of [22], we can see that blossom exists in the trigonometric function space T λ , μ is equivalent to that the space D T λ , μ is a QEC-space. For any λ , μ ∈ ( - 2,1 ] , since the space D T λ , μ is a QEC-space on [ 0 , π / 2 ] , blossom exists in T λ , μ . By Theorem 2.13 of [22], we can also see that the space T λ , μ possesses a normalized Quasi Bernstein-like basis on [ 0 , π / 2 ] for λ , μ ∈ ( - 2,1 ] .
It is worth mentioning that for the special case λ = - 2 or μ = - 2 , F ( t ) may have three distinct zeros on [ 0 , π / 2 ] . For example, for μ = - 2 , C 0 = 3 , C 1 = 0 , and C 2 = - 2 , F ( t ) will vanish at 0 , π / 3 and π / 2 , respectively, which implies that D T λ , μ is not a QEC-space on [ 0 , π / 2 ] for λ = - 2 or μ = - 2 . Thus, from the blossom point of view, T λ , μ is not suited for curve design for the case that λ = - 2 or μ = - 2 .
Theorem 4.
For any λ , μ ∈ ( - 2,1 ] , the cubic T-Bézier basis given in (2) is the optimal normalized totally positive basis of the space T λ , μ .
Proof.
For any λ , μ ∈ ( - 2,1 ] , from the expressions of the four cubic T-Bézier basis functions given in (2), we can easily check the following end-point properties of the basis.
(i) T 0 ( 0 ) = 1 , and T 0 ( t ) vanishes 3 times at π / 2 (counting multiplicities as far as possible up to 3 ).
(ii) T 3 ( π / 2 ) = 1 , and T 3 ( t ) vanishes 3 times at 0 (counting multiplicities as far as possible up to 3 ).
(iii): For i = 1,2 , T i ( t ) vanishes exactly i times at 0 and exactly 3 - i times at π / 2 .
In addition, T i ( t ) ( i = 0,1 , 2,3 ) is strictly positive in ( 0 , π / 2 ) . Thus, by Definition 2.10 of [22], we can see that the cubic T-Bézier basis (2) is precisely the Quasi Bernstein-like basis of the space T λ , μ . And by Theorem 2.18 of [22], we can further conclude that the cubic T-Bézier basis (2) is exactly the optimal normalized totally positive basis of the space T λ , μ .
3. Cubic Trigonometric Bézier Curve
Definition 5.
Given control points P i ( i = 0,1 , 2,3 ) in ... 2 or ... 3 , then [figure omitted; refer to PDF] is called a cubic trigonometric Bézier (T-Bézier) curve with two shape parameters λ and μ .
For any λ , μ ∈ [ - 2,1 ] , from [13], we can see that the cubic T-Bézier basis (2) has the properties of partition of unity and nonnegativity, which implies that the corresponding cubic T-Bézier curve (14) has affine invariance and convex hull property. And for any λ , μ ∈ ( - 2,1 ] , by Theorem 4, we can see that the cubic T-Bézier basis (2) has the optimal totally positive property, which implies that the corresponding cubic T-Bézier curve (14) has the crucial property of variation diminishing; that is, no plane intersects a cubic T-Bézier curve more often than it intersects the corresponding control polygon.
Now we want to develop a new corner cutting algorithm for computing the cubic T-Bézier curve (14). For this purpose, for any t ∈ [ 0 , π / 2 ] , let [figure omitted; refer to PDF] Then, we can rewrite the expression of the cubic T-Bézier curve (14) as follows: [figure omitted; refer to PDF] Furthermore, by setting [figure omitted; refer to PDF] we have [figure omitted; refer to PDF] The expressions (16) and (18) describe a corner cutting algorithm for computing the cubic T-Bézier curve (14). See Figure 1 for an illustration of this new developed algorithm.
Corner cutting algorithm.
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF]
(d) [figure omitted; refer to PDF]
4. Conclusion
The cubic trigonometric Bézier basis with two shape parameters λ and μ given in [13] forms an optimal normalized totally positive basis of the space T λ , μ for λ , μ ∈ ( - 2,1 ] . However, the basis is not suited for curve design for λ = - 2 or μ = - 2 from the blossom point of view. The new developed corner cutting algorithm is an efficient and stable process for computing the corresponding cubic trigonometric Bézier curves.
Acknowledgments
The authors would like to thank the anonymous reviewers for their valuable remarks for improvements. The research is supported by the National Natural Science Foundation of China under Grant nos. 60970097 and 11271376 and Graduate Students Scientific Research Innovation Project of Hunan Province, Grant no. CX2012B111.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
[1] T. Lyche, R. Winther, "A stable recurrence relation for trigonometric B -splines," Journal of Approximation Theory , vol. 25, no. 3, pp. 266-279, 1979.
[2] G. Walz, "Identities for trigonometric B -splines with an application to curve design," BIT Numerical Mathematics , vol. 37, no. 1, pp. 189-201, 1997.
[3] G. Walz, "Trigonometric Bézier and Stancu polynomials over intervals and triangles," Computer Aided Geometric Design , vol. 14, no. 4, pp. 393-397, 1997.
[4] J. W. Zhang, "C-curves: an extension of cubic curves," Computer Aided Geometric Design , vol. 13, no. 3, pp. 199-217, 1996.
[5] J. W. Zhang, "Two different forms of C - B -splines," Computer Aided Geometric Design , vol. 14, no. 1, pp. 31-41, 1997.
[6] J. W. Zhang, F. L. Krause, H. Y. Zhang, "Unifying C-curves and H-curves by extending the calculation to complex numbers," Computer Aided Geometric Design , vol. 22, no. 9, pp. 865-883, 2005.
[7] E. Mainar, J. M. Peña, "A basis of C -Bézier splines with optimal properties," Computer Aided Geometric Design , vol. 19, no. 4, pp. 291-295, 2002.
[8] X. L. Han, "Quadratic trigonometric polynomial curves with a shape parameter," Computer Aided Geometric Design , vol. 19, no. 7, pp. 503-512, 2002.
[9] X. L. Han, "Piecewise quadratic trigonometric polynomial curves," Mathematics of Computation , vol. 72, no. 243, pp. 1369-1377, 2003.
[10] X. L. Han, " C 2 quadratic trigonometric polynomial curves with local bias," Journal of Computational and Applied Mathematics , vol. 180, no. 1, pp. 161-172, 2005.
[11] X. Han, "Quadratic trigonometric polynomial curves concerning local control," Applied Numerical Mathematics , vol. 56, no. 1, pp. 105-115, 2006.
[12] X. L. Han, "Cubic trigonometric polynomial curves with a shape parameter," Computer Aided Geometric Design , vol. 21, no. 6, pp. 535-548, 2004.
[13] X. A. Han, Y. Ma, X. Huang, "The cubic trigonometric Bézier curve with two shape parameters," Applied Mathematics Letters , vol. 22, no. 2, pp. 226-231, 2009.
[14] X. A. Han, X. L. Huang, Y. C. Ma, "Shape analysis of cubic trigonometric Bézier curves with a shape parameter," Applied Mathematics and Computation , vol. 217, no. 6, pp. 2527-2533, 2010.
[15] R. J. Wu, G. H. Peng, "Shape analysis of planar trigonometric Bézier curves with two shape parameters," International Journal of Computer Science , vol. 10, no. 3, pp. 441-447, 2013.
[16] F. Ibraheem, M. Hussain, M. Z. Hussain, A. A. Bhatti, "Positive data visualization using trigonometric function," Journal of Applied Mathematics , vol. 2012, 2012.
[17] U. Bashir, J. Md. Ali, "Data visualization using rational trigonometric spline," Journal of Applied Mathematics , vol. 2013, 2013.
[18] M. Hussain, S. Saleem, " C 1 rational quadratic trigonometric spline," Egyptian Informatics Journal , vol. 14, no. 3, pp. 211-220, 2013.
[19] M. Z. Hussain, M. Hussain, A. Waseem, "Shape-preserving trigonometric functions," Computational and Applied Mathematics , 2013.
[20] M. L. Mazure, "Quasi-Chebyshev splines with connection matrices: application to variable degree polynomial splines," Computer Aided Geometric Design , vol. 18, no. 3, pp. 287-298, 2001.
[21] M. L. Mazure, "On dimension elevation in quasi extended Chebyshev spaces," Numerische Mathematik , vol. 109, no. 3, pp. 459-475, 2008.
[22] M. L. Mazure, "Which spaces for design?," Numerische Mathematik , vol. 110, no. 3, pp. 357-392, 2008.
[23] M. L. Mazure, "On a general new class of quasi Chebyshevian splines," Numerical Algorithms , vol. 58, no. 3, pp. 399-438, 2011.
[24] M. L. Mazure, "Quasi extended Chebyshev spaces and weight functions," Numerische Mathematik , vol. 118, no. 1, pp. 79-108, 2011.
[25] J. M. Carnicer, E. Mainar, J. M. Peña, "Critical length for design purposes and extended Chebyshev spaces," Constructive Approximation , vol. 20, no. 1, pp. 55-71, 2003.
[26] J. M. Carnicer, J. M. Peña, "Shape preserving representations and optimality of the Bernstein basis," Advances in Computational Mathematics , vol. 1, no. 2, pp. 173-196, 1993.
[27] E. Mainar, J. M. Peña, "Corner cutting algorithms associated with optimal shape preserving representations," Computer Aided Geometric Design , vol. 16, no. 9, pp. 883-906, 1999.
[28] L. L. Schumaker Spline Functions: Basic Theory , Cambridge University Press, Cambridge, UK, 2007., 3rd.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2014 Xuli Han and Yuanpeng Zhu. Xuli Han et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
Within the general framework of Quasi Extended Chebyshev space, we prove that the cubic trigonometric Bézier basis with two shape parameters λ and μ given in Han et al. (2009) forms an optimal normalized totally positive basis for λ, μ ∈(- 2,1 ] . Moreover, we show that for λ =- 2 or μ =- 2 the basis is not suited for curve design from the blossom point of view. In order to compute the corresponding cubic trigonometric Bézier curves stably and efficiently, we also develop a new corner cutting algorithm.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer