Shannon wavelet

From HandWiki

In functional analysis, the Shannon wavelet (or sinc wavelets) is a decomposition that is defined by signal analysis by ideal bandpass filters. Shannon wavelet may be either of real or complex type.

Shannon wavelet is not well-localized (noncompact) in the time domain,but its Fourier transform is band-limited (compact support). Hence Shannon wavelet has poor time localization but has good frequency localization. These characteristics are in stark contrast to those of the Haar wavelet. The Haar and sinc systems are Fourier duals of each other.

Definition

Sinc function is the starting point for the definition of the Shannon wavelet.

Scaling function

First, we define the scaling function to be the sinc function.

[math]\displaystyle{ \phi^{\text{(Sha)}}(t) := \frac {\sin \pi t} {\pi t} = \operatorname{sinc}(t). }[/math]

And define the dilated and translated instances to be

[math]\displaystyle{ \phi^n_k(t) := 2^{n/2}\phi^{\text{(Sha)}}(2^n t-k) }[/math]

where the parameter [math]\displaystyle{ n,k }[/math] means the dilation and the translation for the wavelet respectively.

Then we can derive the Fourier transform of the scaling function:

[math]\displaystyle{ \Phi^{\text{(Sha)}}(\omega) = \frac{1}{2\pi}\Pi(\frac{\omega}{2\pi}) = \begin{cases} \frac{1}{2\pi}, & \mbox{if } {|\omega| \le \pi}, \\ 0 & \mbox{if } \mbox{otherwise}. \\ \end{cases} }[/math] where the (normalised) gate function is defined by

[math]\displaystyle{ \Pi ( x):= \begin{cases} 1, & \mbox{if } {|x| \le 1/2}, \\ 0 & \mbox{if } \mbox{otherwise}. \\ \end{cases} }[/math] Also for the dilated and translated instances of scaling function: [math]\displaystyle{ \Phi^n_k(\omega) = \frac{2^{-n/2}}{2\pi}e^{-i\omega(k+1)/2^n}\Pi(\frac{\omega}{2^{n+1}\pi}) }[/math]

Mother wavelet

Use [math]\displaystyle{ \Phi^{\text{(Sha)}} }[/math] and multiresolution approximation we can derive the Fourier transform of the Mother wavelet:

[math]\displaystyle{ \Psi^{\text{(Sha)}}(\omega) = \frac{1}{2\pi}e^{-i\omega} \bigg(\Pi(\frac{\omega}{\pi}-\frac{3}{2})+\Pi(\frac{\omega}{\pi}+\frac{3}{2})\bigg) }[/math]

And the dilated and translated instances:

[math]\displaystyle{ \Psi^n_k(\omega) = \frac{2^{-n/2}}{2\pi}e^{-i\omega(k+1)/2^n} \bigg(\Pi(\frac{\omega}{2^n\pi}-\frac{3}{2})+\Pi(\frac{\omega}{2^n\pi}+\frac{3}{2})\bigg) }[/math]

Then the shannon mother wavelet function and the family of dilated and translated instances can be obtained by the inverse Fourier transform:

[math]\displaystyle{ \psi^{\text{(Sha)}}(t) = \frac{\sin\pi(t-(1/2))-\sin2\pi(t-(1/2))}{\pi(t-1/2)} =\operatorname{sinc}\bigg(t-\frac{1}{2}\bigg)-2\operatorname{sinc}\bigg(2(t-\frac{1}{2})\bigg) }[/math]

[math]\displaystyle{ \psi^n_k(t) = 2^{n/2}\psi^{\text{(Sha)}}(2^nt-k) }[/math]

Property of mother wavelet and scaling function

  • Mother wavelets are orthonormal, namely,

[math]\displaystyle{ \lt \psi^n_k(t), \psi^m_h(t)\gt =\delta^{nm}\delta_{hk}= \begin{cases} 1, & \text{if }h=k \text{ and } n=m\\ 0, & \text{otherwise} \end{cases} }[/math]

  • The translated instances of scaling function at level [math]\displaystyle{ n=0 }[/math] are orthogonal

