Accurate time-frequency estimation in sαs noise via memory-dependent derivative

This letter presents a time-frequency estimation approach based on memory-dependent derivative to obtain accurate spectrograph interpolation information. The memory correlation derivative is the convolution of a time-varying signal with a dynamic weighting function over a past time period with respect to a common derivative. Considering the described method, discrete data from previous times can be derived to estimate the signal values at the current time and to reduce the effect of noise. Fourier transforms with different scales and delay transforms are used as kernel functions to obtain energy-concentrated time-frequency curves with higher resolution and without frequency leakage. Besides, the memory-dependent derivative with adjustable scale factor is used to overcome time-frequency grid mismatches. Furthermore, differing from the phase accumulation manner of conventional time-frequency estimation, ℓ1-norm suppresses the heavy-tailed effect from outliers, thus the robustness of estimator can be enhanced greatly. By suitably choices of scale factor, the estimator can be tuned to exhibit high resolution in targeted regions of the time-frequency spectrum.


Introduction
Extracting the time-frequency trajectory of a non-linear frequency modulated signal from small samples with noise is a foundational problem in statistical signal processing.One challenging task is that the resolution of time-frequency cell is constrained by the time bandwidth product. 1 This indicates that the linear time-frequency analysis (TFA) cannot achieve an arbitrarily high timefrequency resolution at the same time.The other task is that TFA method designed for Gaussian noise performs poorly in the presence of outliers.The noise with heavytailed distribution, like impulsive noise, 2,3 will make the correlation function severely deviated from the theoretical value. 4Therefore, the motivation of this letter is that time-frequency estimator designed for large samples and Gaussian noise perform poorly in the presence of outliers in a small sample.
To improve the performance of TFA methods, many endeavors have been made in the past many years, and plenty of effective methods are proposed.][7] Time-frequency reassignment method, for example, synchrosqueezing transform (SST), belongs to a postprocessing technology, and it can improve the readability of time-frequency representation obviously. 8If the modulated frequency of a signal can be known in advance, the reassignment operator aims to sharpen the time-frequency representation by assigning the average of energy in a domain to the gravity center of these energy contributions.The parameterized TFA method, for example, general parameterized TFA (GPTFA), 9 is inspired by the fact that different analyzed window can achieve the different time-frequency resolution.9][20][21][22] Minimum ' 1 -norm optimization model has found extensive applications in linear parameter estimations.][25][26][27][28] Considering that impulse noise is instantaneous, the average power within a set time is adopted to reduce the effects of instantaneous impulses partially.Memorydependent derivative with adaptive kernel function guarantees exact localization for the time-frequency trajectory, so the SaS noise can be suppressed.][31][32][33] This paper is organized as following.In Section 2, the FM signal model is described.Section 3 discusses non-linear prediction model via memory-dependent derivative.In Section 4, time-frequency estimation with memory-dependent derivative is presented.Section 5 presents some processing results based on numerical simulations.The last section concludes this paper.

FM signal model
An analytic signal x 2 C L with the time-varying instantaneous frequency can be given by: where a 2 R + and u(u) is the amplitude and instantaneous frequency respectively, u stands for time, and i = ffiffiffiffiffiffi ffi À1 p .The conventional linear TFA is established on the assumption of the considered signal being piece-wisely linear frequency modulation in a short time.For example, in the short time ½u À t, u, the instantaneous frequency u(u) can be approximated as a constant, that is, where v 0 is the first-order approximation of the true frequency.
In fact, the window function needs to be short enough for strongly modulated frequency signal, otherwise, the instantaneous frequency is still time-varying.However, small samples within a short windows not only leads to a poor time-frequency resolution, but also sensitive to noise at a low signal-to-noise ratio (SNR) level.Specifically, the impulsive noise whose variance or energy is infinite causes large error for the signal within short window.Thus, this letter proposes an accuracy and robust time-frequency estimation method for non-linear FM signal.

