Cascade algorithm

From HandWiki

In the mathematical topic of wavelet theory, the cascade algorithm is a numerical method for calculating function values of the basic scaling and wavelet functions of a discrete wavelet transform using an iterative algorithm. It starts from values on a coarse sequence of sampling points and produces values for successively more densely spaced sequences of sampling points. Because it applies the same operation over and over to the output of the previous application, it is known as the cascade algorithm.

Successive approximation

The iterative algorithm generates successive approximations to ψ(t) or φ(t) from {h} and {g} filter coefficients. If the algorithm converges to a fixed point, then that fixed point is the basic scaling function or wavelet.

The iterations are defined by

[math]\displaystyle{ \varphi^{(k+1)}(t)=\sum_{n=0}^{N-1} h[n] \sqrt 2 \varphi^{(k)} (2t-n) }[/math]

For the kth iteration, where an initial φ(0)(t) must be given.

The frequency domain estimates of the basic scaling function is given by

[math]\displaystyle{ \Phi^{(k+1)}(\omega)= \frac {1} {\sqrt 2} H\left( \frac {\omega} {2}\right) \Phi^{(k)}\left(\frac {\omega} {2}\right) }[/math]

and the limit can be viewed as an infinite product in the form

[math]\displaystyle{ \Phi^{(\infty)}(\omega)= \prod_{k=1}^{\infty} \frac {1} {\sqrt 2} H\left( \frac {\omega} {2^k}\right) \Phi^{(\infty)}(0). }[/math]

If such a limit exists, the spectrum of the scaling function is

[math]\displaystyle{ \Phi(\omega)= \prod_{k=1}^\infty \frac {1} {\sqrt 2} H\left( \frac {\omega} {2^k}\right) \Phi^{(\infty)}(0) }[/math]

The limit does not depends on the initial shape assume for φ(0)(t). This algorithm converges reliably to φ(t), even if it is discontinuous.

From this scaling function, the wavelet can be generated from

[math]\displaystyle{ \psi(t)= \sum_{n=- \infty}^{\infty} g[n]{\sqrt 2} \varphi^{(k)} (2t-n). }[/math]

Successive approximation can also be derived in the frequency domain.

References