Uniform Approximation and Bernstein Polynomials with Coefficients in the Unit Interval
|
|
- Roderick Harper
- 6 years ago
- Views:
Transcription
1 Unifor Approxiation and Bernstein Polynoials with Coefficients in the Unit Interval Weiang Qian and Marc D. Riedel Electrical and Coputer Engineering, University of Minnesota 200 Union St. S.E. Minneapolis, MN 55455, USA {qianx030, Ivo Rosenberg Matheatics and Statistics University of Montreal C.P. 6128, Succ. Centre-Ville Montreal, Quebec, Canada H3C 3J7 Abstract This paper presents two ain results. The first result pertains to unifor approxiation with Bernstein polynoials. We show that, given a power-for polynoial g, we can obtain a Bernstein polynoial of degree with coefficients that are as close as desired to the corresponding values of g evaluated at the points 0, 1,..., 1, provided that is sufficiently large. The second result pertains to a subset of Bernstein polynoials: those with coefficients that are all in the unit interval. We show that polynoials in this subset ap the open interval 0, 1 into the open interval 0, 1 and ap the points 0 and 1 into the closed interval [0, 1]. The otivation for this wor is our research on probabilistic coputation with digital circuits. Our design ethodology, called stochastic logic, is based on Bernstein polynoials with coefficients that correspond to probability values; accordingly, the coefficients ust be values in the unit interval. The atheatics presented here provide a necessary and sufficient test for deciding whether polynoial operations can be ipleented with stochastic logic. 1 Introduction The Weierstrass approxiation theore is a faous theore in atheatical analysis. It asserts that every continuous function defined on a closed interval can be uniforly approxiated as closely as desired by a polynoial function [1]. The Weierstrass Approxiation Theore: Let f be a continuous function defined on the closed interval [a, b]. For any ɛ > 0, there exists a polynoial function p such that for all x in [a, b], we have fx px < ɛ. The theore can be proved by a transforation with Bernstein polynoials [2]. By a linear substitution, the interval [a, b] can be transfored into the unit interval [0, 1]. Thus, the original stateent of the theore holds if and only if the theore holds for every continuous function f defined on the interval [0, 1]. This wor is supported by a grant fro the Seiconductor Research Corporation s Focus Center Research Progra on Functional Engineered Nano-Architectonics, contract No NT
2 A Bernstein polynoial of degree n is a polynoial expressed in the following for [3]: β,n b,n x, 1 where each β,n, 0, 1,..., n, is a real nuber and n b,n x x 1 x n. 2 The coefficients β,n are called Bernstein coefficients and the polynoials b 0,n x, b 1,n x,..., b n,n x are called Bernstein basis polynoials of degree n. Define the n-th Bernstein polynoial for f to be B n f; x f b,n x. n In 1912, Bernstein showed the following result [4, 5]: The Bernstein Theore: Let f be a continuous function defined on the closed interval [0, 1]. For any ɛ > 0, there exists a positive integer M such that for all x in [0, 1] and integer M, we have fx B f; x < ɛ. Note that the function B f; x is a polynoial on x. Thus, based on the Bernstein Theore, the Weierstrass Approxiation Theore holds. Given a power-for polynoial g of degree n, it is well nown that for any n, g can be uniquely converted into a Bernstein polynoial of degree [6]. Cobining this fact with the Bernstein Theore, we have the following corollary. Corollary 1 Let g be a polynoial of degree n. For any ɛ > 0, there exists a positive integer M n such that for all x in [0, 1] and integer M, we have β, g b, x < ɛ, where β 0,, β 1,,..., β, satisfy gx β, b, x. In the first part of the paper, we prove a stronger result than this: Theore 1 Let g be a polynoial of degree n 0. For any ɛ > 0, there exists a positive integer M n such that for all integers M and 0, 1,...,, we have β, g < ɛ, where β 0,, β 1,,..., β, satisfy gx β, b, x. 2
3 Cobining Theore 1 with the fact that b, x 1, we can easily prove Corollary 1. In the second part of the paper, we consider a subset of Bernstein polynoials: those with coefficients that are all in the unit interval [0, 1]. Definition 1 Define U to be the set of Bernstein polynoials with coefficients that are all in the unit interval [0, 1]: { } U px n 1, 0 β 0,n, β 1,n,..., β n,n 1, such that px β,n b,n x. The question we as is: which polynoials can be converted into Bernstein polynoials in U? Definition 2 Define the set V to be the set of polynoials which are either identically equal to 0 or equal to 1, or ap the open interval 0, 1 into 0, 1 and the points 0 and 1 into the closed interval [0, 1], i.e., V {px px 0, or px 1, or 0 < px < 1, x 0, 1 and 0 p0, p1 1}. We prove that the two sets are equivalent: Theore 2 V U. In what follows, we will refer to a Bernstein polynoial of degree n converted fro a polynoial g as the Bernstein polynoial of degree n of g. When we say that a polynoial is of degree n, we ean that the power-for of the polynoial is of degree n. Exaple 1 Consider the polynoial gx 3x 8x 2 + 6x 3. It aps the open interval 0, 1 into 0, 1 with g0 0, g1 1. Thus, g is in the set V. Based on Theore 2, we have that g is in the set U. We verify this by considering Bernstein polynoials of increasing degree. The Bernstein polynoial of degree 3 of g is Note that here the coefficient β 2,3 2 3 < 0. The Bernstein polynoial of degree 4 of g is Note that here the coefficient β 3,4 1 4 < 0. gx b 1,3 x 2 3 b 2,3x + b 3,3 x. gx 3 4 b 1,4x b 2,4x 1 4 b 3,4x + b 4,4 x. The Bernstein polynoial of degree 5 of gx is Note that here all the coefficients are in [0, 1]. gx 3 5 b 1,5x b 2,5x + b 5,5 x. 3
4 Since the Bernstein polynoial of degree 5 of g satisfies Definition 1, we conclude that g is in the set U. Exaple 2 Consider the polynoial gx 1 4 x + x2. Since g0.5 0, thus g is not in the set V. Based on Theore 2, we have that g is not in the set U. We verify this. By contraposition, suppose that there exist n 1 and 0 β 0,n, β 1,n,..., β n,n 1 such that Since g0.5 0, therefore, gx β,n b,n x. β,n b,n Note that for all 0, 1,..., n, b,n 0.5 > 0. Thus, we have that for all 0, 1,..., n, β,n 0. Therefore, gx 0, which contradicts the original assuption about g. Thus, g is not in the set U. The reainder of the paper is organized as follows. In Section 2, we present soe atheatical preliinaries pertaining to Bernstein polynoials. In Section 3, we prove Theore 1. Based on this theore, in Section 4, we prove Theore 2. Finally, we conclude the paper with a discussion on applications of these theores to our research in probabilistic coputation with digital circuits. 2 Properties of Bernstein Polynoials We list soe pertinent properties of Bernstein polynoials. a The positivity property: For all 0, 1,..., n and all x in [0, 1], we have b The partition of unity property: b,n x 0. 3 The binoial expansion of the left-hand side of the equality x + 1 x n 1 shows that the su of all Bernstein basis polynoials of degree n is the constant 1, i.e., b,n x 1. 4 c Converting power-for coefficients to Bernstein coefficients: The set of Bernstein basis polynoials b 0,n x, b 1,n x,..., b n,n x fors a basis of the vector space of polynoials of real coefficients and degree no ore than n [6]. Each power basis function x j can be uniquely expressed as a linear cobination of the n + 1 Bernstein basis polynoials: x j σ j b,n x, 5 for j 0, 1,..., n. To deterine the eleents of the transforation atrix σ, we write x j x j x + 1 x nj 4
5 and perfor a binoial expansion on the right hand side. This gives j x j n,n x, jb for j 0, 1,..., n. Therefore, we have j σ j { σj n 1, j j for j 0, for j >. 6 Suppose that a power-for polynoial of degree no ore than n is gx a,n x 7 and the Bernstein polynoial of degree n of g is gx β,n b,n x. 8 Substituting Equations 5 and 6 into Equation 7 and coparing the Bernstein coefficients, we have n 1 β,n a j,n σ j a j,n. 9 j j Equation 9 provide a eans for obtaining Bernstein coefficients fro power-for coefficients. d Degree elevation: Based on Equation 2, we have that for all 0, 1,...,, 1 +1 b,+1 x x 1 x 1 b, x, b +1,+1 x x 1 x +1 + x +1 1 x or b, x +1 b,+1 x + +1 b +1,+1 x b,+1x b +1,+1x. 10 Given a power-for polynoial g of degree n, for any n, g can be uniquely converted into a Bernstein polynoial of degree. Suppose that the Bernstein polynoials of degree and +1 degree + 1 of g are β, b, x and β,+1 b,+1 x, respectively. We have β, b, x 5 +1 β,+1 b,+1 x. 11
6 Substituting Equation 10 into the left-hand side of Equation 11 and coparing the Bernstein coefficients, we have β 0,, for 0 β,+1 +1 β 1, β,, for 1 12 β,, for + 1. Equation 12 provides a eans for obtaining the coefficients of the Bernstein polynoial of degree + 1 of g fro the coefficients of the Bernstein polynoial of degree of g. We will call this procedure degree elevation. For convenience, given a Bernstein polynoial gx n β,nb,n x, we can also express it as gx c,n x 1 x n, 13 where c,n n β,n, 14 for 0, 1,..., n. Substituting Equation 14 into Equation 12, we have 3 A Proof of Theore 1 c 0,, for 0 c,+1 c 1, + c,, for 1 c,, for + 1. Suppose that the polynoial g is of degree n. Applying Equation 15 recursively, we can express c, as a linear cobination of c 0,n, c 1,n,..., c n,n. Lea 1 Let g be a polynoial of degree n. For any n, suppose that the Bernstein polynoial of degree of g is gx c, x 1 x. Let c, 0 for all < 0 and all >. Then for all 0, 1,...,, we have c, n n Proof: We prove the lea by induction on n. i 15 c +n+i,n. 16 Base case: For n 0, the right-hand side of Equation 16 reduces to 0 0 c,n c,, so the equation holds. Inductive step: Suppose that Equation 16 holds for soe n and all 0, 1,...,. Consider + 1. Since we assue that c 1, c +1, 0, Equation 15 can be written as c,+1 c 1, + c,, 17 6
7 for all 0,..., + 1. With our convention that c i,n 0 for all i < 0 and i > n, it is easily seen that c 1, 0 n n i c 1+n+i,n, c +1, 0 n n i c +1+n+i,n. Together with the induction hypothesis, we conclude that for all 1, 0,...,, + 1 c, n n Based on Equations 17 and 18, for all 0, 1,..., + 1, we have c,+1 n n i c 1+n+i,n + i c +n+i,n. 18 n n In the first su, we change the suation index to j i 1. We obtain c,+1 n1 j1 n 0 Applying the basic forula c,+1 c +n1,n + n c +n+j,n + j + 1 c +n1,n + n1 r q n1 n r 1 + q n j + 1 n j [ n j + 1 r 1 q + j c +n+j,n n j, we obtain c +n+j,n + c,n c +n+j,n. ] c +n+j,n + +1n + 1 n Thus Equation 16 holds for + 1. By induction, it holds for all. Rear: Equation 16 can be forulated as c, in{,n} iax{0,+n} i n c,n. n c 1+n+i,n. n c i,n, 19 i for all n and 0, 1,...,. Indeed, in Equation 16, first use the basic forula and then change the suation index to j + n + i to obtain c, n n c +n+i,n n i j+n Note that c j,n 0 iplies 0 j n. This yields Equation 19. n c j,n. j r r q r q 7
8 Lea 2 Let n be a positive integer. For all integer, and i such that we have Proof: For siplicity, we define δ n i > n, 0, ax{0, + n} i in{, n}, 20 i 1 ni n2. 21 n i i 1 ni n!!! i! n + i!! ni1 j j j n i. Now 1 i n + i n + 1 i1 i1 i j 1 ni1 j We obtain an upper bound for n i i1 n i 1 i. i j by replacing j in Equation 22 with its least value, 0: 1 ni1 1 i i i ni. i We need the following siple inequality: for real nubers 0 x y 1 and a non-negative integer l, l1 y l x l y x y j x l1j y xl. 23 Fro Equation 20, we obtain 0 i in{, n} and so we can use Equation 23 for 0 x i y 1, l n i 0. We obtain i δ 1 ni n i i i ni ni i i i n i Since 0 1, 0 1, and 0 i n, we obtain i ni i in i i. ni i 22 i in i i 8 in i > n2.
9 Therefore, δ Siilarly, we obtain a lower bound for i 1 ni n i n i the first product and with n i in the second product, obtaining n i1 i 1 ni1 1 i i1 j i j i i n + i ni i i n i > n2. 24 by replacing the index j in Equation 22 with i in 1 i i n + i n + i ni1 ni. 1 i n Thus, proceeding as above, we have i δ 1 ni n i i ni i i n + i ni i n + i [ i ] i i [ ni ni ] n + i ni i i +. i n + i i Due to Equation 20, we have and so we obtain 0 i i 1, 0 n + i n + i 1, δ i i i + i ni n + i ni. 25 n + i Applying Equation 23 twice to the right-hand side of Equation 25, we obtain δ i i + n i i n + i n + i i2 i Fro Equation 20, we have n i2 + n + i. Therefore, δ i 1 ni 0 i 1, 0 n + i 1. n i i2 + n i 2 ni + nn i n2. 26 Equations 24 and 26 together yield Equation 21. Now we give a proof of Theore 1. 9
10 Theore 1 Let g be a polynoial of degree n 0. For any ɛ > 0, there exists a positive integer M n such that for all integer M and 0, 1,...,, we have β, g < ɛ, where β 0,, β 1,,..., β, satisfy that gx β, b, x. Proof: For n 0, g is a constant polynoial. Suppose that gx y, where y is a constant value. We select M 1. Then, for all integers M and all integers 0, 1,...,, we have β, y g. Thus, the theore holds. { } n 2 For n > 0, we select M such that M > ax c i,n, 2n, where the real nubers ɛ c 0,n, c 1,n,..., c n,n satisfy gx c i,n x i 1 x ni. 27 Now consider any M. Since 2n ax { n 2 ɛ } c i,n, 2n < M, we have n > n. Consider the following three cases for. 1. The case where n n. Here ax{0, + n} 0 and in{, n} n. Thus, the suation indices in Equation 19 range fro 0 to n. Therefore, β, c n, i c i,n. 28 Substituting x with in Equation 27, we have g i c i,n 1 ni. 29 By Lea 2, since 0 < n < and 0, Equation 21 holds for all 0 ax{0, + n} i in{, n} n. Thus, by Equations 21, 28, 29 and the well-nown inequality x i x i, we have [ β n i i, g 1 ] ni c i,n n i i 1 ni c i,n n2 c i,n. 10
11 Since n2 ɛ c i,n < M, we have n 2 c i,n < ɛ. 30 Therefore, for all n n, we have β, g < ɛ. 2. The case where 0 < n. Since > 2n, we have + n < n < 0. Thus, ax{0, + n} 0 and in{, n}. Thus, the suation indices in Equation 19 range fro 0 to. Therefore, β, c n, i c i,n. 31 When + 1 i n, we have that i and so i 1 ni i1 1 ni < n n2. 32 By Lea 2, since 0 < n < and 0, Equation 21 holds for all 0 ax{0, + n} i in{, n}. Thus, by Equations 21, 29, 30, 31, 32 and the inequality x i x i, we have β n i i, g c i,n 1 ni c i,n [ n i i 1 ] ni i c i,n 1 ni c i,n i+1 n i i 1 ni c i i,n + 1 ni c i,n n2 c i,n < ɛ. i+1 3. The case where n <. Since > 2n, we have n < n <. Thus, ax{0, + n} + n and in{, n} n. Now, the suation indices in Equation 19 range fro + n to n. Therefore, β, c n, i c i,n. 33 i+n When 0 i + n 1, we have that n i. Thus, i 1 ni 1 i 1 ni1 < n n
12 By Lea 2, since 0 < n < and 0, Equation 21 holds for all + n ax{0, + n} i in{, n} n. Thus, by Equations 21, 29, 30, 33, 34 and the inequality x i x i, we have β n i i, g c i,n 1 ni c i,n i+n [ n ] ni i i i+n n i i+n n2 c i,n < ɛ. i 1 c i,n i 1 ni c i,n + In conclusion, if M, then for all 0, 1,...,, we have β, g < ɛ. 4 A Proof of Theore 2 +n1 +n1 1 ni c i,n i 1 ni c i,n We deonstrate that the sets U and V defined in the introduction see Definitions 1 and 2 are one and the sae. We deonstrate that U V and V U separately. First, we prove the forer the easier one. Then we use Theore 1 to prove the latter. Theore 3 U V. Proof: Let n 1 and β,n 0, for all 0 n. Then the polynoial px Let n 1 and β,n 1, for all 0 n. Then, by Equation 4, the polynoial px Thus 0 U and 1 U. Fro the definition of V, 0 V and 1 V. Now consider any polynoial p U such that p 0 and p 1. 0 β 0,n, β 1,n,..., β n,n 1 such that px β,n b,n x. β,n b,n x 0. β,n b,n x 1. There exist n 1 and Fro Equations 3, 4 and the fact that 0 β 0,n, β 1,n,..., β n,n 1, for all x in [0, 1], we have 0 px β,n b,n x b,n x 1. We further clai that for all x in 0, 1, we ust have 0 < px < 1. By contraposition, we assue that there exists a 0 < x 0 < 1, such that px 0 0 or px 0 1. Since for 0 < x 0 < 1, we have 0 px 0 1, thus px 0 0 or 1. 12
13 We first consider the case that px 0 0. Since 0 < x 0 < 1, it is not hard to see that for all 0, 1,..., n, b,n x 0 > 0. Thus, px 0 0 iplies that for all 0, 1,..., n, β,n 0. In this case, for any real nuber x, px n β,nb,n x 0, which contradicts the assuption that px 0. Siilarly, in the case that px 0 1, we can show that px 1, which contradicts the assuption that px 1. In both cases, we get a contradiction; this proves the clai that for all x in 0, 1, 0 < px < 1. Therefore, for any polynoial p U such that p 0 and p 1, we have p V. Since we showed at the outset that 0 U, 1 U, 0 V and 1 V, thus, for any polynoial p U, we have p V. Therefore, U V. Next we prove the clai that V U. We will first show that each of four possible different categories of polynoials in the set V are in the set U. The different categories are tacled in Theores 4 and 5 and Corollaries 2 and 3. Theore 4 Let g be a polynoial of degree n apping the open interval 0, 1 into 0, 1 with 0 g0, g1 < 1. Then g U. Proof: Since g is continuous on the closed interval [0, 1], it attains its axiu value M g on [0, 1]. Since gx < 1, for all x [0, 1], we have M g < 1. Let ɛ 1 1 M g > 0. By Theore 1, there exists a positive integer M 1 n such that for all integers M 1 and 0, 1,...,, we have β, g < ɛ 1, where β 0,, β 1,,..., β, satisfy that gx β, b, x. Thus, for all M 1 and all 0, 1,...,, β, < g + ɛ 1 M g + 1 M g Denote by r the ultiplicity of 0 as a root of gx where r 0 if g0 0 and by s the ultiplicity of 1 as a root of gx where s 0 if g1 0. We can factorize gx as gx x r 1 x s hx, 36 where hx is a polynoial, satisfying that h0 0 and h1 0. We show that h0 > 0. By the way of contraposition, suppose that h0 0. Since h0 0, we have h0 < 0. By the continuity of the polynoial hx, there exists soe 0 < x < 1, such that hx < 0. Thus, gx x r 1 x s hx < 0. However, gx > 0, for all x 0, 1. Therefore, h0 > 0. Siilarly, we have h1 > 0. gx Since gx > 0 for all x in 0, 1, we have hx x r > 0 for all x in 0, 1. In view of 1 x s the fact that h0 > 0 and h1 > 0, we have hx > 0, for all x in [0, 1]. Since hx is continuous on the closed interval [0, 1], it attains its iniu value h on [0, 1]. Clearly, h > 0. Let ɛ 2 h > 0. By Theore 1, there exists a positive integer M 2 n r s, such that for all integers d M 2 and 0, 1,..., d, we have γ,d h < ɛ 2, where γ 0,d, γ 1,d,..., γ d,d satisfy d that d hx γ,d b,d x
14 Thus, for all d M 2 and all 0, 1,..., d, γ,d > h ɛ 2 h h 0. d Cobining Equations 36 and37, we have gx x r 1 x s hx x r 1 x s d d+r+s γ d,d d+r+s +r d + r + s + r β,d+r+s b,d+r+s x, d γ,d b,d x x r 1 x s x +r 1 x d+s d+r r d d γ,d x 1 x d γ d r,d r b,d+r+s x d+r+s where β,d+r+s are the coefficients of the Bernstein polynoial of degree d + r + s of g and 0, for 0 < r and d + r < d + r + s β,d+r+s γ r,d r d > 0, for r d + r. d+r+s Thus, when d + r + s M 2 + r + s, we have for all 0, 1,...,, β, According to Equations 35 and 38, if we select an 0 ax{m 1, M 2 + r + s}, then gx can be expressed as a Bernstein polynoial of degree 0 : 0 gx β,0 b,0 x, with 0 β,0 1, for all 0, 1,..., 0. Therefore, g U. Theore 5 Let g be a polynoial of degree n apping the open interval 0, 1 into 0, 1 with g0 0 and g1 1. Then g U. Proof: Denote by r the ultiplicity of 0 as a root of gx. We can factorize gx as gx x r hx, 39 where hx is a polynoial satisfying h0 0. By a siilar reasoning as in the proof of Theore 4, we obtain h0 > 0. Since for all x in 0, 1], hx gx x r > 0, we have for all x in [0, 1], hx > 0. Since hx is continuous on the closed interval [0, 1], it attains its iniu value h on [0, 1]. Clearly, h > 0. Let ɛ 1 h > 0. By Theore 1, there exists a positive integer M 1 n r such that for all integers d M 1 and 0, 1,..., d, we have γ,d h < ɛ 1, where γ 0,d, γ 1,d,..., γ d,d satisfy d hx d γ,d b,d x
15 Thus, for all d M 1 and all 0, 1,..., d, γ,d > h ɛ 1 h h 0. d Cobining Equations 39 and 40, we have gx x r hx x r d γ d,d d+r +r d d d γ,d b,d x x r γ,d x 1 x d d+r x +r 1 x d γ d r,d d+r r b,d+r x β,d+r b,d+r x, d + r + r where β,d+r are the coefficients of the Bernstein polynoial of degree d + r of g and 0, for 0 < r β,d+r γ r,d r d > 0, for r d + r. d+r r d+r Thus, when d + r M 1 + r, we have for all 0, 1,...,, Let β, g x 1 gx. 42 Then g aps the open interval 0, 1 into 0, 1 with g 0 1, g 1 0. ultiplicity of 1 as a root of g x. Thus, we can factorize g x as Denote by s the g x 1 x s h x, 43 where h x is a polynoial satisfying that h 1 0. As in the proof of Theore 4, we obtain h 1 > 0. Since for all x in [0, 1, h x g x 1 x s > 0, we have for all x [0, 1], h x > 0. Since h x is continuous on the closed interval [0, 1], it attains its iniu value h on [0, 1]. Clearly, h > 0. Let ɛ 2 h > 0. By Theore 1, there exists a positive integer M 2 n s such that for all integers q M 2 and 0, 1,..., q, we have γ,q h < ɛ 2, where γ0,q q, γ 1,q,..., γ q,q satisfy h x Thus, for all q M 2 and all 0, 1,..., q, γ,q > h q q γ,q b,qx. 44 Cobining Equations 42, 43 and 44, we have ɛ 2 h h 0. gx 1 g x 1 1 x s h x 1 1 x s 1 1 x s q γ,q q x 1 x q 1 15 q γ,q b,qx q γ q,q q+s q + s x 1 x q+s.
16 Further using 4, we obtain q+s gx b,q+s x q γ,q q q+s q+s b,q+s x β,q+s b,q+s x, where the β,q+s s are the coefficients of the Bernstein polynoial of degree q + s of g: 1 γ,q q β,q+s q+s < 1, for 0 q 1, for q < q + s. Thus, when q + s M 2 + s, we have for all 0, 1,...,, β, According to Equations 41 and 45, if we select an 0 ax{m 1 + r, M 2 + s}, then gx can be expressed as a Bernstein polynoial of degree 0 : 0 gx β,0 b,0 x, with 0 β,0 1, for all 0, 1,..., 0, Therefore, g U. Lea 3 If a polynoial p is in the set U, then the polynoial 1 p is also in the set U. Proof: Since p is in the set U, there exist n 1 and 0 β 0,n, β 1,n,..., β n,n 1 such that By Equation 4, we have 1 px b,n x px β,n b,n x. β,n b,n x 1 β,n b,n x γ,n b,n x, where γ,n 1β,n satisfying 0 γ,n 1, for all 0, 1,..., n. Therefore, 1p is in the set U. Corollary 2 Let g be a polynoial of degree n apping the open interval 0, 1 into 0, 1 with 0 < g0, g1 1. Then g U. Proof: Let polynoial h 1 g. Then h aps 0, 1 into 0, 1 with 0 h0, h1 < 1. By Theore 4, h U. By Lea 3, g 1 h is also in the set U. Corollary 3 Let g be a polynoial of degree n apping the open interval 0, 1 into 0, 1 with g0 1 and g1 0. Then g U. 16
17 Proof: Let the polynoial h 1 g. Then h aps 0, 1 into 0, 1 with h0 0, h1 1. By Theore 5, h U. By Lea 3, g 1 h is also in the set U. Cobining Theore 4, Theore 5, Corollary 2 and Corollary 3, we show that V U. Theore 6 V U. Proof: Based on the definition of V, for any polynoial p V, we have one of following five cases. 1. The case where p 0 or p 1. In the proof of Theore 3, we have shown that 0 U and 1 U. Thus p U. 2. The case where p aps the open interval 0, 1 into 0, 1 with 0 p0, p1 < 1. By Theore 4, p U. 3. The case where p aps the open interval 0, 1 into 0, 1 with 0 < p0, p1 1. By Corollary 2, p U. 4. The case where p aps the open interval 0, 1 into 0, 1 with p0 0 and p1 1. By Theore 5, p U. 5. The case where p aps the open interval 0, 1 into 0, 1 with p0 1 and p1 0. By Corollary 3, p U. In suary, for any polynoial p V, we have p U. Thus, V U. Based on Theores 3 and 6, we have proved Theore 2: V U. 5 Discussion We are interested in Bernstein polynoials with coefficients in the unit interval because this concept has applications in the area of digital circuit design. Specifically, the concept is a atheatical prerequisite for a design ethodology that we have been advocating called stochastic logic [7 9]. We provide a brief overview of this application and point the reader to further sources. Stochastic logic ipleents Boolean function with inputs that are rando Boolean variables. A Boolean function f on n variables x 1, x 2,..., x n is a apping f : {0, 1} n {0, 1}. With stochastic logic, the variables x 1, x 2,..., x n are a set of independent rando Boolean variables, i.e., for 1 i n, x i has a certain probability p i 0 p i 1 of being one and a probability 1 p i of being zero. With rando Boolean variables as inputs, the output is also a rando Boolean variable: the function f has a certain probability p o of being one and a probability 1 p o of being zero. If ipleented by digital circuitry, stochastic logic can be viewed as coputation that transfors input probabilities into output probabilities [8]. Given an arbitrary Boolean function f and 17
18 a set of input probabilities p 1, p 2,..., p n that correspond to the probabilities of the input rando Boolean variables being one, the output probability p o is a function on p 1, p 2,..., p n. In fact, we have shown that the general for of the function is a ultivariate polynoial on variables p 1,..., p n with integer coefficients and with the degree of each variable no ore than one [9]. Exaple 3 Consider stochastic logic based on the Boolean function fx 1, x 2, x 3 x 1 x 2 x 1 x 3, where eans logical AND conjunction, eans logical OR disjunction, and eans logical negation. The Boolean function f evaluates to one if and only if the 3-tuple x 1, x 2, x 3 taes values fro the set {0, 0, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1}. The probability of the output being one is p o Prf 1 Pr x 1, x 2, x 3 : x 1, x 2, x 3 {0, 0, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1} Prx 1 0, x 2 0, x Prx 1 0, x 2 1, x Prx 1 1, x 2 1, x Prx 1 1, x 2 1, x 3 1. If x 1, x 2, and x 3 are independent rando Boolean variables with probability p 1, p 2, and p 3 of being one, respectively, then we obtain p o 1 p 1 1 p 2 p p 1 p 2 p 3 + p 1 p 2 1 p 3 + p 1 p 2 p 3 1 p 1 p 3 + p 1 p 2 p 1 p 2 + p 3 p 1 p 3, 46 which confirs that the function coputed by stochastic logic is a ultivariate polynoial on arguents p 1, p 2, and p 3 with integer coefficients and with the degree of each variable no ore than 1. In design probles, we encounter univariate polynoials that have real coefficients and degree greater than 1. Soeties it is possible to ipleent these by setting soe of the probabilities p i to be a coon variable x and the others to be constants. For exaple, if we set p 1 p 3 x and p in Equation 46, then we obtain the polynoial gx 1.75x x 2. With different underlying Boolean functions and different assignents of probability values, we can ipleent any different univariate polynoials. An interesting and yet practical question is: which univariate polynoials can be ipleented by stochastic logic? Define the set W to be the set of univariate polynoials that can be ipleented. We are interested in characterizing the set W. In [9] we showed that U W, i.e., if a polynoial can be expressed as a Bernstein polynoial with all coefficients in the unit interval, then the polynoial can be ipleented by stochastic logic. In this paper, we proved that V U. Thus, we have V W. Further, in [9] we showed that W V, i.e., if a polynoial can be ipleented by stochastic logic, then it is either identically equal to 0 or equal to 1, or it aps the open interval 0, 1 into the open interval 0, 1 and aps the points 0 and 1 into the closed interval [0, 1]. Therefore, we conclude that W V, i.e., a polynoial can be ipleented by stochastic logic if and only if it is either identically equal to 0 or equal to 1, or it aps the open interval 0, 1 into the open interval 0, 1 and aps the points 0 and 1 into the closed interval [0, 1]. This necessary and sufficient conditions allows us to answer the question of whether any given polynoial can be ipleented by stochastic logic. Based on the atheatics, we have proposed a constructive design ethod [9]. An overview of the ethod and its applications in circuit design will appear in a forthcoing Research Highlights article in Counications of the ACM [10]. 18
19 References [1] P. J. Davis, Interpolation and Approxiation. Blaisdell Publishing Copany, [2] G. Lorentz, Bernstein Polynoials. University of Toronto Press, [3] J. Berchtold and A. Bowyer, Robust arithetic for ultivariate Bernstein-for polynoials, Coputer-Aided Design, vol. 32, no. 11, pp , [4] S. N. Bernstein, Déonstration du théorèe de Weierstrass fondée sur le calcul des probabilités, Counications of the Kharov Matheatical Society, vol. 13, pp. 1 2, [5] K. M. Levasseur, A probabilistic proof of the Weierstrass approxiation theore, The Aerican Matheatical Monthly, vol. 91, no. 4, pp , [6] R. Faroui and V. Rajan, On the nuerical condition of polynoials in Bernstein for, Coputer Aided Geoetric Design, vol. 4, no. 3, pp , [7] B. Gaines, Stochastic coputing systes, in Advances in Inforation Systes Science. Plenu, 1969, vol. 2, ch. 2, pp [8] B. Brown and H. Card, Stochastic neural coputation I: Coputational eleents, IEEE Transactions on Coputers, vol. 50, no. 9, pp , [9] W. Qian and M. D. Riedel, The synthesis of robust polynoial arithetic with stochastic logic, in Design Autoation Conference, 2008, pp [10], The synthesis of stochastic logic, Counications of the ACM to appear,
The Weierstrass Approximation Theorem
36 The Weierstrass Approxiation Theore Recall that the fundaental idea underlying the construction of the real nubers is approxiation by the sipler rational nubers. Firstly, nubers are often deterined
More informationCurious Bounds for Floor Function Sums
1 47 6 11 Journal of Integer Sequences, Vol. 1 (018), Article 18.1.8 Curious Bounds for Floor Function Sus Thotsaporn Thanatipanonda and Elaine Wong 1 Science Division Mahidol University International
More informationOn Conditions for Linearity of Optimal Estimation
On Conditions for Linearity of Optial Estiation Erah Akyol, Kuar Viswanatha and Kenneth Rose {eakyol, kuar, rose}@ece.ucsb.edu Departent of Electrical and Coputer Engineering University of California at
More informationNew upper bound for the B-spline basis condition number II. K. Scherer. Institut fur Angewandte Mathematik, Universitat Bonn, Bonn, Germany.
New upper bound for the B-spline basis condition nuber II. A proof of de Boor's 2 -conjecture K. Scherer Institut fur Angewandte Matheati, Universitat Bonn, 535 Bonn, Gerany and A. Yu. Shadrin Coputing
More informationPolygonal Designs: Existence and Construction
Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G
More informationA Bernstein-Markov Theorem for Normed Spaces
A Bernstein-Markov Theore for Nored Spaces Lawrence A. Harris Departent of Matheatics, University of Kentucky Lexington, Kentucky 40506-0027 Abstract Let X and Y be real nored linear spaces and let φ :
More information4 = (0.02) 3 13, = 0.25 because = 25. Simi-
Theore. Let b and be integers greater than. If = (. a a 2 a i ) b,then for any t N, in base (b + t), the fraction has the digital representation = (. a a 2 a i ) b+t, where a i = a i + tk i with k i =
More informationExplicit solution of the polynomial least-squares approximation problem on Chebyshev extrema nodes
Explicit solution of the polynoial least-squares approxiation proble on Chebyshev extrea nodes Alfredo Eisinberg, Giuseppe Fedele Dipartiento di Elettronica Inforatica e Sisteistica, Università degli Studi
More informationPerturbation on Polynomials
Perturbation on Polynoials Isaila Diouf 1, Babacar Diakhaté 1 & Abdoul O Watt 2 1 Départeent Maths-Infos, Université Cheikh Anta Diop, Dakar, Senegal Journal of Matheatics Research; Vol 5, No 3; 2013 ISSN
More informationCharacterization of the Line Complexity of Cellular Automata Generated by Polynomial Transition Rules. Bertrand Stone
Characterization of the Line Coplexity of Cellular Autoata Generated by Polynoial Transition Rules Bertrand Stone Abstract Cellular autoata are discrete dynaical systes which consist of changing patterns
More informationDivisibility of Polynomials over Finite Fields and Combinatorial Applications
Designs, Codes and Cryptography anuscript No. (will be inserted by the editor) Divisibility of Polynoials over Finite Fields and Cobinatorial Applications Daniel Panario Olga Sosnovski Brett Stevens Qiang
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I Contents 1. Preliinaries 2. The ain result 3. The Rieann integral 4. The integral of a nonnegative
More informationM ath. Res. Lett. 15 (2008), no. 2, c International Press 2008 SUM-PRODUCT ESTIMATES VIA DIRECTED EXPANDERS. Van H. Vu. 1.
M ath. Res. Lett. 15 (2008), no. 2, 375 388 c International Press 2008 SUM-PRODUCT ESTIMATES VIA DIRECTED EXPANDERS Van H. Vu Abstract. Let F q be a finite field of order q and P be a polynoial in F q[x
More informationADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE
ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE CHRISTOPHER J. HILLAR Abstract. A long-standing conjecture asserts that the polynoial p(t = Tr(A + tb ] has nonnegative coefficients whenever is
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationa a a a a a a m a b a b
Algebra / Trig Final Exa Study Guide (Fall Seester) Moncada/Dunphy Inforation About the Final Exa The final exa is cuulative, covering Appendix A (A.1-A.5) and Chapter 1. All probles will be ultiple choice
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationReceived: September 24, Research Article. pubs.acs.org/synthbio
wdb ACSJCA JCA..65/W Unicode research.3f (R3.6.i:32 2. alpha 39) 25/7/5 :3: PROD-JCAVA rq_755967 9/3/26 7:52:52 JCA-DEFAULT pubs.acs.org/synthbio Cheical Reaction Networks for Coputing Polynoials 2 Sayed
More informationOn the Communication Complexity of Lipschitzian Optimization for the Coordinated Model of Computation
journal of coplexity 6, 459473 (2000) doi:0.006jco.2000.0544, available online at http:www.idealibrary.co on On the Counication Coplexity of Lipschitzian Optiization for the Coordinated Model of Coputation
More information1 Proof of learning bounds
COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #4 Scribe: Akshay Mittal February 13, 2013 1 Proof of learning bounds For intuition of the following theore, suppose there exists a
More informationSampling How Big a Sample?
C. G. G. Aitken, 1 Ph.D. Sapling How Big a Saple? REFERENCE: Aitken CGG. Sapling how big a saple? J Forensic Sci 1999;44(4):750 760. ABSTRACT: It is thought that, in a consignent of discrete units, a certain
More informationBipartite subgraphs and the smallest eigenvalue
Bipartite subgraphs and the sallest eigenvalue Noga Alon Benny Sudaov Abstract Two results dealing with the relation between the sallest eigenvalue of a graph and its bipartite subgraphs are obtained.
More information#A52 INTEGERS 10 (2010), COMBINATORIAL INTERPRETATIONS OF BINOMIAL COEFFICIENT ANALOGUES RELATED TO LUCAS SEQUENCES
#A5 INTEGERS 10 (010), 697-703 COMBINATORIAL INTERPRETATIONS OF BINOMIAL COEFFICIENT ANALOGUES RELATED TO LUCAS SEQUENCES Bruce E Sagan 1 Departent of Matheatics, Michigan State University, East Lansing,
More informationMax-Product Shepard Approximation Operators
Max-Product Shepard Approxiation Operators Barnabás Bede 1, Hajie Nobuhara 2, János Fodor 3, Kaoru Hirota 2 1 Departent of Mechanical and Syste Engineering, Bánki Donát Faculty of Mechanical Engineering,
More informationList Scheduling and LPT Oliver Braun (09/05/2017)
List Scheduling and LPT Oliver Braun (09/05/207) We investigate the classical scheduling proble P ax where a set of n independent jobs has to be processed on 2 parallel and identical processors (achines)
More informationPattern Recognition and Machine Learning. Learning and Evaluation for Pattern Recognition
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lesson 1 4 October 2017 Outline Learning and Evaluation for Pattern Recognition Notation...2 1. The Pattern Recognition
More informationSolutions 1. Introduction to Coding Theory - Spring 2010 Solutions 1. Exercise 1.1. See Examples 1.2 and 1.11 in the course notes.
Solutions 1 Exercise 1.1. See Exaples 1.2 and 1.11 in the course notes. Exercise 1.2. Observe that the Haing distance of two vectors is the iniu nuber of bit flips required to transfor one into the other.
More informationIn this chapter, we consider several graph-theoretic and probabilistic models
THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions
More informationTight Information-Theoretic Lower Bounds for Welfare Maximization in Combinatorial Auctions
Tight Inforation-Theoretic Lower Bounds for Welfare Maxiization in Cobinatorial Auctions Vahab Mirrokni Jan Vondrák Theory Group, Microsoft Dept of Matheatics Research Princeton University Redond, WA 9805
More informationarxiv: v1 [cs.ds] 3 Feb 2014
arxiv:40.043v [cs.ds] 3 Feb 04 A Bound on the Expected Optiality of Rando Feasible Solutions to Cobinatorial Optiization Probles Evan A. Sultani The Johns Hopins University APL evan@sultani.co http://www.sultani.co/
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationE0 370 Statistical Learning Theory Lecture 5 (Aug 25, 2011)
E0 370 Statistical Learning Theory Lecture 5 Aug 5, 0 Covering Nubers, Pseudo-Diension, and Fat-Shattering Diension Lecturer: Shivani Agarwal Scribe: Shivani Agarwal Introduction So far we have seen how
More informationON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN
More informationComputational and Statistical Learning Theory
Coputational and Statistical Learning Theory Proble sets 5 and 6 Due: Noveber th Please send your solutions to learning-subissions@ttic.edu Notations/Definitions Recall the definition of saple based Radeacher
More informationarxiv: v2 [math.co] 3 Dec 2008
arxiv:0805.2814v2 [ath.co] 3 Dec 2008 Connectivity of the Unifor Rando Intersection Graph Sion R. Blacburn and Stefanie Gere Departent of Matheatics Royal Holloway, University of London Egha, Surrey TW20
More informationComputable Shell Decomposition Bounds
Coputable Shell Decoposition Bounds John Langford TTI-Chicago jcl@cs.cu.edu David McAllester TTI-Chicago dac@autoreason.co Editor: Leslie Pack Kaelbling and David Cohn Abstract Haussler, Kearns, Seung
More informationA symbolic operator approach to several summation formulas for power series II
A sybolic operator approach to several suation forulas for power series II T. X. He, L. C. Hsu 2, and P. J.-S. Shiue 3 Departent of Matheatics and Coputer Science Illinois Wesleyan University Blooington,
More informationAn EGZ generalization for 5 colors
An EGZ generalization for 5 colors David Grynkiewicz and Andrew Schultz July 6, 00 Abstract Let g zs(, k) (g zs(, k + 1)) be the inial integer such that any coloring of the integers fro U U k 1,..., g
More informationOn the summations involving Wigner rotation matrix elements
Journal of Matheatical Cheistry 24 (1998 123 132 123 On the suations involving Wigner rotation atrix eleents Shan-Tao Lai a, Pancracio Palting b, Ying-Nan Chiu b and Harris J. Silverstone c a Vitreous
More information13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices
CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay
More informationTHE AVERAGE NORM OF POLYNOMIALS OF FIXED HEIGHT
THE AVERAGE NORM OF POLYNOMIALS OF FIXED HEIGHT PETER BORWEIN AND KWOK-KWONG STEPHEN CHOI Abstract. Let n be any integer and ( n ) X F n : a i z i : a i, ± i be the set of all polynoials of height and
More informationA1. Find all ordered pairs (a, b) of positive integers for which 1 a + 1 b = 3
A. Find all ordered pairs a, b) of positive integers for which a + b = 3 08. Answer. The six ordered pairs are 009, 08), 08, 009), 009 337, 674) = 35043, 674), 009 346, 673) = 3584, 673), 674, 009 337)
More informationarxiv: v1 [math.nt] 14 Sep 2014
ROTATION REMAINDERS P. JAMESON GRABER, WASHINGTON AND LEE UNIVERSITY 08 arxiv:1409.411v1 [ath.nt] 14 Sep 014 Abstract. We study properties of an array of nubers, called the triangle, in which each row
More information1. INTRODUCTION AND RESULTS
SOME IDENTITIES INVOLVING THE FIBONACCI NUMBERS AND LUCAS NUMBERS Wenpeng Zhang Research Center for Basic Science, Xi an Jiaotong University Xi an Shaanxi, People s Republic of China (Subitted August 00
More informationThe Transactional Nature of Quantum Information
The Transactional Nature of Quantu Inforation Subhash Kak Departent of Coputer Science Oklahoa State University Stillwater, OK 7478 ABSTRACT Inforation, in its counications sense, is a transactional property.
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationThe degree of a typical vertex in generalized random intersection graph models
Discrete Matheatics 306 006 15 165 www.elsevier.co/locate/disc The degree of a typical vertex in generalized rando intersection graph odels Jerzy Jaworski a, Michał Karoński a, Dudley Stark b a Departent
More informationarxiv: v3 [quant-ph] 18 Oct 2017
Self-guaranteed easureent-based quantu coputation Masahito Hayashi 1,, and Michal Hajdušek, 1 Graduate School of Matheatics, Nagoya University, Furocho, Chikusa-ku, Nagoya 464-860, Japan Centre for Quantu
More informationA Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)
1 A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine (1900 words) Contact: Jerry Farlow Dept of Matheatics Univeristy of Maine Orono, ME 04469 Tel (07) 866-3540 Eail: farlow@ath.uaine.edu
More informationComputable Shell Decomposition Bounds
Journal of Machine Learning Research 5 (2004) 529-547 Subitted 1/03; Revised 8/03; Published 5/04 Coputable Shell Decoposition Bounds John Langford David McAllester Toyota Technology Institute at Chicago
More informationTesting Properties of Collections of Distributions
Testing Properties of Collections of Distributions Reut Levi Dana Ron Ronitt Rubinfeld April 9, 0 Abstract We propose a fraework for studying property testing of collections of distributions, where the
More informationResearch Article Rapidly-Converging Series Representations of a Mutual-Information Integral
International Scholarly Research Network ISRN Counications and Networking Volue 11, Article ID 5465, 6 pages doi:1.54/11/5465 Research Article Rapidly-Converging Series Representations of a Mutual-Inforation
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationLost-Sales Problems with Stochastic Lead Times: Convexity Results for Base-Stock Policies
OPERATIONS RESEARCH Vol. 52, No. 5, Septeber October 2004, pp. 795 803 issn 0030-364X eissn 1526-5463 04 5205 0795 infors doi 10.1287/opre.1040.0130 2004 INFORMS TECHNICAL NOTE Lost-Sales Probles with
More informationAbout the definition of parameters and regimes of active two-port networks with variable loads on the basis of projective geometry
About the definition of paraeters and regies of active two-port networks with variable loads on the basis of projective geoetry PENN ALEXANDR nstitute of Electronic Engineering and Nanotechnologies "D
More informationStudy on Markov Alternative Renewal Reward. Process for VLSI Cell Partitioning
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 40, 1949-1960 HIKARI Ltd, www.-hikari.co http://dx.doi.org/10.12988/ia.2013.36142 Study on Markov Alternative Renewal Reward Process for VLSI Cell Partitioning
More informationPrerequisites. We recall: Theorem 2 A subset of a countably innite set is countable.
Prerequisites 1 Set Theory We recall the basic facts about countable and uncountable sets, union and intersection of sets and iages and preiages of functions. 1.1 Countable and uncountable sets We can
More informationOn a Multisection Style Binomial Summation Identity for Fibonacci Numbers
Int J Contep Math Sciences, Vol 9, 04, no 4, 75-86 HIKARI Ltd, www-hiarico http://dxdoiorg/0988/ics0447 On a Multisection Style Binoial Suation Identity for Fibonacci Nubers Bernhard A Moser Software Copetence
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationSupplement to: Subsampling Methods for Persistent Homology
Suppleent to: Subsapling Methods for Persistent Hoology A. Technical results In this section, we present soe technical results that will be used to prove the ain theores. First, we expand the notation
More informationEstimating Parameters for a Gaussian pdf
Pattern Recognition and achine Learning Jaes L. Crowley ENSIAG 3 IS First Seester 00/0 Lesson 5 7 Noveber 00 Contents Estiating Paraeters for a Gaussian pdf Notation... The Pattern Recognition Proble...3
More informationResearch Article Approximate Multidegree Reduction of λ-bézier Curves
Matheatical Probles in Engineering Volue 6 Article ID 87 pages http://dxdoiorg//6/87 Research Article Approxiate Multidegree Reduction of λ-bézier Curves Gang Hu Huanxin Cao and Suxia Zhang Departent of
More informationSymbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm
Acta Polytechnica Hungarica Vol., No., 04 Sybolic Analysis as Universal Tool for Deriving Properties of Non-linear Algoriths Case study of EM Algorith Vladiir Mladenović, Miroslav Lutovac, Dana Porrat
More informationOn Constant Power Water-filling
On Constant Power Water-filling Wei Yu and John M. Cioffi Electrical Engineering Departent Stanford University, Stanford, CA94305, U.S.A. eails: {weiyu,cioffi}@stanford.edu Abstract This paper derives
More information3.8 Three Types of Convergence
3.8 Three Types of Convergence 3.8 Three Types of Convergence 93 Suppose that we are given a sequence functions {f k } k N on a set X and another function f on X. What does it ean for f k to converge to
More informationLectures 8 & 9: The Z-transform.
Lectures 8 & 9: The Z-transfor. 1. Definitions. The Z-transfor is defined as a function series (a series in which each ter is a function of one or ore variables: Z[] where is a C valued function f : N
More informationSolutions of some selected problems of Homework 4
Solutions of soe selected probles of Hoework 4 Sangchul Lee May 7, 2018 Proble 1 Let there be light A professor has two light bulbs in his garage. When both are burned out, they are replaced, and the next
More informationRandomized Recovery for Boolean Compressed Sensing
Randoized Recovery for Boolean Copressed Sensing Mitra Fatei and Martin Vetterli Laboratory of Audiovisual Counication École Polytechnique Fédéral de Lausanne (EPFL) Eail: {itra.fatei, artin.vetterli}@epfl.ch
More informationLinear recurrences and asymptotic behavior of exponential sums of symmetric boolean functions
Linear recurrences and asyptotic behavior of exponential sus of syetric boolean functions Francis N. Castro Departent of Matheatics University of Puerto Rico, San Juan, PR 00931 francis.castro@upr.edu
More informationNew Classes of Positive Semi-Definite Hankel Tensors
Miniax Theory and its Applications Volue 017, No., 1 xxx New Classes of Positive Sei-Definite Hankel Tensors Qun Wang Dept. of Applied Matheatics, The Hong Kong Polytechnic University, Hung Ho, Kowloon,
More informationDescent polynomials. Mohamed Omar Department of Mathematics, Harvey Mudd College, 301 Platt Boulevard, Claremont, CA , USA,
Descent polynoials arxiv:1710.11033v2 [ath.co] 13 Nov 2017 Alexander Diaz-Lopez Departent of Matheatics and Statistics, Villanova University, 800 Lancaster Avenue, Villanova, PA 19085, USA, alexander.diaz-lopez@villanova.edu
More informationPhysics 215 Winter The Density Matrix
Physics 215 Winter 2018 The Density Matrix The quantu space of states is a Hilbert space H. Any state vector ψ H is a pure state. Since any linear cobination of eleents of H are also an eleent of H, it
More informationESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics
ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents
More informationReed-Muller Codes. m r inductive definition. Later, we shall explain how to construct Reed-Muller codes using the Kronecker product.
Coding Theory Massoud Malek Reed-Muller Codes An iportant class of linear block codes rich in algebraic and geoetric structure is the class of Reed-Muller codes, which includes the Extended Haing code.
More informationMath 262A Lecture Notes - Nechiporuk s Theorem
Math 6A Lecture Notes - Nechiporuk s Theore Lecturer: Sa Buss Scribe: Stefan Schneider October, 013 Nechiporuk [1] gives a ethod to derive lower bounds on forula size over the full binary basis B The lower
More informationBernoulli Wavelet Based Numerical Method for Solving Fredholm Integral Equations of the Second Kind
ISSN 746-7659, England, UK Journal of Inforation and Coputing Science Vol., No., 6, pp.-9 Bernoulli Wavelet Based Nuerical Method for Solving Fredhol Integral Equations of the Second Kind S. C. Shiralashetti*,
More informationSTOPPING SIMULATED PATHS EARLY
Proceedings of the 2 Winter Siulation Conference B.A.Peters,J.S.Sith,D.J.Medeiros,andM.W.Rohrer,eds. STOPPING SIMULATED PATHS EARLY Paul Glasseran Graduate School of Business Colubia University New Yor,
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationINDEPENDENT SETS IN HYPERGRAPHS
INDEPENDENT SETS IN HYPERGRAPHS Abstract. Many iportant theores and conjectures in cobinatorics, such as the theore of Szeerédi on arithetic progressions and the Erdős-Stone Theore in extreal graph theory,
More informationThe Euler-Maclaurin Formula and Sums of Powers
DRAFT VOL 79, NO 1, FEBRUARY 26 1 The Euler-Maclaurin Forula and Sus of Powers Michael Z Spivey University of Puget Sound Tacoa, WA 98416 spivey@upsedu Matheaticians have long been intrigued by the su
More informationLeast Squares Fitting of Data
Least Squares Fitting of Data David Eberly, Geoetric Tools, Redond WA 98052 https://www.geoetrictools.co/ This work is licensed under the Creative Coons Attribution 4.0 International License. To view a
More informationTesting equality of variances for multiple univariate normal populations
University of Wollongong Research Online Centre for Statistical & Survey Methodology Working Paper Series Faculty of Engineering and Inforation Sciences 0 esting equality of variances for ultiple univariate
More informationAsynchronous Gossip Algorithms for Stochastic Optimization
Asynchronous Gossip Algoriths for Stochastic Optiization S. Sundhar Ra ECE Dept. University of Illinois Urbana, IL 680 ssrini@illinois.edu A. Nedić IESE Dept. University of Illinois Urbana, IL 680 angelia@illinois.edu
More informationSaddle Points in Random Matrices: Analysis of Knuth Search Algorithms
Saddle Points in Rando Matrices: Analysis of Knuth Search Algoriths Micha Hofri Philippe Jacquet Dept. of Coputer Science INRIA, Doaine de Voluceau - Rice University Rocquencourt - B.P. 05 Houston TX 77005
More informationResearch Article Some Formulae of Products of the Apostol-Bernoulli and Apostol-Euler Polynomials
Discrete Dynaics in Nature and Society Volue 202, Article ID 927953, pages doi:055/202/927953 Research Article Soe Forulae of Products of the Apostol-Bernoulli and Apostol-Euler Polynoials Yuan He and
More informationON SEQUENCES OF NUMBERS IN GENERALIZED ARITHMETIC AND GEOMETRIC PROGRESSIONS
Palestine Journal of Matheatics Vol 4) 05), 70 76 Palestine Polytechnic University-PPU 05 ON SEQUENCES OF NUMBERS IN GENERALIZED ARITHMETIC AND GEOMETRIC PROGRESSIONS Julius Fergy T Rabago Counicated by
More informationA new type of lower bound for the largest eigenvalue of a symmetric matrix
Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box
More informationAnalysis of Polynomial & Rational Functions ( summary )
Analysis of Polynoial & Rational Functions ( suary ) The standard for of a polynoial function is ( ) where each of the nubers are called the coefficients. The polynoial of is said to have degree n, where
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationOn Poset Merging. 1 Introduction. Peter Chen Guoli Ding Steve Seiden. Keywords: Merging, Partial Order, Lower Bounds. AMS Classification: 68W40
On Poset Merging Peter Chen Guoli Ding Steve Seiden Abstract We consider the follow poset erging proble: Let X and Y be two subsets of a partially ordered set S. Given coplete inforation about the ordering
More informationCSE525: Randomized Algorithms and Probabilistic Analysis May 16, Lecture 13
CSE55: Randoied Algoriths and obabilistic Analysis May 6, Lecture Lecturer: Anna Karlin Scribe: Noah Siegel, Jonathan Shi Rando walks and Markov chains This lecture discusses Markov chains, which capture
More informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationThe Fundamental Basis Theorem of Geometry from an algebraic point of view
Journal of Physics: Conference Series PAPER OPEN ACCESS The Fundaental Basis Theore of Geoetry fro an algebraic point of view To cite this article: U Bekbaev 2017 J Phys: Conf Ser 819 012013 View the article
More informationBest Arm Identification: A Unified Approach to Fixed Budget and Fixed Confidence
Best Ar Identification: A Unified Approach to Fixed Budget and Fixed Confidence Victor Gabillon Mohaad Ghavazadeh Alessandro Lazaric INRIA Lille - Nord Europe, Tea SequeL {victor.gabillon,ohaad.ghavazadeh,alessandro.lazaric}@inria.fr
More informationA Low-Complexity Congestion Control and Scheduling Algorithm for Multihop Wireless Networks with Order-Optimal Per-Flow Delay
A Low-Coplexity Congestion Control and Scheduling Algorith for Multihop Wireless Networks with Order-Optial Per-Flow Delay Po-Kai Huang, Xiaojun Lin, and Chih-Chun Wang School of Electrical and Coputer
More informationA Study on B-Spline Wavelets and Wavelet Packets
Applied Matheatics 4 5 3-3 Published Online Noveber 4 in SciRes. http://www.scirp.org/ournal/a http://dx.doi.org/.436/a.4.5987 A Study on B-Spline Wavelets and Wavelet Pacets Sana Khan Mohaad Kaliuddin
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationSupplementary Material for Fast and Provable Algorithms for Spectrally Sparse Signal Reconstruction via Low-Rank Hankel Matrix Completion
Suppleentary Material for Fast and Provable Algoriths for Spectrally Sparse Signal Reconstruction via Low-Ran Hanel Matrix Copletion Jian-Feng Cai Tianing Wang Ke Wei March 1, 017 Abstract We establish
More informationOn Certain C-Test Words for Free Groups
Journal of Algebra 247, 509 540 2002 doi:10.1006 jabr.2001.9001, available online at http: www.idealibrary.co on On Certain C-Test Words for Free Groups Donghi Lee Departent of Matheatics, Uni ersity of
More information