Generating Functions for the q-bernstein Bases

Generating Functions for the q-bernstein Bases Ron Goldman Department of Computer Science, Rice University Houston, Texas 77251 USA Plamen Simeonov Department of Computer and Mathematical Sciences University of Houston-Downtown Houston, Texas 772 USA Yilmaz Simse Department of Mathematics, Adeniz University Antalya, Turey January 19, 214 Abstract We derive explicit formulas for the generating functions of the q-bernstein basis functions in terms of q-exponential functions. Using these explicit formulas, we derive a collection of functional equations for these generating functions which we apply to prove a variety of identities some old and some new for the q-bernstein bases. Keywords and phrases: q-bernstein bases, generating functions, q-exponential functions, q-calculus, q-bézier curves 21 Mathematics Subject Classification: 65Qxx, 33Dxx, 65D17, 26Cxx, 3C1. File name: G-S-S-qGF-January-214. 1 Introduction The Bernstein bases play a central role in the theory of Bézier curves and surfaces. These polynomial curves and surfaces are of fundamental importance in Computer Aided Geometric Design (CAGD) for modeling free-form shapes [6]. Recently Phillips and Oruç [1], [11], [12] advocated using the q-bernstein basis functions as a generalization of the classical Bernstein basis functions for CAGD. Building on this wor of Phillips and Oruç on q-bernstein bases and q-bézier curves as well as on the seminal papers of Ramshaw on blossoming [13], [14], Simeonov et al [15] introduced q-blossoming to investigate the analytic and geometric properties of q- Bernstein bases and q-bézier curves. Simeonov et al also discuss several potential applications of q-bézier curves and surfaces to CAGD. In order to investigate the properties of classical Bernstein bases, Simse [17], [18] defined a novel collection of generating functions to derive many nown and some new identities for the 1

2 classical Bernstein basis functions. The goal of this paper is to construct analogous generating functions for the q-bernstein basis functions. By using an approach similar to that of Simse as well as some nown properties of basic hypergeometric series [2], [5], [8], we shall derive some nown and some new identities and formulas for the q-bernstein basis functions, including the partition of unity property, formulas for representing the monomials, a convolution formula with the monomials, recurrence relations, formulas for q-derivatives, subdivision identities, and a q-marsden identity. We shall also derive for each fixed index two new closed formulas for the sum of the q-bernstein basis functions of degree n scaled by q n. In Section 2 we recall the definition of the classical and the q-bernstein basis functions along with some basic notation from the q-calculus. We also introduce the main generating functions discussed in this paper. In Section 3 we briefly review the definition and some basic properties of the q-exponential function, the q-derivative, and the q-integral, and state the main functional equation for the generating functions defined in Section 2. Section 4 contains most of our main results. Here we use the results from Section 3 concerning the q-exponential function and the q-calculus to establish a collection of functional equations for the generating functions from Section 2, and we apply these functional equations to derive identities for the q-bernstein bases and for their discrete q-derivatives. We also evaluate two discrete integrals for the q-bernstein bases. Section 5 contains several auxiliary definitions and results from the theory of basic hypergeometric series, the proofs of the properties of the q-exponential function listed in Section 3, and the proof of the functional equation in Section 3 for the generating functions from Section 2. In Section 6, we derive for each fixed index two novel closed formulas for the sum of the q-bernstein basis functions of degree n scaled by q n, using a new generating function introduced in Section 5. We close in Section 7 with a brief review of our wor along with a table summarizing the main identities and generating functions derived in this paper. 2 The q-bernstein Bases and their Generating Functions The classical Bernstein basis functions of degree n on the interval [, 1] are defined by ( ) n B n (x) = x (1 x) n, =,..., n. (2.1) Simse and Acigoz [19] and Acigoz and Araci [1] introduced the generating functions g (x, t) = B n (x) tn n! = x t e (1 x)t. (2.2)! n= Using these generating functions, Simse proved many identities for the classical Bernstein basis functions [17], [18]. The q-bernstein basis functions over the interval [, 1] are defined by [ n B n (x; q) = ] q x n 1 j= (1 q j x), =,..., n. (2.3)

