Next Article in Journal
A Cartographic Perspective on the Planetary Geologic Mapping Investigation of Ceres
Next Article in Special Issue
Improved Medium Baseline RTK Positioning Performance Based on BDS/Galileo/GPS Triple-Frequency-Only Observations
Previous Article in Journal
Assessment and Data Fusion of Satellite-Based Precipitation Estimation Products over Ungauged Areas Based on Triple Collocation without In Situ Observations
Previous Article in Special Issue
Assessing the Nonlinear Changes in Global Navigation Satellite System Vertical Time Series with Environmental Loading in Mainland China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Robust Space-Time Adaptive Processing Method for GNSS Receivers in Coherent Signal Environments

1
School of Information and Control Engineering, China University of Mining and Technology, Xuzhou 221116, China
2
School of Instrumentation Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(17), 4212; https://doi.org/10.3390/rs15174212
Submission received: 12 July 2023 / Revised: 24 August 2023 / Accepted: 25 August 2023 / Published: 27 August 2023
(This article belongs to the Special Issue New Progress in GNSS Data Processing Technology and Modeling)

Abstract

:
In the coherent signal environments caused by multipath propagation, the interference suppression performance of the global navigation satellite systems (GNSS) receivers decreases sharply. In this paper, a robust space-time adaptive processing (STAP) method for GNSS receivers is proposed to suppress interferences in coherent signal environments, by using the modified space-time two-dimensional iterative adaptive approach (ST2D-IAA) spectrum estimation. This method applies the IAA algorithm to the ST2D signal model of GNSS receivers, and further modifies the ST2D-IAA algorithm to accurately estimate the power spectrum and noise power simultaneously. The space-time interference-plus-noise covariance matrix (STINCM) is reconstructed by using the estimated power spectrum and noise power in the interference angle region. Based on the reconstructed STINCM, we construct the STAP beamforming optimization problem for the space-time steering vector (STSV) error vector, and further correct the STSV of GNSS signal. Finally, the weight vector of STAP beamforming is calculated by using the reconstructed STINCM and the corrected STSV of GNSS signal. Simulation results show that the proposed method can suppress interferences in coherent signal environments and outperforms the current methods.

1. Introduction

The global navigation satellite systems (GNSS) use the advantages of satellites to provide real-time, high-precision navigation, positioning and timing services, which play a vital role in the economic and military fields [1,2,3]. The current GNSS mainly includes the United States global positioning systems (GPS), China’s Beidou Systems, Europe’s GALILEO systems and Russia’s GLONASS systems. However, the accuracy and reliability of GNSS receivers can be easily affected by all kinds of intentional and unintentional interferences [4]. With the significant increase in the types and quantities of interference, GNSS applications have posed a great threat, especially the rapid increase in cheap jamming devices [5]. Therefore, it is urgent to improve the interference suppression ability of GNSS receivers. An effective method to suppress interferences is to use antenna array processing technology, which can form a nulling in the direction of interference [6].
A large number of interference suppression algorithms for GNSS receivers based on antenna arrays have been proposed. Among these algorithms, space only processing (SOP) and space-time adaptive processing (STAP) are two typical representatives [7]. The SOP algorithms [8] use only one spatial filter to process the observed signals, mainly divided into power inversion algorithms and adaptive beamforming algorithms. The power inversion algorithms [9] do not need to know the array structure and satellite signal information in advance, and utilizes the characteristics that the intensity of the satellite signal is lower than the noise and interferences when it reaches the receivers. However, this kind of algorithms cannot guarantee that the mainlobe of the beampattern always points to the satellite signal. Adaptive beamforming algorithms are proposed to enhance desired satellite signals while suppressing interference signals which are based on the criteria of minimum variance distortionless response (MVDR), minimum power distortionless response (MPDR) and maximum signal-to-interference-plus-noise ratio (SINR). In practice, there may be direction errors, local scattering, wavefront distortion and limited number of sampling snapshots, which lead to the covariance matrix mismatches and steering vector errors. In order to improve the beamforming performance, a large number of robust adaptive beamforming methods have been proposed to eliminate these errors, such as diagonal loading [10], eigenspace-based [11], uncertainty set-based [12], interference-plus-noise covariance matrix (INCM) reconstruction [13,14,15,16] methods. Although these SOP algorithms can effectively suppress multiple interferences, the degrees-of-freedom of interference suppression is limited by the number of antenna array elements.
In order to overcome the shortcomings of the SOP algorithms, the STAP algorithms are proposed to increase the number of interference suppression degrees-of-freedom in GNSS receivers. The STAP algorithms work by placing a finite impulse response (FIR) filter behind each sensor without the need to increase the number of physical array antennas [17,18,19]. The distortionless STAP method can guarantee the distortionless response of GNSS signal and suppress the interference signals [20]. The essence of this method is to ensure the linearity of the phase by constraining the coefficients of the FIR filter. However, the imprecise prior information caused by the direction errors will result in the mismatch of the steering vector, which further deteriorates the performance of STAP method [21]. Therefore, it is of great significance to study a robust STAP method that can deal with steering vector errors. In the vicinity of GNSS receivers, due to the obstruction of buildings or obstacles, GNSS signal will be refracted and reflected, resulting in the coherent interference, which seriously reduces the positioning performance of GNSS receivers [22].
The multipath signal can be used for remote sensing purposes in GNSS interferometric reflectometry, such as the estimation of lake ice thickness in [23]. In general, the GNSS multipath signal is coherent with the direct line-of-sight signal. As a result, the coherent interference causes the spatial spectrum diffusion in the angle region, which reduces the accuracy of the angle estimation algorithms. In order to deal with the coherent interference problems, the spatial smoothing technique [24] resorves the correlation by restoring the rank of the signal covariance matrix. Forward/backward spatial smoothing techniques [25] construct an average forward/backward covariance matrix by dividing the uniform linear array into multiple coincident subarrays, thereby overcoming the defect of rank deficiency. However, this kind of techniques deal with the coherence problems at the cost of the number of effective degrees-of-freedom, and still cannot deal with the coherent interference completely in the low signal-to-noise ratios (SNR) environments. Since the parametric maximum likelihood beamforming technique [26] can handle the coherent interference, it can be considered as an alternative to the methods described above. This method parameterizes all the incident signals, but the derived beamforming optimization problem is a complex high-dimensional expression, which dramatically increases the computational complexity of the system. In [27], a polarization smoothing method is proposed to restore the rank of the matrix by averaging the array data across the polarization domain. Another polarization smoothing method is developed in [28] by performing polarization differential smoothing on the array data, which can work with colored noise. However, both of these methods cannot estimate the polarization parameter. In [29], the 2D direction-of-arrival (DOA) estimation problem of coherent sources is revisited via a polarized uniform rectangular array, where the rank-deficiency of the source matrix is solved by the rearrangement of the array data.
In recent years, various deformation forms of beamformers based on sparse methods [30] have been extensively studied, also known as iterative adaptive approaches (IAA) algorithms, showing good robustness to coherent interference and array perturbations. The IAA algorithms use the weighted minimum variance algorithm to update variables iteratively, and has been widely used in many fields such as multiple-input multiple-output radar imaging, source localization, blood flow velocity estimation, synthetic aperture radar imaging and channel prediction. A generalized IAA [31] is proposed for robust spectral analysis, which utilizes a weighted l p norm to achieve high accuracy and outlier suppression. The IAA-based adaptive beamforming algorithms [32] have robust performance for model errors. The idea of these algorithms is to iteratively estimate the spatial spectrum, which is then used to reconstruct INCM and correct the direction of the desired signal. A virtual uniform linear IAA [33] is proposed for the design of a coprime array robust adaptive beamformer that can directly process snapshots data in the longest virtual uniform linear array. A covariance matrix reconstruction method based on IAA is proposed in [34], which can overcome the direction errors and deal with coherent interference when linear arrays are mutually coupled. So far, there is no literature on the application of IAA algorithm in space-time two-dimensional (ST2D) signal model of GNSS receivers.
In order to deal with the interference suppression of GNSS receivers in coherent signal environments, a robust STAP method based on modified ST2D-IAA spectrum estimation is proposed. First, the IAA algorithm is applied to the ST2D signal model of GNSS receivers, and the ST2D-IAA algorithm is further modified to accurately estimate the power spectrum and the noise power simultaneously. Secondly, the estimated power spectrum and the noise power of the interference angle region are used to reconstruct the space-time interference-plus-noise covariance matrix (STINCM) of GNSS receivers. Then, the STAP beamforming optimization problem is constructed, and the space-time steering vector (STSV) error vector is obtained after solving this optimization problem, which is used to correct the STSV of GNSS signal. Finally, the weight vector of STAP beamforming is calculated. Simulation results show that the proposed method can suppress interferences in coherent signal environments and is robust to model errors.
Notations: The lower-case (upper-case) bold characters are used to represent the vectors (matrices). C denotes the set of the complex numbers. The superscripts ( · ) T , ( · ) H , ( · ) * and ( · ) 1 stand for the transpose, Hermitian transpose, conjugate and matrix inversion, respectively. j denotes the imaginary unit with j 2 = 1 . ⊗ stands for the Kronecker product. · and | · | denote the Euclidean norm and absolute value, respectively. diag { · } denotes the diagonalization operator. E { · } represents the mathematical expectation. I N denotes the N × N identity matrix.

2. ST2D Signal Model of GNSS Receivers

Figure 1 shows a general STAP filter structure for GNSS receiver. Consider a uniform linear array composed of M antenna elements, and the spacing between adjacent elements is d. A temporal filter with N taps and the sampling period T s is connected behind each antenna element. The received signal of each antenna element is down-converted to the baseband by the radio frequency (RF) front end for the analog-to-digital conversion. After passing through the RF front end, the digitized signal is processed by a FIR filter. Finally, the output signal of the STAP filter is sent to each channel of the GNSS receivers for signal acquisition, tracking, and navigation solution. Consider that the direction of GNSS signal reaching the antenna array is θ d , and the directions of G interferences are θ i g ( g = 1 , 2 , , G ) . The received signal vector of the antenna array can be expressed as
x 1 ( t ) = d 1 ( t ) + i 1 ( t ) + n 1 ( t ) = a s ( θ d ) s d ( t ) + g = 1 G a s ( θ i g ) s i g ( t ) + n 1 ( t )
where d 1 ( t ) , i 1 ( t ) , and  n 1 ( t ) are respectively the GNSS signal, interference and noise, a s ( θ d ) and a s ( θ i g ) are respectively the steering vectors of the GNSS signal and the g-th interference, s d ( t ) and s i g ( t ) are respectively the complex waveforms of GNSS signal and the g-th interference. The steering vector a s ( θ ) C M × 1 is defined as
a s ( θ ) = [ 1 , e j 2 π d sin θ / λ , , e j 2 π ( M 1 ) d sin θ / λ ] T
where λ is the signal wavelength.
Stack M N time domain snapshots signals into a vector as
x ( t ) = d ( t ) + i ( t ) + n ( t ) = ( I N a s ( θ d ) ) s d + g = 1 G ( I N a s ( θ i g ) ) s i g + n ( t ) = [ x 11 ( t ) , , x M 1 ( t ) , , x 1 N ( t ) , , x M N ( t ) ] T
where s d = [ s d ( t ) , , s d ( t ( N 1 ) T s ) ] T and s i g = [ s i g ( t ) , , s i g ( t ( N 1 ) T s ) ] T are the complex envelope sets of the GNSS signal and the g-th interference, respectively. x m n ( t ) is the time domain sample at the n-th tap of the temporal filter after the m-th element.
Define the weight vector as w = [ w 11 , , w M 1 , , w 1 N , , w M N ] T C M N × 1 , where w m n is the STAP weight at the n-th tap of the temporal filter after the m-th element. The STAP beamformer output can be written as
y = w H x ( t ) = m = 1 M n = 1 N w m n * x m n ( t )
To generate the weight vector, the well-known Power Inversion (PI) criterion is proposed to minimize the output power of the STAP beamformer. In the PI algorithm, an antenna element is selected as the reference element and the remaining antennas are then viewed as auxiliary elements. The weight vector of the PI algorithm is calculated as
w PI = α R x x 1 c
where α = 1 / ( c H R x x 1 c ) with c = [ 1 , 0 , , 0 ] T C M N × 1 , and  R x x is the auto-correlation matrix of x ( t ) which is denoted as R x x = E { x ( t ) x H ( t ) } .

3. Proposed Robust STAP Method for GNSS Receivers in Coherent Signal Environments

In this section, a robust STAP MVDR beamforming method for coherent signals in GNSS receivers is proposed. We firstly develop a ST2D-IAA algorithm and modify it to obtain the power spectrum of 2D grid points and the noise power. Then we reconstruct the STINCM for GNSS receivers by using the estimated power spectrum of the interference angle region and the noise power. Furthermore, the STSV of GNSS signal is corrected. Finally, the weight vector of the STAP beamforming is calculated based on the reconstructed STINCM and the corrected STSV of GNSS signal.

3.1. ST2D-IAA Spectrum Estimation

Considering the STAP beamformer, the ST2D parameters of impinging signals is denoted as
Φ = { ( θ , f ) | θ [ 0 , π / 2 ] , f ( 0 , f h ] }
where θ represents the pitching angle, f represents the frequency, and  f h represents the cutoff frequency of the low-pass filters in GNSS receivers. The whole space domain Θ = { θ | θ [ 0 , π / 2 ] } can be uniformly discretized as Θ ¯ = { θ 1 , θ 2 , , θ K } by K grid points, and the time domain Ω = { f | f ( 0 , f h ] } can be uniformly discretized as Ω ¯ = { f 1 , f 2 , , f L } by L grid points. Accordingly, the joint space-time domain Φ can be discretized as
Φ ¯ = { ( θ k , f l ) | θ k { θ 1 , θ 2 , , θ K } , f { f 1 , f 2 , , f L } }
by K × L 2D grid points. The STSV of the grid point ( θ k , f l ) is represented as
a ( θ k , f l ) = a s ( θ k ) a t ( f l )
for k = 1 , 2 , , K and l = 1 , 2 , , L . The temporal steering vector a t ( f l ) is defined as
a t ( f l ) = [ 1 , e j 2 π f l / f s , , e j 2 π ( N 1 ) f l / f s ] T
where f s represents the sampling frequency. The space-time steering matrix of all 2D grid points in Φ ¯ is denoted as
A = [ a ( θ 1 , f 1 ) , , a ( θ 1 , f L ) , , a ( θ K , f 1 ) , , a ( θ K , f L ) ]
The complex envelope vector of all 2D grid points in Φ ¯ is written as
s ( t ) = [ s 1 , 1 ( t ) , , s 1 , L ( t ) , , s K , 1 ( t ) , , s K , L ( t ) ] T
The power vector of all 2D grid points in Φ ¯ is given by
p = [ p 1 , 1 , , p 1 , L , , p K , 1 , , p K , L ] T
In order to apply the IAA algorithm to the STAP beamforming of GNSS receivers, we propose a ST2D-IAA algorithm to estimate the power spectrum. Similar to the IAA algorithm, we build the weighted least squares cost function for ST2D-IAA algorithm as
t = 1 T x ( t ) s k , l ( t ) a ( θ k , f l ) Q 1 ( k , l ) 2
where x Q 1 ( k , l ) 2 x H Q 1 ( k , l ) x and the matrix Q ( k , l ) is denoted as
Q ( k , l ) = R p k , l a ( θ k , f l ) a H ( θ k , f l )
with
R = A diag { p } A H
p k , l = 1 T t = 1 T | s k , l ( t ) | 2
Minimizing (13) with respect to s k , l ( t ) yields
s k , l ( t ) = a H ( θ k , f l ) Q 1 ( k , l ) x ( t ) a H ( θ k , f l ) Q 1 ( k , l ) a ( θ k , f l )
Using (14) and the matrix inversion lemma (see, e.g., [35]), (17) can be rewritten as
s k , l ( t ) = a H ( θ k , f l ) R 1 x ( t ) a H ( θ k , f l ) R 1 a ( θ k , f l )
The steps of the ST2D-IAA algorithm are shown in the Appendix A.1. It should be noted that this algorithm has to be performed in an iterative way because it requires R which contains the unknown signal powers. The initialization of this algorithm is implemented using a standard delay-and-sum beamformer. Let p ( i ) denote the power vector at the i-th iteration in the ST2D-IAA algorithm. The iteration will be terminated when the relative change p ( i ) p ( i 1 ) is less than a specified tolerance. Extensive simulations have shown that the number of iterations is less than 12 for the ST2D-IAA algorithm to converge. One can obtain an accurate power spectrum of all 2D grid points from the results of ST2D-IAA algorithm. In addition, the presence of coherent sources has little influence on the performance of ST2D-IAA algorithm.
Although the ST2D-IAA algorithm can obtain accurate power spectrum, it cannot estimate the noise power iteratively. In order to estimate the power spectrum and the noise power simultaneously in the iteration, we further modify the ST2D-IAA algorithm. Based on the weighted least squares method, we modify the optimization problem for estimating s k , l ( t ) as below
min s k , l ( t ) t = 1 T x ( t ) s k , l ( t ) a ( θ k , f l ) R ¯ 1 2
where the matrix R ¯ is denoted as R ¯ = A diag { p ¯ } A H + σ ¯ n 2 I M N with the power vector p ¯ = [ p ¯ 1 , 1 , , p ¯ 1 , L , , p ¯ K , 1 , , p ¯ K , L ] T and the noise power σ ¯ n 2 . By iteratively solving (19), s k , l ( t ) in the ( i + 1 ) -th iteration is estimated as
s ¯ k , l ( i + 1 ) ( t ) = a H ( θ k , f l ) ( R ¯ ( i ) ) 1 x ( t ) a H ( θ k , f l ) ( R ¯ ( i ) ) 1 a ( θ k , f l )
where
R ¯ ( i ) = A diag { p ¯ ( i ) } A H + ( σ ¯ n 2 ) ( i ) I M N
p ¯ ( i ) = [ p ¯ 1 , 1 ( i ) , , p ¯ 1 , L ( i ) , , p ¯ K , 1 ( i ) , , p ¯ K , L ( i ) ] T
p ¯ k , l ( i ) = 1 T t = 1 T | s ¯ k , l ( i ) ( t ) | 2
( σ ¯ n 2 ) ( i ) = 1 M N ( x ( t ) φ ¯ ( i ) A s ¯ ( i ) ( t ) 2
s ¯ ( i ) ( t ) = [ s ¯ 1 , 1 ( i ) ( t ) , , s ¯ 1 , L ( i ) ( t ) , , s ¯ K , 1 ( i ) ( t ) , , s ¯ K , L ( i ) ( t ) ] T
where the scaling coefficient φ ¯ ( i ) is obtained from the following optimization problem
min φ ¯ ( i ) x ( t ) φ ¯ ( i ) A s ¯ ( i ) ( t ) 2
The solution of (26) is given by
φ ¯ ( i ) = ( s ¯ ( i ) ( t ) ) H A H x ( t ) ( s ¯ ( i ) ( t ) ) H A H A s ¯ ( i ) ( t )
The criterion for the termination of the iteration process (20)–(27) is given by p ¯ ( i ) p ¯ ( i 1 ) < ς . Through extensive simulations, the modified ST2D-IAA algorithm will converge after at most 12 iterations. The steps of the modified ST2D-IAA algorithm are shown in Appendix A.2.

3.2. STINCM Reconstruction

Since the powers of interferences are much higher than those of GNSS signal and noise, the ST2D parameters associated with the interferences can be easily estimated through Pauta criterion. The ST2D parameter of the q-th interference can be denoted as
Φ g = { ( θ , f ) | θ ( θ i g , θ i g + U θ ) , f ( f i g , f i g + V f ) }
where θ and f denote the spatial and temporal widths occupied by a certain basic 2D subregion, respectively. Φ g corresponds to the area marked in yellow as shown in Figure 2, which contains the following discretized 2D grid points
Φ ¯ g = { ( θ , f ) | θ { θ i g , 1 , , θ i g , U } , f { f i g , 1 , , f i g , V } } = ( θ i g , 1 , f i g , 1 ) ( θ i g , 1 , f i g , 2 ) ( θ i g , 1 , f i g , V ) ( θ i g , 2 , f i g , 1 ) ( θ i g , 2 , f i g , 2 ) ( θ i g , 2 , f i g , V ) ( θ i g , U , f i g , 1 ) ( θ i g , U , f i g , 2 ) ( θ i g , U , f i g , V )
The ST2D steering matrix corresponding to Φ ¯ g can be denoted as
A ( Φ ¯ g ) = [ a ( θ i g , 1 , f i g , 1 ) , , a ( θ i g , 1 , f i g , V ) , , a ( θ i g , U , f i g , 1 ) , , a ( θ i g , U , f i g , V ) ]
Based on the results of modified ST2D-IAA algorithm, the ST2D spectrum estimation vector corresponding to Φ ¯ g can be written as
p ¯ ( Φ ¯ g ) = [ p ¯ ( θ i g , 1 , f i g , 1 ) , , p ¯ ( θ i g , 1 , f i g , V ) , , p ¯ ( θ i g , U , f i g , 1 ) , , p ¯ ( θ i g , U , f i g , V ) ] T
Considering all of the ST2D steering matrices A ( Φ ¯ g ) , g = 1 , , G and the ST2D spectrum estimation vectors p ¯ ( Φ ¯ g ) , g = 1 , , G , we can construct the space-time interference covariance matrix as
R ^ i = g = 1 G A ( Φ ¯ g ) diag { p ¯ ( Φ ¯ g ) } A H ( Φ ¯ g ) = g = 1 G u = 1 U v = 1 V p ¯ ( θ i g , u , f i g , v ) a ( θ i g , u , f i g , v ) a H ( θ i g , u , f i g , v )
The final noise power σ ¯ n 2 can be obtained from the modified ST2D-IAA algorithm. Thus, the noise covariance matrix can be estimated as
R ^ n = σ ¯ n 2 I M N
By using the constructed space-time interference covariance matrix R ^ i and the noise covariance matrix R ^ n , we reconstruct the STINCM as
R ^ i n = R ^ i + R ^ n = g = 1 G u = 1 U v = 1 V p ¯ ( θ i g , u , f i g , v ) a ( θ i g , u , f i g , v ) a H ( θ i g , u , f i g , v ) + σ ¯ n 2 I M N

3.3. STSV Estimation

In order to suppress the interferences with high powers while remaining a distortionless response for the GNSS signal, we formulate the optimization problem of the STAP beamforming based on the MVDR criterion as
min w w H R i n w subject to N is odd w H a = 1 a = [ 0 , , 0 M ( N 1 ) / 2 , a s ( θ d ) , 0 , , 0 M ( N 1 ) / 2 ] T
where R i n = E { ( i ( t ) + n ( t ) ) ( i ( t ) + n ( t ) ) H } C M N × M N represents the theoretical STINCM. The constraints in (35) are imposed to guarantee the linearity of the STAP beamforming response so that no biases are introduced into the code and carrier phase measurements [12]. The optimal solution of (35) is given by
w o p t = R i n 1 a a H R i n 1 a
Since R i n is unavailable in practice, we replace it with the reconstructed STINCM R ^ i n in (34). Besides, the STSV a is not accurate because it may have STSV error caused by many factors in practical circumstances. By taking the STSV error e into account, the actual STSV can be written as a + e . Based on the reconstructed STINCM R ^ i n and the actual STSV a + e , the STAP beamformer output power can be denoted as
P ^ ( a + e ) = 1 ( a + e ) H R ^ i n 1 ( a + e )
The STSV error e can be estimated by maximizing P ^ ( a + e ) or, equivalently, by minimizing the denominator of P ^ ( a + e ) . The STSV error e can be further decomposed into the orthogonal component e and the parallel component e . Since the parallel component e does not have influence on the STAP beamforming performance, we dismiss e and only take e into computation. The optimization problem for estimating e can be formulated as
min e ( a + e ) H R ^ i n 1 ( a + e ) subject to a H e = 0 ( a + e ) H R ^ i n ( a + e ) a H R ^ i n a
Since R ^ i n 1 0 with the symbol ≻ denoting the positive-definite operation, the convex optimization problem (38) can be easily solved to obtain the solution e ^ . With the estimated e ^ , the STSV a can be corrected as
a ^ = a + e ^
From the last constraint in the optimization problem (35), it can be seen that the corrected STSV a ^ should satisfy a ^ = [ 0 , , 0 M ( N 1 ) / 2 , a ^ s ( θ d ) , 0 , 0 M ( N 1 ) / 2 ] T . It means that a ^ s ( θ d ) can be denoted as a ^ s ( θ d ) = a ^ ( M N M + 2 2 : M ( N + 1 ) 2 ) . The optimal weight vector is obtained as
w ^ = R ^ i n 1 a ^ a ^ H R ^ i n 1 a ^
It should be mentioned that the GNSS signal is generated with the direct-sequence spread-spectrum (DSSS) technique, which shows that it lies in the whole frequency band after passing through low-pass filters in GNSS receivers. In other words, the STSV of the GNSS signal does not have the form a ( θ d , f d ) = a s ( θ d ) a t ( f d ) with a certain frequency f d . Therefore, the constraints of the optimization problem (35) have to be revised for the GNSS signal. The revised optimization problem can be obtained from
min w w H R ^ i n w subject to N is odd w H C = f C = a ^ s ( θ d ) 0 0 0 a ^ s ( θ d ) 0 0 0 a ^ s ( θ d ) M N × N f = [ 0 , , 0 ( N 1 ) / 2 , 1 , 0 , , 0 ( N 1 ) / 2 ] T
By solving (41), the final weight vector can be given by
w ¯ = R ^ i n 1 C [ C H R ^ i n 1 C ] 1 f
Based on the weight vector in (42), the array output carrier-to-noise ratio C / N 0 of the GNSS signal can be denoted as
C / N 0 = B σ d 2 | w H C | 2 w H R ^ i n w
where B denotes the signal bandwidth and σ d 2 denotes the power of the GNSS signal.
The work [36] assesses the performance of receivers by measuring a zero baseline and a short baseline in the field. In the short baseline test, the multipath error is compensated by a set of harmonic functions, which can guarantee the property of unbiasedness of the least-squares estimators. However, The least-squares method provided in [36] may be less agile than the proposed method of this paper in terms of iterations.

4. Simulation Results

Consider a uniform linear array consisting of 10 omnidirectional elements with half wavelength spacing and each element is followed by 8-order FIR filter. The desired GNSS signal is assumed as BeiDou-2 (BD2) signal at B3 band whose carrier frequency is 1268.52 MHz and mainlobe bandwidth is 20.46 MHz. The desired BD2 signal impinges on the array from the direction θ d = 5 . Two interferences whose bandwidths are the same as the desired BD2 signal are respectively incident from the directions θ i 1 = 30 and θ i 2 = 30 and the common INR is 30 dB. Only the interference from the direction 30 is coherent with the desired BD2 signal. In coherent signal environment, the number of samples utilized to calculate the weight vector is chosen as 2400.
Firstly, the spatial-temporal beampatterns corresponding to the proposed method are obtained to observe whether nulls are formed in the directions of interferences. Then, the output C / N 0 performance of the proposed method is mainly compared with the sample matrix inverse method in [37], the diagonal loading in [10], the reconstruct-estimate method in [38]. Although the sampling matrix inverse method is widely used to calculate the weighted vector of STAP beamformers, it can neither suppress the interferences in the coherent signal environments nor be robust to the STSV error caused by the model errors. As for the diagonal loading method, it is a widely used method to deal with space-time covariance matrix errors. The reconstruct-estimate method is also an INCM reconstruction method with STAP structure, which is robust to model errors. In order to verify the robustness of the proposed method in the coherent signal environments and model error conditions, the above methods are selected for comparison. In the reconstruct-estimate method, the angle region of the desired BD2 signal is [ θ d 5 , θ d + 5 ] . In the proposed method of this paper, the two interference angle regions are [ θ i 1 10 , θ i 1 + 10 ] and [ θ i 2 10 , θ i 2 + 10 ] , respectively. The threshold ς is set to 10 3 .
Example 1: Spatial-temporal beampatterns. The 3D and 2D spatial-temporal beampatterns of the proposed method are shown in Figure 3a and Figure 3b respectively. It can be seen that while the entire frequency band at the desired BD2 signal direction forms a mainlobe, nullings are formed in the entire frequency bands of the interference directions 30 and 30 , respectively. This shows that the proposed STAP beamforming method can effectively suppress the interferences in the coherent signal environments, because the proposed method uses the modified ST2D-IAA method to estimate the precise power spectrum.
Example 2: Output C / N 0 performance.Figure 4a shows the relationship curves of the sampling matrix inverse method, diagonal loading method, reconstruct-estimate method and the proposed method with the change of the output C / N 0 versus the input SNR. Note that although the GNSS signal input SNR is usually 20 dB, the spot beam has been or will be equipped with modern GNSS in the future, so the corresponding SNR will be increased to 27 dB, so the SNR range of BD2 signal is [ 30 , 10 ] dB. It can be seen from Figure 4a that the output C / N 0 performance of the proposed method is better than those of the other methods, because the proposed method can suppress the interferences in the coherent signal environments, while the other methods cannot deal with the interferences in the coherent signal environments. Figure 4b shows the relationship curves of the sampling matrix inverse method, diagonal loading method, reconstruct-estimate method and the proposed method with the change of the output C / N 0 versus the number of snapshots. It can be seen that the output C / N 0 performance of the proposed method is still better than the other methods, and the proposed method has good advantages in the output C / N 0 .
Example 3: The scenary of direction error. When the direction error of the desired BD2 signal is 3 , that is, the actual direction of the BD2 signal is 8 . In this case, Figure 5a and Figure 5b respectively show the relationship curves of the change of the output C / N 0 versus the input SNR and the number of snapshots, while Figure 5c shows the relationship curves of the output C / N 0 versus the direction error. It can be seen from these figures that both the proposed method and the reconstruct-estimate method are robust to the signal direction error, because both of them consider the signal STSV error and correct the signal STSV. However, the performance of the diagonal loading method and the sampling matrix inverse method is poor in the case of signal direction error, because these two methods do not consider the signal STSV error, and do not correct the signal STSV. In addition, the proposed method has the best performance of the output C / N 0 , which is because the proposed method not only reconstructs the STINCM, but also corrects the signal STSV.
Example 4: The case of arbitrary STSV error. In the case of arbitrary STSV error of the signal, the actual STSV of the desired BD2 signal can be expressed as a = a ( θ d + δ θ ) + δ a , where the direction error δ θ is fixed at 3 , the STSV error δ a is a random vector with the mean 0 and the variance σ 2 I M N × 1 . The parameter σ 2 is set to 0.64 . Figure 6a and Figure 6b respectively show the relationship curves of output C / N 0 versus the input SNR and versus the number of snapshots under the condition of arbitrary STSV error of the signal. It can be seen that the output C / N 0 of the proposed method is higher than that of the other methods, because the proposed method combines the accurate STINCM and the corrected STSV, which indicates that the proposed method is robust to arbitrary STSV error of signal.
Example 5: The situation of small number of snapshots. We investigate the performance of the tested methods when the number of snapshots is small. The simulated curves of output C / N 0 versus the number of snapshots are displayed in Figure 7, where the number of snapshots ranges from 20 to 200. Obviously, the proposed method has the minimal performance loss, and performs better than the other methods in the situation of small number of snapshots. The reason lies in that the ST2D-IAA spectrum estimation in the proposed method can work with the limited sampling snapshots. As described in [33], the IAA spectrum estimation can even handle the single snapshot in the coprime virtual uniform linear array. Therefore, the performance of the proposed method is robust against the small number of snapshots.

5. Conclusions

In this paper, a robust STAP method based on the modified ST2D-IAA spectrum estimation is proposed to suppress interferences for GNSS receivers in coherent signal environments. By applying the IAA algorithm to the ST2D signal model in GNSS receivers and further modifying the ST2D-IAA algorithm, the power spectrum and noise power can be estimated accurately at the same time. The STINCM of GNSS receivers is reconstructed by selecting the estimated power spectrum and noise power corresponding to the interference angle regions. In order to obtain the error vector of the STSV, the STAP beamforming optimization problem based on the MVDR criterion is constructed to correct the STSV of GNSS signal. The weight vector of STAP beamforming is calculated by combining the reconstructed STINCM and the corrected GNSS signal STSV. Simulation results show that the proposed method can suppress the interferences in the coherent signal environments, and has better robustness than the existing methods under the conditions of model errors. The proposed method has a potential prospect in maintaining the integrity and functionality of GNSS systems against coherent interference in many applications.

Author Contributions

Conceptualization, Z.M. and F.S.; methodology, Z.M.; software, F.S.; validation, F.S.; formal analysis, Z.M.; investigation, Z.M.; resources, F.S.; data curation, F.S.; writing—original draft preparation, Z.M.; writing—review and editing, Z.M.; visualization, Z.M.; supervision, F.S.; project administration, Z.M. and F.S.; funding acquisition, Z.M. and F.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the National Natural Science Foundation of China under Grant 62203447 and Grant 61673128; and in part by the Fundamental Research Funds for the Central Universities under Grant 2022QN1049.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1

The ST2D-IAA algorithm is summarized in Algorithm A1.
Algorithm A1 The ST2D-IAA Algorithm
1:   initialization
2:      for  k = 1 , 2 , , K
3:         for  l = 1 , 2 , , L
4:             p k , l = t = 1 T | a H ( θ k , f l ) x ( t ) | 2 ( a H ( θ k , f l ) a ( θ k , f l ) ) 2 T
5:         end for
6:      end for
7:   repeat
8:       R = A diag { p } A H
9:      for  k = 1 , 2 , , K
10:         for  l = 1 , 2 , , L
11:             s k , l ( t ) = a H ( θ k , f l ) R 1 x ( t ) a H ( θ k , f l ) R 1 a ( θ k , f l ) , t = 1 , 2 , , T
12:             p k , l = 1 T t = 1 T | s k , l ( t ) | 2
13:         end for
14:      end for
15:   until (convergence)

Appendix A.2

The modified ST2D-IAA algorithm is summarized in Algorithm A2.
Algorithm A2 The Modified ST2D-IAA Algorithm
1:   initialization
2:       p ¯ k , l ( 0 ) = t = 1 T | a H ( θ k , f l ) x ( t ) | 2 ( a H ( θ k , f l ) a ( θ k , f l ) ) 2 T k = 1 , , K ;   l = 1 , , L
3:       i = 0
4:   repeat
5:       p ¯ ( i ) = [ p ¯ 1 , 1 ( i ) , , p ¯ 1 , L ( i ) , , p ¯ K , 1 ( i ) , , p ¯ K , L ( i ) ] T
6:       R ¯ ( i ) = A diag { p ¯ ( i ) } A H + ( σ ¯ n 2 ) ( i ) I M N
7:      for  t = 1 , 2 , , T
8:         for  k = 1 , 2 , , K
9:            for  l = 1 , 2 , , L
10:                s ¯ k , l ( i + 1 ) ( t ) = a H ( θ k , f l ) ( R ¯ ( i ) ) 1 x ( t ) a H ( θ k , f l ) ( R ¯ ( i ) ) 1 a ( θ k , f l )
11:            end for
12:         end for
13:          s ¯ ( i + 1 ) ( t ) = [ s ¯ 1 , 1 ( i + 1 ) ( t ) , , s ¯ K , L ( i + 1 ) ( t ) ] T
14:          φ ¯ ( i + 1 ) = ( s ¯ ( i + 1 ) ( t ) ) H A H x ( t ) ( s ¯ ( i + 1 ) ( t ) ) H A H A s ¯ ( i + 1 ) ( t )
15:          ( σ ¯ n 2 ) ( i + 1 ) = 1 M N ( x ( t ) φ ¯ ( i + 1 ) A s ¯ ( i + 1 ) ( t ) ) 2
16:      end for
17:       p ¯ k , l ( i + 1 ) = 1 T t = 1 T | s ¯ k , l ( i + 1 ) ( t ) | 2 k = 1 , , K ;   l = 1 , , L
18:       i = i + 1
19:   until (convergence)

References

  1. Marut, G.; Hadas, T.; Kaplon, J.; Trzcina, E.; Rohm, W. Monitoring the water vapor content at high spatio-temporal resolution using a network of low-cost multi-GNSS receivers. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5804614. [Google Scholar] [CrossRef]
  2. Suzuki, T. GNSS Odometry: Precise trajectory estimation based on carrier phase cycle slip estimation. IEEE Robot. Autom. Lett. 2022, 7, 7319–7326. [Google Scholar] [CrossRef]
  3. Osechas, O.; Fohlmeister, F.; Dautermann, T.; Felux, M. Impact of GNSS-band radio interference on operational avionics. Navig. J. Inst. Navig. 2022, 69, navi.516. [Google Scholar] [CrossRef]
  4. Chen, X.; He, D.; Yan, X.; Yu, W.; Truong, T.K. GNSS interference type recognition with fingerprint spectrum DNN method. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 4745–4760. [Google Scholar] [CrossRef]
  5. Meng, L.; Yang, L.; Yang, W.; Zhang, L. A survey of GNSS spoofing and anti-spoofing technology. Remote Sens. 2022, 14, 4826. [Google Scholar] [CrossRef]
  6. Huang, L.; Lu, Z.; Xiao, Z.; Ren, C.; Song, J.; Li, B. Suppression of jammer multipath in GNSS antenna array receiver. Remote Sens. 2022, 14, 350. [Google Scholar] [CrossRef]
  7. Lu, Z.; Song, J.; Huang, L.; Ren, C.; Xiao, Z.; Li, B. Distortionless 1/2 overlap windowing in frequency domain anti-jamming of satellite navigation receivers. Remote Sens. 2022, 14, 1801. [Google Scholar] [CrossRef]
  8. Zhang, J.; Feng, W.; Yuan, T.; Wang, J.; Sangaiah, A.K. SCSTCF: Spatial-channel selection and temporal regularized correlation filters for visual tracking. Appl. Soft Comput. 2022, 118, 108485. [Google Scholar] [CrossRef]
  9. Zhou, M.; Wang, Q.; He, F.; Meng, J. Impacts of phase noise on the anti-jamming performance of power inversion algorithm. Sensors 2022, 22, 2362. [Google Scholar] [CrossRef]
  10. Cox, H.; Zeskind, R.M.; Owen, M.M. Robust adaptive beamforming. IEEE Trans. Acoust. Speech Signal Process. 1987, 35, 1365–1376. [Google Scholar] [CrossRef]
  11. Zhang, W.; Liu, T.; Yang, G.; Jiang, C.; Hu, Y.; Lan, T.; Zhao, Z. A novel method for improving quality of oblique incidence sounding ionograms based on eigenspace-based beamforming technology and Capon high-resolution range profile. Remote Sens. 2022, 14, 4305. [Google Scholar] [CrossRef]
  12. Huang, Y.; Fu, H.; Vorobyov, S.A.; Luo, Z.Q. Robust adaptive beamforming via worst-case sinr maximization with nonconvex uncertainty sets. IEEE Trans. Signal Process. 2023, 71, 218–232. [Google Scholar] [CrossRef]
  13. Gu, Y.; Leshem, A. Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation. IEEE Trans. Signal Process. 2012, 60, 3881–3884. [Google Scholar]
  14. Meng, Z.; Dong, S.; Shi, X.; Wang, X. Robust beamforming for non-circular signals in uniform linear arrays with unknown mutual couplin. Digit. Signal Process. 2022, 122, 103378. [Google Scholar] [CrossRef]
  15. Li, H.; Geng, J.; Xie, J. Robust adaptive beamforming based on covariance matrix reconstruction with RCB principle. Digit. Signal Process. 2022, 127, 103565. [Google Scholar] [CrossRef]
  16. Meng, Z.; Zhou, W. Robust adaptive beamforming for coprime array with steering vector estimation and covariance matrix reconstruction. IET Commun. 2020, 14, 2749–2758. [Google Scholar] [CrossRef]
  17. Xie, Y.; Chen, F.; Huang, L.; Liu, Z.; Wang, F. Carrier phase bias correction for GNSS space-time array processing using time-delay data. GPS Solut. 2023, 27, 113. [Google Scholar] [CrossRef]
  18. Hao, F.; Li, X.; Wang, W.; Zhao, J. A STAP anti-interference technology with zero phase bias in wireless IoT systems based on high-precision positioning. Front. Phys. 2023, 11, 284. [Google Scholar] [CrossRef]
  19. Stenberg, N.; Axell, E.; Rantakokko, J.; Hendeby, G. Results on GNSS spoofing mitigation using multiple receivers. Navig. J. Inst. Navig. 2022, 69, navi.510. [Google Scholar] [CrossRef]
  20. Wang, Y.; Liu, W.; Huang, L.; Xiao, Z.; Wang, F. Distortionless pseudo-code tracking space-time adaptive processor based on the PI criterion for GNSS receiver. IET Radar Sonar Navig. 2020, 14, 1984–1990. [Google Scholar] [CrossRef]
  21. Brachvogel, M.; Niestroj, M.; Meurer, M.; Hasnain, S.N.; Stephan, R.; Hein, M.A. Space-time adaptive processing as a solution for mitigating interference using spatially-distributed antenna arrays. Navig. J. Inst. Navig. 2023, 70, navi.592. [Google Scholar] [CrossRef]
  22. Chen, X.; Morton, Y. Iterative subspace alternating projection method for GNSS multipath DOA estimation. IET Radar Sonar Navig. 2016, 10, 1260–1269. [Google Scholar] [CrossRef]
  23. Ghiasi, Y.; Duguay, C.R.; Murfitt, J.; van der Sanden, J.J.; Thompson, A.; Drouin, H.; Prévost, C. Application of GNSS interferometric reflectometry for the estimation of lake ice thickness. Remote Sens. 2020, 12, 2721. [Google Scholar] [CrossRef]
  24. Wen, F.; Shi, J.; Zhang, Z. Generalized spatial smoothing in bistatic EMVS-MIMO radar. Signal Process. 2022, 193, 108406. [Google Scholar] [CrossRef]
  25. Shi, J.; Yang, Z.; Liu, Y. On parameter identifiability of diversity-smoothing-based MIMO radar. IEEE Trans. Aerosp. Electron. Syst. 2021, 58, 1660–1675. [Google Scholar] [CrossRef]
  26. Pote, R.R.; Rao, B.D. Maximum likelihood-based gridless doa estimation using structured covariance matrix recovery and sbl with grid refinement. IEEE Trans. Signal Process. 2023, 71, 802–815. [Google Scholar] [CrossRef]
  27. Rahamim, D.; Tabrikian, J.; Shavit, R. Source localization using vector sensor array in a multipath environment. IEEE Trans. Signal Process. 2004, 52, 3096–3103. [Google Scholar] [CrossRef]
  28. He, J.; Jiang, S.; Wang, J.; Liu, Z. Polarization difference smoothing for direction finding of coherent signals. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 469–480. [Google Scholar] [CrossRef]
  29. Zhang, Z.; Wen, F.; Shi, J.; He, J.; Truong, T.K. 2D-DOA estimation for coherent signals via a polarized uniform rectangular array. IEEE Signal Process. Lett. 2023, 30, 893–897. [Google Scholar] [CrossRef]
  30. Park, Y.; Gerstoft, P. Gridless sparse covariance-based beamforming via alternating projections including co-prime arrays. J. Acoust. Soc. Am. 2022, 151, 3828–3837. [Google Scholar] [CrossRef] [PubMed]
  31. Mao, D.; Zhang, Y.; Zhang, Y.; Huo, W.; Pei, J.; Huang, Y. Target fast reconstruction of real aperture radar using data extrapolation-based parallel iterative adaptive approach. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 2258–2269. [Google Scholar] [CrossRef]
  32. Wang, B.; Li, S.; Battistelli, G.; Chisci, L. Fast iterative adaptive approach for indoor localization with distributed 5G small cells. IEEE Wirel. Commun. Lett. 2022, 11, 1980–1984. [Google Scholar] [CrossRef]
  33. Meng, Z.; Zhou, W. Virtual uniform linear iterative adaptive approach for robust adaptive beamforming. In Proceedings of the 2020 39th Chinese Control Conference (CCC), Shenyang, China, 27–29 July 2020; pp. 2980–2985. [Google Scholar]
  34. Yang, X.; Xie, J.; Li, H.; He, Z. Robust adaptive beamforming of coherent signals in the presence of the unknown mutual coupling. IET Commun. 2018, 12, 75–81. [Google Scholar] [CrossRef]
  35. Stoica, P.; Moses, R.L. Spectral Analysis of Signals; Prentice Hall: Upper Saddle River, NJ, USA, 2005. [Google Scholar]
  36. Amiri-Simkooei, A.R.; Tiberius, C.C.J.M. Assessing receiver noise using GPS short baseline time series. GPS Solut. 2007, 11, 21–35. [Google Scholar] [CrossRef]
  37. Capon, J. High-resolution frequency-wavenumber spectrum analysis. Proc. IEEE 1969, 57, 1408–1418. [Google Scholar] [CrossRef]
  38. Li, Z.; Zhang, Y.; Liu, H. A robust STAP method for airborne radar based on clutter covariance matrix reconstruction and steering vector estimation. Digit. Signal Process. 2018, 78, 82–91. [Google Scholar] [CrossRef]
Figure 1. General STAP filter structure for GNSS receiver.
Figure 1. General STAP filter structure for GNSS receiver.
Remotesensing 15 04212 g001
Figure 2. Division of the joint space-time domain for the STAP architecture.
Figure 2. Division of the joint space-time domain for the STAP architecture.
Remotesensing 15 04212 g002
Figure 3. The spatial-temporal beampatterns of the proposed method. (a) 3D spatial-temporal beampattern. (b) 2D spatial-temporal beampattern.
Figure 3. The spatial-temporal beampatterns of the proposed method. (a) 3D spatial-temporal beampattern. (b) 2D spatial-temporal beampattern.
Remotesensing 15 04212 g003
Figure 4. The simulation results in example 2. (a) output C / N 0 versus the input SNR. (b) output C / N 0 versus the number of snapshots.
Figure 4. The simulation results in example 2. (a) output C / N 0 versus the input SNR. (b) output C / N 0 versus the number of snapshots.
Remotesensing 15 04212 g004
Figure 5. The simulation results in example 3. (a) output C / N 0 versus the input SNR. (b) output C / N 0 versus the number of snapshots. (c) output C / N 0 versus the direction error.
Figure 5. The simulation results in example 3. (a) output C / N 0 versus the input SNR. (b) output C / N 0 versus the number of snapshots. (c) output C / N 0 versus the direction error.
Remotesensing 15 04212 g005
Figure 6. The simulation results in example 4. (a) output C / N 0 versus the input SNR. (b) output C / N 0 versus the number of snapshots.
Figure 6. The simulation results in example 4. (a) output C / N 0 versus the input SNR. (b) output C / N 0 versus the number of snapshots.
Remotesensing 15 04212 g006
Figure 7. The simulation results in example 5. Output C / N 0 versus the number of snapshots.
Figure 7. The simulation results in example 5. Output C / N 0 versus the number of snapshots.
Remotesensing 15 04212 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Meng, Z.; Shen, F. Robust Space-Time Adaptive Processing Method for GNSS Receivers in Coherent Signal Environments. Remote Sens. 2023, 15, 4212. https://doi.org/10.3390/rs15174212

AMA Style

Meng Z, Shen F. Robust Space-Time Adaptive Processing Method for GNSS Receivers in Coherent Signal Environments. Remote Sensing. 2023; 15(17):4212. https://doi.org/10.3390/rs15174212

Chicago/Turabian Style

Meng, Zhen, and Feng Shen. 2023. "Robust Space-Time Adaptive Processing Method for GNSS Receivers in Coherent Signal Environments" Remote Sensing 15, no. 17: 4212. https://doi.org/10.3390/rs15174212

APA Style

Meng, Z., & Shen, F. (2023). Robust Space-Time Adaptive Processing Method for GNSS Receivers in Coherent Signal Environments. Remote Sensing, 15(17), 4212. https://doi.org/10.3390/rs15174212

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