Research Article A Family of Even-Point Ternary Approximating Schemes

Transcription

1 International Scholarly Research Network ISRN Applied Mathematics Volume, Article ID 97, pages doi:.5//97 Research Article A Family of Even-Point Ternary Approximating Schemes Abdul Ghaffar and Ghulam Mustafa Department of Mathematics, The Islamia University of Bahawalpur, Bahawalpur, Pakistan Correspondence should be addressed to Ghulam Mustafa, ghulam.mustafa@iub.edu.pk Received 7 July ; Accepted September Academic Editors: L. Ju, F. Lamnabhi-Lagarrigue, G. Psihoyios, H. C. So, and X. Yang Copyright q A. Ghaffar and G. Mustafa. 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 presented a general formula to generate the family of even-point ternary approximating subdivision schemes with a shape parameter for describing curves. Some sufficient conditions for C to C 7 continuity and approximation order for certain ranges of parameter are discussed. The proposed even-point ternary schemes compare remarkably with existing even-point ternary schemes because they are able to generate limit functions with higher smoothness and approximation order. In addition, we measured curvature and torsion that assist the quality of subdivided curves.. Introduction and Preliminaries There are numbers of binary subdivision schemes in the literature. The interest in investigating arities higher than two has been started by Hassan et al.,. Nowadays, we have numbers of ternary schemes introduced by 7 and so forth, But the research communities are still gaining interest in introducing schemes higher than three arities i.e., quaternary, quinary, senary,..., n-ary. Mustafa and Khan introduced a new -point C quaternary approximating subdivision scheme. Lian 9, introduced -, -, 5-, and -point a-ary schemes. Lian also offered m-point and m-point interpolating a-ary schemes for curve design. The -, -,..., -point binary and ternary schemes are very common in the literature. The schemes involving convex combination of more or less than six points at coarse refinement level to insert a new point at next refinement level is introduced by Ko et al. 7. They introduced n - andn -point binary schemes. Zheng et al. investigated ternary interpolatory schemes with an odd number of control points, namely, n -point ternary interpolatory subdivision scheme. They also investigated ternary even symmetric

