A New Approach for the Characterization of Nonstationary Oscillators Using the Wigner-Ville Distribution

Levy, Chagai; Pinchas, Monika; Pinhasi, Yosef

doi:https://doi.org/10.1155/2018/4942938

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusion Appendix Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2018 | Article ID 4942938 | https://doi.org/10.1155/2018/4942938

A New Approach for the Characterization of Nonstationary Oscillators Using the Wigner-Ville Distribution

Chagai Levy,¹Monika Pinchas,¹and Yosef Pinhasi¹

Academic Editor: Akhil Garg

Received28 Feb 2018

Revised09 Jun 2018

Accepted20 Jun 2018

Published11 Jul 2018

Abstract

Oscillators and clocks are affected by physical mechanisms causing amplitude fluctuations, phase noise, and frequency instabilities. The physical properties of the elements composing the oscillator as well as external environmental conditions play a role in the characteristics of the oscillatory signal produced by the device. Such instabilities demonstrate frequency drifts and modulation and spectrum broadening and are observed to be nonstationary processes in nature. Most of tools which are being used to measure and characterize oscillator stability are based on signal processing techniques, assuming time invariance during a temporal window, during which the signal is assumed to be stationary. This paper proposes a new time-frequency metric for the characterization of frequency sources. Our technique is based on the Wigner-Ville distribution, which extends the spectral measures to consist of the temporal nonstationary behavior of the processes affecting the accuracy of the clock. We demonstrate the use of the technique in the characterization of phase errors, frequency offsets, and frequency drift of oscillators.

1. Introduction

In recent years, due to the high data rates required to be transferred in communication networks, the demand for high accuracy satellite navigation systems and the development of high resolution radars and clocks with very high precision are required. Stable oscillators, including atomic ones like rubidium and even cesium atomic, are often used in communication and navigation facilities [1–6]. In addition to the above applications, many other electronic systems require frequency generators. In such systems, time synchronization and an accurate clock frequency are required [7–12].

Every oscillator has its unique clock frequency consisting of a nominal frequency in addition to a normalized frequency difference of Parts Per Million (PPM) from the nominal frequency multiplied by this nominal frequency. In addition, every oscillator has a random jitter . Therefore, the master and the slave clocks may have different frequencies and thus not be synchronized. For example, suppose that the master’s frequency is () and the slave’s frequency is (); thus, their frequencies are not synchronized. In that case, the clock offset generally keeps increasing and may cause the communication link to fail [13]. Therefore, in such systems, frequency synchronization and an accurate clock frequency are required [14].

For example [13], communication networks such as Long Term Evolution Advanced (LTE-A) based cellular networks require base stations and backhaul equipment to maintain time synchronization in order to provide several new features like synchronized transmissions over frequencies from adjacent base stations and interference coordination. An example of an emerging networked application is smart power grid systems which are characterized by a two-way flow of energy and end-to-end communications. The communications are largely machine-to-machine (M2M) in nature [13, 15] and need to share synchronized information in order to improve efficiency and reliability of power delivery. Therefore, in addition to frequency synchronization, time synchronization is also required [16].

Oscillators and clocks are subjected to environmental conditions, internal malfunctions, and inherent physical phenomena causing instabilities in the oscillatory signal being produced [17]. These instabilities are expressed by amplitude and phase fluctuations causing frequency deviations, drift, and spectrum broadening. Instead of being stationary, single frequency, and temporal coherent, the clock signal is observed to be a stochastic process characterized by statistical and spectral features, which are nonstationary time varying [18–21]. Therefore, the design of a set of tools for time domain analysis of a frequency source is therefore needed and significant to investigate.

Several works [22–25] on the time-frequency (TF) analysis of oscillators and atomic clocks have been carried out recently. In [22, 23], the authors presented methods for the detection and identification of atomic clock anomalies via simulation on experimental data. Those methods are based on derivation of the frequency deviation and then sliding Welchs periodogram on the clock noise data.

Another work [24] presents a time-frequency relationship between the Langevin equation and the harmonic oscillator by deriving relationship between the Wigner distribution of the Green’s function of the Langevin equation and of the harmonic oscillator. Their result paves the way for a simplification of the time-frequency representation of differential equations, as well as for a better understanding and filtering of the interference terms.

Most recently, a study of the important role of time-frequency signal representations was presented [25]. This approach was used for enhancement of the performance of global navigation satellite system (GNSS) receivers by using short-time Fourier transform (STFT) and the Wigner distribution as the TF representations.

In order to measure and characterize the oscillator properties, time-frequency signal processing techniques are required. Most often the Allan variance (or the Allan deviation) is employed [26–32]. It calculates the covariance of two adjacent samples of the signal produced by the under test oscillator. The Allan variance is widely used and recommended in international standards [28]. However, the Allan variance does not give the possibility of identifying and interpreting nonstationary effects.

Recently [18–21], an extended version of the Allan variance, the dynamic Allan variance (DAVAR), was proposed for the characterization of atomic clock stability that can vary with time. However, the DAVAR deals only with characterization of the stability of the frequency source and cannot detect phase and frequency offsets relative to a reference source.

The Allan variance is a fundamental tool for characterizing frequency stability of oscillators since it measures the size of frequency fluctuations at different observation intervals, whereas the DAVAR describes the variations with time of these fluctuations. Please note that in this paper we concentrate in ideal oscillators, where no random frequency fluctuations are considered.

In this paper we propose a time-frequency metric for the characterization of oscillators and frequency generators, considering also nonstationary variations. It is based on the Wigner distribution [33, 34], proposed in 1932 for the characterization of quantum fluctuations. In signal processing, it was used by Ville [35], so that it is often called Wigner-Ville distribution [36]. The applications of Wigner distribution are various: analysis of nonstationary signals [37–40], radar signals [41, 42], biomedical signals [43–45], analysis of time varying filters [46, 47], and image processing [48–50]. However, the Wigner-Ville distribution was not yet used in the characterization of nonstationary oscillators such as detection of the phase and frequency offsets relative to a reference source. In this paper, we will introduce both the discrete and continuous forms of the Wigner-Ville distribution for the characterization of a clock signal (nonstationary oscillator) with respect to a reference source and show how this technique is useful for the detection of the phase and frequency offsets relative to a reference source.

A short review of the oscillator clock model and the Wigner-Ville distribution are introduced in Section 2. The proposed approach for the characterization of nonstationary oscillators is presented in Sections 3-6. In Section 7 numerical simulation results are presented including demonstration of phase and frequency deviations. Section 8 summarizes the paper.

2. Theoretical Background

In this section, we will introduce a short theoretical background of the clock model and the Wigner-Ville distribution. In Section 2.1 we will describe briefly the clock model and in Section 2.2 we will present briefly the Wigner-Ville distribution [33].

2.1. The Clock Model

A model for a real clock which includes variations in the amplitude and phase is given by [17, 26]where is the amplitude and is the frequency of the oscillation. and represent random fluctuations observed in amplitude and phase, respectively.

The amplitude fluctuations can be neglected [17, 26], and thus we can write the approximate model [17]:The phase variations result in an instantaneous frequency which is a time dependent process given by [17]In an ideal clock, the phase is constant and can be set to zero [51]; thus, it generates a sinusoidal oscillation of the form [52, 53]:where is the amplitude of the ideal clock.

Physical signals are real, in which case the Wigner-Ville distribution can be applied to either the (real) signal itself, or a complex version of the signal known as the analytic signal [54].

In the following, we define the analytic form of the clock signal asand for an ideal clock, the analytic signal will be

2.2. The Wigner-Ville Distribution

One of the most representative joint representations of time-frequency analysis is the Wigner distribution (WD), which is a quadratic signal representation introduced by Wigner (1932) [34] and later applied by Ville (1948) [35] to signal analysis. A comprehensive discussion of the WD properties can be found in a series of classical articles by Claasen and Mecklenbruker (1980) [55–57]. Originally, the WD has been applied to continuous variables as follows: consider an arbitrary 1D function , . The WD of is given by [33]where “” denotes complex conjugation. By considering the shifting parameter as a variable, (7) represents the Fourier transform (FT) of the product , where denotes the spatial frequency variable and, hence, the WD can be interpreted as the local spectrum of the signal .

The Wigner distribution satisfies many meaningful mathematical properties that are summarized in Appendix A.

The most powerful property is the “Time-Marginal” property and is given bywhere denotes the norm of the signal . Summing up of the energy distribution for all frequencies at a particular time would give the instantaneous energy of the signal. We will use the “Time-Marginal” property in our development in order to detect phase offset of a frequency source with respect to a reference source.

Modern signal processing uses the discrete form of physical signals. Therefore, the practical computation will use the discrete Wigner-Ville distribution form. There are several presentations in the literature for the discrete Wigner-Ville distribution. In Appendix B we present and explain in detail the discrete form of the Wigner-Ville Distribution and some of the main presentations. In this paper we will use the following discrete Wigner-Ville distribution, similar to the one proposed by Eric Chassande-Mottin and Archana Pai [58]:where , denotes the number of samples, and denotes the index numbers for the time and frequency vectors, respectively.

This equation, designed for discrete time or space signals, shows a great similarity with the original continuous version; therefore, it has been selected for our development.

3. A New Metric for Frequency Source Characterization

The system under consideration is presented in Figure 1.

At first we will derive a general expression for as a function of ; then we will develop simplified and closed expression for the Time Marginal of for three different cases.

Case A. In Section 4, we derive the analytic continuous and discrete forms for the Wigner-Ville distribution and its Time-Marginal expression for a constant phase offset case.

Case B. In Section 5, we derive the analytic expressions for the continuous Wigner-Ville distribution and its Time-Marginal expression for the case of constant phase and frequency offsets.

Case C. In Section 6, we derive the analytic expressions for the continuous Wigner-Ville distribution and its Time-Marginal expression for the case of frequency drift.
In this paper, we make use of an analytic signal which is the complex form of the real clock signal. Our choice is based on [54] because of the following reasons. All of the transformations in Cohens class of distributions (of which the Wigner-Ville is one) produce better results when applied to a modified version of the waveform termed the analytic signal, a complex version of the real signal. While the real signal can be used, the analytic signal has several advantages. The most important advantage is due to the fact that the analytic signal does not contain negative frequencies, so its use will reduce the number of cross products. If the real signal is used, then both the positive and negative spectral terms produce cross products.

Theorem 1. For the following assumptions:(1)without loss of generality, ,(2) is the nominal frequency,(3) denotes the continuous complex form of an ideal clock signal,(4) denotes the continuous complex form of the real clock signal.

The Wigner-Ville distribution of the difference is defined as

Proof of Theorem 1. The Wigner-Ville distribution (WD) of is given byBy using the “Summation Property” given in Appendix A (A.2) we obtainNote that for . Then, by using (12) we haveLet us get first a closed form expression for the second component in the right side of (13)In order to get a closed form expression for the first component on the right side of (13) we will use the next definition:where and .
By using the “Convolution Property” given in Appendix A we obtainwhere and are the WDs of and , respectively.
Note that the closed form expression for , where is given in (14) for .
By substituting (14) into (16) we obtainLet us define for simplicity the variable asThus, (17) can be written asNow we turn to the “cross-term” component of which is the third component on the right side of (13):Using (13), (17), and (20), the WD of (12) is given byThis completes our proof of Theorem 1.

4. Constant Phase Offset

As we all know, modern signal processing uses the discrete form of physical signals. Therefore, the practical computation will use the discrete Wigner-Ville distribution form as is introduced later in this subsection. Due to the simplicity of the mathematical development of using the continuous form of the WD upon its discrete one, we will use henceforth in our development process only the continuous form. In order to do so, we need to verify that there is a complete match between the analytic continuous form and the discrete one. Hence, in this subsection we will derive the analytic discrete form for the Wigner-Ville distribution and its Time-Marginal expression for a constant phase offset. In the end of the discrete form development, we will compare the discrete expression against the continuous expression introduced in the last subsection in order to verify a complete match and hence use only the continuous form.

4.1. Continuous Form

In this subsection, we will apply the analytic continuous form of the Wigner-Ville distribution on and then we will derive its Time-Marginal expression for a constant phase offset case.

Theorem 2. For the following assumptions:(1)without loss of generality, ,(2) is the constant phase offset from a reference clock signal,(3) is the nominal frequency,(4) denotes the continuous complex form of an ideal clock signal,(5) denotes the continuous complex form of the real clock signal for the case of constant phase offset.

The Time-Marginal expression of the Wigner-Ville distribution for is defined as

Proof of Theorem 2. For simplicity in the following we denote as .
The complex form of a clock signal with constant phase offset is given bywhere the phase variation considering the constant phase is defined byThe terms presented in (18) are given byBy substituting (25) into (17) we obtainBy substituting (25) into (20) we obtainThen, the real part of isBy substituting (14), (26), and (28) into (13) we obtainLet us derive a closed form expression for the Time Marginal of (29)This completes our proof of Theorem 2.

4.2. Discrete Form

Theorem 3. For the following assumptions:(1)without loss of generality, ,(2) is the constant phase offset from a reference clock signal,(3) is the nominal frequency,(4) denotes the discrete complex form of an ideal clock signal,(5) denotes the discrete complex form of the real clock signal for the case of constant phase offset.

The discrete Time-Marginal expression of the Wigner-Ville distribution for is defined as

Proof of Theorem 3. For discrete time and frequency variables the Wigner distribution is defined as [58]where , denotes the number of samples, and denotes the index numbers for the time and frequency vectors, respectively.
Let us derive a closed form expression for the discrete Wigner-Ville distribution (DWD) of :By using the “Summation Property” given in Appendix A (A.2) we obtainLet us get first a closed form expression for the second component on the right side of (34)Let us derive the first component on the right side of (34)Now we turn to the “cross-term” component of which is the third component on the right side of (34).
The complex form of a discrete clock signal with constant phase offset is given bywhere the constant phase is defined byThe term presented in (18) is given byThen, the cross-term derivation is given byBy substituting (39) into (40) we obtainThen, the real part of isBy substituting (35), (36), and (42) into (34) we obtainFinally, let us drive the closed form expression for the Time Marginal of (43)This completes our proof of Theorem 3.

As can be seen clearly, the discrete expression is the same as the continuous expression (22) introduced in the last subsection. This means that henceforth we can use only the continuous form and it will describe correctly the practical discrete signal processing.

5. Constant Frequency Offset

In this section, we will derive the analytic continuous form for the Wigner-Ville distribution and its Time-Marginal expression for the frequency offset case in addition to the phase offset case.

Theorem 4. For the following assumptions:(1)without loss of generality, ,(2) is the constant phase offset from a reference clock signal,(3) is the constant fractional frequency offset from a reference clock signal,(4) is the nominal frequency,(5) denotes the continuous complex form of an ideal clock signal,(6) denotes the continuous complex form of the real clock signal for the case of constant phase and frequency offsets.

The Time-Marginal expression of the Wigner-Ville distribution for is defined as

Proof of Theorem 4. For simplicity in the following we denote as .
The complex form of a clock signal with constant phase and frequency offsets is given bywhere the phase variation considering the constant phase and frequency offsets is defined bywhere we may defineBy substituting (48) into (18) we obtainIn order to carry out the first term in the right side of (12), the expression in (49) is substituted into (19):Next, in order to carry out the third term in the right side of (12), the expression in (48) is substituted into (20):Then, the real part of (51) isBy substituting (14), (50), and (52) into (13) we obtainNow we turn to derive the Time-Marginal expression of (53)This completes our proof of Theorem 4.

As we can see the result that we get is a periodic signal with a linear period of . In Section 7, we will examine the analytic expression versus the numerical simulation results in order to verify our expression.

6. Frequency Drift

Now we derive the analytic continuous form for the Wigner-Ville distribution and its Time-Marginal expression for the case of frequency drift along with constant phase and frequency offsets.

Theorem 5. For the following assumptions:(1)without loss of generality, ,(2) is the constant phase offset from a reference clock signal,(3) is the constant fractional frequency offset from a reference clock signal,(4) is the nominal frequency,(5) is the linear frequency drift coefficient (),(6) denotes the continuous complex form of an ideal clock signal,(7) denotes the continuous complex form of the real clock signal for the case of frequency drift.

The Time-Marginal expression of the Wigner-Ville distribution for is defined as

Proof of Theorem 5. For simplicity in the following we denote as .
The complex form of a clock signal with frequency drift along with constant phase and frequency offsets is given bywhere the phase is defined byBased on (57), can be written asNow, by substituting (58) into (18) we getNext, by substituting (59) into (17) we obtainSubstituting (58) into (20) we obtainwhere denotes the Fourier transform of .
In the following we derive the outcome of . For that purpose we defineThe Fourier transform of a product is defined by [59]where denotes convolution.
Based on (62), can be written asBased on (63), can be defined asThe Fourier transform of a complex Gaussian is given by [59]Based on (66) and (67) the expression for is given by .
Thus, based on (66) and (67) and that we haveBy substituting (65) and (68) into (64) we obtainWe use the following property [59]:Next, by substituting (69) into (70) we obtainAfter carrying out the integration in (71) we haveBy substituting (72) into (61) we obtainThen, the real part of isBy substituting (14), (60), and (74) into (13) we obtainBased on (8) and (75) we obtainThis completes our proof of Theorem 5.

7. Numerical Simulations

In this section, we will show the simulation results for the continuous and discrete Time-Marginal Wigner-Ville distribution of the differences [] and [] compared to our analytical proposed expressions (22) and (31), respectively. In order to demonstrate different nominal frequencies, we simulate the Time-Marginal Wigner-Ville distribution of the differences [] for three different nominal frequencies.

Figure 2 shows the contour plot of the Wigner-Ville distribution of a reference 10MHz signal , where the sampling frequency is 40MHz.

Figures 3–5 show for three different values of the comparison between the outcomes of the following expressions, (22) and (31), which were derived for the continuous and discrete cases, respectively.

According to Figures 3–5, there are high correlation between the outcomes of (22) and (31).

Figures 6 and 7 show the contour plot of (53), where the upper, lower, and the middle frequencies lines match the first, second, and the third terms of (53), respectively.

Figure 8 shows the contour plot of (75), for , , , and sampling frequency of 40MHz. We can see clearly that the linear chirp line follows the first term of (75). The constant frequency line matches the second term of (75) and the cross-term frequencies (in the middle) correspond to the third term of (75).

8. Summary and Conclusion

In this paper, we have presented a new tool for frequency source characterization based on the Wigner-Ville distribution, whose objective is to represent the variation of clock noise when no stationarity assumption is made. This means that the new tool can be used for a nonstationary phase variation. We presented the analytic terms for the Wigner-Ville distribution and its Time Marginal for constant phase and frequency offsets and frequency drift from a reference source. Simulation results indicate that our derivations hold.

Appendix

A.

In this appendix, we will present many meaningful mathematical properties of the Wigner distribution [33]. The mathematical properties can be summarized as follows:(P1) (the WD is a strictly real-valued function).(P2)If or then or .(P3)For real .(P4) (Time-Marginal property).(P5) ( represents the FT of the function ). ((P4) and (P5) are customarily referred as marginals).(P6)If , then (shift property).(P7)If , then (modulation property).(P8)If , then (convolution property).(P9) (inversion property).(P10)If [60, 61]: (summation property).Then, the Wigner-Ville distribution of is given bywhere and are the Wigner-Ville distributions of and , respectively, and refer to as “auto-terms”.

The cross-term is given by [60, 61]

The WD has an obvious interpretation as defining a local power spectrum at each point of time of the signal. Remarkably, property (P1) implies that the phase of the WD is implicit (phase remains hidden in the real-valued WD function).

Also, (P9) indicates that the WD is a reversible transform and the original signal can be recovered up to a constant factor (Claasen and Mecklenbrauker, 1980) [55–57].

B.

In this Appendix, we present and explain in detail the discrete form of the Wigner distribution. Although the WD was initially defined for continuous variable functions, Claasen and Mecklenbrauker (1980) [55] proposed a first definition for discrete variable functions. However, some attempts to extend definitions of the WD to discrete signals have not yet been completely successful (O’Neill, Flandrin, and Williams, 1998) [62]. For the applications in this paper, the following discrete Wigner distribution, similar to the one proposed by Eric Chassande-Mottin and Archana Pai [58], has been selected:where .

This equation, designed for discrete time or space signals, shows a great similarity with the original continuous version, but it also presents some differences and therefore has been referred to as a pseudo-Wigner distribution (PWD). Simply considered, the PWD involves using finite bounds of integration (i.e., a sliding analysis window). This originates a function which, relative to the true WD, is smoothed with respect to the frequency domain only. One important property preserved within the definition given by (B.1) is the inversion property (P9), which is a greatly desirable feature for the recovering of the original signal. Then, local filtering operations are also possible.

Equation (B.1) represents a PWD limited to a spatial interval . Here and represent the spatial and space-frequency discrete variables, respectively, and is the shifting parameter, which is also discrete. The values of represent the sampling positions in the signal considered as a digital replica of an original continuous representation. In this case, the numerical value of the position in the sequence is equal to the value of the analogical signal at position . Here is recognized as the sampling period and it is the inverse of the sampling frequency. Equation (B.1) originates an -component vector at each position . Furthermore, (B.1) must be interpreted as the discrete Fourier transform (DFT) of product in this interval or window.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

T. Krawinkel, S. Schn, and A. Bauch, “Application of miniaturized atomic clocks in kinematic GNSS single point positioning,” in Proceedings of the 28th European Frequency and Time Forum, EFTF 2014, pp. 97–100, Neuchatel, Switzerland, June 2014.
View at: Google Scholar
Y. Jin, B. Luo, X. Hu, and P. Zheng, “A fast positioning method based on high-precision clock for GNSS receivers,” in Proceedings of the 2015 IEEE International Conference on Information and Automation, ICIA 2015 - In conjunction with 2015 IEEE International Conference on Automation and Logistics, pp. 809–812, Lijiang, China, August 2015.
View at: Google Scholar
P. K. Bhatnagar, “Digital Network Synchronization: Basic Concepts,” in Engineering Networks for Synchronization, CCS 7, and ISDN:Standards, Protocols, Planning and Testing, pp. 6–41, Wiley-IEEE Press, 1997.
View at: Publisher Site | Google Scholar
G. L. Weaver, M. Miranian, and J. F. Garstecki, “Composite USO/CSAC timekeeping system,” in Proceedings of the 2012 IEEE Aerospace Conference, pp. 1–9, Big Sky, MT, USA, March 2012.
View at: Publisher Site | Google Scholar
L. A. Mallette, J. White, and P. Rochat, “Space qualified frequency sources (clocks) for current and future GNSS applications,” in Proceedings of the IEEE/ION Position, Location and Navigation Symposium, PLANS 2010, pp. 903–908, Indian Wells, CA, USA, May 2010.
View at: Publisher Site | Google Scholar
H. Hellwig, “Time and frequency applications,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 40, no. 5, pp. 538–543, 1993.
View at: Publisher Site | Google Scholar
INTERNATIONAL STANDARD, "Timing and Synchronization for Time-Sensitive Applications in Bridged Local Area Networks", ISO/IEC/IEEE 8802-1AS:2014(E), IEEE 2011.
R. Miskinis, D. Smirnov, E. Urba, and B. Dzindzelwta, “Timing and synchronization in mobile telecommunication networks,” in Proceedings of the 2011 Joint conference of the IEEE International Frequency Control and the European Frequency and Time Forum (FCS), pp. 1–5, San Francisco, CA, USA, May 2011.
View at: Publisher Site | Google Scholar
J. M. Cohen, H. E. Moses, and A. Rosenblum, “Accurate Time-Keeping and Accurate Navigation,” in Proceedings of the 1983 IEEE Military Communications Conference, pp. 729-730, Washington, DC, USA, October 1983.
View at: Publisher Site | Google Scholar
X. Guo, K. Liu, and J. Yang, “Adaptively precise time synchronization technique for inter-satellite link equipment of navigation satellite,” in Proceedings of the 2016 2nd IEEE International Conference on Computer and Communications (ICCC), pp. 1844–1848, Chengdu, China, October 2016.
View at: Publisher Site | Google Scholar
E. Re, A. Di Cintio, G. Busca, D. Giunta, and M. Sanchez, “Novel Time Synchronization techniques for Deep Space Probes,” in Proceedings of the 2009 Joint Meeting of the European Frequency and Time Forum (EFTF) and the IEEE International Frequency Control Symposium (FCS), pp. 205–210, Besancon, France, April 2009.
View at: Publisher Site | Google Scholar
F.-C. Chan, J. Mathieu, and P. Boris, “Stochastic modeling of atomic receiver clock for high integrity GPS navigation,” IEEE Transactions on Aerospace and Electronic Systems, vol. 50, no. 3, pp. 1749–1764, 2014.
View at: Publisher Site | Google Scholar
M. Levesque and D. Tipper, “A Survey of Clock Synchronization Over Packet-Switched Networks,” IEEE Communications Surveys & Tutorials, vol. 18, no. 4, pp. 2926–2947, 2016.
View at: Publisher Site | Google Scholar
S. Bregni, Synchronization of Digital Telecommunications Networks, John Wiley & Sons, Ltd, Chichester, UK, 2002.
View at: Publisher Site
K. V. Katsaros, W. Chai, N. Wang, G. Pavlou, H. Bontius, and M. Paolone, “Information-centric networking for machine-to-machine data delivery: A case study in smart grid applications,” IEEE Network, vol. 28, no. 3, pp. 58–64, 2014.
View at: Publisher Site | Google Scholar
A. Vallat and D. Schneuwly, “Clock synchronization in telecommunications via PTP (IEEE 1588),” in Proceedings of the IEEE International Frequency Control Symposium Joint with the 21st European Frequency and Time Forum (FCS '07), pp. 334–341, June 2007.
View at: Publisher Site | Google Scholar
L. Galleani and P. Tavella, “Detection of atomic clock frequency jumps with the Kalman filter,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 59, no. 3, pp. 504–509, 2012.
View at: Publisher Site | Google Scholar
L. Galleani and P. Tavella, “The characterization of clock behavior with the Dynamic Allan Variance,” in Proceedings of the IEEE International Frequency Control Sympposium and PDA Exhibition Jointly with the 17th European Frequency and Time Forum, vol. 2003, pp. 5–8, Tampa, Fl, USA, May 2003.
View at: Publisher Site | Google Scholar
L. Galleani and P. Tavella, “Tracking nonstationarities in clock noises using the dynamic allan variance,” in Proceedings of the 2005 IEEE International Frequency Control Symposium and Exposition, pp. 392–396, Vancouver, Canada, 2005.
View at: Publisher Site | Google Scholar
L. Galleani and P. Tavella, “The dynamic allan variance,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 56, no. 3, pp. 450–464, 2009.
View at: Publisher Site | Google Scholar
I. Sesia, L. Galleani, and P. Tavella, “Application of the dynamic allan variance for the characterization of space clock behavior,” IEEE Transactions on Aerospace and Electronic Systems, vol. 47, no. 2, pp. 884–895, 2011.
View at: Publisher Site | Google Scholar
L. Galleani, “Detection of changes in clock noise using the time-frequency spectrum,” Metrologia, vol. 45, no. 6, pp. S143–S153, 2008.
View at: Publisher Site | Google Scholar
L. Galleani and P. Tavella, “Detection and identification of atomic clock anomalies,” Metrologia, vol. 45, no. 6, pp. S127–S133, 2008.
View at: Publisher Site | Google Scholar
L. Galleani, “A Time-Frequency Relationship Between the Langevin Equation and the Harmonic Oscillator,” in Pseudo-Differential Operators: Groups, Geometry and Applications, pp. 119–131, 2017.
View at: Publisher Site | Google Scholar
M. G. Amin, D. Borio, Y. D. Zhang, and L. Galleani, “Time-Frequency Analysis for GNSSs: From interference mitigation to system monitoring,” IEEE Signal Processing Magazine, vol. 34, no. 5, pp. 85–95, 2017.
View at: Publisher Site | Google Scholar
Keysight Technologies, "Analyzing Frequency Stability in the Frequency and Time Domains", Application Note, Published in USA, August 11, 2014 5991-4797EN.
D. Allan and J. Barnes, “A Modified "Allan Variance" with Increased Oscillator Characterization Ability,” in Proceedings of the Thirty Fifth Annual Frequency Control Symposium, pp. 470–475, May 1981.
View at: Publisher Site | Google Scholar
IEEE Standard Definitions of Physical Quantities for Fundamental Frequency and Time Metrology, IEEE Standard 1139, 1999.
D. W. Allan, “Statistics of atomic frequency standards,” Proceedings of the IEEE, vol. 54, no. 2, pp. 221–230, 1966.
View at: Publisher Site | Google Scholar
J. A. Barnes, A. R. Chi, L. S. Cutler et al., “Characterization of Frequency Stability,” IEEE Transactions on Instrumentation and Measurement, vol. IM-20, no. 2, pp. 105–120, 1971.
View at: Publisher Site | Google Scholar
C. Greenhall, “Does Allan variance determine the spectrum?” in Proceedings of the International Frequency Control Symposium, pp. 358–365, Orlando, FL, USA, 1997.
View at: Publisher Site | Google Scholar
M. Matejcek and M. Sostronek, “Computation and evaluation allan variance results,” in Proceedings of the 2016 New Trends in Signal Processing (NTSP), pp. 1–9, Demanovska dolina, Slovakia, October 2016.
View at: Publisher Site | Google Scholar
L. Cohen, Time-Frequency Analysis, Prentice Hall, 1994.
E. Wigner, “On the quantum correction for thermodynamic equilibrium,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 40, no. 5, pp. 749–759, 1932.
View at: Publisher Site | Google Scholar
J. Ville, “Theorie et applications de la notion de signal analytique,” Cables et Transmission, vol. 2, no. 1, pp. 61–74, 1948.
View at: Google Scholar
W. Martin and P. Flandrin, “Wigner-Ville Spectral Analysis of Nonstationary Processes,” IEEE Transactions on Signal Processing, vol. 33, no. 6, pp. 1461–1470, 1985.
View at: Publisher Site | Google Scholar
R. Torres and E. Torres, “Fractional Fourier analysis of random signals and the notion of alpha-stationarity of the Wigner-Ville distribution,” IEEE Transactions on Signal Processing, vol. 61, no. 6, pp. 1555–1560, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
J. M. O'Toole, M. Mesbah, and B. Boashash, “Improved discrete definition of quadratic time-frequency distributions,” IEEE Transactions on Signal Processing, vol. 58, no. 2, pp. 906–911, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
X. Lv, G. Bi, C. Wan, and M. Xing, “Lvs distribution: principle, implementation, properties, and performance,” IEEE Transactions on Signal Processing, vol. 59, no. 8, pp. 3576–3591, 2011.
View at: Publisher Site | Google Scholar | MathSciNet
P. Flandrin and P. Borgnat, “Time-frequency energy distributions meet compressed sensing,” IEEE Transactions on Signal Processing, vol. 58, no. 6, pp. 2974–2982, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
T. O. Gulum, P. E. Pace, and R. Cristi, “Extraction of polyphase radar modulation parameters using a Wigner-Ville distribution-Radon transform,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '08), pp. 1505–1508, IEEE, Las Vegas, Nev, USA, April 2008.
View at: Publisher Site | Google Scholar
Y. Wu and D. C. Munson Jr., “Multistatic passive radar imaging using the Smoothed Pseudo Wigner-Ville Distribution,” in Proceedings of the IEEE International Conference on Image Processing (ICIP), pp. 604–607, October 2001.
View at: Google Scholar
S. Gabarda, G. Cristobal, J. Martínez-Alajarín, and R. Ruiz-Merino, “Detection of anomalous events in biomedical signals by Wigner analysis and instant-wise Rényi entropy,” in Proceedings of the 14th European Signal Processing Conference, EUSIPCO 2006, pp. 1–5, September 2006.
View at: Google Scholar
F. Akdeniz, I. Kayikiolu, and I. Kaya, “Using Wigner-Ville distribution in ECG arrhythmia detection for telemedicine applications,” in Proceedings of the 2016 39th International Conference on Telecommunications and Signal Processing (TSP), pp. 409–412, Vienna, Austria, June 2016.
View at: Publisher Site | Google Scholar
F. Akdeniz and T. Kayikcioglu, “Wigner-Ville distribution based ECG arrhythmia detection for telemedicine applications,” in Proceedings of the 2016 24th Signal Processing and Communication Application Conference (SIU), pp. 2045–2048, Zonguldak, Turkey, May 2016.
View at: Publisher Site | Google Scholar
S. C. Pei and J. J. Ding, “Fractional Fourier transform, Wigner distribution, and filter design for stationary and nonstationary random processes,” IEEE Transactions on Signal Processing, vol. 58, no. 8, pp. 4079–4092, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
G. F. Boudreaux-Bartels and T. W. Parks, “Time-varying filtering and signal estimation using Wigner distribution synthesis techniques,” Institute of Electrical and Electronics Engineers Transactions on Acoustics, Speech and Signal Processing, vol. ASSP-34, no. 3, pp. 442–451, 1986.
View at: Publisher Site | Google Scholar | MathSciNet
L. Yang, L. Zhao, S. Zhou, and G. Bi, “Sparsity-Driven SAR Imaging for Highly Maneuvering Ground Target by the Combination of Time-Frequency Analysis and Parametric Bayesian Learning,” IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, vol. 10, no. 4, pp. 1453–1465, 2017.
View at: Google Scholar
P. S. Bidari, J. Alirezaie, and J. Tavakkoli, “Shear wave elastography using Wigner-Ville distribution: A simulated multilayer media study,” in Proceedings of the 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBC 2016, pp. 2873–2876, USA, August 2016.
View at: Publisher Site | Google Scholar
D. Pang, P. Zhao, and B. Deng, “Cross-terms suppression in Wigner-Ville distribution based on image processing,” in Proceedings of the 2010 IEEE International Conference on Information and Automation, ICIA 2010, pp. 2168–2171, China, June 2010.
View at: Google Scholar
P. Kartaschoff and H. Brandenberger, Frequency and Time, Burlington, MA: Academic Press, New Jersey, NJ, USA, 1978.
View at: Publisher Site
D. W. Allan, N. Ashby, and C. C. Hodge, "The science of timekeeping" Hewlett Packard Application Note 1289, Palo Alto, Santa Clara, Calif, USA, 1997.
P. Tavella, “Statistical and mathematical tools for atomic clocks,” Metrologia, vol. 45, no. 6, pp. S183–S192, 2008.
View at: Publisher Site | Google Scholar
J. L. Semmlow, Biosignal and Biomedical Image Processing MATLAB based Applications, Robert Wood Johnson Medical School New Brunswick, New Jersey, NJ, USA, 2004.
T. A. C. M. Claasen and W. F. G. Mecklenbruker, “The Wigner distributiona tool for time-frequency signal analysis. I. Continuous-time signals,” Philips Journal of Research, vol. 35, no. 3, pp. 217–250, 1980.
View at: Google Scholar | MathSciNet
T. A. C. M. Claasen and W. F. Mecklenbruker, “The Wigner distributiona tool for time-frequency signal analysis Part II. Discrete-time signals,” Philips Journal of Research, vol. 35, no. 4/5, pp. 276–300, 1980.
View at: Google Scholar | MathSciNet
T. A. C. M. Claasen and W. F. G. Mecklenbruker, “The Wigner distribution A tool for time-frequency analysis, Parts I-III,” Philips Journal of Research, vol. 35, pp. 276–300, 1980.
View at: Google Scholar | MathSciNet
E. Chassande-Mottin and A. Pai, “Discrete time and frequency Wigner-Ville distribution: Moyal's formula and aliasing,” IEEE Signal Processing Letters, vol. 12, no. 7, pp. 508–511, 2005.
View at: Publisher Site | Google Scholar
http://www.thefouriertransform.com/pairs/complexGaussian.php.
M. Amirmazlaghani and H. Amindavar, “Modeling and Denoising Wigner-Ville Distribution,” in Proceedings of the 2009 IEEE 13th Digital Signal Processing Workshop and 5th IEEE Signal Processing Education Workshop, pp. 530–534, Amirkabir University of Technology, Tehran, Iran, January 2009.
View at: Publisher Site | Google Scholar
M. Sandsten, "Time-frequency analysis of non-stationary processes-an introduction", Centre for Mathematical Sciences, 2013.
P. Flandrin, R. Baraniuk, and O. Michel, “Time-frequency complexity and information,” in Proceedings of the ICASSP '94. IEEE International Conference on Acoustics, Speech and Signal Processing, pp. 329–332, Adelaide, SA, Australia, 1994.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2018 Chagai Levy et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1744

Downloads

820

Citations