Low Computational Signal Acquisition for GNSS Receivers Using a Resampling Strategy and Variable Circular Correlation Time

For the objective of essentially decreasing computational complexity and time consumption of signal acquisition, this paper explores a resampling strategy and variable circular correlation time strategy specific to broadband multi-frequency GNSS receivers. In broadband GNSS receivers, the resampling strategy is established to work on conventional acquisition algorithms by resampling the main lobe of received broadband signals with a much lower frequency. Variable circular correlation time is designed to adapt to different signal strength conditions and thereby increase the operation flexibility of GNSS signal acquisition. The acquisition threshold is defined as the ratio of the highest and second highest correlation results in the search space of carrier frequency and code phase. Moreover, computational complexity of signal acquisition is formulated by amounts of multiplication and summation operations in the acquisition process. Comparative experiments and performance analysis are conducted on four sets of real GPS L2C signals with different sampling frequencies. The results indicate that the resampling strategy can effectively decrease computation and time cost by nearly 90–94% with just slight loss of acquisition sensitivity. With circular correlation time varying from 10 ms to 20 ms, the time cost of signal acquisition has increased by about 2.7–5.6% per millisecond, with most satellites acquired successfully.


Introduction
As the most mature Global Navigation Satellite System (GNSS), Global Positioning System (GPS) has occupied the dominated position in many aspects, such as civil application, scientific research, and military defense. Further, the GPS modernization plan has drawn remarkable improvement that adds new civil signals transmitted by parts of GPS satellites for high-accuracy navigating and positioning in complex environments. These new signals possess better performance of self-/cross-correlation, as well as the ability to suppress narrowband interference and correct transmission delay of the ionosphere [1]. The longer code period of new signals, especially GPS L2C signals, indicates that computational complexity and time consumption will become greater for conventional acquisition algorithms [2]. Besides, multimode multi-frequency GNSS receivers are applied in many fields nowadays. Broadband receiving of multi-frequency signals is the popular and feasible technique to minimize the energy consumption and physical size of those receivers [3,4]. However, broadband receiving requires higher sampling frequency, which dramatically increases computation and time consumption of the subsequent signal processing. It brings an inevitable dilemma for signal acquisition.
In principle, signal acquisition is the chief operation to identify the visibility of satellites and provide estimated values of carrier frequency and code phase of received signals. Based on these estimated values, signal tracking can be successfully activated for following baseband operations. Actually, signal acquisition is a two-dimensional search process over carrier frequency and code phase uncertainty [5], by correlating with local replicated ones. Successful acquisitions are affirmed if acquisition decision statistics exceed the acquisition threshold. Conceivably, signal acquisition is quite time-consuming and requires tremendous computation [6]. Meanwhile, computation acts as the bottleneck of GNSS receiver applications, as it determines the lowest sensitivity of the baseband processing. Therefore, requirements for faster, less computationally demanding, and more efficient acquisition approaches are put forward to GNSS receivers in an extensive research field [1, 5,7].
Attributed to the primary convolution and correlation operations in the acquisition process, various sophisticated acquisition algorithms have been investigated to enhance the efficiency of Fast Fourier Transform (FFT) to confront the challenge of computation and time cost. For instance, ideas of sparse FFT [8], fast FFT [9], reduced-size FFT [10], and averaging correlation [11] are attempted to reduce computational complexity when acquiring satellite signals. Specific to n-point FFT operations of parallel search acquisition algorithms, various divisions have presented their respective effects on relieving the burden of computation and time cost, such as radix-2 and radix-4 FFT [12] or by replacing with m sets 2k-point FFT [13]. Besides, multi-step acquisition methods [14][15][16][17][18] also are studied to seek for better performance of signal acquisition on computation and time consumption. Taking advantage of the double-channel structure of modernized GPS signals, the joint data-pilot channel strategy [19] outperforms other individual channel strategies on signal acquisition. In fact, computational complexity and time consumption of signal acquisition are highly dependent on the data size of sampled signals used for convolution and correlation operations. Although extensive study has been conducted, most of them concentrate on reducing instant processing points of FFT. Only a few investigations [20][21][22] have carried out attempts to directly shorten or down-sample the received original signals. The bandpass sampling theory [23], which is a direct expansion of the traditional Shannon sampling theorem, has been widely applied in the fields of communication and radar signal processing [24][25][26]. In the GNSS community, the bandpass sampling theory was adopted to design a direct conversion of GNSS carrier signal [27][28][29][30] in order to reduce the sampling frequency of the direct RF sampling front end of GNSS receivers. However, currently mainstream commercial GNSS receivers generally employ the intermediate frequency (IF) sampling front end. In the design of IF sampling front end for multi-frequency GNSS receivers, its bandwidth is generally much larger than that of the single-frequency front end, in order to incorporate multiple signals located at adjacent frequency points in one channel of the IF front end. Owing to the fact that signals of each frequency point are acquired separately, signal acquisition can be realized within the lower bandwidth. It exactly inspires us to apply the bandpass sampling theory to the efficient and fast signal acquisition of the IF sampling GNSS receivers.
As another important solution contributing to the signal acquisition challenge, massive correlation-related approaches also have been widely investigated. From the perspective of correlation period, Qaisar et al. [31] introduced an interesting chip-wise correlation strategy that accumulated code samples in one period and calculated correlation values for each chip period to reduce the search space and computation. In contrast, Zhu et al. [32] extended the correlation time to an entire duration of navigation data bit, aiming at enhancing the sensitivity of receivers. Further, Zhu et al. [33] proposed a variable time parallel acquisition scheme and different accumulation schemes over multiple code periods to detect weak signals. Jin et al. [34] exploited a fusion acquisition method of multiple correlation strategies to reduce correlation loss caused by sign transition of high data rate signals. Especially, a small increase of computation emerged using the delay-accumulation method [35] and data-pilot optimized combining acquisition methods [36]. Additionally, Ta et al. [37] presented a partial acquisition algorithm using specially-designed matched filters and differential post-correlation techniques to simplify computational complexity of signal acquisition. Although increasing correlation time is an effective way to ensure acquisition accuracy, computation and time consumption generally present higher expense due to the longer correlation time. Further, it is difficult to find out many researches that have done detailed exploration on the relationship between acquisition performance and correlation time, and that have been applied effectively to signal acquisition of GNSS receivers.
With regard to broadband multi-frequency GNSS receivers, we propose a resampling strategy and variable circular correlation time in this paper to decrease computational complexity and time consumption of signal acquisition. Inspired by the bandpass sampling theory, the resampling strategy reduces greatly the data size of the sampled signal with just slight loss of signal quality. It lays the computational foundation for FFT, correlation and convolution operations of signal acquisition. In particular, the resampling strategy is implemented step-by-step to demonstrate the capability of combining with existing FFT-/correlation-based acquisition algorithms. Taking the longer code period into consideration, the relationship between circular correlation time and acquisition performance is concerned for the potential of shortening circular correlation time to adapt to different signal strength conditions and thereby increase the operation flexibility of GNSS signal acquisition. Moreover, the acquisition threshold is defined as the ratio of the highest and second highest correlation results in the search space of carrier frequency and code phase. Computational complexity is formulated by amounts of multiplication and summation operations in the signal acquisition.
This paper is organized as follows: Section 2 introduces the characteristics of the GPS L2C signal as an example, and the framework of GNSS receivers to indicate the crucial role of the proposed strategies. Next, Section 3 elaborates on the proposed resampling strategy and variable circular correlation time, in terms of the principle, realization, coarse and fine acquisition, as well as performance evaluation. Comparative experiments and performance analysis are performed on four sets of real GNSS signals in Section 4, to verify the effectiveness of the resampling strategy and variable circular correlation time. Conclusions are drawn in Section 5.