2 ISRN Applied Mathematics n-point and p-ary approximating subdivision scheme and presented the general ternary even symmetric n-point approximating subdivision rule and design alternative smooth ternary subdivision scheme of higher order. Mustafa and Rehman 5 presented general formulae for the mask of b -point n-ary approximating as well as interpolating subdivision schemes for any integers b andn. This motivates us to present the family of even-point ternary schemes with high smoothness and more degree of freedom for curve design. Proposed schemes not only provide the mask of even-point schemes but also generalize and unify several well-known schemes. Moreover, we measured curvature and torsion that can be used to describe the quality of curve. Also we compared plot of curvature and torsion, obtained by proposed schemes with the other existing schemes. A general compact form of univariate ternary subdivision scheme S, which maps polygon f k {f k i } i Z to a refined polygon f k {f k i } i Z is defined by f k i a j i f k j, i Z,. where the set a {a i : i Z} of coefficients is called the mask at k-th level of refinement. A necessary condition for the uniform convergence of subdivision scheme. is that j a j a a j.. A subdivision scheme is uniformly convergent if for any initial data f {f i : i Z}, there exists a continuous function f such that for any closed interval I R lim sup f k i k i k I ( f i k.. Obviously, f S f, introducing a symbol called Laurent polynomial az i Za i z i,. of the mask a {a i : i Z} which plays an efficient role to analyze the convergence and smoothness of subdivision scheme. From. and. the Laurent polynomial of convergent subdivision scheme satisfies ( a e iπ/ ( a e iπ/, a..5 This condition guarantees existence of a related subdivision scheme for the divided difference of the original control points and the existence of associated Laurent polynomial a z a z z z z az..

3 ISRN Applied Mathematics The subdivision scheme S with Laurent polynomial a z is related to scheme S with Laurent polynomial az by the following theorem. Theorem. see. Let S denote a subdivision scheme with Laurent polynomial az satisfying.5. Then there exists a subdivision scheme S with the property Δf k S Δf k,.7 where f k S k f and Δf k {Δf k i k f k i f k i ; i Z}. Furthermore, S is a uniformly convergent if and only if /S converges uniformly to zero function for all initial data f,inthe sense that ( k lim k S f.. A scheme S satisfying. for all initial data f is termed contractive. The above theorem indicates that for any given scheme S, with mask a satisfying., we can prove the uniform convergence of S by deriving the mask of /S and computing /S i for i,,..., L, where L is the first integer for which /S L <. If such an L exists, then S converges uniformly. Since there are three rules for computing the values at next refinement level, so we define the norm as follows: S max aj, aj, aj, ( L n S max b n,l i L j ; i,,,..., L,.9. where b n,l z L a n( z j, L. j ( ( n z a n z z a n z az, n.. z z z z The paper is organized as follows: the family of even-point ternary approximating scheme and analysis of two even-point ternary schemes are presented in Section. Basic properties of even-point ternary schemes are discussed in Section. Comparison with existing even-point ternary schemes is shown in Section. A few remarks and future work constitute Section 5.

4 ISRN Applied Mathematics. The Even-Point Ternary Approximating Schemes Here we offer a general formula for even-point ternary approximating subdivision schemes with one parameter in the form of Laurent polynomial a λ z ( z z λn λn {( ( ( ( } 5 5 ω ω z ω z z ω,. where ω is shape parameter, λ n,n and depends on n. Despite the fact that we can generate λ-point ternary approximating schemes from. for any λ n, n, however for straightforwardness, we generated and discussed the smoothness of - and -point ternary schemes. The general proof regarding the support of even-point ternary schemes has been presented in this paper... -Point Ternary Scheme From. for n, we get the following: Laurent polynomial for a -point ternary scheme a z {( 7 ω z ( 5 ω ( 9 ω ( z z ( 5 ( 5 ω z ( 7 9ω z 5 5ω z 9 ( 9 5ω z ( 7 9ω z ( ω z ω z ( 5 ω z 7 } z ( ω.. From the above polynomial, we suggest the following -point ternary approximating scheme: f k i [( ( 5 7 5ω f ki 7 9ω f k i f k i [( ( 5 ω f k i 7 ω f k i f k i [( ( 9 7 ω f k i ω f k i From. and., we have a ( ω f k i ω ( f ω ki ω ( 5 f k i ( 9 ( 5 ( 7 9ω {( ( ( z z 9 ω ω z ( ( 7 7 ω z ω z 5 ( ( } ω z z ω 9, ( 5 ω z ( 5 ω ω ( z f k i f k i ], ], ] 5ω f k i. z ω z 7.

5 ISRN Applied Mathematics 5 a {( z z ω z ( 7 ω a z ( ω z ( ω {( a z z ω ( ( 7 7 ω z z ω z ( ( } 7 z z 5 ω z z ω 7, ( z ω z ( z ω z 5, } z z z. ( 5 ω ( 5 ω ( ω. If S m is the scheme corresponding to a m, then for Cm continuity, we require that a m satisfies., which it does and /S m L <, for m,,,. Since by., for 7/ <ω<9/, / <ω</, / <ω</, and 5/ <ω<7/, we have { S max 7 ω 7 5 ω, 7 ω 7 7 ω 7 } ω <, { S max 9 ω 9 7 ω 9 7 ω, } <, { S max ω ω, } <, { S max ω, 5 } ω, 5 ω <..5 In order to obtain a C scheme, we substitute ω / in a z, we have a ( z z z z z.. Using.,weget ( a 5 z z a z z z z { z}..7 If S m is the scheme corresponding to a m, then for Cm continuity, we require that a m satisfy., which it does and /S m L <, for m,,,.

6 ISRN Applied Mathematics For m 5andL, we have S 5 max aj, aj, aj,. max {,, } <. Therefore /S 5 is contractive, that is, /S 5 < and the scheme is C... -Point Ternary Scheme From. for n, we get the following -point ternary approximating scheme: f k i [( 7 79 ( 7 ω f k i ( 59 f k i [ ( w f k i 79 ( f k i [( 79 ( 5 7ω f k i ( 7 5ω f k i 7ω ( w f ki ω f k i f k i ( ω f k i 9ω f k i ( ] 9ω f k i ω f k i ωf k i, ( ( 5 59 ω f k i ω f k i 7ω f k i ( ( ( ω f k i 5ω f ki ω ], f k i ]..9 The continuity of this scheme can be computed in a similar way as we did for the -point ternary scheme. From above discussion, we reach the conclusion shown in Table. Remark... i For n, ω /9u, it becomes -point ternary subdivision scheme ii If we set n, ω / 9w, it becomes -point ternary subdivision scheme. iii In case n, ω / /7u, we get mask of -point ternary subdivision scheme.. Basic Properties of the Schemes In this section, we discuss approximation order and support of basic limit function of evenpoint ternary approximating schemes.

7 ISRN Applied Mathematics 7 Table : The order of continuity of proposed -point and -point ternary schemes for certain ranges of parameter. Scheme Parameter Continuity Scheme Parameter Continuity 7 9 <ω< C 9 9 <ω< C <ω< C <ω< C -point <ω< 5 <ω< 7 ω C 9 C -point 7 C <ω< 75 <ω< 9 <ω< C C C <ω< 5 <ω< 7 ω C 5 C C 7.. Approximation Order Here we only find the approximation order of proposed -point ternary scheme. The approximation order of other even-point ternary schemes can be computed in a similar fashion. Theorem.. A -point ternary approximating scheme has approximating order for ω 5/ < ω<7/ and 5 for ω 5/. A -point ternary approximating scheme has approximating order 7 for ω 5/ <ω< 7/ and for ω /. Proof. We carry out this result by taking our origin the middle of an original span with ordinate..., 5 n, n, n, n, n, 5 n,... If y x n, then we have [ y ]...,a 5 n a n a n a n a 5 5 n a n a n a 5 n a 5 n a n a n a n a n a n a n a n a 5 n a n a n a 5 n. a 5 n a n a n a 5 5 n a n a n a n a 5 n,...,

8 ISRN Applied Mathematics where a /75/ 5ω,a /77/ 9ω,a /79/ ω,a /7/ ω,a 5 /7/ ω,a /75/ ω. If y x, then [ ] y...,,,,,,,,,,... [ ] δy...,,,,,,,,..., [ ] δ y,. where δ represents the differences of the vertices. If y x, then [ y ]..., ω, ω, ω, ω, ω, ω,... 9 ω, 9 9 ω, ω,. Taking further differences, we get δ y. If y x, then [ y ]..., ω, 5 5 ω, ω, ω, 7 7 ω,... 7 ω, ω,, 7 9 ω,. This implies that δ y. If y x, then [ y ]..., ω, ω, ω, 9 7 ω, ω, ω,... 9 ω, ω, ω,.5 By taking further differences, we have [ ] δ 5 y..., 7 ω, 9 7 ω,, 9 7 ω, ω,..., [ ] δ 5 y, at ω 5.. Since by, -point ternary scheme has cubic precision for ω 5/ <ω<7/ and quartic at ω 5/ then -point ternary scheme has approximating order for ω 5/ < ω<7/ and 5 for ω 5/. The proof of other part is similar.

9 ISRN Applied Mathematics a -point scheme b -point scheme Figure : The basic limit functions of proposed schemes at ω /... Support of Basic Limit Function The basic limit functions, ϕ λ S λ f i of proposed λ-point ternary approximating schemes, for λ n, n, are presented in Figures a and b. The following theorem is related to the support of basic limit functions of even point ternary schemes. Theorem.. The basic limit functions ϕ λ of proposed λ -point ternary approximating schemes have support width s λ n /, whereλ n, n, which implies that it vanishes outside the interval λ n/, λ n/. Proof. The support width s of the basic limit functions can be determined by computing howfortheeffect of the nonzero vertex f will propagate along by. As the mask of λ scheme is λ -long sequence by centering it on that vertex, the distances to the last of its left and right nonzero coefficients are equal to λ n andλn, respectively. At the first subdivision step we see that the vertices on the left and right sides of f at λ n / and λ n / are the furthest nonzero new vertices. At each refinement, the distances on both sides are reduced by the factor /. At the next step of the scheme this will propagate along by λn/ /ontheleftandλn/ / on the right. Hence after k subdivision steps the furthest nonzero vertex on the left will be at λ n / / / k λ n / k j /j and on the right will be at λ n / / / k λ n / k j /j. So the total support width is λ n / j /j λ n / j /j λ n/.. Comparison and Application In order to show the performance of the proposed schemes, we compare continuity, support, approximation order, and shape of limit curves. We also discuss curvature and torsion.

10 ISRN Applied Mathematics a -point b -point c -point d -point e -point Figure : Comparison: bold solid continuous curves are generated by proposed -point and -point ternary approximating schemes a Ko et al. 7 and proposed -point ternary schemes b 5, 9 and proposed - point ternary scheme c Hassan and Dodgson and proposed -point ternary schemes d 5, 9 and proposed -point ternary scheme e Khan and mustafa and proposed -point ternary schemes. Table : Comparison of -point ternary schemes. Scheme Type Support Order C n Ternary -point Interpolating 5 Ternary -point Interpolating 5 Ternary -point Interpolating 5 Ternary -point 5, 9 Interpolating 5 Ternary -point 7 Approximating 5.5 Ternary -point Approximating 5 Ternary -point 5 Approximating 5.5 Ternary -point proposed Approximating Comparison of various - and -point ternary schemes are given in Tables and. Figure shows visual comparison of - and -point ternary interpolating and approximating schemes with the proposed ternary approximating schemes... and Torsion The quality of subdivided curves can be assessed quantitatively by measuring curvature and torsion, as functions of cumulative chord length. is the amount by which a geometric object deviates from being flat, or straight in the case of a line, while in the elementary geometry of space curves, torsion measures the rate at which a twisted curve tends to depart from its osculating plane. When the torsion is zero, osculating plane never changes, and we have a plane curve. We used the method described in to determine the curvature and torsion. The appearance of curvature and torsion obtained by our schemes

11 ISRN Applied Mathematics -point [, ].... of [, ] -point [].... of [] a b c d point [5, 9] of [5, 9] -point [7] of [7] e f g h -point [].... of [] -point proposed.... of proposed i j k l Figure : Comparison of existing and proposed -point schemes sampled from control polygon by using the 5th iteration. The results are shown on the left together with their corresponding curvature on the right. Table : Comparison of -point ternary schemes. Scheme Type Support Order C n Ternary -point 5, 9 Interpolating 5 Ternary -point Interpolating Ternary -point proposed Approximating.5 7

12 ISRN Applied Mathematics point [, ] a of [, ] b Torsion.... Torsion of [, ] c point [] of [] Torsion.... Torsion of [] d e f point [5, 9] of [5, 9] Torsion.... Torsion of [5, 9] g h i.... -point [7] of [7] Torsion.... Torsion of [7] j k l Figure : Continued.

13 ISRN Applied Mathematics point proposed m of proposed n Torsion.... Torsion of proposed o Figure : Left: control data sampled from a space curve and four iteration of -point subdivision scheme. The corresponding curvature and torsion are shown in the center and right column, respectively. is measured and compared with the plot or appearance generated by other schemes. The examples in Figure show the results of curvature and torsion, together with their control polygon. The appearance of curvature of the schemes,, 5, 9 is quite similar and plot of curvature of the schemes, 7, is different. Figure shows detailed view of space curve, curvature and torsion generated by existing schemes 5, 7, 9 and proposed schemes. Overall one can observe that the plot of space curve, curvature, and torsion of proposed schemes are more smoothing than other existing schemes. 5. Conclusion The family of even-point approximating schemes for curve design has been established. The - and -point ternary schemes introduced by Zheng et al., are special cases of our proposed even-point approximating schemes. Smoothness and approximation order of proposed - and -point ternary schemes have been discussed. Support of family of even-point ternary schemes has been computed in general. It has been shown that proposed schemes are much better than existing ternary schemes in the sense of smoothness, approximation order, curvature, and torsion. The family of odd-point ternary approximating schemes will be studied in detail in the forthcoming paper. Acknowledgment This work is supported by the Indigenous Ph.D Scholarship Scheme of Higher Education Commission HEC Pakistan. References M. F. Hassan and N. A. Dodgson, Ternary and three-point univariate subdivision schemes, in Curve and Surface Fitting: Saint-Malo,, A. Cohen, J. L. Marrien, and L. L. Schumaker, Eds., pp. 99, Nashboro Press, Brentwood, Calif, USA,. M. F. Hassan, I. P. Ivrissimitzis, N. A. Dodgson, and M. A. Sabin, An interpolating -point C ternary stationary subdivision scheme, Computer Aided Geometric Design, vol. 9, no., pp.,.

14 ISRN Applied Mathematics C. Beccari, G. Casciola, and L. Romani, An interpolating -point C ternary non-stationary subdivision scheme with tension control, Computer Aided Geometric Design, vol., no., pp. 9, 7. C. Beccari, G. Casciola, and L. Romani, Shape controlled interpolatory ternary subdivision, Applied Mathematics and Computation, vol. 5, no., pp. 9 97, 9. 5 G. Deslauriers and S. Dubuc, Symmetric iterative interpolation processes, Constructive Approximation, vol. 5, no., pp. 9, 99. F. Khan and G. Mustafa, Ternary six-point interpolating subdivision scheme, Lobachevskii Mathematics, vol. 9, no., pp. 5,. 7 K. P. Ko, B.-G. Lee, and G. J. Yoon, A ternary -point approximating subdivision scheme, Applied Mathematics and Computation, vol. 9, no., pp. 5 57, 7. G. Mustafa and F. Khan, A new -point C quaternary approximating subdivision scheme, Abstract and Applied Analysis, vol. 9, Article ID 97, pages, 9. 9 J.-A. Lian, On a-ary subdivision for curve design. I. -point and -point interpolatory schemes, Applications and Applied Mathematics, vol., no., pp. 9,. J.-A. Lian, On a subdivision for curve design. II. -point and 5-point interpolatory schemes, Applications and Applied Mathematics, vol., no., pp. 7 7,. J.-a. Lian, On a-ary subdivision for curve design. III. m-point and m -point interpolatory schemes, Applications and Applied Mathematics, vol., no., pp., 9. H. Zheng, M. Hu, and G. Peng, Constructing n--point ternary interpolatory subdivision schemes by using variation of constants, in Proceedings of the International Conference on Computational Intelligence and Software Engineering (CISE 9. H. Zheng, M. Hu, and G. Peng, Ternary even symmetric n-point subdivision, in Proceedings of the International Conference on Computational Intelligence and Software Engineering (CISE 9. H. Zheng, M. Hu, and G. Peng, P-ary subdivision seneralizing B-splines, in Proceedings of the nd International Conference on Computer and Electrical Engineering (ICCEE 9, 9. 5 G. Mustafa and N. A. Rehman, The mask of b -point n-ary subdivision scheme, Computing. Archives for Scientific Computing, vol. 9, no. -, pp.,. M. A. Sabin and N. A. Dodgson, A circle-preserving variant of the four-point subdivision scheme, in Mathematical Methods for Curve and Surfaces. Troms, M. Dhlen, K. Mrken, Schumaker, and L, Eds., Modern Methods in Mathematics, pp. 75, Nashboro Press, Brentwood, Calif, USA, 5.

15 Advances in Operations Research Volume Advances in Decision Sciences Volume Applied Mathematics Algebra Volume Probability and Statistics Volume The Scientific World Journal Volume International Differential Equations Volume Volume Submit your manuscripts at International Advances in Combinatorics Mathematical Physics Volume Complex Analysis Volume International Mathematics and Mathematical Sciences Mathematical Problems in Engineering Mathematics Volume Volume Volume Volume Discrete Mathematics Volume Discrete Dynamics in Nature and Society Function Spaces Abstract and Applied Analysis Volume Volume Volume International Stochastic Analysis Optimization Volume Volume