1. Introduction
The BOC modulation has been adopted in modernized global positioning system (GPS), European Galileo project, and China’s BeiDou-3. BOC signals are characterised by a rectangular subcarrier, which split the signal spectrum around two main lobes at ±
fsc, the subcarrier repetition frequency. It can improve the tracking performance, and enhance spectrum compatibility due to the narrow ACF, and split spectrum [
1]. More detail and characteristics of BOC modulation can be found in literature [
2,
3,
4,
5].
However, the acquisition process of BOC signals becomes more complex, due to the multiple peaks in the ACF envelope, which requires a large number of timing hypotheses to search a given uncertainty window [
6]. Thus the acquisition becomes more computationally expensive and time consuming. Another problem in BOC acquisition is the high Nyquist sampling frequency [
7], which increases the computational burden and power consumption.
Some techniques have been reported to remove the impact of the secondary peaks in the ACF. The subcarrier phase cancellation technique (SPCT) [
8] introduces a quadrature squared subcarrier, called QBOC, to eliminate the ambiguity problem. In the autocorrelation side-peak cancellation technique (ASPeCT) proposed in [
9], subtraction of the cross-correlation between the BOC and PRN signals from the BOC autocorrelation is employed; but the technique is applicable only to sine-BOC(
n,
n). The Pseudo-Correlation Function (PCF) [
10] method constructs a single-peak ACF using a special designed reference signal. However, it will suffer from severe detection performance degradation for high-order BOC signal. The techniques mentioned above can remove the impact of the secondary peaks in the ACF, but they all process the upper and lower sidebands together, resulting in a high sampling frequency and computational complexity.
The main idea behind the BPSK-like method is that the BOC signal can be obtained as the sum of two BPSK signals, located at positive and negative sub-carrier frequencies. The effect of sub-carrier modulation can be removed by using a pair of single-sideband correlators. A receiver may use a single-side band (SSB), either the upper or lower sidebands, or use dual-side band (DSB), where both bands are combined non-coherently. Due to filtering and correlation losses, the BPSK-like methods bring some degradation. The losses are 3dB for SSB and about 0.5 dB for DSB [
11]. However, the BPSK-like method needs two RF channels, mixers and analog-to-digital converters (ADC), which increase the implementation complexity [
11].
In this paper, the band-pass sampling technique [
12,
13,
14,
15,
16] is introduced in BOC signal sampling to remove the secondary peaks, as well as to reduce the sampling frequency. Band-pass sampling is the technique of under-sampling a band-pass signal to achieve a frequency translation via intentional aliasing [
17]. If the sampling frequency is
, then the spectrum of the sampled signal can be obtained by replicating the spectrum of the original signal at multiples of
[
18]. Since a BOC signal can be treated as the sum of two BPSK signals with adjacent center frequencies [
19], the band-pass sampling technique for multiple signals [
14,
17,
20] can be used. Furthermore, taking into account the subcarrier Doppler, a dual-loop based acquisition structure is proposed to acquire the sampled signal.
The remainder of this paper is organized as follows.
Section 2 introduces the dual-sideband model of BOC signals, which is the basis of band-pass sampling; in
Section 3, the band-pass sampling is applied and the sampling parameters are derived. The dual-loop based acquisition structure is also described in detail. The detection probability and computational complexity are analysed in
Section 4. The conclusions are drawn in
Section 5.
2. Dual-Sideband Model of BOC Signals
The BOC(
m,
n) signal can be expressed as:
where
A is the signal amplitude;
is the navigation data;
is the pseudorandom noise (PN) code with code rate
fc =
m × 1.023 MHz;
is the squared subcarrier with frequency
fsc =
n × 1.023 MHz;
and
are the carrier frequency and phase, respectively. Generally, the integration period is shorter than the bit duration. So we will ignore the navigation data, that is,
.
The squared-wave subcarrier can be expanded in Fourier series [
19] (for sine-phased BOC signals):
where
is the Fourier coefficients;
is the subcarrier period.
Equation (2) indicates that the squared-wave subcarrier can be expressed as the sum of a series of sine waves. In practice, both the transmitter and the receiver are band-limited, which means that the high-frequency component of the signals will be filtered out. So in band-limited cases, the squared-wave subcarrier can be approximated by:
where
N is the number of terms preserved by front-end filtering.
For high-order BOC signals, the front-end filter may remove all the harmonic waves. Then the subcarrier can be approximated by a sine wave:
As a result, the BOC(
m,
n) signal can be approximated by
Equation (5) indicates that in band-limited cases, high-order BOC signals can be treated as a dual-sideband signal.
3. BOC Acquisition Based on Band-Pass Sampling
3.1. Problems in BOC Acquisition
As a result of the subcarrier, the ACF of BOC signals has multiple peaks. The squared magnitude ACF is directly relevant to the class of non-coherent energy detection acquisition receivers considered here.
Figure 1a illustrates the normalized squared magnitude ACFs of BOC(15,2.5) and BPSK(2.5). As shown in this figure, the squared magnitude ACFs of BOC(15,2.5) has as many as 23 peaks, while BPSK(2.5) has a unique peak.
This indicates that the step of the searching time bin should be sufficiently small to be able to detect the main peak of the ACF (i.e., we need a higher number of timing hypotheses in order to search a given time-uncertainty window compared to BPSK case) [
6]. Another problem is that the bandwidth of BOC(
m,
n) signals is much wider than BPSK(
n) signals, so that the Nyquist frequency is very high, which results in a higher computational complexity.
3.2. Band-Pass Sampling Technique
As mentioned before, the secondary peaks and the high sampling frequency of high-order BOC signals result in high complexity. Band-pass sampling technique can be used to remove the secondary peaks and reduce computational complexity.
For a single-band signal, the range of
is given by [
12]:
where
and
are the lower and upper frequencies, respectively;
is an positive integer.
The frequency of the sampled signal is [
17]:
where
is the intermediate frequency of the continuous signal; fix(
a) is the truncated portion of argument
a and rem(
a,
b) is the remainder after division of
a by
b.
For a dual-sideband signal, the range of
is:
The signal bands must not overlap in the frequency spectrum of the resultant sampled bandwidth. This can be expressed mathematically for two signals as [
17]:
where
and
are the center frequencies of the sampled signals;
and
are the bandwidth of the two signals, respectively.
3.3. Band-Pass Sampling for High-Order BOC Signals
Since the filtered high-order BOC signal can be treated as a dual-sideband signal, the band-pass sampling technique can be applied. In this study, the signal will be sampled at intermediate frequencies, and the intermediate frequency is an important parameter.
If the intermediate frequency of BOC(
m,
n) is
, then the center frequencies of the two sidebands are:
Define and . To simplify the analysis, we choose the sampling frequency to make sure is even. Then the spectrum transition will be discussed for two cases.
Case a:is odd.
In this case, the frequencies of the sampled signals are:
where
and
.
The spectrum transition is depicted in
Figure 2. Since the signal is real, there is always a negative spectrum.
As shown in
Figure 2, the spectrum marked
is actually from
, and it is an inverse version of
. So, if we treat
and
as the target sampled signal, then
and
are the negative spectrum. For the sake of clarity, we will change the sign of
, then (10) should be rewritten as:
Case b: is even.In this case, the frequencies of the sampled signals are:
The spectrum transition is depicted in
Figure 3.
In summary, the frequencies of the sampled signals are:
where
and:
Consequently, given the center frequencies of sampled sidebands, the sampling frequency and the intermediate frequency can be determined as follows:
where
is the distance between the two sidebands of the sampled signal;
is the intermediate frequency of the sampled signal;
is the bandwidth of either sideband.
Equation (14) can be simplified as:
The sampling frequency must also satisfy (8).
Taking Galileo E1 PRS signal (BOC(15,2.5) modulated) as an example,
Figure 4 illustrates the power spectrum density (PSD) and the ACF of the sampled signal (
,
). The sampling frequency and the intermediate frequency are
and
, respectively.
Figure 4 shows that the PSD and ACF of the sampled signal are similar to that of a BOC(
n,
n) signal and the number of the remaining secondary peaks is reduced to 2, much less than that of BOC(15,2.5).
3.4. Acquisition of Band-Pass Sampled BOC Signals
Taking the velocity into consideration, the sampled signal can be expressed as:
where
fd is carrier Doppler;
fdsc is subcarrier Doppler;
Ts is the sampling interval.
Given
as the relative speed between the satellite and the receiver, the Doppler frequencies of the lower band and upper band are:
Then the Doppler frequencies of the carrier and subcarrier are:
Equation (16) shows that the sampled signal is still a BOC signal, with subcarrier frequency and intermediate frequency .
In the acquisition process, the receiver must determine the carrier frequency offset, the subcarrier frequency offset and the time offset. Then the acquisition turns into a 3-dimensional searching process. However, since the ratio between and always equals the ratio between f0 and fsc; that is, , the acquisition is actually a 2-dimensional searching process.
An acquisition structure based on dual-loop is proposed in
Figure 5.
In
Figure 5,
. The correlators’ outputs are:
where
;
is the ACF of the PN-code;
is the time offset;
and
are the frequencies of the subcarrier numerically controlled oscillator (NCO) and carrier NCO;
and
are the phases of the subcarrier NCO and carrier NCO;
and
are the Doppler errors of the subcarrier and carrier;
and
are the phase errors of the subcarrier and carrier;
,
,
,
are independent zero-mean Gaussian variables.
Then the decision statistic is given by: