Hindawi Publishing Corporation International Journal of Mathematics and Mathematical Sciences Volume 2013, Article ID 893414, 8 pages http://dx.doi.org/10.1155/2013/893414
Research Article Four-Point π-Ary Interpolating Subdivision Schemes Ghulam Mustafa and Robina Bashir Department of Mathematics, The Islamia University of Bahawalpur, Punjab, Bahawalpur 63100, Pakistan Correspondence should be addressed to Ghulam Mustafa;
[email protected] Received 27 July 2013; Revised 27 September 2013; Accepted 27 September 2013 Academic Editor: Palle E. Jorgensen Copyright © 2013 G. Mustafa and R. Bashir. 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. We present an efficient and simple algorithm to generate 4-point n-ary interpolating schemes. Our algorithm is based on three simple steps: second divided differences, determination of position of vertices by using second divided differences, and computation of new vertices. It is observed that 4-point n-ary interpolating schemes generated by completely different frameworks (i.e., Lagrange interpolant and wavelet theory) can also be generated by the proposed algorithm. Furthermore, we have discussed continuity, H¨older regularity, degree of polynomial generation, polynomial reproduction, and approximation order of the schemes.
1. Introduction In general, subdivision schemes can be divided into two categories: approximating and interpolating. For interpolating curve subdivision, new vertices are computed and added to the old polygons for each time of subdivision and the limit curve passes through all the vertices of the original control polygon. Interpolating subdivision schemes are more attractive than approximating schemes in computer aided geometric designs because of their interpolation property. In addition, the interpolation subdivisions are more intuitive to the users. Initial work on interpolating subdivision schemes was started by Dubuc [1]. Later on, Deslauriers and Dubuc [2] have introduced a family of schemes by using Lagrange polynomials indexed by the size of the mask and the arity. In [3], Dyn et al. have studied a family of interpolating schemes with mask size of four. Consequent to this, the research communities are interested in introducing higher arity schemes (i.e., ternary, quaternary, . . . , π-ary) which give better results and less computational cost. Lian [4] has constructed both the 2π-point π-ary for any π β₯ 2 and (2π + 1)-point π-ary for any odd π β₯ 3 interpolatory subdivision schemes for curve design by using wavelet theory. Mustafa and Rehman [5] have presented general formulae for the mask of (2π+4)-point π-ary interpolating and approximating schemes for any integer π β₯ 0 and π β₯ 2. These formulae
provide mask of higher arity schemes and generalize lower arity schemes. Mustafa et al. [6] have presented an explicit formula for the mask of odd points π-ary, for any odd π β₯ 3, interpolating subdivision schemes. In [7], it has been proved that the large support and higher arity schemes may outperform than small support and lower arity schemes. Even though these schemes are not in practice. It has been suggested that the research on large support and higher arity schemes may continue. The multistage approach is very handy to construct subdivision schemes. This idea is variously used by others. Lane and Riesenfeld [8] have presented two fast subdivision algorithms for the evaluation of B-spline and Bernstein curves and surfaces. They have expressed new idea, in which, after a single duplication stage in which the number of control points is doubled by just taking each point twice, a sequence of smoothing operators is applied. Catmull and Clark [9] have used this technique to present the original description of subdivision in which each refinement is expressed in three stages. Zorin and SchrΒ¨oder [10] have considered the construction of an increasing sequence of alternating primal/dual quadrilateral subdivision schemes based on a simple averaging approach. Oswald and SchrΒ¨oder [11] used the same motif to generate families of subdivision schemes. AugsdΒ¨orfer et al. [12] first describe the original 4-point binary subdivision scheme and then apply six variations on the scheme which are obtained by tuning the local stages
2
International Journal of Mathematics and Mathematical Sciences
in various ways, producing some interesting subdivision schemes all of which are improvements on the original. We generalize the same technique to build the family of fourpoint π-ary interpolating subdivision schemes. It is observed that 4-point π-ary interpolating schemes introduced by [2, 4] can also be constructed by our generalized technique, even though these schemes have been constructed by different frameworks. This paper is organized as follows: Section 2 presents some preliminary results. Section 3 consists of multistep algorithm which generates 4-point π-ary interpolating subdivision schemes. Analysis of two schemes is also presented in this section. In Section 4, HΒ¨older regularity, polynomial generation, polynomial reproductions and approximation order of ternary and quaternary subdivision schemes have been discussed. Numerical examples and conclusion are presented in Section 5.
2. Preliminaries A general compact form of univariate π-ary subdivision scheme S [13] which maps polygon ππ = {πππ }πβZ to a refined polygon ππ+1 = {πππ+1 }πβZ is defined by πππ+1 = β πππβπ πππ ,
π β Z,
πβZ
The subdivision scheme π1 with Laurent polynomial π(1) (π§) is related to the scheme π with Laurent polynomial π(π§) by the following theorem. Theorem 1 (see [14]). Let π denote a subdivision scheme with Laurent polynomial π(π§) satisfying (5). Then there exists a subdivision scheme π1 with the property Ξππ = π1 Ξππβ1 ,
π β πππ ); π β Z}. where ππ = ππ π0 and Ξππ = {(Ξππ )π = ππ (ππ+1 Furthermore, π is a uniformly convergent if and only if (1/π)π1 converges uniformly to zero function for all initial data π0 , in the sense that π 1 lim ( π1 ) π0 = 0. πββ π
β πππ = β πππ+1 = β πππ+2 = β
β
β
= β πππ+πβ1 = 1. πβZ
πβZ
πβZ
(1)
{ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ βπββ = max { β σ΅¨σ΅¨σ΅¨σ΅¨πππ σ΅¨σ΅¨σ΅¨σ΅¨ , β σ΅¨σ΅¨σ΅¨σ΅¨πππ+1 σ΅¨σ΅¨σ΅¨σ΅¨ , β σ΅¨σ΅¨σ΅¨σ΅¨πππ+2 σ΅¨σ΅¨σ΅¨σ΅¨ πβZ πβZ {πβZ
σ΅¨ σ΅¨ lim sup σ΅¨σ΅¨σ΅¨πππ β π (πβπ π)σ΅¨σ΅¨σ΅¨σ΅¨ = 0.
π β βπβππ πΌ σ΅¨
(3)
Obviously, π = π π . A symbol called Laurent polynomial π (π§) = β ππ π§π
(4)
πβZ
of the mask π = {ππ : π β Z} plays an efficient role to analyze the convergence and smoothness of subdivision scheme. From (2) and (4) the Laurent polynomial of convergent subdivision scheme satisfies π (πΌππ ) = 0,
π = 1, 2, . . . , π β 1,
π (1) = π,
(5)
where πΌππ = exp((2ππ/π)π) are the πth root of unity. This condition guarantees the existence of a related subdivision scheme for the divided differences of the original control points and the existence of an associated Laurent polynomial π(1) (π§) = ππ§πβ1 (
σ΅©σ΅© πΏσ΅© σ΅©σ΅© { σ΅¨ } σ΅¨σ΅¨ σ΅©σ΅© 1 σ΅¨ ; π = 0, 1, 2, . . . , ππΏ β 1 , σ΅©σ΅©( ππ½ ) σ΅©σ΅©σ΅© = max { β σ΅¨σ΅¨σ΅¨σ΅¨π[π½,πΏ] πΏπσ΅¨ σ΅¨ } π+π σ΅©σ΅© π σ΅©σ΅© σ΅¨ σ΅¨ σ΅© σ΅©β {πβZ } (10) where π[π½,πΏ] (π§) =
β 0
1βπ§ ) π (π§) . 1 β π§π
(6)
(9)
σ΅¨ σ΅¨} , . . . , β σ΅¨σ΅¨σ΅¨σ΅¨πππ+πβ1 σ΅¨σ΅¨σ΅¨σ΅¨} , πβZ }
(2)
A subdivision scheme is uniformly convergent if for any initial data π0 = {ππ0 : π β Z}, there exists a continuous function π such that for any closed interval πΌ β R, it satisfies
(8)
The above theorem indicates that for any given scheme π, with the mask π satisfying (2), we can prove the uniform convergence of π by deriving the mask of (1/π)π1 and computing β((1/π)π1 )π ββ for π = 1, 2, 3 . . . , πΏ, where πΏ is the first integer for which β((1/π)π1 )πΏ ββ < 1. If such an πΏ exists, then π converges uniformly. Since there are βπβ rules for computing the values at the next refinement level, so we define the norm
where the set π = {ππ : π β Z} of coefficients is called the mask at πth level of refinement. A necessary condition for the uniform convergence of subdivision scheme (1) is that πβZ
(7)
ππ½ (π§) = (ππ§πβ1 ( = (ππ§πβ1 (
π 1 πΏβ1 βππ½ (π§π ) , πΏ π π=0
1βπ§ )) ππ½β1 (π§) 1 β π§π 1βπ§ π½ )) π (π§) , 1 β π§π
(11) π½ β©Ύ 1.
Theorem 2 (see [13]). Let π be the subdivision scheme with a characteristic β§-polynomial π(π§) = π ((π§π β 1)/(ππ§πβ1 (π§ β 1))) π(π§), π β β§. If the subdivision scheme ππ , corresponding to the β§-polynomial π(π§), converges uniformly, then πβ π0 β πΆπ (R) for any initial control polygon π0 . Corollary 3 (see [13]). If π is a subdivision scheme of the form above and (1/π)ππ+1 converges uniformly to the zero function for all initial data π0 , then πβ π0 β πΆπ (R) for any initial control polygon π0 .
International Journal of Mathematics and Mathematical Sciences p1 p2 b
z
1 n
2 n
Β·Β·Β·
pnβ1
Β·Β·Β·
nβ1 n
b
c
a
d
z
e
Figure 1: Labeling a sample control polygon. The newly inserted points between old vertices π and π are referred to as π1 , π2 ,. . ., ππβ1 , respectively.
The above Corollary 3 indicates that for any given πary subdivision scheme π, we can prove πβ π0 β πΆπ by first deriving the mask of (1/π)ππ+1 and then computing β((1/π)ππ+1 )π ββ for π = 1, 2, 3, . . . , πΏ (where πΏ is the first integer for which β((1/π)ππ+1 )πΏ ββ < 1). If such an πΏ exists, then πβ π0 β πΆπ . Definition 4 (see [15]). For any subdivision scheme ππ we denote by π = πσΈ (1)/π the corresponding parametric shift and attach the data πππ for π β Z, π β N to the parameter values π‘ππ = π‘0π +
π ππ
with π‘0π = π‘0πβ1 β
π . ππ
(12)
Theorem 5 (see [15]). A convergent subdivision scheme ππ reproduces polynomials of degree π with respect to the parameterizations (12) if and only if πβ1
π(π) (1) = πβ (π β π) , π=0
π = 1, 2, . . . , π β 1
3
π(π) (πΌππ ) = 0,
p2
1 3
2 3
c
d
a
(ii) Determine the position of vertices by using divided differences. In π-ary subdivision scheme each segment is divided into π subsegments at each refinement level. First point is inserted at the position 1/π and second point at the position 2/π and proceeding in the same way the (πβ1)-th point at the position (πβ1)/π. By using divided differences π·π and π·π , we calculate the position of (π β 1)-th newly inserted points between two old vertices π and π by π·ππ = (
πβπ π ) π·π + ( ) π·π , π π
π = 1, 2, 3, . . . , π β 1. (15)
(iii) Computation new vertices. Finally, we calculate positions of new vertices π1 , π2 ,. . ., ππβ1 by using π·π1 , π·π2 ,. . ., π·ππβ1 , respectively, by π·π1 = π2 β 2π1 + π, π·ππ = ππ+1 β 2ππ + ππβ1 ,
where πΌππ = exp((2ππ/π)π), π = 1, 2, . . . , π β 1. Theorem 6 (see [16]). A convergent subdivision scheme ππ that reproduces polynomial ππ (set of polynomials at most degree π) has an approximation order of π + 1.
e
Figure 2: Labeling a sample control polygon. The newly inserted points between old vertices π and π are referred to as π1 and π2 , respectively.
(13)
for π = 0, . . . , π,
p1
(16)
π·ππβ1 = π β 2ππβ1 + ππβ2 , where π = 2, 3, . . . , π β 2. By solving, above set of equations, we get the position of new vertices π1 , π2 , . . . , ππβ1 .
3. Multistep Algorithm We construct 4-point π-ary interpolating subdivision schemes by using three-step algorithm instead of using Lagrange polynomial and wavelets theory, and so forth. These three steps are as follows. (i) Calculate second divided differences.
3.1. Examples. A 4-Point Ternary Interpolating Scheme. In ternary subdivision scheme each segment is divided into three subsegments at each refinement level. One point is inserted at the position 1/3 and another point at the position 2/3 (see Figure 2). For π = 3 in (14), we get second divided differences π·π and π·π at points π and π
At each old vertex compute the second divided difference π·; that. is π·π is the second divided difference at point π and π·π is the second divided difference at point π (see Figure 1): π·π =
π β 2π + π , π2
π β 2π + π π·π = , π2 where π = 3, 4, . . ..
(14)
π·π =
π β 2π + π , 9
(17)
π β 2π + π π·π = . 9 For π = 3 in (15), we get π·ππ =
3βπ π π·π + π·π , 3 3
π = 1, 2.
(18)
4
International Journal of Mathematics and Mathematical Sciences p1
p2
p3
1 4
2 4
3 4
b
c
π4ππ+1 = πππ ,
a
z
4-point quaternary interpolating scheme of Deslauriers and Dubuc
d
e
Figure 3: Labeling a sample control polygon. The newly inserted points between old vertices π and π are referred to as π1 , π2 , and π3 respectively.
π+1 π4π+1 =β
π+1 =β π4π+2 π+1 =β π4π+3
By using (17), we get π·π1 = π·π2 =
(19)
π β 2π + 2π . 27
π·π1 = π2 β 2π1 + π,
π2 =
π + 2π β π·π1 β 2π·π2 3
π2 =
Remark 7. By substituting π β₯ 3 in (14)β(16), we get the mask of 4-point π-ary interpolating scheme generated by two different frameworks, that is, Lagrange interpolation [2] and wavelet theory [4]. In this paper, we propose alternative approach completely different from these approaches. In coming section, we discuss two existing schemes also produced by our framework.
, (21) .
3.2.1. Analysis of 4-Point Ternary Subdivision Scheme. The Laurent polynomial π(π§) for the scheme (23) is π (π§) =
By using (19), we get π1 =
5 π 35 π 105 π 7 π π + π + π β π . 128 πβ1 128 π 128 π+1 128 π+2
3.2. Analysis of Subdivision Schemes. Here we present the analysis of 4-point ternary and quaternary interpolating subdivision schemes. Analysis of other schemes can be done in a similar way.
This implies
3
(24)
(20)
π·π2 = π β 2π2 + π1 .
2π + π β 2π·π1 β π·π2
1 π 9 9 π 1 π β ππ+2 , ππβ1 + πππ + ππ+1 16 16 16 16
This scheme was also reconstructed by [4] in 2009.
2π β 3π + π , 27
For π = 3 in (16), we have
π1 =
7 π 105 π 35 π 5 π π + π + π β π , 128 πβ1 128 π 128 π+1 128 π+2
β5π + 60π + 30π β 4π , 81 β4π + 30π + 60π β 5π . 81
+ 60π§1 + 81 + 60π§β1 + 30π§β2
Using (11) for π = 3, π½ = 1, 2 and πΏ = 1, we get π[1,1] (π§) =
π3ππ+1 = πππ , 5 π 60 30 π 4 π = β ππβ1 + πππ + ππ+1 β ππ+2 , 81 81 81 81
π+1 π3π+2
4 π 30 60 π 5 π = β ππβ1 + πππ + ππ+1 β ππ+2 . 81 81 81 81
(25)
β 5π§β4 β 4π§β5 } .
(22)
Now 4-point ternary scheme can be written as
π+1 π3π+1
1 {β4π§5 β 5π§4 + 30π§2 81
1 4 1 π1 (π§) = β π§5 β π§4 3 81 81 +
5 3 26 2 29 1 26 5 π§ + π§ + π§ + + π§β1 81 81 81 81 81
β
1 β2 4 β3 π§ β π§ , 81 81
(23)
The above scheme was introduced by Deslauriers and Dubuc [2] in 1989 by using Lagrange interpolant. Later on, this scheme was also reconstructed by [4] by using wavelet theory. A 4-Point Quaternary Interpolating Scheme. In quaternary subdivision scheme each segment is divided into four subsegments at each refinement level. First, second and third points are inserted at the positions 1/4, 2/4, and 3/4, respectively (see Figure 3). For π = 4 in (14)β(16), we get the following
π[2,1] (π§) =
1 2 1 4 π (π§) = β π§5 + π§4 + π§3 3 2 27 9 9
(26)
(27)
17 2 1 4 + π§2 + π§1 + β π§β1 . 27 9 9 27 If ππ½ is the scheme corresponding to ππ½ (π§), then by (10) σ΅©σ΅© 1 σ΅©σ΅© σ΅©σ΅© π σ΅©σ΅© = max { β σ΅¨σ΅¨σ΅¨σ΅¨π[π½,1] σ΅¨σ΅¨σ΅¨σ΅¨ : π = 0, 1, 2} , σ΅©σ΅© π½ σ΅©σ΅© } { σ΅¨σ΅¨ π+3π σ΅¨σ΅¨ σ΅© 3 σ΅©β } {πβZ
π½ = 1, 2. (28)
International Journal of Mathematics and Mathematical Sciences Using (9), (26), and (27), we get σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅©σ΅© 1 σ΅©σ΅© σ΅©σ΅© π σ΅©σ΅© = max {σ΅¨σ΅¨σ΅¨ β4 σ΅¨σ΅¨σ΅¨ + σ΅¨σ΅¨σ΅¨ 26 σ΅¨σ΅¨σ΅¨ + σ΅¨σ΅¨σ΅¨ 5 σ΅¨σ΅¨σ΅¨ , σ΅¨σ΅¨σ΅¨ β1 σ΅¨σ΅¨σ΅¨ + σ΅¨σ΅¨σ΅¨ 29 σ΅¨σ΅¨σ΅¨ + σ΅¨σ΅¨σ΅¨ β1 σ΅¨σ΅¨σ΅¨} , σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅©σ΅© 1 σ΅©σ΅© σ΅¨ 81 σ΅¨ σ΅¨ 81 σ΅¨ σ΅¨ 81 σ΅¨ σ΅¨ 81 σ΅¨ σ΅¨ 81 σ΅¨ σ΅¨ 81 σ΅¨ σ΅© 3 σ΅©β σ΅©σ΅© 1 σ΅©σ΅© σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅¨ σ΅©σ΅© π σ΅©σ΅© = max {σ΅¨σ΅¨σ΅¨ β4 σ΅¨σ΅¨σ΅¨ + σ΅¨σ΅¨σ΅¨ 17 σ΅¨σ΅¨σ΅¨ + σ΅¨σ΅¨σ΅¨ β4 σ΅¨σ΅¨σ΅¨ , σ΅¨σ΅¨σ΅¨ 1 σ΅¨σ΅¨σ΅¨ + σ΅¨σ΅¨σ΅¨ 2 σ΅¨σ΅¨σ΅¨} . σ΅©σ΅© 2 σ΅©σ΅© σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅© 3 σ΅©β σ΅¨ 27 σ΅¨ σ΅¨ 27 σ΅¨ σ΅¨ 27 σ΅¨ σ΅¨ 9 σ΅¨ σ΅¨ 9 σ΅¨
Proof. The Laurent polynomial (25) of the scheme (23) can be written as 4
Remark 8. Similarly, we can prove that quaternary 4-point subdivision scheme is πΆ1 .
4. Properties of Subdivision Schemes In this section, we show that how limit curve of 4-point ternary and 4-point quaternary subdivision schemes give response to initial polynomial data. For this we discuss HΒ¨older regularity, degree of polynomial generation, polynomial reproduction and approximation order of the schemes (23) and (24). 4.1. HΒ¨older Regularity. HΒ¨older regularity is an extension of the notion of continuity which gives more information about any scheme. A function π : π
β π
is defined to be regular of order π¦ + πΌ (for π¦ β π0 and 0 < π β€ 1) if it is π¦ time continuously differentiable and ππ¦ is Lipschitz of order πΌ σ΅¨ σ΅¨σ΅¨ (π¦) σ΅¨σ΅¨π (π₯ + β) β π(π¦) (π₯)σ΅¨σ΅¨σ΅¨ β€ π|β|π (30) σ΅¨ σ΅¨ for all π₯ and β in π
and some constant π. According to Dyn and Levin [17] and Rioul [18], the HΒ¨older regularity of subdivision scheme with symbol π(π§) can be computed in the following way. Let π(π§) = π ((1 + π§ + β
β
β
+ π§πβ1 )/π) π(π§), without loss of generality we can assume π0 , . . . , ππ to be the nonzero coefficients of π(π§) and let π΅0 , π΅1 , . . . , π΅π be the π Γ π matrices with elements π, π = 1, . . . , π,
Theorem 9. The HΒ¨older regularity of scheme (23) is π = 4 β log3 (11) = 1.8173.
(29)
As we see β(1/3)π1 ββ < 1, then by Theorem 1 the scheme is πΆ0 . Similarly, β(1/3)π2 ββ < 1 then by Corollary 3 the scheme is πΆ1 .
(π΅π )ππ = ππ+πβππ+π ,
5
π = 0, 1, . . . , π. (31)
Then the HΒ¨older regularity is given by π = π β logπ (π), where π is the joint spectral radius of the matrices π΅0 , π΅1 , . . . , π΅π , that is, π = π (π΅0 , π΅1 , . . . , π΅π ) σ΅©1/π σ΅© = lim sup (max {σ΅©σ΅©σ΅©σ΅©π΅ππ . . . π΅π2 π΅π1 σ΅©σ΅©σ΅©σ΅©β : ππ β {0, 1}}) , πββ max {π (π΅0 ) , . . . , π (π΅π )} β€ π (π΅0 , . . . , π΅π ) σ΅© σ΅© σ΅© σ΅© β€ max {σ΅©σ΅©σ΅©π΅0 σ΅©σ΅©σ΅©β , . . . , σ΅©σ΅©σ΅©π΅π σ΅©σ΅©σ΅©β } . (32) Since π is bounded from below by the spectral radii and from above by the norm of the metrics π΅0 , π΅1 , . . .,π΅π , then σ΅© σ΅© σ΅© σ΅© max {π (π΅0 ) , . . . , π (π΅π )} β€ π β€ max {σ΅©σ΅©σ΅©π΅0 σ΅©σ΅©σ΅©β , . . . , σ΅©σ΅©σ΅©π΅π σ΅©σ΅©σ΅©β } . (33)
1 + π§ + π§2 π (π§) = ( ) π (π§) , 3
(34)
1 (β4 + 11π§ β 4π§2 ) . π§5
(35)
where π (π§) =
From (31) and (34), π0 = β4, π1 = 11, π2 = β4, π = 4, π = 2 and π = 3, thus π = 0, 1, 2, and then π΅0 , π΅1 , and π΅2 are the matrices with elements (π΅0 )ππ = π2+πβ3π , (π΅1 )ππ = π2+πβ3π+1 ,
(36)
(π΅2 )ππ = π2+πβ3π+2 , where π, π = 1, 2. This implies β4 0 π΅0 = ( ), 11 0
11 0 π΅1 = ( ), β4 0
β4 0 π΅2 = ( ) . (37) 0 β4
From (33) and (37) we have max {4, 11, 4} β€ π β€ max {11, 11, 4} .
(38)
Since the largest eigenvalue and the max-norm of the metrics are 11, so π = 4 β log3 (11) = 1.8173.
(39)
Theorem 10. The HΒ¨older regularity of scheme (24) is π = 4 β log4 (24) = 1.7077. Proof. The Laurent polynomial π(π§) of scheme (24) can be written as 4
π (π§) = (
1 + π§ + π§2 + π§3 ) π (π§) , 4
(40)
where π (π§) =
1 (β10 + 24π§ β 10π§2 ) . π§7
(41)
From (31) and (40), π0 = β10, π1 = 24, π2 = β10, π = 4, π = 2 and π = 4, thus π = 0, 1, 2, and then π΅0 , π΅1 , and π΅2 are the matrices with elements (π΅0 )ππ = π2+πβ4π , (π΅1 )ππ = π2+πβ4π+1 , (π΅2 )ππ = π2+πβ4π+2 ,
(42)
6
International Journal of Mathematics and Mathematical Sciences
where π, π = 1, 2. This implies β10 0 π΅0 = ( ), 24 0
24 0 π΅1 = ( ), β10 0
β10 0 π΅2 = ( ). 0 β10 (43)
From (33) and (43) we have max {10, 24, 10} β€ π β€ max {24, 24, 10} .
(44)
Thus the largest eigenvalue and the max-norm of the metrics are 24, so π = 4 β log4 (24) = 1.7077.
(45)
Remark 11. It is generally observed that as we decrease the arity of the scheme the HΒ¨older exponent increases. For example, Cashman et al. [19] have proved that the HΒ¨older exponent for binary scheme is 2 while from above theorems we see that HΒ¨older exponents for ternary and quaternary schemes are 1.8173 and 1.7077, respectively. Trivially, the HΒ¨older exponent approaches to 1 for large arity scheme. 4.2. Polynomial Generation. The generation degree of a subdivision scheme is the maximum degree of polynomials that can potentially be generated by the scheme, provided that the initial data is chosen correctly. Suppose π0 is polynomial of degree π of initial data ππ0 and symbol of the scheme is πβ1 π+1
π (π§) = (1 + π§ + β
β
β
+ π§
)
π (π§) ;
(46)
then the limit curve of the refined data πππ at any level π is polynomial of degree π. So the condition is necessary and sufficient for the scheme being able to generate polynomial of degree π. Theorem 12. The degree of polynomial generation of scheme (23) is 3. Proof. Since the Laurent polynomial π(π§) of the scheme (23) is 2 (3+1)
π (π§) = (1 + π§ + π§ )
π (π§) ,
1 (β4 + 11π§ β 4π§2 ) , (3)4 π§5
π (π§) ,
then degree of polynomial generation of scheme is 3.
π = 1, 2,
(51)
π
for π = 0, . . . , 3, πΌ3 = exp((2ππ/3)π) and π = πσΈ (1)/3. Proof. By taking first derivative of (25) and substituting π§ = 1 in it, we get π(1) (1) = 0.
(52)
This implies that π=
π(1) (1) = 0. 3
(53)
So from (12), the scheme (23) has primal parametrization. For π = 0, π = 1, and from (25), we get π(0) (πΌ31 ) = π (π2ππ/3 ) = 0.
(54)
Similarly, for π = 1, 2 and π = 0, 1, 2, 3 (π denotes the order of derivative) π
(55)
πβ1
(49)
where 1 π (π§) = (β10 + 24π§ β 10π§2 ) , (4)4 π§7
π=0
(48)
Proof. Since the Laurent polynomial of (24) can be written as π (π§) = (1 + π§ + π§ + π§ )
π
π(π) (πΌ3 ) = 0,
By (25), we get π(1) = 3. Also 3ββ1 π=0 (0 β π) = 3, which implies (πβπ). Similarly for π = 1, 2, 3, we can easily that π(1) = 3β0β1 π=0 show that
Theorem 13. The degree of polynomial generation of scheme (24) is 3.
3 (3+1)
πβ1
π(π) (1) = 3β (π β π) ,
π(π) (πΌ3 ) = 0.
then degree of polynomial generation is 3.
2
Theorem 14. A convergent subdivision scheme (23) reproduces polynomials of degree 3 with respect to the parameterizations (12) if and only if
(47)
where π (π§) =
4.3. Polynomial Reproduction and Approximation Order. The polynomial reproduction property has its own importance, as the reproduction property of the polynomials up to a certain degree π implies that the scheme has π + 1 approximation order. Polynomial reproduction of degree π requires polynomial generation of degree π. For this, polynomial reproduction can be made from initial data which has been sampled from some polynomial function. In the view [15], the polynomial reproduction property of the proposed scheme can be obtained after having the parameterizations π given in (12).
(50)
π(π) (1) = 3β (π β π) ,
(56)
π=0
which completes the proof. Since scheme (23) reproduces polynomial of degree 3, so by using Theorem 6, we get the following theorem. Theorem 15. A 4-point ternary interpolating scheme (23) has an approximation order of 4. Proof of the following theorem is similar to the proof of Theorem 14.
International Journal of Mathematics and Mathematical Sciences
(a) 4-point 2-ary
7
(b) 4-point 3-ary
(c) 4-point 4-ary
(d) 4-point 5-ary
(e) 4-point 6-ary
(f) 4-point 7-ary
Figure 4: Comparison of the limit curves generated by proposed 4-point 2-ary, 3-ary, 4-ary, 5-ary, 6-ary, and 7-ary interpolating subdivision schemes at 1st subdivision level.
Theorem 16. A convergent subdivision scheme (24) reproduces polynomials of degree 3 with respect to the parameterizations (12) if and only if πβ1
π(π) (1) = 4β (π β π) , π=0
π
π(π) (πΌ4 ) = 0,
π = 1, 2, 3
(57)
π
for π = 0, . . . , 3, πΌ4 = exp((2ππ/4)π) and π = πσΈ (1)/4.
by using imputation from Lagrange interpolation at four consecutive points [2]; therefore, 3-point ternary and 4-point quaternary schemes have polynomial reproduction of degree 3 and approximation properties are obvious (said by the referee). These schemes also have been generated by [4] using wavelet framework and by our algorithm, so by construction these schemes do no inherit these properties and that is only the reason to include the above theorems.
Again by Theorem 6, we get the following theorem. Theorem 17. A 4-point quaternary interpolating scheme (24) has an approximation order of 4. Remark 18. The considered schemes (i.e., 3-point ternary and 4-point quaternary) are exactly the same as obtained
5. Numerical Examples and Conclusion Six examples are depicted to show the usefulness of 4point 2-ary, 3-ary, 4-ary, 5-ary, 6-ary, and 7-ary interpolating subdivision schemes at 1st subdivision level in Figure 4. In this figure the control polygons are drawn by dotted lines
8 while the subdivision curves are drawn by solid lines. From Figure 4, it is clear that the initial polygon converges rapidly to limit curve as we increase the arity of the subdivision scheme. For many subdivision levels with any of these schemes the limit curves of the π-ary scheme with large π may exhibit sharper singularities at the initial control points compared to the schemes with small π (also mentioned by the referee). But if the initial control points come from noisy source, then π-ary scheme with large π is the better choice. The scheme with small π exhibit overfitting (see [7]). The main purpose to give comparison at first level is to provide the clear visual differences among the refined polygons produced by different schemes. In this paper, we have presented a multistep algorithm which generates 4-point π-ary interpolating subdivision schemes. We have also observed that the 4-point π-ary schemes generated by Lagrange polynomials and wavelet theory can also be generated by proposed multistep algorithm. Some significant properties like HΒ¨older regularity, degree of polynomial generation, degree of polynomial reproduction, and approximation order have been also discussed in this paper.
Acknowledgments This work was supported by Indigenous Ph.D. Scholarship Scheme of HEC Pakistan. The authors would like to thank the referees for their helpful suggestions and comments which show the way to improve this work.
References [1] S. Dubuc, βInterpolation through an iterative scheme,β Journal of Mathematical Analysis and Applications, vol. 114, no. 1, pp. 185β204, 1986. [2] G. Deslauriers and S. Dubuc, βSymmetric iterative interpolation processes,β Constructive Approximation, vol. 5, no. 1, pp. 49β68, 1989. [3] N. Dyn, D. Levin, and J. A. Gregory, βA π-point interpolatory subdivision scheme for curve design,β Computer Aided Geometric Design, vol. 4, no. 4, pp. 257β268, 1987. [4] J.-A. Lian, βOn 2π-ary subdivision for curve design. III. (2π + 1)-point and (2π+4)-point interpolatory schemes,β Applications and Applied Mathematics, vol. 4, no. 1, pp. 434β444, 2009. [5] G. Mustafa and N. A. Rehman, βThe mask of π-point π-ary subdivision scheme,β Computing, vol. 90, no. 1-2, pp. 1β14, 2010. [6] G. Mustafa, J. Deng, P. Ashraf, and N. A. Rehman, βThe mask of odd points π-ary interpolating subdivision scheme,β Journal of Applied Mathematics, vol. 2012, Article ID 205863, 20 pages, 2012. [7] G. Mustafa and I. P. Ivrissimtzis, βModel selection for the Dubuc-Deslauriers family of subdivision scheme,β in Proceedings of the 14th IMA Conference on Mathematics of Surfaces, University of Birmingham, Birmingham, UK, September2013. [8] J. M. Lane and R. F. Riesenfeld, βA theorical development for computer generation and display of piecewise polynomial surfaces,β IEEE Transaction on Pattern Analysis and Machine Intelligence, vol. 2, no. 1, pp. 35β46, 1980.
International Journal of Mathematics and Mathematical Sciences [9] E. Catmull and J. Clark, βRecursively generated B-spline surfaces on arbitrary topological meshes,β Computer Aided Geometric Design, vol. 10, no. 6, pp. 183β188, 1978. [10] D. Zorin and P. SchrΒ¨oder, βA unified framework for primal/dual quadrilateral subdivision schemes,β Computer Aided Geometric Design, vol. 18, no. 5, pp. 429β454, 2001, Subdivision algorithms (Schloss Dagstuhl, 2000). [11] P. Oswald and P. SchrΒ¨oder, βComposite primal/dual sprt(3) subdivision schemes,β Computer Aided Geometric Design, vol. 20, no. 3, pp. 135β164, 2003. [12] U. H. AugsdΒ¨orfer, N. A. Dodgson, and M. A. Sabin, βVariations on the four-point subdivision scheme,β Computer Aided Geometric Design, vol. 27, no. 1, pp. 78β95, 2010. [13] N. Aspert, Non-linear subdivision of univariate signals and Β΄ discrete surfaces [EPFL Thesis], Ecole Polytechnique FΒ΄edΒ΄erale de Lausanne, Lausanne, Switzerland, 2003. [14] M. F. Hassan and N. A. Dodgson, βTernary and three-point univariate subdivision schemes,β in Curve and Surface Fitting, A. Cohen, J. L. Marrien, and L. L. Schumaker, Eds., pp. 199β208, Nashboro Press, Sant-Malo, France, 2002. [15] C. Conti and K. Hormann, βPolynomial reproduction for univariate subdivision schemes of any arity,β Journal of Approximation Theory, vol. 163, no. 4, pp. 413β437, 2011. [16] N. Dyn, Tutorials on Multiresolution in Geometric Modelling, Summer School Lecture Notes, Mathematics and Visualization, Springer, Berlin, Germany, 2002. [17] N. Dyn and D. Levin, βSubdivision schemes in geometric modelling,β Acta Numerica, vol. 11, pp. 73β144, 2002. [18] O. Rioul, βSimple regularity criteria for subdivision schemes,β SIAM Journal on Mathematical Analysis, vol. 23, no. 6, pp. 1544β 1576, 1992. [19] T. J. Cashman, K. Hormann, and U. Reif, βGeneralized LaneRiesenfeld algorithms,β Computer Aided Geometric Design, vol. 30, no. 4, pp. 398β409, 2013.
Advances in
Decision Sciences
Advances in
Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
Mathematical Problems in Engineering Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
International Journal of
Combinatorics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
Journal of Applied Mathematics
Abstract and Applied Analysis Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
Volume 2013
Submit your manuscripts at http://www.hindawi.com Journal of Function Spaces and Applications
Discrete Dynamics in Nature and Society Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
International Journal of Mathematics and Mathematical Sciences
Advances in
Journal of
Operations Research
Probability and Statistics
International Journal of
International Journal of
Stochastic Analysis
Differential Equations Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
ISRN Mathematical Analysis Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
ISRN Discrete Mathematics Volume 2013
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
ISRN Applied Mathematics
ISRN Algebra Volume 2013
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013
ISRN Geometry Volume 2013
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2013