- Research Article
- Open access
- Published:
Robust Time-Frequency Distributions with Complex-Lag Argument
EURASIP Journal on Advances in Signal Processing volume 2010, Article number: 879874 (2010)
Abstract
The robust time-frequency distributions with complex-lag argument are proposed. They can provide an accurate estimation of fast varying instantaneous frequency in the presence of noise with heavy-tailed probability density function. The L-estimate form of this distribution is defined and it includes the L-estimate form of Wigner distribution as a special case. A modification for multicomponent signal representation is proposed, as well. Theoretical considerations are illustrated by the examples.
1. Introduction
Nonstationary signals such as speech, radar, seismic, sonar, and biomedical signals can be found in many practical applications. Due to time-varying spectra of these signals, time-frequency analysis has been used in their analysis. For different types of signals, various time-frequency distributions (TFDs) have been proposed [1–5].
In real applications we deal with signals corrupted by noise. If noise is additive with Gaussian probability density function (pdf), the standard time-frequency distributions represent a maximum likelihood (ML) estimate [6]. However, if the signal is corrupted by noise with heavy-tailed pdf (usually caused by environmental or human-made activities), the standard TFDs produce poor results. Consequently, the robust time-frequency distributions have been introduced [6–12]. The simplest and the most commonly used robust time-frequency representation is the robust short-time Fourier transform (STFT). The marginal median robust STFT has been introduced as an ML estimate of signals with Laplacian noise [8]. This form can also be successfully used for other types of heavy-tailed noises. The L-estimate robust STFT is introduced for signals with a mixture of Gaussian and impulse noises [9]. As in the case of the standard STFT, the main drawback of the robust STFT is a poor time-frequency resolution. In order to improve the time-frequency concentration, the robust forms of the Wigner distribution (WD) have been introduced [9–12]. They can provide an ideal concentration for signals with a linear instantaneous frequency (IF). However, for multicomponent signals the cross-terms appear. The robust S-method that combines good properties of the STFT and the WD has been introduced to provide a cross-terms free representation [12]. However, it cannot provide good concentration for signals with fast varying IF. Thus, the time-frequency distributions with complex-lag argument have been used to estimate nonlinear and fast IF variations [13–20]. Similarly as other TFDs, these distributions provide poor signal representation in the presence of heavy-tailed noise.
In this paper we propose a robust form of the N th-order complex-lag time-frequency distribution (CTD). An arbitrary high concentration can be achieved by increasing the distribution order N. The standard CTD has been defined as convolution of the WD and the Fourier transform of the higher order complex-lag moment, called concentration function (CF) [17, 18]. Similarly, the robust CTD can be obtained as convolution of the robust WD and CF forms. Additionally, a cross-terms free robust complex-lag time-frequency distribution is proposed for multicomponent signals.
The paper is organized as follows. The elementary theory behind linear and quadratic robust time-frequency representations is presented in Section 2. The N th order robust time-frequency distributions with complex-lag argument are proposed in Section 3. The advantages of the proposed distributions are proven through various examples in Section 4. Concluding remarks are given in Section 5.
2. Theoretical Background
2.1. Robust Short-Time Fourier Transform
Consider the noisy signal where is a complex-valued signal corrupted with complex-valued noise The STFT can be obtained as a solution of the following optimization problem:
where F(e) is the loss function, while is the error function:
The number of samples within the window is denoted as while represents complex-valued optimization parameter. In the case of the loss function , the standard STFT is obtained:
The quadratic loss function is the ML estimate of a signal corrupted with Gaussian noise. Hence, the standard STFT is applicable in this case. Nevertheless, when impulse noise is present, the quadratic loss function yields poor results. Therefore, other loss functions should be used. For instance, it has been shown that the loss function exhibits a robust behavior for heavy-tailed noise (e.g., Cauchy noise, Laplacian noise) [6, 11, 21]. The implicit solution of the optimization problem in (1) for requires a computationally demanding iterative procedure. In order to avoid such an iterative procedure, the marginal median estimate, with the loss function , has been introduced in [8]:
If the signal is corrupted by a mixture of Gaussian and impulse noise, the L-estimate robust STFT is used. It is defined as [9]
where the elements and are sorted in nondecreasing order as and , respectively. The coefficients a i are given as
where N s is even, while the parameter α takes values within the range For and the standard STFT and the marginal median STFT are obtained, respectively. Higher value of provides an enhanced reduction of heavy-tailed noise, while smaller value of improves spectral characteristics. Thus, depending on the application, the value of parameter should be chosen to provide good trade-off between these requirements.
The robust spectrogram is obtained as
where STFT h may be STFT M or STFT L .
2.2. Robust Quadratic Time-Frequency Distributions
The previous concept has been extended to the Wigner distribution. The WD can be obtained as a solution of the following optimization problem:
where denotes the complex conjugate of For the loss function the standard WD follows:
The median-based robust WD has been defined in [8] as a solution of optimization problem in (8), for the loss function :
The L-estimate robust WD can be written as [9]
where the elements are sorted in nondecreasing order as , while the coefficientsa i are defined by (6). The marginal median WD follows for =1/2.
The WD can be calculated by using STFT:
where STFT h can be the median-based STFT or the L-estimate STFT.
The robust S-method has been introduced in order to reduce or remove cross-terms in the robust WD. It can be written in the following form [12]:
where P(l) is a frequency domain window with its width equal to For and the robust spectrogram and the robust WD are obtained, respectively. More details about parameter selection can be found in [5].
3. Robust Complex-Lag Time-Frequency Distributions
The time-frequency distributions with complex-lag argument have been introduced for signals with fast varying instantaneous frequency. The general form of the standard N th-order complex-lag time-frequency distribution can be written as [18]
where the standard WD is given by (9), while CF(n,k) is the Fourier transform of the complex-lag signal moment:
where is an even number, representing the distribution order, while the quantity defines the equidistant roots on the unit circle. It has been shown [17–19] that the distribution defined by (14) can provide an arbitrarily high concentration by increasing N. Namely, the complex-lag distributions significantly reduce the spread factor produced by higher phase derivatives. For example, the fourth-order distribution:
contains the number of spread terms that is twice smaller than for the WD (Table 1). Furthermore, the sixth-order distribution is obtained for the roots :
and its spread factor is also given in Table 1. Note that, in general, the first term in the spread factor is of order (N+1).
In order to obtain c(n,m), the signal with complex-lag argument is calculated by using the concept of analytic extension as follows:
where , while X S (k) is the standard Fourier transform. The coordinate is multiplied by . The influence of this term can be such that an additional oversampling (or interpolations) of signal is required.
By analogy with the standard CTD, the robust CTD can be defined as
where RWD is the robust WD, while robust CF (RCF) is obtained as a solution of the optimization problem:
for the loss function The robust WD calculation has been already presented, while the robust CF can be obtained as a solution of nonlinear equation:
The iterative procedure for the nonlinear equation (21) is even more demanding than in the cases of robust STFT and robust WD calculations [6, 11]. Namely, to calculate the robust form of (18) has to be used. It requires an additional iterative procedure. Hence, the CF calculation requires nested iterative procedures, inappropriate for practical realization.
3.1. L-Estimate Form of the Robust CTD
The marginal median and the L-estimate approach can be used to overcome disadvantages of iterative procedure for the robust CTD calculation. In the sequel, only the L-estimate approach is considered, since the marginal median follows as a special case of the L-estimate forms for =1/2. Also, the L-estimates exhibit enhanced performance in the presence of mixture of Gaussian and impulse noise, common to real applications. Thus, the L-estimate approach is used to define the robust CTD.
Having in mind (19), the L-estimate robust CTD can be obtained as a convolution of the L-estimate robust WD and the L-estimate robust CF. By analogy with the robust WD, the L-estimate approach is used for the robust CF calculation, as follows:
where and are the elements of
respectively. They are sorted in nondecreasing order: and . The coefficients are given by (6), while rc(n,m) represent the robust complex-lag signal moment:
The numerical realization is simplified by using the angle and functions. Namely, for a signal in the form , the amplitude modulation terms that may appear in the calculation of and are eliminated, [18]. Also, calculation of signal raised to power j is avoided by using the exponential with function, [14]. The signal with the complex-lag argument is obtained as
where X L (k) represents the L-estimate of Fourier transform:
where and are such that and .
Finally, the L-estimate of complex-lag time-frequency distribution is:
where WDL represents the L-estimate robust WD.
3.2. Robust CTD Form for Multicomponent Signals
Note that the robust CTD form (defined by (27)) can be used for monocomponent signals. However, in the case of multicomponent signal: , the cross-terms appear. Thus, it is necessary to modify (27). The robust S-method, as a cross-terms free distribution, will be used instead of the robust WD, while a modification providing a cross-terms free robust CF should be introduced. In that sense, the signal with complex-lag argument is separately calculated for each component. The component separation is performed by using the robust STFT. Namely, the q th signal component is obtained as
where is the position of the q th signal component maximum in the L-estimate robust STFT. It is assumed that the q-th signal component is of 2W q +1 width; that is, it is within the region [k q (n)-W q , k q (n)+W q ]. Observe that the cross-terms will be avoided, if the distance between signal components is higher than 2W q (see [14] for details). After the first signal component is obtained, the values of STFT L (n,k) within the region [k q (n)-W q , k q (n)+W q ] will be set to 0. Then, this procedure is repeated for other components.
Furthermore, for the q th signal component the complex-lag signal moments and are
For all signal components we have and . By using the L-estimate approach, two corresponding robust CFs can be defined as
where and (sorted in nondescending order) are elements of:
respectively. Similarly, and are elements of:
where and .
The cross-terms free robust CF is:
Finally, the L-estimate robust CTD for multicomponent signals can be written in the form:
The role of window P(l) is same as in the S-method.
4. Examples
Highly nonstationary signals with fast varying instantaneous frequencies are considered. The signals are corrupted with heavy-tailed noise. The standard and the robust forms of CTD are considered and compared with corresponding forms of the Wigner distribution.
Example 1.
Consider a noisy signal:
where is heavy-tailed complex valued noise (cube of Gaussian noise):
where and are mutually independent Gaussian noises (zero mean with variance equal to 1). The time interval with sampling rate is used. The Gaussian window of width is applied in all cases.
The L-estimate forms are calculated by using parameter for all distributions. Namely, this value provides satisfying trade-off between noise reduction and distribution concentration. For a given signal, the standard WD, the L-estimate WD, the standard the marginal median (obtained according to (27) with α=1/2), and the L-estimate are shown in Figure 1. Since the signal has a fast varying IF, both the standard WD (Figure 1(a)) and the L-estimate WD (Figure 1(b)) yield poor results. The fourth-order CTD is introduced to improve concentration. However, as expected, the standard (Figure 1(c)) contains significant disturbances, due to the presence of strong heavy-tailed noise. Furthermore, the robust complex-lag distribution forms are considered. Thus, the results are improved by using the median (Figure 1(d)), but the accuracy of estimation is still not satisfactory. The best result in this example, regarding representation and IF estimation precision, is achieved by using the L-estimate (Figure 1(e)).
The mean squared error (MSE) is used as a quantitative measure of performance, for all distributions:
where is the true IF, while is the estimated IF: . The mean values of MSEs are given in Figure 2 for 100 realizationsof noises. Note that the L-estimate provides the lowest MSE.
In the presence of Gaussian noise, the performance of L-estimate is similar to the performance of the standard The MSEs of IF estimation, calculated as a mean value for 100 realizations of Gaussian noises, are given in Table 2.
Example 2.
This example aims to illustrate how the distribution order has to be increased to achieve concentration improvement. Namely, in the case of signal with IF variations that are faster than in the previous example, for example,
the fourth-order distribution is no longer optimal. Thus, the concentration can be improved by using the sixth-order distribution. The noise is the same as in the previous example. The L-estimate RWD, the L-estimate and L-estimate are shown in Figure 3. The parameter is used in all cases.
As expected, the L-estimate WD is not suitable for analysis. Significant improvements were obtained with the L-estimate while the best results are achieved by using the L-estimate in this example. The MSE of instantaneous frequency estimation is given in Table 3 for 100 realizations of noises. Note that the lowest MSE is obtained for the L-estimate distribution.
Example 3.
Consider a noisy multicomponent signal:
where and represent Gaussian noises. The same parameters for time interval, window, and noise strength are used as in the Example 2. The results for the standard SM, the L-estimate SM, the standard CTDN=4, and the L-estimate RCTDN=4 are shown in Figure 4.
Note that both the standard SM and the standard CTDN=4 cannot provide satisfactory results due to the presence of impulse noise. The L-estimate SM provides good estimation for the component with slow varying instantaneous frequency (Figure 4(b)). However, it cannot follow the fast IF variationsfor the second signal component. The L-estimate RCTDN=4 provides satisfying concentration for both components (Figure 4(d)).
5. Conclusion
The L-estimate-based robust N th-order complex-lag time-frequency distribution has been proposed. It provides an efficient estimation for nonstationary signals corrupted with a mixture of Gaussian and heavy-tailed impulse noise. Additionally, we proposed the modified L-estimate robust CTD form that provides a cross-terms free representation for multicomponent signals.
The L-estimate and standard distribution approaches could be combined in some future work to reduce the calculation complexity. Also, the future research could be focused to generalize the proposed approach to the class of complex-time distributions based on the ambiguity domain [20].
References
Cohen L: Time-Frequency Analysis. Prentice-Hall, Englewood Cliffs, NJ, USA; 1995.
Boashash B: Time-Frequency Analysis and Processing. Elsevier, Amsterdam, The Netherlands; 2003.
Boashash B: Estimating and interpreting the instantaneous frequency of a signal—part 1: fundamentals. Proceedings of the IEEE 1992, 80(4):520-538. 10.1109/5.135376
Hlawatsch F, Boudreaux-Bartels GF: Linear and quadratic time-frequency signal representations. IEEE Signal Processing Magazine 1992, 9(2):21-67. 10.1109/79.127284
Stanković LJ: A method for time-frequency analysis. IEEE Transactions on Signal Processing 1994, 42(1):225-229. 10.1109/78.258146
Katkovnik V: Robust M-periodogram. IEEE Transactions on Signal Processing 1998, 46(11):3104-3109.
Katkovnik V, Djurović I, Stanković LJ: Instantaneous frequency estimation using robust spectrogram with varying window length. AEU—International Journal of Electronics and Communications 2000, 54(4):193-202.
Djurovi I, Katkovnik V, Stanković LJ: Median filter based realizations of the robust time-frequency distributions. Signal Processing 2001, 81(8):1771-1776. 10.1016/S0165-1684(01)00092-5
Djurović I, Stanković LJ, Böhme JF: Robust L-estimation based forms of signal transforms and time-frequency representations. IEEE Transactions on Signal Processing 2003, 51(7):1753-1761. 10.1109/TSP.2003.812739
Katkovnik V, Djurović I, Stanković LJ: Robust time-frequency distributions. In Time-Frequency Signal Analysis and Applications. Edited by: Boashash B. Elsevier, Amsterdam, The Netherlands; 2003:392-399.
Djurović I, Stanković LJ: Robust wigner distribution with application to the instantaneous frequency estimation. IEEE Transactions on Signal Processing 2001, 49(12):2985-2993. 10.1109/78.969507
Djurović I, Stanković LJ, Barkat B: Robust time-frequency distributions based on the robust short time Fourier transform. Annales des Telecommunications/Annals of Telecommunications 2005, 60(5-6):681-697.
Stanković S, Stanković LJ: Introducing time-frequency distribution with a 'complex-time' argument. Electronics Letters 1996, 32(14):1265-1267. 10.1049/el:19960849
Stanković LJ: Time-frequency distributions with complex argument. IEEE Transactions on Signal Processing 2002, 50(3):475-486. 10.1109/78.984717
Morelande M, Senadji B, Boashash B: Complex-lag polynomial Wigner-Ville distribution. Proceedings of the 10th IEEE Annual Conference on Speech and Image Technologies for Computing and Telecommunications (TENCON '97), December 1997, Brisbane, Australia 1: 43-46.
Viswanath G, Sreenivas TV: IF estimation using higher order TFRs. Signal Processing 2002, 82(2):127-132. 10.1016/S0165-1684(01)00168-2
Cornu C, Stanković S, Ioana C, Quinquis A, Stanković LJ: Generalized representation of phase derivatives for regular signals. IEEE Transactions on Signal Processing 2007, 55(10):4831-4838.
Stanković S, Žarić N, Orović I, Ioana C: General form of time-frequency distribution with complex-lag argument. Electronics Letters 2008, 44(11):699-701. 10.1049/el:20080902
Stanković S, Orovic I: Effects of Cauchy integral formula discretization on the precision of IF estimation: unified approach to complex-lag distribution and its L-form. IEEE Signal Processing Letters 2009, 16(4):327-330.
Orović I, Stanković S: A class of highly concentrated time-frequency distributions based on the ambiguity domain representation and complex-lag moment. EURASIP Journal on Advances in Signal Processing 2009, 2009:-9.
Li TAH, Song KAIS: Estimation of the parameters of sinusoidal signals in non-Gaussian noise. IEEE Transactions on Signal Processing 2009, 57(1):62-72.
Acknowledgment
This work is supported by the Ministry of Education and Science of Montenegro.
Author information
Authors and Affiliations
Corresponding author
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
Žarić, N., Orović, I. & Stanković, S. Robust Time-Frequency Distributions with Complex-Lag Argument. EURASIP J. Adv. Signal Process. 2010, 879874 (2010). https://doi.org/10.1155/2010/879874
Received:
Accepted:
Published:
DOI: https://doi.org/10.1155/2010/879874