[math]\displaystyle{ \lt \phi^0_k(t), \phi^0_h(t)\gt =\delta^{kh} }[/math]

  • The translated instances of scaling function at level [math]\displaystyle{ n=0 }[/math] are orthogonal to the mother wavelets

[math]\displaystyle{ \lt \phi^0_k(t), \psi^m_h(t)\gt =0 }[/math]

  • Shannon wavelets has an infinite number of vanishing moments.

Reconstruction of a Function by Shannon Wavelets

Suppose [math]\displaystyle{ f(x)\in L_2(\mathbb{R}) }[/math] such that [math]\displaystyle{ \operatorname{supp}\operatorname{FT}\{f\}\subset[-\pi,\pi] }[/math] and for any dilation and the translation parameter [math]\displaystyle{ n,k }[/math],

[math]\displaystyle{ \Bigg|\int^\infty_{-\infty}f(t)\phi^0_k(t)dt\Bigg|\lt \infty }[/math], [math]\displaystyle{ \Bigg|\int^\infty_{-\infty}f(t)\psi^n_k(t)dt\Bigg|\lt \infty }[/math]

Then

[math]\displaystyle{ f(t)=\sum^\infty_{k=\infty}\alpha_k\phi^0_k(t) }[/math] is uniformly convergent, where [math]\displaystyle{ \alpha_k=f(k) }[/math]

Real Shannon wavelet

Real Shannon wavelet

The Fourier transform of the Shannon mother wavelet is given by:

[math]\displaystyle{ \Psi^{(\operatorname{Sha}) }(w) = \prod \left( \frac {w- 3 \pi /2} {\pi}\right)+\prod \left( \frac {w+ 3 \pi /2} {\pi}\right). }[/math]

where the (normalised) gate function is defined by

[math]\displaystyle{ \prod ( x):= \begin{cases} 1, & \mbox{if } {|x| \le 1/2}, \\ 0 & \mbox{if } \mbox{otherwise}. \\ \end{cases} }[/math]

The analytical expression of the real Shannon wavelet can be found by taking the inverse Fourier transform:

[math]\displaystyle{ \psi^{(\operatorname{Sha}) }(t) = \operatorname{sinc} \left( \frac {t} {2}\right)\cdot \cos \left( \frac {3 \pi t} {2}\right) }[/math]

or alternatively as

[math]\displaystyle{ \psi^{(\operatorname{Sha})}(t)=2 \cdot \operatorname{sinc}(2t)-\operatorname{sinc}(t), }[/math]

where

[math]\displaystyle{ \operatorname{sinc}(t):= \frac {\sin {\pi t}} {\pi t} }[/math]

is the usual sinc function that appears in Shannon sampling theorem.

This wavelet belongs to the [math]\displaystyle{ C^\infty }[/math]-class of differentiability, but it decreases slowly at infinity and has no bounded support, since band-limited signals cannot be time-limited.

The scaling function for the Shannon MRA (or Sinc-MRA) is given by the sample function:

[math]\displaystyle{ \phi^{(Sha)}(t)= \frac {\sin \pi t} {\pi t} = \operatorname{sinc}(t). }[/math]

Complex Shannon wavelet

In the case of complex continuous wavelet, the Shannon wavelet is defined by

[math]\displaystyle{ \psi^{(CSha) }(t)=\operatorname{sinc}(t) \cdot e^{-2\pi i t} }[/math],

References

  • S.G. Mallat, A Wavelet Tour of Signal Processing, Academic Press, 1999, ISBN 0-12-466606-X
  • C.S. Burrus, R.A. Gopinath, H. Guo, Introduction to Wavelets and Wavelet Transforms: A Primer, Prentice-Hall, 1988, ISBN 0-13-489600-9.