Adaptive Gabor representation
Adaptive Gabor representation (AGR) is a Gabor representation of a signal where its variance is adjustable. There's always a trade-off between time resolution and frequency resolution in traditional short-time Fourier transform (STFT). A long window leads to high frequency resolution and low time resolution. On the other hand, high time resolution requires shorter window, with the expense of low frequency resolution. By choosing the proper elementary function for signal with different spectrum structure, adaptive Gabor representation is able to accommodate both narrowband and wideband signal.
Gabor expansion
In 1946, Dennis Gabor suggested that a signal can be represented in two dimensions, with time and frequency coordinates. And the signal can be expanded into a discrete set of Gaussian elementary signals.
Definition
The Gabor expansion of signal s(t) is defined by this formula:
- [math]\displaystyle{ s(t)=\sum_{m=-\infty}^\infty \sum_{n=-\infty}^\infty C_{m,n}h(t-mT)e^{jnt\Omega} }[/math]
where h(t) is the Gaussian elementary function:
- [math]\displaystyle{ h(t)=\left( \frac{\alpha}{\pi} \right)^\frac{1}{4}e^{\left( -\frac{\alpha}{2}t^2 \right)} }[/math]
Once the Gabor elementary function is determined, the Gabor coefficients [math]\displaystyle{ C_{m,n} }[/math]can be obtained by the inner product of s(t) and a dual function [math]\displaystyle{ \gamma(t) }[/math]
- [math]\displaystyle{ C_{m,n}=\int s(t)\gamma^*(t-mT)e^{-jnt\Omega} \, dt. }[/math]
[math]\displaystyle{ T }[/math] and [math]\displaystyle{ \Omega }[/math] denote the sampling steps of time and frequency and satisfy the criteria
- [math]\displaystyle{ T\Omega\leqq2\pi \, }[/math]
Relationship between Gabor representation and Gabor transform
Gabor transform simply computes the Gabor coefficients [math]\displaystyle{ C_{m,n} }[/math] for the signal s(t).
Adaptive expansion
Adaptive signal expansion is defined as
- [math]\displaystyle{ s\left( t \right)=\sum_p B_p h_p(t) }[/math]
where the coefficients [math]\displaystyle{ B_p }[/math] are obtained by the inner product of the signal s(t) and the elementary function [math]\displaystyle{ h_p }[/math]
- [math]\displaystyle{ B_p =\left \langle s,h_p\right \rangle \, }[/math]
Coeffients [math]\displaystyle{ B_p }[/math] represent the similarity between the signal and elementary function.
Adaptive signal decomposition is an iterative operation, aim to find a set of elementary function [math]\displaystyle{ \left\{ h_p(t) \right\} }[/math], which is most similar to the signal's time-frequency structure.
First, start with w=0 and [math]\displaystyle{ s_0\left( t \right)=s\left( t \right) }[/math]. Then find [math]\displaystyle{ h_0\left( t \right) }[/math] which has the maximum inner product with signal [math]\displaystyle{ s_0\left( t \right) }[/math] and
- [math]\displaystyle{ \left| B_p \right|^2 = \max_h \left| \left \langle s_p (t),h_p(t) \right \rangle \right|^2 }[/math]
Second, compute the residual:
- [math]\displaystyle{ s_1\left(t\right) =s_0\left(t\right)-B_0 h_0\left(t\right), }[/math]
and so on. It will comes out a set of residual ([math]\displaystyle{ s_p\left(t\right) }[/math]), projection ([math]\displaystyle{ B_p=\left \langle s_p(t),h_p(t) \right \rangle }[/math]), and elementary function ([math]\displaystyle{ h_p\left(t\right) }[/math]) for each different p. The energy of the residual will vanish if we keep doing the decomposition.
Energy conservation equation
If the elementary equation ([math]\displaystyle{ h_p\left(t\right) }[/math]) is designed to have a unit energy. Then the energy contain in the residual at the pth stage can be determined by the residual at p+1th stage plus ([math]\displaystyle{ B_p }[/math]). That is,
- [math]\displaystyle{ \left \| s_p(t) \right \|^2=\left \| s_{p+1}(t) \right \|^2 + \left| B_p \right|^2, }[/math]
- [math]\displaystyle{ \left \| s(t) \right \|^2=\sum_{p=o}^\infty \left| B_p \right|^2, }[/math]
similar to the Parseval's theorem in Fourier analysis.
The selection of elementary function is the main task in adaptive signal decomposition. It is natural to choose a Gaussian-type function to achieve the lower bound for the inequality:
- [math]\displaystyle{ h_p(t)=\left( \frac{\alpha}{\pi} \right)^\frac{1}{4}e^{ -\frac{\alpha}{2}(t-T_p)^2}e^{jt\Omega_p}, }[/math]
where [math]\displaystyle{ \left(T_p,\Omega_p\right) }[/math] is th mean and [math]\displaystyle{ \alpha_p^{-1} }[/math] is the variance of Gaussian at [math]\displaystyle{ \left(T_p,\Omega_p\right) }[/math]. And
- [math]\displaystyle{ s\left( t \right) = \sum_p B_p h_p(t) = \sum_p B_p\left( \frac{\alpha}{\pi} \right)^\frac{1}{4} e^{-\frac{\alpha}{2}(t-T_p)^2 }e^{jt\Omega_p} }[/math]
is called the adaptive Gabor representation.
Changing the variance value will change the duration of the elementary function (window size), and the center of the elementary function is no longer fixed. By adjusting the center point and variance of the elementary function, we are able to match the signal's local time-frequency feature. The better performance of the adaptation is achieved at the cost of matching process. The trade-off between different window length now become the trade-off between computation time and performance.
See also
References
- M.J. Bastiaans, "Gabor's expansion of a signal into Gaussian elementary signals", Proceedings of the IEEE, vol. 68, Issue:4, pp. 538–539, April 1980
- Shie Qian and Dapang Chen, "Signal Representation using adaptive normalized Gaussian functions," Signal Processing, vol. 42, no.3, pp. 687–694, March 1994
- Qinye Yin, Shie Qian, and Aigang Feng, "A Fast Refinement for Adaptive Gaussian Chirplet Decomposition," IEEE Transactions on Signal Processing, vol. 50, no.6, pp. 1298–1306, June 2002
- Shie Qian, Introduction to Time-Frequency and Wavelet Transforms, Prentice Hall, 2002
 |  |
