Abstract

Oscillators and atomic clocks, as well as lasers and masers, 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 coherence of the oscillatory signal produced by the device. Such instabilities demonstrate frequency drifts, modulation, and spectrum broadening and are observed to be nonstationary processes in nature. Most of the tools which are being used to measure and characterize oscillator stability are based on signal processing techniques, assuming time invariance within a temporal window, during which the signal is assumed to be stationary. This letter proposes a new time-frequency approach for the characterization of frequency sources. Our technique is based on the Wigner–Ville time-frequency distribution, which extends the spectral measures to include the temporal nonstationary behavior of the processes affecting the coherence of the oscillator and the accuracy of the clock. We demonstrate the use of the technique in the characterization of nonstationary phase noise in 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, oscillators, and accurate frequency clocks with very high precision are required. Stable oscillators, including atomic ones such as rubidium and even cesium are often used in communication and navigation facilities [1, 2]. Apart from the above applications, many other electronic systems require frequency generators. In such systems, time synchronization and an accurate frequency are required [3, 4].

Oscillators and clocks, as well as lasers and masers, are subjected to environmental conditions, internal malfunctions, and inherent physical phenomena causing instabilities in the oscillatory signal being produced [5]. These instabilities are expressed by amplitude and phase fluctuations causing frequency deviations, drift, and spectrum broadening. Instead of being stationary coherent, the generated signal is observed to be a stochastic process characterized by statistical and spectral features, which are nonstationary time-varying [6, 7]. Therefore, their coherent features in the time and frequency domain should be investigated.

Phase noise is the frequency domain representation of rapid, short-term, and random fluctuations in the phase of a waveform. Historically, there have been two conflicting yet widely used definitions for phase noise. Some authors define phase noise to be the spectral density of a signal’s phase only [8], while the other definition refers to the phase spectrum (which pairs up with the amplitude spectrum) resulting from the spectral estimation of the signal itself. The IEEE defines phase noise as , where the “phase instability” is given in terms of the power spectral density of a signal’s phase deviation defined for [9]. However, the standard phase noise [9] does not give the possibility of identifying and interpreting nonstationary effects.

In lasers and masers, the radiation line-width limit is calculated via the Schawlow–Townes expression [1012], which assumes a stationary process [1317] although it is not stationary in general.

In this letter, we propose a time-frequency metric for the characterization of oscillators and frequency generators, also considering nonstationary variations. It is based on the Wigner–Ville distribution [18, 19], proposed in 1932 for the characterization of quantum fluctuations. This time-frequency distribution is applicable for the analysis of nonstationary signals [2023], radar signals [24, 25], biomedical signals [2628], analysis of time-varying filters [29, 30], and image processing [3133]. Up to now, the Wigner–Ville distribution was not yet used in the characterization of nonstationary oscillators such as evaluation of the phase noise. In this letter, we introduce an analytic expression of the Wigner–Ville time-frequency distribution for the characterization of the temporal coherency of nonstationary signals produced by classical and quantum oscillators. Our proposed expression can reveal the time-varying frequencies generally seen in the clock error noise under nonstationary conditions and in other platforms such as oscillator under vibrations or shocks.

The analysis is in general considering noise generated by different mechanisms. It is valid also for multinoise sources applied simultaneously to multinodes of the oscillator system. The multinoise sources are combination of two basic classes of noise, additive and parametric, in which their behavior and the underlying physical mechanisms are strikingly different. Each noise caused by an internal behavior of an oscillator or by a near-dc process modulates the carrier in phase.

2. The Oscillator Model

The model for the oscillator includes variation in the phase and given by [5, 34]where is the amplitude, is the central frequency of the oscillation, and represents the random fluctuation in phase. The phase variations is given as a summation of seven terms [35, 36]:where(1) and are the initial phase and frequency offsets, respectively.(2) is defined as the frequency drift coefficient.(3) is a linear variation coefficient added to the frequency drift.(4) are three independent, one-dimensional standard Wiener processes (Brownian motion) [35]. Following [37], the flicker noise is usually omitted from the mathematical model since it is not a rational process. The Wiener process is often referred to as an integral of white noise and is a nonstationary process.(5)The constants , , and represent the magnitude of the white frequency noise, the random walk frequency noise, and the frequency drift, respectively, which are three contributors to the noise in a typical oscillator. The constants , , and are determined by the intensity of each of these respective noise components, resulting from internal electronic noise mechanisms in an oscillator, such as thermal (Johnson) noise, shot noise, and parametric noise, in which a near-dc process modulates the phase of the carrier [38].(6)The frequency drift is a constant [35].

The stochastic term of the phase variations is given by the termwhile the deterministic term of the phase is expressed as shift and drift in the oscillator frequency. Hereafter, we concentrate in the stochastic term (3) only and analyse its time-frequency characteristics using the Wigner–Ville distribution.

The stochastic term (3) is the phase noise causing broadening in the line width.

3. The Wigner–Ville Distribution

In this section, we will present briefly the Wigner–Ville distribution and its spectrum.

One of the most representative joint representations of time-frequency analysis is the Wigner–Ville distribution (WD), which is a quadratic signal representation introduced by Wigner in 1932 [19] and later applied by Ville in 1948 [39] for signal time-frequency analysis.

In the following, we will denote the expectation value of as .

The statistical autocorrelation of is denoted by and defined as follows:where the shifting parameter is denoted by .

The Wigner spectrum is defined as the Fourier transform of the statistical autocorrelation [40]:where denotes the frequency variable.

4. Time-Frequency Analysis of Phase Noise Using Wigner Spectrum

This section proposes a new time-frequency closed form expression for the characterization of nonstationary phase noise in oscillators. Our technique is based on the Wigner spectrum, which extends the spectral measures to consist of the temporal nonstationary behavior of the phase process affecting the accuracy of the clock.

Theorem 1. We derive the Wigner spectrum (5) of the phase variation given in (3) and find that the closed form expression for a nonstationary phase noise is

Proof 1. By substituting (3) into (4), we obtainBecause are two independent, one-dimensional standard Wiener processes, we can write (7) as follows: is a real signal; therefore, , and thus we can write (8) as follows:The standard Wiener process in its general form is a continuous-time Gaussian process on , which starts at zero, has expectation zero for , and has the following covariance function [41, 42]:In our case, the process is the standard Brownian motion, i.e., . By substituting (10) into (9), we obtainNow, we will insert the expectation value into the joint integral argument on the second term of (11):Next, we will preform the expectation value of the second term of (12) by using (10):Finally, we may write (13) in its final form:Then, we will turn to preform the closed form expression of by substituting (14) into (5):whereNext, we will present the closed form expression for by solving the Fourier transform of (16):where denotes the Fourier transform of given in (16), respectively.
Finally, by solving the Fourier transform of (16) and substituting its closed form expressions into (17), we get (6).
This completes our Proof 1.
It is important to note that expression (6) for the Wigner time-frequency spectrum describes the temporal evolution of the phase-noise spectrum instantaneously. The power spectrum of the phase noise is obtained by integration of the resulted Wigner spectrum (6) within a temporal time “window” , an operation which is “smoothing” short-term nonstationary instabilities in the phase variations:A random process is said to be stationary (in a strict-sense) if all its statistical properties are invariant under an arbitrary time shift, that is, the autocorrelation function depends only on the time difference. Therefore, by using the Wiener–Khintchine theorem, the resulted power spectrum of a stationary process is dependent only on the frequency (without any temporal variations).
Please note that the integration of (6) in a finite observation time produced a frequency dependent power spectrum as expected in stochastic processes, which are wide-sense stationary.
Inspection of equation (18) reveals two different stationary noise processes;
The first term describes the white frequency noise which declines by power law while the second term describes the random walk frequency noise which declines by power law [35, 38].
An estimation of the nonstationary characteristics of (6) can be made by fitting the two noise sources expressed in (18) with the power law function for the phase noise spectra [38]. Please note that preforming the estimation gives and , which is exactly the results for and coefficients derived in [35].

5. Numerical Simulations

The Wigner–Ville distribution can be used to identify short-term, nonstationary instabilities in the oscillator phase and follow its spectral evolution in time. Differing from the power spectral density, which is a result of a time assimilation, assuming a stationary behavior within the temporal integration window, the Wigner–Ville distribution is a fundamental time-dependent measure, which can trend phase variations in ultrahigh resolution, depending on the sampling rate only. We employ the Wigner–Ville spectrum to analyse the nonstationary phase noise of a high accuracy Rubidium atomic clocks. Two clock models are used in the experiments: AR13300 [43] and NAC1 [44]. These atomic clocks are designed to generate a harmonic signal at frequency of .

While the NAC1 Rubidium resonator is locked on local Temperature-Compensated Crystal Oscillator (TCXO) and the AR13300 Rubidium resonator is locked on a local Oven-Controlled Crystal Oscillator (OCXO), and thus expected to be more stable (characterized with a lower diffusion parameter ). Table 1 summarizes the atomic clocks characteristics of both models.

The two noise types, white frequency noise and the random walk frequency noise, which are characterized by and , respectively, are those commonly encountered in the measurement of oscillator phase noise. The manufacturer predefined the total power spectrum and then and magnitudes extracted [38].

In Figures 1 and 2, we present a three-dimensional plot of the calculated phase-noise Wigner spectrum of the Rubidium-based TCXO (model NAC1) and the Rubidium-based OCXO (model AR13300) atomic oscillators, respectively. The temporal spectral variations of the OCXO based are revealed to be smaller than that of the NAC1 model. Integration of the Wigner spectrum in time, results in the power spectral density , as shown in Figure 3. A comparison is made between the NAC1 and the AR13300. The last is shown to be more stable as expected.

6. Summary and Conclusion

In this letter, we presented an effective technique for phase noise characterization of oscillators in time-frequency domain. While the common measure of phase noise power spectrum fails to follow nonstationary processes, the Wigner–Ville distribution is shown to demonstrate high-accuracy analysis of short time phase instabilities and a comprehensive description of their evolution in high levels of temporal resolution. It is shown that the Wigner–Ville distribution can be utilized for a nonstationary phase noise characterization that can reveal time-varying spectrum, generally seen in the oscillators and clocks, during fast variating conditions, including vibrations, shocks, and temperature variations. We derived an analytic expression for the Wigner–Ville spectrum for nonstationary phase noise and demonstrate its application on atomic clocks. The technique is applicable for the characterization of different oscillators, frequency generators, and radiation sources including lasers and masers.

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.