The most effective method to overcome GPS jamming problem is to use an adaptive array antenna which has the capability of beamforming or nulling to a certain direction. In this paper, Space Time Adaptive Processing (STAP) beamforming algorithm of four elements array antenna will be designed and the anti-jamming performance will be analyzed. According to the analysis, the signal to noise ratio (SNR) and anti-jamming performance can be enhanced by beamforming algorithm. Also, the time tap effect of STAP algorithm will be analyzed theoretically and verified with array antenna modeling and simulation. Specially, the cautious selection of reference time tap in STAP can prevent the degradation of position accuracy performance.
GPS receivers use very weak satellite signals of 1 pW. Accordingly, they cannot track satellite signals even at weak interference signals, and this has been a major problem (
Parkinson et al. 1996
Kaplan & Hegarty 2006
Rhee & Seo 2013
). For example, 1 Watt portable jammer can interrupt the normal operation of GPS receivers within a 50 km radius. The most effective method for overcoming this problem is to spatially eliminate interference signals using array antennas (
Johnson & Dudgeon 1993
Fante & Vaccaro 2000
Haefner et al. 2003
Kaplan & Hegarty 2006
An array antenna is a system that can control gains depending on incident directions by generating final output through multiplying a number of input signals by a weight vector. If this is applied to GPS anti-jamming devices, the gain in the satellite direction can be conserved while minimizing the gain in the incident direction of a jammer. The performance of an array antenna is basically determined by the number of antenna elements and the arrangement type. In theory, it is known that when the number of antenna elements is N, the number of jammers, to which anti-jamming can be applied, is N-1. Also, there are the arrangement types such as Uniform Linear Array, Uniform Rectangular Array, Uniform Circular Array, and Single Ring Cylindrical Array. In this study, a four-element circular array antenna system, which can process a maximum of three broadband jamming signals, is considered in order to minimize the size, as shown in
. In addition, the distance between antenna elements is a major design factor, and the half wavelength distance principle, which is commonly applied, was used considering mutual coupling and nulling regions (
Johnson & Dudgeon 1993
Digital anti-jamming array antenna system.
In the previous papers, an Space Time Adaptive Processing (STAP) algorithm for array antenna systems, which was extended to time and space domains to extend narrowband anti-jamming number, has been used. In this regard, a pseudorange tracking error, which determines the positioning accuracy of a GPS receiver, occurs. But the cause of this error has not been analyzed (
Carlson et al. 2003
Fante et al. 2004
). In this study, the cause of pseudorange tracking error during the use of an STAP algorithm was analyzed, and the performance analysis of the STAP algorithm which a beamforming algorithm has been added was carried out in terms of signal to noise ratio and changes in the characteristics of satellite signals.
The contents of this paper are as follows. In Chapter 2, an STAP beamforming algorithm was derived, and in Chapter 3 a tap selection technique for preventing pseudorange distortion was proposed. In Chapter 4, the results of the simulation were presented, and conclusions were included in Chapter 5.
2. STAP BEAMFORMING ALGORITHM
- 2.1 Optimal Weight Vector Algorithm For Array Antenna
An array antenna system generates final output by multiplying a number of input signals by a weight vector and adding them together. In this regard, the algorithm for weight vector determination is a constrained optimization problem for the minimization of the objective function,
(w), depending on the constraint condition of the variable, w, as shown in Eq. (1) (
Johnson & Dudgeon 1993
The solution of this constrained optimization problem is one of the stationary points of the Lagrangian, and can be obtained using a Lagrangian multiplier (
Applebaum & Chapman 1976
). First, the Lagrangian of Eq. (1) can be expressed as Eq. (2), and the Lagrangian multiplier,
, of the stationary points of this can be expressed as Eq. (3).
The input value measured at the I/Q channel of the four- element array antenna is defined as x=[x
, and the constraint condition of the weight vector, w, is defined as shown in Eq. (4)
represents the constraint matrix, and c represents the constraining column vector. As the output of the array antenna is
x, the expected value of the output power can be expressed as Eq. (5)
[·] represents the expected value,
], and x
represents the conjugate transpose. Using the variables defined above, the constrained optimization problem for the minimization of the output power can be expressed as Eq. (6).
To obtain the solution of Eq. (6), a Lagrangian multiplier is calculated. Then the weight vector, w
, and the power,
, which are the optimal solution, can be calculated as shown in Eqs. (7) and (8).
The anti-jamming algorithm of the array antenna basically follows the above equations for deriving optimal solution, and the implementation methods are broadly divided into two categories. In the first method, the correlation matrix,
, of an input value is estimated, and the inverse matrix is obtained. Then the weight vector is directly calculated as shown in Eq. (7). In the second method, the weight vector is estimated so that output can be minimized, using a Least Mean Square algorithm without the estimation of the correlation martix,
Lorenzo etal. 2005
- 2.2 STAP Beamforming Algorithm
The STAP input shown in Eq. (9) is used for STAP Algorithm that is extended to time and space domains, which is extended to the time domain by the number of taps T when the input, x=[x
, of the four-element array antenna is used. The remaining processes for calculating the optimal weight vector are identical to the method explained above. For an array antenna with more than four elements, extension can be performed in a similar way.
The degree of freedom of this STAP algorithm is basically obtained by subtracting the number of constraint conditions from the number of free weight vector elements. It is expressed as the following equation, where N-1, which is the degree of freedom of the spatial nulling algorithm that performs spatial filtering, is extended by
(number of taps) times.
For the above STAP algorithm, the degree of freedom is extended by T times. Thus, it is robust to a narrowband jammer, which consumes one degree of freedom, and can counteract
-1narrowband jammers in theory (
Gupta & Ksienski 1983
Berefelt et al. 2003
Vo et al. 2007
The STAP beamforming algorithm uses a Maximum Signal to Interference Ratio method, which calculates the weight vector considering the steering vector, s, in order to conserve the gain and phase in a specific direction (satellite direction) or a specific frequency in Eq. (4) for STAP beamforming algorithm increases the gain in a specif
, is identical to the steering vector,
, and the final optimal weight vector can be expressed as Eq. (11).
The STAP beamforming algorithm increases the gain in a specific direction, and the theoretical gain can be expressed as Eq. (12). In addition, to increase the effect at a specific frequency, thermal noise or jamming signals can be additionally eliminated by applying constrained conditions to a specific tap.
For the above STAP beamforming algorithm, the gain in a specific direction can be achieved by controlling the parts of the steering vector, s, which are related with the array antenna. In addition, it needs to be designed considering the effect of a tap, which is a time delay, and this is analyzed in Chapter 3. Using the STAT beamforming algorithm considering this, a simulation is performed, and the results are presented in Chapter 4.
3. TAP SELECTION IN STAT BEAMFORMING
In this chapter, the distortion characteristics of a pseudorange, which affects the positioning error of a GPS receiver when STAP with an additional filter is applied in the time domain, were analyzed. The tap of STAP is a time delay, and a GPS receiver is very sensitive to a time delay. Especially, for the application of STAP, the effects could be significant depending on the selection of a reference tap. For example, when a satellite and a jammer are arranged as shown in
, and the first tap in
is selected as the reference tap, the front and rear shapes of the Cross Correlation Function (CCF), which is the correlation function of satellite signals, are deformed and distorted as shown in
. In this regard, when the general Early minus Late pseudorange tracking technique was used, a pseudorange error of about 30 m was obtained as shown in
. This error component is problematic because each satellite has a different pseudorange error.
Gain pattern of STAP when SV and jammer exist.
Space Time Adaptive processing structure with taps and 4 antennas.
Distorted CCF of SV signal.
The CCF distortion of satellite signals is due to the nonlinear phase effect of a transversal filter that has time delay taps. The linear phase of the filter represents the case when the group delay, G, of the filter is a constant as shown in Eq. (13).
represents the frequjency variation, and Δ
represents the phase variation.
In the case of a nonlinear phase, each input frequency has a different time delay because a phase variation for a frequency variation is not linearly maintained. Therefore, within a satellite signal bandwidth, each time delay is different depending on the frequency, and the distortion of CCF occurs when correlation functions using replica signals are calculated. To verify this, CCF was obtained by inputting only satellite signals and passing them through a 41st-order transversal filter. When the signals were passed through a nonlinear phase filter as shown in
, the front and rear shapes of CCF were non-symmetric, which resulted in a pseudorange error. In contrast, when the signals were passed through a linear phase filter as shown in
, a time delay occurred, but the front and rear shapes of CCF were symmetric with the same power, which did not result in a pseudorange error when the Early minus Late tracking technique was used.
Nonlinear phase 41 taps FIR filter bode plot (blue: amplitude, red: phase) and non-symmetric CCF.
Linear phase 41 taps FIR filter bode plot (blue: amplitude, red: phase) and symmetric CCF.
From the above discussion, it was found that for a transversal filter with time taps such as the case of an STAP algorithm, the filter should have linear phase characteristics. For an FIR filter with time delay taps to have linear phase characteristics, the filter coefficient should consist of symmetry (
). This also applies to the STAP algorithm, and it indicates that linear phase characteristics can be obtained when the center tap is selected in the process of obtaining a weight vector using a correlation matrix, as shown in the following equations. In other words, the weight vector, w
, is to select a specific column of
, and only the center column of
has symmetry characteristics.
where s represents the vector that located the satellite direction steering elements,
, at the center tap.
shows the CCF correlation function of the satellite obtained after eliminating the jamming signals when the weight vector that selected the center tap was used so that the STAP algorithm could have linear phase characteristics. As shown in the figure, the Early and Late correlation values were symmetric, and thus pseudorange distortion did not occur. On the other hand, there was some time delay, but this is not problematic because it is common to every satellite. In conclusion, the number of taps for the STAP algorithm should be an odd number, and satellite navigation errors can be minimized when the center tap is selected as the reference tap.
CCF of Linear phase STAP filter.
In this chapter, a simulation was performed by locating jamming signals at an azimuth angle of 0˚ and an elevation angle of 10˚, and by locating satellite signals at an azimuth angle of 180˚ and an elevation angle of 50˚ to examine STAT beamforming effects under the condition summarized in
show the gain pattern and phase pattern when a simple power minimization technique without beamforming was applied. Jamming signals were completely eliminated, but a phase shift of -60˚ occurred at an azimuth angle of 180˚ and an elevation angle of 50˚, which was the satellite direction. On the other hand, the satellite signal SNR was calculated to be about 16 dB as shown in
. This was determined by the gain pattern shown in
, and was similar to the input signal SNR (15 dB) of a single antenna in
. In other words, the satellite signal SNR was normal, but the phase value changed depending on the incident direction of satellite signals as shown in the phase pattern. Thus, it would be problematic for systems using a carrier phase.
Gain pattern with anti-jamming of power minimization.
Phase pattern with anti-jamming of power minimization.
SNR with anti-jamming of power minimization.
show the gain pattern and phase pattern for the STAP technique using beamforming. For the gain pattern in
, the gain in the satellite direction was about 6 dB. From this figure, it was found that the maximum gain was maintained in the satellite direction. The phase pattern was calculated as shown in
, and it was maintained at 0˚ in the satellite direction. If beamforming is performed in a different satellite direction, the phase is maintained at 0˚ in the direction, and a phase shift does not occur. The satellite signal SNR was calculated to be 21 dB as shown in
. The value was 5 dB higher than that of the power minimization technique mentioned above, and was 6 dB higher than the input signal of a single antenna in
. This was identical to the gain (6 dB) of the four-element array antenna calculated using Eq. (12). Thus, it was also consistent with the theoretical value.
Gain pattern with anti-jamming of beamforming.
Phase pattern with anti-jamming of beamforming.
SNR with anti-jamming of beamforming.
As shown above, the STAP beamforming technique, which forms a beam in the satellite direction, can prevent a phase shift as well as can enhance satellite signals. This technique conserves signals in the satellite direction, and enhances SNR. Therefore, it is an efficient technique which can improve anti-jamming performance.
On the other hand, the CCF was symmetric as shown in
, which were the results of the STAP algorithm that selected the center tap. Therefore, a pseudorange error did not occur similar to the results of Chapter 3.
In this study, an STAP beamforming algorithm for an array antenna system, which had been extended to time and space domains to increase the narrowband anti-jamming number, was analyzed. The cause of a pseudorange tracking error during the use of the STAP algorithm was theoretically analyzed, and a tap selection technique for preventing this was proposed. In addition, the performance analysis of the STAP algorithm, to which a beamforming algorithm had been added, was carried out in terms of SNR and phase measurements. It was found that the SNR was improved by the beamforming algorithm, and this could enhance antijamming performance. Also, the beamforming algorithm can prevent the shift of carrier phase measurements. Therefore, it could be applied to a reference position survey system or a precise navigation system.
Kihoon Lee is a senior researcher at Agency for Defense Development. He received his B.S. from the Mechanical Engineering Department of POSTECH in 1999. He received his M.S. from the Mechanical Engineering Department of Korean Advanced Institute of Science and Technology in 2001. He has served as a researcher at Agency for Defense Development since 2001. His research focuses on the development of GNSS receiver, Anti-Jamming and Space Based Augmentation System.
Hyoungmin So is a senior researcher of Agency for Defense Development (ADD) in Korea, Republic of. He received B.S. degree in mechanical engineering from Korea Univ. and M.S. and Ph.D. degree in aerospace engineering from Seoul National University (SNU). He worked in the field of GNSS and pseudolite receiver development including SDR and vector tracking loop algorithm in SNU GNSS laboratory. Since 2011, he’s been working for ADD
Kiwon Song is a principal researcher at Agency for Defense Development. He received his B.S. from the Electronical Engineering Department of Seonggyunkwan University in 1984. He received his M.S., Ph.D. degrees from the Electronical Engineering Department of Choongnam University in 1992 and 2002, respectively. He has served as a researcher at Agency for Defense Development since 1984. His research focuses on the satellite navigation system and GNSS signal processing.
Applebaum S. P.
Chapman D. J.
Adaptive Arrays with Main Beam Constraints
IEEE Transactions on Antennas and Propagation
DOI : 10.1109/TAP.1976.1141416
Digital Processing of Signals Theory and Practice
INS/GPS Integration with Adaptive Beamforming
ION GPS/GNSS 2003
Carlson S. G.
Popeck C. A.
Stockmaster M. H.
McDowell C. E
Rockwell Collins’ Flexible Digital Anti-Jam Architecture
ION GPS/GNSS 2003
Fante R. L.
Fitzgibbons M. P.
McDonald K. F.
Effect of Adaptive Array Processing on GPS Signal Crosscorrelation
ION GNSS 17th International Technical Meeting of the Satellite Division
Fante R. L.
Vaccaro J. J.
Wideband Cancellation of Interference in a GPS Receive Array
IEEE Transactions on Aerospace and Electronic systems
DOI : 10.1109/7.845241
Frost O. L.
An Algorithm for Linearly Constrained Adaptive Array Processing
Proceedings of the IEEE
DOI : 10.1109/PROC.1972.8817
Effect of mutual coupling on the performance of adaptive arrays
IEEE Transactions on Antennas and Propagation
DOI : 10.1109/TAP.1983.1143128
G-STAR, Lockheed Martin’s Advanced GPS Anti-Jam Technology
ION 59th Annual Meeting
Adaptive Filter Theory
Johnson D. H.
Dudgeon D. E.
Array Signal Processing Concepts and Techniques
Kaplan E. D.
Hegarty C. J.
Understanding GPS Principles and Applications
Lorenzo D. S.
Adaptive Array Processing for GPS Interference Rejection
Moore T. D.
Analytic study of space-time and spacefrequency adaptive processing for radio frequency interference suppression, PhD Dissertation
The Ohio State University
Parkinson B. W.
The Global Positioning System: Theory and Applications
Rhee J. H.
eLoran Signal Strength and Atmospheric Noise Simulation over Korea
DOI : 10.11003/JKGS.2013.2.2.101
Falchetti C. R.
Morrison A. W.
ADAP: Enhancing GPS Protection for NAVWAR
Proceedings of the 2007 National Technical Meeting of The Institute of Navigation