Next Article in Journal
Review and Evaluation of MAC Protocols for Satellite IoT Systems Using Nanosatellites
Previous Article in Journal
A Novel Four Single-Sideband M-QAM Modulation Scheme Using a Shadow Equalizer for MIMO System Toward 5G Communications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Joint Method Based on Time-Frequency Distribution to Detect Time-Varying Interferences for GNSS Receivers with a Single Antenna

School of Electronic and Information Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(8), 1946; https://doi.org/10.3390/s19081946
Submission received: 14 March 2019 / Revised: 10 April 2019 / Accepted: 22 April 2019 / Published: 25 April 2019
(This article belongs to the Section Remote Sensors)

Abstract

:
In this paper, a joint method combining Hough transform and reassigned smoothed pseudo Wigner-Ville distribution (RSPWVD) is presented to detect time-varying interferences with crossed frequency for a Global Navigation Satellite System (GNSS) receiver with a single antenna. The proposed method can prevent the cross-term interference and detect the time-varying interferences with crossed frequency which cannot be achieved by the classical time-frequency (TF) analysis with the peak detection method. The actual performance of the developed method has been evaluated by experiments with conditions where the real BeiDou system (BDS) B1I signals are corrupted by the simulated chirp interferences. The results of experiments show that the introduced method is effectively able to detect chirp interferences with crossed frequency and provide the same root mean square errors (RMSE) of the parameter estimation for chirp one and the improved initial frequency estimation for chirp two compared with the Hough transform of Wigner-Ville distribution (WVD) when the jamming to noise ratio (JNR) equals or surpasses 4 dB.

1. Introduction

At present, GNSS receivers with a single antenna are threatened by a serious jamming environment where many GNSS receiver failures occurred [1,2,3], limiting the GNSS applications. As a result, the techniques used to detect and mitigate interference effects have become an increasingly important issue and can be divided into the automatic gain control (AGC) method [4], digital signal processing methods [5,6,7,8,9] and receiver methods [10,11] from GNSS signal processing chains.
The AGC acting as an adaptive variable gain amplifier adjusts the input signal level to the analog-to-digital converter (ADC) input range [10,12]. Therefore, the AGC gain variation can be utilized to detect interference, especially continuous wave interference (CWI) and pulsed interference (PI). However, there are not enough effective quantized bits for ADC to realize the gain variation range of AGC for a common GNSS receiver.
From the GNSS receiver, the position accuracy and effective carrier to noise density ratio are commonly adopted to assess the impact of interference on GNSS receivers [13,14]. However, both of them rely on a particular GNSS receiver performance. If the jamming power surpasses the spread spectrum gain of GNSS causing the GNSS receiver not to work, the corresponding algorithm will fail.
For the digital signal processing methods, they include spatial domain methods [5], spatial-temporal domain methods [15,16,17], time domain methods [6], frequency domain methods [7], and time-frequency (TF) domain methods [8,9]. Time domain techniques as well as frequency domain techniques cannot completely describe the nature of time-varying signals.
Spatial domain techniques [5] and spatial-temporal domain methods assume that the GNSS signal angle of arrival (AOA) and the antenna model are known for the receiver. The antenna model is sensitive to AOA estimation error and can cause a high computational complexity. In addition, an antenna array is required, which is difficult to install on small devices.
TF domain methods adopt the classical TF analysis including short-time Fourier transform (STFT), Wigner-Ville distribution (WVD) and RSPWVD [8], which illustrate that jamming signals usually focus their energy in a finite area of the TF plane, while noise power extends over the whole TF plane [9]. However, these TF methods cannot deal with the time-varying interferences with crossed frequency by the peak detection method [18]. In this paper, a joint method combining Hough transform and RSPWVD is depicted to detect chirp interferences with crossed frequency for GNSS receivers with a single antenna. The analytic expression of initial frequency estimation and chirp rate estimation is presented and the double threshold detection is proposed as well. In addition, the effect of sweep period on the estimation of initial frequency and chirp rate has been analyzed.

2. Signal and System Model

The signals received from an antenna of a GNSS receiver pass through the radio frequency (RF) front-end and are down-converted to intermediate frequency (IF). Before the ADC, they can be written as
X I F ( t ) = l = 1 N s S I F , l ( t ) + J I F ( t ) + N ( t )
where N ( t ) is the additive white Gaussian noise term with two-sided power spectral density N 0 / 2 and zero mean, J I F ( t ) is the jamming signal. N s is the number of visible satellites, S I F , l ( t ) is the signal transmitted by the lth satellite, which can be defined
S I F , l ( t ) = 2 A l C l ( t τ l ) D l ( t τ l ) cos [ 2 π ( f I F + f d , l ) t + φ I F , l ]
where,
  • A l is the received GNSS signal power from the lth satellite;
  • C l ( t τ l ) is the pseudorandom noise sequence, and τ l is the code phase delay;
  • D l ( t τ l ) is the navigation data message signal;
  • f I F is GNSS signal intermediate frequency;
  • f d , l is the Doppler-affected frequency;
  • φ I F , l is the carrier phase delay.
The jamming signal J I F ( t ) defined in Equation (1) is redefined
J I F ( t ) = m = 1 L 2 A J , m cos [ 2 π ( f i n s , m ( t ) ) t + θ J , m ]
  • A J , m is the power of the mth jamming signal;
  • f i n s , m ( t ) is the mth jamming instantaneous frequency;
  • θ J , m is the phase delay of the mth jamming signal;
  • L is the number of interferences.
For the linear chirp signal, the instantaneous frequency f i n s , m ( t ) linearly evolves with the time interval [ f b e g i n , f e n d ] . Therefore, the f i n s , m ( t ) is written as follows
f i n s , m ( t ) = f b e g i n + k t ,   0 t T j
where T j is the sweep period, f b e g i n is the initial frequency, and k is the frequency change rate, also named chirp rate, defined as follows
k = f e n d f b e g i n T j = B s w e e p T j
where B s w e e p represents the sweep bandwidth. The mth jamming power to the lth GNSS signal power ratio (JSR) is written as follows
J S R = 10 log 10 { A J , m A l }
The jamming power to noise power ratio (JNR) for the mth jamming signal is written as follows
J N R = 10 log 10 { A J , m N 0 B I F }
where B I F is the RF front-end bandwidth. The analytical expression of the received signal in Equation (1) is defined as follows
X a ( t ) = X I F ( t ) + j X ^ I F ( t )
where the X ^ I F ( t ) is the Hilbert transform of X I F ( t ) ; the use of the analytic signal X a ( t ) can avoid the presence of cross-terms which could be generated by the interaction between positive and negative frequency components [19].

3. The Proposed Method

When it comes to the time-varying interference, the classical time-frequency analysis based on WVD and RSPWVD provides superior performance. However, WVD suffers from cross-term interference seriously when the analytic signals have two or more components [20]. The RSPWVD can reduce the cross-term interference, but it cannot deal with the signals with crossed or overlapped frequency by the peak detection method [18]. An example of RSPWVD with two chirps whose frequencies are crossed is shown as Figure 1. Figure 1a depicts that the peaks of two chirps are clear except the overlapped frequency part. In Figure 1b, it is obvious that the outline of two chirps in the overlapped frequency is blurred so that it is difficult to distinguish which signal the frequency of the overlapped part belongs to. Therefore, the proposed method based on Hough transform of RSPWVD is introduced [21,22,23].

3.1. Hough Transform

Hough transform has been widely used in line detection problems in images. In image I, as shown in Figure 2, the X-Y coordinate origin is at the center of the image with size M × L , then ( x , y ) is satisfied as follows
x = t M 2 y = f L 2
where f and t represent the respective frequency variable and time variable in the TF plane. Equation (9) in polar coordinates can be rewritten
x cos θ + y sin θ = ρ

3.2. The Combination of WVD and Hough Transform

WVD produces an energy distribution concentrated along a straight line. As a result, the problem of detecting the chirp interference can be turned into an issue of detecting the TF plane line and performed by Hough transform, which has been widely used for detecting chirp signal combined with the TF distribution [22,23]. The combination of WVD and Hough transform can be defined as follows
W H x ( f 0 , k ) = X a ( t + τ 2 ) X a ( t τ 2 ) e j 2 π ( f 0 + k t ) τ d τ d t
where X a ( t ) is an analytical signal and X a ( t ) is a complex conjugate of X a ( t ) , f 0 is the initial frequency, k is the chirp rate. Both k and f 0 can be deduced from Figure 2 by a geometric relationship, the results can be written
f 0 = [ L 2 M 2 tan θ + ρ cos θ ] Δ f
k = tan θ Δ f Δ t
where Δ f is the frequency resolution, which equals F s 2 M while Δ t is the time resolution and equal to 1 F s , F s is the sampling rate. The Hough transform of WVD maps the point ( t , f ) of the TF plane line onto the point ( ρ , θ ) in the parameter plane. As a consequence, a peak is formed at the point ( ρ , θ ) of the parameter plane. Once the peak value exceeds a predetermined threshold, it can be determined that there is a chirp interference.

3.3. Double Threshold Detection

This paper adopts the double threshold detection method named primary threshold and secondary threshold [24]. The GNSS signals are buried in thermal noise assumed to be zero mean, independent and identically distributed (IID). When interference is absent, the complex random variable X a ( t ) in Equation (8) is the zero mean and IID. Its magnitude spectrum ϕ ( k , l ) in the TF plane can be written [25]
ϕ ( k , l ) = Re 2 [ X a ( t k , f l ) ] + Im 2 [ X a ( t k , f l ) ]
where k = 1 , , N ;   l = 1 , , N ; N is the length of samples. X a ( t k , f l ) is the result of TF transform for X a ( t ) . The variable | ϕ ( k , l ) | follows a Rayleigh distribution. The primary threshold setting is written as [26]
η = ϕ ( k ) ¯ 4 π ln ( P f a l s e )
where
ϕ ( k ) ¯ = 1 N l = 1 N | ϕ ( k , l ) |
ϕ ( k ) ¯ is the average mean of ϕ ( k , l ) , P f a l s e is the false alarm rate. From Equation (15), the primary threshold can be determined by a predefined P f a l s e . Then, the data in X a ( t k , f l ) are separated to a new set which is higher than the primary threshold and the Hough transform is applied to these data. The secondary threshold in the Hough domain is used to detect the target peak and is based on the primary threshold. The P f a l s e is determined by both the primary threshold and secondary threshold, written as
P f a l s e = f ( η , ε )
where η is the primary threshold and ε is the secondary threshold. When the primary threshold is determined by the initial probability of false alarm, the secondary threshold is set by Monte Carlo simulation with fixed P f a l s e .

3.4. Joint Method Based on Hough Transform of RSPWVD

From the above analysis, not only WVD but also RSPWVD can be combined with Hough transform. However, WVD suffers from severe cross-term interference. As a result, it will affect the result of the Hough transform. Therefore, the Hough transform of RSPWVD is proposed [23]. The general expression of the Hough transform of the time-frequency transform can be written as
T F H X ( f 0 , k ) = T F X ( t , f 0 + k t ) d t
where T F H X ( f 0 , k ) is the TF distribution of the analytical signal X a ( t ) . From Equation (10), the chirp interference in the TF plane is written as polar coordinates as well as Cartesian coordinates. Therefore, Equation (11) can be written as
T F H X ( ρ , θ ) = T F X ( t , f ( ρ , θ ) ) d t
where f ( ρ , θ ) is the instantaneous frequency of the chirp interference in polar coordinates. Equation (18) as well as Equation (19) depicts that the line of the TF plane becomes a peak at point ( ρ , θ ) of the parameter plane by Hough transform. Conversely, a peak in the parameter plane represents a line in the TF plane and can be used to estimate the parameter of the line. As a result, this can be used to detect and identify chirp interferences. The joint method based on the Hough transform of RSPWVD is displayed as Figure 3 and is based on the following steps.
  • map X a ( t ) onto the TF plane by computing its RSPWVD
    X a ( t ) R S P W V D X ( t , f )
  • separate the data of R S P W V D X ( t , f ) into two sets and map the set above the primary threshold onto the Hough transform.
    R S P W V D X ( t , f ) T F H X ( ρ , θ )
  • search the peaks of T F H X ( ρ , θ ) ; if a peak exceeds the secondary threshold, record its point value ( ρ , θ )
  • estimate the parameter ( f 0 , k ) by the point value ( ρ , θ )
    T F H X ( ρ , θ ) T F H X ( f 0 , k )

3.5. Impact on the Acquisition Stage

In the case of a single satellite and a single chirp interference, the input signal X I F ( t ) in Equation (1) enters the ADC without considering the quantization effect and can be rewritten as
X I F ( n ) = 2 A C ( n τ ) D ( n τ ) cos [ 2 π ( f I F + f d ) n T s + φ I F ] + 2 A J cos [ 2 π ( f i n s ( t ) ) n T s + θ J ] + N ( n )
Then, the signal X I F ( n ) is multiplied by two orthogonal sinusoids as well as a local signal replica and is integrated; the results are as follows [27]
S I ( τ , f D ) = 1 N n = 0 N 1 r I ( n ) c ( n τ ) = r I ( τ ) h c ( τ ) S Q ( τ , f D ) = 1 N n = 0 N 1 r Q ( n ) c ( n τ ) = r Q ( τ ) h c ( τ )
where r I ( n ) = X I F ( n ) cos ( 2 π f D n ) , r Q ( n ) = X I F ( n ) sin ( 2 π f D n ) , f D = ( f I F + f d ) T s , ( ) denotes integral operation, N is the length of local code and h c ( τ ) is an equivalent filter. The Cross Ambiguity Function (CAF) is obtained [27]:
S ( τ , f D ) = S I ( τ , f D ) 2 + S Q ( τ , f D ) 2
When the Doppler shift f d and the code delay τ are correctly recovered, the detection probability that the statistical variable S ( τ , f D ) surpasses a fixed threshold β in the GNSS acquisition stage can be defined:
P det ( β ) = P ( S ( τ , f D ) > β )
In order to assess the impact of chirp interference on the detection probability, it is supposed that the GNSS signal and the chirp interference are known. The integration results of GNSS signals in Equation (24) can be written [27]:
S y C / 2 exp { j φ I F }
Similarly, the integration results of chirp interference in Equation (24) can be written [27]:
S J = k 1 A J 2 exp { j 2 π ( f in s + f D ) τ + j θ J + j θ 1 } + k 2 A J 2 exp { j 2 π ( f in s f D ) τ j θ J + j θ 2 }
where θ 1 = { H c ( f i n s + f D ) } and θ 2 = { H c ( f i n s + f D ) } . H c ( f ) is the Fourier transform of h c ( n ) . k 1 and k 2 can be written as
k i 2 = | H c ( ± f i n s + f D ) | 2 = | H c ( f ) | 2 δ ( f ( ± f i n s + f D ) ) d f = G s ( f ) G i ( f ) d f
where δ ( · ) is the delta Dirac, G s ( f ) = | H c ( f + f D ) | 2 and G i ( f ) = δ ( f ± f i n s ) . When the chirp interference appears, the CAF in Equation (25) follows a Rice distribution [27]
S ( τ , f D ) = x σ 2 exp { x 2 + α 2 2 σ 2 } I 0 ( x α σ 2 ) ,   x > 0
where α 2 = | S y + S J | 2 , σ 2 = σ o u t 2 2 , σ o u t 2 is the noise variance of ouput of equivalent filter h c ( n ) and σ o u t 2 = 1 N N 0 B I F , I 0 is the modified Bessel function with first kind and zero order. The detection probability can be defined:
P d ( β ) = β x σ 2 exp { x 2 + α 2 2 σ 2 } I 0 ( x α σ 2 ) d x = Q ( α σ ; β σ )
where Q ( ; ) denotes the Marcum Q function. From Equation (31), the chirp interference seriously affects the detection probability in the GNSS acquisition stage.

4. Results

To test the performance of the introduced method, an experiment under the condition where the BDS B1I signals are corrupted by two simulation chirps is carried out for several scenarios. The main hardware parameters of the down converter and ADC are recorded in Table 1.
First, the JNR is set to 0 dB, 6 dB and 12 dB. The instantaneous frequency of one chirp interference is from 50 MHz to 30 MHz, and the instantaneous frequency of the other chirp interference is from 30 MHz to 50 MHz. The sweep period is 2.56 us and the sample length is 512. The Hough transforms combined by WVD and RSPWVD are shown as Figure 4.
Figure 4a shows that the WVD of two chirps suffers from cross-term interferences seriously due to the interaction of different chirp signal components; the outline of two chirps are not clear. Figure 4b depicts the RSPWVD of two chirps, the outline of which are clear. The cross-terms are less than those in Figure 4a. Figure 4c presents the WVD of chirps with a higher JNR and its energy distribution is clear and suffers from less cross-term interferences compared with Figure 4a. Figure 4d represents a clear RSPWVD which suffers from fewer cross-term interferences compared with that in Figure 4c. Although the power of chirps is 12 dB above that of noise, the WVD of chirps in Figure 4e still suffers from cross-term interference, while in Figure 4f, the outline of RSPWVD is very clear and its energy distribution hardly suffers from cross-term interference.
However, Figure 4 shows that the common frequency (overlapped frequency) which belongs to chirp one as well as chirp two cannot be identified by the peak detection method. As a result, it is difficult to estimate instantaneous frequency in overlapped parts for each chirp. In order to deal with this problem, the experiment with the condition where the Hough transform is combined with WVD and RSPWVD is conducted; the experimental setup is as in Figure 4 and the results are shown in Figure 5.
Figure 5a shows the Hough transform of WVD and there are three peaks in the Hough domain. The peak on the right is the pseudo peak. It is a cross-term, which is integrated by Hough transform in the TF plane and formed into a peak. What is shown in Figure 5b is similar to that in Figure 5a. Although the RSPWVD can reduce cross-term interferences, the power of chirps is as strong as that of noise. As a result, the RSPWVD reallocates the noise energy which is integrated by the Hough transform in the TF plane and is formed into a strong peak. Figure 5c depicts two strong peaks and a weak pseudo peak, while in Figure 5d there are only two strong peaks without a pseudo peak. In Figure 5e,f, there are two peaks without a pseudo peak. Figure 5 depicts that although the WVD suffers from cross-term interferences, the Hough transform can help to reduce the effect of the cross-term interferences when the chirp signals are strong. In addition, the RSPWVD can eliminate the cross-term interference at the expense of its good localization and concentration properties; the combination of RSPWVD and Hough transform can effectively reduce the effect of cross-term interferences.
In order to assess the performance between the Hough transform of WVD and RSPWVD, the root mean square errors (RMSE) of the rate estimation and initial frequency estimation for chirp interference are used and displayed as a function of JNR in Figure 6. The setting parameters of the experiments are the same as in Figure 4. Figure 6a shows that the RMSE of the chirp one rate estimated by the Hough transform of WVD is the same as that estimated by the Hough transform of RSPWVD when JNR ranges from 0 dB to 12 dB. Similarly, the RMSE of the chirp two rate estimated by the Hough transform of WVD is the same as that estimated by the Hough transform of RSPWVD when JNR ranges from 0 dB to 12 dB. In Figure 6a, it is obvious that the Hough transform of RSPWVD as well as that of WVD provides the same RMSE of the rate estimation for each chirp when JNR ranges from 0 dB to 12 dB.
In Figure 6b, the RMSE of the initial frequency for chirp one estimated by the Hough transform of RSPWVD is close to 10 2 when JNR is equal to 0 dB and 2 dB, while the RMSE of the initial frequency for chirp one estimated by the Hough transform of WVD is about 0.5 × 10 2 when JNR is equal to 0 dB and 2 dB. However, when JNR equals or surpasses 4 dB, the Hough transform of RSPWVD and WVD provides the same RMSE of the initial frequency for chirp one. On the other hand, the RMSE of the initial frequency for chirp two estimated by the Hough transform of WVD is close to 0.4 × 10 2 when JNR ranges from 0 dB to 12 dB. The RMSE of the initial frequency for chirp two estimated by the Hough transform of RSPWVD is close to 0.1 × 10 2 when JNR ranges from 0 dB to 12 dB. Figure 6b depicts that the Hough transform of RSPWVD offers a poorer RMSE of the initial frequency estimation for chirp one compared with the Hough transform of WVD when JNR is below 4 dB and provides the same RMSE of the initial frequency estimation as that of WVD for chirp one when JNR equals or surpasses 4 dB. For chirp two, the proposed method provides a better RMSE of the initial frequency estimation compared with the Hough transform of WVD.
In addition, another experiment is performed by the proposed method under the condition where one chirp is adopted and its instantaneous frequency ranges from 30 MHz to 50 MHz. The sample length is 512 and the sweep period is 2.56 us, 5.12 us and 10.24 us. The results of RMSE are shown in Table 2. In Table 2, the estimation of initial frequency is close to 0.5 × 10 2 and the estimation of the chirp rate is about 0.6 × 10 2 when the sweep period is 2.56 us. However, when the sweep period increases to 5.12 us, the estimation of the initial frequency decreases by approximately 1.5 × 10 2 , and the estimation of chirp rate decreases by almost 0.7 × 10 2 . Finally, when the sweep period reaches 10.24 us, the estimation of initial frequency continues to decline by about 1.7 × 10 2 and the estimation of the chirp rate reduces to a level near 0.8 × 10 2 . From Table 2, it is obvious that the estimation accuracy of the initial frequency as well as the chirp rate declines as the sweep period increases.
The computational requirements of the Hough transform of WVD and RSPWVD are shown in Table 3 [18]. N is sample length, M is the length of the time window and L is the length of the frequency window. The computational requirements of the Hough transform of WVD consist of the computational requirements of the Hough transform as well as those of WVD. The output time-frequency matrix of WVD is N × N , the discrete points of θ from 0 to 2 π are K. From Equation (10), the Hough transform needs two complex multiplications and one complex addition. Therefore, the total requirements of the Hough transform are 2 N 2 K complex multiplications and N 2 K complex additions. Similarly, the total requirements of the Hough transform in case of RSPWVD are 2 N M K complex multiplications and N M K complex additions. From Table 3, the requirements of WVD are higher than those of RSPWVD because of N M and N L . As a result, the tatal requirements of the Hough transform of WVD are higher than those of RSPWVD.
Finally, the chirp signal is a continuous wave at any fixed time and can be mitigated by a notch filter based on the second-order direct form structure with two parameters named α and β which are determined by the power of chirp interference and the instantaneous frequency estimated by proposed method, respectively. The receiver operating characteristic (ROC) curves are shown in Figure 7 by Monte Carlo simulations to analyze the performance of the notch filter to mitigate the chirp interference under the condition where the BDS B1I signals are simulated with chirp interferences. The JNR is set to 12 dB and the C / N 0 is set to 43 dB-Hz. The integration time is set to 1 ms. The instantaneous frequency of one chirp ranges from 30 MHz to 50 MHz and the instantaneous frequency of the other ranges from 50 MHz to 30M Hz. The sweep period is 2.56 us.
Figure 7 depicts that the ROC curve in the two chirp scenario presents a worse acquisition performance than that in the one chirp scenario. However, when the notch filter is adopted, the ROC curves present an improved acquisition performance. The results depicts that the proposed method can effectively detect the one or two chirp interferences, especially for those with crossed frequency without a priori knowledge.

5. Conclusions

This paper presents a joint method based on the Hough transform of RSPWVD to detect time-varying interferences with crossed frequency for GNSS receivers with a single antenna. The analytical expression of initial frequency estimation and chirp rate estimation is presented and the double threshold detection is proposed as well.
The performance of the developed method has been assessed by experiments under conditions where the real BDS B1I signals corrupted by the simulated chirp interferences are collected by the GNSS software receiver. The actual performance of experiments has been shown by the quantitative metric RMSE of the parameter estimation. In addition, the effect of the sweep period on the estimation of the initial frequency and chirp rate has been analyzed.
The ROC curves have been used to assess the performance of the notch filter to mitigate chirp interference by Monte Carlo simulations. The results depict that the proposed method can effectively detect and identify the chirp interferences with crossed frequency and provide the same root mean square errors (RMSE) of the parameter estimation for chirp one and the improved initial frequency estimation for chirp two compared with the Hough transform of WVD when JNR equals or surpasses 4 dB. Furthermore, the RSPWVD method itself can provide better performance in reducing cross-term interference and needs less computational requirements compared with the WVD method.

Author Contributions

Q.L. contributed to design the algorithms, performed the experiment, analyzed the results of the experiments and wrote this paper. H.Q. provided valuable feedback and advice during the manuscript modification.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Grabowski, J. Personal privacy jammers: Locating jersey PPDs jamming GBAS safety-of-life signals. GPS World 2012, 23, 28–37. [Google Scholar]
  2. Marcos, E.P.; Konovaltsev, A.; Caizzone, S.; Cuntz, M.; Yinusa, K.; Elmarissi, W.; Meurer, M. Interference and Spoofing Detection for GNSS Maritime Applications using Direction of Arrival and Conformal Antenna Array. In Proceedings of the 31st International Technical Meeting of the Satellite Division of the Institute of Navigation, ION GNSS+ 2018 Conference, Miami, FL, USA, 24–28 September 2018. [Google Scholar]
  3. Son, P.W.; Rhee, J.H.; Seo, J. Novel Multichain-Based Loran Positioning Algorithm for Resilient Navigation. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 666–679. [Google Scholar] [CrossRef]
  4. Motella, B.; Pini, M.; Dovis, F. Investigation on the effect of strong out-of-band signals on global navigation satellite systems receivers. GPS Solut. 2008, 12, 77–86. [Google Scholar] [CrossRef]
  5. Fu, Z.; Hornbostel, A.; Hammesfahr, J.; Konovaltsev, A. Suppression of multipath and jamming signals by digital beamforming for GPS/Galileo applications. GPS Solut. 2003, 6, 257–264. [Google Scholar] [CrossRef]
  6. Palestini, C.; Pedone, R.; Villanti, M.; Corazza, G. Integrated NAV-COM Systems: Assisted Code Acquisition and Interference Mitigation. IEEE Syst. J. 2008, 2, 48–61. [Google Scholar] [CrossRef]
  7. Capozza, P.; Holland, B.; Hopkinson, T.; Landrau, R. A single-chip narrow-band frequency-domain excisor for a Global Positioning System (GPS) receiver. IEEE J. Solid State Circuits 2000, 35, 401–411. [Google Scholar] [CrossRef] [Green Version]
  8. 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] [CrossRef]
  9. Savasta, S.; Presti, L.; Rao, M. Interference Mitigation in GNSS Receivers by a Time-Frequency Approach. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 415–438. [Google Scholar] [CrossRef]
  10. Dovis, F. GNSS Interference Threats and Countermeasures; Artech House: Norwood, MA, USA, 2015. [Google Scholar]
  11. Betz, J.W. Effect of Partial-Band Interference on Receiver Estimation of C/N0: Theory. In Proceedings of the 2001 National Technical Meeting of the Institute of Navigation, Long Beach, CA, USA, 22–24 January 2001. [Google Scholar]
  12. Bastide, F.; Akos, D.; Macabiau, C.; Roturier, B. Automatic Gain Control (AGC) as an Interference Assessment Tool. In Proceedings of the 16th International Technical Meeting of the Satellite Division of The Institute of Navigation, ION GNSS2003, Portland, OR, USA, 9–12 September 2003. [Google Scholar]
  13. Lv, Q.S.; Qin, H.L. A novel algorithm for adaptive notch filter to detect and mitigate the CWI for GNSS receivers. In Proceedings of the ICSIP2018, Shenzhen, China, 13–15 July 2018. [Google Scholar]
  14. Jin, T.; Ren, J. Stability analysis of GPS carrier tracking loops by phase margin approach. GPS Solut. 2013, 17, 423–431. [Google Scholar] [CrossRef]
  15. Chen, Y.H.; Juang, J.C.; Seo, J.; Lo, S.; Akos, D.M.; Lorenzo, D.S.D.; Enge, P. Design and Implementation of Real-Time Software Radio for Anti-Interference GPS/WAAS Sensors. Sensors 2012, 12, 13417–13440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Park, K.; Lee, D.; Seo, J. Dual-polarized GPS antenna array algorithm to adaptively mitigate a large number of interference signals. Aerosp. Sci. Technol. 2018, 78, 387–396. [Google Scholar] [CrossRef]
  17. Zhang, L.; Li, B.; Huang, L.; Kirubarajan, T.; So, H.C. Robust minimum dispersion distortionless response beamforming against fast-moving interferences. Signal Process. 2017, 140, 190–197. [Google Scholar] [CrossRef]
  18. Lv, Q.S.; Qin, H.L. An improved method based on time-frequency distribution to detect time-varying interference for GNSS receivers with single antenna. IEEE Access 2019, 7, 38608–38617. [Google Scholar] [CrossRef]
  19. Hlawatsch, F.; Auger, F. Time-Frequency Analysis; Wiley Press: Hoboken, NJ, USA, 2005. [Google Scholar]
  20. Cohen, L. Time-Frequency Analysis: Theory and Applications; Prentice-Hall Press: Upper Saddle River, NJ, USA, 1995. [Google Scholar]
  21. Hough, P. Method and Means for Recognizing Complex Patterns. U.S. Patent 3069654, 1962. [Google Scholar]
  22. Barbarossa, S. Analysis of Multicomponent LFM Signals by a Combined Wigner-Hough Transform. IEEE Trans. Signal Process. 1995, 43, 1511–1515. [Google Scholar] [CrossRef]
  23. Barbarossa, S.; Scaglione, A. Adaptive time-varying cancellation of wideband interferences in spread-spectrum communications based on time-frequency distributions. IEEE Trans. Signal Process. 1999, 47, 957–965. [Google Scholar] [CrossRef] [Green Version]
  24. Moshe, E. Search Radar Track-Before-Detect Using the Hough Transform. Master’s Thesis, Naval Postgraduate School Monterey, Monterey, CA, USA, 1995. [Google Scholar]
  25. Henttu, P.; Aromaa, S. Consecutive mean excision algorithm. In Proceeding of the IEEE 7th International Symposium on Spread Spectrum Techniques and Application, Prague, Czech Republic, 2–5 September 2002. [Google Scholar]
  26. Vartiainen, J. Concentrated Signal Extraction Using Consecutive Mean Excision Algorithms. Ph.D. Thesis, University of Oulu, Oulu, Finland, 2010. [Google Scholar]
  27. Borio, D. GNSS Acquisition in the Presence of Continuous Wave Interference. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 47–60. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The RSPWVD of two chirps with crossed frequency. (a) Mesh of RSPWVD. (b) Contour of RSPWVD.
Figure 1. The RSPWVD of two chirps with crossed frequency. (a) Mesh of RSPWVD. (b) Contour of RSPWVD.
Sensors 19 01946 g001
Figure 2. Hough transform diagram.
Figure 2. Hough transform diagram.
Sensors 19 01946 g002
Figure 3. Block diagram of the joint method based on the Hough transform of RSPWVD to detect chirp interferences.
Figure 3. Block diagram of the joint method based on the Hough transform of RSPWVD to detect chirp interferences.
Sensors 19 01946 g003
Figure 4. The WVD and RSPWVD of BDS B1I signals with two chirps, C 0 / N is 43 dB-Hz and the JNR is 0 dB, 6 dB and 12 dB. (a) WVD, JNR is 0 dB. (b) RSPWVD, JNR is 0 dB. (c) WVD, JNR is 6 dB. (d) RSPWVD, JNR is 6 dB. (e)WVD, JNR is 12 dB. (f) RSPWVD, JNR is 12 dB.
Figure 4. The WVD and RSPWVD of BDS B1I signals with two chirps, C 0 / N is 43 dB-Hz and the JNR is 0 dB, 6 dB and 12 dB. (a) WVD, JNR is 0 dB. (b) RSPWVD, JNR is 0 dB. (c) WVD, JNR is 6 dB. (d) RSPWVD, JNR is 6 dB. (e)WVD, JNR is 12 dB. (f) RSPWVD, JNR is 12 dB.
Sensors 19 01946 g004
Figure 5. The Hough transform of WVD and RSPWVD with BDS B1I signals corrupted by two chirps, C 0 / N is 43 dB-Hz. (a) Hough transform of WVD, JNR is 0 dB. (b) Hough transform of RSPWVD, JNR is 0 dB. (c) Hough transform of WVD, JNR is 6 dB. (d) Hough transform of RSPWVD, JNR is 6 dB. (e) Hough transform of WVD, JNR is 12 dB. (f) Hough transform of RSPWVD, JNR is 12 dB.
Figure 5. The Hough transform of WVD and RSPWVD with BDS B1I signals corrupted by two chirps, C 0 / N is 43 dB-Hz. (a) Hough transform of WVD, JNR is 0 dB. (b) Hough transform of RSPWVD, JNR is 0 dB. (c) Hough transform of WVD, JNR is 6 dB. (d) Hough transform of RSPWVD, JNR is 6 dB. (e) Hough transform of WVD, JNR is 12 dB. (f) Hough transform of RSPWVD, JNR is 12 dB.
Sensors 19 01946 g005
Figure 6. RMSE of the chirp rate and initial frequency estimation for chirp interferences present in BDS B1I signals by Hough transform combined by WVD and RSPWVD versus JNR. The C 0 / N is 43 dB-Hz. (a) RMSE of chirp rate. (b) RMSE of initial frequency.
Figure 6. RMSE of the chirp rate and initial frequency estimation for chirp interferences present in BDS B1I signals by Hough transform combined by WVD and RSPWVD versus JNR. The C 0 / N is 43 dB-Hz. (a) RMSE of chirp rate. (b) RMSE of initial frequency.
Sensors 19 01946 g006
Figure 7. ROC curves for BDS B1I signals affected by chirp interference. The J N R is 12 dB, the C / N 0 is 43 dB-Hz. The coherent integration time is 1 ms.
Figure 7. ROC curves for BDS B1I signals affected by chirp interference. The J N R is 12 dB, the C / N 0 is 43 dB-Hz. The coherent integration time is 1 ms.
Sensors 19 01946 g007
Table 1. Down-converter digitizer parameters.
Table 1. Down-converter digitizer parameters.
ParameterValue
Bandwidth20 MHz
Intermediate frequency40 MHz
Down Converter Gain60 dB
Dynamic range in Down Converter70 dB
Sampling Rates200 MHz
Bits Per sample14
Table 2. Chirp detection for different sweep periods.
Table 2. Chirp detection for different sweep periods.
Sweep Period
(us)
Estimation of Initial Frequency
RMSE (Normalized)
Estimation of Chirp Rate
RMSE (Normalized)
2.560.0054380.005941
5.120.0150860.006581
10.240.0174630.007522
Table 3. Computational requirements of the Hough transform of WVD and RSPWVD.
Table 3. Computational requirements of the Hough transform of WVD and RSPWVD.
MethodComputational Requirements
WVD 2 N 2 + 2 N 2 log 2 N complex multiplications
4 N 2 log 2 N complex additions
RSPWVD N M ( 2 + L ) + N M / 2 log M complex multiplications
N M ( 4 + L ) + N M log M complex additions
WVD + Hough 2 N 2 + 2 N 2 log 2 N + 2 N 2 K complex multiplications
4 N 2 log 2 N + N 2 K complex additions
RSPWVD + Hough N M ( 2 + L ) + N M / 2 log M + 2 N M K complex multiplications
N M ( 4 + L ) + N M log M + N M K complex additions

Share and Cite

MDPI and ACS Style

Lv, Q.; Qin, H. A Joint Method Based on Time-Frequency Distribution to Detect Time-Varying Interferences for GNSS Receivers with a Single Antenna. Sensors 2019, 19, 1946. https://doi.org/10.3390/s19081946

AMA Style

Lv Q, Qin H. A Joint Method Based on Time-Frequency Distribution to Detect Time-Varying Interferences for GNSS Receivers with a Single Antenna. Sensors. 2019; 19(8):1946. https://doi.org/10.3390/s19081946

Chicago/Turabian Style

Lv, Qingshui, and Honglei Qin. 2019. "A Joint Method Based on Time-Frequency Distribution to Detect Time-Varying Interferences for GNSS Receivers with a Single Antenna" Sensors 19, no. 8: 1946. https://doi.org/10.3390/s19081946

APA Style

Lv, Q., & Qin, H. (2019). A Joint Method Based on Time-Frequency Distribution to Detect Time-Varying Interferences for GNSS Receivers with a Single Antenna. Sensors, 19(8), 1946. https://doi.org/10.3390/s19081946

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop