 Research
 Open Access
 Published:
On uncertainty principle of the local polynomial Fourier transform
EURASIP Journal on Advances in Signal Processing volume 2012, Article number: 120 (2012)
Abstract
In this article, a comprehensive study on uncertainty principle of the local polynomial Fourier transform (LPFT) is presented. It shows that the uncertainty product of the LPFT of an arbitrary order is related to the parameters of the signal and the window function, in addition to the errors of estimating the polynomial coefficients. Important factors that affect resolutions of signal representation, such as the window width, the length of overlap between signal segments, order mismatch and estimation errors of polynomial coefficients, are discussed. The effects of minimizing computational complexities on signal representation by reducing the order of the transform and the overlap length between signal segments are also examined. In terms of the signal concentration, comparisons among the shorttime Fourier transform, the WignerVille distribution and the second order LPFT are presented. The LPFT is shown to be an excellent candidate providing better representations for timevarying signals.
1 Introduction
The Fourier transform is an important tool to obtain signal representation in the frequency domain. However, it is not suitable to deal with timevarying signals that contain frequencies changing with time. To describe the frequency characteristics of such signals, the shorttime Fourier transform (STFT), defined as
is widely used.^{a} The basic assumption of the STFT is that the frequencies in each signal segment are not changed with time so that the frequency variations with time are approximately described by the Fourier transforms of the successive signal segments. Because this assumption is not generally true, however, the resolution of the STFT in the timefrequency domain is often limited.
To overcome this drawback, other timefrequency transforms such as fractional Fourier transform [1, 2] and local polynomial Fourier transform (LPFT) [3] were proposed. As a generalization of the STFT, the kernel of the LPFT uses extra parameters to approximate the phase of the signal into a polynomial form. With these parameters, or the polynomial coefficients, the LPFT can describe the timevarying frequencies with a better accuracy, and therefore the resolution of signal representation in the timefrequency domain can be significantly improved compared with that achieved by the STFT. Because the LPFT is a linear transform, it is free from the cross terms which may exist in some other timefrequency representations. Due to its desirable advantages, the LPFT has found applications in various areas. For example, it has been used to improve radar images of fast maneuvering targets [4, 5]. It has been applied for nonstationary interference suppression in noise radar systems [6], as well as in spread spectrum communications [7], to achieve improved performance compared with that obtained by the STFT. In [8], the LPFTbased beamformer is derived for source localization and tracking in nonstationary environment, which can resolve closely spaced sources when their velocities are sufficiently different. Moreover, the LPFT can be used as an instantaneous frequency estimator, and the corresponding asymptotic covariance matrix and bias of the estimates are studied in [3]. A comprehensive review on the LPFT and its applications can be referred to [9].
In the classical timefrequency domain, the product of time duration and bandwidth cannot be arbitrarily small and is generally described in the form of the uncertainty principle [10]. The uncertainty principle plays an important role in signal analysis [10], and studies on the uncertainty principle of various transforms have been reported in the literature [1, 11–16]. For example, it was reported that the STFT is limited by the uncertainty principle [17]. The local uncertainty principle was demonstrated in [12] by defining local quantities as conditional standard deviations, while the standard formulation of the uncertainty principle is in terms of global standard deviation. It shows that the local uncertainty product of the spectrogram indeed has a lower bound due to the windowing approach but it is not a property of the signal or a fundamental limit. For a large class of bilinear timefrequency distributions, the local uncertainty product is always less than or equal to the global uncertainty product and can be arbitrarily small [12].
It has been observed that the resolution of the LPFT in the timefrequency domain is also influenced by the window length which controls the tradeoff of bias and variance [3, 18]. Therefore, it is believed that there must exist some form of the uncertainty principle based on the LPFT. The uncertainty principle of the second order LPFT was only recently discussed in [19]. It was shown that the uncertainty product of the second order LPFT for the second order PPS is timeindependent when Gaussian window is used to segment the signal. However, a comprehensive study on the uncertainty principle for an arbitrary order LPFT has not been reported in the literature.
Similar to the steps used to derive the uncertainty principle of the STFT [17], this article firstly reports the mathematical derivations of the uncertainty principles of the second, third, fourth order and thereafter M th order LPFTs. It is found that the uncertainty product of an arbitrary order LPFT is related to the parameters of the signal and the window function, as well as the errors of estimating the polynomial coefficients. When the polynomial coefficients of the transform kernel are accurately estimated, the uncertainty product becomes a constant if a Gaussian window is used to segment the input signal. Based on the derived signal duration and bandwidth, we can obtain the uncertainty product. Then we will discuss the effects, in terms of the resolution in the timefrequency domain and the minimization of the required computational complexity, of window width, order mismatch, i.e., using the second order LPFT to process higherorder PPS, and the length of overlap between signal segments. Comparisons are also made on the resolutions of signal representations in the timefrequency domain achieved by using the STFT, the WignerVille distribution (WVD) and the second order LPFT.
The rest of the article is organized as follows. With the similar terminologies used in deriving the uncertainty principle of the STFT in Appendix [17], the uncertainty principles of the LPFTs of various orders are given in Section 2. Section 3 discusses the effects of various issues on obtaining desirable signal representations in the timefrequency domain. Signal concentrations of different timefrequency representations are also compared. Finally conclusions are given in Section 4.
2 Uncertainty principles of the LPFTs
In this section, the definition of the LPFT is given and some important issues that affect the computational complexity are addressed. Then we will give the derivation of the uncertainty principles of different order LPFTs and achieve the general expression of the uncertainty principle for the Mth order LPFT.
2.1 The LPFT
The LPFT is a generalization of the STFT, and the M thorder LPFT is defined as [3, 9]:
where M is the order of the polynomial function, ${\omega}_{1}=\frac{d\mathrm{\Omega}\left(t\right)}{dt}$, ${\omega}_{2}=\frac{{d}^{2}\mathrm{\Omega}\left(t\right)}{d{t}^{2}}$, .., and ${\omega}_{M1}=\frac{{d}^{M1}\mathrm{\Omega}\left(t\right)}{d{t}^{M1}}$ are the polynomial coefficients, and Ω(t) is the instantaneous frequency of the signal.
The LPFT is particularly suited to process the polynomial phase signals (PPSs) with a Gaussian amplitude defined by
where P is the order of the PPSs. To compute the LPFT in (2), the polynomial coefficients ω_{1}, ω_{2},.., ω_{M1}are first estimated and then used to form the polynomial phase. The LPFT is obtained by using an STFT procedure whose input is obtained by modulating s(τ) with $\mathsf{\text{exp}}\left\{j\omega tj{\sum}_{m=2}^{M}\frac{{\omega}_{m1}{\left(\tau t\right)}^{m}}{m!}\right\}$.
A few important issues, that have direct impacts on the resolution of timefrequency representation and the required computational complexity, have to be addressed. It is always desired to accurately estimate ω_{1}, ω_{2},.., and ω_{M1}so that a high resolution of the signal representation in the timefrequency domain can be obtained. The polynomial time frequency transform (PTFT) [20, 21], which is considered as a maximum likelihood estimator, can be used to estimate the polynomial coefficients used in (2). The PTFT is a transform which converts a 1D input signal into a multidimensional output from which the polynomial coefficients are estimated. The maximum number of polynomial coefficients that can be estimated is the same as the order of the PTFT. Figure 1 shows a 2D PTFT of two chirp components. From the location coordinates of the two peaks, the values of ω and ω_{1} for the two chirp components can be approximately achieved.
When M, the order of the LPFT, is large, one practical problem is that the corresponding PTFT requires a heavy computational load even using the fast algorithms reported in [22, 23]. Based on the radix2 fast algorithm for the PTFT [22] and the fast Fourier transform (FFT) algorithm [24–26] for the STFT computation, the number of complex multiplications needed by the LPFT is approximately
where N is the length of the input signal, q is the length of the window (or signal segment), and l, 0 ≤ l ≤ q  1, is the length of overlap between signal segments. In (4), ${C}_{\mathsf{\text{PTFT}}}^{r}\left(q\right)$ is the number of complex multiplications needed by the r th order PTFT of length q, the second term in the square brackets is the number of complex multiplications for an FFT of length q, the third term in the square brackets is the number of complex multiplications for windowing and modulation operations, and the factor outside the square brackets is approximately the number of the segments used for the LPFT. It was reported that for the r th order PTFT of length N, where r = 2 or 3, the required number of complex multiplications are in the order of N^{r} [22]. Because the PTFT requires a heavy computational complexity, especially when N is large, it is important to find various methods to minimize the computational complexity for the parameter estimation by using the PTFT.
Similar to the concept of the STFT which uses Fourier transform to deal with signal segments that contain timevarying frequencies, we can use a smaller number of the polynomial coefficients in the LPFT to approximate the frequency contents of the signal segments. If the second order LPFT is used to process higherorder PPSs, as shown in [27, 28] for example, the number of polynomial coefficients to be estimated is reduced, which directly leads to the reduction of the PTFT order, or equivalently the reduction of the required computational complexity. We will discuss the side effects of this order mismatch, i.e., the order of the LPFT is smaller than the order of the PPSs, on the resolution of the signal representation in the timefrequency domain. Other practical issues, which also exist in the STFT, including the effects of window width and the length of overlap between adjacent signal segments will be discussed. The computational complexity can be significantly reduced if the overlap length can be minimized without obviously degrading the resolution of the timefrequency representation.
2.2 Uncertainty products of the LPFTs
This subsection first presents the uncertainty principles for the lower order LPFTs, i.e., M = 2, 3 and 4. Then, a general expression of the uncertainty principle for the M th order LPFT is deduced.
It is assumed that the input signals are PPSs with the same orders of the LPFTs, i.e., M in (2) is equal to P in (3). The corresponding normalized local signal is
where
is a Gaussian window used to obtain the signal segments and a > 0 is the parameter to control the window width.
With the above defined input signal and the window function, the normalization factor, i.e., the integral in the denominator of (5), becomes
Therefore, the normalized local signal is expressed as
For an arbitrary order PPS, it is easy to derive
which is independent of the order of the PPS. With (32) in Appendix, the mean time, 〈τ〉_{ t }, of the normalized local signal becomes
Similarly, we also obtain
Therefore the time duration of the normalized local signal, which is defined in (33) in Appendix, is expressed as
Since both 〈τ〉_{ t } and ${T}_{t}^{2}$ are related to η_{ t }(τ) instead of η_{ t }(τ), they are independent of the order of the input signal. It means that the signal duration in the time domain is determined by the signal parameter α and window parameter a only.
Let us now compute the mean frequency 〈ω〉_{ t } and the bandwidth ${B}_{t}^{2}$ for the LPFTs of various orders. For the second order LPFT, the local signal is
where s(t), the input of the second order LPFT, is the second order PPS or the chirp signal, obtained by setting P = 2 in (3), and ω_{1} is estimated by using the second order PTFT. In this case, the Fourier transform of the local signal, s_{ t }(τ), becomes the second order LPFT of s(t).
The normalized local signal for P = 2 in (8) at the time t becomes
With the definition in (35), we have
in which the first term is related to the signal parameter a_{0}, the second term is related to ω_{1}, and the last term is a product of the deviation of ω_{1} from a_{1} and the mean time, 〈τ〉_{ t }.
Similarly, by using (37) in Appendix, we have
which shows that the frequency variance is determined by the parameters of the signal and the window in the first term, the square of mean frequency in the second term, and the scaled square of the estimation error of the polynomial parameter ω_{1} in the last term. With the above derivation results, the bandwidth in the domain of the second order LPFT becomes
Therefore
which is related to the parameters of the signal and the window, and the errors of estimating the polynomial parameters ω_{1}.
Let us next consider the uncertainty principle of the third order LPFT. The local signal is
where s(t), the input of the third order LPFT, is the third order PPS obtained by setting P = 3 in (3).
In this case, the local normalized signal becomes
With the same procedure, the mean frequency, defined in (35) in Appendix, is expressed as
and 〈ω^{2}〉_{ t }, defined in (37) in Appendix, becomes
Thus, the bandwidth becomes
Therefore
It is observed that the uncertainty product, ${B}_{t}^{2}{T}_{t}^{2}$, is affected by the errors of estimating the polynomial parameters ω_{1} and ω_{2} in addition to the parameters of the signal and the window function.
By setting P = 4 in (3) and using the fourth order LPFT, we can similarly obtain the corresponding mean frequency 〈ω〉_{ t } and the bandwidth ${B}_{t}^{2}$, which are listed in Table 1 respectively. From the expressions in the table, we can observe the rules of 〈ω〉_{ t } and ${B}_{t}^{2}$ as M, the order of the LPFT, is increased. The mean frequency and the bandwidth for the M thorder LPFT are
The lengthy mathematical derivation of the above two equations is omitted due to the space limit. It is observed that both 〈ω〉_{ t } and ${B}_{t}^{2}$ are related to the moments 〈τ ^{n}〉_{ t }. The expressions of the moments 〈τ ^{n}〉_{ t }, for n = 1, 2, .., 6, which can be calculated from (34), are listed in Table 2.
In general the uncertainty product of an arbitrary order LPFT is related to the parameters of the signal and the window function, as well as the errors of estimating the polynomial coefficients. When the polynomial coefficients are estimated correctly for the M th order LPFT, that is, ${\omega}_{1}=\frac{d\mathrm{\Omega}\left(t\right)}{dt}$, .., ${\omega}_{M1}=\frac{{d}^{M1}\mathrm{\Omega}\left(t\right)}{d{t}^{M1}}$, the mean frequency becomes
which is exactly the instantaneous frequency of the M th order PPS.
Without giving further details, the bandwidth in the domain of the M th order LPFT is derived as ${B}_{t}^{2}=\frac{a+\alpha}{2}$ when the polynomial parameters are estimated correctly, which leads to the uncertainty product to be
Based on the above derived results, we have the following conclusions. When using the Gaussian window to segment the signal, the signal duration of the PPS is only related to the signal parameter α, the window parameter a and is independent of the order. For M > 1, the bandwidth ${B}_{t}^{2}$ is related to the parameters of the signal and the window function, as well as the errors of estimating the polynomial parameters ω_{1}, ω_{2} ..., and ω_{M1}. When the polynomial coefficients are estimated accurately, the uncertainty product of the M th order LPFT, with the same order PPS as the input, is independent of time or frequency, and becomes a constant. When there are errors of estimating the polynomial coefficients, the bandwidth together with the uncertainty product, are polynomial functions of these errors, which directly affect the signal representation in the timefrequency domain.
It should be noted that in this article the uncertainty principle of the LPFT is derived based on the PPSs. For other kind of timevarying signals such as sinusoidal FM signals or parabolic FM signals, we can use the window function to segment the signals to ensure that each segment within the window can be assumed to be the PPS. As long as the signal segments can be modeled as PPSs, the above presented results are still valid. In this way, the derived uncertainty principle of the LPFT can also be extended to general timevarying signals.
3 Discussions
The main objective on studying various issues of the LPFT is to achieve the desired signal resolution in the timefrequency domain with minimum computational costs. This section considers a few important issues, in terms of the uncertainty principles of the LPFTs, that have direct impacts on signal resolution and the required computational complexity.
3.1 Window width effects
The uncertainty relations derived in the previous section show that there exists a tradeoff between the resolutions in the time and frequency domains. With the assumption that the polynomial coefficients used in the LPFT are accurately estimated, the signal duration ${T}_{t}^{2}$ is inversely proportional to the sum of window width a and the signal parameter α, while the signal bandwidth ${B}_{t}^{2}$ is directly proportional to this sum, regardless of the order of the PPS and the LPFT. With different widths of the Gaussian window, Figure 2 shows the signal concentrations in the timefrequency domain achieved by using the second order LPFT to process the PPS which contains chirp components.
It is seen that, as a decreases or the window width increases, the chirp components become more concentrated in the frequency domain, or equivalently, the resolution of the signal representation in the frequency direction is increased. In contrast, the resolution in the time direction decreases. This observation is consistent with the signal duration ${T}_{t}^{2}$ and the bandwidth ${B}_{t}^{2}$ derived in the last Section. For example, decreasing the window parameter a leads to the increase of the signal duration ${T}_{t}^{2}$, or equivalently, the increase of signal spread in time. At the same time, the bandwidth ${B}_{t}^{2}$ decreases as a decreases so that the signal spread in frequency is reduced, that is, the resolution of the signal representation in the frequency direction is increased. Because it is impossible to increase the resolution in both domains at the same time, a compromise has to be made to balance the requirements of the signal representation in the timefrequency domain. When a = 2.5, for example, Figure 2c provides acceptable resolutions for both time and frequency domains, that is, the resolution in the time domain and the frequency domain is high enough to help us clearly separating the closely located chirp components.
3.2 Order mismatch effects
As discussed in the Section 1, it is desired to use lower order LPFTs to deal with higherorder PPSs to significantly reduce the required computational complexity. However, this mismatch between the orders of the LPFT and the PPSs will affect the performance of the signal representation in the timefrequency domain. It is necessary to find out how the signal representation in the timefrequency domain is affected and if possible, how the mismatch effects are to be minimized. Because it is difficult to draw a general conclusion for all possible cases, let us consider the use of the second order LPFT to process the third order PPSs.
By setting ω_{2} = 0 in (17) and (19), the mean frequency and the bandwidth become
When the parameter ω_{1} is estimated correctly, i.e., ω_{1} = a_{1} + a_{2}t, we have
Therefore, the uncertainty product becomes
In (23) or (24), the last term is the effects of the order mismatch, which are directly proportional to the polynomial coefficient, a_{2}. Therefore, the order mismatch may have a small effect when a_{2} is small enough, which is generally true in many practical applications. Another way to reduce the effect of the order mismatch is to minimize the signal duration ${T}_{t}^{2}$ by decreasing the window width or increasing the window parameter a. For example, when a → ∞, 〈ω〉_{ t } in (23) is approaching to ${a}_{0}+{a}_{1}t+\frac{{a}_{2}}{2}{t}^{2}$, which is the instantaneous frequency of the third order PPS. Meanwhile ${B}_{t}^{2}$ in (24) is approaching to $\frac{a+\alpha}{2}$. It means that, when the third order PPS is processed with the second order LPFT, 〈ω〉_{ t } and ${B}_{t}^{2}$ of the second order LPFT approach to those of the third order LPFT as the window width is decreased. The uncertainty product in (25) also approaches to the minimum.
When the signal parameter α approaches to zero, the third order PPS with a Gaussian amplitude becomes a third order PPS with a constant amplitude. Under this situation, the mean frequency and bandwidth become
and the uncertainty product is
which increases with the increase of the polynomial coefficient a_{2} and with the decrease of the window parameter a. Similarly, as a goes to infinity, 〈ω〉_{ t } approaches to the instantaneous frequency, and ${B}_{t}^{2}{T}_{t}^{2}$ approaches to $\frac{1}{4}$ which is the minimum of the uncertainty product.
From the above discussion, it is seen that increasing the window parameter a always helps to achieve a better estimation of the instantaneous frequency 〈ω〉_{ t }. It should be noted that, however, as a goes to infinity, that is as the window becomes narrower in the time domain, the signal bandwidth ${B}_{t}^{2}$, or the spread in frequency in (24) also approaches to infinity, which decreases the resolution in frequency. Therefore, the window width for practical applications should be properly selected to achieve a compromise between the resolutions in both time and frequency. An adaptive procedure for window length selection has been provided in [28]. Similarly, an automated procedure for window width selection is possible and will be discussed in our future study.
3.3 Effects of overlap lengths
It was shown in [28] that when the second order LPFT is used to process the chirp signals, the length of overlap between the adjacent signal segments can be significantly reduced without obviously degrading the signal resolution in the timefrequency domain, which allows us to further minimize the required computational complexity. For example, the required number of complex multiplications in (4) is reduced as l increases. For the third order PPS,
with α = 10^{4}, Figure 3a shows the signal representation obtained by the third order LPFT without any overlap. Compared to the performance with the maximum overlap, i.e., l = q 1, as shown in Figure 3b, no obvious degradation in signal resolution is observed compared to Figure 3a. Because the total computational complexity is directly proportional to the number of signal segments being processed, the LPFT computation without overlap between signal segments makes significant savings on the computational complexity compared to that using the overlap between signal segments.
Let us next consider the effects of using lower order LPFT with different overlap lengths. Figure 4a shows the representation of the third order PPS obtained by using the second order LPFT without any overlap. Obvious unsmoothness is observed at the junctions of two adjacent segments where the values of the estimated chirp rates have sudden changes. This order mismatch effect can be minimized by increasing the length of overlap between the adjacent segments. As shown in Figure 4b, the second order LPFT with the maximum overlap achieves a representation that is comparable to that obtained by using the third order LPFT in Figure 3b.
To process higherorder PPSs with the second order LPFT, another method for improvement on signal representation is to reduce the window length to ensure that each segment can be approximately assumed as a chirp component. Compared with Figure 4a, c shows the improvement on the smoothness of signal representation, which is achieved by reducing the window length from 1/4 to 1/8 of the signal length. It should be noted that, compared with the third order LPFT with the small window length in Figure 4d, the second order LPFT with the same window length in Figure 4c can achieve comparable result and with less computational complexity. It is worth mentioning that for any order LPFTs, the frequency resolution will become degraded as the window length is reduced, as shown in Figure 4c, d. For other higherorder PPSs, the second order LPFT can also be used, as long as the window length is small enough to ensure that each segment within the window can be approximately assumed as the chirp signal.
3.4 Resolution comparisons
Let us now consider the performances achieved by the STFT and the second order LPFT of a chirp signal obtained from (3) with P = 2. Using the STFT, the corresponding mean frequency and bandwidth are [29]
The mean time 〈τ〉_{ t } and the duration ${T}_{t}^{2}$ of the STFT are the same as those of the LPFT. Therefore the uncertainty product of the STFT is
When a_{1} = 0, the input signal becomes one whose frequencies do not change with time and the STFT achieves the lowest bound of the uncertainty product ${B}_{t}^{2}{T}_{t}^{2}=\frac{1}{4}$. Therefore, the STFT can provide good resolution only for signals with constant frequencies. For the chirp signals, however, ${B}_{t}^{2}$ in (29) has a nonzero second term. Therefore, the uncertainty product ${B}_{t}^{2}{T}_{t}^{2}$ of the STFT of timevarying signals must be larger than $\frac{1}{4}$. For the same input signal, if the polynomial coefficient, ω_{1}, of the LPFT is correctly estimated, we achieve the lowest bound of the uncertainty product ${B}_{t}^{2}{T}_{t}^{2}=\frac{1}{4}$, as seen in (22). Therefore the second order LPFT achieves a more concentrated distribution in the timefrequency domain than the STFT.
Another important issue on the LPFT is about the estimation errors on ω_{1}. Although the PTFT is the maximum likelihood estimator, it is still possible to have some estimation errors, especially when the signals are corrupted by heavy noises. Because the second term of (15) is proportional to (a_{1}  ω_{1})^{2}, the uncertain product for the second order LPFT of chirp signals is ${B}_{t}^{2}{T}_{t}^{2}>\frac{1}{4}$ if ω_{1} is not correctly estimated. Comparing (15) and (29) reveals that as long as a_{1}  ω_{1}< a_{1}, i.e., the estimation error of ω_{1} is smaller than a_{1}, the LPFT of the chirp signal achieves a smaller bandwidth than the STFT. Therefore, the second order LPFT can still achieve a more concentrated distribution in the timefrequency domain than the STFT. Furthermore, the resolution performance of the STFT cannot be further improved because a_{1} in (29) or (30) is the signal parameter. In comparison, improvement on the resolution performance of the LPFT can be achieved by minimizing the error of parameter estimation. It should be noted that if the estimation error of ω_{1} is larger than a_{1}, the resolution of the LPFT will become worse than that of the STFT although such a case rarely happens when the PTFT is used to estimate the polynomial coefficients.
Besides the STFT and LPFT, which are linear timefrequency transforms, another bilinear timefrequency transform, known as the WVD, has been often used in many applications [29]. The WVD of the chirp signal is
When is small, the WVD is concentrated on ω = a_{0} + a_{1}t. The associated mean frequency and bandwidth are [29]
It means that the WVD can provide a better frequency resolution for chirp signals as far as α is small enough. However, the WVD suffers from the cross terms for multiplecomponent signals. It is difficult to remove the cross terms without sacrificing the resolution of the signal representation, particularly for nonlinear multiplecomponent chirp signals [29].
In order to clearly show the frequency resolution of the signal, Figure 5 compares signal concentration obtained from the STFT, LPFT and WVD at a particular time instant. Resolution usually means whether we can separate two closely located components in the time and/or frequency domains, while concentration describes how the transforms concentrate along the signal's instantaneous frequency. Resolution and concentration are closely related and high concentration generally contributes to high resolution. For signals with constant frequency components, the STFT can provide the good frequency resolution, as shown in Figure 5. For chirp signals, with the same window length and window parameter a, however, the LPFT achieves a better frequency resolution as seen in Figure 5b, c. Although the WVD also has a better frequency resolution, the cross terms between the two chirp components, as shown in Figure 5d, can be easily mistaken as a valid frequency component. Therefore, the second order LPFT is the best choice for dealing with chirp signals since it is free from the cross terms and also achieves a better resolution than the STFT.
4 Conclusion
In this article, the LPFT is shown to be limited by the uncertainty principle, and the uncertainty relations of various order LPFTs are derived to show the tradeoff between the resolutions of signal representation in the time and frequency domains and the required computational complexity. The uncertainty product of an arbitrary order is determined by the signal parameters, the window function and the errors of estimating the polynomial coefficients. When Gaussian window function is employed to segment the signals, the uncertainty products of the LPFT are time independent under the condition that the polynomial coefficients are accurately estimated. The effects of the window width, the estimation errors, and the order mismatch by using the second order LPFT to process higherorder PPSs are discussed. Resolution comparisons with the STFT and WVD are also provided to show the merits of the LPFT, and the LPFT is found to be a better tool to dealing with timevarying signals.
Appendix
Review on uncertainty principles of the STFT
This appendix provides an example of deriving the uncertainty principle for the STFT. The definitions of some terminologies in this appendix are also used in the main sections of this article. Let us consider the P thorder PPS defined in (3) as the input signal. By multiplying s(t) with a window function h(t), the local signal is defined as s_{ t }(τ) = s(τ)h(t  τ). It is easily seen that the Fourier transform of the local signal s_{ t }(τ) is exactly the STFT of the original signal s(t). The normalized local signal at time t is
which ensures that for any t, ∫η_{ t }(τ)^{2}dτ = 1. Let us review some terminologies based on the normalized local signal η_{ t }(τ) [29] which is a function of τ within a short duration, and its Fourier transform is defined as
to reveal the spectral information around the time instant t. For this normalized local signal, all the relevant quantities, such as the mean time duration, mean frequency and bandwidth, are time dependent [29].
The mean time of the normalized local signal is defined as
and its time duration is defined as
where 〈τ^{2}〉_{ t } is the second order moment of the normalized local signal, defined as
Similarly, the n th order moment of the normalized local signal is defined as
The mean frequency of the normalized local signal is defined as
and its bandwidth is
where
Here the duration ${T}_{t}^{2}$ and bandwidth ${B}_{t}^{2}$ are good measures of the broadness of the local signal in the time and frequency domains, respectively. For example, ${T}_{t}^{2}$ (or ${B}_{t}^{2}$) of a signal means the width of the signal concentrated around the mean time 〈τ〉_{ t } (or mean frequency 〈ω〉_{ t }). The uncertainty product in the STFT domain is generally a function of time, the signal s(t) and the window h(t), and has a lower bound [29]
It should be noted that the focus of this article is on the uncertainty product obtained by multiplying the duration and bandwidth of the local signal. It is important to understand that this uncertainty product places limits on the processing techniques of the windowed transforms. Other kinds of uncertainty products, such as the global uncertainty products, can be referred to [17].
Endnote
^{a}For simplicity in the rest of the article, the integral without limits implies that the integration is from ∞ to ∞.
Abbreviations
 STFT:

