Remote Access Proceedings of the American Mathematical Society
Green Open Access

Proceedings of the American Mathematical Society

ISSN 1088-6826(online) ISSN 0002-9939(print)



Asymptotic analysis of Daubechies polynomials

Authors: Jianhong Shen and Gilbert Strang
Journal: Proc. Amer. Math. Soc. 124 (1996), 3819-3833
MSC (1991): Primary 41A58
MathSciNet review: 1346987
Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: To study wavelets and filter banks of high order, we begin with the zeros of $ {\mathbf {B}}_{p}(y)$. This is the binomial series for $(1-y)^{-p}$, truncated after $p$ terms. Its zeros give the $p-1$ zeros of the Daubechies filter inside the unit circle, by $z+z^{-1} = 2-4y$. The filter has $p$ additional zeros at $z = -1$, and this construction makes it orthogonal and maximally flat. The dilation equation leads to orthogonal wavelets with $p$ vanishing moments. Symmetric biorthogonal wavelets (generally better in image compression) come similarly from a subset of the zeros of $ {\mathbf {B}}_{p}(y)$. We study the asymptotic behavior of these zeros. Matlab shows a remarkable plot for $p = 70$. The zeros approach a limiting curve $|4y(1-y)| = 1$ in the complex plane, which is the circle $|z-z^{-1}| = 2$. All zeros have $|y| \le 1/2$, and the rightmost zeros approach $y = 1/2$ (corresponding to $z= \pm i$ ) with speed $p^{- 1/2}$. The curve $|4y(1-y)| = {(4 \pi p)}^{{1}/{2p}} \, |1-2y|^{ 1/p}$ gives a very accurate approximation for finite $p$. The wide dynamic range in the coefficients of $ {\mathbf {B}}_{p}(y)$ makes the zeros difficult to compute for large $p$. Rescaling $y$ by $4$ allows us to reach $p = 80$ by standard codes.

References [Enhancements On Off] (What's this?)

  • 1. N. Anderson, E.B. Saff and R.S. Varga, On the Eneström-Kakeya Theorem and its sharpness, Linear Algebra Appl. 28 (1979), 5-16. MR 81i:26011
  • 2. I. Daubechies, Ten Lectures on Wavelets, SIAM, Philadelphia, 1992. MR 92e:42025
  • 3. H.E. Fettis, J.C. Caslin and K.R. Cramer, Complex zeros of the error function and of the complementary error function, Math. Comp. 27 (1973), 401-407. MR 48:5333
  • 4. M. Lang and B.C. Frenzel, Polynomial Root Finding, preprint, Rice University, 1994.
  • 5. M. Marden, Geometry of Polynomials, Mathematical Surveys No. 3, AMS, Providence, 1966, p. 137. MR 37:1562
  • 6. Y. Meyer, Wavelets: Algorithms and Applications, SIAM, Philadelphia, 1993. MR 95f:94005
  • 7. G. Strang and T. Nguyen, Wavelets and Filter Banks, Wellesley-Cambridge Press, Wellesley, 1996.
  • 8. G. Szegö, Über eine eigenschaft der exponentialreihe, Sitzungsber. Berlin Math. Ges. 23 (1924), 50-64.
  • 9. R.S. Varga, Scientific Computation on Mathematical Problems and Conjectures, SIAM, Philadelphia, 1990. MR 92b:65012

Similar Articles

Retrieve articles in Proceedings of the American Mathematical Society with MSC (1991): 41A58

Retrieve articles in all journals with MSC (1991): 41A58

Additional Information

Jianhong Shen
Affiliation: Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139

Gilbert Strang
Affiliation: Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139

Received by editor(s): June 25, 1995
Dedicated: Dedicated to Gabor Szegö on the 100th anniversary of his birth
Communicated by: James Glimm
Article copyright: © Copyright 1996 American Mathematical Society

American Mathematical Society