Introduction

Generate a GPS signal by using these three steps.

This figure shows the block diagram of the GPS waveform generator from the configuration parameters.

GPS Signal Structure

The standard [1] describes the transmission of GPS signals on two frequencies: L1 (1575.42 MHz) and L2 (1227.60 MHz). A signal of base frequency 10.23 MHz generates both of these signals. The signal frequency of L1 is 154 x 10.23 MHz = 1575.42 MHz, and the signal frequency of L2 is 120 x 10.23 MHz = 1227.60 MHz. The code-division-multiple-access technique enables you to differentiate between the satellites even though all GPS satellites transmit on the same frequency.

With the evolution of GPS, the signal structure is expanded to improve the navigation performance. For instance, LNAV data is transmitted on the L1-band from the inception of GPS. From GPS block IIR-M onwards, CNAV data is transmitted on the L2-band in addition to the LNAV data that exists on both L1 and L2 bands. This table shows these signal configurations and GPS evolution. "Dc" represents the CNAV data bits and "D" represents the LNAV data bits. Bit operator, $\oplus$, represents the XOR operation.

To select the content for transmitting over the in-phase and quadrature-phase branches, use this table.

IBranchContent = "P(Y) + D"

IBranchContent = 
"P(Y) + D"

QBranchContent = "C/A + D"

QBranchContent = 
"C/A + D"

Many of the transmitted properties in the LNAV and CNAV data are same, but the frame structure is different. The LNAV data is transmitted in 1500 bit-length frames, with each frame consisting of five subframes of 300 bits in each subframe. Because the data rate is 50 bps, transmitting each subframe takes 6 seconds and transmitting each frame takes 30 seconds. Each subframe consists of 10 words with 30 bits (24 data bits and 6 parity bits) in each word. The GPS data contains information regarding the clock and the position of the satellites. This figure shows the frame structure of the LNAV data.

The CNAV data is transmitted continuously in the form of message types. Each message type consists of 300 bits that are transmitted at 25 bps. These bits are passed through a rate-half convolutional-encoder to obtain 600 bits from each message type at 50 bps. Transmitting each message type takes 12 seconds. The standard [1] defines 14 message types in this order: 10, 11, 30, 31, 32, 33, 34, 35, 36, 37, 12, 13, 14, and 15. For a detailed description of each message type and the data transmitted, see IS-GPS-200L Appendix III [1]. The order in which each message type is transmitted is completely arbitrary but is sequenced to provide optimal user experience. In this example, you can choose the order in which these message types are transmitted. This figure shows the CNAV message structure.

The message type shown in this figure has these fields:

Set the ShowVisualizations property to enable the spectrum and correlation plot visualizations. Set the WriteWaveformToFile property to write the complex baseband waveform to a file if needed. This example enables visualizations and disables writing the waveform to a file.

ShowVisualizations =  true;
WriteWaveformToFile = false;

Specify the satellite PRN index as an integer in the range [1,63].

PRNID = 1;

Because generating the GPS waveform for the entire navigation data can take a lot of time and memory, this example demonstrates generating a waveform for only one bit of the navigation data. You can control generating the waveform for a specified number of data bits by using the property NumNavDataBits.

% Set this value to 1 to generate the waveform from the first bit of the
% navigation data
NavDataBitStartIndex = 1321;

% Set this value to control the number of navigation data bits in the
% generated waveform
NumNavDataBits = 1;

GPS Data Initialization

Initialize the data configuration object to generate the CNAV data. You can create a configuration object from the HelperGPSNavigationConfig object. Update the configuration object properties to customize the waveform as required.

cnavConfig = HelperGPSNavigationConfig(SignalType = "CNAV", PRNID = PRNID)

cnavConfig = 
  HelperGPSNavigationConfig with properties:

                           SignalType: "CNAV"
                                PRNID: 1
                         MessageTypes: [4x15 double]
                               HOWTOW: 1
                           L2CPhasing: 0
                         SignalHealth: [3x1 double]
                           WeekNumber: 2149
               GroupDelayDifferential: 0
                 ReferenceTimeOfClock: 0
                  SemiMajorAxisLength: 26560000
            ChangeRateInSemiMajorAxis: 0
                 MeanMotionDifference: 0
           RateOfMeanMotionDifference: 0
                         Eccentricity: 0.0200
                          MeanAnomaly: 0
             ReferenceTimeOfEphemeris: 0
              HarmonicCorrectionTerms: [6x1 double]
                  IntegrityStatusFlag: 0
                    ArgumentOfPerigee: -0.5200
                 RateOfRightAscension: 0
             LongitudeOfAscendingNode: -0.8400
                          Inclination: 0.3000
                      InclinationRate: 0
                              URAEDID: 0
                InterSignalCorrection: [4x1 double]
          ReferenceTimeCEIPropagation: 0
    ReferenceWeekNumberCEIPropagation: 101
                             URANEDID: [3x1 double]
                            AlertFlag: 0
                      AgeOfDataOffset: 0
                      AlmanacFileName: "gpsAlmanac.txt"
                           Ionosphere: [1x1 struct]
                     EarthOrientation: [1x1 struct]
                                  UTC: [1x1 struct]
               DifferentialCorrection: [1x1 struct]
                           TimeOffset: [1x1 struct]
                       ReducedAlmanac: [1x1 struct]
                  TextInMessageType36: 'This content is part of Satellite Communications Toolbox. '
                  TextInMessageType15: 'This content is part of Satellite Communications Toolbox. '

