Navigation Accuracy and Interference Rejection for an Adaptive GPS Antenna Array

Similar documents
Phase Effects Analysis of Patch Antenna CRPAs for JPALS

Mitigation of Continuous and Pulsed Radio Interference with GNSS Antenna Arrays

Miniaturized GPS Antenna Array Technology and Predicted Anti-Jam Performance

Mutual Coupling Estimation for GPS Antenna Arrays in the Presence of Multipath

KINEMATIC TEST RESULTS OF A MINIATURIZED GPS ANTENNA ARRAY WITH DIGITAL BEAMSTEERING ELECTRONICS

TEST RESULTS OF A HIGH GAIN ADVANCED GPS RECEIVER

Test Results of a 7-Element Small Controlled Reception Pattern Antenna

Simulation and Validation of a GPS Antenna Array Concept for JPALS Application

Real-Time Software Receiver Using Massively Parallel

Using GPS to Synthesize A Large Antenna Aperture When The Elements Are Mobile

Small Controlled Reception Pattern Antenna (S-CRPA) Design and Test Results

HIGH GAIN ADVANCED GPS RECEIVER

Galileo GIOVE-A Broadcast E5 Codes and their Application to Acquisition and Tracking

Test Results from a Digital P(Y) Code Beamsteering Receiver for Multipath Minimization Alison Brown and Neil Gerein, NAVSYS Corporation

Phase Center Calibration and Multipath Test Results of a Digital Beam-Steered Antenna Array

Adaptive Array Technology for Navigation in Challenging Signal Environments

Adaptive Antennas. Randy L. Haupt

Null-steering GPS dual-polarised antenna arrays

INTRODUCTION HARDWARE ARCHITECTURE OF RECEIVER

Ionosphere Effects for Wideband GNSS Signals

Towards a Practical Single Element Null Steering Antenna

TEST RESULTS OF A DIGITAL BEAMFORMING GPS RECEIVER FOR MOBILE APPLICATIONS

INTERFERENCE REJECTION OF ADAPTIVE ARRAY ANTENNAS BY USING LMS AND SMI ALGORITHMS

CHAPTER 10 CONCLUSIONS AND FUTURE WORK 10.1 Conclusions

A GPS RECEIVER DESIGNED FOR CARRIER-PHASE TIME TRANSFER

Navigation für herausfordernde Anwendungen Robuste Satellitennavigation für sicherheitskritische Anwendungen

Evaluation of C/N 0 estimators performance for GNSS receivers

GPS ANTENNA WITH METALLIC CONICAL STRUC- TURE FOR ANTI-JAMMING APPLICATIONS

LOW POWER GLOBAL NAVIGATION SATELLITE SYSTEM (GNSS) SIGNAL DETECTION AND PROCESSING

Adaptive Beamforming for Multi-path Mitigation in GPS

Methodology and Case Studies of Signal-in-Space Error Calculation Top-down Meets Bottom-up

Vector tracking loops are a type

A Simulation Tool for Space-time Adaptive Processing in GPS

Satellite Navigation Principle and performance of GPS receivers

GNSS OBSERVABLES. João F. Galera Monico - UNESP Tuesday 12 Sep

The WAAS/L5 Signal for Robust Time Transfer: Adaptive Beamsteering Antennas for Satellite Time Synchronization

Characterization of Signal Deformations for GPS and WAAS Satellites

POWERGPS : A New Family of High Precision GPS Products

The Case for Narrowband Receivers

ADAPTIVE ANTENNAS. TYPES OF BEAMFORMING

Antenna Arrays for Robust GNSS in Challenging Environments Presented by Andriy Konovaltsev

Optimal Pulsing Schemes for Galileo Pseudolite Signals

Prototype Software-based Receiver for Remote Sensing using Reflected GPS Signals. Dinesh Manandhar The University of Tokyo

WELCOME TO. The Role of GNSS Antennas in Mitigating Jamming and Interference. Co Moderator: Lori Dearman, Sr. Webinar Producer

Compass-M1 Broadcast Codes and Their Application to Acquisition and Tracking

Using Frequency Diversity to Improve Measurement Speed Roger Dygert MI Technologies, 1125 Satellite Blvd., Suite 100 Suwanee, GA 30024

Assessing & Mitigation of risks on railways operational scenarios

Improved GPS Carrier Phase Tracking in Difficult Environments Using Vector Tracking Approach

How Effective Are Signal. Quality Monitoring Techniques

Code Generation Scheme and Property Analysis of Broadcast Galileo L1 and E6 Signals

GPS/GNSS Antennas. В. Rama Rao W. Kunysz R. Fante К. McDonald ARTECH HOUSE. BOSTON LONDON artechhouse.com

High Gain Advanced GPS Receiver

Proceedings of Al-Azhar Engineering 7 th International Conference Cairo, April 7-10, 2003.