Characteristics of Satellite Signals
In this paper, we take the GPS L2C signal as an example of GNSS signals. GPS L2C signal is the civil signal, newly added on the existing L2 frequency band and modulated by civil-moderate (CM) and civil-long (CL) codes. The mechanism of the GPS L2C signal is described as follows: Firstly, the navigation message is encoded by forward error correction (FEC) coder and modulated on CM code to comprise the data channel. Meanwhile, without navigation message, CL code acts as the pilot channel by direct sequence spread spectrum (DSSS) method. After that, signals of data and pilot channels are separately added on the baseband signal by time division multiplexing (TDM) method, and further modulated by binary phase shift keying (BPSK) method on carrier wave. Finally, the generated signal is transmitted by satellites. GPS L2C signal received from satellites can be modeled as [4,22]: where A is the signal amplitude and D(t) is the navigation message. C M (t) and C L (t) represent CM and CL codes, respectively. θ is code phase in the received signal. T M denotes the period of CM code and k indicates the number of CM epochs in the current CL epoch, 0 ≤ k ≤ 74. The carrier frequency of the received signal consists of the intermediate frequency f IF and Doppler frequency shift f D , which results from the relative motion between the satellite and the receiver. ϕ refers to the initial carrier phase. Besides, n(t) is the thermal noise in the received signal, denoted as the additive white Gaussian noise. It is noticeable that the signal is just composed of the cosine term on account of the BPSK method.

Framework of GNSS Receivers
Primarily, the signal processing framework of GNSS receivers mainly includes 4 parts: data collection, signal acquisition, signal tracking, and navigation solution, as shown in Figure 1. Firstly, Sensors 2018, 18, 678 4 of 21 the original signal received by the antenna is converted from low-power high-frequency into intermediate-frequency, and then sampled as digital to send to different channels of receivers. Next, visible satellites are identified by signal acquisition. It is also obtained that estimated values of carrier frequency and code phase of acquired satellite signals. Then, based on these estimated values, carrier wave and pseudo-random noise (PRN) code are accurately removed from the received signal by applying phase locked loop (PLL) and delay locked loop (DLL) to signal tracking. Finally, the navigation message is extracted by bit and frame synchronization to calculate position and velocity of the receiver. In this way, satellite signal has been processed by GNSS receivers to accomplish the functions of positioning, navigating, and timing.

Framework of GNSS Receivers
Primarily, the signal processing framework of GNSS receivers mainly includes 4 parts: data collection, signal acquisition, signal tracking, and navigation solution, as shown in Figure 1. Firstly, the original signal received by the antenna is converted from low-power high-frequency into intermediate-frequency, and then sampled as digital to send to different channels of receivers. Next, visible satellites are identified by signal acquisition. It is also obtained that estimated values of carrier frequency and code phase of acquired satellite signals. Then, based on these estimated values, carrier wave and pseudo-random noise (PRN) code are accurately removed from the received signal by applying phase locked loop (PLL) and delay locked loop (DLL) to signal tracking. Finally, the navigation message is extracted by bit and frame synchronization to calculate position and velocity of the receiver. In this way, satellite signal has been processed by GNSS receivers to accomplish the functions of positioning, navigating, and timing. In this paper, the resampling strategy is proposed to decrease the data size of the sampled signal without obvious loss of signal quality, so as to greatly reduce computational complexity and time consumption of signal acquisition. Currently, signal acquisition is fundamentally realized by the serial search algorithm, parallel frequency space search algorithm, and parallel code phase search algorithm. By contrast, the parallel code phase search algorithm is the most efficient one that searches all code phases at the same time for each frequency bin. Hence, we chose the parallel code phase search algorithm on behalf of conventional acquisition algorithms to perform the resampling strategy. Besides, signal acquisition is improved by variable circular correlation time. In this paper, the resampling strategy is proposed to decrease the data size of the sampled signal without obvious loss of signal quality, so as to greatly reduce computational complexity and time consumption of signal acquisition. Currently, signal acquisition is fundamentally realized by the serial search algorithm, parallel frequency space search algorithm, and parallel code phase search algorithm. By contrast, the parallel code phase search algorithm is the most efficient one that searches all code phases at the same time for each frequency bin. Hence, we chose the parallel code phase search algorithm on behalf of conventional acquisition algorithms to perform the resampling strategy. Besides, signal acquisition is improved by variable circular correlation time.

Principle of the Resampling Strategy
Generally, the original signal can be reconstructed from the sampled signal without loss of signal quality if the sampling frequency is higher than twice the upper cutoff frequency of the original signal. However, this principle is not suitable for signal acquisition of multimode multiband GNSS receivers. These receivers can receive multi-frequency satellite signals at the same time. Therefore, a much higher sampling frequency is required to sample the received broadband signal. Obviously, it is rather aimless, low-efficiency, and of low quality to acquire a specific frequency signal by sampling the received broadband signal in this way. Moreover, the higher sampling frequency results in much more computation and time consumption for signal acquisition, because the data size of the sampled signal is much larger. To address this dilemma of balance between the sampling frequency and the greatly increasing computation, we propose the resampling strategy based on the bandpass sampling theory [23].
Theorem 1. Bandpass Sampling Theory: Assume that the frequency band of a time-continuous signal where f l and f u are the lower and upper cutoff frequencies, respectively. B denotes the bandwidth of the signal x(t), B = f u − f l . In order to reconstruct the original signal x(t) without spectrum aliasing from the sampled sequence X[k], the acceptable sampling frequency f s should be in the range of where N = Z − ( f u /B) denotes the maximum integer no more than f u /B. Note that, n = 1 just makes sense for the case that f s ∈ [2 f u , +∞).
The relationship between f s and f u is illustrated in Figure 2. In order to avoid spectrum aliasing when reconstructing the original signal, the sampling frequency f s should be chosen from cyan areas, not gray ones. In Figure 2, the available range of the sampling frequency f s consists of N non-overlapping cyan sections, denoting as Φ n = [2 f u /n, 2 f l /(n − 1)], (n = 1, 2, . . . , N) and f s ∈ ∪ N n=1 Φ n . It is realizable for the sampling frequency to be much lower, such as [2B, 4B], instead of the conventional range ( f s ≥ 2 f u , the cyan n = 1). The parameters, f su , f sl , and f ds will be discussed for the proposed resampling strategy later.

Principle of the Resampling Strategy
Generally, the original signal can be reconstructed from the sampled signal without loss of signal quality if the sampling frequency is higher than twice the upper cutoff frequency of the original signal. However, this principle is not suitable for signal acquisition of multimode multiband GNSS receivers. These receivers can receive multi-frequency satellite signals at the same time. Therefore, a much higher sampling frequency is required to sample the received broadband signal. Obviously, it is rather aimless, low-efficiency, and of low quality to acquire a specific frequency signal by sampling the received broadband signal in this way. Moreover, the higher sampling frequency results in much more computation and time consumption for signal acquisition, because the data size of the sampled signal is much larger. To address this dilemma of balance between the sampling frequency and the greatly increasing computation, we propose the resampling strategy based on the bandpass sampling theory [Error! Reference source not found.].
The relationship between and is illustrated in Figure 2. In order to avoid spectrum aliasing when reconstructing the original signal, the sampling frequency should be chosen from cyan areas, not gray ones. In Figure 2, the available range of the sampling frequency consists of non-overlapping cyan sections, denoting as . It is realizable for the sampling frequency to be much lower, such as [2 , 4 ], instead of the conventional range ( ≥ 2 , the cyan = 1). The parameters, , , and will be discussed for the proposed resampling strategy later.

Realization of the Resampling Strategy for Signal Aquisition
Sampling the time-continuous signal s(t) with the frequency f s , the frequency spectrum of the sampled sequence S(kT s ) is the same as that of the original signal s(t) in the frequency domain, periodically extended with the period f s . Hence, with the appropriate sampling frequency, expanding components of the sampled sequence will not fold with the frequency spectrum of the original signal In this way, the original signal can be recovered from the sampled sequence without loss of signal quality. Based on the bandpass sampling theory, the proposed resampling strategy is realized as Algorithm 1 for signal acquisition of GNSS receivers.

Input:
The received broadband satellite signal;

Output:
Doppler frequency and code phase offset of the received satellite signal.

1.
Filter out side lobes of the original broadband satellite signal; 2.
Update the parameters of the resampling strategy; 3.
Resample the main lobe signal with the resampling frequency; 4.
Restore the resampling acquisition results to that of conventional acquisition algorithms; Return: the acquired Doppler frequency and code phase offset of the received satellite signal.
Note the broadband signal received by multimode multiband GNSS receivers as It includes GPS L1, L2, L5 signals, as well as signals transmitted by other satellite systems, such as Galileo, GLONASS, and BeiDou. We intend to assume that the main lobe refers to the GPS L2C signal, and other signals are regarded as the side lobe signals. A bandpass filter is designed to extract the main lobe signal S L2 (t) from the broadband signal S(t).
where f l and f u are the lower and upper cutoff frequencies of the main lobe signal S L2 (t). f s is the conventional sampling frequency for the broadband signal S(t). f IF and B represent the intermediate frequency and bandwidth of S L2 (t). It is guaranteed that f s 2 f u in order to effectively sample S L2 (t). The bandpass filter is designed as an n-order Hamming-window linear-phase finite impulse response (FIR) filter. The pass band ω p is the vector of the normalized cutoff frequencies. ∆ω f is the margin to ensure the main lobe signal reserved completely.
By employing the designed bandpass filter on the broadband signal, the side lobe signals are filtered out while the main lobe remains. This plays as a foundational block to avoid frequency aliasing and reduce the data size of the sampled signal.
(1) Update the Parameters of the Resampling Strategy A resampling filter is designed to resample the main lobe signal S L2 (t), obtained by the bandpass filter. The lower and upper cutoff frequencies and bandwidth of the resampling filter, f dl , f du and B d , are expressed as Based on the bandpass sampling theory in Theorem 1, the lower and upper boundaries of the acceptable sampling frequency, f sl and f su , are determined as Here, n = 1, 2, . . . , N. And N = Z − ( f du /B d ). The center of the acceptable sampling frequency range [ f sl , f su ] is regard as the resampling frequency f ds : The parameters of the resampling strategy, f sl , f su and f ds , are illustrated as the blue, red and green curves in Figure 2, respectively. The number n is chosen as the maximum integer N to decrease the resampling frequency as much as possible. For that S L2 (t) is the high-frequency narrowband signal, the conventional sampling frequency f s is much higher than the bandwidth ( f s ≥ 2 f u B). As the red curve in Figure 2, the resampling frequency f ds is maintained in the lower frequency range (

2) Resample the Main Lobe Signal with the Resampling Strategy
When applying the resampling filter to the main lobe signal S L2 (t), the sampling frequency and intermediate frequency are updated to f ds and f dIF . Especially, f ds is significantly decreased compared to the conventional one f s . The original intermediate frequency f IF also is equivalent to the much lower one f dIF , which is the remainder of f IF and f ds , as Equation (8). The resampled sequence S d (k) is expressed as Equation (9).
where Z + (x) denotes the minimum integer bigger than x. S s (·) and S d (·) are the sequences obtained by sampling S L2 (t) with f s and f ds , respectively. i is the index of samples in the sequence S s (i), i = 1, 2, . . . , L. L denotes the length of S s (i), while k is that of the resampled sequence S d (k), k = 1, 2, . . . , L d . Due to f ds f s , L d is much shorter than L (L d L). Therefore, S d (k) is extracted from S s (i) to equivalently express the main lobe signal S L2 (t). Though the length is greatly reduced, S d (k) is sufficient to reconstruct the main lobe signal. It is beneficial for reducing the computation and time consumption of signal acquisition.
(3) Restore the Resampling Acquisition Results to that of Conventional Acquisition Algorithms The parameters of the resampling strategy, f ds and f dIF , and the resampled sequence S d (k) are utilized to obtain the acquisition results, including actual carrier frequency f dc , Doppler frequency shift f dD , and code phase offset θ d . In order to be applied to signal tracking in the conventional way, the acquisition results of the resampling strategy are restored equivalently to the carrier frequency f c and code phase offset θ of the conventional acquisition algorithm. In this way, the capability of the proposed resampling strategy to combine with conventional acquisition algorithms has been proven.
With the resampling strategy acting on the original broadband signal, the signal flow chart is illustrated in Figure 3. Remarkably, the data size of the resampled signal is decreased from L to . Therefore, the resampling strategy is capable to greatly reduce computation and time consumption of signal acquisition.
With the resampling strategy acting on the original broadband signal, the signal flow chart is illustrated in Figure 3. Remarkably, the data size of the resampled signal is decreased from to , ( ≪ ). Therefore, the resampling strategy is capable to greatly reduce computation and time consumption of signal acquisition.

Coarse Acquisition with Variable Circular Correlation Time
Conventional acquisition methods usually set the circular correlation time as one period of PRN code. Those are low-efficiency, time-consuming, and of large computation. Especially, code periods of GPS L2C signal are as long as 20 ms for CM code and 1.5 s for CL code. Actually, the shorter circular correlation time could quicken signal acquisition if the received GPS signal is strong enough. For weak signals in noisy environments, the longer circular correlation time ensures the acquisition performance. For the purpose of exploring the effects of different circular correlation times on the performance of signal acquisition, we introduce variable circular correlation time in the coarse acquisition.
The reference signal ( ) of length is extracted from the resampled signal ( ). The zero-padding PRN code ( ) is introduced to calculate circular correlation with the reference signal to obtain accurate carrier frequency and code phase.

Coarse Acquisition with Variable Circular Correlation Time
Conventional acquisition methods usually set the circular correlation time as one period of PRN code. Those are low-efficiency, time-consuming, and of large computation. Especially, code periods of GPS L2C signal are as long as 20 ms for CM code and 1.5 s for CL code. Actually, the shorter circular correlation time could quicken signal acquisition if the received GPS signal is strong enough. For weak signals in noisy environments, the longer circular correlation time ensures the acquisition performance. For the purpose of exploring the effects of different circular correlation times on the performance of signal acquisition, we introduce variable circular correlation time in the coarse acquisition.
The reference signal S r (i) of length L r is extracted from the resampled signal S d (i). The zero-padding PRN code C r (i) is introduced to calculate circular correlation with the reference signal to obtain accurate carrier frequency and code phase. Specifically, P c is the period of PRN code and P x represents variable circular correlation time, P x ∈ (0, P c ]. If P x = P c , the variable circular correlation time method is the same as conventional acquisition algorithms. L r , L c and L x are the lengths of signals with periods P r , P c , and P x . At frequency search step k, the reference signal S r (i) is multiplied with local carrier replica W l (k) of frequency f (k) to remove carrier wave from S r (i) and obtain the baseband signal S w (k) With different code phase offsets, the baseband signal S w (k) is circular correlated with the local zero-padding code C r (p), as illustrated in Figure 4. r x (k, p) is the circular correlation result of S w (k) and C r (p). FFT is adopted to quicken the correlation calculation in the time domain, as Equation (20).
Specifically, is the period of PRN code and represents variable circular correlation time, ∈ (0, ]. If = , the variable circular correlation time method is the same as conventional acquisition algorithms. , and are the lengths of signals with periods , , and . At frequency search step , the reference signal ( ) is multiplied with local carrier replica ( ) of frequency ( ) to remove carrier wave from ( ) and obtain the baseband signal ( ) With different code phase offsets, the baseband signal ( ) is circular correlated with the local zero-padding code ( ), as illustrated in Figure 4.   . Circular correlation results of the baseband signal S w (k) and the local zero-padding code C r (p): (a) code offset between S w (k) and C r (p) is 0 (aligned); (b) code offset between S w (k) and C r (p) is less than (L c − L x ) samples; (c) code offset between S w (k) and C r (p) is more than (L c − L x ) samples.
Correlation magnitude reflects the self-/cross-correlation relationship of the received satellite signal and the local code replica. When they are completely aligned, a much higher correlation peak is gained. In the search space of carrier frequency and code phase, the correlation magnitude M x (k, p) is formulated as where R(·) and I(·) denote the real and imaginary parts of a complex data. By Equations (19)- (21), it is very easy to find out the highest correlation peak in the search space. Denote the corresponding carrier frequency and code phase asf dc andθ d , respectively. In this paper, acquisition threshold is defined as the ratio of the highest and second highest correlation values in the search space of carrier frequency and code phase. The second highest peak is found out in the neighborhood range of the highest one. The left neighborhood range is obtained as where ∆ θ is the one-chip-period range, excluded from the neighborhood range, for it is too close to the highest peak to find out the useful second peak. If the correlation ratio is bigger than the acquisition threshold, it is confirmed that the corresponding satellite is visible and the signal of frequencyf dc and code phaseθ d is acquired successfully.
Note that the frequency spectrum of C r (p) is asymmetrical, such that the former half part occupies most of signal power. Thus, we just use the former to compute circular correlation results. This reduces FFT samples by half, so as to decrease computation of signal acquisition. Additionally, L x is highly related to the amount of FFT samples for signal acquisition. Equation (19) provides a good indication to reduce computation and time cost of signal acquisition by decreasing the sampling frequency f s and circular correlation time P x . That is exactly realized by the proposed resampling strategy and variable circular correlation time.

Fine Acquisition with Pilot Channel
For GPS L2C signals, the phase of CL code is aligned with that of the CM code. The reference signal S f r (i) of length L x is directly extracted from the main lobe signal S L2 (t) with code phaseθ d . Similarly, the CL code is removed from S f r (i) by multiplying with the local CL code replica to get the baseband signal S f x (i). Considering that no navigation message is modulated on the signal of pilot channel, FFT acts as the frequency discriminator to find out the accurate carrier frequency.
where F [S f x (i)] is the FFT result of S f x (i). The amount of FFT samples N L is set as the power of 2 closest to L x to speed up FFT calculation, when the signal length L x is not an exact power of 2.
As F [S f x (i)] is a periodic complex even signal, the frequency spectrum is symmetrical in the frequency range [0, f ds ]. Hence, the one-sided FFT is adopted.
Assume M f x ( m) as the maximum magnitude of F [S f x (i)] in the frequency range [0, f ds ], andm is the corresponding index of F [S f x (i)]. The high-accuracy carrier frequency of the received signal is discriminated as f dc in the fine acquisition: Note that the frequency resolution f ds /N L is inversely proportional to N L , and N L is highly dependent on L x . From this perspective, circular correlation time is positively relative to the resolution of acquired carrier frequency. That is, the longer circular correlation time will ultimately result in higher accurate carrier frequency in the fine acquisition. Besides, circular correlation time also affects computation and time cost of the fine acquisition, for that FFT is utilized as the frequency discriminator.

Performance Evaluation of Signal Acquisition
For signal acquisition, time cost is usually a stochastic variable. Mean value of time cost is concerned as one of the most important measurements to assess the performance of acquisition algorithms. Computation complexity is another crucial evaluation index. Usually, convolution operation in the time domain is converted into multiplication in the frequency domain by FFT to improve the efficiency of signal acquisition algorithms. We choose amounts of multiplication and summation operations to evaluate the computational complexity of signal acquisition.
Equation (25) Note the computations of local code replica generation and correlation ratio calculation as O lc and O pk . And the computations of coarse and fine acquisition are denoted as O ca and O f a , respectively. N SV is the number of satellites in the search list, and N f is the amount of carrier frequency bins.
Therefore, the computational complexity of signal acquisition O total is formulated as (26) in which,

Experimental Platform and Datasets Description
The resampling strategy and variable circular correlation time proposed in this paper are realized by MATLAB software, configured with Intel dual-core 3.5 GHz CPU of i7 − 5930K and 16.0 G RAM. The GNSS receiver is equipped with NT1065 "Nomada" and bladeRF as the RF front end. They can simultaneously receive various GNSS satellite signals, including GPS (L1, L2, L3, L5), GLONASS (E1, E5a, E5b, E6), BeiDou (B1, B2, B3), Galileo, IRNSS, QZSS. For comparison, experiments and performance analysis are conducted on 4 sets real GNSS signals collected at the same observation station. Different configurations are adopted, such as set-up status of the resampling strategy (on/off) and circular correlation time (10-20 ms). Table 1 presents the data type, intermediate frequency, and conventional sampling frequency of experimental datasets. Besides, basic parameters of the experimental GNSS signals are provided here: the bandwidth of main lobe is 2.046 MHz, the basis carrier frequency is 1227.6 MHz, the code frequency is 511.5 KHz, the range of frequency search band is ±5 KHz, and the bandwidth of frequency bins is 100 Hz. Acquisition threshold is set as correlation ratio R c = 1.5.

