GPS HDL Acquisition and Tracking Using C/A Code

Model Overview

The model that you use in this example mainly consists of acquisition and tracking modules. The input to the model is a GPS baseband waveform with a sample rate of 32.768 Msps. This baseband waveform contains multiple satellite waveforms with Doppler offsets, code phase offsets, and additive white Gaussian noise (AWGN) noise in it. The model upsamples the GPS waveform by a factor of six. This increases the clock from 32.768 MHz to 196.6 MHz to achieve faster acquisition. The model supports sample rates 8.192 Msps, 16.384 Msps, and 32.768 Msps. The example is explained with the sample rate of 32.768 Msps.

The Acquisition block decimates the GPS waveform from 32.768 Msps to 4.096 Msps and operates at this sampling rate. The acquisition block detects the satellites and generates coarse estimates for them. You can search specified pseudorandom noise identification numbers (PRNIDs) in the received waveform, instead of searching all the 32 PRNIDs.

The Tracking block uses the GPS waveform at 32.768 Msps, the PRNID of the detected satellites, and the coarse estimates to track the satellites. At least four satellites are required for GPS position estimation, so tracking does not start until the acquisition process detects four satellites. The Tracking block contains eight instances of tracking logic to simultaneously track eight satellites. The Tracking block uses phase-locked loop, frequency-locked loop, and delay-locked loop to recover the pi/2-BPSK modulated LNAV symbols of each detected satellite. This figure shows an overview of the example model.

File Structure

This example uses one Simulink model, two functions, and two scripts.

Model Interface

This figure shows the structure of the gpshdlAcquisitionTracking model.

Model Inputs

Model Outputs:

The output ports are vectors of length 8 because the example tracks eight satellites simultaneously.

Input Configuration

Double-click the Input Configuration subsystem to configure the transmitter parameters and channel impairments.

Model Structure

This figure shows the Acquisition and Tracking subsystem. This subsystem consists of the Acquisition, Timing Adjust, and Tracking subsystems.

Acquisition Subsystem

The Acquisition subsystem accepts the GPS waveform at a sampling rate of 32.768 Msps, decimates it to 4.096 Msps, and stores one millisecond duration of the decimated waveform. The subsystem selects coarse Doppler frequencies sequentially from -10 kHz to 10 kHz in steps of 1 kHz, generates local carrier waveforms at these frequencies, compensates for these frequencies in the decimated waveform, and outputs carrier-wiped-off waveform. The subsystem converts the carrier-wiped-off waveform into the frequency domain using a 4096-point fast Fourier transform (FFT). The subsystem fetches the frequency-domain C/A code from a look-up table (LUT) and correlates it with the waveform to find the correlation peak. When the peak is greater than a dynamic threshold, the subsystem detects the satellite with the C/A code. The corresponding Doppler frequency and C/A code phase are the coarse estimates of the satellite. The subsystem performs this frequency-domain correlation for four satellites in parallel and for eight sequential searches to finish searching all 32 GPS satellites. The subsystem then sorts and selects the eight detected satellites with the strongest correlation peaks. The subsystem also generates a 1 ms epoch signal, which asserts, for every 1 ms, to use that signal for timing adjust.

The Acquisition subsystem contains these main subsystems:

Timing Adjust Subsystem

The Timing Adjust subsystem accepts the GPS waveform, the detected satellite PRNIDs, coarse Doppler frequencies, coarse code phase offsets, and the 1 ms epoch signal from the Acquisition subsystem. After the Acquisition subsystem finishes detecting the satellites, the Timing Adjust subsystem outputs valid signal for the Tracking subsystem starting from the instant the 1 ms epoch signal asserts.

Tracking Subsystem

The Tracking subsystem tracks the satellites detected through acquisition. The subsystem accepts the time-adjusted waveform at 32.768 Msps, the detected satellite PRNIDs, coarse Doppler offsets, and code phase offsets. The subsystem uses the coarse Doppler estimate and phase estimate to generate a local carrier using an NCO, removes the Doppler offset from the waveform, and returns a carrier-wiped-off waveform. The subsystem uses the detected PRNID and the coarse code phase offset to fetch replica C/A code from an LUT. The subsystem generates early, prompt, and late versions of this C/A code to correlate with the carrier-wiped-off waveform. The subsystem integrates these correlated outputs for every 1 millisecond (predetection integration time) and returns integrated samples after every millisecond. The integrated prompt outputs are the LNAV symbols of the receiver. The subsystem uses the integrated early and late outputs to estimate the delay error and the integrated prompt output to estimate the frequency and phase errors. The subsystem filters these errors using loop filters and feeds the filtered values to the NCO and C/A code LUT to aid local carrier generation and replica C/A code generation. This tracking logic applies to a single satellite. Similarly, the Tracking subsystem generates eight instances of this tracking logic and tracks eight satellites simultaneously.

The Tracking subsystem contains the Tracking Core subsystem that carries out the tracking logic. The Tracking Core subsystem contains these main subsystems:

Run Model

Open the gpshdlAcquisitionTracking model and double-click the Input Configuration subsystem to change the transmitter configuration and channel impairments. Run the model.

Note: It might take around 25 minutes to complete the simulation.

    Input PRNID    Detected PRNID
    ___________    ______________

        11               11      
        18               18      
        22               22      
        23               23      

    Input code phase offset    Estimated code phase offset
    _______________________    ___________________________

            300.34                       300.35           
            312.88                       312.86           
            587.21                       587.21           
            425.89                       425.86           

    Input doppler offset    Estimated doppler offset
    ____________________    ________________________

            3289                     3265.4         
            1568                     1606.6         
            5856                     5871.7         
            7796                     7838.9         

Generate HDL Code

To generate the HDL code, you must have an HDL Coder™ license. Use makehdl and makehdltb functions to generate HDL code and a HDL test bench for the Acquisition and Tracking subsystem.

The resulting HDL code is synthesized for a Xilinx® Zynq®-7000 ZC706 evaluation board. This table shows the post place and route resource utilization for sample rate of 32.768 Msps. The maximum frequency of operation is 202.76 MHz.

       Resources       Usage
    _______________    _____

    Slice LUT          68297
    Slice Registers    73005
    RAMB36             82   
    RAMB18             50   
    DSP48              233  

Appendix

This example uses these helper files:

References

[1]. Kaplan, Elliot D., and C. Hegarty, eds., _Understanding GPS/GNSS: Principles and Applications. Third edition. GNSS Technology and Applications Series. Boston; London; Artech House, 2017.

[2]. IS-GPS-200, Rev: L. "NAVSTAR GPS Space Segment/Navigation User Segment Interfaces." GPS Enterprise Space & Missile Systems Center (SMC) - LAAFB. https://www.gps.gov/technical/icwg/IS-GPS-200L.pdf

[3]. Ward, P.W. "GPS Receiver Search Techniques." In Proceedings of Position, Location and Navigation Symposium - PLANS '96, 604 - 11. Atlanta, GA, USA: IEEE, 1996. https://doi.org/10.1109/PLANS.1996.509134.