Research Article Assessment of Measurement Distortions in GNSS Antenna Array Space-Time Processing

A Survey on SQM for Sat-Nav Systems

Antenna Measurements using Modulated Signals

3D-Map Aided Multipath Mitigation for Urban GNSS Positioning

MITIGATION OF SIGNAL BIASES INTRODUCED BY CONTROLLED RECEPTION PATTERN ANTENNAS IN A HIGH INTEGRITY CARRIER PHASE DIFFERENTIAL GPS SYSTEM

MAKING TRANSIENT ANTENNA MEASUREMENTS

GAJET, a DRDC Evaluation Testbed for Navigation Electronic Warfare. Michel Clénet

Adaptive beamforming using pipelined transform domain filters

Utilizing Batch Processing for GNSS Signal Tracking

The Influence of Multipath on the Positioning Error

UNIT 1 - introduction to GPS

A Slope-Based Multipath Estimation Technique for Mitigating Short-Delay Multipath in GNSS Receivers

Some Notes on Beamforming.

The Case for Recording IF Data for GNSS Signal Forensic Analysis Using a SDR

Analysis on GNSS Receiver with the Principles of Signal and Information

Signal Quality Checks For Multipath Detection in GNSS

Every GNSS receiver processes

Using a Sky Projection to Evaluate Pseudorange Multipath and to Improve the Differential Pseudorange Position

GPS Receiver Autonomous Interference Detection

3 RANGE INCREASE OF ADAPTIVE AND PHASED ARRAYS IN THE PRESENCE OF INTERFERERS

Orion-S GPS Receiver Software Validation

GPS software receiver implementations

Receiving the L2C Signal with Namuru GPS L1 Receiver

Performance and Jamming Test Results of a Digital Beamforming GPS Receiver

Unambiguous BOC Acquisition in Galileo Signal

Chapter 5. Clock Offset Due to Antenna Rotation

PERFORMANCE ASSESSMENT OF MAXIMUM LIKELIHOOD IN THE DETECTION OF CARRIER INTERFERENCE CORRUPTED GPS DATA IN MOBILE HANDSETS

University of Bristol - Explore Bristol Research. Link to published version (if available): /VTCF

GNSS for Landing Systems and Carrier Smoothing Techniques Christoph Günther, Patrick Henkel

Evaluation of L2C Observations and Limitations

Adaptive Antenna Array Processing for GPS Receivers

UHF Phased Array Ground Stations for Cubesat Applications

BIOGRAPHY ABSTRACT. This paper will present the design of the dual-frequency L1/L2 S-CRPA and the measurement results of the antenna elements.

It is well known that GNSS signals

K.NARSING RAO(08R31A0425) DEPT OF ELECTRONICS & COMMUNICATION ENGINEERING (NOVH).

Survey on Interference Mitigation via Adaptive Array Processing in GPS

Monitoring Station for GNSS and SBAS

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

GPS receivers built for various

Nominal Signal Deformations: Limits on GPS Range Accuracy

B SCITEQ. Transceiver and System Design for Digital Communications. Scott R. Bullock, P.E. Third Edition. SciTech Publishing, Inc.

Analysis of Processing Parameters of GPS Signal Acquisition Scheme

Demonstrations of Multi-Constellation Advanced RAIM for Vertical Guidance using GPS and GLONASS Signals

Decoding Galileo and Compass

Advanced Communication Systems -Wireless Communication Technology

Cancellation of Space-Based Interference in Radio Telescopes 1. Lou Nigra 2. Department of Astronomy University of Wisconsin Madison, Wisconsin

Smart antenna technology

Transcription:

Navigation Accuracy and Interference Rejection for an Adaptive GPS Antenna Array David S. De Lorenzo, Stanford University Jason Rife, Stanford University Per Enge, Stanford University Dennis M. Akos, University of Colorado BIOGRAPHY David De Lorenzo is a member of the Stanford University GPS Laboratory, where he is pursuing a Ph.D. degree in Aeronautics and Astronautics. He received a Master of Science degree in Mechanical Engineering from the University of California, Davis, in 199. David has worked previously for Lockheed Martin and for the Intel Corporation. Dr. Jason Rife is a Research Associate studying JPALS and the Local Area Augmentation System (LAAS) at Stanford University. After receiving his B.S. in Mechanical and Aerospace Engineering from Cornell University (199), he spent one year working in the turbine aerodynamics group of the commercial engine division of Pratt & Whitney. He resumed his studies at Stanford and earned M.S. (1999) and Ph.D. () degrees in Mechanical Engineering. His thesis work focused on positioning and navigation technologies for underwater robots. Dr. Per Enge is a Professor of Aeronautics and Astronautics at Stanford University, where he is the Kleiner-Perkins, Mayfield, Sequoia Capital Professor in the School of Engineering. He directs the GPS Research Laboratory, which develops satellite navigation systems based on the Global Positioning System. Dr. Enge has received the Kepler, Thurlow, and Burka Awards from the Institute of Navigation for his work. He is a Fellow of the Institute of Navigation and the Institute of Electrical and Electronics Engineers. Dr. Dennis M. Akos is an Assistant Professor with the Aerospace Engineering Science Department at the University of Colorado at Boulder. He also has served as a faculty member with the Luleå Technical University, Sweden, and as a Research Associate in the GPS Laboratory at Stanford University. Dr. Akos completed the Ph.D. degree in Electrical Engineering at Ohio University within the Avionics Engineering Center. ABSTRACT In the Joint Precision Approach and Landing System (JPALS) both the differential GPS reference station and the airborne user will employ Controlled Reception Pattern Array (CRPA) antennas. In a high precision system such as JPALS, CRPA antennas may suffer phase center biases that introduce significant integrity risk. These phase center biases result from both hardware design and from algorithm selection (beamsteering and/or nullforming). This study shows that there is a clear tradeoff between radio frequency interference (RFI) rejection and the introduction of biases in the pseudorange and carrier-phase navigation outputs from a space-time adaptive processor (STAP) GPS receiver. Deterministic corrections based either on single-element or array calibration (and implemented as a line-of-sight-based lookup table) will reduce pseudorange and carrier-phase biases in the tracking output. For the STAP algorithms and patch-element-based antenna array considered here, the carrier-phase bias residuals are on the order of -1 and the pseudorange bias residuals are in the 1 s of cm. While the carrier-phase residuals are likely tolerable for high-integrity carrier-phase-differential integer resolution, the code-phase residuals are troubling and will need further work in regards either to algorithm development, to antenna design improvements, or to both. INTRODUCTION The Joint Precision and Approach Landing System (JPALS) is a United States Navy and Air Force program to provide local-area augmentation to the on-board GPS navigation solution for pilots on approach to aircraft carrier, fixed-base, and tactical airfields. Sea-based JPALS provides carrier-phase-differential navigation, and

Correlation Output x 1 5 3 PRN 1 PRN PRN 3 PRN PRN 5 PRN PRN 7 PRN 8 PRN 9 PRN 1 PRN 11 Test satellite ideal correlation peak Phase 3 1 Antenna Phase vs. Frequency Response, w.r.t. Look Direction PRN 1 PRN PRN 3 PRN PRN 5 PRN PRN 7 PRN 8 PRN 9 PRN 1 1-1 -1.5-1 -.5.5 1 1.5 Code Offset (chips) - 15 155 1 15 Frequency (MHz) Figure 1. Antenna gain/phase response varies as a function of incoming signal azimuth, elevation, and frequency L1 microstrip patch antenna []. Left: correlation peak distortion. Right: carrier-phase bias. as such there is the requirement to maintain high-integrity GPS measurements while simultaneously rejecting radio frequency interference (RFI) and multipath. To increase the available C/No, the baseline JPALS architecture includes a multi-element antenna array. Adaptive beamsteering and nullforming are being studied to further improve interference rejection. Unlike the response of an ideal isotropic element, that of a single-element fixed reception pattern antenna (FRPA) varies as a function of incoming signal azimuth, elevation, and frequency. This distortion leads to biases in codephase (pseudorange) and carrier-phase, and possibly to attenuation of the incoming signal power (Figure 1 and Table 1). Bias errors, such as those shown in Table 1 for microstrip patch antennas developed at Stanford University [1], may be sufficient to make integer determination and carrier-phase-differential navigation problematic. In the GPS solution, the mean bias in pseudorange across all incoming signals is assigned to the user clock offset term; however, the residual pseudorange bias after this correction is still on the order of m. The carrier-phase bias needs to be a small fraction of a carrier cycle, which clearly is not the case for the microstrip patch antenna elements characterized here. A multi-element antenna array may introduce additional biases in the pseudorange and carrier-phase estimates. Not only is this due to mutual electronic coupling between elements that does not exist for stand-alone antennas, but also it is a consequence of distortion that may be introduced by the spatial and temporal weighting in forming the composite array output signal. Table 1. Antenna response and signal attenuation L1 microstrip patch antenna []. Pseudorange and carrier-phase biases for the center antenna of a 7-element array, as a function of incoming signal line-of-sight; isotropic signal power of db-hz. PRN Incoming Signal Pseudorange Carrier-Phase C/No Line-of-Sight Bias (m) Bias (deg) (db-hz) Az El 1 -. 85 37.5 3 3-1.9 3.8 3-1.5-7 37. 9 5 -.3-83 38.1 5 1 -.8-15 3.3 15 -. 15 38.7 7 1 7 -. -19 38.9 8 -.5 11 33. 9 7 3-3. 7 35.9 1 3 8 -.7 7 38.8

In order to increase the probability of correctly fixing integer ambiguities during carrier-phase-differential navigation, some form of bias mitigation is necessary. This could be implemented, for example, by frequencydomain equalization or by a line-of-sight-based bias calibration look-up table []. Equalization could take the form of appropriate filtering of the incoming satellite signals to undo the distortion caused by the antenna gain and phase response, and therefore would apply to any subsequent method of calculating antenna weights. However, since the antenna response is a strong function of the incoming signal line-of-sight, the filter coefficients would likewise be line-of-sight dependent, requiring a massive database of filter coefficient terms. An alternative approach would be to apply deterministic pseudorange and carrier-phase corrections to the tracking output quantities themselves the corrections would be stored in a look-up table and tagged according to the incoming signal line-of-sight. The look-up table approach applies corrections for only one set of antenna weights. Because this mitigation method is simpler to implement and verify, it is considered preferable to frequencydomain equalization and is the method studied in the remainder of this report. In non-adaptive spatial-only beamsteering the beamsteering weights are calculated deterministically based on satellite position and array orientation. Consequently, biases may be eliminated with a precise look-up table, namely pseudorange and carrier-phase corrections tagged to incoming signal line-of-sight (Figure ). The use of multiple antennas inherently improves C/No relative to the single antenna case. However, the drawbacks of a deterministic controlled reception pattern antenna array (CRPA) include high sidelobes and uncontrolled nulls that may not effectively reject incoming interference or signal multipath. Adaptive processing for noise-rejection or powerminimization allows the automatic suppression of narrowband interference [3-] but may exacerbate bias errors compared to deterministic weighting. The reason for the modified biases is that adaptive weighting of single-element distortion from the slave antennas shifts the pseudorange and carrier-phase measurements. Even greater interference rejection, particularly of wideband sources, may be realized by incorporating temporal filtering in the array processing e.g., with a tapped delay line antenna array [5-9]. However, adding time taps to allow space-time adaptive processing (STAP) yields a finite impulse response filter that may further distort the spread-spectrum GPS ranging signal [1-11]. This paper evaluates the trade-offs between pseudorange and carrier-phase bias errors and RFI rejection in adaptive, multi-antenna GPS receivers. To this end, array performance was evaluated with regards to both code- and carrier-phase estimation errors as well as tracking output signal-to-noise ratio. Adaptive algorithms that were considered in this study include: 1) A least-mean-square (LMS) error approach to weight vector determination that seeks to Deterministic Biases Deterministic Corrections Front-end Code & carrier wipeoff Correlators Tracking Loops Bias Compensation Code & carrier wipeoff Receiver Signal Processing Bias Calibration Look-up Table (code & carrier az/el-dependent corrections) PVT Processing Corrected estimates of code & carrier phase Figure. Pseudorange and carrier-phase bias compensation.

carrier wipeoff Weight Control Algorithm code wipeoff Correlators integrate & dump Tracking Loops sin Early FFT-based Coherent Acquisition cos Carrier NCO reference signal Prompt Code NCO Late Weight control processing supports FRPA & CRPA antennas and adaptive LMS & Applebaum arrays C/No Estimator Pseudorange & Carrier-phase Bias Estimator Figure 3. Software receiver block diagram. minimize the difference between the actual array output and the value of a pilot or reference signal this approach can also be termed blindadaptive beamsteering/nullforming, since no knowledge of incoming signal line-of-sight or array geometry is required [1-13]. ) A steering-vector-based (Applebaum) approach that constrains the array response in the direction of the arriving signal (neglecting for the moment the effects of antenna distortion in the calculation of the weight vector) and then seeks to minimize the remaining array output power [1-15]. These algorithms were developed as part of a softwarebased, all-in-view GNSS receiver that is capable of tracking either simulated or actual satellite signals (Figure 3). In conjunction with this software receiver, a signal simulation environment was created that allows generation of GPS C/A-code or P-code signals incident on a multi-element array, as well as injection of narrow-band or wide-band interference. An important feature of this simulator is the ability to incorporate electronic models of the gain and phase response of the various antennas in the array as a function of incident signal arrival direction (Figure ). Using this signal simulation environment and software receiver, it is possible to show the advantages of adaptive arrays with respect to RFI rejection, and to quantify their pseudorange and carrier-phase bias residuals after suitable look-up table corrections are applied. Deterministic corrections based either on single-element or array calibration (and implemented as a line-of-sight-based lookup table) will reduce pseudorange and carrier-phase biases in the tracking output. For the STAP algorithms and patch-element-based antenna array considered here, the carrier-phase bias residuals are on the order of -1 and the pseudorange bias residuals are in the 1 s of cm. While carrier-phase residuals are likely tolerable for highintegrity carrier-phase-differential integer resolution, the code-phase residuals are troubling and will need further work in regards either to algorithm development, to antenna design improvements, or to both. METHODOLOGY Antennas introduce biases in the estimates of code-phase and carrier-phase due to their variation in response as a function of incoming signal azimuth, elevation, and frequency. Multi-element antenna arrays, while increasing C/No and nulling/rejecting RFI, may compound these distortion-induced biases. If the weights themselves are deterministic functions of satellite and array geometry, then the biases are likewise deterministic. However, if the weights are calculated adaptively, then the biases are determined by the signal environment, STAP algorithm, and receiver tracking implementation. Therefore, in order to evaluate pseudorange and carrierphase corrections, or the bias residuals after deterministic compensation, it is necessary to evaluate receiver

C/A & P code 1.3 & 1.3 Mchip/sec Gain/phase data for each look direction (1 S/V constellation) and antenna (7 rectangular patches) Bandlimited WGN & CW interference Front-end A/D A/D Front-end 1. 1.8 Gain/phase response for PRN #1 [Azimuth =, Elevation = ] Antenna Gain Response for PRN #1 [Az = deg, El = deg] Antenna Phase Response for PRN #1 [Az = deg, El = deg] 5 Antenna #1 Antenna #1 Antenna # Antenna # Antenna #3 Antenna #3 Antenna # Antenna # Antenna #5 Front-end BW ~MHz Antenna #5 Antenna # 15 Antenna # Antenna #7 Antenna #7 Gain. Phase 1. 5 Front-end BW ~MHz. 18 15 15 15 15 158 1 1 1 1 18 Frequency (MHz) -5 15 15 155 1 15 17 Frequency (MHz) Figure. Software-based signal simulator models satellite signal reception, correlated interference, non-correlated, Gaussian noise, and antenna-induced signal distortion. performance in the context of realistic GPS tracking scenarios. This dictates an end-to-end study that encompasses antenna characterization, signal/rfi generation, adaptive weight computation, and signal tracking. Then, the receiver estimates of the tracking observables (code-phase, carrier-phase, and C/No) may be compared to the truth values from the signal simulator. With all variables under direct control, this allows isolation and estimation of the parameters of interest. A software signal generator, as shown in Figure, was developed to produce P-code-like signals mixed-down from the GPS-L1 frequency to a suitable intermediate frequency for sampling (f S = 8MHz, f IF = MHz, frontend bandwidth of MHz). The spreading-code modulation used the P-code generator and a chipping rate of 1.3 Mchips/sec, but short-cycled after the first millisecond to make implementation of the tracking code more straightforward. Time-offsets based on incoming signal line-of-sight and array geometry were applied to advance or retard the received signal with respect to the array physical center; the array comprised a center master element and a hexagonal pattern of slave elements with λ/ baselines. A standard constellation of 1 satellites was used (Figure 5), with locations selected to correspond to previously-identified extremes in the code/carrier bias maps. The ratio of signal to white Gaussian noise (WGN) was set to deliver a C/No of db-hz for an isotropic receiving element; this noise is not correlated across antenna elements. In addition, there is wideband interference incident on the array this interference is correlated in space/time between the various elements of the array. The RFI emitters are placed at or near the horizon, and the power level can be varied to achieve the desired jammer-to-signal (J/S) power ratio. The signal generator code also simulated signal distortion introduced by the multi-element antenna array. The GPS signals were received either by an array of isotropic antenna elements or by a 7-element array whose gain and phase response characteristics were as determined by previous electronic simulation models and anechoic chamber measurements []. The method of incorporating the gain and phase characteristics is by frequency-domain convolution with the generated satellite signals. Finally, real-valued samples are stored to disk at -bit resolution, after passing through an automatic gain-control stage to ensure good dynamic range performance. The generated signals were then fed as inputs into a software-based GNSS receiver (Figure 3). This receiver uses standard PLL and DLL tracking loops, and has the

PRN Azimuth Elevation 1 3 3 3 9 5 5 1 RFI Azimuth Elevation 15 1 7 1 7 5 8 3 1 9 7 3 5 1 3 8 5 5 Figure 5. Satellite and RFI constellation. 7 1 capability of tracking either GPS C/A code, pseudo P- code originating from the previously described signal simulator, or Galileo signals recorded from the GIOVE-A satellite either on the L1, E5, or E bands. The receiver is a flexible, all-in-view, MATLAB implementation, which incorporates some organizational and datastructure elements from other publicly-available code bases [1]. Distinguishing this receiver implementation is the inclusion of array-processing modules, either for deterministic beamsteering or for adaptive beamsteering and nullforming using space/time adaptive processing. Since the tracking outputs can be compared directly with the known code-phase and carrier-phase input values from the signal simulator, it is possible to determine precisely the corresponding pseudorange and carrierphase biases due to antenna-induced distortion effects and array processing. Likewise, the software receiver can estimate the C/No ratio for the received signal, which allows quantification of the noise/rfi performance of the various array processing approaches (Figure ). In determining pseudorange and carrier-phase biases, it is important to suppress the effects of noise in the output signal to obtain a precise estimate. Because the pseudorange standard deviation shown in Figure is in the 1 s of centimeters, substantial filtering would be needed to attenuate the noise effects. Instead, a two-step process is implemented to first track the signals in the presence of WGN (noise uncorrelated across elements) and RFI (correlated interference) while updating the antenna weights based on the selected adaptation scheme (Figure 7a), and then second to track the signals without WGN or RFI, while playing back the stored weight vector from the previous step at each processing epoch (Figure 7b). In this way, bias values may be computed directly from a short-duration simulation (Figure 8). RESULTS The first important result from the signal simulation and tracking process is the characterization of the RFI-free pseudorange and carrier-phase biases and C/No performance of the non-ideal single-element FRPA, the 7- element deterministic CRPA, and the Averages over 1-satellite constellation Pseudorange bias (m) Carrier-phase bias (deg) C/No (db-hz) Table. Uncompensated code-phase and carrier-phase biases. Uncompensated pseudorange and carrier-phase biases; isotropic signal power of db-hz. Single-element FRPA 7-element deterministic CRPA Blind-adaptive STAP beam/null steering (LMS-based) Steering-vector STAP beam/null steering (Applebaum-based).13 1.88 1.9 1.88 9 9 9 9 37. 5.7 5.7 5.7

Code Phase Bias (m) 1 8 - - Code Phase Biases PRN 1 PRN PRN 3 PRN PRN 5 PRN PRN 7 PRN 8 PRN 9 PRN 1 PRN 11 Carrier Phase Bias (deg) 15 1 5-5 Carrier Phase Biases C/No (db-hz) 55 5 5 35 3 Filtered C/No -1 5 - -8-15 -1 1 3 5-1 3 5 15 1 3 5 Figure. Software receiver tracking output and noise/bias estimation; Applebaum-based STAP, C/No = db-hz plus six RFI sources at J/S = 3dB. beamsteering/nullforming adaptive arrays (Table ). The bias levels are clearly unacceptable for carrier-phasedifferential integer resolution and navigation. For the antenna array studied, there is an average code-phase bias due to correlation peak distortion of approximately m (standard deviation ~7cm). The corresponding carrierphase biases are essentially uniformly distributed from - 18 to +18 degrees. As expected, the multi-element antennas achieved better C/No than the single-element FRPA. Carrier-to-noise ratio was evaluated taking non-ideal antenna effects into account (Figure ). For this RFI-free case, the thermal noise level was set so that an ideal isotropic antenna would achieve a C/No of db-hz. The attenuation of signal energy due to non-ideal gain response of the single- Signal Simulator Software Receiver GPS signals + Antenna distortion Weight Control Algorithm Receiver tracking RFI sources + WGN Store weights Estimate C/No Constellation & array geometry Standard PLL & DLL P-code sampled at 8 MHz FRPA, CRPA, & STAP modules Antenna response from HFSS Spatial & temporal constraints Figure 7a. Methodology for estimating code/carrier biases and C/No.

Signal Simulator Software Receiver GPS signals + Antenna distortion Weight Vector Multiply Receiver tracking RFI sources + WGN Stored weights Estimate code & carrier biases Constellation & array geometry Standard PLL & DLL P-code sampled at 8 MHz FRPA, CRPA, & STAP modules Antenna response from HFSS Spatial & temporal constraints Figure 7b. Methodology for estimating code/carrier biases and C/No. element FRPA results in tracked C/No of about 3 db-hz below the db-hz level that would be produced for an isotropic antenna. All of the multi-element arrays yield ~8.5 db-hz improvement in C/No [1*log1(7) = 8.5 db] with respect to the single-element FRPA. The response of the adaptive arrays closely matches the deterministic case, since in the absence of RFI the converged (steady-state) adaptive weight vector solutions approach that of the deterministic CRPA. RFI-free simulations also revealed that different types of look-up tables should be used for different CRPA Step 1: Track signals containing RFI & noise Save weight vector for STAP algorithms Determine noise performance: C/No Step : Track noise-free signals Use stored weight vector for STAP algorithms Determine code & carrier biases Code Phase Bias (m) Code Phase Biases Carrier Phase Biases 1 PRN 1 PRN PRN 3 8 PRN PRN 5 PRN PRN 7 PRN 8 PRN 9 PRN 1 PRN 11 - - - -8 Carrier Phase Bias (deg) 15 1 5-5 -1-15 C/No (db-hz) 55 5 5 35 3 5 Filtered C/No Code Phase Bias (m) Code Phase Biases Carrier Phase Biases Filtered C/No 1 PRN 1 PRN PRN 3 8 PRN PRN 5 PRN PRN 7 PRN 8 PRN 9 PRN 1 PRN 11 - - - -8 Carrier Phase Bias (deg) 15 1 5-5 -1-15 C/No (db-hz) 55 5 5 35 3 5 C/No estimate from step 1-1 1 3 5-1 3 5 15 1 3 5-1 1 3 5-1 3 5 15 1 3 5 Figure 8. Two-step process for estimating C/No and code/carrier biases in presence of noise/rfi. This example shows 7-antenna, Applebaum-based space-time adaptive processing (with a single time-tap), C/No of db-hz, plus six RFI sources at J/S of 3dB.

Bias residuals w.r.t. single-element FRPA calibration Bias residuals w.r.t. 7-element CRPA calibration Averages over 1-satellite constellation Table 3. Bias residuals after compensation the best calibration option depends on STAP algorithm. Pseudorange and carrier-phase biases; isotropic signal power of db-hz. Pseudorange bias (m) Carrier-phase bias (deg) Pseudorange bias (m) Carrier-phase bias (deg) Blind-adaptive STAP beam/null steering (LMS-based) Steering-vector STAP beam/null steering (Applebaum-based).3.8.5 9... 9.3. algorithms. As discussed previously, a look-up table of pseudorange and carrier-phase biases could be calibrated for either the single-element FRPA or the deterministic CRPA, reducing the antenna-induced navigation biases to zero. Either of these two look-up table compensation strategies could be applied for use with adaptive weighting algorithms. Because these look-up strategies are tuned for specific weight vector combinations, however, some calibration error is expected when a static look-up table is used with an adaptive algorithm. Of the FRPA and CRPA deterministic look-up corrections, the better compensation option depends on the adaptive algorithm (the blue-shaded cells in Table 3). For LMS-based blind/adaptive beamsteering/nullforming, pseudorange bias corrections should come from the deterministic CRPA calibration, while the carrier-phase corrections should come from the single-element FRPA calibration. The reason for the near-zero carrier-phase bias with respect to the single-element FRPA is that the weight adjustment algorithm seeks to steer all desirable incoming signal energy to the in-phase channel of the master element. For Applebaum-based steering-vector beamsteering/nullforming, both pseudorange and carrierphase corrections should come from the deterministic CRPA calibration. The reason for the effectiveness of this compensation scheme in the RFI-free case is that the converged weight vector is identical to that of the deterministic CRPA. By comparing high-rfi simulations (J/S = 5dB for each of six wideband RFI emitters) to the RFI-free case, it is possible to define the trade-offs between bias and noise performance for various adaptive weighting schemes (Figure 9). In Figure 9, the horizontal axis shows the improvement in C/No for the adaptive algorithms with respect to the deterministic CRPA for the RFI-free and high-rfi conditions, while the vertical axis shows the increase in pseudorange and carrier-phase bias after bias compensation. It should be noted that more noise rejection and lower biases are to be desired this corresponds to movement in Figure 9 down and to the right. While in the RFI-free case the C/No performance of the adaptive arrays matches that of the deterministic CRPA, the presence of RFI causes the adaptive arrays to clearly outperform the deterministic CRPA. This is demonstrated by the C/No improvement for the STAP algorithms with respect to the deterministic CRPA when including six 5dB RFI emitters (i.e., J/S for each emitter is 5 db). The increase in C/No with respect to the deterministic CRPA comes at the expense of greater pseudorange and carrier-phase biases; this is shown in Figure 9 as the increase in pseudorange and carrier-phase bias residual after compensation for the six RFI emitter case. Neither the LMS-based nor Applebaum-based algorithms clearly outperform to other in all respects. Figure 9 clearly shows that steering-vector-based STAP (the Applebaum approach) yields better carrier-phase performance than blind-adaptive beamsteering and nullforming (the LMS approach) in the presence of RFI when calibration data are available. However, blindadaptive STAP has better RFI rejection capability than the steering-vector-based method for the scenario herein described, i.e. 1dB vs. 8dB C/No improvement with respect to the deterministic CRPA, or a db advantage. Both STAP approaches produce comparable code-phase bias residuals. CONCLUSION This paper has evaluated the trade-offs between pseudorange and carrier-phase bias errors and RFI rejection in adaptive, multi-antenna GPS receivers. Deterministic corrections based either on single-element