Performance Analysis of the Resampling Strategy
With emphasis on effects of the proposed strategy on signal acquisition, we execute the resampling strategy and the conventional algorithm (without the resampling strategy) on all datasets. Experimental results are discussed in terms of acquisition effectiveness, sensitivity, computation and time cost. Figure 5 presents the characteristics of the received broadband signal in the frequency and time domains, as well as the amplitude distribution of the received signal in Dataset 1. By applying the bandpass filter to the received broadband signal, the power of the side lobe signals was filtered out and the main lobe signal was retained, as shown in Figure 6. The bandwidth and intermediate frequency of the main lobe signal were the same as that of the received broadband signal, B d = 2.046 MHz and f IF = 7.4 MHz. Besides, the sampling frequency of the main lobe signal was the conventional one, f s = 53 MHz. Applying the resampling strategy to the main lobe signal, the resampled signal was obtained, as presented in Figure 7. Specially, the resampled signal had the same bandwidth of the main lobe signal, while the intermediate frequency changed from f IF = 7.4 MHz to f dIF = 1.43 MHz. Meanwhile, as a result of the resampling strategy, the sampling frequency was greatly reduced from the conventional one f s = 53 MHz to f ds = 5.97 MHz, to reduce the data size of the resampled sequences used in FFT, IFFT, and correlation calculation. Here, circular correlation time of signal acquisition was fixed as L x = 20 ms.          The acquisition results of Dataset 1 are listed in Table 2. There were 7 satellites ( PRN 5,6,12,17,24,25,and 29) acquired successfully. It also has proved the capability of the proposed strategy to adapt to the signals with different carrier-to-noise ratios. Comparing experimental results of the resampling strategy and the conventional method (without the resampling strategy), it can be found that the two algorithms demonstrated almost the same acquisition performances, including acquired satellites, carrier frequency, Doppler frequency shift, and code phase offset. There is no obvious loss of acquisition performance for the resampling strategy. In addition, it is remarkable that the correlation magnitude of each acquired satellite decreased to about 10% as a result of the resampling strategy, shown in the bold columns in Table 2. This can be qualitatively attributed to the reduced-size samples for circular correlation by the resampling strategy, see Equation (10). Moreover, the correlation ratio of signal acquisition reflects the relationship of the satellite signal and local PRN replica, acting as the acquisition threshold in this paper. From the bold rows listed in Table 2, the correlation magnitude of Satellite PRN12 was a little smaller than that of Satellite PRN6 when acquired with the resampling strategy. Nevertheless, the correlation ratio of Satellite PRN12 was still the biggest among all acquired satellites. The acquisition results of the resampling strategy kept consistent with that of the conventional acquisition algorithm in Table 2. It proves that setting the correlation ratio as acquisition threshold can enhance the acquisition performance, especially for cases: (1) correlation magnitudes are too small in the search space due to noise in the received broadband signal; (2) there exists more than one value which is higher than the conventional correlation threshold (always a constant of correlation magnitude).
Satellites acquired by the resampling strategy and the conventional acquisition algorithm were completely the same, as illustrated in Figure 8. This has proved the effectiveness of the proposed resampling strategy for signal acquisition. Furthermore, for acquired (cyan and bright green bars) and not acquired (dark blue and blackish green bars) satellites, the correlation magnitude of the resampling strategy (the left Y axis) was about 10% that of the conventional one (the right Y axis), corresponding to results in Table 2. The acquisition results of Dataset 1 are listed in Table 2. There were 7 satellites ( 5, 6, 12, 17, 24, 25, and 29) acquired successfully. It also has proved the capability of the proposed strategy to adapt to the signals with different carrier-to-noise ratios. Comparing experimental results of the resampling strategy and the conventional method (without the resampling strategy), it can be found that the two algorithms demonstrated almost the same acquisition performances, including acquired satellites, carrier frequency, Doppler frequency shift, and code phase offset. There is no obvious loss of acquisition performance for the resampling strategy. In addition, it is remarkable that the correlation magnitude of each acquired satellite decreased to about 10% as a result of the resampling strategy, shown in the bold columns in Table 2. This can be qualitatively attributed to the reduced-size samples for circular correlation by the resampling strategy, see Equation (10).
Moreover, the correlation ratio of signal acquisition reflects the relationship of the satellite signal and local PRN replica, acting as the acquisition threshold in this paper. From the bold rows listed in Table 2, the correlation magnitude of Satellite PRN12 was a little smaller than that of Satellite PRN6 when acquired with the resampling strategy. Nevertheless, the correlation ratio of Satellite PRN12 was still the biggest among all acquired satellites. The acquisition results of the resampling strategy kept consistent with that of the conventional acquisition algorithm in Table 2. It proves that setting the correlation ratio as acquisition threshold can enhance the acquisition performance, especially for cases: (1) correlation magnitudes are too small in the search space due to noise in the received broadband signal; (2) there exists more than one value which is higher than the conventional correlation threshold (always a constant of correlation magnitude).
Satellites acquired by the resampling strategy and the conventional acquisition algorithm were completely the same, as illustrated in Figure 8. This has proved the effectiveness of the proposed resampling strategy for signal acquisition. Furthermore, for acquired (cyan and bright green bars) and not acquired (dark blue and blackish green bars) satellites, the correlation magnitude of the resampling strategy (the left Y axis) was about 10% that of the conventional one (the right Y axis), corresponding to results in Table 2.  Doppler frequency shift, code phase offset, and correlation ratio of acquired satellites in Table 2 can be demonstrated clearly in the search space of frequency and code phase, for instance, the acquisition result of Satellite PRN12 shown in Figure 9. The correlation peak occurred at the point of frequencŷ f = 7.39853 MHz and code phaseθ = 193394 (samples). The highest peak was much bigger (R c > 7) than other correlation magnitudes in the search space. Therefore, Satellite PRN 12 was acquired successfully, the carrier frequency is 7.39853 MHz and code phase is 0.365π (1866 chips). The acquisition results were adopted as estimated parameters for the subsequent signal tracking. Doppler frequency shift, code phase offset, and correlation ratio of acquired satellites in Table 2 can be demonstrated clearly in the search space of frequency and code phase, for instance, the acquisition result of Satellite PRN12 shown in Figure 9. The correlation peak occurred at the point of frequency = 7.39853 MHz and code phase = 193394 (samples). The highest peak was much bigger (R > 7) than other correlation magnitudes in the search space. Therefore, Satellite PRN 12 was acquired successfully, the carrier frequency is 7.39853 MHz and code phase is 0.365π (1866 chips). The acquisition results were adopted as estimated parameters for the subsequent signal tracking. Aimed at investigating the sensitivity of the resampling strategy under weak signal conditions, extensive Monte Carlo simulations were performed by employing a high-confidence software IF signal simulator [38]. The carrier-to-noise ratio was set to vary from 25 to 45 dB-Hz with a step of1 dB-Hz and circular correlation time of signal acquisition was set as = 20 ms. The experiments of signal acquisition were repeated for 200 times for each carrier-to-noise ratio. Detection probabilities are evaluated and compared for signal acquisition without/with the resampling strategy. From the experimental results of sensitivity illustrated in Figure 10, it can be found that the signal could be acquired successfully with the resampling strategy when the carrier-to-noise ratio was bigger than 32 dB-Hz, and detection probability had reached 100% with the carrier-to-noise ratio of 38 dB-Hz. Compared to the detection probability without the resampling strategy (the blue plot), the sensitivity of the resampling strategy had a loss of nearly 1 dB. It can be attributed to the bandpass sampling filter that filters out some powers of the side lobe signals from the received broadband signal. However, the significant computational cost reduction obtained can make up for the sensitivity loss, as longer circular correlation time and non-coherent integration can be adopted to increase the sensitivity. Aimed at investigating the sensitivity of the resampling strategy under weak signal conditions, extensive Monte Carlo simulations were performed by employing a high-confidence software IF signal simulator [38]. The carrier-to-noise ratio was set to vary from 25 to 45 dB-Hz with a step of 1 dB-Hz and circular correlation time of signal acquisition was set as L x = 20 ms. The experiments of signal acquisition were repeated for 200 times for each carrier-to-noise ratio. Detection probabilities are evaluated and compared for signal acquisition without/with the resampling strategy. From the experimental results of sensitivity illustrated in Figure 10, it can be found that the signal could be acquired successfully with the resampling strategy when the carrier-to-noise ratio was bigger than 32 dB-Hz, and detection probability had reached 100% with the carrier-to-noise ratio of 38 dB-Hz. Compared to the detection probability without the resampling strategy (the blue plot), the sensitivity of the resampling strategy had a loss of nearly 1 dB. It can be attributed to the bandpass sampling filter that filters out some powers of the side lobe signals from the received broadband signal. However, the significant computational cost reduction obtained can make up for the sensitivity loss, as longer circular correlation time and non-coherent integration can be adopted to increase the sensitivity.
In order to validate the effects of the proposed strategy on computation and time cost of signal acquisition, comparative experiments without/with the resampling strategy were conducted on all the datasets. Experimental results are presented in Table 3 and Figure 11. Meanwhile, the circular correlation time was fixed as L x = 20 ms to eliminate the effects of variable circular correlation time. With respect to the comparison in Table 3, the resampling frequency was reduced to about 10% of the conventional sampling frequency by applying the resampling strategy to Dataset 1, 2, and 3. In Figure 11, computations (multiplication and summation operations) of all the datasets are linearly related to the amount of frequency bins and the data size of the sampled signal. The amount of frequency bins was set as a constant, based on the width of frequency search range and the bandwidth of each frequency bin. And the data size of the sampled signal was exactly determined by the resampling frequency. Resulting from the reduced sampling frequency, the computation cost of signal acquisition with the resampling strategy was greatly decreased, nearly 90-94% of that without the proposed strategy, as shown in Table 3 and Figure 11. Besides, the received broadband signal in Dataset 4 was an 8-bit complex signal with negative intermediate frequency. The bandpass filter of the proposed resampling strategy was removed since it was not applicable to extract the main lobe signal of Dataset 4. Additionally, the conventional sampling frequency of Dataset 4 was low enough. Thus, the resampling frequency reduced to about 59%, not as nearly 10% for other datasets. Correspondingly, the computation reduction of 43% was gained by applying the resampling strategy to signal acquisition of Dataset 4. Meanwhile, for all experimental datasets, the time cost of signal acquisition with the resampling strategy reduced to about 7.7 ∼ 65% of that with the conventional acquisition algorithm, accompanying with the reduced computation. In order to validate the effects of the proposed strategy on computation and time cost of signal acquisition, comparative experiments without/with the resampling strategy were conducted on all the datasets. Experimental results are presented in Table 3 and Figure 11. Meanwhile, the circular correlation time was fixed as = 20 ms to eliminate the effects of variable circular correlation time. With respect to the comparison in Table 3, the resampling frequency was reduced to about 10% of the conventional sampling frequency by applying the resampling strategy to Dataset 1, 2, and 3. In Figure 11, computations (multiplication and summation operations) of all the datasets are linearly related to the amount of frequency bins and the data size of the sampled signal. The amount of frequency bins was set as a constant, based on the width of frequency search range and the bandwidth of each frequency bin. And the data size of the sampled signal was exactly determined by the resampling frequency. Resulting from the reduced sampling frequency, the computation cost of signal acquisition with the resampling strategy was greatly decreased, nearly 90-94% of that without the proposed strategy, as shown in Table 3 and Figure 11. Besides, the received broadband signal in Dataset 4 was an 8-bit complex signal with negative intermediate frequency.
The bandpass filter of the proposed resampling strategy was removed since it was not applicable to extract the main lobe signal of Dataset 4. Additionally, the conventional sampling frequency of Dataset 4 was low enough. Thus, the resampling frequency reduced to about 59%, not as nearly 10% for other datasets. Correspondingly, the computation reduction of 43% was gained by applying the resampling strategy to signal acquisition of Dataset 4. Meanwhile, for all experimental datasets, the time cost of signal acquisition with the resampling strategy reduced to about 7.7~65% of that with the conventional acquisition algorithm, accompanying with the reduced computation. Table 3. Sampling frequency, computation, and time cost of signal acquisition without/with the resampling strategy for the received broadband signals of different experimental datasets. Figure 10. Sensitivity of the resampling strategy for weak signals. The red plot is the detection probability of signal acquisition with the resampling strategy, and the blue one is that of the conventional acquisition algorithm (without the resampling strategy). Experimental results in Table 3 and Figure 11 prove that the adoption of the resampling strategy directly resulted in the significant reduction on the resampled signals, which were used for FFT, correlation, and convolution operations. It makes the dominated contribution to the improved performance of computation and time cost of signal acquisition. Moreover, it is concluded that the resampling strategy can demonstrate better performance for the received broadband signal of higher conventional sampling frequency. Experimental results in Table 3 and Figure 11 prove that the adoption of the resampling strategy directly resulted in the significant reduction on the resampled signals, which were used for FFT, correlation, and convolution operations. It makes the dominated contribution to the improved performance of computation and time cost of signal acquisition. Moreover, it is concluded that the resampling strategy can demonstrate better performance for the received broadband signal of higher conventional sampling frequency.

