An Improved Time-Frequency Analysis Method in Interference Detection for GNSS Receivers
Abstract
: In this paper, an improved joint time-frequency (TF) analysis method based on a reassigned smoothed pseudo Wigner–Ville distribution (RSPWVD) has been proposed in interference detection for Global Navigation Satellite System (GNSS) receivers. In the RSPWVD, the two-dimensional low-pass filtering smoothing function is introduced to eliminate the cross-terms present in the quadratic TF distribution, and at the same time, the reassignment method is adopted to improve the TF concentration properties of the auto-terms of the signal components. This proposed interference detection method is evaluated by experiments on GPS L1 signals in the disturbing scenarios compared to the state-of-the-art interference detection approaches. The analysis results show that the proposed interference detection technique effectively overcomes the cross-terms problem and also preserves good TF localization properties, which has been proven to be effective and valid to enhance the interference detection performance of the GNSS receivers, particularly in the jamming environments.1. Introduction
The demand for safety-critical applications (e.g., civil aviation, aircraft landing) using GNSSs has gained extensive and increased interest in recent years. With the advent of the new GNSSs, such as the European Union's Galileo system [1], the USA's modernized GPS [2] and China's Beidou/Compass [3], new radio navigation signals will be broadcast, and more attention has been devoted to the design of the signal structure, which is expected to increase the accuracy, availability, integrity and continuity of service, especially in the field of the safety of life (SOL) applications (e.g., the accuracy needed during the landing of an aircraft). High precision positioning and reliable SOL services represent the main challenges for the upcoming satellite navigation systems. The presence of disturbing signals, such as spurious, harmonic and electromagnetic interferences, will result in serious performance degradation for GNSS receivers.
Radio frequency interference (RFI) is one of the biggest threats for satellite navigation systems. Although the satellite navigation system has a certain capability to be immune from interference, since the direct sequence spread spectrum (DSSS) is utilized, due to the low power of the received GNSS signals, the presence of intentional or unintentional disturbing signals, such as spurious, harmonic and electromagnetic interferences, will result in serious performance degradation for GNSS receivers. Among all of the different error sources that can potentially corrupt GNSS signals, RFI is particularly harmful, since, in some cases, it cannot be mitigated by a simple correlation process [4]. The jamming environment is threatening for satellite navigation systems. Many systems rely on the transmission of radio frequency (RF) energy in the L-band. For example, the European Galileo E5a and E5b radio bands, located within 1164–1214 MHz, occupy frequencies already allocated for aeronautical radio navigation services (ARNS), such as tactical air navigation (TACAN), distance measuring equipment (DME) and secondary surveillance radar (SSR) [5]. In addition, other RF transmissions, such as continuous wave (CW) signals originating from the European Digital Video Broadcast Terrestrial (DVB-T) service, can be considered as the main threat for the GNSS signals, since they appear in Radio Navigation Satellite System (RNSS) frequency bands. The presence of RFI and other channel impairments can heavily degrade the reception of useful GNSS signals, which results in poor navigation accuracy or complete loss of tracking for GNSS receivers.
Currently, interference detection and mitigation (ID&M) have become very important issues for GNSS applications. Interference can be detected and mitigated through various means. On the hardware side, specialized instrumentation, such as choke rings or active beam-forming antennas, can be used to suppress interference and improve the reception of line of sight (LOS) satellite signals. The main drawbacks of these antenna techniques are their requirements for particular hardware configurations and their computational complexity [6–8]. In the GNSS receiver, a special front-end architecture design adopts pulse blanker or automatic gain control (AGC) to reduce interference. For example, a specific AGC and analog-to-digital converter (ADC) design performing digital pulse blanking has been implemented in the GPS L5 receiver [9]. In addition, a nonuniform ADC controlled by a digital AGC can be adaptively adjusted based on the interference power strength to obtain an optimized conversion gain of the ADC (i.e., the SNR loss in ADC) and improved receiver performance [10,11]. However, these approaches require digital access to the feedback control of the AGC device, which is not a common output of generic analog front-ends used in the GNSS receivers.
In the literature, several interference mitigation techniques have been proposed and investigated, and each of them differs according to the operate domain (time, frequency or space). These techniques can be classified according to the specific processing domain. Considering time domain techniques, temporal filtering can be usually adopted. This method can be implemented at the digital intermediate frequency (IF) level after the ADC in the front-end of the GNSS receiver, which is effective only against narrow band RFI sources, because a wide band interference cannot be easily discriminated from the thermal noise by adopting temporal filtering. Frequency domain techniques are generally based on spectral estimation of the incoming signal, which is obtained by applying signal processing techniques, such as the discrete Fourier transform (DFT). These frequency domain techniques are typically performed by comparing the spectrum of the received signal with a theoretical threshold, which is usually determined according to a statistical model representing the received signal [12].
Recently, the research topic on transformed domain techniques (e.g., TF transform) has obtained increasing attention in the ID&M for GNSS receivers [13–16]. These techniques allow one to observe the received GNSS signals in a joint domain. In many cases, interferences may appear for a limited time and present a very variable behavior in frequency. In comparison to the GNSS signal, interferences are extremely different in terms of time and frequency characteristics. The presence of an interfering signal is limited to a region of the two-dimensional TF plane, and the adoption of the TF analysis is allowed to detect different types of disturbing signals. TF representations (TFRs) map a one-dimensional time signal into a two-dimensional function of time and frequency [17], which have found significant applications in non-stationary signal analysis. An interference mitigation technique based on the TFR approach has been described [18], where the TFR of the received GNSS signal is obtained by performing an orthogonal-like Gabor expansion on the samples at the output of the ADC of the GNSS receiver front-end. Another class of ID&M algorithm is reported aiming at obtaining a representation of the received signal in a different domain by making use of the time-scale transformation, which can be performed by means of wavelet transform [19].
Interfering signals are usually concentrated in a limited area of the TF plane, while noise is spread over the entire plane. The TF analysis approaches are very appealing as countermeasures in the detection of a large variety of interfering signals for GNSS receivers. There are different tools representing TF distributions, and the commonly-used TFRs in interference detection for GNSS applications include spectrogram, Wigner–Ville distribution (WVD) [13,16] and Choi–Williams transform (CWT) [16]. The spectrogram and WVD have been considered in the interference detection for GNSS applications [13]. The spectrogram approach presents the TF resolution trade-off problems according to the uncertainty principle, providing poor TF localization properties. In order to overcome the TF resolution trade-off problems of the spectrogram, WVD has been used in interference detection for GNSS receivers [13,16]. WVD is well known, since it provides nearly the best TF resolution among all of the TF distributions and also satisfies a large number of good properties [20], but it presents very severe cross-interfering terms without any physical meaning between true signal components (auto-terms), due to the interaction of different frequency components [16].
Most of the previous works on interference detection for GNSS receivers based on TF analysis are concerned with various TFR plots. Nevertheless, the quadratic TF distribution is usually a biased estimator for signal instantaneous frequency due to the presence of cross-term problems or makes a trade-off of temporal and frequency resolution due to the limitation of the uncertainty principle. In order to overcome or attenuate the cross-interfering terms present in the the quadratic TF distributions, the Choi–Williams distribution has been proposed to detect the sweep interference for GNSS receivers [16]. Appropriate image processing techniques can be also used for detection and parameter estimation of chirp signals by line detection in an image [21]. In addition, several kernel design methods have been proposed to mitigate the cross-term effect [22–24]. Unfortunately, these techniques need heavy computational complexity when they are applied in a real-time context. A reasonable method is to to introduce a window function in the time domain to reduce the undesired cross-term effects; therefore, the concept of the pseudo Wigner–Ville distribution (PWVD) is educed [20]. Consequently, the window function in PWVD can partially suppress the cross-terms to some ex0tent; the disadvantage of the filtering window operation is the degradation of the resolution, particularly in the frequency domain.
In PWVD, the time window operation is equal to frequency filtering in WVD, which can reduce the number of cross-interfering terms by suppressing the interferences between signal components sufficiently separated in time. In order to obtain a better readable result, the cross-terms between components in the frequency domain should be also minimized. Thus, an additional window function is added in order to perform a smoothing in time independently of the frequency smoothing. Therefore, the smoothed version of PWVD, namely the smoothed pseudo Wigner–Ville distribution (SPWVD), can be obtained [25]. The SPWVD is characterized by a separable kernel, which allows the time and frequency smoothing to be adjusted independently, which becomes one of the most versatile of Cohen's class TF distributions.
The smoothing windows can be adopted to reduce the cross-terms significantly; unfortunately, the SPWVD method also smears localized components, leading to less accurate localization of the signal auto-components in the TF plane compared to the WVD approach. Therefore, a reassignment method can be advantageously applied to improve TF localization properties in SPWVD [25,26]. In this way, the reassigned smoothed pseudo Wigner–Ville distribution (RSPWVD) can be obtained. The RSPWVD method is used to compensate for faults in mapping the TF energy distribution by relocating the value of the neighboring energy to the gravity center rather than the geometric center. In the RSPWVD method, by the adoption of the two-dimensional low-pass filtering smoothing function, the cross-term artifacts present in the quadratic TF distribution can be efficiently eliminated; meanwhile, by the use of the reassignment, the TF localization and aggregation properties of the auto-terms of the signal can be significantly improved.
In this paper, an improved TF analysis method by adopting RSPWVD has been proposed in interference detection for GNSS receivers. To the best of our knowledge, this interference detection technique based on joint TF analysis by adopting RSPWVD in the interference detection units for GNSS receivers is new. The performance of the proposed method has been deeply evaluated in comparison with the existing TF analysis approaches. Different localization properties and cross-term effects in the TF plane have been well investigated and compared among the aforementioned TF distributions adopted in interference detection for GNSS receivers.
In order to prove the effectiveness of the proposed TF analysis method by adopting RSPWVD in interference detection for GNSS receivers, an experiment is accomplished in the GPS L1 signal, which is characterized in additive white Gaussian noise (AWGN) corrupted by linearly-modulated sweep interference (chirp disturbance). The analysis results show that the proposed joint TF analysis by using RSPWVD eliminates the cross-terms significantly and preserves the high resolution of time and frequency in the TF plane at the same time. This developed improved TF analysis technique by adopting RSPWVD in interference detection makes the spectral characteristic of the interfering term sharply distinguishable among the received GNSS signal, which provides improved readability in the TF plane and enhanced detection performance for GNSS receivers with respect to the state-of-the-art TF analysis approaches.
2. Signal and System Model
The signal at the input of a GNSS receiver, in a noisy environment with RFI, can be written as:
Ai is the amplitude of the i-th useful signal;
τi is the code phase delay introduced by the transmission channel;
ci(t − τi) is the pseudo random noise (PRN) code sequence, which is assumed to take a value in the set {−1,1};
di(t − τi) is the bit stream of the navigation message, binary phase-shift keying (BPSK) modulated, including satellite data; and each binary unit is called a bit;
fd,i is the Doppler frequency shift affecting the i-th useful signal, and φRF,i is the initial carrier phase offset;
fRF is the carrier frequency, and it depends on the GNSS signal band under analysis; in the case of the GPS L1 signal, fRF = fL1 = 1575.42 MHz.
In general, the disturbing term ηRF(t) can be expressed as:
The interfering signal jRF(t) can assume different forms depending on the jammer that generates it. Several efforts have been devoted to the analysis and characterization of civilian GNSS jammers; despite significant differences, the transmitted jamming signal is usually frequency modulated with an almost constant amplitude. In this paper, the interference term jRF(t) is assumed to be in the class of sweep interference (linear chirp). Sweep interference is one of the main classes of the interfering signals, and its corresponding time-domain function represented by sinusoids can be written as follows:
In a linear chirp, the instantaneous frequency finst(t) of the interfering signal evolves linearly with time over the interval [fRF+∆f0, fRF+∆f1], where fRF is the GNSS signal center frequency. Therefore, the instantaneous frequency finst(t) can be expressed as:
The input signal yRF(t) defined in Equation (1) is received by the receiver antenna, down-converted and filtered by the receiver front-end. Then, the received signal before the analog to digital (A/D) conversion can be written as:
Considering the interference term j[t], the mean power of the sweep interference can be defined as:
The jammer-to-noise ratio (JNR) is defined as follows:
In order to avoid the cross-terms resulting from the interaction between the positive and negative frequency parts of the spectrum, the analytic representation of the received signal is adopted, provided as follows:
3. Time Frequency Transforms
The classical method for analyzing a signal with time-varying frequency content is to split the time-domain signal into many segments. The signal to be transformed is multiplied by a window function, which is nonzero for only a short period of time, and then, take the Fourier transform of each segment as the window slid along the time axis, resulting in a two-dimensional representation of the signal. This is known as the short-time Fourier transform (STFT) operation, which is the most widely-used method for analyzing non-stationary signals. Additionally, simply, in the continuous time case, it is defined as:
In order to understand the time properties at a particular frequency, the definition of the STFT can also be expressed in the frequency domain by manipulating Equation (12), obtaining the following result:
In the STFT analysis, one intends to achieve both high time and frequency resolution if possible. However, the resolution in the time domain is limited by the width of the window function h(t); similarly, the resolution in the frequency domain is limited by the width of the frequency window H(ω). As a result, the choice of a window to represent the signal by its spectrogram imposes a compromise between the conversation of temporal localization and that of frequency localization. This compromise is due to Heisenberg's uncertainty principle, which states that the window width in time and the window width in frequency are inversely proportional to each other. Therefore, choosing a small time window leads to good resolution in time and necessarily implies poor resolution in frequency; conversely, a long time window yields poor time resolution, but good frequency resolution. The length of the window plays a fundamental role in this TF compromise.
The squared magnitude of the STFT, denoted by SPEC(t, ω), is called spectrogram, which can be written as follows:
In order to avoid the TF resolution trade-off problem of the spectrogram, WVD is adopted in interference detection for GNSS receivers. In the WVD, a time-dependent instantaneous auto-correlation function is chosen as:
The WVD of ya(t) is then defined as the Fourier transform of this time-dependent instantaneous auto-correlation function [27], written as follows:
In addition, this class of bilinear (or quadratic) TF distributions can be most easily understood in terms of the ambiguity function. If the inverse Fourier transform of the instantaneous auto-correlation function R(t, τ) is taken with respect to t instead of τ, the ambiguity function can be obtained as follows:
The ambiguity function can be used to monitor the disturbing effect in the received GNSS signals.
The WVD has a number of desirable properties that make it a good indicator of how the energy of the signal can be viewed as a function of time and frequency. First, the WVD of any signal is always real. Second, it satisfies the time marginal condition:
That is, by summing the TF distribution over all frequencies, the instantaneous energy of the signal at a particular time instant can be obtained. Similarly, the WVD also satisfies the frequency marginal condition:
In this case, by summing the TF distribution over all time, the power spectrum of the signal at a particular frequency can be obtained.
Although WVD has many good properties and provides nearly the best resolution among all of the TF techniques, its main drawback comes from undesirable cross-term interference. The WVD is said to be bilinear, because the analyzed signal enters twice in its calculation. Consider the signal y(t) = y1(t) + y2(t), where y(t), y1(t) and y2(t) are analytic. Expanding the instantaneous auto-correlation function of y(t), we can obtain:
Thus, the WVD of the sum of two signals is not the sum of their corresponding WVDs, but also of their XWVDs. This means that the spectrum energy density of the sum of two signals does not reduce to the sum of the individual densities (unless the signals are spectrally disjoint). If y1(t) and y2(t) are mono-component signals, WVDy1(t,ω) and WVDy2(t,ω) are the auto-terms, while 2Re{WVDy1y2(t, ω)} is a cross-term.
As a result, if a signal contains more than one component, in the TF plane, its WVD suffers from spurious features containing cross-terms that occur halfway between each pair of auto-terms. The magnitude of these oscillatory cross-terms can be twice as large as the auto-terms, and they do not possess any physical meaning. As an example, the time-domain signal, which contains four Gaussian components, is shown in Figure 1. The WVD of this signal is provided in Figure 2a, and correspondingly, the contour of the computed WVD is presented in Figure 2b. From Figure 2, there exist four peaks, which respectively denote the corresponding auto-terms of the four Gaussian components in the TF plane; in addition, the WVD also presents six cross-terms, which occur between each pair of Gaussian signal components, and among them, two cross-terms overlap in the diagonal intersection point of the rectangle connected by the four vertices (i.e., four Gaussian components in the joint TF plane). These extra cross-terms have large oscillating amplitudes due to the interaction of the different signal components. The magnitudes of the oscillatory cross-terms can be twice as large as the auto-terms, but in the diagonal intersection position, the intensity of the corresponding cross-term can be four-times as large as the auto-terms, since two cross-terms overlap in this intersection point.
In Figure 2b, the serious cross-terms apparently present in the regions of the TF plane where we expect no energy at all, which make proper interpretation impossible. This is the main drawback of the WVD approach.
A feasible method to depress the effect of the cross-terms is to introduce windows in WVD; therefore, the concept of PWVD is educed. The WVD weighs equally all times of the future and past for a given time, but in the practical calculation of the distribution for a time instant t, we may concentrate the properties of the signal near the time of interest rather than from minus to plus infinity. Therefore, if we want to emphasize the signal around time t, a window function h(τ) can be multiplied with the instantaneous ambiguity function in order to define the PWVD, which is the WVD windowed in the time direction [20], written as follows:
Since the WVD is highly nonlocal, the window function is adopted in the PWVD to make it local. The overlapped window function contributes the cross-term suppression for multi-component signals. The operation of time-window overlay is equal to the frequency filtering in WVD; consequently, the interference between time-shifted signals is usually attenuated. An example of PWVD of the four Gaussian components signal is provided in Figure 3a, and correspondingly, the contour of the computed PWVD is presented in Figure 3b. From the results, it is easy to find that only two cross-terms remain in the TF plane, and the number of the extra cross-terms decreases from six to two, since the the cross-interfering term between each pair of time-shifted signal components is attenuated. In addition, from Figure 3b, it is clear to observe that the TF localization properties of the auto-terms of the analyzed signal are deteriorated in comparison with the WVD approach. In the PWVD, the cross-interfering terms have been partially mitigated by the adoption of the window function. The time window operation is equal to frequency filtering in WVD, which can reduce the number of cross-interfering terms by suppressing the interferences between signal components sufficiently separated in time. However, these advantages are achieved at the price of a blurring of the auto-terms of the signal and a loss of many desirable theoretical properties.
4. Interference Detection Based on Joint TF Analysis by Using RSPWVD
Although the WVD has many good properties and provides nearly the best resolution among all of the TF analysis techniques, due to its intrinsic quadratic nature, it suffers from cross-term interference when it is applied to multi-component signals. This drawback severely hinders the usefulness of the WVD for detecting RFI characteristics in the TF plane. In Section 3, the PWVD has been used to suppress the cross-terms for multi-component signals; however, many desirable properties of the WVD, such as marginals and instantaneous frequency, are annihilated, and the TF concentration property is also attenuated. In order to eliminate the cross-terms present in the quadratic TF distribution and to preserve good time and frequency resolution at the same time, a reassigned SPWVD method has been proposed in interference detection for GNSS receivers.
4.1. SPWVD
The unsatisfactory results obtained with the existing TF distributions justify the search for better tools; one way of achieving this is to start from the general form of quadratic representations. All of these existing TF distributions could be written in a generalized form, which can be used to facilitate the design of desirable TF transforms. This class of transform is known as Cohen's class [20], and the definition of Cohen's class of bilinear (or quadratic) TF distributions can be written as follows:
Cohen's class can be also rewritten as the double convolution of the WVD of the signal ya(t) and a two-dimensional smoothing function, provided as follows:
If a separable smoothing function is considered, it can be written as the product of windows from both time and frequency domains:
Since the cross-terms with the WVD are strongly oscillating, the most effective way of removing cross-term interference is to apply two-dimensional low-pass filtering in the ambiguity domain. The resulting two-dimensional convolution of the WVD in Equation (25) defines the smoothed version of PWVD, that is SPWVD, which can be written as follows:
Therefore, an independent and progressive control can be applied to the WVD in both time and frequency directions. The independency of h(t) and g(t) makes SPWVD more flexible to reduce the cross-terms present in WVD.
4.2. Reassignment Method
Compared with WVD, SPWVD can be used to effectively depress the influence of the cross-interfering terms of a multi-component signal, but its TF concentration and localization properties decrease somewhat. In order to improve the TF aggregation properties in SPWVD, a reassignment method has been considered.
From Equation (25), we can know that the two-dimensional smoothing function Π(t − s, ω − θ) determines a certain TF region at the neighborhood nearby the point (t, ω), inside which a weighted average of the WVD Wya(s, θ) of the signal ya(t) is performed. However, these mean values may not be symmetrically distributed around a certain point (t, ω), which is the geometrical center of the TF region. Consequently, the point (t, ω) is not truly representative for such a region. In contrast, the energy gravity center of such a region is more approximate to represent the local energy distribution of the analyzed signal. The local energy distribution Π(t − s, ω − θ)C(s, θ; Π) of Cohen's class TF distribution can be assumed as the distribution of mass, and it is better to assign the total mass to the gravity center rather than to its geometrical center.
In this paper, the reassignment method is considered to relocate each value of Cohen's class distribution C(t, ω) at any point (t, ω) to another point (t̂, ω̂), which is the gravity center of the signal's energy distribution around the point (t, ω). Then, the reassigned Cohen's class TF distribution can be defined as follows:
The RSPWVD RSPWVD(t′, ω′; g, h) can be used to eliminate the cross-term artifacts inherent in the quadratic TF distributions, which presents good resolution in both the time and frequency domains; therefore, in this work, the joint TF analysis based on RSPWVD has been firstly proposed in interference detection for GNSS receivers, which is clearly illustrated in Figure 4. The double convolution of the WVD Wya(t, ω) of the analytical signal and the two-dimensional smoothing function Π(t, ω) is performed to obtain the SPWVD(t, ω; g, h); after reassignment, the cross-term-free TF distribution RSPWVD(t′, ω′; g, h) has good readability, which is suitable for multi-component signal analysis. The proposed joint TF analysis by using RSPWVD makes the spectral characteristic of the interfering signal sharply distinguishable among the received GNSS signal; therefore, the instantaneous frequency of the disturbing term in the received GNSS signals can be effectively estimated by detecting the peaks of the proposed joint TF distribution.
5. Performance Evaluation
In this section, the performance of the proposed interference detection method based on joint TF analysis by using RSPWVD is analyzed. In particular, this proposed new algorithm, which adopts RSPWVD to detect sweep interference in GNSS receivers, is compared with the conventional interference detection approaches in the disturbing scenario.
The mentioned interference detection approaches are tested on real GPS data collected by using the GNSS software receiver developed at Beihang University [28]. The scheme of the test is reported in Figure 5, while an image of the experimental setup adopted for collecting the GPS data corrupted by sweep interference is depicted in Figure 6. The real GPS samples are collected by using the GNSS software receiver connected to the Trimble Zephyr Geodetic 2 antenna placed on the roof of the new main building at Beihang University in an open-sky static condition; the developed software interferer is adopted for generating the sweep interfering signal for GNSS applications. The generated sweep interference is added to the GPS samples collected by the GNSS receiver front-end.
In the experiment, the scenario adopted for the test is characterized by the setting parameters provided in Table 1, representing the real GPS L1 signal in zero mean Gaussian noise corrupted by a constant amplitude linearly frequency modulated interference (linear chirp), which is commonly considered as a test bench in interference detection for GNSS receivers.
In the experiment settings, the JNR value of the sweep interference is set to be −1 dB; in a linear chirp, the instantaneous frequency finst of the interfering signal evolves linearly with time over the interval [fL1 + ∆f0, fL1 + ∆f1], where fL1 is the GPS L1 signal center frequency, ∆f0 = +7 MHz and ∆f1 = −9 MHz.
In Figure 7a, the ambiguity function of the GPS L1 signal without interference is shown as a spike. In the case with the presence of sweep interference, the ambiguity function of the interfered GPS L1 signal is depicted in Figure 7b, where the disturbing effect can be clearly observed. The adoption of the ambiguity function of the received interfering signal is beneficial to better monitor the interference contribution in the received GNSS signals.
In Figure 8, the spectrogram of the GPS L1 signal with sweep interference is depicted, where the Hamming window function is chosen. In Figure 8a, the window size is 63 samples, where the disturbing term more or less emerges in the TF plane, but very poor TF localization properties are obtained with this approach. In order to evaluate the TF characteristics of the spectrogram by increasing the length of the analysis window, the spectrogram of the interfered GPS L1 signal is presented in Figure 8b, where the window length is increased to 127 samples. It is clear that a long time window leads to improved frequency resolution and inevitably yields poor resolution in time. The spectrogram is nonlinear, but this nonlinearity results from the operation of squared magnitude and, therefore, does not lead to undesirable cross-terms present in the WVD. In practice, the spectrogram approach cannot be used to provide the instantaneous frequency estimate for the interfering signal, although sometimes, a good approximation can be achieved.
In Figure 9a, the WVD of the GPS L1 signal with the sweep interference is depicted, and correspondingly, the contour of the computed WVD is presented in Figure 9b. The interfering term presents a linear behavior in frequency, which is well localized in the restricted portion of the TF plane. The price that is paid for high TF resolution with the WVD approach is the undesirable cross-terms. Due to the interaction of the different signal components, the presence of the serious cross-terms is apparently observed in Figure 9b, which makes interpretation quite difficult and brings much error for the estimation of the instantaneous frequency of the interfering signal.
In Figure 10a, the PWVD of the GPS L1 signal with the sweep interference is depicted, and correspondingly, the contour of the computed PWVD is presented in Figure 10b. It is clear that although the sweep interfering term presents a linear behavior in the TF plane, the PWVD presents less accurate TF localization precision, since the TF concentration property of the PWVD is attenuated. In addition, the cross-terms can be still observed in the TF plane by using the PWVD approach.
In order to mitigate the cross-interfering terms present in the WVD, a two-dimensional smoothing function is used to smooth out the interfering signal in both the time and frequency domains. In Figure 11a, the SPWVD of the interfered GPS L1 signal is provided, where the peaks of the SPWVD are clearly observed, denoting the sweep interference contribution. The corresponding contour of the SPWVD is provided in Figure 11b, and it is easy to find that the undesired sweep interference shows a linear behavior in frequency, which is localized in the linear portion of the TF plane. From Figure 11b, it is very clear that the undesired cross-interfering terms are partially mitigated by the adoption of the two-dimensional smoothing function in the SPWVD. Although the adoption of the smoothing window in SPWVD is beneficial to suppress the cross-terms in the TF plane, it smears localized components, leading to less accurate TF localization of the auto-terms of the signal, as compared to the WVD approach. The disadvantage of this filtering window operation is that it limits the original excellent TF resolution features.
Therefore, in order to remove the cross-terms present in the quadratic TF distribution and to preserve the good TF energy concentration properties at the same time, an improved TF distribution by adopting RSPWVD has been proposed in interference detection for GNSS receivers. This proposed joint TF distribution of the interfered GPS L1 signal is depicted in Figure 12a, and the corresponding contour of the calculated TF distribution of RSPWVD is provided in Figure 12b. From the results, it is clear that the undesired sweep interfering term presents a very strict linear behavior in the TF plane with the proposed joint TF distribution based on RSPWVD, which provides representations that are easy to interpret in interference detection for GNSS receivers. This proposed improved TF analysis technique effectively eliminates the cross-term artifacts present in the quadratic TF distributions, which shows good readability in the TF plane. In the reassigned method, the smoothing window is moved from the geometrical center to the energy gravity center of the TF distribution. Therefore, with the proposed RSPWVD method, the disturbing interference is highly localized in the restricted region of the TF plane, and a much improved TF localization property can be achieved in comparison to the SPWVD approach. The proposed joint TF distribution based on RSPWVD has been proven to be valid and effective in interference detection for GNSS receivers in jamming environments.
This proposed joint TF analysis technique based on RSPWVD effectively solves the cross-term problem of the bilinear TF distribution and keeps good time and frequency resolution in the TF plane at the same time, which has been proven to be effective to adopt in interference detection for GNSS receivers. With this developed interference detection technique, the interfering signal can be correctly identified, and its instantaneous frequency can be precisely estimated. In the following interference excision/mitigation unit (anti-jamming device), the estimated instantaneous frequency of the interfering signal will be further used to control the coefficients of an excision filter (notch filter) that adaptively removes the disturbing interference signal.
6. Conclusions
In this paper, an improved TF analysis method based on RSPWVD has been proposed to detect sweep interference for GNSS receivers. In order to prove the advantages and effectiveness of the developed technique, a comprehensive performance comparison has been carried out compared with the existed TF analysis approaches. The experiments have been performed on GPS L1 signals in the disturbing scenario in order to support the theoretical analysis among the aforementioned TF distributions adopted in interference detection for GNSS receivers.
From the analysis results, the spectrogram approach presents the TF resolution trade-off problems and provides poor localization properties in the TF plane; the interference detection based on WVD presents a severe cross-interfering term problem due to the interaction of different frequency components, which makes proper disturbance interpretation impossible. Then, the SPWVD method is adopted to partially mitigate the cross-terms present in the bilinear TF distribution, since a two-dimensional smoothing function is added in both the time and frequency domains, but the disadvantage of using filtering windows is the degradation of time and frequency resolution.
In order to eliminate the cross-term problem and preserve good time and frequency resolution in the TF plane at the same time, an improved joint TF distribution by using RSPWVD has been firstly proposed in interference detection for GNSS receivers. The interference detection based on joint TF distribution by using RSPWVD efficiently combines the removal of the cross-interfering terms provided by a two-dimensional smoothing kernel function and an increased TF concentration of the auto-terms of the signal components achieved by the reassignment. From the analysis results, the proposed interference detection method based on RSPWVD successfully overcomes the cross-term problem and, meanwhile, presents good TF localization properties, which provide much improved interference detection performance in comparison with the existing TF analysis approaches.
The developed interference detection technique based on joint TF analysis by adopting RSPWVD has been proven to be suitable to improve the interference detection performance in the jamming environments, which is promising for adoption in the anti-interference GNSS receiver design for civil aviation and military purposes, particularly in disturbing environments.
Acknowledgements
This research work is supported by the Scientific Research Foundation for Returned Scholars, Ministry of Education of China.
Author Contributions
Kewen Sun conceived the work, designed the algorithms, performed the experiment, analyzed the experiment results and wrote the manuscript. Tian Jin helped during the phase of GPS data collection. Dongkai Yang commented on the work.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Galileo OS SIS ICD (Open Service Signal-In-Space Interface Control); Issue 1.1; European Union: September; 2010.
- Systems Enginneing & Integration Interface Specification IS-GPS-705D, Navstar GPS Space Segment/User Segment L5 Interface; Global Positioning Systems Directorate: Washington, DC, USA; 24; September; 2013.
- BeiDou Navigation Satellite System Signal in Space Interface Control Document Open Service Signal (Version 2.0); BDS-SIS-ICD-2.0; China Satellite Navigation Office: Beijing, China; December; 2013.
- Kaplan, E.; Hegarty, C. Understanding GPS: Principles and Applications, 2nd ed.; Artech House Press: Norwood, MA, USA, 2006. [Google Scholar]
- Bastide, F.; Chatre, E.; Macabiau, C.; Roturier, B. GPS L5 and Galileo E5a/E5b signal-to-noise density ratio degradation due to DME/TACAN signals: Simulations and theoretical derivation. Proceedings of the 2004 National Technical Meeting of The Institute of Navigation, San Diego, CA, USA, 26–28 January 2004.
- Zhang, Y.D.; Amin, M.G. Anti-Jamming GPS Receiver With Reduced Phase Distortions. IEEE Signal Process. Lett. 2012, 19, 635–638. [Google Scholar]
- Zhang, Y.D.; Amin, M.G. Array Processing for Nonstationary Interference Suppression in DS/SS Communications Using Subspace Projection Techniques. IEEE Trans. Signal Process. 2001, 49, 3005–3014. [Google Scholar]
- De Lamare, R.C.; Sampaio-Neto, R. Reduced-Rank Space-Time Adaptive Interference Suppression With Joint Iterative Least Squares Algorithms for Spread-Spectrum Systems. IEEE Trans. Veh. Technol. 2010, 59, 1217–1228. [Google Scholar]
- Bastide, F.; Akos, D.; Macabiau, C.; Roturier, B. Automatic Gain Control (AGC) as an Interference Assessment Tool. Proceedings of the 16th International Technical Meeting of the Satellite Division of The Institute of Navigation, Portland, OR, USA, 9–12 September, 2003.
- Amoroso, F.; Bricker, J.L. Performance of the Adaptive A/D Converter in Combined CW and Gaussian Interference. IEEE Trans. Commun. 1986, COM-34, 209–213. [Google Scholar]
- Amoroso, F. Adaptive A/D Converter to Suppress CW Interference in DSPN Spread-Spectrum Communications. IEEE Trans. Commun. 1983, 31, 1117–1123. [Google Scholar]
- Kuriger, G.; Grant, H.; Cartwright, A.; Heirman, D. Investigation of Spurious Emissions From Cellular Phones and the Possible Effect on Aircraft Navigation Equipment. IEEE Trans. Electromagn. Compat. 2003, 45, 281–292. [Google Scholar]
- Borio, D.; Camoriano, L.; Savasta, S.; Lo Presti, L. Time-Frequency Excision for GNSS Application. IEEE Syst. J. 2008, 2, 27–37. [Google Scholar]
- Ouyang, X.; Amin, M.G. Short-Time Fourier Transform Receiver for Nonstationary Interference Excision in Direct Sequence Spread Spectrum Communications. IEEE Trans. Signal Process. 2001, 49, 851–863. [Google Scholar]
- Lach, S.R.; Amin, M.G.; Lindsey, A.R. Broadband Interference Excision for Software-Radio Spread-Spectrum Communications Using Time-Frequency Distribution Synthesis. IEEE J. Sel. Areas Commun. 1999, 17, 704–714. [Google Scholar]
- Sun, K.; Liu, W.; Xu, H.; Yang, D. Interference Detection Based on Time-Frequency Analysis for GNSS Receivers. Proceedings of the 26th International Technical Meeting of the Satellite Division of the Institute of Navigation, ION GNSS + 2013 Conference, Nashville, TN, USA, 16–20 September 2013.
- Hlawatsch, F.; Boudreaux-Bartels, G.F. Linear and Quadratic Time-Frequency Signal Representations. IEEE Signal Process. Mag. 1992, 9, 21–67. [Google Scholar]
- Savasta, S.; Lo Presti, L.; Rao, M. Interference Mitigation in GNSS Receivers by a Time-Frequency Approach. IEEE Trans. Aerosp. Electron. Syst. 1992, 49, 415–438. [Google Scholar]
- Paonni, M.; Jang, J.; Eissfeller, B.; Wallner, S.; Avila-Rodriguez, J.A.; Samson, J.; Fernandez, F.A. Wavelets and Notch Filtering: Innovative Digital Processing Techniques for Mitigating RF Interference, Inside; GNSS; Gibbons Media & Research LLC: Eugene, OR, USA, 2011; pp. 54–62. [Google Scholar]
- Cohen, L. Time-Frequency Analysis: Theory and Applications; Prentice-Hall Press: Upper Saddle River, NJ, USA, 1995. [Google Scholar]
- Barbarossa, S. Analysis of Multicomponent LFM Signals by a Combined Wigner-Hough Transform. IEEE Trans. Signal Process. 1995, 43, 1511–1515. [Google Scholar]
- Ristic, B.; Boashash, B. Kernel Design for Time-Frequency Signal Analysis using the Radon Transform. IEEE Trans. Signal Process. 1993, 41, 1996–2008. [Google Scholar]
- Stankovic, L. Auto-term Representation by the Reduced Interference Distributions: A Procedure for Kernel Design. IEEE Trans. Signal Process. 1996, 44, 1557–1563. [Google Scholar]
- Jones, D.L.; Baraniuk, R.G. An Adaptive Optimal-Kernel Time-Frequency Representation. IEEE Trans. Signal Process. 1995, 43, 2361–2371. [Google Scholar]
- Auger, F.; Flandrin, P. Improving the Readability of Time-Frequency and Time-Scale Representations by the Reassignment Method. IEEE Trans. Signal Process. 1995, 43, 1068–1089. [Google Scholar]
- Richard, C.; Lengelle, R. Joint Recursive Implementation of Time-Frequency Representations and Their Modified Version by the Reassignment Method. Signal Process. 1997, 60, 163–179. [Google Scholar]
- Grochenig, K. Foundations of Time-Frequency Analysis; Birkhauser Press: Boston, MA, USA, 2001. [Google Scholar]
- Jin, T.; Ren, J. Stability Analysis of GPS Carrier Tracking Loops by Phase Margin Approach. GPS Solut. 2013, 17, 423–431. [Google Scholar]
Parameter | Value |
---|---|
Carrier-to-noise ratio, C/N0 | 46 dB-Hz |
Sampling frequency, fs | 24 MHz |
Intermediate frequency, fIF | 40.42 MHz |
Code length | 1023 chips |
Sweep period | 1 ms |
Spectrogram analysis window | Hamming |
© 2015 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license ( http://creativecommons.org/licenses/by/4.0/.
Share and Cite
Sun, K.; Jin, T.; Yang, D. An Improved Time-Frequency Analysis Method in Interference Detection for GNSS Receivers. Sensors 2015, 15, 9404-9426. https://doi.org/10.3390/s150409404
Sun K, Jin T, Yang D. An Improved Time-Frequency Analysis Method in Interference Detection for GNSS Receivers. Sensors. 2015; 15(4):9404-9426. https://doi.org/10.3390/s150409404
Chicago/Turabian StyleSun, Kewen, Tian Jin, and Dongkai Yang. 2015. "An Improved Time-Frequency Analysis Method in Interference Detection for GNSS Receivers" Sensors 15, no. 4: 9404-9426. https://doi.org/10.3390/s150409404
APA StyleSun, K., Jin, T., & Yang, D. (2015). An Improved Time-Frequency Analysis Method in Interference Detection for GNSS Receivers. Sensors, 15(4), 9404-9426. https://doi.org/10.3390/s150409404