Bias Residuals - 7-element Array with Single-element and Array Calibration Bias Residuals - 7-element Array with Single-element and Array Calibration Code-phase Bias (m).5..3..1 no RFI LMS, single time-tap LMS, 3 time-taps Applebaum, single time-tap Applebaum, 3 time taps six 5dB RFI sources 8 1 1 C/No Improvement vs. Deterministic CRPA (db-hz) Carrier-phase Bias (deg) 1 8 no RFI LMS, single time-tap LMS, 3 time-taps Applebaum, single time-tap Applebaum, 3 time taps Figure 9. Bias residuals and C/No improvement for STAP arrays with respect to a 7-element deterministic beamsteering array. six 5dB RFI sources 8 1 1 C/No Improvement vs. Deterministic CRPA (db-hz) or array calibration (and implemented as a line-of-sightbased lookup table) reduce pseudorange and carrier-phase biases in the tracking output. For the STAP algorithms and patch-element-based antenna array considered here, the carrier-phase bias residuals are on the order of -1 and the pseudorange bias residuals are in the 1 s of cm. While carrier-phase residuals are likely tolerable for highintegrity carrier-phase-differential integer resolution, the code-phase residuals are troubling and will need further work in regards either to algorithm development, to antenna design improvements, or to both. There is significant future relevance to the work described here, not only for military landing systems but also for civilian receivers operating with next-generation GNSS signals. This is due to the facts that pseudorange and carrier-phase biases are a strong function of the spreading-code bandwidth and that military and nextgeneration signals have higher code chipping rates than the current GPS C/A code rate of 1.3 Mchips/sec. In addition, the increased attention paid to denial-of-service due to unintentional or harmful RFI makes STAP algorithms desirable for safety-of-life applications. ACKNOWLEDGMENTS The authors gratefully acknowledge the support of the JPALS Program Office, and the Naval Air Warfare Center Aircraft Division through contract N1-5-C- 8. REFERENCES 1. Kim, U.S.,, Mitigation of Signal Biases Introduced by Controlled Reception Pattern Antennas in a High Integrity Carrier Phase Differential GPS System, PhD dissertation defense, Stanford University.. Kim, U.S., 5, Analysis of Carrier Phase and Group Delay Biases Introduced by CRPA Hardware, Proc. ION GNSS 5. 3. Gecan, A. and Zoltowski, M., 1995, Power Minimization Techniques for GPS Null Steering Antenna, Proc. ION GPS 1995, pp. 81-88.. Nicholson, B.W., Upton, D.M., Cotterill, S., Marchese, J., Upadhyay, T., and Vander Velde, W.E., 1998, Computer Simulations of Digital Beam Forming Adaptive Antennae for GPS Interference Mitigation, Proc. ION NTM 1998, pp. 355-3. 5. Hatke, G.F., 1998, Adaptive Array Processing for Wideband Nulling in GPS Systems, Proc. Thirty- Second Asilomar Conf. on Signals, Systems & Computers, vol., pp. 133-133.. Fante, R.L. and Vaccaro, J.J.,, Wideband Cancellation of Interference in a GPS Receive Array, IEEE Trans. on Aerospace and Electronic Systems, vol. 3, no., pp. 59-5. 7. Gupta, I.J. and Moore, T.D., 1, Space-Frequency Adaptive Processing (SFAP) for Interference Suppression in GPS Receivers, Proc. ION NTM 1, pp. 377-385. 8. Moore, T.D.,, Analytic Study of Space-Time and Space-Frequency Adaptive Processing for Radio Frequency Interference Suppression, PhD. Dissertation, The Ohio State University.

9. Hatke, G.F. and Phuong, T.T.,, Design and Test of a GPS Adaptive Antenna Array Processor: The Multipath Adaptive Multibeam Array (MAMBA) Processor, Lincoln Laboratory Project Report GPS-1, Massachusetts Institute of Technology. 1. Fante, R.L., Fitzgibbons, M.P., and McDonald, K.F.,, Effect of Adaptive Array Processing on GPS Signal Crosscorrelation, Proc. ION GNSS, pp. 579-583. 11. McGraw, G.A., Ryan Young, S.Y., and Reichenauer, K.,, Evaluation of GPS Anti-Jam System Effects on Pseudorange and Carrier Phase Measurements for Precision Approach and Landing, Proc. ION GNSS, pp. 7-751. 1. Widrow, B., Mantey, P.E., Griffiths, L.J., and Goode, B.B., 197, Adaptive Antenna Systems, Proc. IEEE, vol. 55, no. 1, pp. 13-159. 13. Widrow, B. and Stearns, S.D., 1985, Adaptive Signal Processing, Pearson Education, Inc. 1. Applebaum, S.P., 197, Adaptive Arrays, IEEE Trans. on Antennas and Propagation, vol., no. 5, pp 585-598. 15. Compton, R.T., 1988, Adaptive Antennas, Prentice- Hall, Inc. 1. Borre, K., Akos, D.M., Bertelsen, N., Rinder, P., and Jensen, S.H., A Software-Defined GPS and Galileo Receiver: A Single-Frequency Approach, Birkhauser, in press.