Next Article in Journal
Exiting Inflation with a Smooth Scale Factor
Next Article in Special Issue
Curve-Surface Pairs on Embedded Surfaces and Involute D-Scroll of the Curve-Surface Pair in E3
Previous Article in Journal
Residual Tsallis Entropy and Record Values: Some New Insights
Previous Article in Special Issue
Proposed Theorems on the Lifts of Kenmotsu Manifolds Admitting a Non-Symmetric Non-Metric Connection (NSNMC) in the Tangent Bundle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Shape Preserving Class of Two-Frequency Trigonometric B-Spline Curves

1
School of Mathematics, National University of Colombia, Medellín Campus, Medellín 4309511, Colombia
2
Department of Applied Mathematics, University Research Institute of Mathematics and Its Applications (IUMA), University of Zaragoza, 50009 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(11), 2041; https://doi.org/10.3390/sym15112041
Submission received: 9 October 2023 / Revised: 4 November 2023 / Accepted: 8 November 2023 / Published: 10 November 2023
(This article belongs to the Special Issue Symmetry/Asymmetry: Differential Geometry and Its Applications)

Abstract

:
This paper proposes a new approach to define two frequency trigonometric spline curves with interesting shape preserving properties. This construction requires the normalized B-basis of the space U 4 ( I α ) = span { 1 , cos t , sin t , cos 2 t , sin 2 t } defined on compact intervals I α = [ 0 , α ] , where α is a global shape parameter. It will be shown that the normalized B-basis can be regarded as the equivalent in the trigonometric space U 4 ( I α ) to the Bernstein polynomial basis and shares its well-known symmetry properties. In fact, the normalized B-basis functions converge to the Bernstein polynomials as α 0 . As a consequence, the convergence of the obtained piecewise trigonometric curves to uniform quartic B-Spline curves will be also shown. The proposed trigonometric spline curves can be used for CAM design, trajectory-generation, data fitting on the sphere and even to define new algebraic-trigonometric Pythagorean-Hodograph curves and their piecewise counterparts allowing the resolution of C ( 3 Hermite interpolation problems.

1. Introduction

The definition of new spaces with more flexibility than polynomials, but with the same nice structural properties, such as variation diminishing, containment in the convex-hull, affine invariance, or tangency to the control-polygon at the endpoints, is an interesting trend in Computer-Aided Geometric Design (CAGD).
For a suitable basis ( b 0 , , b n ) of a given space of functions defined on I R , γ ( t ) = i = 0 n P i b i ( t ) , t I , provides a parametric representation of curves, where the coefficients are points in a given R d determining a polygon P 0 P n , which is called the control polygon of γ .
A system ( b 0 , , b n ) of functions defined on I is normalized if i = 0 n b i ( t ) = 1 for all t I . On the other hand, a totally positive (TP) basis is a basis whose collocation matrices have nonnegative minors. We say that a basis provides a shape preserving representation if the shape of any parametric curve γ imitates the shape of its control polygon. Normalized and totally positive (NTP) bases provide shape preserving representations (cf. [1,2]). In [3], it was proved that a space of functions with NTP bases always has an optimal shape preserving basis that it is called the normalized B-basis and satisfies that the matrix of change of basis of any NTP basis with respect to the normalized B-basis is totally positive and stochastic. Roughly speaking, any parametric curve more faithfully imitates the shape of its control polygon with respect to the normalized B-basis than using other shape preserving representation.
The Bernstein polynomials of a given degree n on a compact interval [ a , b ] are defined by
b i n ( t ) = n i ( t a ) i ( b t ) n i / ( b a ) n , t [ a , b ] , i = 0 , , n .
It is well-known that the Bernstein basis ( b 0 n , , b n n ) is the normalized B-basis of the space of polynomials of degree at most n on the considered interval. The B-spline basis is the normalized B-basis in the case of the space of spline polynomials.
Trigonometric functions are usually considered for the representation of closed curves and periodical functions. The space of n order trigonometric polynomials,
U 2 n : = span { 1 , cos t , sin t , cos ( 2 t ) , sin ( 2 t ) , , cos ( n t ) , sin ( n t ) } ,
is a classical space of functions with applications in CAGD in order to represent approximately parametric curves arising in civil engineering and robotics, among other fields. Traditionally these curves were approximated by polynomial and rational functions.
Due to the oscillation properties of trigonometric functions, the representation of curves using trigonometric spaces is not shape preserving for very long intervals. However, using a domain of a length α < π , the space U 2 n admits a normalized B-basis and solves several shortcomings of the algebraic polynomials. Consequently, it provides an alternative to the rational model for several purposes (see [4]). In the literature one can find many papers introducing shape parameters for a more flexible design of trigonometric curves (see [5,6,7,8,9,10,11,12]). For trigonometric polynomials in U 2 n , the length α of the parameter domain can also be considered as a shape parameter with a tension-like effect in the parametric curves generated by the corresponding normalized B-bases. Trigonometric spline functions were introduced in [13] (see also [14]) and the recurrence relation for the trigonometric B-splines of arbitrary order was obtained in [15]. In [16,17], interested readers can find a novel technique based on the collocation finite element method for the numerical resolution of the wave equation in which trigonometric cubic B-splines are used as approximate functions.
Spline spaces provide more flexibility for curve design and then common curve and surface representations in CAGD use piecewise polynomials. Integral recurrence formulae for B-splines have been often used in the past. The definition of B-spline functions as a divided difference of a truncated power function and the Hermite–Gennochi formula lead to integral recursions. However, trigonometric polynomial spaces do not admit an integral construction. For this reason, this paper introduces NTP bases of U 2 n , such that their derivatives at the ends of the interval of [ 0 , α ] allow us the definition of regular piecewise trigonometric curves. In [18], one-frequency trigonometric spline curves are defined in the space U 2 . This paper proposes and analyzes piecewise trigonometric curves in U 4 providing two different frequencies for the design.
Polynomial Pythagorean-Hodograph (PH) curves have been widely analyzed (see [19,20,21,22,23,24,25] and their references). Polynomial PH curves are defined with Bernstein bases and possess a closed-form polynomial representation of their arc lengths, as well as an exact rational parameterization of their offset curves. New Pythagorean-Hodograph (PH) B-spline curves were proposed in [21]. The results in this paper will be considered in the near future to define new algebraic-trigonometric PH curves and their piecewise counterparts for the resolution of C ( 3 Hermite interpolation problems. Moreover, the proposed trigonometric spline curves can also be used for CAM design, trajectory generation or data fitting on the sphere where the use of trigonometric splines provide better results than the conventional polynomial counterpart (see Chapter 12 of [1]).
The paper is organized as follows. Section 2 provides the normalized B-basis of the space U 4 ( I α ) defined on a given interval with a length less than π . A corner cutting algorithm for the evaluation and subdivision of curves generated in U 4 ( I α ) is also provided. This algorithm allows us to generate curves through a control polygon in a similar way to the Bézier case [26,27]. In Section 3, we introduce NTP bases of U 4 defined on appropriate intervals [ 0 , α ] , and such that their derivatives at the ends of the interval allow us the definition of regular spline spaces and their bases called normalized T 4 B-spline basis. The matrices relating the normalized B-bases of U 4 and the spline spaces are also derived and their total positivity analyzed. Piecewise trigonometric curves with nice properties are also defined. Section 4 shows the convergence of the introduced trigonometric curves to polynomial B-spline curves as the parameter α 0 . Finally, Section 5 summarizes the conclusions and future work.

2. The Normalized B-Basis of the Space of Two Frequency Trigonometric Functions

Trigonometric polynomial spaces U 2 n are invariant under translations and reflexions and, for this reason, they are usually analyzed on compact intervals I α : = [ 0 , α ] .
We shall consider the design of parametric trigonometric curves whose components are functions in the space
U 4 ( I α ) = span { 1 , cos ( t ) , sin t , cos ( 2 t ) , sin ( 2 t ) } , t I α .
In Section 4.2 of [28], the normalized B-basis ( B 0 4 , , B 4 4 ) of U 4 ( I α ) is given by:
B 0 4 ( t ) : = ( 1 cos ( α t ) ) 2 ( 1 cos α ) 2 , B 1 4 ( t ) : = 2 ( 1 cos ( α t ) ) ( cos t + cos ( α t ) cos α 1 ) ( 1 cos α ) 2 , B 2 4 ( t ) : = 2 ( 1 cos ( α t ) ) ( 1 cos t ) + ( cos t + cos ( α t ) cos α 1 ) 2 ( 1 cos α ) 2 , B 3 4 ( t ) : = 2 ( 1 cos t ) ( cos t + cos ( α t ) cos α 1 ) ( 1 cos α ) 2 , B 4 4 ( t ) : = ( 1 cos t ) 2 ( 1 cos α ) 2 ,
for t I α and 0 < α < π . Note that the normalized B-basis can be also defined as follows:
B 0 4 ( t ) = 1 sin 4 α / 2 sin 4 ( α t ) / 2 , B 4 4 ( t ) = B 0 4 ( α t ) , B 1 4 ( t ) = 4 cos α / 2 sin 4 α / 2 sin 3 ( α t ) / 2 sin t / 2 , B 3 4 ( t ) = B 1 4 ( α t ) , B 2 4 ( t ) = 2 ( 1 + 2 cos 2 α / 2 ) sin 4 α / 2 sin 2 ( α t ) / 2 sin 2 t / 2 .
The basis functions in (1), equivalently in (2), can be regarded as the equivalent to the Bernstein polynomial basis in U 4 ( I α ) due to the following properties:
  • Symmetry:
    B i 4 ( t ) = B 4 i 4 ( α t ) , i = 0 , , 4 , t I α .
  • Positivity:
    B i 4 ( t ) 0 , i = 0 , , 4 , t I α .
  • Partition of unity:
    i = 0 4 B i 4 ( t ) = 1 , t I α .
Figure 1 shows the graphs of the normalized B-basis of U 4 ( I α ) for α = π / 3 .
Let us note that the parameter α not only controls the length of the domain but also can be considered as a shape parameter with a tension-like effect in the parametric curves generated by the normalized B-basis of U 4 . The following result will explain the behavior as α 0 . In order to prevent the normalized B-bases from losing their domain intervals I α , we shall use a reparametrization so that the new functions are defined on [ 0 , 1 ] allowing the parameter α 0 .
Lemma 1.
The normalized B-basis of U 4 ( I α ) converges uniformly to the Bernstein basis ( b 0 4 , , b 4 4 ) of degree 4 on the interval [ 0 , 1 ] whenever α 0 .
Proof. 
Using the change t = α τ , and developing by the Taylor expansion at τ = 0 , we have
b 4 4 ( τ ) B 4 4 ( α τ ) = 4 ( 1 cos α ) 2 α 4 4 ( 1 cos α ) 2 τ 4 + 2 α 6 ( 16 cos ( 2 α ξ ) cos ( α ξ ) ) 6 ! ( cos α 1 ) 2 τ 6 , τ [ 0 , 1 ] ,
for ξ [ 0 , τ ] . Then
b 4 4 ( τ ) B 4 4 ( α τ ) 4 ( 1 cos α ) 2 α 4 4 ( 1 cos α ) 2 + α 6 ( cos α 1 ) 2 , τ [ 0 , 1 ] ,
and, since
lim α 0 4 ( 1 cos α ) 2 α 4 4 ( 1 cos α ) 2 = 0 , lim α 0 α 6 ( cos α 1 ) 2 = 0 ,
we derive
lim α 0 max 0 τ 1 | b 4 4 ( τ ) B 4 4 ( α τ ) | = 0 .
By definition, b 0 4 ( τ ) = b 4 4 ( 1 τ ) , B 0 4 ( α τ ) = B ˜ 4 4 ( α ( 1 τ ) ) and therefore, using (4), we deduce that
lim α 0 max 0 τ 1 | b 0 4 ( τ ) B 0 4 ( α τ ) | = 0 .
Using the Taylor expansion and a similar reasoning, we can also write
lim α 0 max 0 τ 1 | b 3 4 ( τ ) B 3 4 ( α τ ) | = 0 , lim α 0 max 0 τ 1 | b 1 4 ( τ ) B 1 4 ( α τ ) | = 0
Finally, from the normalization property of the bases, 1 = i = 0 4 b i 4 ( τ ) = i = 0 4 B i 4 ( α τ ) , τ [ 0 , 1 ] , and taking into account Formulaes (4)–(6), we conclude
lim α 0 max 0 τ 1 | b 2 4 ( t ) B ˜ 2 4 ( t ) | = 0 .
Definition 1.
Let 0 < α < π and d N . Given p i R d , i = 0 , , 4 , we say that the parametric trigonometric curve
p 4 ( t ) : = i = 0 4 p i B i 4 ( t ) , t I α ,
is a T 4 -curve.
Let us observe that, using the symmetry property (3), we can write
p 4 ( α t ) = i = 0 4 p 4 i B i 4 ( t ) , t I α ,
and derive that T 4 -curves possess a symmetry similar to that of Bézier curves. Furthermore, since T 4 -curves are curves expressed in terms of a normalized B-basis, they admit a de Casteljau-type algorithm, which is a corner cutting algorithm providing evaluation and subdivision.
Given a T 4 -curve (8), and a parameter value t [ 0 , α ] , the B-algorithm provides certain values λ i k ( t ) that define the intermediate points p i k ( t ) of a de Casteljau-like algorithm:
p i k + 1 ( t ) = ( 1 λ i k ( t ) ) p i k ( t ) + λ i k ( t ) p i + 1 k ( t ) ,
where k = 0 , 1 , 2 , 3 and i = 0 , , 3 k . In fact, beginning with p i 0 = p i , this algorithm yields the final value
p 0 4 = p ( t ) .
In addition, the two segments in which the parameter t divides the curve have control points given by { p 0 0 , p 0 1 , p 0 2 , p 0 3 , p 0 4 } and { p 0 4 , p 1 3 , p 2 3 , p 3 1 , p 4 0 } , respectively. Here we provide compact expressions of λ i k ( t ) for an arbitrary t [ 0 , α ] :
λ 0 0 ( t ) = sin ( t / 2 ) sin ( α / 2 ) cos ( α / 2 ) cos ( t / 2 ) , λ 3 0 ( t ) = 1 λ 0 0 ( α t ) , λ 1 0 ( t ) = sin ( t / 2 ) sin ( α / 2 ) 1 + 2 cos 2 ( α / 2 ) cos ( ( α t ) / 2 ) + 2 cos ( α / 2 ) cos ( t / 2 ) , λ 2 0 ( t ) = 1 λ 1 0 ( α t ) , λ 0 1 ( t ) = sin ( t / 2 ) sin ( α / 2 ) cos ( ( α t ) / 2 ) + 2 cos ( α / 2 ) cos ( t / 2 ) 1 + 2 cos 2 ( t / 2 ) , λ 2 1 ( t ) = 1 λ 0 1 ( α t ) , λ 1 1 ( t ) = sin ( t / 2 ) sin ( α / 2 ) 3 cos ( t / 2 ) + 2 sin ( α / 2 ) sin ( ( α t ) / 2 ) cos ( α / 2 ) + 2 cos ( t / 2 ) cos ( ( α t ) / 2 ) , λ 0 2 ( t ) = sin ( t / 2 ) sin ( α / 2 ) cos ( α / 2 ) + 2 cos ( t / 2 ) cos ( ( α t ) / 2 ) 3 cos ( t / 2 ) , λ 1 2 ( t ) = 1 λ 0 2 ( α t ) , λ 0 3 ( t ) = sin ( t / 2 ) sin ( α / 2 ) cos ( ( α t ) / 2 ) .
Taking into account Lemma 1, whenever α 0 , the T 4 -curve (8) degenerates to an integral Bézier curve with control points p i , i = 0 , , 4 , and the coefficients λ i k ( t ) of the algorithm degenerate to t / α (see Figure 2). In other words, the B-algorithm reduces to the standard de Casteljau algorithm. This property can be also checked by introducing the Taylor expansion of the functions λ i k and taking limits.

3. Two Frequency Trigonometric Spline Curves

Now, we shall introduce NTP bases of U 4 ( I α ) such that their derivatives up to the third order at the ends of the interval of [ 0 , α ] allow us the definition of regular piecewise curves. We shall use the following notation:
C p , q : = p + q cos α , C p , q , r : = p + q cos α + r cos 2 α ,
for given p , q , r N .
Let us define the following systems of functions on I α :
( N 0 4 , , N 4 4 ) : = ( B 0 4 , , B 4 4 ) A 4 , ( N 0 j , 4 , , N 4 j , 4 ) : = ( B 0 4 , , B 4 4 ) A j , 4 , j = 0 , 1 , 2 ,
where ( B 0 4 , , B 4 4 ) is the normalized B-basis of U 4 ( I α ) (see (1) or (2)),
A 4 : = 12 C 1 , 1 0 0 0 0 0 6 C 1 , 1 0 0 0 0 0 2 C 2 , 1 0 0 0 0 0 6 C 1 , 1 0 0 0 0 0 12 C 1 , 1 1 1 C 5 , 6 C 5 , 6 1 0 0 2 C 1 , 1 C 3 , 4 1 0 0 1 2 C 1 , 1 1 0 0 1 C 3 , 4 2 C 1 , 1 0 0 1 C 5 , 6 C 5 , 6 1 ,
and the matrices A j , 4 , for j = 0 , 1 , 2 , are defined as follows:
A 0 , 4 : = D 0 1 0 0 0 0 0 1 0 0 0 0 3 C 1 , 2 0 0 0 3 C 3 , 4 2 C 0 , 1 0 0 3 C 1 , 2 C 7 , 18 , 12 C 1 , 10 , 12 C 1 , 2 , D 0 : = diag ( 1 , 1 , 2 C 2 , 1 , 6 C 1 , 1 , 12 C 1 , 1 C 1 , 2 ) 1 ,
A 1 , 4 : = D 1 3 C 1 , 2 C 7 , 18 , 12 C 1 , 10 , 12 C 1 , 2 0 0 4 C 1 , 1 2 C 1 , 8 , 8 C 1 , 2 0 0 2 C 1 , 1 C 1 , 6 , 4 C 1 , 2 0 0 1 C 1 , 4 C 1 , 2 0 0 2 C 1 , 1 4 C 1 , 3 C 1 , 1 C 5 , 6 C 1 , 2 C 1 , 2 , D 1 : = diag ( 12 C 1 , 1 C 1 , 2 , 6 C 1 , 1 C 1 , 2 , 2 C 2 , 1 C 1 , 2 , 3 C 1 , 2 , 12 C 1 , 1 C 1 , 2 ) 1 ,
A 2 , 4 : = D 2 2 C 1 , 1 4 C 1 , 3 C 1 , 1 C 5 , 6 C 1 , 2 C 1 , 2 0 0 2 C 1 , 1 C 3 , 4 1 0 0 1 2 C 1 , 1 1 0 0 1 C 3 , 4 2 C 1 , 1 0 0 1 C 5 , 6 C 5 , 6 1 , D 2 : = diag ( 12 C 1 , 1 C 1 , 2 , 6 C 1 , 1 , 2 C 2 , 1 , 6 C 1 , 1 , 12 C 1 , 1 ) 1 .
The following result proves that the systems (10) are all NTP bases of U 4 ( I α ) .
Theorem 1.
For 0 < α < π / 2 , the systems ( N 0 4 , , N 4 4 ) and ( N 0 j , 4 , , N 4 j , 4 ) , j = 0 , 1 , 2 , defined in Formula (10) are NTP basis of U 4 ( I α ) .
Proof. 
It can be easily checked that the matrix A 4 in (11), as well as the matrices A j , 4 , j = 0 , 1 , 2 , in (12), (13) and (14), respectively, are stochastic for 0 < α < π / 2 . Moreover, since for 0 < α < π / 2
det ( A 4 ) = 1 5184 C 0 , 1 C 1 , 2 2 C 1 , 1 4 C 2 , 1 0 ,
and
det ( A 0 , 4 ) = 1 72 C 0 , 1 C 1 , 2 C 1 , 1 2 C 2 , 1 0 , det ( A 1 , 4 ) = 1 864 C 0 , 1 C 1 , 2 C 1 , 1 3 C 2 , 1 0 , det ( A 2 , 4 ) = 1 2592 C 0 , 1 C 1 , 2 C 1 , 1 3 C 2 , 1 0 ,
we conclude that they are also nonsingular matrices and the introduced systems are all bases of U 4 ( I α ) .
Since ( B 0 4 , , B 4 4 ) is the normalized B-basis of U 4 ( I α ) (see [4]), by Corollary 3.9 (iv) of [3], it remains to prove that the matrices are TP. Let us recall that Neville elimination is an alternative procedure to Gaussian elimination and has been used to characterize TP and STP matrices. Applying the Neville elimination to the matrix A 4 in (11), we can factorize A 4 as follows:
A 4 = F 3 , 4 F 2 , 4 F 1 , 4 G 1 , 4 G 2 , 4 G 3 , 4 ,
with
F 3 , 4 : = ( 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 1 / 2 1 / 2 ) ,   F 2 , 4 : = ϕ 2 , 4 ( 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 C 2 , 1 C 1 , 2 0 0 0 0 2 / 3 1 / 3 ) ,   F 1 , 4 : = ϕ 1 , 4 ( 1 0 0 0 0 0 1 0 0 0 0 3 C 1 , 2 0 0 0 0 C 1 , 1 C 0 , 1 0 0 0 0 C 1 , 2 1 ) ,
and
ϕ 2 , 4 : = diag ( 1 , 1 , 1 , 3 C 1 , 1 , 1 ) 1 , ϕ 1 , 4 : = diag ( 1 , 1 , 2 C 2 , 1 , C 1 , 2 , 2 C 1 , 1 ) 1 .
Additionally,
G 1 , 4 : = ψ 1 , 4 1 C 11 , 12 0 0 0 0 C 11 , 12 2 C 2 , 3 0 0 0 0 2 C 2 , 3 C 5 , 6 0 0 0 0 1 0 0 0 0 0 1 , G 2 , 4 : = ψ 2 , 4 1 0 0 0 0 0 C 5 , 6 6 C 1 , 1 0 0 0 0 3 C 1 , 2 1 0 0 0 0 1 0 0 0 0 0 1 , G 3 , 4 : = ψ 3 , 4 1 0 0 0 0 0 1 0 0 0 0 0 C 5 , 6 1 0 0 0 0 1 0 0 0 0 0 1 ,
with
ψ 1 , 4 : = diag ( 12 C 1 , 1 , C 15 , 18 , C 5 , 6 , 1 , 1 ) 1 , ψ 2 , 4 : = diag ( 1 , C 11 , 12 , 2 C 2 , 3 , 1 , 1 ) 1 , ψ 3 , 4 : = diag ( 1 , 1 , 6 C 1 , 1 , 1 , 1 ) 1 .
It can be easily checked that, for 0 < α < π / 2 , the entries of the matrix factors in (15) are nonnegative and consequently, the bidiagonal matrices F k , 4 and G k , 4 , k = 1 , 2 , 3 are TP. Taking into account that, by Theorem 3.1 of [29], the product of TP matrices is a TP matrix, we conclude that A 4 is TP.
Now, we provide the bidiagonal factorization obtained by applying Neville elimination to the matrices A j , 4 , j = 0 , 1 , 2 . It can be checked that
A 0 , 4 = F 3 , 4 F 2 , 4 F 1 , 4 ,
where F 3 , 4 , F 2 , 4 and F 1 , 4 are the nonsingular, stochastic bidiagonal lower triangular matrices in the factorization (15). Furthermore,
A 1 , 4 = F 3 , 4 F 2 , 4 F 1 , 4 G 1 , 4 ( 1 ) G 2 , 4 ( 1 ) G 3 , 4 ( 1 ) ,
where F 3 , 4 , F 2 , 4 and F 1 , 4 are the matrices in the factorization (15) and
G 1 , 4 ( 1 ) : = ψ 1 , 4 ( 1 ) 1 C 3 , 4 0 0 0 0 2 C 1 , 1 C 3 , 4 ( C 1 , 2 ) 2 0 0 0 0 3 ( C 1 , 2 ) 2 2 C 0 , 1 0 0 0 0 1 0 0 0 0 0 1 , G 2 , 4 ( 1 ) : = ψ 2 , 4 ( 1 ) 1 0 0 0 0 0 C 7 , 18 , 12 2 C 1 , 6 , 6 0 0 0 0 C 1 , 6 , 6 C 7 , 18 , 12 0 0 0 0 1 0 0 0 0 0 1 , G 3 , 4 ( 1 ) : = ψ 3 , 4 ( 1 ) 1 0 0 0 0 0 1 0 0 0 0 0 C 1 , 10 , 12 C 1 , 2 0 0 0 0 1 0 0 0 0 0 1 .
with
ψ 1 , 4 ( 1 ) : = diag ( 4 C 1 , 1 , C 7 , 18 , 12 , C 3 , 14 , 12 , 1 , 1 ) 1 , ψ 3 , 4 ( 1 ) : = diag ( 1 , 1 , 2 C 1 , 6 , 6 , 1 , 1 ) 1 , ψ 2 , 4 ( 1 ) : = diag ( 1 , 3 C 1 , 2 C 3 , 4 , 3 C 1 , 1 C 1 , 2 C 1 , 10 , 12 , 1 , 1 ) 1 .
Finally,
A 2 , 4 = F 3 , 4 F 2 , 4 F 1 , 4 G 1 , 4 ( 2 ) G 2 , 4 ( 2 ) G 3 , 4 ( 2 ) ,
where F 3 , 4 , F 2 , 4 and F 1 , 4 are the matrices in the factorization (15) and
G 1 , 4 ( 2 ) : = ψ 1 , 4 ( 2 ) 1 C 5 , 12 0 0 0 0 C 5 , 12 C 1 , 6 0 0 0 0 C 1 , 2 C 1 , 6 2 C 0 , 1 0 0 0 0 1 0 0 0 0 0 1 , G 2 , 4 ( 2 ) : = ψ 2 , 4 ( 2 ) 1 0 0 0 0 0 2 C 1 , 3 3 C 1 , 2 0 0 0 0 3 C 1 , 10 , 12 2 C 1 , 3 0 0 0 0 1 0 0 0 0 0 1 , G 3 , 4 ( 2 ) : = ψ 3 , 4 ( 2 ) 1 0 0 0 0 0 1 0 0 0 0 0 C 5 , 6 1 0 0 0 0 1 0 0 0 0 0 1 ,
with
ψ 1 , 4 ( 2 ) : = diag ( 6 C 1 , 2 , 6 C 1 , 3 , C 1 , 10 , 12 , 1 , 1 ) 1 , ψ 2 , 4 ( 2 ) : = diag ( 1 , C 5 , 12 , C 5 , 36 , 36 , 1 , 1 ) 1 , ψ 3 , 4 ( 2 ) : = diag ( 1 , 1 , 6 C 1 , 1 , 1 , 1 ) 1 .
It can be easily checked that the bidiagonal matrices of the factorization (16), (17) and (18) are TP for 0 < α < π / 2 and the result follows. □
We shall see that the NTP bases in (10) allow us to define T 4 -B-spline bases and T 4 -B-spline curves, that is, piecewise functions and curves, respectively, on U 4 ( I α ) . In the following result, the derivatives up to the third order of the functions of these NTP bases are provided.
Lemma 2.
The functions N i 4 , N i j , 4 , j = 0 , 1 , 2 , i = 0 , , 4 , of the NTP bases of U 4 ( I α ) defined in (10) satisfy the following properties (see Figure 3):
1.
i = 0 4 N i 4 ( t ) = 1 , i = 0 4 N i j , 4 ( t ) = 1 , t I α .
2.
( N 0 4 ) ( k ) ( α ) = 0 , ( N 0 j , 4 ) ( k ) ( α ) = 0 , ( N 4 4 ) ( k ) ( 0 ) = 0 , ( N 4 j , 4 ) ( k ) ( 0 ) = 0 for j = 0 , 1 , 2 .
3.
( N i 4 ) ( k ) ( 0 ) = ( N i + 1 4 ) ( k ) ( α ) , k = 0 , , 3 .
4.
( N i 4 ) ( k ) ( 0 ) = ( N i + 1 2 , 4 ) ( k ) ( α ) , k = 0 , , 3 .
5.
( N i 2 , 4 ) ( k ) ( 0 ) = ( N i + 1 1 , 4 ) ( k ) ( α ) , k = 0 , , 3 .
6.
( N i 1 , 4 ) ( k ) ( 0 ) = ( N i + 1 0 , 4 ) ( k ) ( α ) , k = 0 , , 3 .
Given p N , p 4 , we are going to consider the equally spaced partition
π : = { u i } i = 0 p + 5 = { i α } i = 0 p + 5 ,
as well as the partition
μ : = { u i } i = 0 p + 5 ,
with
0 = u 0 = = u 4 < u 5 < < u p < u p + 1 = = u p + 5 ,
and
u k = ( k 4 ) α , k = 4 , , p + 1 .
Now, we define piecewise functions on either partition π or μ . When considering the equally spaced partitions π , trigonometric B-spline bases will be obtained for the shape-preserving representation of closed curves. On the other hand, using the piecewise functions on partitions μ , clamped trigonometric curves satisfying the tangency to the control-polygon will be defined.
For any i = 0 , , p , let
N i , 4 ( u ) : = N j 4 u u i + 4 j , u [ u i + 4 j , u i + 5 j ) , j = 0 , , 4 , 0 , else .
Moreover, for the partition μ , we consider the following piecewise functions,
N ˜ i , 4 ( u ) : = N j i j , 4 u u i + 4 j , u [ u i + 4 j , u i + 5 j ) , 0 j i , 0 , else ,
for i = 0 , , 3 , with the convention N 0 3 , 4 : = N 0 4 , and
N ˜ i , 4 ( u ) : = N j p ( i + 4 j ) , 4 u u i + 4 j , u [ u i + 4 j , u i + 5 j ) , i + 4 p j 4 , 0 , else ,
for i = p 3 , , p , with the convention N 4 3 , 4 : = N 4 4 . In Figure 4 the piecewise functions N i , 4 ( u ) , N ˜ i , 4 ( u ) , i = 0 , , 3 are depicted.
The following properties of the above introduced piecewise functions can be deduced from their definition and taking into account Theorem 1 and Lemma 2.
Proposition 1.
The functions N i , 4 , i = 0 , , p , defined in (21), and the functions N ˜ i , 4 , i = 0 , , 3 and i = p 3 , , p , defined in (22) and (23) satisfy the following properties:
1.
All the mentioned functions are piecewise trigonometric functions of the space U 4 ( I α ) .
2.
The functions N i , 4 ( u ) are symmetrical with respect to the middle of their supports and they can be obtained by translation, i.e.,
N i , 4 ( u ) = N 0 , 4 ( u u i ) .
3.
N i , 4 ( u ) > 0 and N ˜ i , 4 ( u ) > 0 for u ( u i , u i + 5 ) and all applicable indices i. In fact, N i , 4 and N ˜ i , 4 have minimal support [ u i , u i + 5 ] .
4.
N i , 4 ( u ) | [ u l , u l + 1 ] 0 for i = l 4 , , l .
5.
On the partitions π and μ ,
i = 0 p N i , 4 ( u ) = 1 , u [ u 4 , u p + 1 ] .
6.
On the partitions π and μ , the functions N i , 4 ( u ) for i = 0 , , p and N ˜ i , 4 ( u ) for i = 0 , , 3 and i = p 3 , , p are C j m u l t ( u k ) -continuous, where m u l t ( u k ) is the multiplicity of the knot u k in the support of the respective function.
Due to the analogy to the well-known polynomial B-splines, we will say that the introduced piecewise functions are T 4 -B-splines.
Definition 2.
Given p N , p 4 , and partitions π and μ from (19) and (20), respectively. We say that ( N 0 , 4 , π , , N p , 4 , π ) with
N i , 4 , π ( u ) : = N i , 4 ( u ) , i = 0 , , p ,
and N i , 4 defined in (21), is the normalized T 4 -B-spline basis over the partition π (see Figure 5a for an illustration). On the other hand, we say that ( N 0 , 4 , μ , , N p , 4 , μ ) with
N i , 4 , μ ( u ) : = N i , 4 ( u ) , i = 0 , , p ,
and N i , 4 = N ˜ i , 4 defined in (22), for i = 0 , , 3 , and defined in (23), for i = p 3 , , p , is the normalized T 4 -B-splines over the partition μ (see Figure 5b for an illustration).
The following properties guarantee shape preserving properties of normalized T 4 -B-spline systems for knot partitions μ and π , with 0 < α < π / 2 .
Proposition 2.
Given 0 < α < π / 2 , the T 4 -B-spline functions N i , 4 , π and N i , 4 , μ , i = 0 , , p , introduced in Definition 2, satisfy the following properties:
1.
( N 0 , 4 , π , , N p , 4 , π ) is an NTP basis of the generated space U of piecewise functions on U 4 ( I α ) defined on [ u 4 , u p + 1 ] .
2.
( N 0 , 4 , μ , , N p , 4 , μ ) is the normalized B-basis of C ( 3 functions on [ u 0 , u p + 5 ] .
Proof. 
Let ( N 0 4 , , N 4 4 ) and ( N 0 j , 4 , , N 4 j , 4 ) , j = 0 , 1 , 2 , be the NTP bases of U 4 ( I α ) defined in (10). Given u [ u i , u i + 1 ) ,
N i 4 + j , π ( u ) = N j 4 u u i , j = 0 , , 4 ,
for i = 4 , , p , and j = 0 , , 4 , we have
N i 4 + j , μ ( u ) = N j i 4 , 4 u u i , i = 4 , 5 , 6 , N j 4 u u i , i = 7 , , p 3 , N j p i , 4 u u i , i = p 2 , p 1 , p .
Consequently, taking into account Theorem 1, we deduce that the restrictions to [ u i , u i + 1 ) of ( N 0 , π , , N p , π ) and the restrictions to [ u i , u i + 1 ) of ( N 0 , μ , , N p , μ ) are NTP bases of U 4 ( I α ) for 0 < α < π / 2 . Then, it can be deduced that ( N 0 , π , , N p , π ) and ( N 0 , μ , , N p , μ ) are NTP bases of the corresponding generated space S of piecewise functions defined on [ u 4 , u p + 1 ] .
In addition, the basis ( N 0 , μ , , N p , μ ) also satisfies
lim t u j 4 + ( N k , μ ( t ) / N j , μ ( t ) ) = 0 , lim t u k + 1 , μ ( N j , μ ( t ) / N k , μ ( t ) ) = 0 ,
whenever 0 j < k p . Then, by Theorem 3.2 of Chapter 4 of [1], ( N 0 , μ , , N p , μ ) is the normalized B-basis of U . □
The previous result implies that the T 4 -B-spline basis has optimal shape preserving properties (see [3] and Chapter 4 of [1]) and it also has optimal stability properties for the evaluation (cf. Chapter 5 of [1]) Now, we can define the corresponding piecewise trigonometric curves.
Definition 3.
Given d , p N , p 4 , let s i R d , i = 0 , , p , and partitions π and μ from (19) and (20), respectively. The parametric curve defined by
s ( u ) : = i = 0 p s i N i , 4 , π ( u ) , u [ u 4 , u p + 1 ) ,
is called T 4 -B-spline curve with respect to the partition π and control points s 0 , , s p . In particular, if p = m , with m N , m 4 , then we say that s ( u ) is an open T 4 -B-Spline curve. If p = m + 4 and s m + 1 + i = s i , i = 0 , , 3 , we say that s ( u ) is a closed T 4 -B-Spline curve. See Figure 6a for an illustration of a closed T 4 -B-Spline curve.
The parametric curve defined by
s ( u ) : = i = 0 p s i N i , 4 , μ ( u ) , u [ u 4 , u p + 1 ) ,
is called clamped T 4 -B-spline curve with respect to the partition μ and control points s 0 , , s p . See Figure 6b for an illustration of a clamped T 4 -B-Spline curve.
T 4 -B-Spline curves satisfy the following properties.
Proposition 3.
The T 4 -B-Spline curves described in Definition 2 satisfy:
1.
The relation between a T 4 -B-Spline curve and its control points is affinely invariant.
2.
Any T 4 -B-Spline curve s ( u ) is locally controlled, i.e., moving a control point s l only modifies the curve for u [ u l , u l + 5 ) , moreover for τ = π or τ = μ we have
s ( u ) | u [ u l , u l + 1 ) = i = l 4 l s i N i , 4 , τ ( u ) ,
and the curve s ( u ) lies in the convex hull of its control points s i , i = l 4 , , l .
3.
The T 4 -B-Spline curves are monotonicity preserving: the curve has the same monotonicity as the monotone control points.
4.
The length of a T 4 -B-Spline curve is bounded above by the length of its control polygon.
5.
If the control polygon of a T 4 -B-Spline curve is planar and convex, then the T 4 -B-Spline curve is also planar and convex.
6.
The T 4 -B-Spline curve never crosses a hyperplane more often than does the control polygon.
7.
Clamped T 4 -B-Spline curves have end point and end tangent interpolation properties:
s ( 0 ) = s 0 , s ( u m + 1 ) = s m , s ( 0 ) = cot α / 2 ( s 1 s 0 ) , s ( u m + 1 ) = cot α / 2 ( s m s m 1 ) .
Finally, let us note that Definition 1, relations (10), the corner cutting algorithm for T 4 -curves in (9) and the matrix factorizations (11)–(14) result in a corner cutting algorithm for T 4 -B-Spline curves analogous to the one for T 2 -B-Spline curves detailed in [18].

4. Convergence of T 4 -B-Spline Curves to Quartic Polynomial B-Spline Curves

In this section we prove the convergence, when α 0 , of T 4 -B-Spline curves to polynomial B-spline curves.
Theorem 2.
Let π and μ partitions described in (19) and (20), respectively. When α 0 , the T 4 -B-Spline curve (24) and the clamped T 4 -B-Spline curve (25), with respect to π and μ , respectively, and control points s 0 , , s p approaches uniformly to the quartic polynomial B-spline curve with knot vector π and μ , respectively and control points s 0 , , s p .
Proof. 
Let us observe that, by (26), for τ = π or τ = μ , we can write
s ( u ) | u [ u l , u l + 1 ) = i = l 4 l s i N i , 4 , τ ( u ) = j = 0 4 s l + j 4 N j 4 ( u u l ) .
Let τ : = ( u u l ) / α for reparameterizing each segment curve on the interval 0 τ 1 . By Lemma 1, as α 0 , the function B i 4 ( α τ ) approaches uniformly the Bernstein polynomial b i 4 ( τ ) , 0 τ 1 , for all i = 0 , , 4 . It can be easily checked that that the matrix A 4 in (11) satisfies
lim α 0 A 4 = A , A : = 1 / 24 11 / 24 11 / 24 1 / 24 0 0 1 / 3 7 / 12 1 / 12 0 0 1 / 6 2 / 3 1 / 6 0 0 1 / 12 7 / 12 1 / 3 0 0 1 / 24 11 / 24 11 / 24 1 / 24 ,
and therefore for τ = π we obtain
lim α 0 s ( u ) | u [ u l , u l + 1 ) = lim α 0 B 0 4 ( α τ ) , , B 4 4 ( α τ ) A 4 s l 4 s l = b 0 4 ( τ ) , , b 4 4 ( τ ) A s l 4 s l ,
which is the matrix form of a uniform B-spline curve of degree 4.
For τ = μ , we can follow a similar reasoning taking into account that the matrices of () satisfy
lim α 0 A 0 , α = A 0 , A 0 = 1 0 0 0 0 0 1 0 0 0 0 1 / 2 1 / 2 0 0 0 1 / 4 7 / 12 1 / 6 0 0 1 / 8 37 / 72 23 / 72 1 / 24 ,
lim α 0 A 1 , α = A 1 , A 1 = 1 / 8 37 / 72 23 / 72 1 / 24 0 0 4 / 9 17 / 36 1 / 12 0 0 2 / 9 11 / 18 1 / 6 0 0 1 / 9 5 / 9 1 / 3 0 0 1 / 18 4 / 9 11 / 24 1 / 24 ,
and
lim α 0 A 2 , α = A 2 , A 2 = 1 / 18 4 / 9 11 / 24 1 / 24 0 0 1 / 3 7 / 12 1 / 12 0 0 1 / 6 2 / 3 1 / 6 0 0 1 / 12 7 / 12 1 / 3 0 0 1 / 24 11 / 24 11 / 24 1 / 24 .

5. Conclusions and Future Work

We have proposed two frequency trigonometric spline bases with shape preserving properties associated to uniform knot vectors. The corresponding parametric trigonometric spline curves have been also described. It is also shown that these curves share many properties of polynomial spline curves. In fact, they converge to uniform quartic B-spline curves.
There is some worthwhile work to study further. We want to extend the bases to knot vectors with multiple knots and to investigate whether it is possible to construct new Algebraic-Trigonometric Pythagorean-Hodograph B-Splines curves taking into account the results from [21].

Author Contributions

Conceptualization, G.A., E.M., J.M.P. and B.R.; methodology, G.A., E.M., J.M.P. and B.R.; investigation, G.A., E.M., J.M.P. and B.R.; writing—original draft, G.A., E.M., J.M.P. and B.R.; writing—review and editing, G.A., E.M., J.M.P. and B.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by the CICYT under Research Grant CICYT BFM2000–1253. This research was also partially supported by Spanish research grants PGC2018-096321-B-I00 (MCIU/AEI) and RED2022-134176-T (MCI/AEI) and by Gobierno de Aragón (E41_23R).

Data Availability Statement

The data and codes used in this work are available under request. The code of the experimentation was developed in Matlab R2022b. All experiments were ran on a Apple M1 Pro chip with 10-Core CPU and 32 GB RAM.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peña, J.M. (Ed.) Shape Preserving Representations in Computer-Aided Geometric Design; Nova Science Publishers Commack: New York, NY, USA, 1999. [Google Scholar]
  2. Goodman, T.N.T. Shape preserving representations, In Mathematical Methods in CAGD; Lyche, T., Schumaker, L.L., Eds.; Academic Press: Boston, MA, USA, 1989; pp. 333–351. [Google Scholar]
  3. Carnicer, J.M.; Peña, J.M. Totally positive bases for shape preserving curve design and optimality of B-splines. Comput. Aided Geom. Des. 1994, 11, 635–656. [Google Scholar] [CrossRef]
  4. Sánchez-Reyes, J. Harmonic rational Bézier curves, p-Bézier curves and trigonometric polynomials. Comput. Aided Geom. Des. 1998, 9, 909–924. [Google Scholar] [CrossRef]
  5. Han, X.L.; Zhu, Y.P. Curve construction based on five trigonometric blending functions. Bit Numer. Math. 2012, 52, 953–979. [Google Scholar] [CrossRef]
  6. Wang, K.; Zhang, G.C.; Nhon, N.T. New trigonometric basis possessing denominator shape parameters. Math. Probl. Eng. 2018. [Google Scholar] [CrossRef]
  7. Han, X.L.; Ma, Y.C.; Huang, X.L. The cubic trigonometric Bézier curve with two shape parameters. Appl. Math. Lett. 2008, 22, 226–231. [Google Scholar] [CrossRef]
  8. Han, X.L. Cubic trigonometric polynomial curves with a shape parameter. Comput. Aided Geom. Des. 2004, 21, 535–548. [Google Scholar] [CrossRef]
  9. Wu, X.Q.; Han, X.L.; Luo, S. Quadratic trigonometric polynomial Bézier curves with a shape parameter. Eng. Graph. 2008, 1, 82–87. [Google Scholar]
  10. Majeed, A.; Abbas, M.; Qayyum, F.; Miura, K.T.; Misro, M.Y.; Nazir, T. Geometric modeling using new cubic trigonometric B-spline functions with shape parameter. Mathematics 2020, 8, 2102. [Google Scholar] [CrossRef]
  11. Zhu, Y.P.; Han, X.L. New trigonometric basis possessing exponential shape parameters. J. Comput. Math. 2015, 33, 642–684. [Google Scholar]
  12. Zhu, Y.P.; Liu, Z. A class of trigonometric Bernstein-type basis functions with four shape parameters. Math. Probl. Eng. 2019. [Google Scholar] [CrossRef]
  13. Schoenberg, I.J. On trigonometric spline interpolation. J. Math. Mech. 1964, 13, 795–825. [Google Scholar]
  14. Schumaker, L.L. Spline Functions: Basic Theory; John Wiley and Sons: New York, NY, USA, 1981. [Google Scholar]
  15. Lyche, T.; Winther, R. A stable recurrence relation for trigonometric B-splines. J. Approx. Theory 1979, 25, 266–279. [Google Scholar] [CrossRef]
  16. Yağmurlu, N.M.; Karakaş, A.S. Numerical solutions of the equal width equation by trigonometric cubic B-spline collocation method based on Rubin–Graves type linearization. Numer. Methods Partial Differ. Equ. 2020, 36, 1170–1183. [Google Scholar] [CrossRef]
  17. Yağmurlu, N.M.; Karakaş, A.S. A novel perspective for simulations of the MEW equation by trigonometric cubic B-spline collocation method based on Rubin-Graves type linearization. Comput. Methods Differ. Equ. 2022, 10, 1046–1058. [Google Scholar]
  18. Albrecht, G.; Mainar, E.; Peña, J.M.; Rubio, B. A New Class of Trigonometric B-Spline Curves. Symmetry 2023, 15, 1551. [Google Scholar] [CrossRef]
  19. Farouki, R.T. The conformal map of the hodograph plane. Comput. Aided Geom. Des. 1994, 11, 363–390. [Google Scholar] [CrossRef]
  20. Sˇir, Z.; Jüttler, B. Spatial Pythagorean Hodograph Quintics and the Approximation of Pipe Surfaces. In Mathematics of Surfaces XI. Lecture Notes in Computer Science; Martin, R., Bez, H., Sabin, M., Eds.; Springer: Berlin/Heidelberg, Germany, 2005; Volume 3604. [Google Scholar]
  21. Albrecht, G.; Beccari, C.V.; Canonne, J.C.; Romani, L. Planar Pythagorean-Hodograph B-Spline curves. Comput. Aided Geom. Des. 2017, 57, 57–77. [Google Scholar] [CrossRef]
  22. Farouki, R.T.; Giannelli, C.; Sestini, A. Identification and “reverse engineering” of Pythagorean-hodograph curves. Comput. Aided Geom. Des. 2015, 34, 21–36. [Google Scholar] [CrossRef]
  23. Farouki, R.T. Pythagorean-Hodograph Curves: Algebra and Geometry Inseparable; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  24. Farouki, R.T.; Giannelli, C.; Sestini, A. Local modification of Pythagorean-hodograph quintic spline curves using the B-spline form. Adv. Comput. Math. 2016, 42, 199–225. [Google Scholar] [CrossRef]
  25. Pelosi, F.; Sampoli, M.L.; Farouki, R.T.; Manni, C. A control polygon scheme for design of planar PH quintic spline curves. Comput. Aided Geom. Des. 2007, 24, 28–52. [Google Scholar] [CrossRef]
  26. Farin, G. Curves and Surfaces for Computer Aided Geometric Design, 4th ed.; Academic Press: San Diego, CA, USA, 1997. [Google Scholar]
  27. Hoschek, J.; Lasser, D. Fundamentals of Computer Aided Geometric Design; AKPeters: Wellesley, MA, USA, 1993. [Google Scholar]
  28. Mainar, E.; Peña, J.M.; Sánchez-Reyes, J. Shape preserving alternatives to the rational Bézier model. Comput. Aided Geom. Des. 2001, 18, 37–60. [Google Scholar] [CrossRef]
  29. Ando, T. Totally positive matrices. Linear Algebra Appl. 1987, 90, 165–219. [Google Scholar] [CrossRef]
Figure 1. Normalized B-basis of U 4 ( I α ) for α = π / 3 .
Figure 1. Normalized B-basis of U 4 ( I α ) for α = π / 3 .
Symmetry 15 02041 g001
Figure 2. T 4 -curves (8) associated to a control polygon with several parameters α , as well as the Bézier curve corresponding to the mentioned control polygon.
Figure 2. T 4 -curves (8) associated to a control polygon with several parameters α , as well as the Bézier curve corresponding to the mentioned control polygon.
Symmetry 15 02041 g002
Figure 3. NTP bases defined in (10) for α = π 3 . (a) ( N 0 0 , 4 , , N 4 0 , 4 ) , (b) ( N 0 1 , 4 , , N 4 1 , 4 ) , (c) ( N 0 2 , 4 , , N 4 2 , 4 ) , (d) ( N 0 4 , , N 4 4 ) .
Figure 3. NTP bases defined in (10) for α = π 3 . (a) ( N 0 0 , 4 , , N 4 0 , 4 ) , (b) ( N 0 1 , 4 , , N 4 1 , 4 ) , (c) ( N 0 2 , 4 , , N 4 2 , 4 ) , (d) ( N 0 4 , , N 4 4 ) .
Symmetry 15 02041 g003
Figure 4. Trigonometric piecewise functions (a) N ˜ 0 , 4 ( u ) , (b) N ˜ 1 , 4 ( u ) , (c) N ˜ 2 , 4 ( u ) , (d) N ˜ 3 , 4 ( u ) and (e) N i , 4 ( u ) , for α = π 3 .
Figure 4. Trigonometric piecewise functions (a) N ˜ 0 , 4 ( u ) , (b) N ˜ 1 , 4 ( u ) , (c) N ˜ 2 , 4 ( u ) , (d) N ˜ 3 , 4 ( u ) and (e) N i , 4 ( u ) , for α = π 3 .
Symmetry 15 02041 g004
Figure 5. Normalized T 4 -B-splines over a partition π (a), and μ (b) for α = π 3 .
Figure 5. Normalized T 4 -B-splines over a partition π (a), and μ (b) for α = π 3 .
Symmetry 15 02041 g005
Figure 6. T 4 -B-splines curves (closed in (a), clamped in (b)) associated to a control polygon with several parameters α , as well as the B-spline curve corresponding to the mentioned control polygon.
Figure 6. T 4 -B-splines curves (closed in (a), clamped in (b)) associated to a control polygon with several parameters α , as well as the B-spline curve corresponding to the mentioned control polygon.
Symmetry 15 02041 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Albrecht, G.; Mainar, E.; Peña, J.M.; Rubio, B. A Shape Preserving Class of Two-Frequency Trigonometric B-Spline Curves. Symmetry 2023, 15, 2041. https://doi.org/10.3390/sym15112041

AMA Style

Albrecht G, Mainar E, Peña JM, Rubio B. A Shape Preserving Class of Two-Frequency Trigonometric B-Spline Curves. Symmetry. 2023; 15(11):2041. https://doi.org/10.3390/sym15112041

Chicago/Turabian Style

Albrecht, Gudrun, Esmeralda Mainar, Juan Manuel Peña, and Beatriz Rubio. 2023. "A Shape Preserving Class of Two-Frequency Trigonometric B-Spline Curves" Symmetry 15, no. 11: 2041. https://doi.org/10.3390/sym15112041

APA Style

Albrecht, G., Mainar, E., Peña, J. M., & Rubio, B. (2023). A Shape Preserving Class of Two-Frequency Trigonometric B-Spline Curves. Symmetry, 15(11), 2041. https://doi.org/10.3390/sym15112041

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop