Matching Pursuit Method

A general family of time-frequency atoms can be generated by scaling, translating and modulating a single window function g(t). For any scale s > 0, frequency modulation and translation u, we denote = (s, u,) and define

If g(t) is even, which is generally the case, g(t) is centered at the abscissa u. Its energy is mostly concentrated in a neighborhood of u, whose size is proportional to s. The Fourier transform of g(t) can be expressed by

Since is even, is centered at the frequency =. Its energy is concentrated in a neighborhood of , whose size is proportional to 1/s.

In this study we use software implementation of matching pursuit algorithm developed by Mallat and Zhang (1993). We use dictionary composed of Gabor functions supplemented with a canonical basis of discrete Dirac functions and the discrete Fourier basis of cosine and sine functions. The discrete Gabor time-frequency atom is defined by

where constant K_s normalizes the discrete norm of g_s to 1, and , for 0 < p < N and 0 < k < N. This dictionary provides a richer collection of atomic waveforms which are located on a much finer grid in time-frequency space than wavelet and cosine packet tables.

We compute a linear expansion of signal f over a set of atoms selected from the dictionary, in order to best match its inner structures. After m iterations, a matching pursuit decomposes a signal into

where R^mf is the residual vector after m iterations, and <f, g> denotes inner product of functions f and g. The Matching Pursuit algorithm at each step selects atom g_n(t) , for which inner product <Rⁿf, g_n> is largest.

To illustrate the decomposition into the time-frequency atoms we compute energy density defined by

where Wg_n(t,) is the Wigner distribution of atom g_n(t,). Unlike the Wigner and the Cohen class distributions, it does not include cross terms. Energy distribution of atoms from our dictionary are visible as horizontal lines for cosine functions, vertical lines for Dirac , or ellipses with axes proportional to time and frequency spread of Gabor functions (Figure 1).