In theory, clocks of all satellites must be synchronized, which implies that all GPS satellite clocks must show the same time at a given moment in time. In practice, deterministic satellite clock error characteristics of bias, drift, and aging as well as satellite implementation characteristics of group delay bias and mean differential group delay exist. These errors deviate the satellite clocks from the GPS system time.

GPS satellites revolve around the Earth in an elliptical orbit, with Earth at one of the focal points of the ellipse. A set of orbital parameters that accurately defines a satellite position in this elliptical orbit is called ephemeris. Each GPS satellite transmits its own ephemeris data on subframes 2 and 3 for LNAV data and on message types 10 and 11 for CNAV data. This figure shows five orbital parameters for a satellite in the Earth-centered Earth-fixed (ECEF) coordinate system. This figure does not show an actual GPS satellite, and the figure is for illustration only.

Per Kepler's second law of planetary motion, the angular velocity (rate of change of true anomaly) is different at different locations in the orbit. You can define the mean anomaly, whose rate of change is constant over the entire orbit of the satellite. In GPS ephemeris parameters, rather than specifying true anomaly, mean anomaly is specified (from which true anomaly can be found). IS-GPS-200L Table 20-IV [1] specifies the algorithms that relate mean anomaly and true anomaly.

The eccentricity of the ellipse also defines the orbit. Eccentricity gives a measure of deviation of the elliptical orbit from the circular shape.

A similar set of properties exist for the LNAV data. Create a configuration object to store the LNAV data.

lnavConfig = HelperGPSNavigationConfig(SignalType = "LNAV", PRNID = PRNID)

lnavConfig = 
  HelperGPSNavigationConfig with properties:

                       SignalType: "LNAV"
                            PRNID: 1
                     FrameIndices: [1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25]
                       TLMMessage: 0
                           HOWTOW: 1
                    AntiSpoofFlag: 0
                        CodesOnL2: "P-code"
                      L2PDataFlag: 0
                         SVHealth: 0
                 IssueOfDataClock: 0
                            URAID: 0
                       WeekNumber: 2149
           GroupDelayDifferential: 0
    SVClockCorrectionCoefficients: [3x1 double]
             ReferenceTimeOfClock: 0
              SemiMajorAxisLength: 26560000
             MeanMotionDifference: 0
                  FitIntervalFlag: 0
                     Eccentricity: 0.0200
                      MeanAnomaly: 0
         ReferenceTimeOfEphemeris: 0
          HarmonicCorrectionTerms: [6x1 double]
             IssueOfDataEphemeris: 0
              IntegrityStatusFlag: 0
                ArgumentOfPerigee: -0.5200
             RateOfRightAscension: 0
         LongitudeOfAscendingNode: -0.8400
                      Inclination: 0.3000
                  InclinationRate: 0
                        AlertFlag: 0
                  AgeOfDataOffset: 0
        NMCTAvailabilityIndicator: 0
                          NMCTERD: [30x1 double]
                  AlmanacFileName: "gpsAlmanac.txt"
                       Ionosphere: [1x1 struct]
                              UTC: [1x1 struct]
                      TextMessage: 'This content is part of Satellite Communications Toolbox. Thank you. '

GPS Signal Generation

To generate a GPS signal at the baseband, follow these steps.

Based on the configuration, generate the CNAV data.

cnavData = HelperGPSNAVDataEncode(cnavConfig);

Pass the CNAV data through the convolutional encoder.

% Initialize the trellis for convolutional encoder
trellis = poly2trellis(7,["1+x+x^2+x^3+x^6" "1+x^2+x^3+x^5+x^6"]);
cenc = comm.ConvolutionalEncoder(TrellisStructure = trellis, ...
    TerminationMethod = "Continuous");
encodedCNAVData = cenc(cnavData);

Based on the configuration, generate the LNAV data.

lnavData = HelperGPSNAVDataEncode(lnavConfig);

Specify all of the required properties for waveform genertion.

CLCodeResetIdx = 75; % CL-code spans over 75 data bits before resetting
numBBSamplesPerDataBit = 204600;
CLCodeIdx = mod(NavDataBitStartIndex-1,CLCodeResetIdx);
IQContent = [IBranchContent,QBranchContent];
pgen = gpsPCode(PRNID = PRNID, InitialTime = ...
    lnavConfig.ReferenceTimeOfEphemeris, ...
    OutputCodeLength = numBBSamplesPerDataBit);
% Pre-initialize the baseband waveform for speed
gpsBBWaveform = zeros(numBBSamplesPerDataBit*NumNavDataBits,1);