Here [ n] denotes the q-binomial coefficient defined by q where and for Z +, [] q is the q-integer defined by [ n ] = [n] q!, =,..., n (2.4) q [] q![n ] q! [] q! = [1] q [2] q [] q, q 1, [] q! := 1 (2.5) 3 [] q = 1 + q + + q 1 = { 1 q 1 q, q 1,, q = 1. (2.6) We shall have occasions to consider both [] q and [] 1/q. It follows from (2.5) and (2.6) that [] q = q 1 [] 1/q (2.7) and [] q! = q ( 2) []1/q!. (2.8) The q-bernstein basis functions were introduced by Phillips [11] and Oruç [1], who advocated their use in CAGD. q-bernstein basis functions on arbitrary intervals were introduced in [9], and later in [15] as blending functions for q-bézier curves over arbitrary intervals. Other types of q-bernstein polynomials have been studied, see for example [19] for a type used in interpolation theory. The classical Bernstein polynomials B n (x) are obtained as the limiting case q 1 of the q-bernstein polynomials in [19]. In analogy with the generating functions g (x, t) in equation (2.2) for the classical Bernstein basis functions, we define generating functions for the q-bernstein basis functions by setting G (x, t; q) = n= B n tn (x; q) [n] q!. (2.9) We shall see shortly that the series for G (x, t; q) converges for all t < 1/ 1 q if q < 1, and for all t C if q > 1 or q = 1. Therefore for any fixed value of q >, there is an interval for t in which the series for G (x, t; q) converges. 3 The q-exponential Function and the q-calculus Just as the generating functions g (x, t) in equation (2.2) for the classical Bernstein basis functions are closely related to the classical exponential function e x, so too the generating functions G (x, t; q) for the q-bernstein basis functions defined in equation (2.9) are closely related to the q-exponential function E q (x). The q-exponential function E q (x) is defined by [4] E q (x) = n= x n [n] q!. (3.1)

4 This function is closely related to the standard q-exponential functions-see Section 5. By the ratio test this series for E q (x) converges and therefore E q (x) is well defined for all x < 1/ 1 q if q < 1, and for all x C if q > 1 or q = 1. Therefore for any fixed value of q >, there is an interval in which the series for E q (x) converges. The continuous derivatives of the classical exponential functions are easy to compute. Similarly the discrete q-derivatives of the q-exponential function have an analogous simple form. To understand what we mean by the discrete q-derivatives of the q-exponential function, we need to tae a quic detour into the q-calculus. 3.1 The q-calculus The discrete q-derivative of a function f(x) is defined by [2], [5], [8] D q,x f(x) = f(qx) f(x) (q 1)x, q 1. Notice that if f (x) exists, then f (x) = lim q 1 D q,x f(x). The q-derivative is a linear operator. In addition there is a simple product rule and chain rule for the q-derivative [8] and for any constant c D q,x (u(x)v(x)) = D q,x (u(x))v(x) + u(qx)d q,x (v(x)) (3.2) D q,x (f(cx)) = c D q,y f(y) y=cx. (3.3) These two formulas are easy to derive directly from the definition of the q-derivative. For polynomials the q-derivatives are easy to compute. Indeed, it follows easily from the definition of the q-derivative and equation (2.7) that D q,x (x ) = [] q x 1, (3.4) D 1/q,x (x ) = [] q (x/q) 1. (3.5) Just as there is a discrete form of differentiation, there is also a discrete form of integration. The discrete q-integral of a function f(x) is defined by [2], [5], [8] t f(x) d q x = t(1 q) f(tq j )q j, (3.6) provided that the infinite sum converges. Since the q-integral is an infinite Riemann sum with the division points in geometric progression, the q-integral converges to the classical integral as q 1 for functions that are continuous on [, t]. The q-integral has many properties analogous to the properties of the classical integral. For example, the q-integral is a linear operator. Moreover if f(x) is continuous at, then j= In addition, by (3.6) t t D q,x f(x) d q x = f(t) f(). (3.7) f(cx) d q x = 1 c ct f(x) d q x. (3.8)

5 For polynomials, q-integrals are easy to compute. Indeed, it follows easily from equations (3.7), (3.4), and (3.5) that t t x n d q x = x n d 1/q x = tn+1 [n + 1] q, (3.9) tn+1 = qn t n+1. (3.1) [n + 1] 1/q [n + 1] q 3.2 Properties of the q-exponential Function Returning to the q-exponential functions, we now have the following simple formulas for the q- and the 1/q-derivatives: D q,x E q (x) = E q (x), (3.11) D q,x (1/E q (x)) = 1/E q (qx), (3.12) D 1/q,x (1/E q (x)) = 1/E q (x). (3.13) Equation (3.11) follows directly from (3.1) and (3.4). Equations (3.12) and (3.13) will be proved in Section 5. In Section 5 we shall also derive the following identities: E 1/q (x) = 1/E q ( x), (3.14) E q (qx) = (1 (1 q)x)e q (x). (3.15) We close this section with an explicit formula for the generating functions G (x, t; q). Theorem 3.1. Let max{ q, x } < 1 and t < 1/ 1 q. Then Proof. See Section 5. G (x, t; q) = (xt) [] q! E q (t) E q (xt). (3.16) 4 Identities for the q-bernstein Basis Functions from Functional Equations for their Generating Functions In this section we will use our generating functions to derive properties of the q-bernstein basis functions. We shall proceed in the following fashion. In each case we first derive a functional equation satisfied by the generating functions. We then extract an identity for the q-bernstein basis functions by comparing the coefficients of t n on both sides of this functional equation. All of the identities in this section can also be verified using standard summation formulas for basic hypergeometric series. The formulas for the basis functions are typically valid for all values of q, even though there are constraints on where (for what values of t and q) the functional equations for the generating functions are valid. The reason these identities hold for the q-bernstein basis functions for all values of q is that the identities for the functional equations hold for all values

6 of x [, 1/(1 q)] and all values of t [, 1] for < q < 1. Thus the identities for the q-bernstein basis functions are valid for all values of x [, 1/(1 q)] for < q < 1. But the q-bernstein basis functions B n (x; q) are also polynomials in q, and polynomial identities that are valid over an interval are valid everywhere. Thus it is enough to derive identities for the functional equations for < q < 1 and x [, 1/(1 q)]. Our functional equations hold for all values of x [, 1/ 1 q ] and all values of t [, 1] for q < 1. Using Theorem 3.1 and equation (3.1), we immediately obtain the following result. Theorem 4.1. Let max{ q, x, z } < 1 and t < 1/ 1 q. Then = G (x, t; q)z = E q(t) E q (xt) E q(zxt). (4.1) In particular, G (x, t; q) = E q (t). (4.2) = From Theorem 4.1, we can easily show that the q-bernstein bases form a partition of unity. Corollary 4.2. n B n (x; q) = 1. (4.3) = Proof. Changing the order of summation in (4.2) yields n= t n [n] q! n B n (x; q) = E q(t) = = n= t n [n] q!. Now (4.3) follows by equating the coefficients of t n [n] q! on both sides of this equation. A convolution identity between the monomials and the q-bernstein basis functions can also be derived from the functional equation (3.16) in Theorem 3.1. Corollary 4.3. The q-bernstein basis functions satisfy the convolution identity Proof. By (3.16) we have n [ n ] [ n ] l q xl B n l (x; q) = q x. (4.4) l= E q (xt)g (x, t; q) = (xt) [] q! E q(t). Expanding E q (xt), G (x, t; q), and E q (t) in this equation using (3.1) and (2.9), and then equating the coefficients of t n /[n] q! on both sides of the resulting equation yields (4.4).

7 Theorem 4.4. Let max{ q, x, z } < 1 and t < 1/ 1 q. Then [ ] G (x, t; q)z j = (xt)j [j] q! j [ ] j =j =j q q E q (t) E q (xt) E q(zxt). (4.5) G (x, t; q) = (xt)j [j] q! E q(t). (4.6) Proof. Applying D j q,z to (4.1), then applying (3.4), (3.11), and (3.3), and finally dividing both sides by [j] q!, we obtain (4.5). Setting z = 1 in (4.5) yields (4.6). Corollary 4.5. n =j [ ] B n (x; q) = j q [ ] n x j. (4.7) j q Proof. We use (2.9) and (3.1) to expand both sides of (4.6) in powers of t, and then we equate the coefficients of t n /[n] q! on both sides of the resulting equation to arrive at (4.7). Corollary 4.6. j ( 1) n q (n 2) n= [ ] j B n n (x; q) = 2) q ( 1/q [ ] j x j. (4.8) 1/q Proof. From (3.16) and (3.14) we get G (x, t; q)e 1/q ( t) = (xt) [] q! E 1/q( xt). We expand the functions on both sides of this equation using (2.9) and (3.1) to obtain n= B n (x; q) tn [n] q! l= ( t) l [l] 1/q! = (xt) [] q! Equating the coefficients of t j on both sides of this equation yields j n= B n j n ( 1) (x; q) [n] q![j n] 1/q! = l= ( xt) l [l] 1/q!. ( 1)j x j [] q![j ] 1/q!. The last equation reduces to (4.8) after multiplying both sides by ( 1) j [j] 1/q! and applying (2.8). The partial q-derivatives of the generating functions also contain information about the q-bernstein basis functions. We begin with two PDE s for the q-derivatives of the generating functions with respect to x. These PDE s lead to formulas for the q-derivatives of the q- Bernstein basis functions. Theorem 4.7. Let max{ q, x } < 1 and t < 1/ 1 q. Then D q,x G (x, t; q) = tq +1 G 1 (qx, t; q) tq G (qx, t; q), (4.9) D 1/q,x G (x, t; q) = tq +1 G 1 (x, t; q) tq G (x, t; q). (4.1)

8 Proof. From (3.16), (3.2) (with u(x) = 1/E q (xt) and v(x) = x ), (3.12), (3.3), and (3.4), it follows that { } { D q,x G (x, t; q) = t E q (t) x D q,x = t E q (t) x t [] q! E q (xt) [] q! E q (qxt) + [] qx 1 } E q (qxt) which proves (4.9). The proof of (4.1) is similar. Again using (3.16), (3.2) (with u(x) = x and v(x) = 1/E q (xt)), (3.13), (3.3), and (3.5), we obtain { } { D 1/q,x G (x, t; q) = t E q (t) x D [] q! 1/q,x = t E q (t) []q (x/q) 1 } (x/q) t E q (xt) [] q! E q (xt) E q (xt) which yields (4.1). From definition (2.9) of the generating functions G (x, t; q) and formulas (4.9) and (4.1), we can derive differentiation formulas for the q-bernstein basis functions by equating the coefficients of t n on both sides of these equations. Corollary 4.8. D q,x B n (x; q) = [n] q q ( qb 1 n 1 (qx; q) Bn 1 (qx; q) ), (4.11) D 1/q,x B n (x; q) = [n] qq ( qb n 1 1 (x; q) Bn 1 (x; q) ). (4.12) Formula (4.12) is a special case of formula (2.2) in [7]. As a direct consequence of formulas (4.11) (4.12), we evaluate two discrete integrals for the q-bernstein basis functions. Corollary 4.9. If q < 1, then If q > 1, then B n (x; q) d q x = B n (x; q) d 1/qx = { q +1 [n+1] q, =,..., n 1, 1 [n+1] q, = n. (4.13) q [n + 1] q, =,..., n. (4.14) Proof. First we prove (4.13). To begin, observe that by (2.3) and (3.9) B n n(x; q) d q x = x n d q x = 1 [n + 1] q. (4.15) Also by q-integrating (4.11) (with n replaced by n + 1) from to 1/q and applying (3.7) on the left hand side and (3.8) on the right hand side, we get B n+1 (1/q; q) B n+1 (; q) = [n + 1] q q Notice that by (2.3) B n+1 (; q) =, 1, B n+1 (1/q; q) = ( B n 1 (x; q) q 1 B n (x; q)) d q x. (4.16) { [n + 1]q q n (1 1/q), = n,, < n.

9 Therefore when = n, (4.16) yields [n + 1] q q n (1 1/q) = [n + 1] q q n ( B n n 1 (x; q) q 1 Bn(x; n q) ) d q x which together with (4.15) implies Bn n 1 (x; q) d qx = q n /[n + 1] q. When < < n, (4.16) becomes B 1 n (x; q) d qx = q 1 B n(x; q) d qx. Iterating this relation n times yields B 1 n (x; q) d qx = q (n ) Bn 1(x; n q) d q x = q /[n + 1] q. Next, we prove (4.14). To begin, observe that by (2.3) and (3.1) B n n (x; q) d 1/qx = x n d 1/q x = qn [n + 1] q. (4.17) Now we 1/q-integrate (4.12) (with n replaced by n + 1) from to 1 and apply (3.7) and (2.3). When < < n + 1, we get ( qb n 1 (x; q) B n (x; q)) d 1/q x = q [n + 1] q ( B n+1 (1; q) B n+1 (; q) ) =, hence Bn 1 (x; q) d 1/qx = q 1 Bn (x; q) d 1/qx. Iterating this relation n +1 times yields B 1 n (x; q) d 1/qx = q (n +1) Bn(x; n q) d 1/q x = q 1 [n + 1] q. Next we derive two more PDE s for the generating functions G (x, t; q), this time by taing the partial q-derivatives with respect to t. These new PDE s for the generating functions G (x, t; q) lead to recurrences for the q-bernstein basis functions. Theorem 4.1. Let max{ q, x } < 1 and t < 1 q. Then D q,t G (x, t; q) = xg 1 (x, t; q) + q G (x, t; q) xg (x, qt; q), (4.18) D q,t G (x, t; q) = xq +1 G 1 (x, qt; q) + G (x, t; q) xq G (x, qt; q). (4.19) Proof. First we prove (4.18). We multiply both sides of equation (3.16) by E q (xt) and then we apply D q,t to the resulting equation. Using the product rule (3.2) with u(t) = G (x, t; q) and v(t) = E q (xt) on the left hand side and with u(t) = t and v(t) = E q (t) on the right hand side yields D q,t (G (x, t; q))e q (xt) + G (x, qt; q)xe q (xt) = x ( ) [] q t 1 E q (t) + (qt) E q (t) [] q!

1 where we also used (3.3), (3.11), and (3.4). Dividing both sides of this equation by E q (xt) and applying (3.16), we get (4.18). To establish (4.19), we apply D q,t to equation (3.16) and use the product rule (3.2) twice on the right hand side, first with u(t) = E q (t)/e q (xt) and v(t) = t, and then with u(t) = 1/E q (xt) and v(t) = E q (t). We then get = (xt) [] q! D q,t G (x, t; q) = x [] q! = (xt) [] q! { ( t Eq (t) D q,t E q (xt) { xe q(t) E q (qxt) + E q(t) E q (qxt) { x(1 (1 q)t)e q(t) E q (qxt) = (xt) [] q! { xe q(qt) E q (qxt) + E q(t) E q (xt) ) + E } q(qt) E q (qxt) [] qt 1 } + xq +1 G 1 (x, qt; q) + (1 (1 q)xt)e q(t) E q (qxt) } + xq +1 G 1 (x, qt; q) } + xq +1 G 1 (x, qt; q) = xq G (x, qt; q) + G (x, t; q) + xq +1 G 1 (x, qt; q), (4.2) where we also used (3.4) in the first line, (3.12), (3.11), (3.3), and (3.16) to get the second line, (3.15) to get the fourth line, and (3.16) to get the fifth line. Corollary 4.11. The q-bernstein basis functions satisfy the recurrence relations B n+1 (x; q) = xb 1(x; n q) + q (1 xq n )B n (x; q), (4.21) B n+1 (x; q) = xq n +1 B 1 n (x; q) + (1 xqn )B n (x; q). (4.22) Proof. We use (2.9) and (3.4) to expand both sides of (4.18) in powers of t, and then we equate the coefficients of t n /[n] q! on both sides of the resulting equation to obtain (4.21). Similarly, we use (2.9) and (3.4) to expand both sides of (4.19) in powers of t, and then we equate the coefficients of t n /[n] q! on both sides of the resulting equation to obtain (4.22). Next, we establish another functional equation for the generating functions G (x, t; q) and we use this equation to derive a subdivision relation for the q-bernstein basis functions. Theorem 4.12. Let max{ q, x, y } < 1 and t < 1/ 1 q. Then the generating functions G (x, t; q) satisfy the functional equations for all l, s. In particular, G (xy, t; q) = [l] q![s] q! x l s y s t l s G l (x, t; q)g s (y, xt; q) (4.23) [] q! G (xy, t; q) = G (x, t; q)g (y, xt; q). (4.24) Proof. Formula (4.23) follows immediately from (3.16) and the identity E q (t) E q (xyt) = E q(t) E q (xt) E q (xt) E q (yxt). Formula (4.24) follows from (4.23) by letting l = and s =.

11 Corollary 4.13. The q-bernstein basis functions satisfy the subdivision relation Proof. By (2.9) and (2.3) we have Also by (2.9) we have G (x, t; q) = n B n (xy; q) = n= j= B n n j(x; q)b n j (y; q). (4.25) B n tn (x; q) [n] q! = n 1 (1 xq l ) tn [n] n= q!. l= G (xy, t; q) = n= B n tn (xy; q) [n] q! and G (y, xt; q) = n= B n (y; q)xn t n [n] q!. Substituting the last three expressions into (4.24) and equating the coefficients of t n on both sides of the resulting equation, we obtain B n n (xy; q) j 1 = [n] q! j= By (2.3), this equation reduces to (4.25). l= (1 xql ) [j] q! B n j (y; q)x n j. [n j] q! Formula (4.25) was first established in [3] using a matrix based technique, and later in [15, Proposition 5.4] using a method based on q-blossoming. Corollary 4.14. n j= j 1 Bj n (x; q) (1 tq l ) = l= n 1 (1 txq l ). (4.26) l= Proof. Formula (4.26) follows immediately from (4.25) and (2.3) by setting = and y = t in (4.25) and then replacing j by n j. A more general version of Corollary 4.14 for q-bernstein basis functions on arbitrary intervals can be derived using the same method or a method based on q-blossoming [15]. The next result is a q-version of Marsden s identity on the interval [, 1]. Corollary 4.15. (q-marsden s Identity on the Interval [, 1]). n 1 (t xq l ) = l= n ( 1) j q (j 2) Bn j (x; q)bn n j (t; 1/q) [ ]. (4.27) n j j= 1/q

12 Proof. In equation (4.26) we replace t by 1/t and multiply both sides of the resulting equation by t n, thus obtaining n 1 n j 1 (t xq l ) = Bj n (x; q)t n j (t q l ). By (2.3) we have j 1 t n j l= l= j= j 1 (t q l ) = ( 1) j q (j 2) t n j (1 tq l ) = ( 1) j q (j 2) Bn n j (t; ] 1/q). l= Combining the last two equations yields formula (4.27). This q-marsden identity is a special case of a more general q-marsden identity over arbitrary intervals established in [15, Proposition 5.1] using a method based on q-blossoming. We close this section with a very general subdivision relation. Corollary 4.16. The q-bernstein basis functions satisfy the subdivision relations l= [ n j 1/q for all l, s. B n (xy; q) = [l] q![s] q![n] q! [] q! n+l y s m=l B m l (x; q) [m] q! Bs n+l+s m (y; q)x n m [n + l + s m] q! (4.28) t Proof. If l+s, formula (4.28) is obtained by equating the coefficients of n [n] q! on both sides of (4.23). If l + s, formula (4.28) follows after multiplying (4.23) by t l+s and equating the coefficients of tn+l+s [n] q! on both sides of the resulting equation. 5 Proofs of Properties of the q-exponential Functions We are now going to prove the assertions we made without proof in Section 3 about the q- exponential function E q (x) and about the generating functions G (x, t; q). In order to avoid excessive length in some of the formulas in our proofs, here we introduce a much more concise notation. We begin with the q-shifted factorials defined by [2], [5], [8] (a; q) n = n (1 aq j 1 ), (a; q) := 1, j=1 and Notice that (a; q) = (1 aq j 1 ), q < 1. j=1 (a; q) (aq ; q) n = (a; q) n, (5.1) (a; q) (aq ; q) = (a; q). (5.2)

13 The q-bernstein basis functions can now be rewritten more concisely in terms of the q- shifted factorials. First the q-binomial coefficients defined in (2.4) can be written as [ n ] q = (q; q) n (q; q) (q; q) n = (q+1 ; q) n (q; q) n, =,..., n. (5.3) Then, by (2.3) and (5.3) the q-bernstein basis functions defined in (2.3) can be expressed as [ n B n (x; q) = ] q x (x; q) n = (q+1 ; q) n x (x; q) n, n. (5.4) (q; q) n We shall also use the identities (A; q) n = (A; q) n ( A) q (n 2) ( n 2 ) (q n+1 /A; q), (5.5) (q n+1 /A; q) n = ( 1/A) n q (n 2) (A; q)n. (5.6) These identities follow directly from the definition of the q-shifted factorials. In order to analyze the q-exponential function E q (x) which appears in the explicit formula (3.16) for the generating functions G (x, t; q), we now introduce the standard q-exponential functions e q (x) and E q (x) defined for q < 1 by [2], [5], [8] and e q (x) = n= E q (x) = x n (q; q) n, x < 1, (5.7) n= q (n 2) (q; q) n x n. (5.8) We shall see shortly (Proposition 5.2) that the three q-exponential functions E q (x), e q (x), and E q (x) are closely related. We shall use the basic hypergeometric series notation [2], [5], [8] where rφ s (a 1,..., a r ; b 1,..., b s ; q, z) = n= (a 1,..., a r ; q) n = (a 1,..., a r ; q) n (q, b 1,..., b s ; q) n z n ( ( 1) n q (n 2) ) s+1 r, (5.9) r (a j ; q) n. j=1 In order to analyze the q-exponential functions E q (x), e q (x), and E q (x), we will need the q-binomial formula [2, Theorem 1.2.1], [8, Theorem 12.2.5]. Let q < 1. Then 1φ (a; ; q, z) = n= where z < 1 or a = q n for some integer n. (a; q) n (q; q) n z n = (az; q) (z; q), (5.1) With this preparation, we are now ready to prove some identities for the q-exponential functions.

14 Lemma 5.1. (Euler s Theorem, [8, Theorem 12.2.6]) 1 e q (x) =, (x; q) (5.11) E q (x) = ( x; q). (5.12) Proof. Setting a = in (5.1) yields (5.11), and replacing z by x/a in (5.1) and letting a yields (5.12). Proposition 5.2. Proof. Equation (5.13) follows from (3.1) and (5.7). E q (x) = e q ((1 q)x), q < 1, (5.13) E 1/q (x) = E q ((1 q)x), < q < 1. (5.14) To prove equation (5.14), notice that for < q < 1, the function E 1/q (x) defined by the series in (3.1) is an entire function of x. From (3.1), (2.8), and (5.8) it follows that E 1/q (x) = n= x n [n] 1/q! = n= q (n 2) x n [n] q! = E q ((1 q)x). Next we prove equations (3.12) (3.15). Below we recall these four equations: D q,x (1/E q (x)) = 1/E q (qx), (5.15) D 1/q,x (1/E q (x)) = 1/E q (x), (5.16) E 1/q (x) = 1/E q ( x), (5.17) E q (qx) = (1 (1 q)x)e q (x). (5.18) Proof of Equations (5.15) and (5.16). By (5.13) and (5.11) we get Equation (5.15) is proved similarly. D 1/q,x (1/E q (x)) = ((1 q)x/q; q) ((1 q)x; q) (1/q 1)x = ((1 q)x; q) [(1 (1 q)x/q) 1] (1 q)x/q = 1 E q (x). Proof of Equations (5.17) and (5.18). To prove equation (5.17), observe that from equations (5.14), (5.12), (5.11), and (5.13) it follows that E 1/q (x) = E q ((1 q)x) = ( (1 q)x; q) = 1/e q ( (1 q)x) = 1/E q ( x). Equation (5.18) is an immediate consequence of (5.13) and (5.11).

15 To derive the explicit formula (3.16) for the generating functions G (x, t; q), next we introduce a new family of generating functions defined by G (x, t, c; q) = B n t n (x; q). (5.19) (c; q) n n= Theorem 5.3. Let max{ q, x, t, c } < 1. Then G (x, t, c; q) = (xt) (c; q) 2φ 1 (x, q +1 ; cq ; q, t). (5.2) Proof. This formula follows from equations (5.4), (5.1), and (5.9). Now we are ready to derive our explicit formula for the generating functions G (x, t; q) that was stated without proof in Theorem 3.1. Proof of Theorem 3.1. From (2.9), (5.19), (5.2), (5.1), (5.7), (5.11), and (5.13), we derive G (t, x; q) = G (x, (1 q)t, q; q) = = (xt) [] q! ((1 q)xt; q) = (xt) ((1 q)t; q) [] q! ((1 q)xt) 1φ (x; ; q, (1 q)t) (q; q) e q ((1 q)t) e q ((1 q)xt) = (xt) [] q! E q (t) E q (xt). 6 Two New Summation Formulas In this section we are going to derive two q-analogues of the following summation formula for the classical Bernstein polynomials defined in (2.1): B n (x) = 1 x n= (6.1) for < x < 1. Equation (6.1) can be verified by direct summation, by integration, or by induction on. To derive these new identities, we need Heine s transformation formula [8, Theorem 12.5.1]: for z < 1 and b < 1. 2φ 1 (a, b; c;q,z) = (b, az; q) (c, z; q) 2φ 1 (c/b, z; az;q,b), (6.2) Corollary 6.1. n= for < x < 1 and q < 1. B n (x; q)q n = q x ( 1 (x; q) l= x l (q; q) l ) (6.3)

16 Proof. Applying (5.19), (5.2), Heine s transformation (6.2) (with a = q +1, b = x, c =, and z = q), (5.9), (5.2), (5.1), (5.7), and (5.11), we derive which yields (6.3). B n (x; q)q n = G (x, q, ; q) = (qx) 2φ 1 (q +1, x; ; q,q) n= = (qx) (q+2, x; q) 2φ 1 (, q; q +2 ; q, x) = (qx) (x; q) (q; q) (q; q) +1 = q (x; q) x l=+1 ( x l = q (x; q) (q; q) l x 1 (x; q) l= l= x l (q +2 ; q) l ) x l (q; q) l Corollary 6.2. ( ) B n(x; q)qn = ( x) 1 (x; q) [] q! D 1/q,x x n= for < x < 1 and q < 1. (6.4) Proof. The following differentiation formula is easily verified by induction on n: D 1/q,x ((x; q) n) = ( 1) [n] q! [n ] q! (x; q) n. (6.5) By the first equality in (5.4), (2.4), and (6.5), we get = x [] q! n= n= B n (x; q)qn = x [] q! q n [n] q! [n ] q! (x; q) n n= q n ( 1) D 1/q,x ((x; q) n) = ( x) [] q! D 1/q,x ( ) (x; q) n q n. n= (6.6) Notice that in the last expression on the second line we can start the sum from n = instead of from n =, since D1/q,x annihilates polynomials of degree less than. By (5.9) we have (x; q) n q n = 2 φ 1 (q, x; ; q,q) = (x; q) 1 q 2 φ 1 (, q; q 2 ; q, x) n= = (x; q) 1 q l= x l (q 2 = (x; q) ; q) l x ( 1 1 (x; q) = (x; q) x x l+1 (q; q) l+1 l= ), (6.7) where to get the second equality we used Heine s transformation (6.2) (with a = q, b = x, c =, and z = q), and to get the last equality we used (5.7) and (5.11). Formula (6.4) now follows from (6.6) and (6.7).

17 7 Summary We have derived a suite of identities for the q-bernstein basis functions [ n B n (x; q) = ] q x n 1 j= (1 q j x), =,..., n from a collection of functional equations for their generating functions G (x, t; q) = n= B n tn (x; q). [n] q Our method in each case is first to derive a functional equation satisfied by the generating functions, and then to extract an identity for the q-bernstein basis functions by equating the coefficients of t n on both sides of this functional equation. Similar identities and proofs hold for the h-bernstein basis functions [16] using analogues of the h-exponential functions, h-derivatives, and h-integrals. Since we have derived many identities, we close here with a list of these formulas in Table 1 for easy future reference.

18 Formula q-bernstein Identity Functional Equation Basis Functions B n (x; q) = h n iq x (x; q) n G (x, t; q) = (xt) Eq(t) []q! Eq(xt) Partition of Unity Convolution Identity Monomials I Monomials II Discrete q-derivative Discrete 1/q-Derivative First Recurrence Second Recurrence Subdivision Formula I Subdivision Formula II Polynomial Identity q-marsden s Identity Summation Formula I Summation Formula II q-integration Formula 1/q-Integration Formula P n= B n P (x; q) = 1 = G (x, t; q) = E q(t) P h i n n j= j q xj B n j h (x; q) = n iq x E q(xt)g (x, t; q) = (xt) []q! Eq(t) P h i h i n=j j q Bn (x; q) = n P h j q xj =j j iq G (x, t; q) = (xt)j [j]q! Eq(t) P j n= ( 1)n q n h i 2 j n 1/q Bn (x; q) = q h i G (x, t; q)e 1/q ( t) = (xt) []q! E 1/q ( xt) 2 j 1/q xj D q,xb n (x; q) = [n] qq qb n 1 D q,xg (x, t; q) 1 (qx; q) Bn 1 (qx; q) = tq +1 G 1 (qx, t; q) tq G (qx, t; q) D 1/q,x B n (x; q) = [n] qq qb n 1 D 1/q,x G (x, t; q) 1 (x; q) Bn 1 (x; q) = tq +1 G 1 (x, t; q) tq G (x, t; q) B n+1 (x; q) D q,tg (x, t; q) = xg 1 (x, t; q) = xb 1 n (x; q) + q (1 xq n )B n (x; q) +q G (x, t; q) xg (x, qt; q) B n+1 (x; q) D q,tg (x, t; q) = xq +1 G 1 (x, qt; q) = xq n +1 B 1 n (x; q) + (1 xqn )B n (x; q) +G (x, t; q) xq G (x, qt; q) B n (xy; q) = P n j= Bn n j (x; q)bn j (y; q) G (xy, t; q) = G (x, t; q)g (y, xt; q) B n [l]q![s]q![n]q! (xy; q) = y s []q! B n+l+s j G (xy, t; q) = [l]q![s]q! x l s y s t l s []q! G l (x, t; q)g s(y, xt; q) P n+l x n j B j l (x;q) s (y;q) j=l [j]q! [n+l+s j]q! P n= B n (x; q)(t; q) = (tx; q) n G (xy, t; q) = G (x, t; q)g (y, xt; q) Q n 1 l= (t xql ) = P j B n j= ( 1) q j n 2 n j (t;1/q)b j n h i (x;q) n j P 1/q n= B n (x; q)qn = q P x 1 (x; q) l= x l «(q;q) P l n= B n (x; q)qn = ( x) 1 (x;q) []q! D 1/q,x x R 1 B n (x; q) dqx = q(1 δ n, )+1 R 1 B n (x; q) d 1/q x = [n+1]q q [n+1]q G (xy, t; q) = G (x, t; q)g (y, xt; q) G (x, t, c; q) = (xt) (c;q) 2φ 1 (x, q +1 ; cq ; q, t) G (x, t, c; q) = (xt) (c;q) 2φ 1 (x, q +1 ; cq ; q, t) D q,xb n (x; q) = [n] qq qb n 1 1 (qx; q) Bn 1 (qx; q) D 1/q,x B n (x; q) = [n] qq qb n 1 1 (x; q) Bn 1 (x; q) Table 1: Identities for the q-bernstein basis functions and the corresponding functional equations for their generating functions. Acnowlegment. The research of Yilmaz Simse was supported by the Scientific Research Project Administration of Adeniz University. References [1] M. Acigoz, S. Araci, On the generating function for Bernstein polynomials, American Institute of Physics CP1281, 21. [2] G. E. Andrews, R. A. Asey, and R. Roy, Special Functions, Cambridge University Press, Cambridge, 1999. [3] C. Disibuyu, H. Oruç, A generalization of rational Bernstein-Bézier curves, BIT 47 (27) 313-323. [4] H. Exton, q-hypergeometric Functions and Applications, New Yor: Halstead Press, Chichester: Ellis Horwood, 1983.

19 [5] G. Gasper, M. Rahman, Basic Hypergeometric Series, Cambridge University Press, Cambridge, 24. [6] R. Goldman, Pyramid Algorithms, A Dynamic Programming Approach to Curves and Surfaces for Geometric Modeling, The Morgan Kaufmann Series in Computer Graphics and Geometric Modeling, Elsevier Science, 23. [7] R. Goldman, P. Simeonov, Formulas and algorithms for quantum differentiation of quantum Bernstein bases and quantum Bézier curves based on quantum blossoming, Graphical Models, 74(6)(212) 326 334. [8] M. E. H. Ismail, Classical and Quantum Orthogonal Polynomials in one Variable, Cambridge University Press, Cambridge, 25. [9] S. Lewanowicz, P. Woźny, Generalized Bernstein polynomials, BIT, 44(24) 63 78. [1] H. Oruç, Generalized Bernstein Polynomials and Total Positivity, Ph.D. Thesis, School of Mathematical and Computational Sciences, University of St. Andrews, 1998. [11] G. M. Phillips, Bernstein polynomials based on the q-integers, The heritage of P. L. Chebyshev: a Festschrift in honor of the 7th birthday of T. J. Rivlin, Annals of Numerical Math. 1-4, 511 518, 1997. [12] G. M. Phillips, A survey of results on the q-bernstein polynomials, IMA J. Numer. Analysis 3(21) 277-288. [13] L. Ramshaw, Blossoming: A Connect-the-Dots Approach to Splines, Digital Equipment Corp., Systems Research Center, Technical Report no. 19, 1987. [14] L. Ramshaw, Blossoms are polar forms, Computer Aided Geometric Design, 6(1989) 323-358. [15] P. Simeonov, V. Zafiris, R. Goldman, q-blossoming: A new approach to algorithms and identities for q-bernstein bases and q-bézier curves, Journal of Approximation Theory 164(212), 77 14. [16] P. Simeonov, V. Zafiris, R. Goldman, h-blossoming: A new approach to algorithms and identities for h-bernstein bases and h-bézier curves, Computer Aided Geometric Design, 28(9)(211), 549-565. [17] Y. Simse, Generating functions for the Bernstein type polynomials: a new approach to deriving identities and applications for these polynomials, to appear in Hacettepe Journal of Mathematics and Statistics. [18] Y. Simse, Functional equations from generating functions: a novel approach to deriving identities for the Bernstein basis functions, Fixed Point Theory and Applications, 1:8(213), 13 pages. [19] Y. Simse, M. Acigoz, A new generating function of (q ) Bernstein-type polynomials and their interpolation function, Abstr. Appl. Anal., (21), 12 pages, Article ID 76995, DOI: 1.1155/21/76995.

Ron Goldman: rng@rice.edu Plamen Simeonov: simeonovp@uhd.edu Yilmaz Simse: ysimse@adeniz.edu.tr 2