Performance of Variable Circular Correlation Time
Variable circular correlation times are adopted for the coarse signal acquisition of all the datasets in the experiments to explore the relationship between circular correlation time and time cost of signal acquisition, as well as acquisition effectiveness.
For all the datasets, the number of satellites acquired without/with the resampling strategy is shown in Figure 12 when circular correlation time was varying from 10 ms to 20 ms . The experimental results are exhibited as the comparative red and blue plots to illustrate the effects of the resampling strategy. There are several cases, marked as green circles in Figure 12c,d, in which the number of satellites acquired with the resampling strategy was fewer than that with the conventional acquisition algorithm. However, there were still enough acquired satellites for navigation solution in these cases. Except that, the majority of experimental results with the resampling strategy were the same with that of the conventional one. This means that acquisition effectiveness of the resampling strategy can primarily reach that of the conventional acquisition algorithm, which cost about 10 times computation and time of the resampling strategy. Further, the number of satellites acquired without/with the resampling strategy nearly remained unchanged when circular correlation time was increasing from 12 ms to 20 ms. This indicates that there is

Performance of Variable Circular Correlation Time
Variable circular correlation times are adopted for the coarse signal acquisition of all the datasets in the experiments to explore the relationship between circular correlation time and time cost of signal acquisition, as well as acquisition effectiveness.
For all the datasets, the number of satellites acquired without/with the resampling strategy is shown in Figure 12 when circular correlation time was varying from 10 ms to 20 ms. The experimental results are exhibited as the comparative red and blue plots to illustrate the effects of the resampling strategy. There are several cases, marked as green circles in Figure 12c,d, in which the number of satellites acquired with the resampling strategy was fewer than that with the conventional acquisition algorithm. However, there were still enough acquired satellites for navigation solution in these cases. Except that, the majority of experimental results with the resampling strategy were the same with that of the conventional one. This means that acquisition effectiveness of the resampling strategy can primarily reach that of the conventional acquisition algorithm, which cost about 10 times computation and time of the resampling strategy. Further, the number of satellites acquired without/with the resampling strategy nearly remained unchanged when circular correlation time was increasing from 12 ms to 20 ms. This indicates that there is potential to shorten the circular correlation time so as to reduce computation and time cost, without obvious loss of acquisition performance. Whereas, too short circular correlation time reasonably leads to the incomplete acquisition of visible satellites, as cyan circles in Figure 12a,c. potential to shorten the circular correlation time so as to reduce computation and time cost, without obvious loss of acquisition performance. Whereas, too short circular correlation time reasonably leads to the incomplete acquisition of visible satellites, as cyan circles in Figure 12a,c. The time cost performance of variable circular correlation time on all experimental datasets is presented in Figure 13, when acquiring without/with the resampling strategy. Noticeably, when circular correlation time was varying from 10 ms to 20 ms, the time cost of signal acquisition with the resampling strategy (red bars) were much less than that of the conventional algorithm (blue bars), nearly decreased by 60-90% for Datasets 1, 2, and 3, while 5-35% for Dataset 4. This is attributed to the fact that the bandpass filter of the resampling strategy is not suitable for the signal with negative intermediate frequency in Dataset 4. The experimental results in Figure 13 illustrate the capability of the resampling strategy to decrease the time cost of signal acquisition. In addition, the blue and red plots in Figure 13 reflected the changes of time cost when the circular correlation time was varying. The time cost of signal acquisition without the resampling strategy seemed to increase plainly. Nevertheless, by linearly fitting these values, the time cost exhibited an increase of 3.6%, 3.6%, 3.8%, and 4.7% per millisecond for the 4 datasets, respectively, when the circular correlation time was varying from 10 ms to 20 ms. By contrast, the time cost for the resampling strategy looked stable, despite the existence of some outliers. Actually, the time cost of signal acquisition increased by 2.7%, 2.7%, 5.7%, and 2.7% per millisecond for the 4 datasets, respectively.
These experimental results indicate that a relatively shorter circular correlation time is beneficial for signal acquisition with both the proposed resampling strategy and the conventional acquisition algorithm, leading to less computation and time cost without obvious loss of acquisition performance. The time cost performance of variable circular correlation time on all experimental datasets is presented in Figure 13, when acquiring without/with the resampling strategy. Noticeably, when circular correlation time was varying from 10 ms to 20 ms, the time cost of signal acquisition with the resampling strategy (red bars) were much less than that of the conventional algorithm (blue bars), nearly decreased by 60-90% for Datasets 1, 2, and 3, while 5-35% for Dataset 4. This is attributed to the fact that the bandpass filter of the resampling strategy is not suitable for the signal with negative intermediate frequency in Dataset 4. The experimental results in Figure 13 illustrate the capability of the resampling strategy to decrease the time cost of signal acquisition. In addition, the blue and red plots in Figure 13 reflected the changes of time cost when the circular correlation time was varying. The time cost of signal acquisition without the resampling strategy seemed to increase plainly. Nevertheless, by linearly fitting these values, the time cost exhibited an increase of 3.6%, 3.6%, 3.8%, and 4.7% per millisecond for the 4 datasets, respectively, when the circular correlation time was varying from 10 ms to 20 ms. By contrast, the time cost for the resampling strategy looked stable, despite the existence of some outliers. Actually, the time cost of signal acquisition increased by 2.7%, 2.7%, 5.7%, and 2.7% per millisecond for the 4 datasets, respectively.
These experimental results indicate that a relatively shorter circular correlation time is beneficial for signal acquisition with both the proposed resampling strategy and the conventional acquisition algorithm, leading to less computation and time cost without obvious loss of acquisition performance. Sensors 2018, 18, x FOR PEER REVIEW 19 of 21 Figure 13. Time cost of signal acquisition with/without the resampling strategy for variable circular correlation time. Red bars and plot represent acquisition results with the resampling strategy, while blue ones are for that of the conventional acquisition algorithm.

