The signal of the GPS contains three parts: carrier signals (L1 and L2), spread spectrum sequences (C/A codes, P codes or Y codes) and navigation data (D codes) [4]. The local code is represented as a form of C/A code multiplied by an RF signal. The radio frequency signal RF is a complex number, and the local code is multiplied by the RF code and the C/A code [5]. Assuming that the L1 frequency (1575.42 MHz) is down converted to 21.25 MHz, digitized at a rate of 5 MHz to obtain an output frequency of 1.25 MHz. The frequency range of the capture program is (1250 ± 10) kHz, step 1 kHz, and a total of 21 frequency components. The local code can be expressed as:

In the above equation, the subscript s represents the satellite number, is the C/A code of the satellite S, = 1250 − 10, 1250 − 9, ・・・, 1250 + 10 kHz. The local signal must also be digitized at a 5 MHz sampling rate to produce 5000 data points.

The above 21 sets of data (only corresponding to a satellite C/A code) re- presents 21 frequencies at intervals of 1 kHz. If the C/A code and the frequency of the local signal are correct, and the C/A code phase is aligned, the output is maximized. The flow chart of the acquisition process is as shown in Figure 1:

The operation of the group is as follows: the digitized input signal and the locally generated signal are multiplied by a point (5000 points in total). The square of the product is added to the square of the imaginary part and the square root is the magnitude of the output frequency. Every 200 ns carries on the above operation to the data, after the input data shift 5000 points, has completed the search to the 1ms data, altogether outputs 5000 magnitude values. 21 groups of local signals output 105,000 (5000 × 21) amplitude values. In these amplitude values, the maximum frequency component beyond the threshold corresponds to the Doppler frequency. The maximum value corresponding to the data shift K is the starting position of the C/A code [6].

When the starting position of the C/A code of the GPS signal is 1234 and the carrier frequency is 1,255,123 Hz, setting the signal-to-noise ratio to −10 dB, the simulation results of the acquisition of the No. 1 satellite signal are shown in Figures 2-4 :

This chapter briefly introduces the basic principle and implementation of GPS signal acquisition based on time domain correlation method. Through the simulation results, we can see that this method can realize GPS signal acquisition under certain SNR conditions. However, from the principle and the actual simulation we can see that the method has a lot of shortcomings in the calculation time, because it is in the time domain point-by-point correlation, the calculation is larger. Therefore, it is necessary to propose a faster GPS signal acquisition method, that is, the following GPS signal acquisition method based on FFT cyclic correlation.

The cyclic correlation DFT of the two sequences is equal to the conjugate of one of the sequences DFT multiplied by another sequence DFT [7]. The cyclic correlation method can be used for acquisition, and the computation time can be effectively reduced by fast Fourier transform. The specific steps are as follows:

1) FFT operation is performed on the input data x(n) of 1 ms to obtain X(k). n = k = 0, 1, 2, ・・・, 4999.

2) Take X(k) complex conjugate, get X(k)*.

3) Generate 21 local codes l_si(n), where i = 1, 2, ・・・, 21.

4) Take FFT for l_si(n) and get L_si(k).

5) X(k)* and L_si(k) are multiplied by point to obtain R_si(k).

6) Take R_si(k) inverse FFT transform, get the time domain signal r_si(n), take the absolute value, get |r_si(n)|. A total of 105,000 |r_si(n)|.

7) Give the starting position of the C/A code at a resolution of 200 ns and the carrier frequency at a frequency resolution of 1 kHz at the nth position and the maximum value of the ith frequency |r_si(n)|.

When the starting position of the C/A code of the GPS signal is 1234 and the carrier frequency is 1,255,123 Hz, setting the signal-to-noise ratio to −10 dB, the simulation results of the acquisition of the No. 6 satellite signal are shown in Figures 5-7.

This chapter gives the GPS acquisition method based on FFT cyclic correlation, introduces its basic principle and implementation process, and then gives the simulation results. In the code phase, the fast acquisition algorithm based on FFT requires only one FFT and one IFFT operation to determine the correlation value for all samples in a cycle of C/A codes. Rather than the serial search acquisition method, the code phase is gradually sliding to find the correlation value. And it greatly reduces the amount of calculation.

The fast acquisition algorithm based on FFT in the carrier frequency using a direct search method, but it requires only once FFT and once IFFT operation in the code phase to determine the correlation value for all samples in a cycle of C/A codes. Rather than the serial search acquisition method, the code phase is gradually sliding to find the correlation value.

Assuming that using a C/A code period of data for acquisition, the sampling points is N. When using the traditional serial search acquisition algorithm, each time of the correlation results need N times multiplication and N-1 times addition. To calculate all the N code phase, a total of N² need to multiply. Using the fast acquisition algorithm based on FFT, the total amount of computation required is about 2 times the amount of FFT calculation, that is, times multiplication, reducing the amount of computation [8].

According to the actual simulation results of MATLAB, the average time of GPS signal acquisition method based on time domain correlation method is 21.124146 seconds, while the average time of GPS signal acquisition method based on FFT cyclic correlation method is only 0.261626 seconds.

According to the fact that GPS receiver receives the GPS satellite signal in actual situation, the received power typical value of C/A code is between −162.5 - 154.5 dBW, and the SNR is between −35.5 - −27.5 dB. Therefore, set the SNR range of simulation between −36 - 17 dB, and the starting position of the C/A code of the GPS signal is 1234 and the carrier frequency is 1,255,123 Hz. The simulation results of two acquisition methods of the No. 6 satellite signal are shown in Figure 8.

As can be seen from Figure 9, the time delay error tends to decrease with the increase of SNR in the given SNR range. In the case of small SNR, because the SNR has a great impact on the GPS signal acquisition results, the results are not stable, and the curve has a floating.

It can be seen from the figure that the results obtained by using the FFT method are similar to those of the time domain correlation method. The time delay error tends to decrease with the increase of SNR in the given SNR range. In the case of small SNR, because the SNR has a great impact on the GPS signal acquisition results, the results are not stable, and the curve has a floating. However, it can be seen that the results obtained using the FFT method are more ideal than the time domain correlation method, and the floating is relatively small.