The Wavelet Digest Homepage
Return to the homepage
Search the complete Wavelet Digest database
Help about the Wavelet Digest mailing list
About the Wavelet Digest
The Digest The Community
 Latest Issue  Back Issues  Events  Gallery
The Wavelet Digest
   -> Volume 6, Issue 4


Answer: Smoothest Scaling Function (WD 6.3 #21)
 
images/spacer.gifimages/spacer.gif Reply into Digest
Previous :: Next  
Author Message
unser@helix.nih.gov (Michael Unser)
Guest





PostPosted: Thu Mar 06, 1997 4:02 pm    
Subject: Answer: Smoothest Scaling Function (WD 6.3 #21)
Reply with quote

#18 Answer: Smoothest Scaling Function (WD 6.3 #21)

>What is the smoothest scaling function, w(x), solving a
>'wavelet equation' of type:
>(1) w(x) = c(0)w(2x) + c(1)w(2x-1) + c(2)w(2x-2) +...+ c(m)w(2x-m),
>where the coefficients c(i) satisfy the 'usual conditions':
>(2) sum c(i) over odd i = 1,
>(3) sum c(i) over even i = 1,
>and where w is normalized to have, say, integral = 1.

First, the smoothness of a scaling function cannot be greater than p
(the order of the function = its ability to reproduce polynomials of
degree p-1 = corresponds to the factor (1+z)^p in the refinement
filter C(z)=(1+z)^p*Q(z)).

Second, it is easy to see that the B-spline of order p (or degree p-1)
is the shortest scaling function of order p, because Q(z)=1.

Third, B-splines of order p have a Sobolev smoothness index of p-1/2
(cf. Strang's book). However, B-splines are not the smoothest scaling
functions of order p! For instance, one can take Q(z)=(1+EPSILON+z),
which for EPSILON>0 sufficiently small (but non-zero !) can achieve
the maximum smoothness s_max=p. But this also shows that the B-splines
are the smoothest scaling for a refinement filter of a given length:
the example above with EPSILON=0 is a B-spline of order p+1 which has
a smoothness s_max=p+1/2 --- this is 1/2 better than anything else of
the same length !

Michael Unser

BEIP, Building 13/3N17
National Institutes of Health
Bethesda, MD 20892-5766, USA

Email : Unser@helix.nih.gov
All times are GMT + 1 Hour
Page 1 of 1

 
Jump to: 
 


disclaimer - webmaster@wavelet.org
Powered by phpBB

This page was created in 0.025918 seconds : 18 queries executed : GZIP compression disabled