Conclusions
The resampling strategy and variable circular correlation time are proposed to decrease computational complexity and time consumption of signal acquisition for GNSS receivers. The resampling strategy is inspired by the bandpass sampling theory and is applicable for broadband multi-frequency signals of widely-used multimode multiband GNSS receivers. The realization of the proposed resampling strategy is described in detail for combination with conventional acquisition algorithms. Specific to the much longer code period of satellite signals, variable circular correlation time is investigated to explore the potential of shortening the circular correlation time without obvious loss of acquisition performance. Moreover, the signal of the pilot channel is utilized to acquire the high-accuracy carrier frequency in the fine acquisition process. The acquisition threshold is set as the ratio of the highest and second highest correlation results in the search space of carrier frequency and code phase. Besides, the sensitivity of the resampling strategy to weak signals is investigated by extensive Monte Carlo simulations, and the computational complexity of signal acquisition is formulated by amounts of multiplication and summation operations to evaluate the efficiency of the signal acquisition algorithms.
Four sets of real GPS L2C signals are applied by comparative experiments to verify the effectiveness and efficiency of the proposed resampling strategy and variable circular correlation time. It has been proven that the resampling strategy has effectively decreased the computation and time cost of signal acquisition by nearly 90-94% without obvious loss of acquisition performance, although the sensitivity of the resampling strategy had a loss of nearly 1 dB and detection probability reached 100% when carrier-to-noise ratio was up to 38 dB-Hz. The greater decrease is achieved for the received broadband signal with the higher intermediate frequency and conventional sampling frequency. Besides, with circular correlation time varying from 10 ms to 20 ms, the time cost of signal acquisition increased by about 2.7-5.6% per millisecond, and the number of acquired satellites has no obvious changes for the cases with the resampling strategy and the conventional acquisition algorithm.

Conclusions
The resampling strategy and variable circular correlation time are proposed to decrease computational complexity and time consumption of signal acquisition for GNSS receivers. The resampling strategy is inspired by the bandpass sampling theory and is applicable for broadband multi-frequency signals of widely-used multimode multiband GNSS receivers. The realization of the proposed resampling strategy is described in detail for combination with conventional acquisition algorithms. Specific to the much longer code period of satellite signals, variable circular correlation time is investigated to explore the potential of shortening the circular correlation time without obvious loss of acquisition performance. Moreover, the signal of the pilot channel is utilized to acquire the high-accuracy carrier frequency in the fine acquisition process. The acquisition threshold is set as the ratio of the highest and second highest correlation results in the search space of carrier frequency and code phase. Besides, the sensitivity of the resampling strategy to weak signals is investigated by extensive Monte Carlo simulations, and the computational complexity of signal acquisition is formulated by amounts of multiplication and summation operations to evaluate the efficiency of the signal acquisition algorithms.
Four sets of real GPS L2C signals are applied by comparative experiments to verify the effectiveness and efficiency of the proposed resampling strategy and variable circular correlation time. It has been proven that the resampling strategy has effectively decreased the computation and time cost of signal acquisition by nearly 90-94% without obvious loss of acquisition performance, although the sensitivity of the resampling strategy had a loss of nearly 1 dB and detection probability reached 100% when carrier-to-noise ratio was up to 38 dB-Hz. The greater decrease is achieved for the received broadband signal with the higher intermediate frequency and conventional sampling frequency. Besides, with circular correlation time varying from 10 ms to 20 ms, the time cost of signal acquisition increased by about 2.7-5.6% per millisecond, and the number of acquired satellites has no obvious changes for the cases with the resampling strategy and the conventional acquisition algorithm.