Create a file into which the waveform is written.

if WriteWaveformToFile == 1
    bbWriter = comm.BasebandFileWriter("Waveform.bb",10.23e6,0);
end

Independently process each navigation data bit in a loop.

for iDataBit = 1:NumNavDataBits
    dataBitIdx = iDataBit+NavDataBitStartIndex-1;
    bbSamplesIndices = ((iDataBit-1)*numBBSamplesPerDataBit+1): ...
        (iDataBit*numBBSamplesPerDataBit);
    gpsBBWaveform(bbSamplesIndices) = HelperGPSBasebandWaveform(IQContent,pgen,PRNID, ...
        CLCodeIdx,lnavData(dataBitIdx),encodedCNAVData(dataBitIdx));
    CLCodeIdx = mod(CLCodeIdx+1,CLCodeResetIdx);
    if WriteWaveformToFile == 1
        bbWriter(gpsBBWaveform(bbSamplesIndices));
    end
end

Close the file if it is opened.

if WriteWaveformToFile == 1
    release(bbWriter);
end

Signal Visualization

Plot autocorrelation of the C/A-code and visualize the spectrum of the GPS signals.

if ShowVisualizations

Autocorrelation of the ranging code sequence is near-zero except at zero delay, and crosscorrelation of two different sequences is near-zero. Because the C/A-code is periodic with a period of 1023 bits, autocorrelation has a peak for a delay of every 1023 bits. Calculate and plot the autocorrelation of GPS spreading code.

    % Because P-code is 10 times faster than C/A-code or L2 CM-/L2 CL-code,
    % initialise down sample factor to 10
    downsampleFactor = 10;
    IBranchData = real(gpsBBWaveform);
    QBranchData = imag(gpsBBWaveform(1:downsampleFactor:end));
    lags = (-1023:1023).';
    plot(lags,xcorr(real(QBranchData(1:1023)),1023))
    grid on
    xlabel("Number of Samples Delayed")
    ylabel("Autocorrelation")
    title("Autocorrelation of GPS Spreading Code")

Compare the power spectral density of in-phase branch and quadrature-branch signals. This spectrum plot shows that the P-code is wider when P-code is used.

    repeatFactor = 40;
    % Repeat the generated BPSK signal of C/A-code to see the adjacent bands spectrum
    QBranchUpsampled = repmat(QBranchData(:).',repeatFactor,1);
    QBranchUpsampled = QBranchUpsampled(:);
    % Repeat the generated BPSK signal of in-phase component to see the
    % adjacent bands spectrum. Repeat the in-phase branch samples ten times less
    % as every sample in quadrature-branch corresponds to 10 samples in in-phase branch
    IBranchUpsampled = repmat(IBranchData(:).',repeatFactor/10,1);
    IBranchUpsampled = real(IBranchUpsampled(:));
    iqScope = spectrumAnalyzer(SampleRate = 1.023e6*repeatFactor, ...
        SpectrumType = "Power density", ...
        SpectrumUnits = "dBW", ...
        YLimits = [-130, -50], Title = ...
        "Comparison of Power Spectral Density of GPS baseband I and Q Signals", ...
        ShowLegend = true, ChannelNames = ...
        ["Q-branch spectrum with content: " + QBranchContent, ...
        "I-branch spectrum with content: " + IBranchContent]);

    iqScope([QBranchUpsampled,IBranchUpsampled]);

Plot the signal power spectral density at the baseband. To observe the adjacent band spectrum for the GPS signal, repeat the signal at the baseband.

    repeatFactor = 4;
    % Repeat the generated BPSK signal to see the adjacent bands spectrum
    updata = repmat(gpsBBWaveform(:).',repeatFactor,1);
    updata = updata(:);
    bbscope = spectrumAnalyzer(SampleRate = 10*1.023e6*repeatFactor, ...
        SpectrumType = "Power density", ...
        SpectrumUnits = "dBW", ...
        YLimits = [-120,-50], ...
        Title = "Power Spectral Density of Complex Baseband GPS Signal");
    bbscope(updata);
end

Further Exploration

This example uses a configuration object to generate GPS data bits and the navigation signal in the baseband. You can replace the property values in this configuration object and observe how the GPS data is generated. You can also change the ephemeris parameters with an existing real data set. Additionally, you can specify your own almanac file. If you are using your own almanac file, the week number in the almanac file and the week number in the configuration object must match.

Further, this example shows how to generate a GPS baseband waveform, which can be extended to generate an intermediate frequency (IF) waveform from the baseband waveform by multiplying a cosine signal on the in-phase branch and a sine signal on the quadrature-branch.

Additionally, this example shows how to generate a GPS waveform from one satellite, which can be combined along with multiple satellite PRN codes to get an integrated signal.

Appendix

This example uses these data and helper files:

Bibliography

[1] IS-GPS-200, Rev: N. NAVSTAR GPS Space Segment/Navigation User Segment Interfaces. Aug 22, 2022; Code Ident: 66RP1.