Cascaded Multitype Interferences Suppression Method Using Sparse Representation and Array Processing for GNSS Receiver

Interference suppression techniques have been intensively studied in nearly two decades due to their importance formaintaining the integrity and functionality of global navigation satellite system (GNSS). However, the interference suppression method applicable for the complex receiving environment in which there are multitype interfering signals has not been considered in most of the researches. To deal with this problembetter, a cascadedmultitype interferences suppressionmethod using sparse representation and array processing is proposed. In the first stage, according to the sparsity of the narrowband and modulated wideband interference signals, a novel parallel multichannel signal interference suppression method based on matching pursuit (MP) algorithm and a design strategy for the overcomplete dictionary are proposed tomitigate the interferences with sparse features.Then, the minimum power distortionless response (MPDR) beamformer is employed in the second stage to suppress the residuary interferences (such as Gaussian noise interferences). Compared with existing algorithms, the proposed method can not only effectively suppress the interference arriving from the samedirectionwith the desired signal and increase theDegree of Freedom (DoF) of the array antenna, but also introduce no distortion into the navigation signal. The effectiveness of the proposed method is illustrated by theoretical analysis and several simulation results.


Introduction
The global navigation satellite system (GNSS) has been widely used in military and civilian applications; however, the problems caused by the vulnerability of satellite navigation signals are increasingly serious.Therefore, interference suppression techniques are necessary to ensure the reliability, accuracy, and continuity of GNSS services.According to the characteristics of the GNSS, there are two ways to improve the capability of radio frequency (RF) interference suppression [1].On the one hand, we can improve the design of navigation satellites to enhance the interference suppression performance of GNSS, such as increasing the signal transmission power, optimizing the structure of the signal; but these methods are too complex, and the cost of time and materials is too huge.On the other hand, we can improve antijamming performance of GNSS receivers.These methods are gaining significant attention due to its effectiveness for interference from multiple, strong interferers plus multipath suppression.
The interference suppression techniques for GNSS receivers include time-domain processing, transform-domain processing, spatial processing, and spatial-time processing [1,2].The methods based on time-domain and transform-domain processing have been widely studied (e.g., [3,4]).Although they are effective for suppressing narrowband interference signals, their performance degrades when they deal with wideband interferers such as Gaussian interferers or when the interfering signal changes rapidly in time-domain and frequency domain (e.g., swept continuous wave interference).
On the contrary, the interference suppression techniques using an antenna array can effectively suppress both narrowband and wideband interferences regardless of their time and frequency characteristics.One of the most effective spacebased processing methods has been referred to as the minimum variance distortionless response (MVDR) beamformer [5].The MVDR beamformer has a distortionless response for the desired signal while rejecting all interfering signals arriving from other directions.Antenna array processing in 2 International Journal of Antennas and Propagation GNSS applications has been mostly centered on interference suppression [6].Reference [7] drew the attention on utilizing minimum power distortionless response (MPDR) beamformer to suppress interfering signals whose power is significantly higher than that of the GPS signals.However, spatial processing methods suffer from hardware complexity, and cannot effectively suppress interference whose direction of arrival (DOA) is close to the GNSS signal.To deal with these problems, techniques employing both time (or frequency) and spatial domain processing have been of great interest since [8] was published.These methods combine spatial and temporal filters to suppress more narrowband interference by increasing the Degree of Freedom (DoF) of the array without physically increasing the number of antenna array elements.Besides the superior advantages of space-time processing, these methods can introduce serious distortion and biases into the GNSS signal and pseudorange measurements [9].To reduce this distortion, many methods have been proposed, but these methods increase the complexity of the algorithm and reduce the DoF of interference suppression [10,11].
In addition, with the rapid development of jamming technology and the increasingly complex electromagnetic environment, there are multiple types of interferences existing simultaneously rather than a single type of interference in the environment.Although the interference suppression methods mentioned in the prior context improve the performance of GNSS receivers, they still face the following problems when dealing with multiple types of interferences: (1) the cost of hardware or space is huge; (2) the methods employing both time/frequency and spatial domain processing may introduce biases and distortion into the GNSS signals.
To solve the above problems, [12,13] drew the attention on cascaded interference suppression methods based on frequency domain and spatial domain.These methods are "simple cascade structure" of the conventional algorithms, and one of the obvious shortcomings is that when the power of wideband interference is higher than that of the narrowband interference, the interference mitigation methods in the prestage are not able to detect and cancel the narrowband interference.
In this paper, a cascaded multitype interferences suppression method using sparse representation and array processing for GNSS receivers is proposed.Firstly, the signal sparse decomposition theory [14] is introduced into array signal processing, which is gaining significant attention and has been successfully applied to many fields, such as clutter and jamming suppression for airborne radar [15] and signal detection [16].A novel overcomplete dictionary composed of linear frequency modulation atoms is designed.And a parallel multichannel signal interference suppression method based on matching pursuit (MP) algorithm is proposed, which can effectively detect and suppress the narrowband and modulated wideband interferences even when they fall in the other wideband interferences (e.g., Gaussian jammers).Secondly, the MPDR beamformer is employed in the second stage to mitigate the residual interference, which enhances the steering gain of the DOA of interest.Meanwhile, the null gain is set at DOA of interference.Finally, in order to verify the effectiveness of the proposed method, it is compared to the well-known MPDR beamformer [7] and the distortionless space-time processor proposed in [10].

Signal Model
The signal received by the GNSS receiver is the aggregate of the satellite navigation signal, interfering signals, and the thermal noise.Without loss of generality and for the sake of simplicity, assume that they are independent of each other and there is only one GNSS signal.Considering an arbitrary antenna array with  elements, the complex baseband representation of the received signal can be presented as where b  ,  = 1, 2, . . ., , indicates the steering vector of th interfering signal and   ,  = 1, 2, . . ., , is the th interfering signal. represents the GNSS signal and n is the receiver white Gaussian noise vector.And a is the steering of the GNSS signal where  is the wavelength of the signal; z  ,  = 1, 2, . . ., , is a 3 × 1 unit vector pointing to the th antenna element; d is a 3 × 1 unit vector representing the DOA of the GNSS signal; "(⋅)  " denotes transpose.

The Proposed Method
In the coexistence of multiple types of interfering signals, we can suppress them according to their different characteristics in different dimensions.In view of that, a cascaded multitype interferences suppression method using sparse representation and array processing for GNSS is proposed, and the structure of the proposed method is as shown in Figure 1.
In the first stage, in order to save the spatial DoF of the antenna array, the received signals are sparsely decomposed to suppress interference signals with sparse features according to the sparse characteristic of narrowband and wideband modulation interfering signals.To make the MP-based sparse decomposition be suitable for multitype interferences suppression scenario, we redesign the overcomplete dictionary in accordance with the characteristics of the interference signal, and the necessary conditions for interfering signals detection are analyzed, especially for the scenario that the energy of the interference signals with sparse feature is smaller than that of the wideband Gauss interferences.Then on the base of these, a parallel multichannel sparse decomposition based on MP and the termination condition is proposed.
In the second stage, the MPDR beamformer is employed to suppress Gaussian noise interferences and other residuary interferences by utilizing the spatial DoF of the antenna array.theory was proposed.And with the flexibility of the redundant dictionary and the adaptability to signals, it has gained more and more attentions.MP is a kind of iterative "greedy" algorithm.At each iteration, the best matching atom that is the most similar to the residual signal is selected and regarded as one of the components of the sparse representation.

Signal Sparse Decomposition
Let  ∈  be the signal to be decomposed, where  is a Hilbert space, and let  be the overcomplete dictionary defined by  = (g   )   ∈ Γ, where Γ is the parameter set of the atoms and ‖g   ‖ = 1 is named an atom.And the basic steps of MP are as follows.
(1) Choose an atom g  0 from  which best matches the ; that is, Then the signal can be decomposed into covariates on the best atom g  0 and the residual; that is, (2) Decompose the residual iteratively, for the ( + 1)th iteration (3) If termination condition is reached after the  steps decomposition, the signal can be expressed as Since ∑ −1 =0 ⟨  , g   ⟩g   ≫   , the signal can be decomposed into a sum of dictionary elements: Formula ( 6) and  ≪   , where   is the size of , represent the key idea of sparse representation.It should be noted that, with the decomposition, namely, as  is increasing, the residual    is gradually decreasing, until it disappears.

The Proposed Atoms of the Overcomplete Dictionary.
In order to facilitate the analysis, only one single-tone or LFM interfering signal is considered in the formulations below.And assume that there are  wide band Gaussian noise interfering signals.Because the strength of navigation satellite signals is usually very weak, there is no need to be concerned about obtaining signal-free samples.Then, for  snapshots, at the th antenna, the signal in quadrature  path can be written as where  is the thermal noise with zero mean and   variance; and   is the th snapshot data, expressed as International Journal of Antennas and Propagation in which   is the linear modulation frequency ratio;   represents the fixed frequency;   is the phase; ,  = 1, 2, . . .,  represents the number of samples;  is sampling period; p is the energy of the interfering signal; and   (⋅) is Gaussian interfering signals with zero mean and   variance.
According to the characters of interference signals, select a series of linear frequency modulated signals as atoms in the overcomplete dictionary, which could be expressed as where   is the normalized coefficient; represents the search range of the linear modulation frequency ratio;   ∈ [ 0 − /2,  0 + /2] is fixed frequency parameters and  and  0 represent the bandwidth and the center frequency of the received signal, respectively;   ∈ (−,] is phase parameters.The values of   ,   , and   form the parameter set Γ and are selected at regular intervals according to their respective reasonable ranges and the number of atoms.The more the number of atoms is, the higher the decomposition accuracy is.

Analysis of Interference Detection
Performance.The signal to be decomposed contains only one signal that can be sparse in the overcomplete dictionary.Then ( 8) can be rewritten as where Then, ℵ can be treated as a Gaussian noise with zero mean and  variance The best matching atom is selected according to formula (3) Assume that the closest matching atom, cos(2(   +   ) +   ), has been known.Then where Δ, Δ, and Δ represent the linear frequency modulation ration spacing, frequency spacing, and phase spacing, respectively.And since  is large enough, we can suppose that the normalized coefficient of each atom,   , is equivalent to each other in the following analysis.In general signal processing,  satisfies the Nyquist sampling theorem or the bandpass sampling theorem and  is large enough, so the latter is much less than the former in the polynomial; then In order to further understand the nature of R  , it is analyzed in the following three cases.
(1) Let Δ = Δ = 0, and   = 0 and   = 0; then Obviously, R  is a periodic function over Δ, and the period is 1/.When Δ = / or Δ = −/ − 2  , the maximum value of R  can be obtained When the sampling rate is 8 times as much as   , the relation between R  and Δ is as shown in Figure 2.
where R  is a periodic function over Δ and the period is 2.The relation between R  and Δ is as shown in Figure 3.The maximum value of R  is (3) Let Δ = Δ = 0, and   = 0 and   = 0; then where R  is a periodic function over Δ and the period is 1/ Accordingly, R  , R  , and R  do not have periodicity within the range defined by overcomplete dictionary.In other words, R  does have only maximum value in the search range of the overcomplete dictionary.From the above analysis, it is shown that there is only one best matching atom in the overcomplete dictionary for each sparse representation of the interference signal.And the maximum value of R  is the inner product of the signal to be decomposed and the best matching atom In addition, since G is the Gauss white noise with zero mean and  variance, then Therefore, in order to effectively detect the interfering signal, the number of samples  should satisfy that is, 3.1.4.Design Strategy of the Hierarchical Adaptive Overcomplete Dictionary.Although signal sparse decomposition based on MP has better performance with the increase of the number of atoms in the overcomplete dictionary, the improvement of the performance for such ways is at the cost of the increase of complexity.In order to reduce the amount of computation while maintaining the accuracy, we propose a design strategy of the hierarchical adaptive overcomplete dictionary.In step 1, an overcomplete dictionary is generated according to the prior information and the signal is decomposed using MP algorithm.In step 2, the parameters of the best atom obtained in step 1 serve as the new prior information to generate the next level overcomplete dictionary; repeat steps 1 and 2 until the overcomplete dictionary which satisfies the accuracy requirement is obtained.
Taking the fixed frequency parameter as an example, assume that the bandwidth of the desired signal is 2 MHz, while the decomposition accuracy is 0.0001 MHz.The number of fixed frequency parameters in the overcomplete dictionary is 2 × 10 4 .However, conduct the three-layer adaptive overcomplete dictionary in the following steps.
Step 1. Decompose the received signal with 0.02 MHZ precision in 2 MHz bandwidth, and get the primary best frequency parameter  1 .
Step 2. Generate new frequency parameters dictionary with 0.001 MHz precision at 0.02 MHz bandwidth centered in  1 , then decompose the received signal with the new dictionary, and get the best frequency parameter  2 .
Step 3. Generate new frequency parameters dictionary with 0.0001 MHz precision at 0.001 MHz bandwidth centered in  2 , then decompose the received signal, and get the frequency parameter meeting the accuracy requirements.
The total number of fixed frequency parameters in the 3stage overcomplete dictionary is 130.The computation of the latter is 0.65 percent of the former.

The Best Atom Search Strategy for Multichannel Signal.
Because the GNSS receiver used in the proposed method is equipped with the array antenna, there are multichannel signals to deal with.If these signals are decomposed independently, the calculation is very great.Fortunately, array signal processing theory demonstrates that there is a potential correlation among the signals of each channel.In particular, they have the same frequency, similar amplitude, and distinct phases.Therefore, we can extend the concept of single channel signal sparse decomposition to the multichannel signal by using the correlation among the signals of each channel.The detailed implementation steps are as follows.
Step 1. Divide the overcomplete dictionary into  subovercomplete dictionary, and the division principle is expressed as where  = 1, 2, . . ., ;  is the number of antenna; Step 2. Each  signal at th antenna corresponds to a subdictionary; and use the conventional MP algorithm to select the subbest matching atom, respectively, where x   is the  signal vector at th antenna. Step Then, the retrieval interference from the best atom in each channel can be expressed as And the  signal in each channel can be written as where   is the number of iterations.
International Journal of Antennas and Propagation

The Terminate Condition of MP.
The conventional terminative condition of the MP algorithm is to determine whether the number of iterations or the energy of the residual signal meets the requirements.However, in complicated electromagnetic environment, since the number of signals that can be sparse representation is unknown, the number of iterations cannot be preset; and when the energy of the signals which could be sparse representation is smaller than that of the other interferences, the residual signal energy may have little difference with the original signal energy, so the existing terminative condition is not able to guarantee the effectiveness of the algorithm.Hence, a ratio principle is introduced to judge whether the ratio of inner product of the best matching atom obtained in the  step and the signal to be decomposed to the sum of residual signal meets the threshold.It can be expressed mathematically as where sum(⋅) represents the sum function.Formulas (26) and ( 27) show that when the signal to be decomposed contains the interference signal that can be sparse representation in the overcomplete dictionary, inevitably there is a best matched atom making  > 1.

Parallel Multichannel Signal Interference Suppression
Method Based on MP.Based on the conclusions above, a multichannel signal interference suppression method based on sparse decomposition is proposed.The basic steps are as follows.
Step 1. Algorithm initialization: set up the overcomplete dictionary classification layer, search accuracy, and the threshold value of the termination conditions.
Step 2. Generate the global overcomplete dictionary .Step 3. Use parallel multichannel signal sparse decomposition method to decompose array signals.

International Journal of Antennas and Propagation
Step 4. Judge the accuracy of decomposition; if it reaches the preset accuracy, then take the next step; otherwise, regenerate global overcomplete dictionary according to the information acquired, and iterate to step 3.
Step 5. Shift the phase of the best atoms 90 ∘ to get the best atoms of the corresponding  signal; then the  signal can be written as Step 6. Determine whether the termination condition is met; if not, iterate to step 2; otherwise, take the next step.
Step 7. Output the residual signal of each channel, which can be written as where  = 1, 2, . . ., .

MPDR Beamformer.
In the second stage, by applying the spatial filtering, the beam of receiver antenna arrays is pointed towards the GNSS satellite and away from interferers to protect the GNSS signal and reject interferences.In GNSS applications, MPDR beamformer is one of the powerful approaches available to suppress interfering signals while maintaining desired signals due to its effectiveness for interference suppression without considering the structure and direction of the interfering signals.And the optimization problem for the MPDR beamformer can be expressed as where w represents the array weight vector; a is defined by formula (2); R is the spatial covariance matrix of the residual signals obtained by stage 1, which can be expressed by where X = [(x  1 )  (x  2 )  ⋅ ⋅ ⋅ (x   )  ]; "" denotes conjugate transpose; M is the size of snapshot data.Then the optimal weight vector is

Simulation Results and Analysis
In order to verify the effectiveness of the proposed method, three simulations have been conducted.In all simulations, a linear half-wavelength space antenna array with 5 elements has been considered.A navigation signal operates on 1.023 MHz with bandwidth of 2 MHz and C/A code rate of 1.023 MHz, whose incident angle is 80 ∘ and Doppler frequency is 2 KHz.The analog baseband signal with SNR = −20 dB is sampled at 16.328 MHz, and the number of snapshots is cut to 16328 points.The parameters of interference signals are shown in Table 1, and the parameters involved in the proposed algorithm are shown in Table 2.According to the parameters in Table 2, we can obtain the computation complexity of the MP stage for each iteration.And the computational complexities of the proposed method and the conventional algorithm are given in Table 3.    6 shows the relationship between terminate condition and the iteration number of the MP stage.It can be found that the values of  in the first three iterations are greater than the termination threshold, and it consists with the simulation conditions.Figure 7 shows the time-domain waveforms of the original interfering signals (Interferences 1, 2, and 4) and the residual signals in theory; Figure 8 shows the time-domain waveforms of the estimated interfering signals and the residual signal by the proposed method.Comparing Figures 7 and 8, it can be shown that the proposed method in this paper can exactly and effectively estimate the interfering signal.And the errors between these two signals can be measured with the normalized mean square error (NMSE), which are shown in Table 4. Accordingly, we can assume that there are not multitone and LFM interfering signals in the residual signal.

Simulation 2.
In this simulation, the proposed method is compared to the well-known space-only MPDR (S-MPDR) beamformer and three simulation scenarios are considered.In scenario 1, Interferences 3, 4, 5, and 6 are used.In other words, there is not any interfering signals with the same direction as the GNSS signal and the number of interferences is less than that of antenna elements.In scenario 2, Interferences 1, 4, 5, and 6 are adopted.It means that there is one interfering signal with the same direction as the GNSS signal.In scenario 3, Interferences 2, 3, 4, 5, and 6 are used.It implies that the number of interferences is equal to that of antenna elements.
The normalized correlation peaks after interference suppression are shown in Figures 9-14 for each scenario, which are used to measure the performance of signal acquisition after interference suppression by these two methods.For scenario 1, Figures 9 and 10 show that the two methods can effectively suppress the interferences.However, Figures 11 and  13 show that the MPDR beamformer has failed to capture the GNSS signal for both scenario 2 and scenario 3. Figures 12 and  14 show that the proposed method can effectively suppress the interferences and does not distort the normalized correlation peak of the GNSS signal even there is one interfering signal with the same direction as the GNSS signal and when the number of interference is more than the DoF of the antenna array.

Simulation 3.
In this simulation, the proposed method is compared to the well-known space-time MPDR (ST-MPDR) beamformer [7] and the distortionless space-time adaptive processor (DST-MPDR) [10].The number of time delay taps is 7.And three simulation scenarios are conducted.In scenario 1, Interferences 1, 2, 4, 5, 6, and 7 are considered, which means that there is one single-tone interfering signal with the same direction as the GNSS signal and the number of wideband interferences is less than that of antenna elements.In scenario 2, Interferences 1, 2, 4, 5, 6, 7, and 8 are used.In other words, the number of wideband interferences is more than that of antenna elements.In scenario 3, Interferences 1, 2, 5, 6, 7, and 9 are used.It means that there is one wideband interfering signal with the same direction as the GNSS signal.
The normalized correlation peaks after interference suppression by the three methods for each scenario are shown in Figures 15-23, which are used to measure the performance of interference suppression.For scenario 1: Figure 15 shows that the space-time MPDR processing method causes the distortion and shift of the correlation peak; Figure 16 shows that although the distortionless space-time adaptive processor introduces code biases into the GNSS signal, it does not distort the correlation peak of the GNSS signal; Figure 17 shows the proposed algorithm does not affect the shape and the position of the correlation peaks.For both scenario 2 and scenario 3: Figures 18 and 19 and Figures 21 and 22 show that space-time MPDR processing and the distortionless spacetime adaptive processor fail to suppress the interference and it is failed to capture the GNSS signal; Figures 20 and 23 show that the proposed method is still effective since the modulated wideband can be canceled by sparse decomposition in the first stage.

Conclusion
Based on the sparsity of the interfering signals and the advantage of the spatial processing, a novel cascaded multitype interferences suppression method using sparse representation and array processing is proposed and examined in this paper.Firstly, the parallel multichannel signal interference suppression method based on MP is proposed, which can save the spatial DoF of the antenna array by effectively detecting and canceling the narrowband and modulated wideband interference even when they vanish into the Gaussian interferences.Then, the MPDR beamformer is employed to suppress the residuary interferences (such as Gaussian noise interferences) by utilizing the spatial DoF of the antenna array.Numerical simulations show that the proposed method not only can suppress more interferences, but also does not affect the shape and position of the correlation peaks.Compared with the space-only MPDR International Journal of Antennas and Propagation  beamformer, space-time MPDR beamformer, and the distortionless space-time adaptive processor, the proposed method is able to deal with multitypes interferences more effectively and can suppress the interference with the same direction as the desired signal.Therefore, the proposed method is able to effectively improve the interference suppression ability of GNSS receiver while reducing the cost of space and hardware.

Figure 1 :
Figure 1: Block diagram of the proposed method.

Figure 7 :
Figure 7: Time-domain waveform of original interfering signals.

Figure 10 :Figure 11 :
Figure 10: Correlation peaks after interference suppression by the proposed method for scenario 1.

Figure 12 :Figure 13 :
Figure 12: Correlation peaks after interference suppression by the proposed method for scenario 2.

Figure 14 :Figure 15 :
Figure 14: Correlation peaks after interference suppression by the proposed method for scenario 3.
of fixed frequency in the first layer (corresponding range) 0.01 MHz (1.023 ± 1 MHz) The decomposition accuracy of fixed frequency in the second layer (corresponding range) 0.001 MHz ( 1 ± 0.01 MHz) The decomposition accuracy of fixed frequency in the third layer (corresponding range) 0.0001 MHz ( 2 ± 0.0005 MHz) The decomposition accuracy of phase in the first layer (corresponding range) 0.1 ((−, ]) The decomposition accuracy of phase in the second layer (corresponding range) 0.01 ( 1 ± 0.05) The decomposition accuracy of phase in the third layer (corresponding range) 0.001 ( 1 ± 0.005) The decomposition accuracy of linear modulation frequency rate in the first layer (corresponding range) 10 7 (0 ± 10 9 ) The decomposition accuracy of linear modulation frequency rate in the second layer (corresponding range) 10 6 ( 1 ± 5 × 10 6 ) The decomposition accuracy of linear modulation frequency rate in the third layer (corresponding range) 10 5 ( 2 ± 5 × 10 5 )

Figure 20 :
Figure 20: Correlation peaks after interference suppression by the proposed method for scenario 2.
According to the frequency parameters of the global best matching atom, calculate the unique phase and amplitude parameters of the signals in each channel     ⟨x   , g  (  ,   , )⟩     ⟨x   , g  (  ,   , )⟩     .
3. Compare the inner products of subbest matching atoms and corresponding signals, and select the global closest matching atom g   (  ,   , ), which satisfies the following condition: max     ⟨x   , g  sub ⟩     .(31)Step 4.

Table 2 :
The parameters involved in the proposed algorithm.

Table 3 :
Computational complexities of the proposed method and the conventional MP.