shorttime Fourier transform
 LPFT:

local polynomial Fourier transform
 WVD:

WignerVille distribution
 PPSs:

polynomial phase signals
 PTFT:

polynomial time frequency transform
 FFT:

fast Fourier transform.
References
 1.
Ozaktas H, Zalevsky Z, Kutay M: The Fractional Fourier Transform With Applications in Optics and Signal Processing. Wiley, Chichester; 2001.
 2.
Zhang F, Chen YQ, Bi G: Adaptive harmonic fractional Fourier transform. IEEE Signal Process Lett 1999, 6(11):281283.
 3.
Katkovnik V: A new form of Fourier transform for timefrequency estimation. Signal Process 1995, 47(2):187200. 10.1016/01651684(95)001077
 4.
Djurovic I, Thayaparan T, Stankovic L: Adaptive local polynomial Fourier transform in ISAR. EURASIP J Appl Signal Process 2006, 2006: 115.
 5.
Li X, Bi G, Ju Y: Quantitative SNR analysis for ISAR imaging using local polynomial Fourier transform. IEEE Trans Aerosp Electron Syst 2009, 45(3):12411248.
 6.
Dakovic M, Thayaparan T, Djukanovic S, Stankovic L: Timefrequencybased nonstationary interference suppression for noise radar systems. IET Radar Sonar Navig 2008, 2(4):306314. 10.1049/ietrsn:20070137
 7.
Stankovic L, Djukanovic S: Order adaptive local polynomial FT based interference rejection in spread spectrum communication systems. IEEE Trans Instrum Meas 2005, 54(6):21562162. 10.1109/TIM.2005.858567
 8.
Katkovnic V, Gershman A: A local polynomial approximation based beamforming for source localization and tracking in nonstationary environments. IEEE Signal Process Lett 2000, 7(1):35.
 9.
Li X, Bi G, Stankovic S, Zoubir A: Local polynomial Fourier transform: a review on recent development and applications. Signal Process 2011, 91: 13701393. 10.1016/j.sigpro.2010.09.003
 10.
Gabor D: Theory of communication. J Inst Electr Eng 1946, 93(26):429457.
 11.
Korn P: Some uncertainty principles for timefrequency transforms of the Cohen class. IEEE Trans Signal Process 2005, 53(2):523527.
 12.
Loughlin P, Cohen L: The uncertainty principle: global, local, or both? IEEE Trans Signal Process 52(5):12181227.
 13.
Shinde S, Gadre V: An uncertainty principle for real signals in the fractional Fourier transform domain. IEEE Trans Signal Process 2001, 49(11):25452548. 10.1109/78.960402
 14.
Mustard D: Uncertainty principle invariant under fractional Fourier transform. J Austr Math Soc Ser B 1991, 33: 180191. 10.1017/S0334270000006986
 15.
Sharma K, Joshi S: Uncertainty principle for real signals in the linear canonical transform domain. IEEE Trans Signal Process 2008, 56(7):26772683.
 16.
Xu G, Wang X, Xu X: The logarithmic, Heisenberg's and shorttime uncertainty principles associated with fractional Fourier transform. Signal Process 2009, 89(3):339343. 10.1016/j.sigpro.2008.09.002
 17.
Cohen L: The uncertainty principles of the shorttime Fourier transform. Proc SPIE 1995, 2563: 8090.
 18.
Katkovnikm V: Discretetime local polynomial approximation of the instantaneous frequency. IEEE Trans Signal Process 1998, 46(10):26262637. 10.1109/78.720366
 19.
Li X, Bi G: Uncertainty principle of the secondorder LPFT. In IEEE International Symposium on Circuits and Systems (ISCAS'2009). Taipei, Taiwan; 2009:325328.
 20.
Ikram M, AbedMeraim K, Hua Y: Fast quadratic phase transform for estimating the parameters of multicomponent chirp signals. Dig Signal Process 1997, 7(2):127135. 10.1006/dspr.1997.0286
 21.
Xia X: Discrete chirpFourier transform and its applications to chirp rate estimation. IEEE Trans Signal Process 2000, 48(11):31223133. 10.1109/78.875469
 22.
Bi G, Wei Y: Radix2 DIF fast algorithms for polynomial time frequency transforms. IEEE Trans Aerosp Electron Syst 2006, 42(4):15401546.
 23.
Ju Y, Bi G: Generalized fast algorithms for the polynomial timefrequency transform. IEEE Trans Signal Process 2007, 55(10):49074915.
 24.
Bi G, Chen Y: Fast DFT algorithms for length N = q * 2 ^{m} . IEEE Trans Circ Syst IIAnalogu Dig Signal Process 1998, 45(6):685690. 10.1109/82.686687
 25.
Bi G, Chen Y, Zeng Y: Fast algorithms for generalized discrete Hartley transform of composite sequence lengths. IEEE Trans Circ Syst IIAnalogu Dig Signal Process 2000, 47(9):893901. 10.1109/82.868457
 26.
Bi G, Chen Y: Fast generalized DFT and DHT algorithms. Signal Process 1998, 65(3):383390. 10.1016/S01651684(97)00234X
 27.
Djurovic I: Robust adaptive local polynomial Fourier transform. IEEE Signal Process Lett 11(2):201204.
 28.
Wei Y, Bi G: Efficient analysis of timevarying muliticomponent signals with modified LPTFT. EURASIP J Appl Signal Process 2005, 2005(1):12611268.
 29.
Cohen L: TimeFrequency Analysis. Prentice Hall, Englewood Cliffs, NJ; 1995.
Acknowledgements
This research work was supported by the National Natural Science sFoundation of China (No. 61102164), and the Scientific Research Fund of Hangzhou Normal University (No. 2011QDL021).
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Li, X., Bi, G. & Li, S. On uncertainty principle of the local polynomial Fourier transform. EURASIP J. Adv. Signal Process. 2012, 120 (2012). https://doi.org/10.1186/168761802012120
Received:
Accepted:
Published:
Keywords
 uncertainty principle
 shorttime Fourier transform
 local polynomial Fourier transform
 WignerVille distribution