Non-linear prediction model via memorydependent derivative
In this section, the memory-dependent derivative as a predict model of FM signal is derived from the integral of Taylor series at expansion points over an interval.Linear prediction is a mathematical operation where the current value of a linear frequency signal is estimated as a weighted function of previous samples.The most common representation is given by where v stands for time, x(t) is the predicted signal value, x(v) is the previous observed values, and h(t À v) is the kernel function.
Most prediction method can model linear frequency signal, but fails for FM signal.One of the simplest and effective specifications is to make x(t) piecewise linear.The Taylor series of FM signal converges for all t in the interval (v À s, v + s) (s.0 is infinitesimal) and the sum of the series is equal to x(t).
In this case, the first-order Taylor series expansion as a linear function approximates x(t) within interval (v À s, v + s), that is, Additionally, the second and higher order are ignored because the derivative amplifies the noise inevitably and the resulting inaccuracy worsens for higher derivatives. 10ext, as expansion point v slides in the short time interval v 2 ½t À t, t, we can obtain a series of linear functions converged to x(t).Liking the linear prediction method as shown in equation ( 4), the current value x(t) can be estimated by using the weighted sum of firstorder Taylor series expanded in the previous period v 2 ½t À t, t.The prediction value is an average of linear superposition of a series of local linear convergence values as follows Inserting equation ( 1) and x 0 (v) ' iv 0 x(v) into equation ( 5), it can be rewritten as which indicates that x(t) can be estimated by its previous values within ½t À t, t.
In fact, v 0 in equation ( 6) is unknown for the problem of instantaneous frequency estimation.Thus, , a kernel function h should be chosen properly to present the time-frequency characteristics of x.Enlightened by these, we have introduced the concept of memory-dependent derivative to predict the non-linear FM signal in a distinct way.Generalizing the kernel function of equation ( 6) to the general form, the memory-dependent derivative D t ½Á 11 can be defined as follows which is an integral form of x 0 (v) with a kernel function h(t À v) on a slipping interval ½t À t, t.A memorydependent derivative is characterized by its kernel function which can be chosen according to the frequency function of signal.

Time-frequency estimation with memorydependent derivative
The kernel function determines the accuracy of the prediction model for FM signal, so a memory-dependent derivative is characterized by its kernel function.A fundamental issue is to design an appropriate kernel function which improve accuracy of time-frequency estimation result.Given a rectangular window g(t À v) to truncate the signal to instead of the integral interval ½t À t, t, the spectrum expression of D t ½x(t) can be obtained through Fourier delay transformation as follows where F Á f g and Ã denotes the Fourier transform and convolution operator, respectively.According to the definition of Fourier transform, we define these Fourier transform pairs as follows Substituting equation ( 9) into equation ( 8), it can be rewritten as Traditionally, the TFA can discretize the continuous space to a finite number of time-frequency cells presented the time-frequency trajectory.This discretized model is simple and easy to handle analytically, but it inevitably brings the drawback that the restriction of uncertainty principle causes the energy dispersion.Especially, in the case of short window size, the conventional TF representation is not sufficient to generate a more energy concentrated results.
Utilizing the Fourier scale transformation, an kernel function with tunable grids is designed to generate the more energy-concentrated result by introducing the scale coefficient r.0.The kernel function is given by where v is estimated frequency.Its Fourier transform can be expressed as Obviously, the frequency domain grid can be refined for r.1.By using the scale transform, the basis mismatch problem can be resolved in exact location of instantaneous frequency.For signals with wellseparated frequencies, we show the r is roughly proportional to the number of frequencies, up to polylogarithmic factors.Instantaneous frequency estimation is possible even though kernel function is not inherently continuous at all and does not satisfy any sort of restricted isometry conditions.Inserting equation ( 12) into equation ( 10), we can obtain as follows Then, equation ( 13) is given by The energy of result in the time-frequency space divided into many small grid is well-concentrated, and it is not restricted by the Heisenberg uncertainty principle because it is based on an inner product operator with the refined kernel function.Furthermore, because a peak at the true frequency at v = v and some null for the side-lobes is presented at non-targeted cells, the estimated frequency v within window function can be calculated the derivative of TF x (v, vjr) with respect to v as follows As the window function g(t À v) moves, residuals of time-frequency trajectories in the whole time-frequency spectrum can be readily formulated as The ' 0 or ' 2 norm is sensitive to outlier observations.This sensitivity is the principal reason for exploring alternative norms.Procedures for the proposed method that involve the norm have been developed to increase robustness.The ' 1 norm has been applied in numerous variations of frequency estimation.5][36][37][38] Therefore, this equation ( 16) proposes an ' 1 norm procedure based on the efficient calculation of the optimal solution of the ' 1 norm best-fit hyperplane problem.Additionally, a proper scale factor can ensure a good concentration of energy in the time-frequency presentation.Thus, the optimal scale factor r opt can be determined by calculating the derivative of TF x (v, vjr) with respect to r as follows Due to both of TF x (v, vjr) and r are positive, the optimal scale factor r opt is obtained by solving equation ( 17) as follows Optimal scaling factor r opt is determined by the frequency deviation.Specifically, the estimated frequency is equal to the true, that is, v = v, at r opt !+ '.Considered an FM signal x(t) and a short and sliding rectangular window g(t À v), the estimated instantaneous frequency v at time v is obtained by minimizing equation ( 16) as follows From a statistical perspective, '1-norm is robust to impulsive noise. 12Equation ( 19) not only effectively mitigates the effect of artifacts caused by the noise, but also produces a time-frequency estimation for the nonlinear FM signal with an excellent concentration.

Numerical simulations
We compare the performance of the memorydependent derivative algorithm with several classical algorithms, and they are assessed by tests on generalized sinusoidal FM (GSFM) signal.The signal is given as The sample time is T = 1 S and the sample frequency is fs = 5000 Hz.4][15][16] Additionally, the average relative error of the estimated instantaneous frequency û½m with respect to the true instantaneous frequency is defined as follows Figure 1 displays the performance of SST, 8 GPTFA 9 and the proposed algorithm at a = 1:6 and GSNR = 25 dB (GSNR means generalized signal-tonoise ratio).The window size is 0.02 S. The frequency modulation rate of the FM signal can be estimated by minimizing equation (21).As shown in Figure 1(b), the energy of the signal would be scattered over a relatively wide frequency band, so their estimated IF trajectories cannot show themselves off very clearly.These estimated IF trajectories fluctuate around the true IF, even deviate from the true IF significantly at some moments in Figure 1(c) and (d).However, from the results, we can see that the proposed algorithm clearly outperform GPTFA and SST, especially at the moment of impulses appearance.Time-frequency estimation generated by the proposed algorithm is more excellent concentration because its time-frequency cell is adjustable adaptively.The relative error of GPTFA, SST and the proposed algorithm are 1.63%, 10.9% and 0.54%, respectively.This shows that the proposed method is suitable for heavy tailed and low GSNR environment.Next, Figure 2 illustrates the average relative error with the sample time for SST, GPTFA, and the proposed algorithm for 200 Monte-Carlo simulations.As shown in Figure 2(a), j(T ) by the proposed algorithm is very small and is insensitive to the SaS noise with various degrees of heavy tailed effect, that is, the characteristic parameter a 2 ½0:2, 2 every 0.2 at the fixed GSNR = 5 dB.Although j(T ) by GPTFA is insensitive as well to a, its error is much higher than the proposed algorithm.Besides, the estimated results of SST have larger bias as the characteristic parameter a increases.For non-Gaussian noise with high GSNR, results of the estimated IF by the proposed method and SST are similar for the whole data, but the proposed method is more accurate under a .1.2.For the different GSNR, that is, GSNR 2 [210, 10] dB every 2 dB at the fixed a = 1:6, the Figure 2(b) reveals that the proposed algorithm has a better estimation performance than the other two methods whether GSNR is high or low.The proposed method is superior to other two algorithms, which are sensitive to the low GSNR SaS noise.And the SST performance is seriously degraded under the condition of low GSNR.There is no lower bound for the estimation performance shown in Figure 2 due to the complex probability density of the alphastable noise.It is difficult to determine the algorithm's lower limit for frequency estimation.However, the error pattern is shown in Figure 2.

Conclusion
This letter presents the time-frequency estimation algorithm based on memory-dependent derivative for the exact and robust instantaneous frequency estimation of nonlinear FM signal.Unlike previous work in timefrequency estimation, the instantaneous frequency is not assumed to lie on a grid, but can assume any value in the frequency domain.Furthermore, the phase accumulation used by conventional method is converted to '1-norm to guarantee robust frequency localization in SaS noise.Simulation results demonstrate that the algorithm performs well at various levels of the heavytailed noise, and its estimation accuracy is high at low GSNR levels.